Agricultural Frontiers, Health Care, and Population Size Impact the Recovery Patterns of Brazilian Indigenous Nations

After centuries of decline, the populations of indigenous nations in Brazil began to increase in the 1970s. Population Ecology theory predicts that population size affects the dynamics of three basic recovery patterns: intra-specific cooperation (a positive effect of population size); competition (a negative effect); and exponential growth (no effect of population size). We analyzed the dynamics and recent history of 25 Brazilian indigenous populations using a cross-sectional approach to understand how exogenous and cultural variables (e.g., deforestation, diet richness) interact with population levels. We found that population size, access to health care, the extent of indigenous territories, and degree of deforestation impact the recovery of indigenous population levels.


Introduction
American indigenous populations suffered drastic declines and widespread acculturation following the arrival and expanding settlement of Europeans from the late fifteenth century onward (Ribeiro 1956). Conservative estimates of indigenous pre-contact Brazilian populations report around three million people with 1000 distinct languages (Melatti 1993). After 500 years of contact, the indigenous populations are now estimated at only slightly more than 800,000 individuals with approximately 180 languages persisting (IBGE 2010; S1 Text).
In the last century, Latin America experienced a political shift towards the protection of indigenous rights (McSweeney and Shahna 2005) notably in Brazil during the 1950s and 1960s, peaking with the creation of the Parque Indígena do Xingu (PIX) in 1961. After 22 years of dictatorship, the Brazilian Constitution of 1988 explicitly guaranteed the rights of indigenous nations to a territory and access to basic health care. Since then, Indigenous Territories (ITs) have expanded further and currently occupy 13% of the country (Instituto Socioambiental 2011). Indigenous communities began to receive better health assistance and government subsidies that reduced mortality rates and allowed demographic recovery rates with an average of approximately 3% per year (estimated from at least 50 indigenous nations), a value much higher than the national rates (Melatti 1993;Ricardo 2000;Pagliaro et al. 2005;Instituto Socioambiental 2011;IBGE 2014). This growth has led to many nations regaining pre-contact population levels, allowing them to maintain their social and cultural integrity (Pozzobon 1994;S2-S4 Text).
Nonetheless, Brazilian indigenous nations continue to endure higher disease burdens and mortality levels as well as lower life expectancies than non-indigenous populations (Junior 1985(Junior , 1989Junior and Santos 1991;Fagundes et al. 2004;Orellana et al. 2006). ITs are subject to habitat destruction ( Fig. 1) and land conflicts arising from illegal mining and logging activities and expansion of agricultural frontiers, particularly for ranching and large-scale cultivation enterprises 1 (Villas Bôas and Villas Bôas 1994;Flowers 1994;Junior et al. 2002;Baruzzi 2005;Fearnside 2005;Paula 2005; Barbosa et al. 2008;Filho and de Souza 2009;Lemos and Silva 2011). These developments serve to undermine any benefits to indigenous communities brought by government designation of their territories and improved access to health care (Ribeiro 1970; S5 and S6 Text; Fig. 1).
Population Ecology theory adopts principles from Dynamic Systems theory, which identifies both endogenous and exogenous processes that could affect the dynamics of any variable (Royama 1992;Berryman 1999;Turchin 2001;Hopfenberg 2003). Endogenous processes occur when a variable affects directly or indirectly (through another variable) its own rate of change and dynamics, as in the predator-prey interaction. An exogenous process occurs when a variable affects another, but the effect does not feed back, e.g., climate variables, natural disasters, the introduction of invasive species, and wars could affect population dynamics (Lima and Berryman 2011;Lima 2014;Lima et al. 2016).
Population Ecology Theory dictates that under a favorable and constant environment any population may grow exponentially (Ginzburg 1986;Royama 1992;Berryman 1999). In this case, the recovery rate is constant and independent of population size. There are two basic deviations from exponential growth: (1) when population size enhances environmental suitability such as through cooperative activities (intraspecific cooperation) in a stressful environment (Allee, 1932;Berryman 1999;Luděk et al. 2006); and (2) when population size has a negative effect on environmental suitability (intra-specific competition) such as by reducing per capita available resources (Malthus 1798;Verhulst 1838;Allee, 1932;Royama 1992;Berryman 1999;Turchin 2001;Hopfenberg 2003;Luděk et al. 2006). In summary, adding individuals to a population may either enhance or deplete environmental suitability (Malthus 1798;Verhulst 1838;Royama 1992;Berryman 1999;Turchin 2001;Hopfenberg 2003). Moreover, endogenous and exogenous variables can interact, influencing the dynamics of human populations (Lima and Berryman 2011;Lima 2014;Lima et al. 2016).
These ecological principles and models have been used to explore and explain the human dynamics of pre-and postindustrial periods ranging from small non-hierarchical societies to large empires (Royama 1992;Berryman 1999;Turchin 2001;Lima and Berryman 2011). By taking into account the effects of limiting factors and potential feedback loops (see below), here we aim to correlate recent historical events with the demography of 25 Brazilian indigenous nations and to explore the interaction between exogenous variables and population sizes to explain their observed population growth rates from the perspective of Population Ecology Theory (Royama 1992;Berryman 1999;Turchin 2001;Lima and Berryman 2011).

Methods
We analyzed the dynamics of 25 Brazilian indigenous nations (Table S1). We selected nations with at least five population size records within the period 1930 to 2010/12. We focused on the relationship among the population per capita rate of change (R), population size, and external variables. R can be estimated as the natural logarithmic ratio between consecutive population sizes or as the natural logarithmic difference between natality and mortality rates in a closed population (Royama 1992;Berryman 1999).
Using principles from System Dynamics Theory and biological interactions, population ecologists have developed a framework to explore the contributions of distinct drivers to population dynamics. R measures the average (i.e., per capita) fitness of the population, i.e., how the average individual performs in the population. Hence, R can be used as a proxy of per capita life-standard. Under a favorable environment and with sufficient resources individuals have more children and survive for longer periods (R is positive) than under unfavorable environmental conditions (R is negative) (Royama 1992;Berryman 1999).
The endogenous view holds that as the population grows individuals may become more efficient in food acquisition, finding sexual partners (i.e., changes in sex ratio), as well as in protecting themselves against predators (Allee, 1932;Berryman 1999;Luděk et al. 2006;Gregory et al. 2010). Population growth may allow for better labor organization (Boserup 1965(Boserup , 1981, reducing the costs and increasing the benefits of economic and social activities that are important to the social and cultural persistence of any indigenous population. Population growth improves individuals' life-standards, and thus accelerates its own rate of increase (R). It marks a positive feedback process biologically known as intra-specific cooperation (Berryman 1999) (S1 Fig. and S7 Text).
However, demographic cooperation is unstable; it cannot persist indefinitely. The positive demographic effects of accumulated cooperative activities eventually reach a maximum for a given environment. Life-quality decelerates and enters a new phase where population growth is independent of population size. These populations exhibit an exponential growth pattern, which is also unstable since a population cannot grow indefinitely (Royama 1992;Berryman and Turchin 2001;Ginzburg 1986;Turchin 2001; S1 Fig.).
Nevertheless, as the population increases, individuals may face per capita reductions in resource supply (e.g., space or food), initiating a decline in the population per capita rate of change (R). This process is known as intra-specific competition, and these populations display a logistic growth pattern (Ginzburg 1986;Berryman 1999;Berryman and Turchin 2001;Lima and Berryman 2011;S1 Fig. and S7 Text).
An exogenous process occurs when a variable affects another without, however, being affected by the changes it causes (Royama 1992;Berryman and Turchin 2001), e.g., climate variables, political changes, habitat destruction, the introduction of diseases and vaccines (Lima 2009;Lima and Berryman 2011).
We can summarize a general R-function for any population (Royama 1992;Berryman 1999): R is an undefined function of endogenous (N t−d where t-d is the time lag) and exogenous processes.

Endogenous Processes (N t−d )
We aimed to identify nations with similar R behavior with population size (intra-specific cooperation, competition, and exponential growth). Due to the low number of observations per nation, we first visually explored the relationship between R and population size for the 25 nations. Depending on whether R declined, increased, or stayed relatively constant with population size, we classified the nations as intra-specific competition, cooperation, and exponential growth patterns, respectively. To support our classification, we first analyzed Mucaji Yanomami time series, which covers 60 years of annual data. Then, we grouped nations that exhibited the same behavior and fitted the following discrete models (Verhulst 1838;Gregory et al. 2010) that capture the intra-population principles: Eq. 2, Eq. 3 and Eq. 4 capture, respectively, the intraspecific competition, cooperation, and both behaviors. R is the population per capita rate of change, R max represents the maximum per capita rate of change, N t − 1 is the population size at a previous observation, K is the carrying capacity, i.e., the stable population size determined by the supply of resources in the environment and the caloric demands of the individuals to survive and reproduce. Above K, individuals suffer from food scarcity and mortality is higher than natality, resulting in population decline. Below K, there is sufficient resource surplus and population grows. Q is a non-linear term that reflects the intensity of intra-specific competition around K. E denotes an unstable equilibrium population size, where, if higher, population always grows and, if lower, it tends towards extinction; T represents how fast a population deviates from the unstable equilibrium size (Berryman 1999).
We also explored the patterns of natality, mortality, and migrations rates in order to identify which contributed to the observed R behavior for the Kawaiweté, Ikpeng, Kamaiurá, and Kisêdjê, the nations for which there is available data. We used the Partial Least Squared Regression (PLSR), with R levels as the response variable and the demographic components and population size as exploratory variables. PLSR is an extension of multiple regression analyses and is useful when the number of exploratory variables is higher than or similar to the number of observations. PLSR groups non-independent exploratory variables on orthogonal linear gradients maximizing the covariance between scores and loadings of exploratory and response variables (Naes and Martens 1985;Bates and Watts, 1988;Tobias 2003;Carrascal et al. 2009).
The behavior of R in function of population size will suggest the actions of the cooperation or competition (or the lack of it), and the PLSR analysis (Kawaiweté, Ikpeng, Kamaiurá, and Kisêdjê) will suggest the demographic paths of the effects of endogenous and exogenous processes.

Exogenous Variables
Endogenous processes may interact with external variables to affect the dynamics of a population (Berryman 1999). Theory posits that external variables are determinant if they also affect the demographic components (natality, mortality, and migration) of the population per capita rate of change (Berryman 1999). Thus we reviewed the recent history of the 25 nations from the literature and websites (<http://pib.Socioambiental .org/pt>ref) looking for variables and processes that could affect R's demographic components. We performed a multivariate crosssectional analysis to explore the effects of potential exogenous variables on the dynamics. Based on our review and available data, we can rewrite the R-function as follows: Here, N t−d refers to initial population size during recovery period, which was used as a proxy of the contribution of population size to the dynamics. The threats can be divided into extractive (e.g., illegal logging, mining) and land pressures (e.g., illegal invasions and settlements) and may be acting as negative exogenous variables (S2, S3 and S4 Text and Fig. 1). The creation of new villages may represent a cost to individuals (e.g., land clearing, construction of houses, migration) but may also ease potential problems of population growth (e.g., local resource depletion) by distributing individuals throughout the territory (S1 Table). Diet items are the variety of hunted, gathered, and cultivated species used for subsistence. The number of biomes found inside the ITs (e.g., semi-arid, dense-ombrophilous, semi-decidualombrophilous, among others) (Instituto Socioambiental 2011) was used as a proxy of ecosystem diversity invoking the diversity-equilibrium theory (Odum 1953;MacArthur 1955;Elton, 1958;Chu and Tai 2001). Briefly, the theory postulates that food-webs with more species are more resistant and resilient than less diverse systems when faced by external perturbations. We could posit that indigenous populations living in contact with more biomes would be more resistant and resilient when faced by habitat loss, for example through deforestation. We used the number of governmental and nongovernmental (NGO) projects as a proxy of political awareness about the indigenous rights.
Combining the data set (Eq.5) generated in our review and the data available in the ISA website, we built a Pearson's correlation matrix to investigate the relationship among exogenous variables. We also used the PLSR with a cross-sectional approach to reduce data complexity and ascertain which variables might be more important to explain the R trend. We used the average rate of change of R (dR), which captures the general R trend as the response variable and the data cited above as the exploratory variables (S1 Fig.). We started with the entire set of variables, and step-by-step removed exploratory variables with no or small weights on the PLSR components. We present only the final PLSR results. As noted, multivariate analyses like PLSR are usually used to reduce the complex universe of measured variables. PLSR is an exploratory analysis, not a process of hypothesis testing. Hence the results from the multivariate cross-sectional analysis should be considered as hypotheses. All statistical analyses were performed in the R environment (R Core Team 2016). We fitted Eqs. 3-4 using the nls library in R through nonlinear regression analysis (Mevik and Wehrens 2007). We interpret the statistical results under the scope of Population Ecology Theory and the recent history of indigenous nations (S1 Table and S3 Text).

Endogenous Processes (N t-d )
We observed a variety of recovery patterns. Some nations showed a general competition or cooperation recovery, while others presented an exponential growth. We found nations with evidence of an initial cooperation growth and, at higher population size, competition growth. Three nations presented a distinct sequence of behaviors. Gavião Pkopje, Karitiana, and Krikati showed an initial competition phase and a later cooperation growth (Figs. 2,3,4,5,6,7,and 8).
The dynamics of the Mucaji Yanomami is the best time series that we have and reflect most of the behaviors observed in the others nations, describing an initial cooperation phase, followed by periods of exponential and competition growth (Fig. 2). The cooperation-competition model (Eq. 3) explained nearly 50% of the observed data from 1930 to 1989. When we increased the time intervals between observations, reducing the noise, the R 2 (i.e., the coefficient of determination) of the model increased (Table 1).
We classified nations based on R behavior in function of population size. We then pooled all demographic data per group. Nations grouped as cooperation, competition, and exponential growth respectively showed a general positive, negative, and independent R trend with population size (Figs. 2, 3, 4, 5, 6, 7, 8 and Table 1). Regardless of the behavior displayed, we found that R was related to more changes in natality than mortality levels (Table 2 and Fig. 7). In all cases, mortality was at lower levels and generally declined gradually (Fig. 7).

Exogenous Variables
The correlation matrix revealed some potential collinearities. Most of the variables were positively correlated to the numbers of indigenous territories and vegetation types (Table 3). The  (Table 4).   (Table 4).

Discussion
Based on our data from the recent history of 25 Brazilian indigenous nations we propose a hypothesis concerning the interaction between the effects of population size and exogenous variables from the perspective of Population Ecology Theory. The population recovery patterns of the Brazilian indigenous nations suggest the existence of intra-specific cooperation and competition processes (or exponential growth) and exogenous variables. The designation of ITs and improved access to health care have reduced mortality levels, leading to natality-driven population per capita rate of change. At the same time, the expansion of agricultural frontiers has been degrading environmental suitability for population growth (Ribeiro 1970; Villas Bôas and Villas Bôas 1994; Instituto Socioambiental 2011; S3 Text).
Together with health care, officially designated ITs are central to the social and cultural persistence of indigenous populations. The correlation matrix showed that the number of ITs correlates with most other variables. Increasing the number of ITs per nation raises the possibility of confronting new threats and new biomes (positive correlations), and allows higher initial population sizes and the potential for the creation of more villages.
The PLSR results point to the importance of deforestation, population size, and diet to nations' demography. We focus on the results from the first component of the analysis and suggest that the number of biomes may be a proxy for deforestation. Most nations analyzed in this study are from the Central-Western region, which is a transitional area between the dense-ombrophilous forest (Amazon) and other biomes (Brando et al. 2013), presents the highest diversity of the biomes. While this high richness of biomes could invoke the diversity-equilibrium theory (Odum 1953;MacArthur 1955;Elton, 1958;Chu and Tai 2001), in our case of Brazilian indigenous nations the predictions of the diversity-equilibrium theory did not Fig. 6 The R functions of the nations with both accelerating and declining R behaviors. The initial decline trend of the Krikati, Karitiana and Gavião Pykopjê can be related to the costs prior the land homologation apply (see Methods). First, the dense-ombrophilous forest (the Amazon in the Northern region) present more species than all the other biomes together (Brando et al. 2013). Secondly, the Central-Western region has been suffering from the advance of deforestation for a longer period than the Northern region (Fearnside 2005;Paula 2005;Salomão et al. 2007; Barbosa et al. 2008;Brando et al. 2013). The region with more vegetation types is also where agricultural frontiers have been expanding during the last decades. Therefore, we suggest that vegetation types and habitat degradation correlate geographically and may be a proxy of habitat impoverishment. Hence, the PLSR first component can be interpreted as a combination of deforestation and population size effects on dR.
The Mucaji Yanomami dynamics provide the best evidence of population size effects from our sample of nations. They live in the Northern region (on the frontier with Venezuela) and are therefore more distant from the negative exogenous variables than the other nations we analyzed, located in the Central-Western region. Hence, we can interpret the Mucaji time series as being driven mainly by the effects of population size on the population per capita rate of change (R). Their population started to recover prior to the introduction of health care however, and the initial cooperation recovery phase can be related to sex ratio stabilization and changes in their reproductive regime as they moved from polyandry to monogamy, which is a demographic Allee process with cultural changes (Early and Peters, 1990;Albert 1999).
Habitat degradation could interact with population to explain the dynamics of nations from Central-Western and Northeastern regions. Deforestation reduces availability of big game, fish stocks, and both timber and a range of non-timber forest products (NTFPs) (Ribeiro 1970; Vi l l a s B ô a s a n d Vi l l a s B ô a s 1 9 9 4 ; I n s t i t u t o Socioambiental 2011). We propose that as the advance of agricultural frontiers degrades environmental sustainability, it may increase the importance of cooperative activities, which in turn impact demography (Lima and Berryman 2011;S5 Text). For example, in the PIX, almost all nations with more villages near the agricultural frontiers showed a general R positive trend (Fig. 9). Fig. 7 Indigenous nations' available natality and mortality rates. Shows that independent of the behavior, R is closer related to natality than mortality rates. All nations are from PIX During the R acceleration phase, population growth together with cooperative activities may improve environments previously degraded by deforestation (S8 and S9 Text). Population growth may diminish individual costs of foraging (hunting and gathering), crop cultivation, house construction, and caring for children (Ribeiro 1970;Berryman 1999;Smith and Winterhalder, 2003;Hamilton et al. 2009).
However, this allows for the potential existence of a minimum population size (E) below which the population becomes susceptible to extinction (Berryman 1999). Indigenous nations have two E, one biological and other cultural (Pozzobon 1994). Any human population needs a minimal population size to function culturally as a group. Some activities both social and related to subsistence (e.g., group hunting or fishing) require a minimum number of individuals involved (Ribeiro 1970;Pozzobon 1994),. Depending on the degree of resource stock reductions, habitat destruction may also remove activities related to subsistence even when there are enough individuals to perform them (Flowers 1994), pushing the cultural and numerical E to higher population sizes, increasing the odds of extinction during negative exogenous events (e.g., epidemics, drought) within populations with demographic evidence of the intra-cooperation process. These nations should invest in population size to counteract approaching the numerical and cultural E. Investing in population size to counter the negative effects of deforestation may also allow better territory defense and attract more media attention, which may translate into increased and improved governmental assistance and protection of nations´human rights, as occurred for the Karitiana, Krikatí, and Gavião Pykopjê (Ladeira, 1991;Ladeira and Azanha 1996, 2005a, 2005b. Our results also captured the potential negative impact of deforestation on diet richness leading to compromised levels of nutrition, mentioned in nearly all the nationsć hronicles (Junior 1985(Junior , 1989Junior and Santos 1991;Fagundes et al. 2004;Orellana et al. 2006; S1 Table; S2, S4, S8 Text). This suggests that habitat destruction could exacerbate the decline of the per capita resource availability as population grows, as suggested by recent Xavante history and dynamics (S9 Text). Habitat destruction could result in a lower carrying capacity (K), also increasing the odds of extinction during negative events. Deforestation increases dependence on family farm  production (Villas Bôas and Villas Bôas 1994;Flowers 1994;Junior et al. 2002;Instituto Socioambiental 2011;Nilsson and Fearnside, 2011;S9 Text). Within Brazilian indigenous communities, agriculture is mainly women's responsibility (Ribeiro 1970;da Silva, 1999), so that any increase in agricultural production could potentially produce a negative trade-off between family size and overall demographic fitness (Kaplan 1996;Hill and Kaplan 1999). The intra-specific competition phenomenon implies that eventually all nations will suffer from the negative effects of population growth (Berryman 1999). The ultimate goal of indigenous nations and Brazilian governments should be to create conditions that enable nations to attain high and sustainable equilibrium population size (K). Indigenous nations with demographic evidence of intra-specific competition should invest in the creation of more villages to reduce the negative effects of local resource depletion. Diluting the population could allow higher and sustainable K while the Brazilian government continues to invest in the expansion of agricultural frontiers. The Brazilian government should invest in health care improvements, the creation of more ITs and villages to ease the negative effects of population growth.
Political and economic interests are the main impediments to the creation of new ITs and the extension of existing ones, as evidenced by the lack of impact of governmental and NGO projects on the demography of indigenous communities. International economic growth (e.g., China) increases demand for Brazilian natural commodities, pushing mono-agricultural frontiers closer to indigenous territories (IBGE 2014). As a significant proportion of Brazilian GDP, agribusiness has a strong political representation in the Brazilian congress, and its supporters are generally opposed to indigenous rights. Variables are: number of territories (N ITs), total territory size (Territory, Km 2 ), initial population size (Population), number of villages (Villages), extractive pressure (Extractive), land pressure (Land), gathering diet items (Gather), agricultural items (Agriculture), hunting items (Hunt), government and non-government projects (Gov-NGO). Significant results in bold (P < .05) Consequently, the number of officially designated ITs has been declining drastically during the recent governments, while deforestation in favor of agriculture has been advancing (Instituo Socioambiental 2011). The platform of the recently elected president Jair Bolsonaro seems to represent a continuation of this trend, so it seems likely that indigenous communities will continue to suffer as a consequence of the belief of non-indigenous society that they represent a hindrance to economic growth (S10 Text). Brazilian GDP estimates do not directly include the benefits of ecosystem services from natural reserves (IBGE 2014). Yet habitat degradation and the loss of natural services do impact actual GDP. For example, deforestation and urbanization in the Amazon region increased the frequency of malaria outbreaks and consequently governmental costs on health care (Amaral et al. 2001;Loioloa et al. 2002). Also, the Amazon forest influences rainfall regimes through humidity to the southeastern region, which has the most populous and important industrial centers in Brazil. Therefore, rain deficits, a product of deforestation, may impact energy costs and industrial productivity. In 2017, the government sanctioned a new law (No. 13.493/17) that incorporates costs of ecosystem services into GDP estimates, but how these estimates will be calculated and when the law will be implemented are still under debate. If and when it is implemented, the law may change the actual paradigm regarding hindrances to economic growth currently faced by ITs.
The United Nations Sustainable Development Goals (United Nations 2017) include the target of BPromot[ing] mechanisms for raising capacity for effective climate change-related planning and management in least developed countries and small island developing States, including focusing on women, youth and local and marginalized communities.^Brazil has been and seems to be set to continue to follow policies contrary to such efforts to reduce the impacts of climate change as it continues to fuel the expansion of agricultural frontiers. The loss of nearly 800 distinct languages during the last centuries represents not only a demographic genocide, but also the loss of the indigenous human capital in the form of local environmental knowledge (LEK) acquired over millennia, including knowledge of food and medicinal plant species as well as the natural cycles of ecosystems and the services they provide. 2 Also, LEK can be important to reforestation programs, as observed today in PIX where indigenous communities collect the seeds for a project aiming to rehabilitate the degraded headwaters of Xingu's hydrographic basin. It is also notable that smallscale agriculture as practiced in indigenous communities allows for greater diversity than commercial monocultures and is less affected by pests, reducing the demand for and use of pesticides (Instituto Socioambiental 2011).
The demographic recovery of indigenous populations in Brazil allowed better territory protection and surveillance as well as intra-and inter-nation social mobilization to pressure the government to guarantee land rights that potentially impact the demographic recovery pattern. For example, the population growth of the Karitiana, Krikatí, and Gavião Pykopjê nations allowed better internal organization and also attracted media attention that ultimately led to guaranteed land rights, changing the R decline to increase. The inter-nations' mobilizations together with the efforts of pro-indigenous NGOs and universities act to counteract the agribusiness agenda of the government (see https://pib.socioambiental.org/pt/Lista_de_ organiza%C3%A7%C3%B5es_de_apoio_aos_povos_ind% C3%ADgenas).
Education may be the best way to construct a future social awareness in Brazilian society as well as the international community of the current situation and potential contributions of indigenous nations, and most particularly to counteract the secular paradigm that indigenous nations are a threat to economic growth. The new national curricular base in Brazil (approved by the congress and the senate in 2017) clearly includes under History the valorization of afro-and indigenous cultures (Law n°10.639/200349 and Law n°11.645/200850). If successfully implemented, these new guidelines might change social and political awareness of indigenous rights established in the 1988 Constitution (S10 Text).

Conclusion
Most of the population predictions of human dynamics are deterministic or stochastic projections of current 2 For example, the Kawaiweté recognize more than 200 species of plants that can be cultivated for agricultural purposes (Instituto Socioambiental 2011). Variables with low or no weights on PLS components were removed from the analysis and are not shown. W1, W2 and W3 represent the loading weights of each variable in the first, second and third PLS regression components. R 2 % is the proportion of the variance in the response variable accounted for by each component of the PLS regression. All predictor variables were standardized trends based either on exponential or faster than exponential growth models (Galor and Weil 2000;Lutz et al. 2006). These approaches are problematic because they are not based on general principles of population dynamics. Consequently, most of the human population dynamics forecasts derive from inappropriate models because they do not take into account the effects of limiting factors and potential feedback loops (Royama 1992;Berryman 1999). By exploring the dynamics and the recent history of indigenous nations with the proper ecological principles we were able to correlate exogenous variables and processes (such as the introduction of health care and advances of agriculture frontiers) with population size. The ecological principles helped to disentangle the complex universe of variables in which indigenous populations live. The framework encourages hypothesizing and testing processes that could affect population dynamics and would be readily applicable across a range of human societies.