Species Distribution Modeling (SDM) is used to derive spatially explicit predictions of environmental suitability for a species (plant or animal). This is typically achieved by classifying grid cells (also called “background locations”) according to the degree in which they are suitable for a species in the study region (Elith & Leathwick 2009; Guisan & Thuiller 2005; Miller 2010; Phillips, Anderson & Schapire 2006). Among the available tools for SDM, Maximum Entropy Modeling (Maxent) is one of the most widely used approaches in ecology for predicting habitat suitability from presence-only records (Phillips, Dudík & Schapire 2021). This method involves using a collection of known species locations as sample data, along with relevant environmental factors to model the distribution of that species within a known geographical extent (Phillips, Dudík & Schapire 2021). Having the ability to derive predictive models from presence-only records is particularly relevant in archaeology, as archaeologists often know where some, but not all, sites are located, and they rarely know for certain where past sites are absent (Howey, Palace & McMichael 2016). Maxent has been broadly used by archaeologists to predict archaeological site locations in different environments and time-periods around the world (Conolly et al. 2012; Franklin et al. 2015; Galletti et al. 2013; Howey, Palace & McMichael 2016; Howey et al. 2020; McMichael et al. 2014; Muttaquin, Heru & Susilo 2019; Noviello et al. 2018; Politis et al. 2011, 2019). When compared to traditional predictive models, such as GIS-based multiparametric spatial analysis (i.e., logistic regression), Maxent has been shown to perform well for archaeology predictive models (Noviello et al. 2018; Wachtel et al. 2018). The Maxent modeling provides an initial platform that is flexible and can be run anytime there is additional input data, such as new archaeology site locations or refined environmental layers.
Few archaeological predictive models have been built for the Pampas region of Argentina. One of the reasons for this is related to the properties of the current landscape, which is largely spatially homogenous, temporally dynamic with variations on a year-to-year basis due to extreme weather conditions (e.g., flooding), and highly affected by modern farming activities. Apart from the hills (sierras), the environmental characteristics of the Pampas grasslands are “less than optimal” for predictive models, which work best in heterogeneous landscapes that are relativity stable over long periods. Although the study area is not ideal, refined high resolution and freely available geo-environmental resource data are available. These, combined with powerful geostatistical software are providing the opportunity to build predictive models in less-than-optimal environments, such as forests with high tree coverage, and wetland areas (Jones 2017). At the local and sub-regional scale, archaeological predictive modeling has been used along the north coast of the Pampas region (Eugenio, Aldazabal & Macchi 2013; Macchi 2010), and in different sectors of the Tandilia hills (Mariano et al. 2014; Mazzia & Gómez, 2013). At the regional scale, a Maxent model was built to infer paleoenvironmental conditions during pre-historic times using multiple environmental variables and the presence data of three native fauna species (i.e., Lama guanicoe, Ozotoceros bezoarticus and Blastocerus dichotomus) (Politis et al. 2011). For the southern Pampas and northeast Patagonia, prediction maps were built to study the distribution of lithic artifacts using geostatistical interpolation tools to discuss issues related with the archaeological lithic landscapes (Barrientos, Catella & Oliva 2015). Finally, just north of the Pampas region, in the Northeast of Argentina, a Maxent predictive model was built for archaeological sites of the Goya-Malabrigo archaeological entity (Politis et al. 2019). This last study, which uses geographically referenced sites for the occurrence data (i.e., species), is the closest in terms of methodology to the present one. All of these studies demonstrated the utility of geospatial software to quantify large archaeological databases containing multiple variables, and to integrate different sources of information, allowing a better understanding of pre-historic hunter-gatherer groups from the Pampas region.
The southern Pampas offered abundant and diverse resources for hunter-gatherer groups, including –but not limited to- large sized game in the inter-hill Pampas (or grasslands), shelter and good quality stones for making tools in the hills (Tandilia and Ventania ranges), and marine faunal resources and secondary lithic raw materials along the Atlantic coast. Evidence of consumption of marine resources has been found in sites at least 50 km from the coast (Politis, Scabuzzo & Tykot 2009); however, groups from the Pampas region mainly consumed terrestrial herbivores, such as the guanaco (Lama guanicoe) and the Pampas deer (Ozotoceros bezoarticus) (Gutiérrez & Martínez, 2008; Martínez et al. 2016). Archaeological data suggests that hunter-gatherers occupied diverse environments within the annual cycle of mobility, and maintained social networks over large geographic areas, which facilitated access to distant resources and information (Armentano, Martínez & Gutiérrez 2007; Barros 2012; Barros, Messineo & Colantonio 2015; Bayón & Flegenheimer 2004; Martínez et al. 2015; Mazzanti 2006; Politis 2008). For the later part of the Holocene, in relation to a growing increase in population density, data suggests a reduction of residential mobility; though, a hunter-gatherer way of life persisted until the Spanish conquest.
In the immense extent of the southern Pampas, where the ground surface is covered by vegetation and soil formation predominates, with thousands of shallow lakes, multiple levels of streams, and variable changes in elevations; archaeological surveying is challenging. Traditional survey strategies are usually driven by informants, such as landowners or farmers who come across isolated remains or sites. The riverbanks, streams, and shallow lakes are particularly useful for exposing in-situ archaeological material. During wet periods, the riverbanks and streams erosion expose older deposits, giving archaeologists a window into the buried record, and guiding most archaeological pedestrian survey when no prior information is provided. Along the hills, survey strategies are similar and mainly oriented to the exploration of caves or rock shelters where there is a better preservation of sedimentary deposits. Although animal bone does not preserve well in these sites; most of them present exceptional records of stone tool technology (Flegenheimer, Mazzia & Weitzel 2015; Martínez 2017).
The following paper presents the results of a SDM for hunter-gatherer archaeological sites in the southern Pampas, using the geographical location of archaeological sites as the occurrence data (species). The objective is to build a probability distribution of pre-historic hunter-gatherer sites that helps advance pedestrian surveys towards sectors of the landscape with greater archaeological potential (favorable characteristics for human settlement in the past). The SDM will also help characterize and discuss site locations and regional distribution patterns of hunter-gatherer occupations; and may also serve as a potential model for future heritage management projects. Research on mobility and land use is of great importance for understanding the way of life of the hunter-gatherer societies that inhabited the Pampas region. If human behavior is modeled with respect to the natural and social environment, we can learn a lot about how people interacted and moved in space by observing the relationships between the archaeological record and the environment.
One of the major challenges of SDM studies is to appropriately determine the extent of the study area. The study area should be objectively determined to cover accessible areas by the species within its known range, allowing for an assortment of environmental variation and extremes occupied by the species (Barve et al. 2011; Kantner 2008; Sánchez-Fernández, Lobo & Hernández-Manrique 2011). It should also be constrained enough not to introduce a spatial sampling bias. The SDM presented here is designed to capture hunter-gatherer sites from the southern Pampas region of Argentina; specifically, it will be centered on the inter-hill Pampas eco-complex and will capture neighboring areas such as the Sierra Tandilia and Sierra Ventania Hills, as well as parts of the Flooding Pampas, Chained Lakes, and Sandy Pampas; for a total area of approximately 124,740 km2 (not including the ocean) (Figure 1).
The executable Maxent software (v.3.4.4) (Phillips, Dudík & Schapire 2021) has several modifiable parameters that influence the model performance, and ultimately the effectiveness of the predictive model (e.g., Barbet-Massin et al. 2012; Elith, Kearney & Phillips 2010; Fourcade et al. 2014; Merow, Smith & Silander 2013; Morales, Fernández & Baca-González 2017; Phillips et al. 2017). There are several discussions as to how to objectively set these parameters; ultimately, we need to have a clear understanding of the specific dataset (species) and knowledge of the spatial distribution of sampling effort (Araújo & Guisan 2006; Pearce & Boyce 2006). Here, I adjusted three parameters to optimize the predictive accuracy: (1) bias correction, (2) feature combination, and (3) the regularization multiplier. The remaining parameters (see Table 1 and Data Availability section) were left with the default software values (Phillips & Dudík 2008). For example, for the replicate type, “cross-validation” was used, which has been shown to work better with small data sets and helps with spatial autocorrelations and overfitting (Merckx et al. 2011; Phillips 2017). In order to use this re-sampling procedure, the value of “random test percentage” must be set to zero.
|Feature selection:||Linear, Quadratic, Product, Hinge.|
|Variable assessment:||Create response curves, Do jackknife.|
|Random test percentage:||0|
|Regulation multiplier:||Tested at 0.25, 0.5, 0.75, 1.0, 1.25|
|Max number of background points:||10,000|
|Replicates:||1 (30 with optimal model)|
|Replicated run type:||Crossvalidate|
When using presence-only occurrence data in SDMs, sampling bias is a serious issue (see discussion in Phillips et al. 2009), especially in archaeology. The only way to remove sampling bias is through study design and/or modelling of that bias (Yackulic et al. 2013: 238). Few archaeological surveys are made using planned random survey coverage strategies (Banning 2021; Orton 2000). Without clear records on survey effort across the landscape, we cannot distinguish between areas that are environmentally unsuitable and those that are under-sampled (Elith, Kearney & Phillips 2010: Appendix S1). Constraints and opportunities of survey design result in many areas that are overrepresented, or alternatively, underrepresented because of accessibility (e.g., distance to towns, roads, land ownership, etc.) or research bias (projects dedicated to studying specific areas of the landscape, well known hotspots, etc.). To help correct sampling bias, Maxent has three main parameter settings. First, the program discards redundant records that occur in a single grid cell; second, 10,000 pseudo-absence sample background points are generated (see Barbet-Massin et al. 2012); and third, an optional bias grid file can be added to the model. Using the location of known archaeological sites, I decided to build a bias grid file in QGIS (QGIS Development Team, 2021), following the methods proposed by Elith, Kearney & Phillips (2010). The map was derived using kernel density estimation, and set to a maximum radius of 10 km around each site. The choice of that radius is based on information of hunter-gatherers daily foraging trips (see: Kelly 1995; Lee 1979; Morgan 2008). This is admittedly variable for different regions and groups, as logistical trips greater than 10 km were likely to occur in some circumstances; however, this method provides a reasonable estimate for the purpose of this exploratory model. The result is a continuous surface layer with values ranging between 1 thru 12.
To avoid overfitting (i.e., matching the input data too closely), Maxent allows the user to adjust the feature classes (also called types) used to build the model (linear, quadratic, product, threshold, and hinge features). The default settings allow certain feature classes to be used based on the sample size of occurrence points. Allowing more feature classes enables a more flexible fit to the observed data; however, higher flexibility has been shown to result is an overfitted model (Hastie, Tibshirani & Friedman 2009; Muscarella et al. 2014; Peterson et al. 2011; Phillips & Dudík 2008). The regularization multiplier (β) also affects how focused the output distribution is. A small regularization multiplier will result in a more localized output (Phillips et al., 2006; Phillips & Dudík 2008). Essentially, adjusting both the regularization multiplier and feature classes will provide the best predictive performance. The only way to know which to use requires multiple runs. Here, I ran all possible combinations of feature classes (n = 15), excluding threshold (which is not recommended for smaller samples), against five regularization multiplier increments (0.25, 0.5, 0.75, 1.0 and 1.25). This led to 75 unique combinations. The optimal parameters were selected based on the lowest Akaike Information Criterion (AICc) with a correction for small sample size value (Muscarella et al. 2014; Warren & Seifert 2011), which I calculated using the executable software ENMTools (Warren, Glor & Turelli 2010; Warren & Seifert 2011). The optimal parameters correspond to the combination of feature classes and a regularization multiplier with the smallest AICc value. A low AICc value estimates the relative importance of variables, the suitability of habitat, and the transferability of the model to different periods (Warren & Seifert, 2011); the latter of which is particularly relevant in archaeology. The AICc value makes a useful criterion to objectively choose appropriate settings for Maxent, and it has also been shown to help reduce overfitting issues and complexity (Li et al. 2020); however, it is not fully related with the predictive accuracy of the SDM (see discussion in Velasco & González-Salazar 2019).
Once the optimal parameters were determined, I made a final run of the model in Maxent with response curves and a jackknife variable response to display and determine how each environmental variable contributed to the model (marginal response) and performed on its own (corresponding response). I ran 30 replicates with cross-validation to bring the average of the outputs to a suitable final model (Morales, Fernández & Baca-González 2017; Spiers et al. 2018). Each replicate randomly removed between 2 to 3 presence locations for testing. To help characterize the model performance, I used the Area Under the receiver operator characteristic Curve (AUC) evaluation metric (Warren & Seifert 2011). The AUC is the probability that a randomly chosen presence site will be ranked above a randomly chosen absence site (Merow, Smith & Silander 2013; Phillips and Dudík 2008). The values range from 0 to 1. A random ranking has on average an AUC of 0.5, and a perfect ranking achieves the best possible AUC of 1.0. Models with values above 0.75 are considered potentially useful (Elith 2000).
Georeferenced archaeological sites (n = 66) were used as the presence-only occurrence data (Figure 2 and Data Availability section). The samples used in this model included all pre-contact hunter-gatherer sites dated between 12,240 and 470 14C years BP. Only sites with radiocarbon dates were used in this model. While site density clearly fluctuates with time, given the low sample size for some time periods, and the mixed nature of much of the regional archaeological record; it was appropriate to begin to understand the patterning from a location-based perspective rather than a period specific approach (Conolly 2018: 197). In four cases, there were no published coordinates, just location descriptions or maps, which I georeferenced in Google Earth (Google Earth Inc. 2020). It was imperative to maintain these sites, not only because they are significant archaeology findings, but also because it is important to maximize the sample size for the SDM (van Proosdij et al. 2016; Wisz et al. 2008). From the total number of sites (n = 76), 10 sites were removed because they share a single cell with a neighboring site, which is part of the bias correction method as described above. Deciding which neighboring site to remove was a random selection.
Determining which environmental layers to use in any SDM is dependent on spatial scale, time-period, known environmental changes (datasets), and the species under investigation. Bioclimatic variables, such as seasonality (e.g., annual range in temperature and precipitation) and extreme or limiting environmental factors (e.g., temperature of the coldest and warmest month) are key to determine the distribution of many plants and animals (Fick & Hijmans 2017); however, humans are less susceptible to these variables, particularly at a smaller scale of analysis such as the southern Pampas, where humans where likely able to adapt to ongoing changes in their environment (Bettinger 1981; Kelly 1983; Yellen 1977). In the Pampas region, there is evidence of climate shifts during the last ca. 12,000 14C years BP (known date for the initial arrival of humans). The most dramatic changes occurred during the Late Pleistocene, when dry and cold weather and lower precipitation levels (ca. 100 mm) resulted in a sea-level drop of ca. 60 m below present day (Cione, Tonni, & Soibelzon 2003; Ponce et al. 2011; Quattrocchio et al. 2008; Zárate & Blasi 1991). During the Early to Middle Holocene, a warm period (called Hypsithermal), which lasted from ca. 8000 to 4500 14C years BP., brought with it more humid conditions accompanied by marine ingression and flooding (Aguirre & Whatley, 1995; Bonadonna, Leone, & Zanchetta 1999; Cruz, Prado & Arroyo-Cabrales 2021; Grill et al. 2007; Iriondo, Brunetto & Kröhling 2009; Mancini et al. 2005; Prado & Alberdi 1999; Prieto 2000; Tonni, Cione & Figini 1999; Zárate, Espinosa & Ferrero 1998). After this, the climate eventually stabilized.
High-resolution paleo-environmental data in the study area is fragmentary, and in some cases controversial (see references above). Most datasets are localized models, limited to specific sites or localities. Ideally, all environmental layers should be based on models of the past climate (e.g., paleo-SDM, see Franklin et al. 2015); however, these datasets are often geographically sparse and may not intersect with the time of interest. Projecting this data onto a larger scale can untimely have negative consequences on the effectiveness of the predictive model. Due to the dynamic formation processes in the landscape (both past and present), chronological resolution is also typically poor, as the archaeological record is frequently time-averaged palimpsests of multiple occupations.
Taking into consideration the above information, I decided to focus on some physical attributes of the landscape that remained stable over extended periods of time (Table 2, Figure 3). These include elevation, slope, aspect, permanent water sources, topographic wetness index (TWI) (an index in hydrological analysis that describes the tendency of an area to accumulate water) (Mattivi et al. 2019), and toolstone sources (areas where hunter-gatherers went to collect raw material for making tools). Elevation, slope, aspect, and the TWI were processed using a 3 arc second (~90-meter) void filled digital elevation model (DEM) in QGIS. While higher resolution DEMs are available for the region (e.g., 1 arc), the 90-meter DEM model is considered the most appropriate here primarily because finer scale DEMs tend to capture non-essential features (noise), such as buildings or roads which can skew the results.
|ENVIRONMENTAL LAYER||SOURCE (TYPE)||GIS PROCESSING||RASTER VALUE|
|Digital Elevation Model (DEM)||STRM 3 arc second void-filled (raster).||Raw data. Clipped to extent and pixel size, and converted to .asc.||Meters above sea level.|
|Slope||Derived from DEM (raster).||r.slope.aspect||Degrees|
|Aspect||Derived from DEM (raster).||r.slope.aspect||Degrees|
|Topographic Wetness Index (TWI)||Derived from DEM (raster).||r.watershed||TWI|
|Watercourse||IGN (vector).||Euclidean Distance (r.grow.distance).||Normalized matrix (0 to 1).|
|Bodies of water||IGN (vector).||Euclidean Distance (r.grow.distance).||Normalized matrix (0 to 1).|
|Toolstone source||SEGMAR (vector).||Euclidean Distance (r.grow.distance).||Normalized matrix (0 to 1).|
The permanent water sources (watercourses and bodies of water, not including the ocean) were derived from vectorized layers provided by the Instituto Geográfico Nacional (IGN), Argentina. The watercourses include permanent rivers, arms, streams, ravines, and estuaries. The bodies of water include both lakes and shallow lakes. Essentially, the size, shape and formation of these water sources changed throughout the Holocene (Iriondo & Garcia 1993; Johnson et al. 2012; Stutz et al. 2010; Stutz & Prieto, 2003; Zárate et al. 2000). For example, along the Quequén Grande River, evidence suggests that in some areas during the Early and Middle Holocene, a channel simply did not exist; just a broad, marshy, low-energy environment (e.g., paleo-swamps) represented by interconnected longitudinal pools of flowing water (Martínez & Gutiérrez 2018; Zárate et al. 2000). Apart from some dry pulses (see discussion in Martínez & Gutiérrez 2018), paleo-environmental data indicates that throughout the Holocene, water was available in these locations (in one form or another).
Toolstone sources are highly localized and heterogeneously distributed across the southern Pampas (Bayón et al. 1999; Messineo & Barros 2015). Lithic resources were available both in the Tandilia and Ventania Hills, along the Atlantic coast, and in some smaller isolated outcrops in the inter-hill Pampas. Diverse rock types were suitable for tool manufacture (e.g., silicified dolomite, chert, quartz, rhyolite, granite, sandstone, and costal cobbles), but it was the fine-grained orthoquartzite from the Sierras Bayas Group Formation, located primarily in the central sector of the Tandilia Hills, that was the regionally predominant raw material for stone tools throughout the Holocene (Bayón, Flegenheimer & Pupio 2006; Colombo 2011; Flegenheimer, Mazzia & Weitzel 2015; Flegenheimer et al. 1996; Politis, Messineo & Kaufmann 2004). For the predictive model, I decided to include the distance to this fine-grained orthoquartzite, as the primary toolstone source for hunter-gatherers in the area. To determine the limits of the Sierras Bayas Group Formation, I georeferenced and vectorized in QGIS, multiple map sources of the Tandilia Hill geology (Cingolani 2011; Rapela et al. 2007, 2011), including a web-based layer provided by the Instituto de Geologia y Recursos Minerales (SEGMAR). I also used the 90-meter DEM to help define the spatial limits of the formation.
Environmental layers with collinearity will affect the interpretation of the results. To correct this parameter, some authors suggest eliminating highly correlated variables (see Perkins-Taylor & Frey 2020). Multicollinearity between predictors was investigated using ENMTools Correlation Pearson correlation coefficient (Warren, Glor & Turelli 2010).
The Maxent predictive model was evaluated using the average test AUC output, the estimates of relative contributions of the environmental variables, and the marginal and corresponding response curves. The average test AUC output is calculated over the replicate runs (30), and provides the most comprehensive descriptive analysis for interpreting the results, including the standard deviation, which is used to evaluate the model performance. All of these datasets are reported in the Maxent output files.
To help further evaluate the Maxent model performance, I used information from previously completed pedestrian transects in a section of the inter-hill Pampas (Figure 4). Transects were purposely directed towards environmental areas with watercourses, bodies of water, and ease of access (Rafuse & Massigoge 2018). Currently, the total lineal distance covered by the transects is approximately 81 km. This work has resulted in the discovery of various isolated artifacts, and artifact clusters or surface sites (i.e., assemblages of at least five remains in approximately 10 m2). To avoid bias, any transects which were in the same grid cell as previously known sites (occurrence data used in the Maxent model) where not included in this evaluation. While the Maxent prediction model is not specifically designed to locate small scale isolated materials or clusters of archaeological materials, this field evidence is what generally leads to sub-surface exploration (e.g., shovel tests) and the eventual finding of buried archaeological sites (i.e., occurrence data).
Archaeological sites from the study area show a dispersed trend (Nearest neighbor index = 72146, z score = 1121277). The average distance between sites is 163 km, with a minimum distance of 21 m and a maximum distance of 465 km. Many of the sites are located within a range of 89–223 km from each other. The density of sites is 0.0005 sites per km2. The highest density per km2 is in the inter-hill Pampas (n = 25; 38%; 0.0009 per km2); while the greater number (n) of sites is in the hills (n = 35; 53%; 0.0007 per km2). However, if we considered Tandilia and Ventania as two separate areas, the highest density per km2 is recorded in the Tandilia range (n = 25; 0.0019 per km2). This range shows some possible clusters in the southern sector (as shown in the Bias grid map, Figure 3). These clusters have been interpreted by some researchers as a result of hunter-gatherer decisions based on geographic features and resource availability (Mazzanti & Bonnat 2013; Politis et al. 2004); while others suggest geological bias created by a greater preservation of archaeological materials in caves and rock shelters (Martínez et al. 2015). Research bias may also be a factor, as this sector of the hill has been under systematic investigation since the 1990s (Mazzanti 1993, 1997). The low-lying Pampas is the eco-complex with the lowest density of sites (9%; 0.0002 per km2). While not an independent eco-complex, it is important to mention that along the Atlantic coast, there are 10 sites located within 5 km of the ocean.
The study area presents a mostly gradual change in elevation, from 0 m along the coast to a maximum height of 1041 m in the Ventania range. The sites contain elevation values from 9–533 m; with an average of 147.3 m (Table 3). Around half of the sites (n = 34; 51%) have an elevation between 39–208 m (Figure 5), and all of the sites over 200 m are located in the hills.
|ENVIORMENTAL LAYER||TOTAL||HILLS||INTER-HILL PAMPAS||LOW-LYING PAMPAS|
|TWI (index value)||31.2||9.4||4.6||23.1||8.4||4.6||31.2||10.7||5.5||20||10.2||6.2|
|Watercourse (distance from site in km)||34.8||3.1||0.01||16.2||2.8||0.06||9.7||2.4||0.01||34.8||7.6||0.32|
|Bodies of water (distance from site in km)||37.2||8.7||6.1||37.2||11.6||6.1||15.9||4.9||6.2||12.5||6.9||4.1|
|Toolstone (distance from site in km)||248||87.2||0||248||83.2||0||154||91.4||37.2||125||95.8||85.6|
The average slope angle of sites is 3.6°, with a minimum of 0.1° and a maximum of 15.2°. The majority of sites (n = 51; 77%) present slopes between 1.2° to 4.8°. All of the sites above 4.8° correspond to caves and rock shelters located in the hills. Areas with low slope angle are prone to accumulate water, which could be linked to a high topographic wetness index (TWI). The TWI values for the study area range from 3.5 to 32 (mean = 12; sd = 6.2). Half of the sites (50.7%) have values between 6.3 and 9.1, which suggests that sites are mostly located on dry land. Of the 17 sites below the TWI index of 6.3; 14 (82%) correspond to caves and rock shelters in the hills.
Aspect ranged from 38° to 356°, with an average of 190°. Half of the sites (n = 34; 51%) have an aspect between 109° to 266°, which is roughly equal to a northern direction (in the southern hemisphere). There is no clear pattern in the distribution of sites and aspect, as sites in the hills have a similar range of aspect as those in the inter-hill and low-lying Pampas. There are also no outliers in this variable (Figure 5), which suggests a homogenous range in aspect.
The average distance to the watercourses is 3.1 km. The minimum value is slightly overestimated, as the distances were calculated from the center vertex of the watercourse, not the banks or edges, which is variable and dependent on water flow and accumulation during wet and dry seasons. This however, is not an issue with the model because the minimum scale for this analysis is 90 m, so all sites located less than 90 m from the watercourses have a cell value equal to 1. There is at least one outlier from the low-lying Pampas, at 34.8 km.
The average distance to permanent bodies of water is 8.7 km, with a minimum distance of 6.1 m and maximum distance of 37 km. Half of the sites (n = 34; 51%) are between 3.6–12.9 km. With the exception of one site located along the coast, all sites at distances greater than 12.9 km are located in the hills. Similar to watercourses, the minimum value is slightly overestimated, as the polygon vertices are irregular; however, once again this is not an issue for sites under the 90 m cell size (n = 4; 6%).
Finally, the average distance to the primary toolstones sources is 87 km, with a minimum of 0 m (6 sites located directly in the polygon layer) and a maximum of 248 km. The majority of the sites are located between 37–105 km of tool stones sources. Only 9 sites (all in the Tandilia range) are located under the 10 km foraging range from these sources.
From the 75 model runs, the average AICc value was 1725.39 (sd = 322.59). The parameter with the lowest AICc value (1544.83) had a regularization multiplier of 0.75 and used the feature combination of just linear and quadratic. There was no one particular parameter or regulation multiplier that outperformed, however the β = 0.50 had some of the highest values, and the hinge feature on its own did not perform well (Figure 6).
Using the parameters indicated above (β = 0.75; feature classes = LQ), 30 model replicates were run. The average test AUC for the 30 runs was 0.857 (sd = 0.128), meaning that this is a potentially useful model. In other words, the model shows moderate performance, which is able to discriminate between random points and the environment associated with site locations (Franklin 2010; Galletti et al. 2013). All seven of the environmental variables were uncorrelated (r < 0.99) (see Data Availability section). The environmental variable with highest percent contribution gain were watercourses (55.3%), followed by slope (29.7%) (Table 4). The remaining variables contributed to 5.7% or less. The permutation importance (i.e., the measure of the change in model fit when values of a given covariate are randomly shuffled among sites) show qualitatively similar results; with the exception of bodies of water, which increases from 5.4% to 10.6%; suggesting that bodies of water may play an important role in the distribution of sites. Fundamentally, results suggest that areas close to water, in particular watercourses, are the most suitable locations to detect new archaeological sites.
|VARIABLE||PERCENT CONTRIBUTION||PERMUTATION IMPORTANCE|
|Bodies of water||5.4||10.6|
Figure 7 shows the mean response curves of the 30 replicate Maxent runs (red) and the mean +/– one standard deviation (blue) of the three variables with the highest permutation importance (watercourses, slope, and bodies of water). These curves show how each environmental variable influences the prediction of the hunter-gatherer site location. The x-axis represents the grid cell values of the environmental variable (e.g., distance in meters), and the y-axis gives the probability of occurrence (cloglog output) on a scale from 0 (low probability) to 1 (high probability). The distance to watercourses shows an abrupt curve, in both the marginal and corresponding response curves, staring at grid cell values of 0.60 (with a less than 10% probability). Grid cells with a value of 0.95 to 1.00 (which can be roughly translated to a distance 40 meters or less from the watercourse) have 80% or more probability of sites. This shows that suitable distance for hunter-gatherer archaeology sites is less than 40 meters from watercourses.
The response curve for slope suggests higher probably of sites occurring between 10° to 30°. However, as mentioned above, almost 80% of the presence-only sites are located on slopes under 5°, and no sites are located on slopes greater than 15.2°. Looking at the response curve created using only the corresponding variable (slope), we see how site probably begins to decrease after 20°. Therefore, this corresponding model appears to more accurately describe the response of sites to slope.
A less abrupt curve is shown for bodies of water, with a general correlation between grid cell values and probability. For example, sites with grid cell values between 0.80 and 0.90 (which can be roughly translated to a distance 10 km from the body of water), have a site probability of close to 85%. This is consistent with the average distance to bodies of water for the study area, which is 8.7 km. When using only the corresponding variable, these values decrease (e.g., cell values of 0.80 = 60% probability). Both curves show how this variable plays a lesser role in the contribution to the model.
The Maxent model generated a map using the mean and standard deviation of the 30 output grids (Figure 8). Visual analysis once again suggests that watercourses are the primary contributor to the model. The average AUC value for the individual sites across the study area was 0.74 (sd = 0.25). The majority (n = 40; 59.7%) contained AUC values greater than 0.80 (Figure 9). Of these, 4 sites (all corresponding to sites in the hills; specifically, caves/rock shelters) had AUC values of 1.00; 22 sites (13 in the hills, and 9 in the inter-hill Pampas) had AUC values between 0.90 and 0.99; and 13 sites (7 in the hills, 6 in inter-hill Pampas) had AUC values between 0.80 and 0.89. A Wilcox rank sum test was unable to reject the null hypothesis using the AUC values from the hills and the inter-hill Pampas (W = 523, p-value = 0.2876; r = 0.137; see Data Availability section for Rscript and raw data). The environmental region with the lowest AUC values was the low-lying Pampas, with a maximum value of just 0.60. Sample size is likely an issue in this region.
The pedestrian transects that were completed previous to the Maxent model included a total of 104 grid cells. Only 11 grid cells contained archaeological material in the form of isolated remains (4 grid cells) or clusters (7 grid cells). The mean AUC for the grid cells with archaeological material was 0.68. The remaining 93 grid cells, which contained no archaeological material had a similar mean AUC value (0.66). Over half of these grid cells (n = 57; 55%) had AUC values between 0.60 and 0.80; and just 14 grid cells (13%) had AUC values above 0.80. This suggests the probability of finding archaeological materials in this survey area is closer to random. Looking at some specific transects, the model did successfully predict the location of archaeological materials along grid cells near watercourses; as well as locations with no archaeological materials (see example in Figure 10). In other words, the model does appear to be useful at finer scales of analysis (localized areas), identifying environmental variables such as bodies of water with low average probabilities.
Despite a non-optimal regional environment for predictive models, and a challenging archaeological record (e.g., non-representative sampling), the Maxent model produced a reliable exploratory model for the study area. This is promising not only for archaeology predictive models, but also for the Maxent model performance in this type of environment (Jones 2017; Noviello et al. 2018). The results show that the key environmental factors influencing the habitat distribution of hunter-gatherer archaeological sites in this region are watercourses; specifically, classified grid cells which are close or less than 40 meters from this variable. Additionally, sites with slope between 5° and 20°, and that are less than 10 km from a permanent body of water are also promising. What seems to be less relevant to site distribution are the distances to high quality toolstones, ground wetness, aspect, and elevation.
The influence of watercourses in archaeological site distribution from the southern Pampas has two possible, but not mutually exclusive, explanations: (1) hunter-gatherer settlement choices, and (2) research bias. First, hunter-gatherers logically occupied areas near watercourses for resource availability, including fresh water, and associated animal and plants resources. The watercourses also likely acted as natural travel routes (waterway networks) for groups. In many cases, these watercourses are interconnected with permanent bodies of water, and both created the structure for a foraging territory (Conolly 2018; Kelly 1995). In regards to hunting strategies, the water attracted different mammals, like the guanaco (the most prominent game recovered in archaeological sites in the region), providing hunter-gatherers the opportunity to stalk or trap their prey using these natural features of the landscape (Kaufmann et al. 2021). Many archaeological sites from the inter-hill Pampas near water show evidence of repeated and frequent use, and are interpreted as persistent places (sensu Schlanger 1992: 92) (e.g., Arroyo Seco 2, Paso Mayor, and Las Toscas 5; Bayón et al. 2010; Massigoge et al. 2021; Politis et al. 2016). While bodies of water did not show high relative contributions to the Maxent model, they were clearly a vital resource. One positive characteristic that could be taken from the model is its effectiveness to identify bodies of water with low probability, as shown in the field validation results. This is potentially an effective method to deciding which bodies of water to survey, which is particularly challenging in the inter-hill Pampas, considering the immense number of shallow lakes.
A second possibility for the abundance of sites near watercourses is bias related. While the bias grid used in the Maxent model was designed to counter this variable, the fact that pedestrian survey along watercourses is one of the most widely used methods in the region for the detection of archaeological sites is a good reason for the abundance of sites associated with this variable. One option would be to create a bias grid that used not only the location of known archaeological sites, but also the location of watercourses. This can be tested in future models; however, the expected result will reduce the contribution of water resources, which as mentioned above, was an important variable in hunter-gatherer ways of life. To avoid this bias, a better understanding of the paleo landscape; in particular the paleochannel network during different periods, may be a possible solution.
Slope appears to be relevant for the distribution of archaeological sites, although this seems to play a minor role in the inter-hill and low-lying Pampas. Sites in the hills are naturally located on sloped terrain, and thus show higher slope gradient. From a hunter-gatherer perspective, settlement areas (habitable terrain) and travel routes (the cost of movement) were constrained by slope (Llobera 2000). The slope of the land can also impact behavior of grazing animals and act as a natural feature for trapping prey (similar to watercourses). According to the Maxent model, archaeologists working in the southern Pampas should focus on classified grid cells with less than a 20° slope gradient. Slope is directly related to the flatness of the landscape, and hunter-gatherers will prefer dry terrain for settlement areas. The Maxent model showed a very low percentage contribution to the topographic wetness. TWI has been shown to perform badly in a flat landscape and areas characterized by gentle slope (Mattivi et al. 2019). The flat landscape of the inter-hill Pampas may also have reduced the contribution of the elevation variable. Regardless of the fact that there are a greater number of sites in the hills, elevation appears to be less relevant when predicting site locations in the southern Pampas. This result is somewhat inconsistent with the Maxent predictive model by Politis et al. (2019) in the Lower Paraná River. Their research suggested that elevation is the most relevant environmental variable. Since the environment in that region is mostly fluvial, settlement elevation is considered a relevant part of an adaptation to these types of environments (Politis et al. 2019). Conversely, the Maxent model for the southern Pampas region suggests that hunter-gatherers preferred lower elevations, which, as some authors propose would be better for domestic activities (Mazzanti 2003). Furthermore, most locations along the hills are easily accessible, and social practices were not necessarily related to altitude (Mazzia & Flegenheimer 2012).
Regarding the distance to the high-quality raw materials, the low contribution of this variable is consistent with the results of previous regional archaeological studies, which suggests that hunter-gatherers from the Pampas region kept themselves supplied with artifacts through diverse provisioning strategies or used alternative raw materials from outcrops located near the sites (Messineo & Barros 2015). If we consider the daily foraging trip of ca. 10 km, and the known location of archaeological sites, grassland hunter-gatherers at present day conditions would have had good access to fresh watercourses (average distance = 3.1 km) or bodies of water (average distance = 8.7 km); but the distance to high-quality SBG orthoquartzite would not be accessible in most areas (average distance = 87 km). The average distance between archaeological sites in the southern Pampas is currently 163 km, meaning that hunter-gatherers accessed raw materials through either provisioning places, or provisioning individuals of long-distance exchange networks (Armentano, Martínez & Gutiérrez 2007; Barros 2012; Flegenheimer, Mazzia & Weitzel 2015; Martínez & Mackie 2003).
One last factor to consider when interpreting the results of the Maxent model is the detectability issues. On the one hand, portions of the geographic range where the population is most dynamic (low densities, high turnover, etc.) are precisely where non-detection is most likely to occur (Yackulic et al. 2013). Regions of rapid change in occupancy with corresponding low-detection probabilities are expected to occur in archaeology. Hunter-gatherer groups form the southern Pampas region were likely highly mobile, small sized family groups, and occupied short-term camps and hunting stands; therefore, non-detection is probable in most locations. The occurrence of greater seasonal or interannual fluctuations in water availability in smaller basins could have contributed to a shorter duration and high turnover of human occupations (Massigoge 2009; Massigoge & Pal, 2011). On the other hand, some watercourses have, and continue to undergo significant fluvial processes that negatively affect the preservation and visibility of archaeological deposits, which help explain in part the lack of buried archaeological material in certain sectors of the landscape (Dubois 2006; Favier-Dubois, Massigoge & Messineo 2017).
The Maxent model is a flexible cost-effective way to identify areas most likely to reveal the presence of archaeological sites. Initial results of the model for the southern Pampas are promising, with a moderate performance. When complete random sampling is not feasible, this information can be used by archaeologists to focus their efforts on areas with greater probabilities of site discovery. In other words, this method should be considered a complement, and not a replacement to traditional archaeological survey methods. Random sampling strategies should not be discarded, as this is needed to improve the overall model performance. Archaeologists must also reinforce a standardized field transect recording method. Recording track survey data with information of occurrence data (i.e., positive and negative grid cells) will help to improve field validation. The Maxent model can also be used to complement heritage cultural management strategies.
With new environmental and archaeological data regularly being gathered and updated, the model needs to be constantly adjusted. Paleoenvironmental data will help improve the performance of the model, in particular, geographic mapping of ancient watercourses. In order to capture the range of variability in the use of the environment and settlement strategies, future analysis should develop a period specific approach (although this may reduce sample size significantly), in addition to site classifications (e.g., base camp vs hunting camps; mortuary vs non-mortuary sites; etc.), and the correlation between distance to major watercourses and the proportions of materials. This will help to test hypothesis and explain archaeological patterns along the watercourses (e.g., waterway-based network models; see Conolly, 2018). The model should also attempt to consider non-environmental variables. Many times, human settlement selection is based on cultural factors, such as symbology or taboo, or social factors, such as the presence of other groups. For this reason, we should be cautious and modest in our expectations for inference when using predictive modelling based exclusively on environmental variables (Grøn 2018; Wheatley 2004). Currently, site databases for the southern Pampas are inconsistent, where some site records are very detailed and provide data from sampling strategies, whereas other site records are basic qualitative descriptions. Sample size is still geographically sparse for certain areas (e.g., north of the Tandilia hills for instance), which makes interpretations problematic.
In conclusion, the site and environmental distributions described in this paper, along with a new Maxent predictive model, allowed us to start to identify key landscape variables critical to hunter-gatherer settlement decisions in the southern Pampas. This information can ultimately lead to the detection of new archaeology sites. This is an important advance in the regional analysis, and for similar environments around the world (i.e., grassland hunter-gatherer settings). This work also reinforced the need for using standardized protocols when running and analyzing distribution models. In order for the model to improve we need to continue to use these protocols to enhance transparency, reproducibility, evaluation and reuse in the research (Merow, Smith & Silander 2013; Yackulic et al. 2013; Zurell et al. 2020: 1273). We also need to continue to study the uncertainty of the model through more calibration (adjustment of the input parameters), sensitivity analysis (examination of the effects of different input parameters), and model validation (evaluation of how well the model provides reproducible results that are close to reality) (Brouwer Burg 2017; Brouwer Burg, Peeters & Lovis 2016).
Datasets related to this article can be found at https://github.com/drafuse/maxent-grassland-archaeology, hosted at https://github.com/.
Survey transect data from the southern Pampas was processed in collaboration with Agustina Massigoge. Thanks to Agustina for her review and corrections of the manuscript, and to James Conolly and Jamie Kass for their help with Maxent. Finally, thanks to the two external journal reviewers for their comments and suggestions.
This work was funded by the Unidad Ejecutora INCUAPA-CONICET (Facultad de Ciencias Sociales, UNICEN), and the Agencia Nacional de Promoción Científica y Tecnológica (ANPCyT), under Grant numbers PICT-2015-2777 and PICT-2017-0433.
The author has no competing interests to declare.
Aguirre, ML and Whatley, RC. 1995. Late Quaternary marginal marine deposits and palaeoenvironments from northeastern Buenos Aires Province, Argentina: A review. Quaternary Science Reviews, 14(3): 223–254. DOI: https://doi.org/10.1016/0277-3791(95)00009-E
Araújo, MB and Guisan, A. 2006. Five (or so) challenges for species distribution modelling. Journal of Biogeography, 33(10): 1677–1688. DOI: https://doi.org/10.1111/j.1365-2699.2006.01584.x
Armentano, G, Martínez, GA and Gutiérrez, MA. 2007. Revisión del sitio Paso Otero 5: aspectos tecnológicos y fuentes de aprovisionamiento. In Bayón, C, Pupio, A, González, MI, Flegenheimer, N and Frere, M (eds.). Arqueología En Las Pampas, 535–548. Tomo II. Sociedad Argentina de Antropología.
Banning, EB. 2021. Sampled to Death? The Rise and Fall of Probability Sampling in Archaeology. American Antiquity, 86(1): 43–60. DOI: https://doi.org/10.1017/aaq.2020.39
Barbet-Massin, M, Jiguet, F, Albert, CH and Thuiller, W. 2012. Selecting pseudo-absences for species distribution models: How, where and how many? Methods in Ecology and Evolution, 3(2): 327–338. DOI: https://doi.org/10.1111/j.2041-210X.2011.00172.x
Barrientos, G, Catella, L and Oliva, F. 2015. The Spatial Structure of Lithic Landscapes: the Late Holocene Record of East-Central Argentina as a Case Study. Journal of Archaeological Method and Theory, 22: 1151–1192. DOI: https://doi.org/10.1007/s10816-014-9220-0
Barros, MP. 2012. La circulación y el modo de explotación de las rocas en el sitio Paso Otero 4 durante el Holoceno temprano y medio (partido de Necochea, provincia de Buenos Aires, Argentina). In Boëda, E, Farias, M and Lourdeau, A (eds.), Peuplement et Modalités d’occupation de l’Amérique Du Sud: L’apport de La Technologie Lithique/Povoamento e Modalidades de Ocupação Humana Na América Do Sul: A Contribuição Da Tecnologia Lítica, 255–278. Couverture P. Binant pour @rchéo-éditions.com.
Barros, MP, Messineo, PG and Colantonio, MJ. 2015. Chert quarries and workshops in the Humid Pampa sub-region: New contributions on exploitation techniques and circulation through study of chaînes opératoires. Quaternary International, 375: 99–112. DOI: https://doi.org/10.1016/j.quaint.2014.07.017
Barve, N, Barve, V, Jiménez-Valverde, A, Lira-Noriega, A, Maher, SP, Peterson, AT, Soberón, J and Villalobos, F. 2011. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecological Modelling, 222(11): 1810–1819. DOI: https://doi.org/10.1016/j.ecolmodel.2011.02.011
Bayón, C and Flegenheimer, N. 2004. Cambio de planes a través del tiempo para el traslado de roca en la pampa bonaerense. Estudios Atacameños, 28: 59–70. DOI: https://doi.org/10.4067/S0718-10432004002800006
Bayón, C, Flegenheimer, N and Pupio, A. 2006. Planes sociales en el abastecimiento y traslado de roca en la pampa bonaerense en el Holoceno temprano y tardío. Relaciones de la Sociedad Argentina de Antropología, XXXI: 19–45.
Bayón, C, Flegenheimer, N, Valente, M and Pupio, A. 1999. Dime cómo eres y te diré de dónde vienes: procedencia de rocas cuarcíticas en la Región Pampeana. Relaciones de la Sociedad Argentina de Antropología, XXIV: 187–217.
Bettinger, RL. 1981. Explanatory/Predictive Models of Hunter-Gatherer Adaptation. In: Schiffer, MB (ed.), Advances in Archaeological Method and Theory, 157–223. San Diego: Academic Press. DOI: https://doi.org/10.1016/B978-0-12-624180-8.50009-3
Bonadonna, FP, Leone, G and Zanchetta, G. 1999. Stable isotope analyses on the last 30 ka molluscan fauna from Pampa grassland, Bonaerense region, Argentina. Palaeogeography, Palaeoclimatology, Palaeoecology, 153(1): 289–308. DOI: https://doi.org/10.1016/S0031-0182(99)00063-2
Brouwer Burg, M. 2017. It must be right, GIS told me so! Questioning the infallibility of GIS as a methodological tool. Journal of Archaeological Science, 84: 115–120. DOI: https://doi.org/10.1016/j.jas.2017.05.010
Brouwer Burg, M, Peeters, H and Lovis, WA. 2016. Uncertainty and Sensitivity Analysis in Archaeological Computational Modeling. Springer. DOI: https://doi.org/10.1007/978-3-319-27833-9
Cingolani, CA. 2011. The Tandilia system of Argentina as a southern extension of the Río de la Plata craton: An overview. International Journal of Earth Sciences, 100(2): 221–242. DOI: https://doi.org/10.1007/s00531-010-0611-5
Cione, AL, Tonni, EP and Soibelzon, LH. 2003. The Broken Zig-Zag: Late Cenozoic large mammal and tortoise extinction in South America. Revista del Museo Argentino de Ciencias Naturales, 5(1): 1–19. DOI: https://doi.org/10.22179/REVMACN.5.26
Colombo, M. 2011. El área de abastecimiento de las ortocuarcitas del grupo Sierras Bayas y las posibles técnicas para su obtención entre los cazadores y recolectores pampeanos. Intersecciones en Antropología, 12: 231–243.
Conolly, J, Manning, K, Colledge, S, Dobney, K and Shennan, S. 2012. Species distribution modelling of ancient cattle from early Neolithic sites in SW Asia and Europe. The Holocene, 22(9): 997–1010. DOI: https://doi.org/10.1177/0959683612437871
Cruz, JA, Prado, JL and Arroyo-Cabrales, J. 2021. The mutual ecogeographical range and paleoclimatic reconstruction during the Late Pleistocene-Holocene in the Pampas (Argentina) using meso and microvertebrate fossils. Holocene, 31: 983–992. DOI: https://doi.org/10.1177/0959683621994652
Dubois, CMF. 2006. Dinámica fluvial, paleoambientes y ocupaciones humanas en la localidad arqueológica Paso Otero, río Quequén Grande, Provincia de Buenos Aires Fluvial dynamics, palaeoenvironment and human occupations in the Paso Otero archaeological locality, Quequén Gran. Intersecciones en antropología, 7: 109–127.
Elith, J. 2000. Quantitative Methods for Modeling Species Habitat: Comparative Performance and an Application to Australian Plants. In: Ferson, S and Burgman, M (eds.), Quantitative Methods for Conservation Biology, 39–58. New York: Springer-Verlag. DOI: https://doi.org/10.1007/0-387-22648-6_4
Elith, J, Kearney, M and Phillips, S. 2010. The art of modelling range-shifting species. Methods in Ecology and Evolution, 1(4): 330–342. DOI: https://doi.org/10.1111/j.2041-210X.2010.00036.x
Elith, J and Leathwick, JR. 2009. Species distribution models: Ecological explanation and prediction across space and time. Annual Review of Ecology, Evolution, and Systematics, 40: 677–697. DOI: https://doi.org/10.1146/annurev.ecolsys.110308.120159
Favier-Dubois, CM, Massigoge, A and Messineo, PG. 2017. El Holoceno medio en valles fluviales del sudeste pampeano: ¿escasez de sitios o de unidades portadoras? Una perspectiva geoarqueologica. Revista del Museo de Antropología, 10(2): 19–34. DOI: https://doi.org/10.31048/1852.4826.v10.n2.15414
Fick, SE and Hijmans, RJ. 2017. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology, 37(12): 4302–4315. DOI: https://doi.org/10.1002/joc.5086
Flegenheimer, N, Mazzia, N and Weitzel, C. 2015. Landscape and Rocks in the East-Central Portion of the Tandilia Range (Buenos Aires Province, Argentina). PaleoAmerica, 1(1): 163–180. DOI: https://doi.org/10.1179/2055556315Z.00000000017
Fourcade, Y, Engler, JO, Rödder, D and Secondi, J. 2014. Mapping Species Distributions with MAXENT Using a Geographically Biased Sample of Presence Data: A Performance Assessment of Methods for Correcting Sampling Bias. PLoS ONE, 9(5): e97122. DOI: https://doi.org/10.1371/journal.pone.0097122
Franklin, J. 2010. Mapping Species Distributions. Cambridge: Cambridge University Press. DOI: https://doi.org/10.1017/CBO9780511810602
Franklin, J, Potts, AJ, Fisher, EC, Cowling, RM and Marean, CW. 2015. Paleodistribution modeling in archaeology and paleoanthropology. Quaternary Science Reviews, 110: 1–14. DOI: https://doi.org/10.1016/j.quascirev.2014.12.015
Galletti, CS, Ridder, E, Falconer, SE and Fall, PL. 2013. Maxent modeling of ancient and modern agricultural terraces in the Troodos foothills, Cyprus. Applied Geography, 39(May 2015): 46–56. DOI: https://doi.org/10.1016/j.apgeog.2012.11.020
Google Earth Inc. 2020. Google Earth Pro. https://www.google.com/earth/. 22.214.171.12486 (64-bit).
Grill, S, Borromei, A, Martínez, GA, Gutiérrez, MA, Cornou, ME and Olivera, D. 2007. Palynofacial analysis in alkaline soils and paleoenvironmental implications: The Paso Otero 5 archaeological site (Necochea district, Buenos Aires province, Argentina). Journal of South American Earth Sciences, 24(1): 34–47. DOI: https://doi.org/10.1016/j.jsames.2007.02.010
Grøn, O. 2018. Some problems with modelling the positions of prehistoric hunter-gatherer settlements on the basis of landscape topography. Journal of Archaeological Science: Reports, 20: 192–199. DOI: https://doi.org/10.1016/j.jasrep.2018.04.034
Guisan, A and Thuiller, W. 2005. Predicting species distribution: Offering more than simple habitat models. Ecology Letters, 8(9): 993–1009. DOI: https://doi.org/10.1111/j.1461-0248.2005.00792.x
Gutiérrez, MA and Martínez, GA. 2008. Trends in the faunal human exploitation during the Late Pleistocene and Early Holocene in the Pampean region (Argentina). Quaternary International, 191(1): 53–68. DOI: https://doi.org/10.1016/j.quaint.2007.09.024
Hastie, T, Tibshirani, R and Friedman, J. 2009. The Elements of Statistical Learning. Data Mining, Inference, and Prediction. New York, NY: Springer. DOI: https://doi.org/10.1007/978-0-387-84858-7
Howey, MCL, Palace, MW and McMichael, CH. 2016. Geospatial modeling approach to monument construction using Michigan from A.D. 1000-1600 as a case study. Proceedings of the National Academy of Sciences of the United States of America, 113(27): 7443–7448. DOI: https://doi.org/10.1073/pnas.1603450113
Howey, MCL, Sullivan, FB Burg, MB and Palace, MW. 2020. Remotely Sensed Big Data and Iterative Approaches to Cultural Feature Detection and Past Landscape Process Analysis. Journal of Field Archaeology, 45(sup1): S27–S38. Taylor & Francis. DOI: https://doi.org/10.1080/00934690.2020.1713435
Iriondo, M, Brunetto, E and Kröhling, D. 2009. Historical climatic extremes as indicators for typical scenarios of Holocene climatic periods in the Pampean plain. Palaeogeography, Palaeoclimatology, Palaeoecology, 283(3–4): 107–119. DOI: https://doi.org/10.1016/j.palaeo.2009.09.005
Iriondo, MH and Garcia, NO. 1993. Climatic variations in the Argentine plains during the last 18,000 years. Palaeogeography, Palaeoclimatology, Palaeoecology, 101(3–4): 209–220. DOI: https://doi.org/10.1016/0031-0182(93)90013-9
Johnson, E, Holliday, VT, Martínez, G, Gutiérrez, M and Politis, GG. 2012. Geochronology and Landscape Development Along the Middle Río Quequén Grande at the Paso Otero Locality, Pampa Interserrana, Argentina. Geoarchaeology, 27(4): 300–323. DOI: https://doi.org/10.1002/gea.21410
Jones, EE. 2017. Significance and context in GIS-based spatial archaeology: A case study from Southeastern North America. Journal of Archaeological Science, 84: 54–62. DOI: https://doi.org/10.1016/j.jas.2017.05.009
Kantner, J. 2008. The archaeology of regions: From discrete analytical toolkit to ubiquitous spatial perspective. Journal of Archaeological Research, 16(1): 37–81. DOI: https://doi.org/10.1007/s10814-007-9017-8
Kelly, RL. 1983. Hunter-Gatherer Mobility Strategies. Journal of Anthropological Research, 39(3): 277–306. DOI: https://doi.org/10.1086/jar.39.3.3629672
Macchi, D. 2010. Modelos predictivos en arqueología: un ejemplo práctico (partidos de General Lavalle y de la Costa, pcia. de Buenos Aires). In: Bárcena, R and Chiavazza, H (eds.), Arqueología Argentina En El Bicentenario de La Revolución de Mayo. XVII Congreso Nacional de Arqueología Argentina, 1787–1792. Tomo 5. Mendoza: Facultad de Filosofía y Letras, Universidad Nacional de Cuyo.
Mancini, MV, Paez, MM, Prieto, AR, Stutz, S, Tonello, M and Vilanova, I. 2005. Mid-Holocene climatic variability reconstruction from pollen records (32°–52°S, Argentina). Quaternary International, 132(1): 47–59. DOI: https://doi.org/10.1016/j.quaint.2004.07.013
Mariano, CI, Endere, ML, Pedrotta, V and Mariano, M. 2014. Anatomía de un Sistema de Información Geográfica (SIG) para el patrimonio arqueológico del centro de la provincia de Buenos Aires. Comechingonia. Revista de Arqueología, 225–247. DOI: https://doi.org/10.37603/2250.7728.v18.n2.18162
Martínez, G and Gutiérrez, MA. 2018. Early Holocene water well in the Pampas of Argentina: Human responses to water shortage events. The Holocene. DOI: https://doi.org/10.1177/0959683618804643
Martínez, G, Gutiérrez, MA, Messineo, PG, Kaufmann, CA and Rafuse, DJ. 2016. Subsistence strategies in Argentina during the late Pleistocene and early Holocene. Quaternary Science Reviews, 144: 51–65. DOI: https://doi.org/10.1016/j.quascirev.2016.05.014
Martínez, GA and Mackie, Q. 2003. Late Holocene human occupation of the Quequen Grande River valley bottom: settlement systems and an example of a built environment in the Argentine Pampas. Before Farming, 4(1): 1–27. DOI: https://doi.org/10.3828/bfarm.2003.4.1
Martínez, GA, Prates, L, Flensborg, GA, Stoessel, L, Alcaráz, AP and Bayala, PD. 2015. Radiocarbon trends in the Pampean region (Argentina). Biases and demographic patterns during the final Late Pleistocene and Holocene. Quaternary International, 356: 89–110. DOI: https://doi.org/10.1016/j.quaint.2014.09.056
Massigoge, A. 2009. Arqueología de los cazadores-recolectores del sudeste de la región Pampeana: una perspectiva tafonómica. Unpublished Ph.D thesis. Universidad Nacional de La Plata, La Plata, Argentina. DOI: https://doi.org/10.35537/10915/22949
Massigoge, A and Pal, N. 2011. Producción y uso de artefactos líticos en contextos cazadores-recolectores del Holoceno tardío del área Interserrana (Argentina): análisis integral de la diversidad tecno-morfológica y funcional. Revista Española de Antropología Americana, 41(1): 51–73. DOI: https://doi.org/10.5209/rev_REAA.2011.v41.n1.3
Massigoge, A, Rafuse, DJ, Rodriguez, JM, Torino, RF and Favier-Dubois, CM. 2021. Technology, subsistence and mobility of Middle to Late Holocene hunter-gatherers in the Southern Pampas: new archaeological data from Las Toscas shallow lake (Argentina). Archaeological and Anthropological Sciences, 13(69). DOI: https://doi.org/10.1007/s12520-021-01314-3
Mattivi, P, Franci, F, Lambertini, A and Bitelli, G. 2019. TWI computation: a comparison of different open source GISs. Open Geospatial Data, Software and Standards, 4(1). Open Geospatial Data, Software and Standards. DOI: https://doi.org/10.1186/s40965-019-0066-y
Mazzanti, DL. 2003. Human settlements in caves and rockshelters during the Pleistocene-Holocene Transition in the Eastern Tandilia Range, Pampean Region, Argentina. In: Miotti, L, Salemme, M and Flegenheimer, N (eds.), Where the South Winds Blow. Ancient Evidence of Paleo South Americans, 57–63. Texas A&M University, Collage Station: Center for the Study of the First Americans.
Mazzanti, DL. 2006. La constitución de territorios sociales durante el Holoceno Tardío. El caso de las sierras orientales de Tandilia, Argentina. Relaciones de la Sociedad Argentina de Antropología, XXXI: 277–300.
Mazzanti, DL and Bonnat, GF. 2013. Paisajes Arqueológicos Y Cazadores-Recolectores De La Transición Pleistoceno-Holoceno. Análisis De Las Cuencas De Ocupación En Tandilia Oriental, Provincia De Buenos Aires, Argentina. Relaciones de la Sociedad Argentina de Antropología, XXXVIII(2): 521–541.
Mazzia, N and Flegenheimer, N. 2012. Early settlers and their places in the Tandilia range (Pampean Region, Argentina). In: Miotti, LL, Salemme, M, Flegenheimer, N and Goebel, T (eds.), Southbound. Late Pleistocene Peopling of Latin America, 105–110. Texas A&M University, Texas: Current Research in the Pleistocene. Center for the Study of the First Americans.
Mazzia, N and Gómez, JC. 2013. GIS and Landscape Archaeology: A Case of Study in the Argentine Pampas. International Journal of Heritage in the Digital Era, 2(4): 527–546. DOI: https://doi.org/10.1260/2047-49126.96.36.1997
McMichael, CH, Palace, MW, Bush, MB, Braswell, B, Hagen, S, Neves, EG, Silman, MR, Tamanaha, EK and Czarnecki, C. 2014. Predicting pre-Columbian anthropogenic soils in Amazonia. Proceedings of the Royal Society B: Biological Sciences, 281(1777). DOI: https://doi.org/10.1098/rspb.2013.2475
Merckx, B, Steyaert, M, Vanreusel, A, Vincx, M and Vanaverbeke, J. 2011. Null models reveal preferential sampling, spatial autocorrelation and overfitting in habitat suitability modelling. Ecological Modelling, 222: 588–597. DOI: https://doi.org/10.1016/j.ecolmodel.2010.11.016
Merow, C, Smith, MJ and Silander, JA. 2013. A practical guide to MaxEnt for modeling species’ distributions: What it does, and why inputs and settings matter. Ecography, 36(10): 1058–1069. DOI: https://doi.org/10.1111/j.1600-0587.2013.07872.x
Messineo, PG and Barros, MP. 2015. Lithic Raw Materials and Modes of Exploitation In Quarries and Workshops from the Center of the Pampa Grasslands of Argentina. Lithic Technology, 40(1): 3–20. DOI: https://doi.org/10.1179/2051618514Y.0000000007
Miller, J. 2010. Species Distribution Modeling. Geography Compass, 4(6): 490–509. DOI: https://doi.org/10.1111/j.1749-8198.2010.00351.x
Morales, NS, Fernández, IC and Baca-González, V. 2017. MaxEnt’s parameter configuration and small samples: Are we paying attention to recommendations? A systematic review. PeerJ, 2017(3): 1–16. DOI: https://doi.org/10.7717/peerj.3093
Morgan, C. 2008. Reconstructing prehistoric hunter-gatherer foraging radii: a case study from California’s southern Sierra Nevada. Journal of Archaeological Science, 35(2): 247–258. DOI: https://doi.org/10.1016/j.jas.2007.02.025
Muscarella, R, Galante, PJ, Soley-Guardia, M, Boria, RA, Kass, JM, Uriarte, M and Anderson, RP. 2014. ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for Maxent ecological niche models. Methods in Ecology and Evolution, 5(11): 1198–1205. DOI: https://doi.org/10.1111/2041-210X.12261
Muttaquin, LA, Heru, SH and Susilo, B. 2019. MaxEnt (Maximum Entropy) model for predicting prehistoric cave sites in Karst area of Gunung Sewu, Gunung Kidul, Yogyakarta. In: Wibowo, SB, Rimba, AB, Aziz, A, Phinn, S, Sri Sumantyo, JT, Widyasamratri, H and Arjasakusuma, S (eds.), Sixth Geoinformation Science Symposium, 3. SPIE. DOI: https://doi.org/10.1117/12.2543522
Noviello, M, Cafarelli, B, Calculli, C, Sarris, A and Mairota, P. 2018. Investigating the distribution of archaeological sites: Multiparametric vs probability models and potentials for remote sensing data. Applied Geography, 95: 34–44. DOI: https://doi.org/10.1016/j.apgeog.2018.04.005
Orton, C. 2000. Sampling in Archaeology. Cambridge University Press. DOI: https://doi.org/10.1017/CBO9781139163996
Pearce, JL and Boyce, MS. 2006. Modelling distribution and abundance with presence-only data. Journal of Applied Ecology, 43(3): 405–412. DOI: https://doi.org/10.1111/j.1365-2664.2005.01112.x
Perkins-Taylor, IE and Frey, JK. 2020. Predicting the distribution of a rare chipmunk (Neotamias quadrivittatus oscuraensis): Comparing MaxEnt and occupancy models. Journal of Mammalogy, 101(4): 1035–1048. DOI: https://doi.org/10.1093/jmammal/gyaa057
Peterson, AT, Soberón, J, Pearson, RG, Anderson, RP, Martínez-Meyer, E, Nakamura, M and Araújo, MB. 2011. Ecological Niches and Geographic Distributions (MPB-49). Princeton University Press. DOI: https://doi.org/10.23943/princeton/9780691136868.001.0001
Phillips, SJ. 2017. A Brief Tutorial on Maxent. Available from url: http://biodiversityinformatics.amnh.org/open_source/maxent/.
Phillips, SJ, Anderson, RP, Dudík, M, Schapire, RE and Blair, ME. 2017. Opening the black box: an open-source release of Maxent. Ecography, 40(7): 887–893. DOI: https://doi.org/10.1111/ecog.03049
Phillips, SJ, Anderson, RP and Schapire, RE. 2006. Maximum entropy modeling of species geographic distributions. Ecological Modelling, 190(3–4): 231–259. DOI: https://doi.org/10.1016/j.ecolmodel.2005.03.026
Phillips, SJ and Dudík, M. 2008. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography, 31(2): 161–175. DOI: https://doi.org/10.1111/j.0906-7590.2008.5203.x
Phillips, SJ, Dudík, M, Elith, J, Graham, CH, Lehmann, A, Leathwick, J and Ferrier, S. 2009. Sample selection bias and presence-only distribution models: Implications for background and pseudo-absence data. Ecological Applications, 19(1): 181–197. DOI: https://doi.org/10.1890/07-2153.1
Phillips, SJ, Dudík, M and Schapire, RE. 2021. Maxent software for modeling species niches and distributions (Vesion 3.4.1). Available from url: http://biodiversityinformatics.amnh.org/open_source/maxent/.
Politis, G, Domic, A, Bonomo, M and Capriles, JM. 2019. Modelamiento de los parámetros ambientales de la distribución de sitios goya-malabrigo mediante maxent. Instituto de Arqueología, Lingüística y Folklore “Dr. Pablo Cabrera”, 13(May): 289–304.
Politis, GG. 2008. The Pampas and Campos of South America. In: Silverman, H and Isbell, W (eds.). Handbook of South American Archaeology, 235–260. New York, NY: Springer New York. DOI: https://doi.org/10.1007/978-0-387-74907-5_14
Politis, GG, Gutiérrez, MA, Rafuse, DJ and Blasi, A. 2016. The Arrival of Homo sapiens into the Southern Cone at 14,000 Years Ago. PLOS ONE, 11(9): e0162870. DOI: https://doi.org/10.1371/journal.pone.0162870
Politis, GG, Messineo, PG and Kaufmann, CA. 2004. El poblamiento temprano de las llanuras pampeanas de Argentina y Uruguay. Complutum, 15: 207–224. Available at: http://revistas.ucm.es/index.php/CMPL/article/view/30767.
Politis, GG, Prates, L, Merino, ML and Tognelli, MF. 2011. Distribution parameters of guanaco (Lama guanicoe), pampas deer (Ozotoceros bezoarticus) and marsh deer (Blastocerus dichotomus) in Central Argentina: Archaeological and paleoenvironmental implications. Journal of Archaeological Science, 38(7): 1405–1416. DOI: https://doi.org/10.1016/j.jas.2011.01.013
Politis, GG, Scabuzzo, C and Tykot, RH. 2009. An approach to pre-hispanic diets in the Pampas during the Early/Middle Holocene. International Journal of Osteoarchaeology, 19(2): 266–280. DOI: https://doi.org/10.1002/oa.1053
Ponce, JF, Rabassa, JO, Coronato, A and Borromei, AM. 2011. Palaeogeographical evolution of the Atlantic coast of Pampa and Patagonia from the last glacial maximum to the Middle Holocene. Biological Journal of the Linnean Society, 103(2): 363–379. DOI: https://doi.org/10.1111/j.1095-8312.2011.01653.x
Prado, JL and Alberdi, MT. 1999. The mammalian record and climatic change over the last 30,000 years in the Pampean Region, Argentina. Quaternary International, 57: 165–174. DOI: https://doi.org/10.1016/S1040-6182(98)00057-3
Prieto, AR. 2000. Vegetational history of the Late glacial–Holocene transition in the grasslands of eastern Argentina. Palaeogeography, Palaeoclimatology, Palaeoecology, 157(3): 167–188. DOI: https://doi.org/10.1016/S0031-0182(99)00163-7
QGIS Development Team. 2021. QGIS Geographic Information System. Available at: http://qgis.osgeo.org.
Quattrocchio, ME, Borromei, AM, Deschamps, CM, et al. 2008. Landscape evolution and climate changes in the Late Pleistocene–Holocene, southern Pampa (Argentina): Evidence from palynology, mammals and sedimentology. Quaternary International, 181(1): 123–138. DOI: https://doi.org/10.1016/j.quaint.2007.02.018
Rafuse, DJ and Massigoge, A. 2018. In search for early peopling evidence: archaeological survey in the upper Quequén Grande River basin (Pampas region, Argentina). 9° Simposio Internacional El Hombre Temprano en América: la gente y sus lugres, 182–184. Necochea, Argentina, 2018. DOI: https://doi.org/10.6067/XCV8447818
Rapela, CW, Fanning, CM, Casquet, C, Pankhurst, RJ, Spalletti, L, Poiré, D and Baldo, EG. 2011. The Rio de la Plata craton and the adjoining Pan-African/brasiliano terranes: Their origins and incorporation into south-west Gondwana. Gondwana Research, 20(4): 673–690. DOI: https://doi.org/10.1016/j.gr.2011.05.001
Rapela, CW, Pankhurst, RJ, Casquet, C, Fanning, CM, Baldo, EG, González-Casado, JM, Galindo, C and Dahlquist, J. 2007. The Río de la Plata craton and the assembly of SW Gondwana. Earth-Science Reviews, 83(1–2): 49–82. DOI: https://doi.org/10.1016/j.earscirev.2007.03.004
Sánchez-Fernández, D, Lobo, JM and Hernández-Manrique, OL. 2011. Species distribution models that do not incorporate global data misrepresent potential distributions: a case study using Iberian diving beetles. Diversity and Distributions, 17(1): 163–171. DOI: https://doi.org/10.1111/j.1472-4642.2010.00716.x
Schlanger, SH. 1992. Recognizing Persistent Places in Anasazi Settlement Systems. In: Rossignol, J and Wandsnider, L (eds.). Space, Time, and Archaeological Landscapes, 91–112. Boston, MA: Springer US. DOI: https://doi.org/10.1007/978-1-4899-2450-6_5
Spiers, JA, Oatham, MP, Rostant, LV and Farrell, AD. 2018. Applying species distribution modelling to improving conservation based decisions: A gap analysis of trinidad and tobago’s endemic vascular plants. Biodiversity and Conservation, 27: 2931–2949. DOI: https://doi.org/10.1007/s10531-018-1578-y
Stutz, S, Borel, CM, Fontana, SL, del Puerto, L, Inda, H, García-Rodriguez, F and Tonello, MS. 2010. Late Holocene climate and environment of the SE Pampa grasslands, Argentina, inferred from biological indicators in shallow, freshwater Lake Nahuel Rucá. Journal of Paleolimnology, 44(3): 761–775. DOI: https://doi.org/10.1007/s10933-010-9450-4
Stutz, S and Prieto, AR. 2003. Modern pollen and vegetation relationships in Mar Chiquita coastal lagoon area, southeastern Pampa grasslands, Argentina. Review of Palaeobotany and Palynology, 126(3–4): 183–195. DOI: https://doi.org/10.1016/S0034-6667(03)00084-8
Tonni, EP, Cione, AL and Figini, AJ. 1999. Predominance of arid climates indicated by mammals in the pampas of Argentina during the Late Pleistocene and Holocene. Palaeogeography, Palaeoclimatology, Palaeoecology, 147(3–4): 257–281. DOI: https://doi.org/10.1016/S0031-0182(98)00140-0
van Proosdij, ASJ, Sosef, MSM, Wieringa, JJ and Raes, N. 2016. Minimum required number of specimen records to develop accurate species distribution models. Ecography, 39(6): 542–552. DOI: https://doi.org/10.1111/ecog.01509
Velasco, JA and González-Salazar, C. 2019. Akaike information criterion should not be a “test” of geographical prediction accuracy in ecological niche modelling. Ecological Informatics, 51: 25–32. DOI: https://doi.org/10.1016/j.ecoinf.2019.02.005
Wachtel, I, Zidon, R, Garti, S and Shelach-Lavi, G. 2018. Predictive modeling for archaeological site locations: Comparing logistic regression and maximal entropy in north Israel and north-east China. Journal of Archaeological Science, 92: 28–36. DOI: https://doi.org/10.1016/j.jas.2018.02.001
Warren, DL, Glor, RE and Turelli, M. 2010. ENMTools: a toolbox for comparative studies of environmental niche models. 1(September 2009): 607–611. DOI: https://doi.org/10.1111/j.1600-0587.2009.06142.x
Warren, DL and Seifert, SN. 2011. Ecological niche modeling in Maxent: The importance of model complexity and the performance of model selection criteria. Ecological Applications, 21(2): 335–342. DOI: https://doi.org/10.1890/10-1171.1
Wheatley, D. 2004. Making space for an archaeology of place. Internet Archaeology, 15(15). http://intarch.ac.uk/journal/issue15/wheatley_inde. DOI: https://doi.org/10.11141/ia.15.10
Wisz, MS, Hijmans, RJ, Li, J, Peterson, AT, Graham, CH and Guisan, A. 2008. Effects of sample size on the performance of species distribution models. Diversity and Distributions, 14(5): 763–773. DOI: https://doi.org/10.1111/j.1472-4642.2008.00482.x
Yackulic, CB, Chandler, R, Zipkin, EF, Royle, JA, Nichols, JD, Campbell Grant, EH and Veran, S. 2013. Presence-only modelling using MAXENT: When can we trust the inferences? Methods in Ecology and Evolution, 4(3): 236–243. DOI: https://doi.org/10.1111/2041-210x.12004
Zárate, M and Blasi, A. 1991. Late pleistocene and holocene loess deposits of the Southeastern Buenos Aires province, Argentina. GeoJournal, 24(2): 211–220. DOI: https://doi.org/10.1007/BF00186018
Zárate, M, Kemp, RA, Espinosa, M and Ferrero, L. 2000. Pedosedimentary and palaeoenvironmental significance of a Holocene alluvial sequence in the southern Pampas, Argentina. The Holocene, 10(4): 481–488. DOI: https://doi.org/10.1191/095968300669846317
Zurell, D, Franklin, J, König, C, Bouchet, PJ, Dormann, CF, Elith, J, Fandos, G, Feng, X, Guillera-Arroita, G, Guisan, A, Lahoz-Monfort, JJ, Leitão, PJ, Park, DS, Peterson, AT, Rapacciuolo, G, Schmatz, DR, Schröder, B, Serra-Diaz, JM, Thuiller, W, Yates, KL, Zimmermann, NE and Merow, C. 2020. A standard protocol for reporting species distribution models. Ecography, 43(9): 1261–1277. DOI: https://doi.org/10.1111/ecog.04960