Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Modelling the Pleistocene colonisation of Eastern Mediterranean islandscapes

  • Theodora Moutsiou ,

    Contributed equally to this work with: Theodora Moutsiou, Christian Reepmeyer, Athos Agapiou

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing

    tmouts01@ucy.ac.cy

    Affiliations Archaeological Research Unit, University of Cyprus, Nicosia, Cyprus, ARC Centre of Excellence for Australian Biodiversity and Heritage, College of Arts, Society and Education, James Cook University, Cairns, Australia, Department of Archaeology, Max Planck Institute for the Science of Human History, Jena, Germany

  • Christian Reepmeyer ,

    Contributed equally to this work with: Theodora Moutsiou, Christian Reepmeyer, Athos Agapiou

    Roles Data curation, Investigation, Writing – review & editing

    Affiliations Archaeological Research Unit, University of Cyprus, Nicosia, Cyprus, ARC Centre of Excellence for Australian Biodiversity and Heritage, College of Arts, Society and Education, James Cook University, Cairns, Australia, Department of Archaeology, Max Planck Institute for the Science of Human History, Jena, Germany, Science and Technology in Archaeology and Culture Research Center, The Cyprus Institute, Nicosia, Cyprus

  • Vasiliki Kassianidou ,

    Roles Funding acquisition, Supervision, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliation Archaeological Research Unit, University of Cyprus, Nicosia, Cyprus

  • Zomenia Zomeni ,

    Roles Data curation, Methodology, Resources, Validation, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliation Geological Survey Department, Lefkosia, Cyprus

  • Athos Agapiou

    Contributed equally to this work with: Theodora Moutsiou, Christian Reepmeyer, Athos Agapiou

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – review & editing

    Affiliations Department of Civil Engineering and Geomatics, Faculty of Engineering and Technology, Cyprus University of Technology, Limassol, Cyprus, Eratosthenes Centre of Excellence, Limassol, Cyprus

Abstract

Predictive models have become an integral part of archaeological research, particularly in the discovery of new archaeological sites. In this paper, we apply predictive modeling to map high potential Pleistocene archaeological locales on the island of Cyprus in the Eastern Mediterranean. The model delineates landscape characteristics that denote areas with high potential to unearth Pleistocene archaeology while at the same time highlighting localities that should be excluded. The predictive model was employed in surface surveys to systematically access high probability locales on Cyprus. A number of newly identified localities suggests that the true density of mobile hunter-gatherer sites on Cyprus is seriously underestimated in current narratives. By adding new data to this modest corpus of early insular sites, we are able to contribute to debates regarding island colonisation and the role of coastal environments in human dispersals to new territories.

Introduction

Recent advances in the study of island colonisation, which showed the scale of hunter-gatherer interaction with maritime resources and island landscapes [16], require us to revisit the question of a potential Pleistocene human presence on Cyprus. Unequivocal evidence of material culture of this period is currently largely absent from Cyprus, but recent research has shown that an earlier colonisation of the island is viable [7]. Moreover, the current absence of evidence might relate to missing research focus on earlier time periods rather than representing an archaeological reality. This paper presents a methodological approach that combines traditional archaeological practices and technological applications from the field of geoinformatics as a tool for the investigation of a potential Pleistocene colonisation of the island of Cyprus. Although our emphasis is the Late Pleistocene (last 125 ka), we envisage that such an approach will be a useful means in the identification of archaeological sites of other periods too, particularly in areas, where no archaeological sites have been found so far.

Here we report on a study designed to systematically investigate high potential localities for Late Pleistocene human presence either for purposes of habitation or exploitation of subsistence resources (e.g. food, raw materials). The geospatial modelling of target areas is presented here functioning as the starting point for a large-scale research-based field survey aimed at collecting empirical archaeological data that will feedback into the model. A simple set of assumptions regarding Palaeolithic decision-making underpin this study; fundamentally, we assume that although limited concrete evidence for Pleistocene archaeology exists on the island, multiple Pleistocene mobile hunter-gatherer groups could have entered Cyprus, either on short-term forays or longer stays. In addition to hunting medium-sized game and using good quality chert for stone tool making, we assume that they were familiar with coastal and marine resources, also present on Cyprus. This assumption stems from recent studies that document a generalised subsistence system, including coastal foods, for early humans [8, 9]. Thirdly, we assume that Pleistocene hunter-gatherers could have used all regions of the island, including mountainous areas. We base this assumption on the fact that Terminal Pleistocene and Early Holocene sites on Cyprus with hunter-gatherer-forager elements are found located not only coastally but also inland and upland, such as Vretsia Roudias [10] and Krittou Marottou Ais Giorkis [11].

Predicting locales of past human activity

Predictive modeling is a method used to identify archaeological site locations on the basis of observed patterns and assumptions about human behaviour [1215]. It was initially developed in the mid-1970s and mainly applied within the context of governmental land management and Cultural Heritage management [12, 16, 17]. Early predictive models have focused on the application of multivariate statistical techniques to discuss site and settlement patterns by extrapolating from sample to population [12]. More recently, digital predictive methodologies combining Geographical Information Systems (GIS) and geostatistical methods have been used to estimate the probability of encountering archaeological sites outside of areas where they have already been found in the past. Spatial analysis and predictive modeling using GIS have been subjected to fierce debate regarding the suitability of such approaches in archaeological research. Predictive models have been heavily criticised as leading to environmental determinism in archaeology [1821] and often considered as unable to predict the location of all archaeological sites [22, 23]. Despite its shortcomings, predictive modeling has proved a valuable tool for the rescue of archaeological remains and data that would have otherwise been lost due to modern development. Beyond Cultural Heritage management applications, predictive modeling has already been successfully exploited in spatial planning of archaeological research [24, 25]. Predictive models have resulted in particularly useful insights especially in areas where there is an extreme paucity of data that could be used as a baseline for expanding archaeological work or in cases where extensive areas needed to be covered with traditional surveys [2628]. In such instances, predictive models can assist in the determination of which regions are more likely to produce archaeological localities and which regions are less likely to do so.

In our study, the focus is predominantly on low-density lithic scatters—the main material remains of prehistoric hunter-gatherer activity. Such archaeological signatures are largely ephemeral with low visibility and, hence, particularly difficult to identify in the landscape (see [7] for a detailed discussion). It is of paramount importance, then, to pay a closer look at the insular landscape in order to understand how Pleistocene hunter-gatherers would have inhabited and used their environment. Here we combine environmental and socio-cultural factors, focusing on features of the landscape that would have had socio-cultural significance for early hunter-gatherers, such as raw material sources.

Regional setting

Cyprus is the third largest island in the Mediterranean extending over 9,251 km2 and at distances around 30–70 km (depending on marine transgressions) from the nearest coastline. It is also a resource-rich island, encompassing good quality lithic sources and an ecosystem that could sustain hunter-fisher-gatherer groups [7]. Despite all this, however, Cyprus has remained at the periphery of early island colonisation discussions mainly due to an ongoing assumption of its marginality due to open water gaps functioning as barriers to pre-Holocene maritime crossings (see discussion in [29]). Global interest in early seafaring and maritime colonisation in recent years [1, 2, 30] has led to significant advances in the thematic and there is tentative evidence that Homo sapiens as well as pre-sapiens species were capable of maritime crossings and successfully colonised insular environments in various parts of the world (see [27, 31, 32] for some eastern Mediterranean examples).

Archaeological background

Akrotiri Aetokremnos [33, 34] in southern Cyprus is currently the oldest securely dated site on the island with dates that place its occupation at the Terminal Pleistocene (ca. 12 ka). The coastal sites of Nissi Beach and Akamas Aspros [35], and Vretsia Roudias [10] and Ayia Varvara Aspokremmos [36, 37] further inland, add to the number of likely Early Holocene sites with hunter-gatherer/forager attributes. The first permanent settlements, for example Ayios Tychonas Klimonas [38] also appear during the Early Holocene (ca. 10 ka) followed by a Neolithic settlement expansion well into the Middle Holocene (ca. 8 ka).

Beyond these excavations, only small collections of surface lithic finds of a tentative Palaeolithic character [7] have been documented from various parts of the island during the 1960s-1980s (Fig 1). Stockton ascribes an Upper Palaeolithic character to a lithic assemblage, including cores, flakes and formal tools (86 lithics), which he collected from four different localities at the Kyrenia range, on the northern coast of Cyprus, near the area of Khrysokava during a short visit to the island in 1967. In 1968, during a brief visit to the island, Vita-Finzi [39] collected a small assemblage of lithic artefacts from a fossil beach near the mouth of the Maroni River at Zygi (Larnaka District), southeast Cyprus. The artefacts were assigned to the Middle Palaeolithic, based on their typo-technological characteristics and geological context. The archaeological reconnaissance of the Khrysochou Bay area project [40] and the Tremithos Valley Project [41] targeting the Khrysochou drainage near Polis (Paphos District), in northwest Cyprus, and the Tremithos Valley (Larnaka District), southern Cyprus, respectively, documented amongst other material culture, lithic collections with typo-technological characteristics of a potential Palaeolithic age. Finally, in 2016 Strasser et al. [42] ascribed a Lower Palaeolithic age to a biface, a surface find which was collected by S. Swiny in 1992 near the Aceramic Neolithic site of Kholetria Ortos (Pafos District), southwest Cyprus [43], after examining a three-dimensional printed replica of the artefact. The significance of these finds in demonstrating a Pleistocene human presence in Cyprus, either in the form of ephemeral visits or a more permanent habitation, is, nevertheless, limited. Surface finds are unreliable age markers, especially when dealing with extremely small numbers. Additional data are needed to add to this modest corpus of information and fast-track research into the Pleistocene past of the island.

thumbnail
Fig 1. Terminal Pleistocene sites and legacy Palaeolithic surface lithic assemblages from Cyprus.

Map depicting known Terminal Pleistocene archaeological sites and localities where legacy Palaeolithic lithic assemblages (with typo-technological characteristics ranging from late Middle to Late Pleistocene) have been collected across Cyprus used in the development of the predictive model. Compiled in ArcGIS with Digital Elevation Model (DEM) data from the Land and Surveys Department of the Republic of Cyprus. Original copyright with the authors.

https://doi.org/10.1371/journal.pone.0258370.g001

Materials and methods

A key problem in predictive modeling is the availability of representative archaeological data that can be used either as input, or as a test for a set of expectations [44]. Cyprus has a limited dataset related to Pleistocene archaeology. To counterbalance this paucity, we have utilised data from better-studied regions within the Eastern Mediterranean region. Essentially, we correlate the locations of known sites with landscape and topographic features and use this information to model locations of other unknown sites in an area. The list of sites used in the analysis is included in S1 Table. Although the actual number of archaeological sites in these areas might be underestimated, we consider these as representative site samples sufficient for our model. Limitations of this approach are discussed in the appropriate section below.

Variables

Here we consider parameters on the Cypriot landscape that would have been attractive to early humans: freshwater sources, raw material sources and topographical features. Freshwater sources apart from being vital to human survival would have attracted fauna also necessary for human subsistence; sources of good quality raw materials are essential for the manufacture of a reliable tool-kit. Topographical features of relevance include slope, aspect and elevation. Mobile hunter-gatherers of the Late Pleistocene would have been familiar with coastal resources and marine foods [9] but not restricted to such environments even when on islands as often assumed. In order to account for this, we also consider distance from the (modern) coast and fresh water sources, such as rivers and springs. We considered all these variables initially independently from each other to determine whether any given variable is preferentially selected. Once it became clear that we do not have any distinguishable patterns, we used these variables together in developing our model to avoid biases in estimates of importance (see also [12]). The model was developed using various methods: Boolean geometry, proximity estimation and similarity analysis. We performed our GIS analyses based on two spatial input-data variables: local and regional.

Regional-based data analysis.

We collected data from various published sources [4548], including material from surveys and excavated horizons from mainland Greece, the Ionian and Aegean Seas, Turkey, Israel, Syria, Jordan and Lebanon (S1 Table). Data and maps from these reports were georeferenced and then digitised in a GIS environment for further analysis.

Local-based data analysis.

On a local-based data analysis, which encompasses the full extent of the island of Cyprus, we concentrated on the localities of the documented legacy collections of a potential Palaeolithic character. Furthermore, taking into consideration that hunter-gatherers exploited Cyprus at the end of the Pleistocene, we opted to include known Epipalaeolithic sites dating to the Terminal Pleistocene in our analyses (Fig 1).

Data

The geospatial analyses concerning the broader Eastern Mediterranean relied upon global Digital Elevation Models (DEM) from the Shuttle Radar Topography Mission (SRTM) Digital Elevation Database. For the geospatial analyses concerning Cyprus we make use of the DEM, acquired from the Department of Land and Surveys of Cyprus (DLS). In addition to the DEMs, other geographical data such as major rivers in the Eastern Mediterranean region and the drainage maps of Cyprus have been used.

The model

The predictive model was created using a combination of layers for elevation, slope, aspect, distance from rivers/springs and the coast, as well as lithology and general geology on the local scale. In choosing these parameters, we consider the natural environment but also cognitive and behavioural aspects of Pleistocene hunter-gatherers. As such, our approach follows a theoretically informed data-driven modeling. We expect our analyses to result in a robust model for the investigation of Pleistocene localities in Cyprus but there remain some limitations. The most obvious one is inherent in any GIS approach, namely that in the application of geospatial tools, we are estimating past landscapes with features we can only quantify in the present (for a discussion, see the ‘complex topography concept’ [49, 50]). Moreover, remotely sensed data can only record the present surface and very shallow subsurface topography.

In coastal regions the onset and offset of 60+ ice ages, which characterise the Quaternary led to extreme sea-level changes (as much as 120 m sea-level drop during ice ages and as much as 6m sea-level rise during interglacial periods, both in relation to present sea-level), resulted in substantial impacts to near-shore landscape morphology. Phenomena such as tectonics, isostacy and changes in erosion patterns have all contributed to the Quaternary being a geologic period with extreme landscape changes, at least as it pertains to humans, their evolution and behaviour in relation to the environment [51, 52].

In the first iteration of the model, our input-data were based on a regional scale (Fig 2) and, particularly, on topographic variables in order to determine whether any clear patterns regarding the location of Pleistocene sites within the broader Eastern Mediterranean emerge. For this, we analysed over 100 archaeological localities from the circum-Cyprus region, specifically looking at their elevation, slope, aspect and distance from rivers and the sea. The attributes of all the analysed Eastern Mediterranean sites fed back to the predictive model used for Cyprus.

thumbnail
Fig 2. Eastern Mediterranean Pleistocene sites.

Map showing Pleistocene archaeological sites in the circum-Cyprus Eastern Mediterranean region. A detailed list of the sites and references can be found in S1 Table. Original copyright with the authors.

https://doi.org/10.1371/journal.pone.0258370.g002

In the application of the model to Cyprus, we tested the results of the regional model on a local scale and incorporated information from legacy Palaeolithic surface collections and known hunter-gatherer, Epipalaeolithic, sites on Cyprus. Fifteen sites were plotted as shown in Fig 1 and Table 1. Their environmental attributes where then extracted and analysed. At this stage, we also incorporated additional variables, specifically smaller fresh water sources (springs), lithology and geology. In detail, the distance from coastline, rivers, elevation, aspect and slope was examined, as well as lithology and geological attributes. Initially, each parameter was evaluated separately to investigate their importance for the creation of a predictive model. However, this did not provide any clear patterns, and therefore a multi-variable analysis was processed using Boolean algebra. In addition, an alternative approach was followed: a probability map, based on one-class classification analysis, using the Spectral Angle Mapper (SAM) classifier was performed. The results eventually provided areas where there is a higher probability to share similar attributes with the ones of the fifteen earliest known archaeological locales of Cyprus. At the last stage, we also tried to limit the areas of future investigations taking into consideration the Quaternary (terrestrial) geology of the island ignoring the Holocene and focusing on the Pleistocene, including chert-bearing formations that would have provided raw materials for tool manufacture (see maps below).

thumbnail
Table 1. Main topographic features of Terminal Pleistocene and legacy Palaeolithic surface lithic assemblages from Cyprus.

https://doi.org/10.1371/journal.pone.0258370.t001

Field methodology

Field methodology followed a classic reconnaissance survey approach with a small team of experts targeting high probability areas as delineated by the predictive model. Permits for this fieldwork were acquired in written form by the Department of Antiquities of the Republic of Cyprus. The primary objective of the field surveys was to systematically search specific regions for remains of human activity and to collect surface archaeological samples to the validity of the predictive model and to demonstrate function and date of human activity. Between September and October 2020, the research team conducted systematic by-foot surface surveys, focusing on the areas/zones delineated by the model as ‘hotspots’ and especially targeting specific topographic features (see ‘Local’ Results section) in these zones, such as rockshelters, riverbanks and sediments of a known Pleistocene age, as more likely to preserve in situ material. All investigated areas/tracks were documented via tracklogs with GPS points and detailed soil and landscape descriptions (e.g. average surface visibility, present land use, and vegetation) recorded for localities that preserved cultural remains. Surface sampling was conducted in order to determine the function and age of the artefacts and the use of different parts of the area in different periods. A spatially well-defined concentration of lithic cultural materials (stone tools) was designated a findspot (lithic scatter). In areas of high interest, sediment depth was investigated through probing and augering in order to record points for possible future test excavations. At the end of the research, all the new data that have been georeferenced were collated in density/distribution maps.

Results

Regional

The main aim of our model is to delineate a zone of high probability to unearth Pleistocene archaeological remains on the island of Cyprus that is as precise as possible. We started by modeling the main characteristics of known Pleistocene localities within the Eastern Mediterranean, expecting to define potential topographical patterns that may be used in delineating a probability zones on Cyprus. Fig 3A–3D summarises the main topographic patterns of the region. Slope, aspect and elevation do not appear to be significant factors in the location of Pleistocene sites in the Eastern Mediterranean. Aspect is particularly non-informative and elevation similarly so. Two thirds of the sites are at an elevation that reaches or exceeds the highest mountain ranges of Cyprus. Slope is somewhat more informative as it suggests a preference for locations with a slope component <20°, with the majority of the sites found at <15°. Regarding distances from water, the majority of the sites appear to be within 10 km (maximum distance) from either a river or the coastline. Sites are predominantly located closer to a river source than the coast. Sites found at a great distance from rivers are probably supplied with fresh water by smaller water sources not included in the available dataset. Overall, the regional scale analyses highlight the significance of freshwater proximity and the secondary nature of topographic features, especially aspect, in site location preferences.

thumbnail
Fig 3. Main topographic patterns of Eastern Mediterranean Pleistocene archaeological sites.

Graphs (a = distance from coast/rivers, b = elevation, c = slope, d = aspect) summarising the main topographic patterns of known Pleistocene archaeological sites across the Eastern Mediterranean.

https://doi.org/10.1371/journal.pone.0258370.g003

Local

In our next step we focused on Cyprus itself exploring patterns of Terminal Pleistocene archaeological sites and surface assemblages with typo-technological characteristics of a probable Palaeolithic age known from across the island (Table 2). Here too, aspect was found not to be informative. With regards to elevation most of the sites are located <350 m with 46.6% of these occurrences ranging between 0–100 m. Slope largely reflects the patterns observed in the regional model with a preference for locations with a slope component <20°, although here the majority of sites is found at 0–10°. In terms of distances from the coast we note that sites range from <100 m to > 5000 m from the modern coastline, with almost half of the sites (46.6%) located up to 500 m from the coast. Considering distance from freshwater sources, our analysis shows that 60% of the sites are found between 500–5000 m from a river source, with only 13.3% of the sites located within 500 m. It is known that tectonic and/or climatic events can alter river courses over millennia; however, here we have concentrated on major rivers with a well-established Quaternary evolutionary history. It is possible that the existence of smaller tributaries or perennial/seasonal streams that were not included in the analysis skews the actual distance from freshwater sources for these early sites. Finally, if we consider the maximum distance of these early sites from either the coastline or rivers and springs, we note that most of these archaeological localities are located within 5 to 10 km from either water source (Figs 46).

thumbnail
Fig 4. Prediction map of high probability areas for Pleistocene archaeological sites on Cyprus.

Prediction map delineating high probability areas for the location of potential Pleistocene archaeological sites on Cyprus based on attractive topographic attributes. The dark red dots represent known early archaeological sites on Cyprus (see Fig 1 for site names). Note how all but one archaeological sites fall within the high probability zone (green area) as defined by the predictive model. Original copyright with the authors.

https://doi.org/10.1371/journal.pone.0258370.g004

thumbnail
Fig 5. Probability map based on one-class classification analysis.

Probability map based on one-class classification analysis, using the Spectral Angle Mapper (SAM) classifier. The map illustrates areas where there is a higher probability to share similar attributes with the ones of the fifteen earliest known archaeological sites of Cyprus. The dark red dots represent known early archaeological sites on Cyprus (see Fig 1 for site names). Original copyright with the authors.

https://doi.org/10.1371/journal.pone.0258370.g005

thumbnail
Fig 6. Prediction analysis final output.

The map illustrates the final output of the predictive analyses, taking into account the combination of slope, elevation, aspect, distance from coast and rivers results as well as Pleistocene terrestrial geology and raw material availability. The map delineates the smallest possible area of high potential to preserve Pleistocene archaeology based on the currently available data (green zone). The dark red dots indicate the location of known early archaeological sites on Cyprus (see Fig 1 for site names). Original copyright with the authors.

https://doi.org/10.1371/journal.pone.0258370.g006

thumbnail
Table 2. Analysis outcomes of the main topographic patterns of Terminal Pleistocene and legacy Palaeolithic surface lithic assemblages from Cyprus.

https://doi.org/10.1371/journal.pone.0258370.t002

We acknowledge that in these analyses we rely on the modern coastline and, presumably, archaeological localities may currently be under water as a result of post-glacial sea level changes. High-resolution bathymetric data are not readily available for the submerged landscape of Cyprus, limiting our understanding of the island’s underwater relief and pre-Holocene land extension. However, considering that our research timeframe is the last 125 ka, we should revisit the issue of presently submerged landscapes that were terrestrial landscapes during glacial periods. About 110 ka and then again 90 ka sea level had dropped as much as 50 m below present sea level. About 65 ka sea level had dropped as much as 90 m below present sea level, but the most dramatic drop is during the Last Glacial Maximum (about 22 ka) when sea level dropped as much as 120 m below present sea level. These presently submerged terrains consist of a shelf that can be fairly narrow or wide, depending on the shelf topography, the steeper the submerged coast topography, the more narrow the zone of interest. This zone is 1.5 km wide in the more precipitous areas and as much as 8 km wide in the more shallow shelfs. Most of the onshore analysed sites are found at a distance of around 10 km from the modern coast we assume that distance from the coast was not a controlling factor on location preferences for Pleistocene occupation of Cyprus. In their study of the Last Glacial Maximum in the Western Mediterranean, Barton et al. [53, 54] note that there is a general tendency for assemblages accumulating in base camps to predominate in a zone approximately 40–60 km from the contemporaneous coast, while ephemeral camps are located both seaward of the base camps, in a zone 0–15 km from the coast, and further inland, at 100–125 km from the coast. Thus, even though marine transgressions at the end of the Pleistocene may have altered distance of sites to the coast on Cyprus, some sites are expected to be located not directly associated with coastlines, but in the interior of the island.

Our analyses show that distance from a water source would have been an important parameter for site selection localities on behalf of Pleistocene hunter-gatherers. Indeed, on Cyprus, one can find the earliest securely dated water wells in the world (Early Holocene), which suggests that early populations on the island were aware of their availability and had good knowledge of how to best exploit the fresh water sources of the island [55]. Our knowledge of the spatio-temporal distribution of paleohydrological resources during episodes of climatic amelioration (increased/decreased rainfall) remains limited and for our analyses, we rely on modern day data. Despite these constraints and considering that permanent water sources require particular lithologies, hydrogeological environments and geomorphic settings [5659], it is reasonable to expect that the distribution of freshwater sources in the past is likely broadly analogous with their modern distribution [6064].

Another key element in site location choices would have been lithic raw material availability. Lithic resources play an essential role in Palaeolithic lifeways as they provide the raw materials for stone tool production for the fulfillment of various subsistence-related everyday tasks. On Cyprus, chert is the main raw material resource with performance characteristics–‘knappability’, sharpness and durability–and fracture mechanics ideal for tool manufacture [65]. We included lithology as a parameter in our model, focusing particularly on chert-bearing geological formations (Fig 7). Most of the Cypriot sites used in the model appear to be located within these resource-rich zones, although it is not presently possible to estimate the exact distances to the various sources. Rarely, a site is found outside a chert-bearing formation; this may indicate that raw material acquisition could also take place at a distance, possibly from secondary sources. For example, Akrotiri Aetokremnos is located at a distance from any of the region’s chert-bearing geological formations but Kouris River, one of the major rivers of Cyprus, is located only a short distance from the site, suggesting hunter-gatherers could have procured material not necessarily from the rock outcrops but from rich river deposits, such as the active channel and flood plain deposits.

thumbnail
Fig 7. Chert-bearing geological formations of Cyprus.

Original copyright with the authors.

https://doi.org/10.1371/journal.pone.0258370.g007

Field survey results

A preliminary research-based field survey has been conducted as a spatially explicit island-scale assessment of the opportunities for Pleistocene human occupation of Cyprus. The survey used initial outcomes of the GIS model and targeted appropriate landscape features in an endeavor to collect field archaeological data that will feed back into the model. Although the GIS model considers the island as a whole, field surveys are only possible in the southern part of the island, which is currently under the effective control of the Republic of Cyprus. Over two seasons, key localities–as determined by the geospatial predictive model–across the island were located and investigated via intensive foot survey. Reconnaissance surveys targeted both the coastal zone and the interior of the island, as delineated by the GIS model, while also focusing on specific topographic features, such as rockshelters, rivers and Pleistocene-aged deposits. Surface lithic scatters as well as individual artefacts (isolated finds) were identified and recorded in various localities (Fig 8) confirming the expectations of the geospatial model. Based on a preliminary study of their morphological-typological characteristics, these chert stone tools likely range chronologically from at least the Epipalaeolithic to the Neolithic (Terminal Pleistocene–Early Holocene). The artefacts, mainly flakes and blades, are made from different types of chert. The majority of the identified archaeological localities have been found on a coastal setting, reflecting similar patterns to those observed from the Cypriot legacy collections and excavated Terminal Pleistocene sites.

thumbnail
Fig 8. Archaeological localities identified during field surveys based on the predictive model outcomes.

Archaeological localities identified during field surveys superimposed over predictive map (green area). Insert shows the location of surveyed region on the island of Cyprus. Original copyright with the authors.

https://doi.org/10.1371/journal.pone.0258370.g008

Another important outcome of the field survey is the open-air character of the documented localities. Caves are not a common feature of the Cypriot landscape but several small rockshelters were noted during fieldwork. The team investigated those as they are assumed to preserve intact deposits. No evidence of prehistoric exploitation of those landscape features was recorded on any of the studied rockshelters. This outcome is an important negative confirmation of the predictive model highlighting its validity since most of the rockshelter sites are located on slopes >20% and are, hence, outside the model’s high probability zone. The identification of several lithic scatters of at least a Terminal Pleistocene age is significant as it adds to a modest corpus of data from Cyprus for this period. These sites begin to illustrate the true density of likely hunter-gatherer/forager sites on the island and show that mobile groups were exploiting and perhaps residing on Cyprus for several millennia, potentially earlier than the current archaeological record allows for.

Discussion: Modelling the Pleistocene exploitation of Cyprus

Decision-making relies on the active engagement of humans with their environments and includes any choice among alternative courses of action that usually involve a set of feasible alternatives/options [66]. For Pleistocene hunter-gatherers exploiting a new territory, crucial decisions will involve subsistence strategies, resource acquisition and home base locale preferences at the very basic level, which in turn affect social interactions. Forward planning and advanced organisational capabilities enabled Pleistocene hunter-gatherers, for example, to successfully colonise northern latitudes during full interglacial conditions in the Middle Pleistocene [67]. Planning ahead provides a degree of cognitive plasticity that must not be underestimated. Whatever the difficulties of adapting to a new territory or environmental context, the fact that Pleistocene hunter-gatherers elsewhere were able to flourish in all climatic conditions and habitats argues for a high degree of resourcefulness [6769].

The island of Cyprus has traditionally remained isolated from discussions concerning the pre-Holocene colonisation of islands. This is primarily a result of persisting narratives of the Mediterranean islands as marginal environments to a hunter-fisher-gatherer lifestyle. Archaeological discoveries of the last 15–20 years are slowly bringing to light evidence of the earliest phases of the Neolithic on the island. Unsurprisingly, however, this ‘marginality’ reasoning has discouraged archaeologists from actively searching for earlier, Pleistocene, sites on Cyprus. Using a geospatial predictive modeling approach, and keeping in mind the remarkable plasticity observed in hunter-gatherer societies (e.g. [70] and references therein), this paper proposes a systematic methodology for the study of the island’s Pleistocene exploitation. The geospatial model developed considered aspects of topography, such as elevation and slope, as well as environmental parameters, including fresh water and raw material sources. Although separately these parameters were not found to be of particular significance in site location, in combination they delineated areas, which were more attractive to Pleistocene exploitation. These areas are characterised by elevations <350 m and usually range between 0–100 m, slopes <20%, and mainly 0–10°, distance from the coast <500 m and freshwater sources 500–5000 m.

Proximity to water, especially freshwater, is important in the selection of site locations by hunter-gatherers on Cyprus and elsewhere. A recent study from Australia [59] demonstrated the importance of water in Palaeolithic societies by showing that 84% of archaeological sites >30 ka old are within 20 km of modern permanent water. In the Arabian subcontinent, approximately 45% of Lower and Middle Palaeolithic sites are located within 10 km of drainages and 68% and 58% respectively within 20 km of drainages [71]. A progressive drop in the number of recorded archaeological sites occurs as distance from major drainage features increases [72]. Cyprus today experiences a semi-arid/arid climate but several major permanent rivers transcend its territory. Numerous springs further contribute to fresh water accessibility. Although the island’s Pleistocene hydrological network is not well understood, the distribution of freshwater sources in the past is likely broadly comparable with their modern distribution. Permanent water points, such as rivers and springs, require particular lithologies, hydrogeological environments and geomorphic settings [56, 57]. Thus, while it is possible that smaller water sources (e.g. small springs, rockholes) may activate or deactivate over time, their general location as well as the location of larger water points can be safely assumed to be relatively constant. Apart from being vital to human survival, freshwater sources are important landscape features for human subsistence as they are rich with aquatic plants and prey species such as fish, shellfish and turtles and they also attract faunal species, thus, contributing to energy intake as locales for finding hunting prey and carcass scavenging [7375].

The coast would have also contributed to the subsistence supply for Pleistocene hunter-gatherers. Our model showed that distance to the coast functioned as an important parameter in site location preferences. This outcome adds to an existing body of evidence regarding the exploitation of the Cypriot coastline by Epipalaeolithic groups as attested in sites such as Akrotiri Aetokremnos [33, 76], Nissi Beach and Akamas Aspros [35] on southern Cyprus. The faunal assemblage from Akrotiri Aetokremnos, the earliest known Terminal Pleistocene site on Cyprus (12 ka), provides evidence for the consumption of shellfish and small quantities of fish, alongside pygmy hippos and elephants [33, 76].

The distribution, abundance and quality of raw material sources for stone tool manufacture conditions how these resources can be exploited [7780]. An examination of the geology of Cyprus shows an abundance of chert sources across the island [81, 82]; the Lapithos chalks with cherts border the northern Kyrenian range, while the larger Lefkara formation lies within the sedimentary succession that encircles the Troodos Mountains (see also [83] for a detailed description of the formations). Surveys in central and western Cyprus have identified over fifty chert and jasper sources, ranging in quality, located along the pillow lava and sedimentary boundary of the Troodos foothills [36, 37, 84]. Raw material procurement while undertaking subsistence activities minimises direct costs [85] and the abundance of sources in these localities would have facilitated an embedded strategy in these areas, allowing Pleistocene humans to carry out other activities in the knowledge that a chert source was in the vicinity whenever new tools were required. Our model shows that although immediate access to a good quality source of stone tool raw material was not a determining factor in site location, distance from such a resource influenced site preferences. All documented archaeological localities are found within geological zones with chert outcrops.

While field investigations have confirmed the power of our predictive analyses to successfully identify areas that preserve cultural material of a Pleistocene age, it has also highlighted zones of low probability to unearth prehistoric archaeology. High elevations, over 350–500 m, can be excluded from investigations targeting a pre-Holocene occupation of Cyprus. The availability of Pleistocene fauna that would have been available for hunting, namely pygmy hippos and elephants, at low elevations, supports this predictive outcome. We discussed elsewhere in this paper the plasticity of the human species to successfully adapt to resource availability at a range of environmental circumstances. By excluding certain areas from future investigations, our model does not deny Pleistocene hunter-gatherers on Cyprus the ability to occupy and explore the full extent of the island. Rather, it suggests localities that would have been preferentially selected for repeated visits and are, thus, more likely to preserve material evidence of prehistoric choices.

Movement is a primary means for humans to read and experience the landscape and connect with one another [24, 86, 87]. The manifestation of evidence of movement on the landscape largely depends on the frequency and duration of the routes used, possible intentionality of making them visible, and, their chance conservation through time [87]. The selection of a location for resource exploitation or settlement is a reflection of decisions that are influenced by the affordances of the natural environment [88, 89] but also by personal preference and social factors [90]. In the case of Pleistocene Cyprus, we are called to make sense of hunter-gatherer adaptation and human-insular environment relationship by looking for low density, dispersed lithic scatters that are notoriously difficult to detect in the field.

Conclusions

GIS-based analyses were performed on Pleistocene Cyprus as a heuristic model applied to determine likely hotspots of early human activity on the island. Topography alone fails to provide sufficient predictive markers; a variety of factors, primarily freshwater sources, lithic raw material availability and the distribution of floral and faunal resources, as well as social interactions, all played key and dynamically varying roles in these processes throughout the Pleistocene. The results of the analyses can be used predictively to identify areas worthy of investigation for their potential to yield previously unrecognised archaeological sites and delimit the pre-Holocene colonisation of the eastern Mediterranean island. Preliminary field surveys confirm the validity of the model and provide new data concerning early human presence on the island. We view the model we have presented here as a useful first step in attempting to approach the question of when and how early humans reached Cyprus and contribute to discussions of the potential characteristics of hominin lifeways during the Pleistocene (see also [7]). Further work and, especially, high-resolution palaeoenvironmental and palaeoclimatic records are required to refine these initial outcomes and excavations to ground truth the predictions. We are hopeful that our study will help set the foundations for targeted and multi-disciplinary Pleistocene research in Cyprus and, alongside other recent advances in archaeological thought, contribute to the island colonisation thematic in a part of the world, where this topic remains, peculiarly, largely unexplored.

Supporting information

S1 Table. Summary table of archaeological sites of a Pleistocene age documented in the Eastern Mediterranean.

https://doi.org/10.1371/journal.pone.0258370.s001

(XLSX)

Acknowledgments

Dr Theodora Moutsiou is the Principal Investigator of the project Water Routes in Human Island Dispersals: Modeling the Pleistocene Exploitation of Cyprus (PLEICY, 2018–2020), which is hosted at the Archaeological Research Unit of the University of Cyprus and supervised by Professor Vasiliki Kassianidou. The authors wish to dedicate this article to Dr Carole McCartney (†) who has tirelessly worked to bring the early prehistory of Cyprus to the forefront of archaeological investigations on the island and beyond.

References

  1. 1. Simmons AH. Stone Age sailors: Paleolithic seafaring in the Mediterranean. Walnut Creek, CA: Left Coast; 2014.
  2. 2. Leppard TP. The evolution of modern behaviour and its implications for maritime dispersal during the Palaeolithic. Camb Archaeol J. 2015; 25(4): 829–846.
  3. 3. van den Bergh GD, Li B, Brumm A, Grün R, Yurnaldi D, Moore MW et al. Earliest hominin occupation of Sulawesi, Indonesia. Nature 2016; 529 (7585): 208–211. pmid:26762458
  4. 4. Clarkson C, Jacobs Z, Marwick B, Fullagar R, Wallis L, Smith M. et al. 2017. Human occupation of northern Australia by 65,000 years ago. Nature 2017; 547: 306–310. pmid:28726833
  5. 5. O’Connor S, Louys J, Kealy S, Samper Carro SC. Hominin Dispersal and Settlement East of Huxley’s Line. Curr Anthropol. 2017; 58: S567–S582.
  6. 6. Louys J, Braje TJ, Chang C-H, Cosgrove R, Fitzpatrick SM, Fujita M. et al. No evidence for widespread island extinctions after Pleistocene hominin arrival. Proc Nat Acad Sci USA 2021; 118 (20): https://www.pnas.org/content/118/20/e2023005118 pmid:33941645
  7. 7. Moutsiou T. Climate, environment and cognition in the colonisation of the Eastern Mediterranean islands during the Pleistocene. Quat Int. 2021; 577: 1–14. https://doi.org/10.1016/j.quaint.2020.09.012
  8. 8. Henshilwood C, Sealy J, Yates R, Cruz-Uribe K, Goldberg P, Grine F, et al. Blombos Cave, Southern Cape, South Africa: Preliminary Report on the 1992–1999, Excavations of the Middle Stone Age Levels 2001; 28: 421–448.
  9. 9. Klein RG, Avery G, Cruz-Uribe K, Halkett D, Parkington J, Steele T. et al. The Ysterfontein 1 Middle Stone Age site, South Africa, and early human exploitation of coastal resources. Proc Natl Acad Sci USA 2004; 101: 5708–5715. pmid:15007171
  10. 10. Efstratiou N, McCartney C, Karkanas P, Kyriakou D. An upland early site in the Troodos Mountain. Report of the Department of Antiquities, Cyprus; 2010: 1–26.
  11. 11. Simmons AH, DiBenedetto K., Keach L. Neolithic Kritou Marottou-Ais Giorkis, Cyprus–Living in the Uplands. BASOR 2018; 379: 171–195.
  12. 12. Kohler TA, Parker SC. Predictive models for archaeological resource location. In: Schiffer MB, editor. Advances in archaeological method and theory Vol. 9. New York: Academic; 1986. pp. 397–452.
  13. 13. Kvamme KL. Development and testing of quantitative models. In: Judge WJ, Sebastian L, editors. Quantifying the present and predicting the past: Theory, method, and application of archaeological predictive modeling. Denver: U.S. Department of the Interior, Bureau of Land Management; 1988. pp. 325–428.
  14. 14. Kvamme KL. The fundamental principles and practice of predictive archaeological modeling. In: Voorrips A, editor. Mathematics and Information Science in Archaeology: A flexible framework. Bonn: Holos Verlag; 1990. pp. 257–295.
  15. 15. Verhagen P, Kamermans H, van Leusen M, Deeben J, Hallewas D, Zoetbrood P. First thoughts on the incorporation of cultural variables into predictive modeling. In: Verhagen JWHP. Case studies in archaeological predictive modelling. Doctoral Thesis, Archaeological Studies Leiden University 14. Leiden: Leiden University Press; 2007. pp. 203–210. Available from: https://scholarlypublications.universiteitleiden.nl/handle/1887/11863
  16. 16. Balla A, Pavlogeorgatos G, Tsiafakis D, Pavlidis G. Recent advances in archaeological predictive modeling for archeological research and cultural heritage management. Mediterr Archaeol Archaeom 2014; 14 (4): 143–153.
  17. 17. Verhagen JWHP. Case studies in archaeological predictive modelling. Doctoral Thesis, Archaeological Studies Leiden University 14. Leiden: Leiden University Press; 2007. Available from: https://scholarlypublications.universiteitleiden.nl/handle/1887/11863
  18. 18. Lock GS, Stancic Z, editors. Archaeology and Geographical Information Systems: A European perspective. London: Taylor and Francis; 1995.
  19. 19. Wheatley D, Gillings M. Spatial technology and archaeology: The archaeological applications of GIS. London: Taylor and Francis; 2002.
  20. 20. van Leusen M, Kamermans H, editors. Predictive Modelling for Archaeological Heritage Management: A research agenda. Amersfoort: Rijksdienst voor het Oudheidkundig Bodemonderzoek; 2005. https://doi.org/10.1039/b512135c pmid:16482351
  21. 21. Wansleeben M, Verhart LBM. Geographical Information Systems. Methodical progress and theoretical decline? Archaeological Dialogues 1997; 4: 53–70.
  22. 22. Verhagen P, Whitley TG. Integrating archaeological theory and predictive modeling: A live report from the scene. J Archaeol Method Theory 2012; 19: 49–100.
  23. 23. Altschul JH, Sebastian L, Heidelberg K. Predictive modeling in the military. Similar goals, divergent paths. Rio Rancho: SRI Foundation; 2004.
  24. 24. Ashmore W, Knapp BA. Archaeologies of landscape: Contemporary perspectives. USA: Blackwell Publishers Ltd; 1999.
  25. 25. Aldenderfer M, Maschner DG, editors. Anthropology, space and Geographic Information Systems. Oxford: Oxford University Press; 1999.
  26. 26. Phoca-Cosmetatou N, Rabett RJ. Pleistocene island occupation in the Mediterranean: Insights from a tied-biome approach to glacial refugia. In: Boyle K, Rabett RJ, Hunt CO, editors. Living in the landscape: Essays in honour of Graeme Barker. (McDonald Institute Monograph Series). Cambridge: McDonald Institute for Archaeological Research; 2014. pp. 83–108.
  27. 27. Strasser T, Panagopoulou E, Runnels C, Murray PM, Thompson NC, Karkanas P, et al. Stone Age seafaring in the Mediterranean: Evidence from the Plakias region for Lower Palaeolithic and Mesolithic habitation of Crete. Hesperia 2010; 79: 145–190.
  28. 28. Runnels C, Panagopoulou E, Murray P, Tsartsidou G, Allen S, Mullen K, et al. A Mesolithic Landscape in Greece: Testing a Site-Location Model in the Argolid at Kandia. J Mediterr Archaeol. 2005; 18 (2): 259–285.
  29. 29. Simmons AH. Faunal extinction in an island society. Pygmy hippopotamus hunters of Cyprus. Interdisciplinary Contributions to Archaeology. New York: Kluwer Academic/Plenum Publishers; 1999.
  30. 30. Bird MI, Condie SA, O’Connor S, O’Grady D, Reepmeyer C, Ulm S, et al. Early human settlement of Sahul was not an accident. Sci Rep. 2019; 9: 8220. pmid:31209234
  31. 31. Kopaka K, Matzanas C. Paleolithic industries from the island of Gavdos, near neighbour to Crete in Greece. Antiquity Project Gallery 2009; 321.
  32. 32. Runnels C, DiGregorio C, Wegmann KW, Gallen SF, Strasser T, Panagopoulou E. Lower Palaeolithic artifacts from Plakias Crete: implications for hominin dispersals. In: Ammerman AJ, Davis T, editors. Island Archaeology and the origins of seafaring in the Eastern Mediterranean, Eurasian Prehistory 2013; 11 (1–2): 129–152.
  33. 33. Simmons AH. Preliminary report of the interdisciplinary excavations of Akrotiri-Aetokremnos (site E): 1987, 1988, 1990. Report of the Department of Antiquities, Cyprus 1991: 7–14.
  34. 34. Simmons AH. Of tiny Hippos, Large Cows and Early Colonists in Cyprus, J. Mediterr Archaeol. 1998; 11 (2): 232–241.
  35. 35. Ammerman AJ, Flourentzos P, McCartney C, Noller J, Sorabji D. Two new early sites on Cyprus. Report of the Department of Antiquities, Cyprus 2006: 1–22.
  36. 36. McCartney C, Croft P, Manning SW, Rosendahl S. Preliminary report on the 2008 EENC excavations at Agia Varvara Asprokremmos and regional field survey. Report of the Department of Antiquities, Cyprus 2009: 1–16.
  37. 37. McCartney C, Manning SW, Stewart ST. Agia Varvara-Asprokremmos 2009: Excavations, radio-carbon dating and geo-chemical analysis of chert sources. Report of the Department of Antiquities, Cyprus, 2010: 77–94.
  38. 38. Vigne J-D, Briois F, Guilaine J. Klimonas, the oldest Pre-Pottery Neolithic village in Cyprus. In: Astruc L, McCartney C, Briois F, Kassianidou V, editors. Near Eastern lithic technologies on the move. Interactions and contexts in Neolithic traditions. 8th international conference on PPN chipped and ground stone industries of the Near East, Nicosia, November 23rd-27th 2016. Studies in Mediterranean Archaeology Vol. CL. Astrom Editions; 2019. pp. 3–11.
  39. 39. Vita-Finzi C. Palaeolithic finds from Cyprus? Proc Prehist Soc. 1963; 39: 453–454.
  40. 40. Adovasio JM, Fry GF, Gunn JD, Maslowski RF. Prehistoric and historic settlement patterns of western Cyprus with a discussion of Cypriote Neolithic stone tool technology. World Archaeol. 1975; 6 (3): 339–364.
  41. 41. Baudou E, Engelmark R. The Tremithos Valley Project. A preliminary report for 1981–1982. Report of the Department of Antiquities, Cyprus, 1983: 1–8.
  42. 42. Strasser TF, Runnels C, Vita-Finzi C. A possible Palaeolithic hand axe from Cyprus. Antiquity 2016; 90: 350.
  43. 43. Simmons AH. Preliminary report on multidisciplinary investigations at Neolithic Kholetria-Ortos, Paphos District. Report of the Department of Antiquities, Cyprus, 1996: 29–44.
  44. 44. Verhagen P. Prospection strategies and archaeological predictive modelling. In: van Leusen M, Kamermans H, editors. Predictive modelling for archaeological heritage management: A research agenda. (NAR 29). Amersfoort: Rijksdienst voor het Oudheidkundig Bodemonderzoek; 2005. pp. 109–122.
  45. 45. Bar-Yosef O, Belfer-Cohen A. From Africa to Eurasia: Early dispersals. Quat Int. 2001; 75: 19–28.
  46. 46. Galanidou N. Advances in the Palaeolithic and Mesolithic archaeology of Greece for the new millennium. Pharos 2014; 20 (1): 1–40.
  47. 47. Kadowaki S, Omori T, Nishiaki Y. Variability in Early Ahmarian lithic technology and its implications for the model of a Levantine origin of the Protoaurignacian. J Hum Evol. 2015; 82: 67–87. pmid:25924809
  48. 48. Dinçer B. The Lower Paleolithic in Turkey: Anatolia and hominin dispersals out of Africa. In: Harvati K, Roksandic M, editors. Paleoanthropology of the Balkans and Anatolia: Human evolution and its context. Vertebrate Paleobiology and Paleoanthropology Series, Springer; 2016. pp. 213–228.
  49. 49. King G, Bailey G. Tectonics and human evolution. Antiquity 2006; 80 (308): 265–286.
  50. 50. Bailey GN, King GCP. Dynamic landscapes and human dispersal patterns: tectonics, coastlines, and the reconstruction of human habitats. Quat Sci Rev. 2011; 30 (11–12): 1533–1553.
  51. 51. Sakellariou D, Galanidou N. Pleistocene submerged landscapes and Palaeolithic archaeology in the tectonically active Aegean region. In: Harff J, Bailey G, Lüth F, editors. Geology and Archaeology: Submerged landscapes of the continental shelf. Geological Society of London, Special Publications 411; 2016. pp. 145–178.
  52. 52. Bailey G, Galanidou N, Peeters H, Jöns H, Mennenga M. The Archaeology of Europe’s drowned landscapes. Coastal Research Library 35. SpringerOpen; 2020.
  53. 53. Barton CM, Villaverde V, Zilhão J, Aura JE, Garcia O, Badal E. In glacial environments beyond glacial terrains: Human eco-dynamics in late Pleistocene Mediterranean Iberia. Quat Int. 2013; 318: 53–68. https://doi.org/10.1016/j.quaint.2013.05.007.
  54. 54. Barton CM, Aura Tortosa JE, Garcia-Puchol O, Riel-Salvatore JG, Gauthier N, Vadillo Conesa M, et al. Risk and resilience in the late glacial: A case study from the western Mediterranean. Quat Sci Rev. 2018; 184: 68–84. https://doi.org/10.1016/j.quascirev.2017.09.015
  55. 55. Peltenburg E. East Mediterranean water wells of the 9th–7th millennium BC. In: Klimscha F, Eichmann R, Schuler C, Fahlbusch H, editors. Wasserwirtschaftliche innovationen im archäologischen kontext von den prähistorischen anfängen bis zu den metropolen der antike, Rahden/Westf.: Marie Leidorf GmbH; 2012. pp. 69–82.
  56. 56. Morton SR, Short J, Barker RD. Refugia for biological diversity in arid and semiarid Australia. Department of the Environment Water Heritage and the Arts Biodiversity Series Paper 4; 1995.
  57. 57. Wischusen JD, Fifield LK, Cresswell RG. Hydrogeology of Palm Valley, central Australia; a Pleistocene flora refuge? J Hydrol. (Amst) 2004; 293 (1–4): 20–46.
  58. 58. Davis J, Pavlova A, Thompson R, Sunnucks P. Evolutionary refugia and ecological refuges: Key concepts for conserving Australian arid zone freshwater biodiversity under climate change. Glob Change Biol. 2013; 19 (7): 1970–1984. pmid:23526791
  59. 59. Bird MI, O’Grady D, Ulm S. Humans, water, and the colonization of Australia. Proc Natl Acad Sci USA 2016; 113 (41): 11477–11482. pmid:27671630
  60. 60. Tullström NHO. Hydrogeological map of Cyprus. Geological Survey Department of Cyprus; 1970.
  61. 61. Enzel Y, Ely LL, Mishra S. High-resolution Holocene environmental changes in the Thar Desert, Northwestern India. Science 1999; 284: 125–128. pmid:10102808
  62. 62. Kar A, Singhvi AK, Rajaguru SN, Juyal N, Thomas JV, Banerjee D, et al. Reconstruction of late Quaternary environment of the lower Luni Plains, Thar Desert, India. J Quat Sci. 2001; 16: 61–68.
  63. 63. Srivastava P, Singh IB, Sharma M, Singhvi AK. Luminescence chronology and Late Quaternary geomorphic history of the Ganga plain, India. Palaeogeogr Palaeoclimatol Palaeoecol. 2003; 197: 15–41.
  64. 64. Bhandari S, Maurya DM, Chamyal LS. Late Pleistocene alluvial plain sedimentation in Lower Narmada Valley, Western India: Palaeoenvironmental implications. J Asian Earth Sci. 2005; 24: 433–444.
  65. 65. Andrefsky W. Lithics: Macroscopic approaches to analysis. Cambridge: Cambridge University Press; 2000. https://doi.org/10.1097/00003246-200005000-00008 pmid:10834670
  66. 66. Mithen SJ. Thoughtful foragers: A study of prehistoric decision making. Cambridge: Cambridge University Press; 1990.
  67. 67. Roebroeks W, Conard NJ, van Kolfschoten T. Dense Forests, Cold Steppes, and the Palaeolithic Settlement of Northern Europe. Curr Anthropol. 1992; 33 (5): 551–586.
  68. 68. Roebroeks W, Kolen J, Rensink E. Planning depth, anticipation and the organization of Middle Palaeolithic technology: The archaic natives meet Eve’s descendants. Helinium 1988; 28 (1): 17–34.
  69. 69. Kuhn SL. On planning and curated technologies in the Middle Paleolithic. J Anthropol Res. 1992; 48 (3): 185–214.
  70. 70. Roberts P, Stewart BA. Defining the ‘generalist specialist’ niche for Pleistocene Homo sapiens. Nat Hum Behav. 2018; 2(8): 542–550. pmid:31209320
  71. 71. Breeze PS, Drake NA, Groucutt HS, Parton A, Jennings RP, White TS, et al. Remote sensing and GIS techniques for reconstructing Arabian palaeohydrology and identifying archaeological sites. Quat Int. 2015; 382: 98–119. https://doi.org/10.1016/j.quaint.2015.01.022
  72. 72. Field JS, Petraglia MD, Lahr MM. The Southern dispersal hypothesis and the South Asian archaeological record: Examination of dispersal routes through GIS analysis. J Anthropol Archeol. 2007; 26: 88–108.
  73. 73. Thrash I, Theron GK, Bothma JP. Dry season herbivore densities around drinking troughs in the Kruger National Park. J Arid Environ. 1995; 29: 213–219.
  74. 74. Redfern JV, Grant R, Biggs H, Getz WM. Surface-water constraints on herbivore foraging in the Kruger National Park, South Africa. Ecology 2003; 84: 2092–2107.
  75. 75. Valeix M, Loveridge AJ, Davidson Z, Madzikanda H, Fritz H, Macdonald DW. How key habitat features influence large terrestrial carnivore movements: Waterholes and African lions in a semi-arid savanna of north-western Zimbabwe. Landscape Ecol. 2010; 25: 337–351.
  76. 76. Simmons AH. Akrotiri-Aetokremnos (Cyprus) 20 years later: An assessment of its significance. In: Ammerman AJ, Davis TW, editors. Island Archaeology and the origins of seafaring in the Eastern Mediterranean. Eurasian Prehistory 2013; 10 (1–2): 139–155.
  77. 77. Andrefsky W. Raw-Material Availability and the Organization of Technology. Am Antiq. 1994; 59 (1): 21–34.
  78. 78. Bamforth DB. Technological efficiency and tool curation. Am Antiq. 1986; 51: 38–50.
  79. 79. Bamforth DB. Settlement, Raw Material, and Lithic Procurement in the Central Mojave Desert. J Anthropol Archaeol. 1990; 9: 70–104.
  80. 80. Bamforth DB. Quarries in Context: A Regional Perspective on Lithic Procurement. In: Arnold JE, editor. Stone tool procurement and distribution in California prehistory. Perspectives in California Archaeology. Vol. 2. Los Angeles: Institute of Archaeology, University of California; 1992. pp. 131–150.
  81. 81. Constantinou G, Panagides I. Cyprus and geology: Science, environment, culture. Nicosia: Bank of Cyprus Cultural Foundation; 2013.
  82. 82. Geological Survey Department of Cyprus. Geological map of Cyprus; 1995.
  83. 83. McCartney C. Prehistoric chert use in Cyprus. In: Delage C, editor. Chert availability and prehistoric exploitation in the Near East. Oxford: British Archaeological Reports International Series 1615; 2007. pp. 304–326.
  84. 84. Stewart ST. Hill and vale: Understanding prehistoric lithic use in Cyprus. In: Peltenburg E, Wasse A, editors. Neolithic revolution. New perspectives on Southwest Asia in light of recent discoveries on Cyprus. Oxford: Oxbow Books (Levant Supplementary Series, Volume 1); 2004. pp. 123–132.
  85. 85. Duke C. Quantifying Palaeolithic landscapes: Computer approaches to terrain analysis and visualization. In: Doerr M, Sarris A, editors. The digital heritage of Archaeology. CAA2002. Computer Applications and Quantitative Methods in Archaeology. Proceedings of the 30th CAA Conference, Heraklion, Crete, April 2002; 2003. pp.138-145.
  86. 86. Balée W. The research program of historical ecology. Annu Rev Anthropol. 2006; 35: 75–98.
  87. 87. Milheira RG, De Souza JG, Iriarte J. Water, movement and landscape ordering: A GIS-based analysis for understanding the mobility system of late Holocene mound-builders in southern Brazil. J Archaeol Sci. 2019; 111: 105014.
  88. 88. Llobera M. Exploring the topography of mind: GIS, social space and archaeology. Antiquity 1996; 70 (269): 612–622.
  89. 89. Gillings M. Landscape phenomenology, GIS and the role of affordance. J Archaeol Method Theory 2012; 19: 601–611.
  90. 90. Podobnikar T, Veljanovski T, Stanèiè Z, Oštir K. Archaeological predictive modeling in Cultural Resource Management. In: Konečný M, editor. GI in EUROPE: Integrative—Interoperable–Interactive. Proceedings of 4th Agile Conference on Geographic Information Science, 2001, Brno, Czech Republic; 2001. pp. 535–544.