Next Article in Journal
The Role of El Niño in Driving Drought Conditions over the Last 2000 Years in Thailand
Previous Article in Journal
Influences of West Pacific Sea Surface Temperature on Covarying Eurasian Droughts Since the Little Ice Age
Previous Article in Special Issue
Sparse Radiocarbon Data Confound Culture-Climate Links in Late Pre-Columbian Amazonia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Between Foragers and Farmers: Climate Change and Human Strategies in Northwestern Patagonia

1
Instituto de Evolución, Ecología Histórica y Ambiente (IDEVEA-CONICET & UTN FRSR), Av. Urquizas 350, San Rafael, Mendoza 5600, Argentina
2
Departamento de Historia, Facultad de Filosofía y Letras, Universidad Nacional de Cuyo, Mendoza 5500, Argentina
3
Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales (IANIGLA-CONICET), Mendoza 5500, Argentina
4
Departamento de Historia, Facultad de Historia, Facultad de Filosofía y Letras, Universidad Nacional de Cuyo, Mendoza 5500, Argentina
*
Author to whom correspondence should be addressed.
Quaternary 2020, 3(2), 17; https://doi.org/10.3390/quat3020017
Submission received: 21 December 2019 / Revised: 29 April 2020 / Accepted: 13 May 2020 / Published: 17 June 2020

Abstract

:
In this paper we explore how changes in human strategies are differentially modulated by climate in a border area between hunter-gatherers and farmers. We analyze multiple proxies: radiocarbon summed probability distributions (SPDs), stable C and N isotopes, and zooarchaeological data from northwestern Patagonia. Based on these proxies, we discuss aspects of human population, subsistence, and dietary dynamics in relation to long-term climatic trends marked by variation in the Southern Annular Mode (SAM). Our results indicate that the farming frontier in northwestern Patagonia was dynamic in both time and space. We show how changes in temperature and precipitation over the last 1000 years cal BP have influenced the use of domestic plants and the hunting of highest-ranked wild animals, whereas no significant changes in human population size occurred. During the SAM positive phase between 900 and 550 years cal BP, warmer and drier summers are associated with an increase in C4 resource consumption (maize). After 550 years cal BP, when the SAM changes to the negative phase, wetter and cooler summer conditions are related to a change in diet focused on wild resources, especially meat. Over the past 1000 years, there was a non-significant change in the population based on the SPD.

1. Introduction

For most archaeologists, farming was an evolutionary, progressive, and irreversible process [1,2,3,4]. From this perspective, farming’s origin and expansion led to population increases and occupation of new areas. Farming expanded in a continuous process and only stopped when the environment made farms impossible [5,6], defining a stable frontier that persisted and that would have persisted in time. In historical times, some of these borders were occupied by foragers, but this last scene is masked and omitted in the models [2,7]. Recent research demonstrated that the border area among hunter-gatherers and farming societies never was the case for a lineal non-reversible expansion [8,9,10]. Both subsistence strategies were reversible, and even complementary on incipient farming. Therefore, our objective is to understand the stability of agricultural borders. Under which conditions did the border area perpetuate itself? Was the border unstable and variable over long periods of time (e.g., millennia, centuries)? This paper evaluates the degree to which climatic variations could modulate a retraction in the farming borders. We will focus on northwestern Patagonia, where researcher proposed a pre-Hispanic farming border [11,12].
In the northern sector of the region, the land consisted of enclaves of groups with an agricultural component in their subsistence, beginning ca. 2300 years BP [13,14,15]. The rest of northwestern Patagonia was the land of the hunter-gatherers [11]. In historic chronicles, only hunter-gatherers are mentioned within southern Mendoza [16,17]. However, mostly inferred by domestic plant macrorests, during pre-Hispanic times, farmers might have occupied terrains farther south than the farmer/hunter-gatherer border observed during historical times. If this assertion is right, there was a retraction of farming after it reached its maximum expansion. This contradiction between the historical and archaeological records cast doubt on the idea of a stable long-term border.
In northwestern Patagonia, previous research focused on the last 2000 years cal BP, proposing and discussing an intensification process that developed from the middle of the late Holocene [18]. Most of this research focused on a multi-millennial time scale, assuming few changes after the intensification process started to develop [18,19,20,21]. This lack of fine-grain research on archaeological trends after 2300 years cal BP implicitly generated a static view of this period, implying that little changes occurred in the human-environment system. To assess how stable the farming/hunter-gatherer border was in northwestern Patagonia, we compared the temporal trends in the SAM (Southern Annular Mode) with human population dynamics (based on radiocarbon summed probability distributions (SPDs)), human diet (stable isotope ratios in human bones, namely (δ13C and δ15N), and subsistence (zooarchaeological record).

1.1. Climate and Environment in Northwestern Patagonia during the Last Millennium

This paper focuses on southern Mendoza, in northwestern Patagonia, where we will concentrate on two desert ecosystems: Monte Desert and Patagonia Desert (Figure 1; [22,23,24]). Given the biogeographical similarity between the Altoandino and Patagonia deserts and the low number of archeological samples available for the former, this work includes both groups of samples under the Patagonia Desert unit. The Monte Desert occupies the lowland plain between 200 and 1400 meters above sea level (MASL). It is a temperate region, with summer-dominant rainfall (100–350 mm per year) [22,24,25,26]. The Patagonia Desert is a cold, dry, and wind-swept environment. In the study area, it is restricted to a piedmont fringe bordering the Andean cordillera between 1000 and 2500 MASL (Figure 1). At 1400 MASL, when the Monte Desert limits with the Patagonia Desert, annual precipitation averages 300 mm.
The Southern Annular Mode (SAM) or Antarctic Oscillation is the main mode of climate variability in the middle and high latitudes of the Southern Hemisphere (SH), associated with approximately 30% of the precipitation and temperature variability in the austral hemisphere [27,28]. Its spatial structure comprises pressure-synchronous anomalies of opposite signs between the mid (45° S) and high (60° S) latitudes of the Southern Hemisphere. When pressures are anomalously high around Antarctica, but low at mid-latitudes, the SAM is defined as in its negative phase and vice versa. These pressure changes in the SH modulates the latitudinal position and the strength of the storm tracks associated with the Westerly winds. As a result, many extratropical regions in South America show significant changes in temperature and precipitation relative to the positive or negative phases of the SAM [27,28].

1.2. Human Population Dynamic and Climate in Northwestern Patagonia: Background and Expectations

Hunter-gatherers settled in northwestern Patagonia since the late Pleistocene–early Holocene [30,31,32,33,34]. During the late Holocene, people occupied new areas, including the most inhospitable regions of La Payunia and the high-altitude Andes [35,36]. The intensity of land use increased in terms of both the diversity and density of archaeological assemblages during the last 2300 years cal BP [31,37,38,39].
Subsistence patterns changed over time, showing an increase in the range of exploited taxa [40,41,42]. The guanaco, a wild South American camelid, is the largest animal prey available in this area, and the best ranked in terms of energetic returns [42,43,44,45]. However, there is a decline in the abundance and significance of guanaco after ca. 2000 years cal BP [21,41]. While camelids are regular components of Monte archaeofaunas, smaller prey, such as armadillos, fish, birds, and medium/big size rodents, are also frequent [20,35,41,45]. Rheid eggs are common in both deserts, but they seem to have been more significant in the Monte Desert [46]. Archaeobotanical evidence indicates that algarrobo (Prosopis spp.) and chañar (Geoffroea spp.) were the most common taxa consumed in the Monte Desert, whereas molle (Schinus polygamous) was the most frequently eaten in the Patagonia Desert [47].
Based on the spatial–temporal distribution of domesticated plants, the Atuel and Diamante rivers (Figure 1) were considered as the southern boundary since ca. 2300 years BP, with “low-level food production systems” [11,15,48,49,50]. South of these valleys, the archaeological record includes only scattered evidence of domesticated plants in the vicinity of Cerro Nevado (Zea mays and Cucurbita spp.) [11,15,35,47]. Although microfossils of maize may have also been found in the Patagonian region of northern Neuquén [51] the only unequivocal domesticates known from the Patagonian Desert comes from macro remains specimens of bottle gourds (Lagenaria spp.), and are not quantitatively significant [15,52].
From an evolutionary and ecological perspective, there is a relationship between resource productivity and human population stability [2]. The origin of agriculture is explained as a human response to the unbalance between human demography and resources availability [4,50]. The stability of human populations results from interactions between climate, resources, internal population processes, life history characteristics, and social processes [2,4,53]. In terms of the population dynamic, Zahid et al. [54] argues that “…radiocarbon record shows that transitioning farming societies experienced the same rate of growth as contemporaneous foraging societies. The same rate of growth measured for populations dwelling in a range of environments and practicing a variety of subsistence strategies suggests that the global climate and/or endogenous biological factors, not adaptability to local environment or subsistence practices, regulated the long-term growth of the human population during most of the Holocene.” In addition, Bettinger [2] affirms that pre-industrial agricultural populations did not grow faster than hunter-gatherer populations. This implies that the differences in adaptation, environment, and time are mild influential determinants of population growth. If climate impacted on the resources, the tradeoff between farming or hunter-gathering will change in organizational variables to maintain optimal responses in the demography-resource equation.
In the frame of the diet-breath model [55,56], we accept wild and domestic plants as low-ranked resources for northwestern Patagonia [21,43]. On the other hand, we consider the faunas, mostly guanacos, armadillos, and ostrich eggs, as high-rank resources [43]. In a scenario where population increase leads to pressure on resources, we expect a change in diet, shifting towards greater reliance on plants, including domesticated plants [21]. The diet-breadth model predicts that guanacos and armadillos, the highest ranked items, will always be in the diet. Low-rank resources, such as algarrobo and maize, should be included in the diet only when highest-rank items are scarce or in decline. Two alternative responses under domesticated use can be modeled: low and high yield density of cultigens [57]. Under a low density of domesticates, they will be a small part of the diet and a non-significant increase in human population is expected; but, in a high density cultigen context, we expect domesticates to form a significant part of the diet, with a subsequent increase in human population density.
When climatic change affects resources productivity, we can expect changes in the human resources balance and, consequently, in their diet and population size [55,56,58]. If climatic change affected farming, diminishing the productivity, different human responses can occur [4]. If domesticates are a significant part of the human diet (more than 40% or 50%), a drop in demography could be expected. Freeman et al. [53] proposed that a reversion between maize farmers and wild seed gatherers can be a response when the domestic plant component is less than 40%. Under this scenario no significant population change can be expected. In addition, if climate increase wild resources productivity, or a human demography drop reduces pressure on resources, we would expect a focus on the highest-rank wild resources.

2. Materials and Methods

2.1. SAM Reconstruction

Paleoclimatic records for our study region are very scarce. In order to infer past climate variations in southern Mendoza, we employed a recently developed multi-proxy-based reconstruction of summer (December–January) SAM over the past 1000 years cal BP [59]. To establish the relationship between the SAM and regional climate, spatial correlation patterns were estimated across southern South America between Climatic Research Unit (CRU) temperature and precipitation, and the SAM instrumental index over the common period 1957–2016 [27]. The November–January spatial correlation patterns over South America south of 10° S latitude between temperature and precipitation were estimated using the facilities provided by the KNMI (https://climexp.knmi.nl/start.cgi). During late spring–summer (November–January), climate strongly influences the primary productivity of the desert ecosystems in southern Mendoza–northern Patagonia [60].

2.2. Radiocarbon Dates as Data: Summed Probability Distribution (SPD)

We monitored human population dynamics using the radiocarbon summed probability distributions (SPDs), a tool often used to reconstruct prehistoric populations [30,61,62]. This analysis is based on the “date as data” framework proposed by Rick [63]. The SPD as proxy for human population dynamics has been subject to criticism, including sampling errors, time-dependent taphonomic loss, spatial and temporal differences in site-to-population ratios, and spatial and temporal differences in sampling intensity [54,64,65,66,67,68,69,70,71,72,73]. Crema et al. [61] provide details on how the SPD limitations should be considered for obtaining consistent results.
We calculated SPD using RCarbon, an R software package for the analysis of large collections of radiocarbon dates [74]. Observed SPDs are fitted to different demographic models [74]. Random calendar years dates were generated via Monte-Carlo simulation and subsequently “back-calibrated” into 14C dates to generate null SPDs [61,75]. This approach allows the detection of statistically significant local deviations (i.e., portions of the SPD showing smaller or greater SPD than the null model) as well as the estimation of a global significance test [76]. A non-parametric extension of the hypothesis-testing approach enables the statistical comparison of two or more sets of 14C dates. The null-hypothesis is the equality of the SPDs.
This paper focuses on the last 1000 years cal BP, in which taphonomic loss can be considered as a minor factor affecting our database. Most of the dates are from samples processed in radiocarbon laboratory during the last 20 years (references in Table S1) and from research programs oriented to the study of long-term cultural process with no special temporal focus [15]. We can remark a possible bias against post-Hispanic dates as far as the Hispanic material culture have archaeological markers to date the last 500 years. If so, we can observe fewer dates for the last 500 years. But it is not the case in northwest Patagonia, where the Hispanic material culture is weakly recorded after 100 or 150 years BP [77]. We built the SPD using different spatial scales: northwestern Patagonia, and secondly divided this area into the Monte Desert and Patagonia Desert. We compared the shape of the SPDs of the two deserts with the permutation test described above, using the same calibration procedure and with 10,000 iterations. To avoid the edge effect, we used a radiocarbon database for the last 1400 years cal BP (n = 128; archaeological sites n = 78; Table S1). The radiocarbon dates come from Monte Desert (n = 60; archaeological sites n = 33) and Patagonia Desert (n = 68; archaeological sites n = 45). We included the 14C results obtained from charcoal (n = 78), plants/maize remains (n = 7), bone collagen (n = 52), leather (n = 1), and material without information. The archaeological sites were open air sites (n = 53), caves, rockshelters (n = 21), and sites without a type information (n = 4). Most of the archaeological contexts are human burials loci (n = 34), some are domestic multi-activity (n = 37), others are specific activities sites (n = 5), and few have no information about their context (n = 2). In order to reduce the effect of bias we binned the 14C dates based on the clustering of the mean 14C years BP by archaeological site, using a threshold of 50 years.
This approach is robust to inter-regional differences in sample size, as the comparison is based on the “shape” of the SPDs and not on differences in their absolute magnitudes [74]. We first assessed whether the SPD for the northwestern Patagonia area showed statistically relevant fluctuations when compared against the uniform and the exponential null models, following the procedure described in Timpson et al. [68], using 10,000 Monte-Carlo simulations, and calibrated with the SHCal13 curve. We used “uniform” in order to detect any possible change of the SPD as significant. On the other hand, we use “exponential” assuming it helps to consider taphonomic loss [66]. To detect differences in the population dynamics between the Monte and Patagonia deserts, we ran the permTest function in Rcarbon [74,76,78].

2.3. Bone Collagen Stable Isotope Carbon and Nitrogen

To reconstruct the human diet, we employed stable isotopes from bone collagen and bone structural carbonate [79]. Our data comes from previously published studies of this region [31,80,81] and new samples provided in Table S2. Extraction of bone collagen and structural carbonate was performed in the Museo de Historia Natural de San Rafael [80,82]. The analyses were carried out at the Stable Isotope Ratio Facility for Environmental Research (SIRFER) at the University of Utah, Stable Isotope Facility Wyoming University and Laboratorio de Isótopos Estables para Ciencias Ambientales (LIECA).
Collagen δ13C and δ15N were determined by flash combustion to produce CO2 and N2. In SIRFER the resulting gases were analyzed using a Finnigan Delta Plus isotope ratio mass spectrometer (Finnigan, Bremen, Germany) coupled with a Carlo Erba Model 1110 elemental analyzer (Carlo Erba, Milan, Italy) through a CONFLO III open split interface (Finnigan, Bremen, Germany). Stable carbon and nitrogen isotopic compositions were calibrated relative to the VPDB and AIR scales using USGS40 and USGS41. Both stable isotope measurements and the sample weight percent carbon and nitrogen were obtained from a single sample combustion. Analytic precision is 0.1‰ for carbon isotope ratios (VPDB) and 0.2‰ for nitrogen (AIR). In LIECA, we followed a similar procedure using a Thermo Scientific DELTA V Advantage continuous flow isotope ratio mass spectrometer coupled via ConFlo IV to an Elementar Analyzer Flash 2000. Stable carbon and nitrogen isotopic compositions were calibrated relative to the V-PDB and AIR scales using USGS-40 and USGS-41a. Measurement uncertainty was monitored using in-house collagen standards with well-characterized isotopic compositions: Caffeine LIECA 17 (δ13C = −33.02‰, δ15N = −2.02‰), SRM-14 polar beer bone collagen (δ13C = −13.66‰, δ15N = +21.52‰), and bone collagen LIECA-17 (δ13C = −18.16‰, δ15N = +11.07‰). Precision (u(Rw)) was determined as ±0.17‰ for both δ13C and δ15N on the basis of repeated measurements of calibration standards, check standards, and sample replicates. C:N ratios are in the range of 3.1 to 3.6, indicating generally good preservation [83]. The bone carbonate was sent to the University of Wyoming Stable Isotope Facility (Gas Bench online with a Finnigan Delta Plus XP) for analysis. The ratios of stable isotopes for carbon and oxygen were determined there and reported in values related to the V-PDB standard (analytical precision 0.3‰).
As a consequence of the fractionation between plant tissues and the consumer, bone collagen was 5‰ more enriched than the dietary sources [84]. Bone collagen from individuals with a diet mostly composed of C3 plants exhibit δ13C values between −24‰ and −19‰; individuals with a diet mostly focused on C4 plants have δ13C values between −8.5‰ and −6.5‰ [85]. Carbon isotope signatures on bone collagen reflect mostly the protein component of the diet; they are preferentially routed from dietary sources and incorporated into body tissues [86,87]. Fernandes et al. [88] indicate that 75% of the carbon in bone collagen is derived from dietary protein; the other 25% is derived from carbohydrates and fats. Using a 9.5‰ offset diet–bone carbonate, we expect that carbonate δ13C signals for 100% C3 and 100% C4 consumers would be around −20 to −14‰ and −4‰ to −2‰, respectively.
Nitrogen isotope ratios reflect diets mediated by the metabolic pathways by which nitrogen is fixed in various tissues [86]. The diet is the primary determinant of isotope values. Most animals exhibit an approximately 3‰ shift in δ15N values relative to their diet by the “trophic enrichment factor” (TEF) [89,90]. Domestic plant nitrogen values vary with soil characteristics and nitrogen cycling, as do the nitrogen values of many other wild plants.
This paper adopts the following fractionation values for analysis (resources–consumer bone collagen): δ13C = 1.0‰ and δ15N = 4.6‰ for animals [91,92,93]; δ13C = 3.1‰ and δ15N = 4‰ for plants [91,92]. For δ13Ccarb the offsets used were δ13C = 8.0‰ and δ15N = 4.0‰ for animals [91,92] and δ13C = 11.1 and δ15N = 4.6‰ for plants. Modern samples were corrected for the Suess effect by adding 1.5‰ to their δ13C values [94,95].
The human bone stable isotope values database includes 43 individuals (Table S2), and only adults and juveniles are included. All samples’ C:N ratios were between the accepted threshold [83]. Radiocarbon dates were calibrated using the OxCal/SHCal13 curve and reported as median ages plus or minus a 68.7‰ probability [96]. Few bone samples have chronology inferred from an associated radiocarbon date. In order to improve the explanation of human bone stable isotope values (δ13C and δ15N) in terms of the diet, we employed regional stable isotopes data on resources [25,80,97]. Using these resources as a reference framework, we explore in a more solid way the meaning of variation in δ13C and δ15N.
In addition, we explored the distribution of stable isotopes in human bone to define a spatial pattern and to discuss variations in our expectations. Q-GIS was used to generate an “isoscape” [98,99] base on δ13Ccarb and δ15N from human bone using the temporal frame defined by SAM. We created a raster surface model using the inverse distance weighted (IDW) interpolator model [100,101]. Interpolation was performed with a power parameter (p) equal to one, between 12 and 20 neighbors within a circular window. The maximum and minimum values on the interpolated surface occurs only on those sample dots with available data. For this paper, we used the Las Palmas 2.18 Q-GIS Geographic Information System, version 2.18.

2.4. Artiodactyl Index

We investigated variations in faunal consumption by measuring the variability in consumption of the best-ranked prey, the guanaco. To test if resource depression occurred in the last millennium, we use the artiodactyl index (AI), which measures the relative importance of artiodactyls in prehistoric human diets [102,103,104,105]. For each zooarchaeological assemblage, we analyzed the taxonomic structure and composition [106]. Each bone specimen was classified according to a taxonomic category and the Number of Identified Specimens (NISP) for each taxon was tallied. We estimated a set of different indices of taxonomic structure and similarity for each assemblage, such as taxonomic richness, diversity index, and evenness [45,107]. The AI considers the amount of specimens identified as guanaco, and estimates its abundance over the rest of the animal preys. It offers values between 0 and 1. An index equal to 1 indicates the dominance of the guanaco in the assemblage. Values close to 0 indicate little presence of this taxon, and as a consequence a broader diet including other taxa. This index is calculated by the formula:
AI = ΣNISP Artiodactyl/(Σ NISP Artiodactyl + NISP all other taxa)
We defined temporal analytical units for each zooarchaeological assemblage following a chrono-stratigraphic criterion. The oldest 14C date defined the lower limit, and the youngest 14C date defined the upper limit of each zooarchaeological assemblage. We considered the whole assemblage as a single temporary unit when the archaeological sites have a single date. In cases with intrusive features in their stratigraphy (i.e., Rincón del Atuel-1) [108], we considered the dated features as an independent assemblage. For sites that have more than one 14C date, which statistically overlapped (i.e., Risco de los Indios), we obtained an average value and we assigned to the assemblage a single date (Table S3) The statistic tests were performed using Past 4.01 [109]. Temporal trends in the archaeological proxies were analyzed by plotting smooth Loess curves using the GGplot2/R [109,110].

3. Results

3.1. SAM Reconstruction for the Last 1000 Years

Based on instrumental records, the SAM positive phase associates with warmer November–January temperatures over Patagonia. These anomalies are stronger on the Atlantic coast but expands to the southern Mendoza lowlands (Patagonian and Monte deserts; Figure 2A). On the other hand, the spatial patterns of correlation between SAM and precipitation from November to January show that the SAM positive phase relates to negative deviations in precipitation over the Pacific coast at 35°–38° S, reaching southern Mendoza with less intensity (Figure 2B).
The SAM reconstruction shows a period of predominant positive anomalies from the beginning of the record, 950 years cal BP, until about 550 years cal BP (SAM positive phase). After that, the negative phase persistently prevailed until the ca. 150 years cal BP, reaching extreme negative departures around 470 years cal BP (Figure 3; SAM negative phase). The persistence of the SAM positive phase implies the occurrence of warmer late spring–early summer temperatures in southern Mendoza during the first centuries of the past millennium. On the contrary, more negative SAM values prevailed between 550 and 450 years cal BP and persisted continuously, although less intense since ca. 150 years cal BP. This interval is associated with dominant below-average temperatures across the region. We assumed that the first centuries of the last millennium were relatively warm and dry, since the SAM positive phase relates to below-average rainfalls. Conversely, from the early 550 to 100 years cal BP, the climate during the late spring–summer season was predominantly colder and humid. This wetter, cooler period, which began in the region towards 550 years cal BP and persisted until approximately 100 years cal BP, was concurrent with glacial advances documented for the Andes in southern Mendoza [111] and northern Patagonia [112,113]. This cold and relatively humid period in the region, synchronous with the Little Ice Age in the North Hemisphere, contrasts with the warmer conditions recorded during the first three centuries of the last millennium. This warmer interval has been identified as the Medieval Warm Period in the Northern Hemisphere. Its occurrence in South America is still under discussion [114].

3.2. Human Populations in Northwestern Patagonia: Trends in Radiocarbon SPD

The empirical SPD for the last millennium shows variation between 0.05 and 0.11 (Figure 4). We compared the empirical SPD with those from the exponential and uniform models and found that the empirical SPD is not different from model simulations, indicating a global p = 0.6 in both cases. There are no statistically significant empirical SPD deviations from the uniform and exponential null models. The null models take the noise of the sampling and calibration processes into account. This means that while there are no changes in the empirical SPDs through time; these changes could occur due to random noise caused by sampling and calibration. By the width of the grey bar (the null models), we consider that this is product of the small sample sizes. Based on these results, we assumed that this empirical SPD represents a relatively stable population through time.
The permutation test shows that Monte and Patagonia SPDs follow the critical envelopes of the null model (Figure 5). For the interval 700 to 300 years cal BP, we noted that increases in the empirical curve for the Monte were not related to variations in the Patagonia SPD. However, after 350 years cal BP, the Monte empirical SPD declines whereas the Patagonia SPD increases (Figure 5). This opposite pattern in SPD could indicate that these deserts did not have the same human demographic response during the last 1000 years cal BP, but still the difference is not enough to be statistically significant.

3.3. Stable Isotopes (C and N) in Human Bone Collagen

We explored the stable isotope results obtained from 43 human bone samples—bone collagen and bone structural carbonate (Figure 6) (see Materials and Methods; Table S2). The northwestern Patagonia human bones have a δ13Ccoll mean of −16.0‰ (SD = 2.5; max. = −10.0‰, min. = −19.0‰). The δ13Ccoll range indicates a notable variation (9‰), direct or indirect, in C3/C4 plant consumption. For northwestern Patagonia, Gil et al. [80] proposed as the extreme to calculate diet on herbivorous bone collagen the following values: δ13C −21.4% (100% C3 in diet) and −7.5 (100% C4 in diet). Based on these values, human bone δ13Ccoll from the northwestern Patagonia mean diet composition (direct or indirect) of C4 plants is between 86% and 18% (mean 40%). The δ15N mean is 10.3 ± 1.4‰ and ranges between 7.2‰ to 13.2‰ (Figure 7; Table 1 and Table S2), showing a significant variation. Based on the mean of 5.2 ± 1.2 ‰ for northwestern Patagonia guanaco, which could be considered as the herbivorous regional trophic base line [80], the mean difference of 5.1‰ suggests a human diet focused on meat and not in vegetables. The regional human bone δ13Ccarb mean value is −11.0% ± 2.9‰, ranging between −5.1‰ and −14.9‰ (Figure 7; Table 1). It indicates a difference of 5‰ with δ13Ccoll; the correlation between δ13Ccarb and δ13Ccoll is 0.9.
Based on the chronological definition of the SAM positive (SAM+; 950 to 550 years cal BP) and negative (SAM-; 550 to ca. 150 years cal BP) phases, we used a GIS interpolation method to estimate differences in stable isotope on human bones over northwestern Patagonia. Nine human bones from Monte Desert and nine from Patagonia Desert composed the SAM positive phase samples. Human bones assigned to the SAM negative phase included nine and sixteen samples from the Monte and Patagonia deserts, respectively. Table 1 shows the basic statistics for δ13Ccoll, δ13Ccarb, and δ15N from human bones according to the desert environments and SAM phases. Figure 6 displays the relationships between δ13Ccoll and δ15N according to the deserts and SAM phases. In the Monte, two groups are recognized. The first Monte group shows δ13Ccoll between −10‰ and −15‰ and δ15N between 10‰ and 9‰. The second Monte group shows δ13Ccoll between −15.3‰ and −19% and higher δ15N between 10‰ and 13‰. A third group composed of human bones from the Patagonia Desert includes mostly (ca. 90% of the samples) individuals with lower δ13Ccoll values between −15.5‰ and −19.5‰, and highly variable δ15N between 8‰ to 13‰. Inferred from these results, the first group corresponds to individuals with higher consumption of C4 in their diets, possibly maize, than the other two groups. The last two groups have a minor, if any proportion of C4. During the SAM positive phase, the human bone samples were mostly (66%) conformed by the first Monte group. In contrast, during the SAM negative phase, the human bones include all individuals from the third group and most of the second group (Figure 6).
Figure 7 presents the relationship between human bone collagen stable isotopes with regional plants and animals stable isotopes (details in Gil et al. [97]). We grouped all data according to desert type and SAM phase. Individuals from the Patagonia Desert have higher δ15N values than plants, indicating a diet largely based on meat. In addition, most individuals (with the exception of two) in the Patagonia Desert show lower δ13C than C4 plants and cacti (Figure 7). We observe two individual groups in the Monte Desert samples: one closely related to guanaco/vizcacha and armadillo resources, and another to maize/cacti. This suggests different diets between the groups. A higher consumption of cacti and maize by the second group could be associated with higher human δ13C values. However, the presence of cacti is not significant in the archaeobotanical record nor in the assessments of the isotopic mixing model diet [97]. We suggest that maize was the primarily resource responsible for enriching these human diets.
Figure 8 compares the modelled continuous spatial distribution of δ13Ccarb and δ15N contrasting the SAM positive and negative phases. The interpolation of δ13Ccarb indicates a mean highest value during the SAM positive phase, with most of the area between −7‰ and −8.5‰. In the SAM negative phase, the individuals show values predominantly between −12‰ and −14‰. In addition, the δ15N also show contrasting spatial patterns. Lowest δ15N values, between 7‰ and 8‰, occurred during the SAM positive phase. On the contrary, highest δ15N values, between 11‰ and 13‰, were recorded during the SAM negative phase and cover most of the area. Combining both stable isotopes patterns, we infer a more plant-focused diet with stronger focus on C4 resources during the SAM positive than the negative phase.

3.4. Artiodactyl Index

The Artiodactyl Index (AI) shows no significant changes through time. The relationship between time (years cal BP) and AI is not significant (r = −0.2; p = 0.3). However, the regional pattern shows that guanaco dominated faunal assemblages after 550 years cal BP, concurrent with a dominant SAM negative phase (AI > 0.9). In contrast, AI ranges between 0.05 and 0.9 during the previous SAM positive phase (Figure 9). The Monte Desert shows lower AI values in comparison with the Patagonia Desert. It is important to highlight that all samples dating to SAM negative phase are from the Patagonia Desert. SAM positive phase zooarchaeological samples from the Monte Desert are taxonomically richer and more diverse than the Patagonia Desert samples [45]. The guanaco-dominant diet in the Patagonia Desert during the SAM negative phase reflects access to the highest-ranked resources. The observed trend in the Patagonia Desert zooarchaeological assemblages explains the highest values of δ15N in human samples from the same period noted in Figure 6 and Figure 9. At the same time, the highest variation and lowest AI during the SAM positive phase are concordant with lower δ15N and enriched δ13C, indicating a more variable and broader diet, which included a greater focus on C4 resources, such as maize, during the SAM positive phase.

4. Discussion

After 2300 years cal BP, the human population size in northwestern Patagonia changed its low growth previous trajectory and reached levels never seen before [15]. Cultigens were not used homogenously in northwestern Patagonia. In addition, probably their use also shifted across the following 2300 years cal BP [9]. Maize productivity is highly variable as a consequence of inter- and intra-annual variation in temperature and precipitation [115,116,117,118,119]. Variation in productivity affects its energetic value and order in the resources rank. Extremes temperatures and moisture have negative effects. Low temperatures during the seedling-emergence period cause crop damage by late and early frosts. Heat stress during the critical flowering period can reduce or eliminate yields. Demerit et al. [117] and Benson et al. [115] proposed the use of the corn heat unit or growing degree days (GDD) to explore this temperature constriction. Their scheme assumes that a plant matures by a daily accumulation of heat units. From this perspective, changes in temperature can affect the geographic limit and the productivity of corn by reducing the number of GDD’s accumulated per day and the length of the growing season. Either or both processes can lead to reduced yields and even total crop failure.
During SAM positive phase (from ca. 950 to ca. 550 years cal BP), lower summer precipitation but higher temperatures prevailed in northwestern Patagonia, likely increasing the GGDs. Isotopic markers of human diet indicate temporal variation in stable carbon and nitrogen isotopes (Figure 9D–F). Fitted Loess curves show a δ13Ccoll increase, ranging between −12‰ and −13‰ over the SAM positive phase (Figure 9E). The higher δ13Ccoll responds to an increase in the consumption of C4 resources, including maize (Figure 7). Both deserts, Monte and Patagonia, have similar trends but their means are different. Monte shows higher δ13C than Patagonia, which in turn presents more stable values than Monte (Figure 9D,F). Similarly, in both deserts, bone collagen δ15N averages in ca. 10‰ (Figure 9D). Since the SAM positive phase is associated with higher summer temperatures and lower precipitation, higher δ15N values should be expected in resources to this SAM positive phase [120,121,122]. Recent analyses from guanaco bone collagen show increased δ15N between 1200 and 600 years cal BP in the Monte Desert, but not in the Patagonia Desert. The increased δ15N was associated with more extreme aridity conditions in the Monte [123]. If this is the case, the herbivorous δ15N base line increased from ca. 5‰ (Figure 7) to 9‰ in the Monte Desert, but not in the Patagonia Desert (Figure 4, [123]). We assume that human diet in the Monte Desert, between 900 and 600 years cal BP, had a larger component of plants than in the Patagonia Desert, despite people occupying both deserts had similar δ15N values.
The long-term SAM negative phase, between 550 years cal BP to 100 years cal BP, is associated with increased summer precipitation and lower summer temperatures. After 550 years cal BP, the human isotopic diet changed both in δ13C and δ15N (Figure 9D–F). In the Monte Desert there is a drop in δ13Ccoll, from −13‰ to −17‰. A similar trend in δ13Ccarb is associated with increase in δ15N, from ca. 9.5‰ to ca. 12‰. The decrease in δ13C suggests a drop in C4 consumption; the increase in δ15N supports a boost in meat consumption. In the Patagonia Desert, the δ13Ccoll is stable and similar to those recorded during the SAM positive phase, from −13‰ to −17‰, with a similar trend in δ13Ccarb. In contrast, there is an increase in δ15N, from ca. 9‰ to ca. 11.5‰. According to De Fina et al. [124], maize in Mendoza shows a decrease in the production of mature grains when the mean temperature of the warmest month is lower than 19.5˚C, and maize production fails completely when summers shorten. We assume that during the SAM negative phase, summers were shorter and impacted negatively on maize production. At the same time, wetter and cooler summers could have increased primary productivity. This implies an increase of wild plants and therefore in animal biomass.
The cooler and wetter summer could have induced a wild fauna drop in δ15N due to changes in animal physiological responses to climate [120]. If so, the recorded changes in human bone collagen δ15N is contrary to our expectation based on climate. It is a change in the human diet as far as δ15N increased instead to decline as expected by the SAM negative phase. In such case, this increase in δ15N indicates a diet more focus on meat than in plants.
Variations in agricultural strategies are influenced by the effects of local environmental conditions on marginal gains, for time spent in maize production, and the abundance of high-ranked wild resources [57]. They are not influenced by fluctuations in harvest yields. We expect more agricultural investment with decreasing opportunities to exploit highest-ranked items, such as guanaco in our region [43]. We expect less time spent on farming if there is an increase in alternative economic opportunities [125]. The decline in δ13C to −18‰ indicates a reduction in maize consumption, the more ubiquitous domesticated plant in the region. This drop in maize could mean less engagement in farming.
Climatic changes associated with the SAM phases could have triggered the recorded shifts in diets. Cool summers in northwestern Patagonia during the SAM negative phase represent a major risk for farming, particularly based on maize production [119]. However, the increase in summer precipitations under the same negative phase could have facilitated the increase in primary and secondary productivity. Other things being equal, if maize farming was more expensive, and at the same time wild resources increased, we can expect a change in the diet. This change implies a retraction of the farming frontier to the north and an internal change in the population dynamics.
The regional population dynamic did not change after the drop in maize consumption; on the contrary, it remained stable throughout the last millennium. Low-level maize production was probably implemented as a supplement to a diet based on wild resources that sustained the highest population levels observed in the last millennium. At the same time, regular to low consumption of maize indicates a highly variable human diet. Furthermore, we sustain that the decline in maize production was likely associated with an increase in wild resource productivity. If this were not the case, and the wild resources did not increase, a drop in population should have been recorded, a scenario that did not happen.
Figure 4 shows a drop in the empirical northwestern Patagonia SPD after 300 years cal BP. However, this drop is not significant when compared with the uniform and empirical null models. Spaniards arrived in Mendoza in 1551, and Mendoza city was founded in 1561 [126]. Peripheral areas, such as southern Mendoza, were occupied effectively by Europeans at the beginning of the 19th century. A change in human subsistence has been interpreted as a response to a declined population caused by European diseases and pathogens [107,127,128,129]. This post 300 years cal BP drop could be significant based on historic information. If this was the case, the SAM negative phase plus a drop in human demography could impulse the same pattern: a return to a diet based on high-ranked wild resources.
Our results confirm that the northwestern Patagonia farming border was dynamic through time and space. Interactions between climate variations, ecosystem productivity, and human organization could be invoked to explain past changes in the diet of the population that occupied these territories. The increase in primary productivity triggered during the SAM negative phase could impel a significant reorganization of the human/environment system after 550 years cal BP, in which human population size was maintained despite changes in subsistence strategies and diet. The scenario post-550 years cal BP shows a shift from a system that relied more on domestic plants, and perhaps to some degree on agriculture, towards a system with emphasis on hunting and gathering, as noted during historical times by chroniclers. In an area where the environmental conditions make farming a high-cost activity—a “low rank, low yield” in terms of Winterhalder and Kennett [57]—any change in the system can impact human decision-making. Therefore, the balance between resources can shift from a system more concentrated on low-level food production to another focused on hunter-gathering and wild resources exploitation. This would imply significant flexibility in adaptation strategies that are extremely sensitive to climate changes.

Supplementary Materials

The following are available online at https://www.mdpi.com/2571-550X/3/2/17/s1, Table S1: Radiocarbon database Northwestern Patagonia (South Mendoza), Table S2: Human samples information from northwestern Patagonia with stable isotope values, references, contextual and chronological information; Table S3: Artiodactyla Index northwestern Patagonia.

Author Contributions

Conceptualization, A.F.G., R.V., C.O., F.R.F., C.C.A., E.A.P., and G.N.; methodology, R.V., G.N. and A.F.G.; software, R.V., G.N. and A.F.G.; validation, F.R.F., C.C.A., and C.O.; formal analysis, A.F.G., R.V., and C.C.A.; investigation, A.F.G., R.V., C.O., F.R.F. and E.A.P.; data curation, E.A.P. and C.C.A.; writing—original draft preparation, A.F.G., R.V., and F.R.F.; writing—review and editing, C.C.A., C.O., E.A.P. and G.N.; visualization, C.C.A. and E.A.P.; supervision, A.F.G.; project administration, A.F.G.; funding acquisition, A.F.G. and R.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by ANPCyT (grant PICT-2016/2667) and CONICET (grant PIP-11220150100342).

Acknowledgments

This research was funded by CONICET and ANPCYT. This research is part of the Past Global Changes (PAGES) PalEOclimate and the Peopling of the Earth (PEOPLE3K) network. We thank to LIECA lab (Ings. Gisela Quiroga and Armando Dauverné) for stables isotopes measurement and to Andrew Ugan for English review and comments about ideas and data presented in this paper. We acknowledge the use of imagery provided by services from NASA’s Global Imagery Browse Services (GIBS), part of NASA’s Earth Observing System Data and Information System (EOSDIS). Eugenia Ferrero (IANIGLA-CONICET) helped with preparation of Figure 3. This paper born as part of the PEOPLE3K workshop in Los Reyunos (May 2018), we thank to the participant by stimulate debates and the data exploration. A preliminary version of this paper was presented in INQUA Meeting, Dublin 2019. Thank to Erick Robinson and Jacob Freeman for exchange and comments about this topic. We thank to Museo de Historia Natural de San Rafael and Universidad Tecnológica Nacional Regional San Rafael. To editor of this especial volume, Encarni Montoya and Bronwen S. Whitney for invitation. We strongly highlight the contribution of four reviewers for their comments and criticism that helped to improve this paper. R.V. and AG were partially supported by BNP Paribas through the THEMES project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fagan, B. The Little Ice Age: How Climate Made History 1300-1850; Basic Books: Hachette, UK, 2000. [Google Scholar]
  2. Bettinger, R.L. Orderly Anarchy: Sociopolitical Evolution in Aboriginal California; University of California Press: Berkeley, CA, USA, 2015; Volume 8. [Google Scholar]
  3. Brothwell, D.R.; Cohen, M.N. The Food Crisis in Prehistory. Overpopulation and the Origins of Agriculture. Popul. Stud. 1978, 32, 593. [Google Scholar] [CrossRef]
  4. Shennan, S. The First Farmers of Europe; Cambridge University Press (CUP): Cambridge, UK, 2018. [Google Scholar]
  5. Kroeber, A.L. Cultural and Natural Areas of Native North America; University of California Press: Berkeley, CA, USA, 1939; Volume 38. [Google Scholar]
  6. Sauer, C. American Agriculture Origins: A Consideration of Nature and Culture. In Essays in Anthropology Presented to A. L. Kroeber in Celebration of His Sixtieth Birthday; University of California Press: Berkeley, CA, USA, 1936. [Google Scholar]
  7. Svizzero, S.; Tisdell, C.A. The Persistence of Hunting and Gathering Economies. Soc. Evol. Hist. 2015, 14, 3–26. [Google Scholar]
  8. Freeman, J. Alternative adaptive regimes for integrating foraging and farming activities. J. Archaeol. Sci. 2012, 39, 3008–3017. [Google Scholar] [CrossRef]
  9. Gil, A.F.; Menendez, L.P.; Atencio, J.P.; Peralta, E.A.; Neme, G.A.; Ugan, A. Estrategias Humanas, Estabilidad Y Cambio en la Frontera Agricola Sur Americana. Lat. Am. Antiq. 2017, 29, 6–26. [Google Scholar] [CrossRef] [Green Version]
  10. Simms, S.R. Ancient Peoples of the Great Basin and Colorado Plateau; Left Coast Press: Walnut Creek, CA, USA, 2016. [Google Scholar]
  11. Lagiglia, H. Actas del V Congreso Nacional de Arqueología Argentina; Universidad Nacional de San Juan: San Juan, Argentina, 1982; Volume 1, pp. 231–252. [Google Scholar]
  12. Gil, A.F. Cultígenos Prehispánicos En El Sur de Mendoza. Relac. Soc. Argent. Antropol. 1998, 22, 295–318. [Google Scholar]
  13. Johnson, A.; Gil, A.; Neme, G.; Freeman, J. Hierarchical method using ethnographic data sets to guide archaeological research: Testing models of plant intensification and maize use in Central Western Argentina. J. Anthr. Archaeol. 2015, 38, 52–58. [Google Scholar] [CrossRef]
  14. Lagiglia, H. Dinámica Cultural En El Centro Oeste y Sus Relaciones Con Áreas Aledañas Argentinas y Chilenas. In Actas del VII Congreso de Arqueología Chilena. Tomo II; Sociedad Chilena de Arqueologia: Santiago de Chile, Chile, 1977; pp. 531–560. [Google Scholar]
  15. Gil, A.F.; Giardina, M.A.; Neme, G.A.; Ugan, A. Demografía Humana e Incorporación de Cultígenos En El Centro Occidente Argentino: Explorando Tendencias En Las Fechas Radiocarbónicas; Universidad Complutense de Madrid: Madrid, Spain, 2015; Volume 44. [Google Scholar] [CrossRef]
  16. Lagiglia, H.A. Primeros Contactos Hispanico-Indígenas de Mendoza:(La Arqueología Histórica y Su Periodificación); Museo de Historia Natural de San Rafael: San Rafael, Argentina, 1983. [Google Scholar]
  17. Prieto, M. Formación y Consolidación de Una Sociedad de Frontera En Un Área Marginal Del Reino de Chile: La Provincia de Cuyo En El Siglo XVII. An. Arqueol. Etnol. 2000, 52, 17–366. [Google Scholar]
  18. Neme, G.A. Un Enfoque Regional En Cazadores-Recolectores Del Oeste Argentino: El Potencial de La Ecología Humana. In Perspectivas acutales en Arqueología Argentina; Dunken: Buenos Aires, Argentina, 2009; pp. 306–326. [Google Scholar]
  19. Bernal, V.; Gonzalez, P.N.; Gordón, F.; Perez, S.I. Exploring Dietary Patterns in the Southernmost Limit of Prehispanic Agriculture in America by Using Bayesian Stable Isotope Mixing Models. Curr. Anthr. 2016, 57, 230–239. [Google Scholar] [CrossRef] [Green Version]
  20. Corbat, M.; Zangrando, A.F.J.; Gil, A.F.; Chiavazza, H.D. Explotacion de Peces e intensificación en ambientes áridos: Comparando el registro en humedales del centro-occidente de argentina. Lat. Am. Antiq. 2017, 28, 196–212. [Google Scholar] [CrossRef]
  21. Neme, G. Cazadores-recolectores de altura en los Andes meridionales: El alto valle del río Atuel, Argentina; BAR Publishing: Oxford, UK, 2007; Volume 1591. [Google Scholar]
  22. Abraham, E.; Del Valle, H.; Roig, F.; Torres, L.M.; Ares, J.; Coronato, F.; Godagnone, R. Overview of the geography of the Monte Desert biome (Argentina). J. Arid Environ. 2009, 73, 144–153. [Google Scholar] [CrossRef]
  23. Cabrera, A.L. Regiones Fitogeográficas Argentinas; Acme: Buenos Aires, Argentina, 1976. [Google Scholar]
  24. Roig, F.; Roig-Juñent, S.; Corbalán, V. Biogeography of the Monte Desert. J. Arid Environ. 2009, 73, 164–172. [Google Scholar] [CrossRef]
  25. Otaola, C.; Ugan, A.; Gil, A. Environmental diversity and stable isotope variation in faunas: Implications for human diet reconstruction in Argentine mid-latitude deserts. J. Archaeol. Sci. Rep. 2018, 20, 57–71. [Google Scholar] [CrossRef]
  26. Labraga, J.C.; Villalba, R. Climate in the Monte Desert: Past trends, present conditions, and future projections. J. Arid Environ. 2009, 73, 154–163. [Google Scholar] [CrossRef]
  27. Marshall, G.J. Trends in the Southern Annular Mode from Observations and Reanalyses. J. Clim. 2003, 16, 4134–4143. [Google Scholar] [CrossRef]
  28. Gillet, C.; Quétin, P. Effect of temperature changes on the reproductive cycle of roach in Lake Geneva from 1983 to 2001. J. Fish Boil. 2006, 69, 518–534. [Google Scholar] [CrossRef]
  29. Earthdata. Available online: https://earthdata.nasa.gov (accessed on 15 April 2020).
  30. Perez, S.I.; Gonzalez, P.N.; Bernal, V. Past population dynamics in Northwest Patagonia: An estimation using molecular and radiocarbon data. J. Archaeol. Sci. 2016, 65, 154–160. [Google Scholar] [CrossRef]
  31. Giardina, M.; Corbat, M.; Otaola, C.; Salgán, L.; Ugan, A.; Neme, G.; Gil, A. Recursos y Dietas Humanas en Laguna Llancanelo (Mendoza, Nordpatagonia): Una Discusión Isotópica del Registro Arqueológico. Magallania (Punta Arenas) 2014, 42, 111–131. [Google Scholar] [CrossRef] [Green Version]
  32. Tripaldi, A.; Zárate, M.; Neme, G.; Gil, A.; Giardina, M.; Salgán, M.L. Archaeological site formation processes in northwestern Patagonia, Mendoza Province, Argentina. Geoarchaeology 2017, 32, 605–621. [Google Scholar] [CrossRef]
  33. Barberena, R.; Borrazzo, K.; Rughini, A.A.; Romero, G.; Pompei, M.P.; Llano, C.; De Porras, M.E.; Durán, V.; Stern, C.R.; Re, A.; et al. Perspectivas arqueológicas para Patagonia Septentrional: Sitio Cueva Huenul 1 (Provincia del Neuquén, Argentina). Magallania (Punta Arenas) 2015, 43, 137–163. [Google Scholar] [CrossRef] [Green Version]
  34. Neme, G.; Gil, A.; Otaola, C.; Giardina, M. Resource Exploitation and Human Mobility: Trends in the Archaeofaunal and Isotopic Record from Central Western Argentina. Int. J. Osteoarchaeol. 2013, 25, 866–876. [Google Scholar] [CrossRef]
  35. Gil, A.F. Arqueología de La Payunia (Mendoza, Argentina): El poblamiento humano en los márgenes de la agricultura. Payunia (Mendoza, Argentina) 2006, 1477. [Google Scholar] [CrossRef]
  36. Neme, G.; Gil, A.; Durán, V. Late Holocene in Southern Mendoza (Northwestern Patagonia) Radiocarbon Pattern and Human Occupation. Before Farming 2005, 2005, 1–18. [Google Scholar] [CrossRef]
  37. Martínez, G. Arqueología y Pobladores Antiguos de La Cuenca Del Río Colorado; Serie Aportes al Desarrollo Nacional de la Fundacion ArgenINTA: Buenos Aires, Argentina, 2015. [Google Scholar]
  38. Berón, M. Integración de Evidencias Para Evaluar Dinámica y Circulación de Poblaciones En Las Fronteras Del Río Colorado; CEQUA: Punta Arena, Chile, 2007. [Google Scholar]
  39. Franchetti, F. Hunter-Gatherer Adaptation in the Deserts of Northern Patagonia; University of Pittsburgh: Pittsburgh, PA, USA, 2019. [Google Scholar]
  40. Wolverton, S.; Otaola, C.; Neme, G.; Giardina, M.; Gil, A. Patch Choice, Landscape Ecology, and Foraging Efficiency: The Zooarchaeology of Late Holocene Foragers in Western Argentina. J. Ethnobiol. 2015, 35, 499–518. [Google Scholar] [CrossRef]
  41. Rindel, D. Explorando La Variabilidad En El Registro Zooarqueológico de La Provincia Del Neuquén: Tendencias Cronológicas y Patrones de Uso Antrópico; Gordon, F., Barberena, R., Bernal, V., Eds.; Aspha: Buenos Aires, Argentina, 2017. [Google Scholar]
  42. Otaola, C.; Giardina, M.; Corbat, M.; Fernández, F. Zooarqueología en el sur de Mendoza: Integrando Perspectivas en un Marco Biogeográfico; Neme, G., Gil, A., Eds.; Sociedad Argentina de Antropología: Buenos Aires, Argentina, 2012. [Google Scholar]
  43. Corbat, M.; Gil, A.; Zangrando, A.F. Ranking de Recursos y El Rol de Los Peces En Las Dietas Humanas Del Centro Occidente Argentino. In MC Salemme et al. (Comps.), Libro de Resúmenes del IV Congreso Nacional de Zooarqueología Argentina; Congreso Nacional de Zooarqueologia: Buenos Aires, Argentina, 2016; p. 92. Available online: https://www.researchgate.net/publication/309459807_IV_Congreso_Nacional_de_Zooarqueologia_Argentina_Libro_de_Resumenes (accessed on 15 January 2020).
  44. Neme, G.A.; Gil, A.F. Faunal exploitation and agricultural transitions in the South American agricultural limit. Int. J. Osteoarchaeol. 2008, 18, 293–306. [Google Scholar] [CrossRef]
  45. Otaola, C.; Wolverton, S.; Giardina, M.; Neme, G. Geographic scale and zooarchaeological analysis of Late Holocene foraging adaptations in western Argentina. J. Archaeol. Sci. 2015, 55, 16–25. [Google Scholar] [CrossRef]
  46. Giardina, M.A. El Aprovechamiento de La Avifauna Entre Las Sociedades Cazadoras-Recolectoras Del Sur de Mendoza: Un Enfoque Arqueozoológico; Universidad Nacional de La Plata: La Plata Argentina, Argentina, 2010. [Google Scholar]
  47. Llano, C. On optimal use of a patchy environment: Archaeobotany in the Argentinean Andes (Argentina). J. Archaeol. Sci. 2015, 54, 182–192. [Google Scholar] [CrossRef]
  48. Smith, B.D. Low-Level Food Production. J. Archaeol. Res. 2001, 9, 1–43. [Google Scholar] [CrossRef]
  49. Neme, G. El Indígeno and High-Altitude Human Occupation in the Southern Andes, Mendoza (Argentina). Lat. Am. Antiq. 2016, 27, 96–114. [Google Scholar] [CrossRef]
  50. Winterhalder, B.; Kennett, D.J. Behavioral Ecology and the Transition to Agriculture; University of California Press: Berkeley, CA, USA, 2006; Volume 1. [Google Scholar]
  51. Lema, V.S.; Della Negra, C.; Bernal, V. Explotación de recursos vegetales silvestres y domesticados en Neuquén: Implicancias del hallazgo de restos de maíz y algarrobo en artefactos de molienda del holoceno tardío. Magallania (Punta Arenas) 2012, 40, 229–247. [Google Scholar] [CrossRef] [Green Version]
  52. Llano, C.; Barberena, R. Explotación de Especies Vegetales En La Patagonia Septentrional: El Registro Arqueobotánico de Cueva Huenul 1 (Provincia de Neuquén, Argentina). Darwiniana 2013, 1, 5–19. [Google Scholar]
  53. Freeman, J.; Peeples, M.; Anderies, J.M. Toward a theory of non-linear transitions from foraging to farming. J. Anthr. Archaeol. 2015, 40, 109–122. [Google Scholar] [CrossRef]
  54. Zahid, H.J.; Robinson, E.; Kelly, R.L. Agriculture, population growth, and statistical analysis of the radiocarbon record. Proc. Natl. Acad. Sci. USA 2015, 113, 931–935. [Google Scholar] [CrossRef] [Green Version]
  55. Simms, S.R. Behavioral Ecology and Hunter-Gatherer Foraging: An example from the Great Basin; BAR Publishing: Oxford, UK, 1987; Volume 381. [Google Scholar]
  56. Bettinger, R.L. Hunter-Gatherer Foraging; Eliot Werner Publications: New York, NY, USA, 2009. [Google Scholar]
  57. Winterhalder, B.; Kennett, D.J. Seven Behavioral Ecology Reasons for the Persistence of Foragers with Cultivars. In Cowboy Ecologist: Essays in Honor of Robert L. Bettinger; Delacorte, M., Terry, J.L., Eds.; CARD (Center for Archaeological Research at Davis): Davis, CA, USA, 2020; pp. 93–110. [Google Scholar]
  58. Lupo, K. Evolutionary Foraging Models in Zooarchaeological Analysis: Recent Applications and Future Challenges. J. Archaeol. Res. 2007, 15, 143–189. [Google Scholar] [CrossRef]
  59. Dätwyler, C.; Neukom, R.; Abram, N.J.; Gallant, A.J.E.; Grosjean, M.; Jacques-Coper, M.; Karoly, D.J.; Villalba, R. Teleconnection stationarity, variability and trends of the Southern Annular Mode (SAM) during the last millennium. Clim. Dyn. 2017, 51, 2321–2339. [Google Scholar] [CrossRef]
  60. Bianchi, E.; Villalba, R.; Solarte, A. NDVI Spatio-temporal Patterns and Climatic Controls over Northern Patagonia. Ecosystem 2019, 23, 84–97. [Google Scholar] [CrossRef]
  61. Crema, E.R.; Bevan, A.; Shennan, S. Spatio-temporal approaches to archaeological radiocarbon dates. J. Archaeol. Sci. 2017, 87, 1–9. [Google Scholar] [CrossRef] [Green Version]
  62. Robinson, E.; Nicholson, C.; Kelly, R.L. The Importance of Spatial Data to Open-Access National Archaeological Databases and the Development of Paleodemography Research. Adv. Archaeol. Pr. 2019, 7, 395–408. [Google Scholar] [CrossRef] [Green Version]
  63. Rick, J.W. Dates as Data: An Examination of the Peruvian Preceramic Radiocarbon Record. Am. Antiq. 1987, 52, 55. [Google Scholar] [CrossRef]
  64. Bevan, A.; Colledge, S.; Fuller, D.Q.; Fyfe, R.; Shennan, S.; Stevens, C.J. Holocene fluctuations in human population demonstrate repeated links to food production and climate. Proc. Natl. Acad. Sci. USA 2017, 114, E10524–E10531. [Google Scholar] [CrossRef] [Green Version]
  65. Freeman, J.; A Byers, D.; Robinson, E.; Kelly, R.L. Culture Process and the Interpretation of Radiocarbon Data. Radiocarbon 2017, 60, 453–467. [Google Scholar] [CrossRef] [Green Version]
  66. Shennan, S.; Downey, S.S.; Timpson, A.; Edinborough, K.; Colledge, S.; Kerig, T.; Manning, K.; Thomas, M.G. Regional population collapse followed initial agriculture booms in mid-Holocene Europe. Nat. Commun. 2013, 4, 2486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Surovell, T.A.; Brantingham, P.J. A note on the use of temporal frequency distributions in studies of prehistoric demography. J. Archaeol. Sci. 2007, 34, 1868–1877. [Google Scholar] [CrossRef]
  68. Timpson, A.; Colledge, S.; Crema, E.R.; Edinborough, K.; Kerig, T.; Manning, K.; Thomas, M.G.; Shennan, S. Reconstructing regional population fluctuations in the European Neolithic using radiocarbon dates: A new case-study using an improved method. J. Archaeol. Sci. 2014, 52, 549–557. [Google Scholar] [CrossRef] [Green Version]
  69. Gordón, F.; Béguelin, M.; Rindel, D.; Della Negra, C.; Hajduk, A.; Vazquez, R.; Cobos, V.; Perez, S.I.; Bernal, V. Estructura Espacial y Dinámica Temporal de La Ocupación Humana de Neuquén (Patagonia Argentina) Durante El Pleistoceno Final-Holoceno. Intersecc. Antropol. 2019, 20, 93–105. [Google Scholar]
  70. Kelly, R.L.; Surovell, T.A.; Shuman, B.N.; Smith, G.M. A continuous climatic impact on Holocene human population in the Rocky Mountains. Proc. Natl. Acad. Sci. USA 2012, 110, 443–447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Peros, M.C.; Munoz, S.E.; Gajewski, K.; Viau, A. Prehistoric demography of North America inferred from radiocarbon data. J. Archaeol. Sci. 2010, 37, 656–664. [Google Scholar] [CrossRef]
  72. Riede, F. Climate and Demography in Early Prehistory: Using Calibrated 14 C Dates as Population Proxies. Hum. Boil. 2009, 81, 309–337. [Google Scholar] [CrossRef] [Green Version]
  73. Surovell, T.A.; Finley, J.B.; Smith, G.M.; Brantingham, P.J.; Kelly, R.L. Correcting temporal frequency distributions for taphonomic bias. J. Archaeol. Sci. 2009, 36, 1715–1724. [Google Scholar] [CrossRef]
  74. Bevan, A.; Crema, E.R. Rcarbon v1. 2.0: Methods for Calibrating and Analy. Available online: https://CRAN.R-project.org/package=rcarbon (accessed on 31 March 2020).
  75. Manning, K.; Timpson, A. The demographic response to Holocene climate change in the Sahara. Quat. Sci. Rev. 2014, 101, 28–35. [Google Scholar] [CrossRef] [Green Version]
  76. Crema, E.R.; Habu, J.; Kobayashi, K.; Madella, M. Summed Probability Distribution of 14C Dates Suggests Regional Divergences in the Population Dynamics of the Jomon Period in Eastern Japan. PLoS ONE 2016, 11, e0154809. [Google Scholar] [CrossRef] [Green Version]
  77. Durán, V. Nuevas Consideraciones Sobre La Problemática Arqueológica Del Río Grande (Malargüe, Mendoza); Gil, A.F., Neme, G.A., Eds.; Sociedad Argentina de Antropología: Buenos Aires, Argentina, 2002. [Google Scholar]
  78. Roberts, N.; Woodbridge, J.; Bevan, A.; Palmisano, A.; Shennan, S.; Asouti, E. Human responses and non-responses to climatic variations during the last Glacial-Interglacial transition in the eastern Mediterranean. Quat. Sci. Rev. 2018, 184, 47–67. [Google Scholar] [CrossRef]
  79. Katzenberg, M.A. Stable Isotope Analysis: A Tool for Studying Past Diet, Demography, and Life History. In Biological Anthropology of the Human Skeleton; Wiley: Hoboken, NJ, USA, 2008; pp. 411–441. [Google Scholar]
  80. Gil, A.F.; Ugan, A.; Neme, G.A. More Carnivorous than Vegetarian: Isotopic Perspectives on Human Diets in Late Holocene Northwestern Patagonia. J. Archaeol. Sci. Rep 2019. in review. [Google Scholar]
  81. Gil, A.; Neme, G.; Tykot, R.H. Stable isotopes and human diet in central western Argentina. J. Archaeol. Sci. 2011, 38, 1395–1404. [Google Scholar] [CrossRef]
  82. Coltrain, J.B.; Janetski, J.C.; Carlyle, S.W. The Stable and Radio-Isotope Chemistry of Eastern Basketmaker and Pueblo Groups in the Four Corners Region of the American Southwest. In Histories of Maize; Elsevier BV: Amsterdam, The Netherlands, 2006; pp. 275–287. [Google Scholar]
  83. Ambrose, S.H. Preparation and characterization of bone and tooth collagen for isotopic analysis. J. Archaeol. Sci. 1990, 17, 431–451. [Google Scholar] [CrossRef]
  84. Hedges, R.E.M. On bone collagen?apatite-carbonate isotopic relationships. Int. J. Osteoarchaeol. 2003, 13, 66–79. [Google Scholar] [CrossRef]
  85. Coltrain, J.B.; Leavitt, S.W. Climate and Diet in Fremont Prehistory: Economic Variability and Abandonment of Maize Agriculture in the Great Salt Lake Basin. Am. Antiq. 2002, 67, 453–485. [Google Scholar] [CrossRef]
  86. Ambrose, S.H.; Norr, L. Experimental Evidence for the Relationship of the Carbon Isotope Ratios of Whole Diet and Dietary Protein to Those of Bone Collagen and Carbonate; Springer Science and Business Media LLC: Berlin, Germany, 1993; pp. 1–37. [Google Scholar]
  87. Tieszen, L.L.; Fagre, T. Effect of Diet Quality and Composition on the Isotopic Composition of Respiratory CO2, Bone Collagen, Bioapatite, and Soft Tissues; Springer Science and Business Media LLC: Berlin, Germany, 1993; pp. 121–155. [Google Scholar]
  88. Fernandes, R.; Nadeau, M.-J.; Grootes, P.M. Macronutrient-based model for dietary carbon routing in bone collagen and bioapatite. Archaeol. Anthr. Sci. 2012, 4, 291–301. [Google Scholar] [CrossRef]
  89. Ambrose, S.H. Controlled Diet and Climate Experiments on Nitrogen Isotope Ratios of Rats. In Advances in Archaeological and Museum Science; Springer Science and Business Media LLC: Berlin, Germany, 2006; Volume 5, pp. 243–259. [Google Scholar]
  90. Hedges, R.E.; Reynard, L.M. Nitrogen isotopes and the trophic level of humans in archaeology. J. Archaeol. Sci. 2007, 34, 1240–1251. [Google Scholar] [CrossRef]
  91. Bocherens, H.; Drucker, D.G. Trophic level isotopic enrichment of carbon and nitrogen in bone collagen: Case studies from recent and ancient terrestrial ecosystems. Int. J. Osteoarchaeol. 2003, 13, 46–53. [Google Scholar] [CrossRef]
  92. Tykot, R.H.; Falabella, F.; Planella, M.T.; Aspillaga, E.; Sanhueza, L.; Becker, C. Stable isotopes and archaeology in central Chile: Methodological insights and interpretative problems for dietary reconstruction. Int. J. Osteoarchaeol. 2009, 19, 156–170. [Google Scholar] [CrossRef]
  93. Froehle, A.; Kellner, C.; Schoeninger, M. FOCUS: Effect of diet and protein source on carbon stable isotope ratios in collagen: Follow up to. J. Archaeol. Sci. 2010, 37, 2662–2670. [Google Scholar] [CrossRef]
  94. Keeling, C.D. The Suess effect: 13Carbon-14Carbon interrelations. Environ. Int. 1979, 2, 229–300. [Google Scholar] [CrossRef]
  95. Athfield, N.B.; Green, R.C.; Craig, J.; McFadgen, B.; Bickler, S. Influence of marine sources on 14 C ages: Isotopic data from Watom Island, Papua New Guinea inhumations and pig teeth in light of new dietary standards. J. R. Soc. N. Z. 2008, 38, 1–23. [Google Scholar] [CrossRef] [Green Version]
  96. Telford, R.J.; Heegaard, E.; Birks, H.J.B. The intercept is a poor estimate of a calibrated radiocarbon age. Holocene 2004, 14, 296–298. [Google Scholar] [CrossRef] [Green Version]
  97. Gil, A.; Ugan, A.; Otaola, C.; Neme, G.; Giardina, M.; Menendez, L.P. Variation in camelid δ13C and δ15N values in relation to geography and climate: Holocene patterns and archaeological implications in central western Argentina. J. Archaeol. Sci. 2016, 66, 7–20. [Google Scholar] [CrossRef]
  98. Ascough, P.; Mainland, I.; Newton, A. From Isoscapes to Farmscapes: Introduction to the Special Issue. Environ. Archaeol. 2018, 23, 1–4. [Google Scholar] [CrossRef] [Green Version]
  99. Bowen, G.J. Isoscapes: Spatial Pattern in Isotopic Biogeochemistry. Annu. Rev. Earth Planet. Sci. 2010, 38, 161–187. [Google Scholar] [CrossRef] [Green Version]
  100. Wang, S.; Fan, J.; Zhong, H.; Li, Y.; Zhu, H.; Qiao, Y.; Zhang, H. A multi-factor weighted regression approach for estimating the spatial distribution of soil organic carbon in grasslands. Catena 2019, 174, 248–258. [Google Scholar] [CrossRef]
  101. Barrientos, G.; Catella, L.; Morales, N.S. A journey into the landscape of past feeding habits: Mapping geographic variations in the isotope (δ15N) -inferred trophic position of prehistoric human populations. Quat. Int. 2020. [Google Scholar] [CrossRef]
  102. Janetski, J.C. Fremont Hunting and Resource Intensification in the Eastern Great Basin. J. Archaeol. Sci. 1997, 24, 1075–1088. [Google Scholar] [CrossRef] [Green Version]
  103. Bayham, F.E. Factors Influencing the Archaic Pattern of Animal Exploitation. KIVA 1979, 44, 219–235. [Google Scholar] [CrossRef]
  104. Broughton, J.M. Late Holocene Resource Intensification in the Sacramento Valley, California: The Vertebrate Evidence. J. Archaeol. Sci. 1994, 21, 501–514. [Google Scholar] [CrossRef]
  105. Broughton, J.M. Declines in Mammalian Foraging Efficiency during the Late Holocene, San Francisco Bay, California. J. Anthr. Archaeol. 1994, 13, 371–401. [Google Scholar] [CrossRef]
  106. Lyman, R.L. Quantitative Paleozoology by R. Lee Lyman; Cambridge University Press (CUP): Cambridge, UK, 2008. [Google Scholar]
  107. Neme, G.A.; Gil, A.F.; Abbona, C.C.; Otaola, C.; Nagaoka, L.; Giardina, M.; Wolverton, S.; Johnson, J. Exploring the Environmental Rebound in Northwest Patagonia: Zooarchaeological, Stable Isotope, Radiocarbon Frequency and DNA Trends; Book Unpublished; 2020. [Google Scholar]
  108. Dieguez, S.; Gil, A.; Neme, G.; Zárate, M.; De Francesco, C.; Strasser, E. Cronoestratigrafía Del Sitio Rincón Del Atuel-1 (San Rafael, Mendoza): Formación Del Sitio y Ocupación Humana. Intersecc. Antropol. 2004, 5, 71–80. [Google Scholar]
  109. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontol. Electron. 2001, 4, 9. [Google Scholar]
  110. Baxter, M.J. Statistics in Archaeology; Arnold: London, UK, 2003. [Google Scholar]
  111. Espizua, L.E.; Pitte, P. The Little Ice Age glacier advance in the Central Andes (35° S), Argentina. Palaeogeogr. Palaeoclim. Palaeoecol. 2009, 281, 345–350. [Google Scholar] [CrossRef]
  112. Villalba, R. Climatic Fluctuations in Northern Patagonia during the Last 1000 Years as Inferred from Tree-Ring Records. Quat. Res. 1990, 34, 346–360. [Google Scholar] [CrossRef]
  113. Masiokas, M.H.; Rivera, A.; Espizua, L.E.; Villalba, R.; Delgado, S.; Aravena, J.C. Glacier fluctuations in extratropical South America during the past 1000 years. Palaeogeogr. Palaeoclim. Palaeoecol. 2009, 281, 242–268. [Google Scholar] [CrossRef]
  114. IPCC. The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2007; Volume 996, p. 2007. [Google Scholar]
  115. Benson, L. Factors Controlling Pre-Columbian and Early Historic Maize Productivity in the American Southwest, Part 1: The Southern Colorado Plateau and Rio Grande Regions. J. Archaeol. Method Theory 2010, 18, 1–60. [Google Scholar] [CrossRef] [Green Version]
  116. Claessens, L.; Antle, J.M.; Stoorvogel, J.J.; Valdivia, R.; Thornton, P.; Herrero, M. A method for evaluating climate change adaptation strategies for small-scale farmers using survey, experimental and modeled data. Agric. Syst. 2012, 111, 85–95. [Google Scholar] [CrossRef] [Green Version]
  117. Demeritt, D. Agriculture, Climate, and Cultural Adaptation in the Prehistoric Northeast. Archaeol. East. N. Am. 1991, 19, 183–202. [Google Scholar]
  118. Maddonni, G. Ángel Analysis of the climatic constraints to maize production in the current agricultural region of Argentina—A probabilistic approach. Theor. Appl. Clim. 2011, 107, 325–345. [Google Scholar] [CrossRef]
  119. Gil, A.; Villalba, R.; Ugan, A.; Cortegoso, V.; Neme, G.; Michieli, C.T.; Novellino, P.; Durán, V. Isotopic evidence on human bone for declining maize consumption during the little ice age in central western Argentina. J. Archaeol. Sci. 2014, 49, 213–227. [Google Scholar] [CrossRef]
  120. Heaton, T.H.E.; Vogel, J.C.; Von La Chevallerie, G.; Collett, G. Climatic influence on the isotopic composition of bone nitrogen. Nature 1986, 322, 822–823. [Google Scholar] [CrossRef]
  121. Sealy, J.; Van Der Merwe, N.J.; Thorp, J.A.; Lanham, J.L. Nitrogen isotopic ecology in southern Africa: Implications for environmental and dietary tracing. Geochim. Cosmochim. Acta 1987, 51, 2707–2717. [Google Scholar] [CrossRef]
  122. Ambrose, S.H.; Deniro, M.J. Climate and Habitat Reconstruction Using Stable Carbon and Nitrogen Isotope Ratios of Collagen in Prehistoric Herbivore Teeth from Kenya. Quat. Res. 1989, 31, 407–422. [Google Scholar] [CrossRef]
  123. Gil, A.; Otaola, C.; Neme, G.; Peralta, E.A.; Abbona, C.; Quiroga, G.; Dauverné, A.; Seitz, V.P. Lama guanicoe bone collagen stable isotope (C and N) indicate climatic and ecological variation during Holocene in Northwest Patagonia. Quat. Int. 2019. [Google Scholar] [CrossRef]
  124. De Fina, A.L.; Giannetto, F.; Richard, A.E.; Sabella, L. Difusión Geográfica de Cultivos Índices En La Provincia de Mendoza y Sus Causas. INTA Inst. Suelos Agrotec. Mendoza 1964, 83, 2–87. [Google Scholar]
  125. Barlow, K.R. Predicting Maize Agriculture among the Fremont: An Economic Comparison of Farming and Foraging in the American Southwest. Am. Antiq. 2002, 67, 65–88. [Google Scholar] [CrossRef]
  126. Michieli, C.T. Antigua Historia de Cuyo (Siglos XVI a XVIII); Ansilta: San Juan, Argentina, 1994. [Google Scholar]
  127. Denevan, W.M. The Pristine Myth: The Landscape of the Americas in 1492. Ann. Assoc. Am. Geogr. 1992, 82, 369–385. [Google Scholar] [CrossRef]
  128. Denevan, W.M. After 1492: Nature Rebounds. Geogr. Rev. 2016, 106, 381–398. [Google Scholar] [CrossRef]
  129. O’Fallon, B.D.; Fehren-Schmitz, L. Native Americans experienced a strong population bottleneck coincident with European contact. Proc. Natl. Acad. Sci. USA 2011, 108, 20444–20448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Study area with archaeological sites mentioned in the text. Archaeological sites references: 1: Laguna El Diamante; 2: Risco de Los Indios; 3: La Hedionda; 4: El Indígeno; 5: Gruta El Carrizalito; 6: Dique 25 de Mayo; 7: Los Peuquenes; 8: Los Coroneles; 9: Puesto Agua del Zapallo; 10: Los Reyunos; 11: Alero Montiel; 12: Loma del Eje; 13: El Durazno; 14: Rincón del Atuel 1; 15: Gruta del Indio; 16: MDA.114.116; 17: Cerro Negro; 18: Real del Padre 1; 19: Agua Buena; 20: Arroyo Malo 1; 21: La Olla; 22: El Bosquecillo 5; 23: El Bosquecillo 3; 24: Respolar; 25: El Nihuil; 26: El Sosneado; 27: Tierras Blancas; 28: Arbolito; 29: El Nihuil 1; 30: Las Ramada; 31: Arbolito-6; 32: Ojo de Agua; 33: Puesto Jaque 2; 34: Gendarmería Nacional 5; 35: Cueva A. Colorado; 36: Los Leones 5; 37: Agua de los Caballos; 38: Zanjón El Morado; 39: Agua de la Mula; 40: Cerro Mesa; 41: Cueva Zanjón del Buitre; 42: Arroyo Panchino; 43: Puesto Ortubia; 44: LLancanelo 22; 45: Norte de Llancanelo; 46: A Negro de Pincheira; 47: Bajada de Las Tropas; 48: Bajada de las Tropas; 49: Ea. Las Chacras; 50: Llancanelo; 51: Los Menucos; 52: Llancanelo 17; 53: Cochico; 54: Llancanelo 50; 55: Pozos de Carapacho; 56: Manqui Malal; 57: Cañada de las Vizcachas; 58: Casa de Piedra El Chequenco; 59: Ciénaga de Borbarán; 60: El Manzano; 61: Alero Puesto Carrasco; 62: A. Mechenquil; 63: El Alambrado; 64: El Payén; 65: La Peligrosa 2; 66: La Matancilla; 67: Agua del Toro; 68: Agua de Perez 8B; 69: Agua de Perez 1; 70: Médano R6; 71: Corcovo 1; 72: Cueva Palulo; 73: Campos Las Julias; 74: Carmonina; 75: Gruta Puesto Las Tinajas; 76: Los Gallegos 1; 77: Cueva Salamanca; 78: Cupertino; 79: Gruta del Manzano; 80: Cañón del Atuel; 81:Arroyo de los Jilgueros; 82: Camping Familiar Cristiano; 83: Campos Las Julias; 84: La Hedionda; 85: Puesto Aisol; 86: Laguna El Sosneado. Map base [29].
Figure 1. Study area with archaeological sites mentioned in the text. Archaeological sites references: 1: Laguna El Diamante; 2: Risco de Los Indios; 3: La Hedionda; 4: El Indígeno; 5: Gruta El Carrizalito; 6: Dique 25 de Mayo; 7: Los Peuquenes; 8: Los Coroneles; 9: Puesto Agua del Zapallo; 10: Los Reyunos; 11: Alero Montiel; 12: Loma del Eje; 13: El Durazno; 14: Rincón del Atuel 1; 15: Gruta del Indio; 16: MDA.114.116; 17: Cerro Negro; 18: Real del Padre 1; 19: Agua Buena; 20: Arroyo Malo 1; 21: La Olla; 22: El Bosquecillo 5; 23: El Bosquecillo 3; 24: Respolar; 25: El Nihuil; 26: El Sosneado; 27: Tierras Blancas; 28: Arbolito; 29: El Nihuil 1; 30: Las Ramada; 31: Arbolito-6; 32: Ojo de Agua; 33: Puesto Jaque 2; 34: Gendarmería Nacional 5; 35: Cueva A. Colorado; 36: Los Leones 5; 37: Agua de los Caballos; 38: Zanjón El Morado; 39: Agua de la Mula; 40: Cerro Mesa; 41: Cueva Zanjón del Buitre; 42: Arroyo Panchino; 43: Puesto Ortubia; 44: LLancanelo 22; 45: Norte de Llancanelo; 46: A Negro de Pincheira; 47: Bajada de Las Tropas; 48: Bajada de las Tropas; 49: Ea. Las Chacras; 50: Llancanelo; 51: Los Menucos; 52: Llancanelo 17; 53: Cochico; 54: Llancanelo 50; 55: Pozos de Carapacho; 56: Manqui Malal; 57: Cañada de las Vizcachas; 58: Casa de Piedra El Chequenco; 59: Ciénaga de Borbarán; 60: El Manzano; 61: Alero Puesto Carrasco; 62: A. Mechenquil; 63: El Alambrado; 64: El Payén; 65: La Peligrosa 2; 66: La Matancilla; 67: Agua del Toro; 68: Agua de Perez 8B; 69: Agua de Perez 1; 70: Médano R6; 71: Corcovo 1; 72: Cueva Palulo; 73: Campos Las Julias; 74: Carmonina; 75: Gruta Puesto Las Tinajas; 76: Los Gallegos 1; 77: Cueva Salamanca; 78: Cupertino; 79: Gruta del Manzano; 80: Cañón del Atuel; 81:Arroyo de los Jilgueros; 82: Camping Familiar Cristiano; 83: Campos Las Julias; 84: La Hedionda; 85: Puesto Aisol; 86: Laguna El Sosneado. Map base [29].
Quaternary 03 00017 g001
Figure 2. (A) Spatial correlations between the Southern Annular Mode (SAM) and seasonal (November–January) temperatures from CRU version TS4.01 (1957–2006). (B) Spatial correlations between SAM and seasonal (November–January) precipitation from CRU version TS4.01 (1957–2006).
Figure 2. (A) Spatial correlations between the Southern Annular Mode (SAM) and seasonal (November–January) temperatures from CRU version TS4.01 (1957–2006). (B) Spatial correlations between SAM and seasonal (November–January) precipitation from CRU version TS4.01 (1957–2006).
Quaternary 03 00017 g002
Figure 3. Reconstruction based on annual-resolution proxy records of the summer (DJF) Southern Annular Mode (SAM) or Antarctic Oscillation over the past 1000 years. The observed SAM record is shown in light blue. The reconstruction (thin blue line) is also shown in a smoothed version (violet line) using a digital filter (64-year spline). During the interval 950-550 years cal BP, associated with the Medieval Warm Period, the SAM oscillated around its long-term mean value. On the contrary, its values were mostly negative between 550 and 100 years cal BP (particularly 550-450 years cal BP), concurrent with the Little Ice Age.
Figure 3. Reconstruction based on annual-resolution proxy records of the summer (DJF) Southern Annular Mode (SAM) or Antarctic Oscillation over the past 1000 years. The observed SAM record is shown in light blue. The reconstruction (thin blue line) is also shown in a smoothed version (violet line) using a digital filter (64-year spline). During the interval 950-550 years cal BP, associated with the Medieval Warm Period, the SAM oscillated around its long-term mean value. On the contrary, its values were mostly negative between 550 and 100 years cal BP (particularly 550-450 years cal BP), concurrent with the Little Ice Age.
Quaternary 03 00017 g003
Figure 4. Summed probability distributions (SPDs) for northwestern Patagonia (southern Mendoza) and their fitted hypothetical null models (exponential and uniform). The black lines represent the empirical SPD while the grey areas are the 95% confidence interval for the simulated null models.
Figure 4. Summed probability distributions (SPDs) for northwestern Patagonia (southern Mendoza) and their fitted hypothetical null models (exponential and uniform). The black lines represent the empirical SPD while the grey areas are the 95% confidence interval for the simulated null models.
Quaternary 03 00017 g004
Figure 5. Radiocarbon SPDs for the Patagonia and Monte deserts over the past millennium.
Figure 5. Radiocarbon SPDs for the Patagonia and Monte deserts over the past millennium.
Quaternary 03 00017 g005
Figure 6. Individual variations in human bone collagen δ13Ccoll15N according to desert (Monte/Patagonia) and SAM phase (positive/negative).
Figure 6. Individual variations in human bone collagen δ13Ccoll15N according to desert (Monte/Patagonia) and SAM phase (positive/negative).
Quaternary 03 00017 g006
Figure 7. Regional resource trends in δ13C and δ15N for the Monte and Patagonia deserts (data base from Gil et al. [80]) and their relationships with human bone collagen isotopes (Table S2).
Figure 7. Regional resource trends in δ13C and δ15N for the Monte and Patagonia deserts (data base from Gil et al. [80]) and their relationships with human bone collagen isotopes (Table S2).
Quaternary 03 00017 g007
Figure 8. Northwestern Patagonia–southern Mendoza isoscapes (IDW interpolation) segregated by SAM phase (see details in the text).
Figure 8. Northwestern Patagonia–southern Mendoza isoscapes (IDW interpolation) segregated by SAM phase (see details in the text).
Quaternary 03 00017 g008
Figure 9. Multiproxy temporal analysis over the past millennium. (A) SAM variations from Figure 3; (B) radiocarbon SPD based on Figure 4; (C) Temporal variations for the Artiodactyl Index (AI); (D) long-term δ15N variations using the LOESS filter (individual data segregated by desert); (E) δ13Ccoll temporal trend using the LOESS filter (individual data segregated by desert); (F) δ13Ccarb temporal trend using the LOESS filter (individual data segregated by desert).
Figure 9. Multiproxy temporal analysis over the past millennium. (A) SAM variations from Figure 3; (B) radiocarbon SPD based on Figure 4; (C) Temporal variations for the Artiodactyl Index (AI); (D) long-term δ15N variations using the LOESS filter (individual data segregated by desert); (E) δ13Ccoll temporal trend using the LOESS filter (individual data segregated by desert); (F) δ13Ccarb temporal trend using the LOESS filter (individual data segregated by desert).
Quaternary 03 00017 g009
Table 1. Basic statistic in δ13Ccoll, δ13Ccarb, and δ15N from human bones in southern Mendoza. Results are arranged by deserts (Monte and Patagonia) and SAM phases (positive and negative; see details in main text).
Table 1. Basic statistic in δ13Ccoll, δ13Ccarb, and δ15N from human bones in southern Mendoza. Results are arranged by deserts (Monte and Patagonia) and SAM phases (positive and negative; see details in main text).
Descriptive StatisticsSAM Positive PhaseSAM Negative Phase
MONTEPATAGONIAMONTEPATAGONIA
δ13Ccollδ13Ccarbδ15Nδ13Ccolδ13Ccarbδ15Nδ13Ccolδ13Ccarbδ15Nδ13Ccolδ13Ccarbδ15N
N999999999161616
Min−7.1−12.39.6−19.4−14.57.2−8.9−13.58.99.4−14.97.4
Max−12.1−6.512.0−13.7−7.611.7−10.1−5.112.9−14.7−8.213.2
Mean−13.5−8.410.3−17.6−12.99.5−14.4−8.810.2−17.6−12.810.9
SD1.71.80.91.82.11.72.62.61.41.41.91.4
Median−12.9−8.010.2−18.3−13.68.9−14.5−7.99.9−17.9−13.110.9

Share and Cite

MDPI and ACS Style

Gil, A.F.; Villalba, R.; Franchetti, F.R.; Otaola, C.; Abbona, C.C.; Peralta, E.A.; Neme, G. Between Foragers and Farmers: Climate Change and Human Strategies in Northwestern Patagonia. Quaternary 2020, 3, 17. https://doi.org/10.3390/quat3020017

AMA Style

Gil AF, Villalba R, Franchetti FR, Otaola C, Abbona CC, Peralta EA, Neme G. Between Foragers and Farmers: Climate Change and Human Strategies in Northwestern Patagonia. Quaternary. 2020; 3(2):17. https://doi.org/10.3390/quat3020017

Chicago/Turabian Style

Gil, Adolfo F., Ricardo Villalba, Fernando R. Franchetti, Clara Otaola, Cinthia C. Abbona, Eva A. Peralta, and Gustavo Neme. 2020. "Between Foragers and Farmers: Climate Change and Human Strategies in Northwestern Patagonia" Quaternary 3, no. 2: 17. https://doi.org/10.3390/quat3020017

APA Style

Gil, A. F., Villalba, R., Franchetti, F. R., Otaola, C., Abbona, C. C., Peralta, E. A., & Neme, G. (2020). Between Foragers and Farmers: Climate Change and Human Strategies in Northwestern Patagonia. Quaternary, 3(2), 17. https://doi.org/10.3390/quat3020017

Article Metrics

Back to TopTop