Alexander Koch: Earth System Impacts of the European arrival and Great Dying in the Americas after 1492

* Science Direct: Earth System Impacts of the European arrival and Great Dying in the Americas after 1492 [PDF: Copy PDF]; The Guardian: European colonization of America’s killed so many it cooled Earth’s climate [SQ Copy] » Genghis Khan Reforestation Effect: Julia Pongratz: Genghis Khan the GREEN: Invader killed so many people that carbon levels plummeted.

Earth system impacts of the European arrival and Great Dying in the Americas after 1492

Alexander Koch,  Chris Brierley, Mark M. Maslin, Simon L. Lewis, Geography, University College London, London, WC1E6BT, United Kingdom; School of Geography, University of Leeds, Leeds, LS29JT, United Kingdom

01 March 2019, Science Direct.



Human impacts prior to the Industrial Revolution are not well constrained. We investigate whether the decline in global atmospheric CO2 concentration by 7e10 ppm in the late 1500s and early 1600s which globally lowered surface air temperatures by 0.15+C, were generated by natural forcing or were a result of the large-scale depopulation of the Americas after European arrival, subsequent land use change and secondary succession. We quantitatively review the evidence for (i) the pre-Columbian population size, (ii) their per capita land use, (iii) the post-1492 population loss, (iv) the resulting carbon uptake of the abandoned anthropogenic landscapes, and then compare these to potential natural drivers of global carbon declines of 7e10 ppm. From 119 published regional population estimates we calculate a pre-1492 CE population of 60.5 million (interquartile range, IQR 44.8-78.2 million), utilizing 1.04 ha land per capita (IQR 0.98-1.11). European epidemics removed 90% (IQR 87-92%) of the indigenous population over the next century. This resulted in secondary succession of 55.8 Mha (IQR 39.0-78.4 Mha) of abandoned land, sequestering 7.4 Pg C (IQR 4.9-10.8 Pg C), equivalent to a decline in atmospheric CO2 of 3.5 ppm (IQR 2.3-5.1 ppm CO2). Accounting for carbon cycle feedbacks plus LUC outside the Americas gives a total 5 ppm CO2 additional uptake into the land surface in the 1500s compared to the 1400s, 47-67% of the atmospheric CO2 decline. Furthermore, we show that the global carbon budget of the 1500s cannot be balanced until large-scale vegetation regeneration in the Americas is included. The Great Dying of the Indigenous Peoples of the Americas resulted in a human-driven global impact on the Earth System in the two centuries prior to the Industrial Revolution. © 2019 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license


1. Introduction

The arrival of Europeans in the Americas in 1492 CE marks the onset of disease epidemics resulting in the loss of the majority of indigenous people living in the Americas over the subsequent century (Berlinguer, 1993; Cook, 1998; Crosby, 1972, 1976; Nunn and Qian, 2010). Indigenous land use was widespread before European arrival, particularly in Mexico, Central America, Bolivia and the Andes where terraced fields and irrigated agriculture was practised (e.g. Abrams and Nowacki, 2008; Chepstow-Lusty and Jonsson, 2000; Heckenberger et al., 2003; Hunter and Sluyter, 2015; Whitmore and Turner, 1992), and across Amazonia where diverse pre-Columbian land uses left its traces in the composition of contemporary Amazon forests (Clement et al., 2015; Levis et al., 2017; Maezumi et al., 2018a; Watling et al., 2018). Thus the Great Dying of the Indigenous Peoples of the Americas (Wolf, 1982) after 1492 CE likely led to a reduction in land use. Fields and fallow areas then underwent secondary succession and in many cases increased carbon stocks as they reverted towards similar prior states e with local, regional and potentially global consequences for the Earth System. The uptake of carbon on the abandoned anthropogenic lands after European contact may have been large enough to impact the atmospheric CO2 record (Dull et al., 2010; Faust et al., 2006; Lewis and Maslin, 2015; Nevle and Bird, 2008; Nevle et al., 2011; Ruddiman, 2005). Furthermore, at the same time high-resolution Antarctic ice-core records of atmospheric CO2 concentration show an anomalously large decline of ~7e10 ppm (Ahn et al., 2012; MacFarling Meure et al., 2006) beginning in the 1500s with a minimum in the early 1600s (Fig. 1). Isotope analysis shows that the anomaly was driven by an increase in the terrestrial carbon sink (Fig. 1B, Bauska et al., 2015; Francey et al., 1999; Trudinger et al., 1999). Hence, the carbon uptake that is thought to have occurred following the arrival of epidemics in the Americas may have reduced atmospheric CO2 levels and led to a decline in radiative forcing that may then have contributed to the coldest part of the Little Ice Age (Faust et al., 2006; Neukom et al., 2014).

Fig. 1. Global atmospheric CO2 time series. (A) The CO2 concentrations recorded in two Antarctic ice cores: Law Dome (grey, MacFarling Meure et al., 2006) and West Antarctic Ice Sheet (WAIS) Divide (blue, Ahn et al., 2012). (B) Carbon isotopic ratios recorded in CO2 from the WAIS Divide ice core (black, Bauska et al., 2015) showing an increased terrestrial carbon uptake over the 16th century (B). Yellow box is the span of the major indigenous depopulation event (1520e1700 CE). Loess smoothed lines for visual aid. Error bars are the 1-s standard deviation. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

However, modelling efforts attempting to reproduce the CO2 decline from land use change (LUC) have so far been inconclusive, some showing an impact (Kaplan et al., 2011), while others do not (Pongratz et al., 2011). An alternative approach, using a forcing- fingerprint analysis, rules out internal variability of the Earth System, changes in the total solar irradiance (TSI), and volcanic forcing, as drivers of the change (Schurer et al., 2013). This study finds that the decline in atmospheric CO2 contributed to the observed cooling, which is consistent with uptake following secondary succession in the Americas.

We therefore test the hypothesis that human actions impacted CO2 levels and global surface air temperatures in the 16th and 17th centuries, just before the Industrial Revolution, by examining the evidence at each stage of its primary assumptions. First, the existence of a sufficiently large indigenous population in the Americas before 1492 CE. Second, the impact of anthropogenic land use on the environment. Third, the population decline estimates following the arrival and spread of European diseases. Fourth, the magnitude, extent and timing of the subsequent successional regrowth of vegetation and resulting carbon sequestration. We therefore compile all existing estimates from a comprehensive literature survey. These records are then cross-combined and sampled using two different statistical methods (Fast Fourier and Monte Carlo approaches) to obtain revised estimates of the 1492 CE population, per capita land use, mortality from European-contact diseases, and the subsequent uptake of carbon on the abandoned lands. For each term we also compute uncertainty estimates. We then combine these estimates to show the most likely impact of the Great Dying of the Indigenous Peoples of the Americas on atmospheric CO2 levels. Finally, we place these new estimates of carbon uptake in the context of other changes to the 16th century carbon cycle.

2. The population of the Americas in 1492

The first population groups to arrive in North America between 20,000 and 15,000 years ago were of east Asian- and north Eurasian ancestry (e.g. Amick, 2017; Moreno-Mayar et al., 2018; Pedersen et al., 2016). The adoption of a sedentary, agricultural way of life in the Americas began 10,000e8000 BP (Piperno and Dillehay, 2008; Piperno et al., 2009). Large, complex civilizations emerged in North, Central and South America, further increasing population density, with abundant evidence for a large population living in the Americas prior to European arrival (e.g. Bolt and van Zanden, 2014; Canuto et al., 2018; Clement et al., 2015; Cook, 1998; Denevan, 1992a). However, as the epidemics spread, often ahead of the European explorers, pre-European population estimates were never formally documented in colonial censuses (Diamond, 1997; Dobyns, 1993). Hence, Native American populations were only documented in the decades after European contact, necessitating the use of indirect methods to estimate the pre-1492 population of the Americas. Here, we summarize all existing methodologies applied to reconstruct pre-contract indigenous population, and extract published regional estimates to calculate a revised hemispheric population estimate with uncertainty ranges.

2.1. Approaches to estimating the indigenous population at European arrival

Historical anthropologists and archaeologists estimate indigenous population numbers at contact by utilizing a range of documentary evidence, including sizes of armies, number of adult males, census data, tribute records, numbers of buildings, depopulation ratios, and historic clerical chronicles such as the number of baptisms and number of deaths in a community (Borah and Cook, 1960, 1969; Dobyns, 1966, 1983). The conversion from such evidence into population estimates has several limitations. Records of the sizes of armies and the sizes of settlements from just after Spanish arrival may be prone to miscounting and exaggeration (Denevan, 1992a) while the conversion from tributes to population numbers requires contentious assumptions to be made (Henige, 1998; Zambardino, 1980). These include the proportion of a population excluded from tributes, average tribute paid per house, number of people per house, conversion from goods into monetary value, spatial homogeneity of population structure and analogies to present day population structures (Borah and Cook, 1960; Sanders, 1976; Sanders et al., 1979).

Colonial census estimates are generally considered more reliable (Denevan, 1992a), but do not capture population levels before European contact (Borah and Cook, 1960). Accepting such censuses as representative of pre-contact populations led to studies from the early and mid-20th century giving very low 1492 CE population estimates (e.g. Kroeber, 1939; Rosenblat, 1954). Some estimate of the population loss on contact is necessary to provide more robust population numbers.

Archaeological data can also be used to estimate pre-contact population size, but again involves several assumptions, such as the number of houses in a settlement that are occupied at one time (Schacht, 1984). But in contrast to historical documents, new archaeological sites are regularly discovered and new and innovative multidisciplinary approaches are being applied to the data, adding new constraints on regional pre-contact population sizes (e.g. Canuto et al., 2018; Goldberg et al., 2016; Liebmann et al., 2016; Storey, 2012). Population densities can be inferred from the number of archaeological sites (Heckenberger et al., 2003), the number of archaeological features such as pottery, charcoal and fertile anthropogenic soils known as Amazonian Dark Earth (ADEs Heckenberger et al., 1999; Meggers, 2001). Expanding on this, the “habitat density method” projects well-studied local population densities to other regions with similar environmental characteristics to calculate the total population under the assumption that population densities in different locations with the same habitat are similar (e.g. Denevan, 1970, 2003; Newson, 1996). Thus, all 1492 CE population estimates require data to be combined with assumptions to arrive at estimates.

2.2. New pre-Columbian population estimates

We summarize the existing evidence for pre-Columbian populations using seven geographical regions: the Caribbean, Mexico, Central America, the Inca Territory, Amazonia and contiguous forested area, North America and the Rest of the Americas. The Inca Territory at contact encompassed Peru, Bolivia, Ecuador, Colombia, Chile and parts of north-western Argentina; the Rest of Americas is composed of Venezuela, Uruguay, Paraguay and the rest of Argentina. We complied 129 estimates from 82 studies (Caribbean, n ¼ 18; Mexico, n ¼ 17; Central America, n ¼ 23; Inca Territory, n ¼ 26; Amazonia, n ¼ 20; North America, n ¼ 9; Rest of Americas, n ¼ 6; total hemisphere, n ¼ 10; Table 1). Below we discuss the range of estimates for each region and use them to calculate a total population estimate of the Americas in 1492 CE based on all existing data available at the time of this review.

Caribbean. The first contact between Europeans and Amerindians occurred when Columbus arrived in 1492 CE on the Caribbean island of Hispaniola (what is now Haiti and the Dominican Republic). The primary source for most modern estimates of the contact population is a contemporary eyewitness report by the Dominican friar Bartolome de las Casas who reported a Hispaniola  population of 4 million at the time of contact (Denevan, 1992a). However, studies based on this same report range from 60,000 that neglects contemporaneous reports of post-contact deaths (Verlinden, 1973, in Henige, 1978), to acknowledging some degree of depopulation with initial populations estimated at between 100,000 (Amiama, 1959; Rosenblat, 1976) and 8 million (Cook and Borah, 1971, Table S1). The 8 million estimate has drawn criticism for the worker-to-population ratio used, the assumption of a logarithmic population decline during 1492e1496 CE, and the inclusion of a controversial tribute count from 1496 CE (Henige, 1978). Most estimates are between 300,000 and 500,000 people in the Caribbean before European contact (Cordova, 1968  ; Dobyns, 1966; Morison, 1948; Moya Pons, 1979; Williams, 1970).

Mexico. Europeans arrived in what is now Mexico in 1519 CE. Yet, the first comprehensive census only took place in 1568 CE, and tribute records are only available as documentary evidence from the late 1540s (Cook and Borah, 1960; Whitmore, 1991). For central Mexico (the most populous region in Mexico) a population of 25.3 million in 1519 CE has been estimated from tribute records (Borah and Cook, 1969), which would make it one of the most populous regions in the world at the time. Zambardino (1980) highlights the uncertainties in Borah and Cook’s approach and arrived at a population between 2.2 and 28 million. Dobyns (1966, 1983) suggests the highest numbers, up to 52 million for central Mexico, based on an extrapolation of depopulation ratios from the Valley of Mexico. These values are higher than some estimates for all of Mexico. For the other region with a substantial pre-Columbian population, the Yucatan peninsula, estimates range between 2.3 million (historic documents, Lange, 1971) and 13 million (extrapolation from number of houses, Morley, 1968). Considering all of Mexico, the lowest estimates, which either outright reject (Rosenblat, 1954) or question the accuracy (Sanders, 1972) of the early tribute records, are between 4.5 and 6 million. Intermediate estimates are ~16 million, based on a linear extrapolation from results of a modelling study incorporating a combination of population, agricultural and epidemiological models for the Valley of Mexico (Fig. 2, near Valle de Mezquital, Whitmore, 1991) and 17.2 million as a best-estimate based on a synthesis of previous estimates by (Denevan, 1992a). Pre-contact population estimates for central Mexico and Yucatan combined, which is considered representative for all of Mexico, range from less than 3 million to over 52 million with many at around 20 million (Table S1).

Central America. Population numbers for Central America (Belize, Costa Rica, El Salvador, Guatemala, Honduras, Nicaragua, Panama) are based on sums of regional population estimates (Lovell and Swezey, 1982), reports on indigenous army sizes, tributes and carrying capacities (Fowler, 1988) as well as depopulation ratios (Dobyns, 1966). Estimates range from 0.8 million, based on the extrapolation of census data, omitting any significant population changes prior to the census date (Rosenblat, 1954) to 10.8e13.5 million based on depopulation ratios from better-documented Mexico (Dobyns, 1966), similar to a recent 7e11 million estimate for the Late Classic Maya period (650e800 CE, Canuto et al., 2018). Most estimates range between 4.75 million and 6 million (Sapper, 1924; Denevan, 1992a; Lovell and Lutz, 1995, Table S1).

Inca Territory. At its height, just before conquest in 1533 CE, the Inca Empire encompassed Peru, Bolivia, Ecuador, southern Colombia, Chile and parts of north-western Argentina (Fig. 2). The lowest estimate for the pre-Columbian Inca empire population is two million, based on a reconstructing depopulation rates from tribute records (Shea, 1976). The most comprehensive study on population dynamics of the Inca heartlands of Peru and Bolivia (Cook, 1981), is 9 million (3.3e14.2 million range) at the time of contact, derived from carrying capacity calculations, mortality figures and rates of change between colonial censuses throughout time. Subsequent studies tend to favour the higher end of these estimates (9e11.7 million, Montenegro et al., 1990; Denevan, 1992a). For the other regions of the Inca realm populations are suggested to be relatively small, between 300,000 and 1.5 million in each of Bolivia, Ecuador, Chile and Argentina, plus ~3 million in Colombia (Denevan, 1992a; Etter, 2000; Ocampo, 1997; Villamarín and Villamarín, 2000), with a range of 850,000 (Rosenblat, 1954) to 5 million (Colmenares, 2007). Overall Inca Territory estimates range from 4.1 million to 43.8 million with a likely population of around 20 million, based on the sum of the most widely accepted figures for each of the regions (Table S1).

Large: Table 1. Population and land use at 1500 CE and 1600 CE in P08 (Pongratz et al., 2008a), HYDE 3.1 (Klein Goldewijk et al., 2010), KK10 (Kaplan et al., 2011), and this study. Also shown are their implied net global carbon uptakes as published in Pongratz et al. (2011) (for P08) and Kaplan et al. (2011) (for KK10 and HYDE 3.1) and the calculated net global carbon uptake from this study.
Fig. 2. Regions known to have been affected by disease outbreaks by 1600 CE and pre-Columbian land use. Disease outbreak data derives from Dobyns (1993) and Cook (1998). Land use extent is based on high densities of archaeological sites, earthworks and Amazonian Dark Earths occurrence (Clement et al., 2015; De Souza et al., 2018; McMichael et al., 2014; Whitmore and Turner, 1992). The locations of sites and records referred to in the text and later figures are also shown. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

Amazonia. The vast drainage basin of the Amazon and contiguous forested areas are assumed to be regions of relatively low population density (Denevan, 1992a). Based on archaeological excavations, as historical documents are sparse, larger populations are assumed to have primarily settled in the lowland and lower upland drainage area of the Amazon and its tributaries (Clement et al., 2015; Denevan, 1992a; Roosevelt et al., 1991) as well as at the edges of savannas in Bolivia and Brazil (Erickson, 2000; Heckenberger et al., 2003) and along the northeastern coastline of South America (Rostain, 2008; Schaan, 2008, see Fig. 2). Population studies give a range between 1 million (Rosenblat, 1954) and 20 million (Oyuela-Caycedo in Denevan, 2012) living in Amazonia. Estimates include 1.5e2 million based on an average of present-day densities (Meggers, 1992), 3.2 million based on tribe-by-tribe counts (Hemming, 1995), 5.5 million extrapolated from eastern Ecuador (Newson, 1996), and from 5.1 to 20 million based on the extrapolation of studied population densities along floodplains and known regions of ADEs scaled to an estimated total extent of ADEs (Clement et al., 2015; Denevan, 1970, 1992b, 2012, 2014; Smith, 1999; Woods et al., 2013, see Fig. 2). The large range of the last approach is due to (i) uncertainty in the extent of ADEs; (ii) the proportion occupied at contact, and (iii) settlement size uncertainty (Woods et al., 2013). Recent findings, such as the early domestication of Amazon rice ~4000 BP (Hilbert et al., 2017; Watling et al., 2018), indicate larger populations, with most recent estimates ranging between 8 and 20 million people before Europeans arrived (Table S1).

North America. Europeans began to colonize North America (defined here as the United States of America and Canada) after Central and South America, thus regional and continent wide estimates are largely based on archaeological evidence, tribe-by-tribe counts and environmental carrying capacities (Milner and Chaplin, 2010, see Table S1). The lower range of population estimates for North America lies between 900,000 and 2.4 million, based on tribe-by-tribe counts for the period 1600 CE to mid-1800 CE (Kroeber, 1939; Ubelaker, 1976, 2006). The highest estimate of 18 million, established from analysing environmental carrying capacities (Dobyns, 1983), has been criticized for its assumptions on food acquisition strategies (Milner and Chaplin, 2010). More recent estimates derived from geospatial interpolation of archaeological sites range between 2.8 million and 5.7 million (Milner and Chaplin, 2010). These intermediate figures are supported by a recent comprehensive regional-scale archaeological study (Liebmann et al., 2016).

Rest of the Americas. Not included in the other regions are estimates for Venezuela with 600,000e1.5 million, Uruguay and Paraguay, estimated together as 285,000e1.1 million, and Argentina with 300,000e500,000 people (Denevan, 1992a; Pyle, 1976; Rosenblat, 1954). All low estimates are based on census data (Rosenblat, 1954) while the higher bounds are based on precensus eye-witness accounts for Argentina (Pyle, 1976) and a synthesis of previous estimates for Venezuela, Uruguay and Paraguay (Denevan, 1992a). The total estimate for the remainder of the Americas is between 1.2 and 3.1 million.

We then synthesize the 119 published population estimates, including the minimum and maximum estimate when given (Caribbean, n ¼ 18; Mexico, n ¼ 17; Central America, n ¼ 23; Inca Territory, n ¼ 26; Amazonia, n ¼ 20; North America, n ¼ 9; Rest of Americas, n ¼ 6; Table 1) into a revised estimate. We included all the prior studies and did not make any judgement on their relative quality. To obtain regional population estimates, we first calculated estimates by country. For sub-national population data, we crosscombine each sub-national estimate from one part of a country with each sub-national estimate from another part of the country and calculate their sums to obtain national estimates (e.g. each estimate for Yucatan is paired with each estimate for Central Mexico, the sum of those forms estimates for Mexico). These estimates are then added to the existing estimates for their country (e.g. Mexico). National estimates within a region are crosscombined and their sums form a regional estimate for each of the seven regions: Caribbean, Mexico, Central America, Inca Territory, Amazonia and North America, and the Rest of the Americas (Table S1). Next, cross-combining and taking the sums of these regional estimates (2 1028 combinations) gives a hemispherewide population frequency distribution, with the higher occurrence rate of similar results reflecting higher frequencies in the distribution (Fig. 3A). Our estimate of the number of people living in the Americas in 1492 CE is 60.5 million, with an interquartile range (IQR) of 44.8-78.2 million, using the FFT approach (Fig. 4A). An alternate Monte-Carlo approach to combining the data results in a similar best estimate of 64.1 million (IQR 48.4-82.3 million; Table S2). When sufficient data is available we adopt the crosscombination approach for the rest of the variables we estimate later in the paper. The programming code for both approaches in available in the supplement.

These results are significantly higher than the often-quoted hemisphere 1492 CE population of around 40-41 million, based on Clark’s 1967 literature review (Biraben, 1979; Durand, 1977; Henige, 1998; Klein Goldewijk et al., 2010, 2011; Pongratz et al., 2008b), but lower than the Dobyns (1966, 1983) estimate of 112 million. Our results are more similar to Denevan (1992a) 54 million, incorporating another 25 years of research findings after Clark’s 1967 review, and the 60 million estimate used in the LUC dataset KK10 (Kaplan et al., 2011; Krumhardt, 2010, Table 1). For comparison, 60.5 million people living in the Americas in 1492 is lower than the population of Europe at the beginning of the 16th century, at 70-88 million (Durand, 1977) and that of China (including Mongolia) at about 100 million (Krumhardt, 2010), both larger populations in a smaller area (Europe 10 million km2 , China and Mongolia 11 million km2 , Central and South America, 18 million km2 ). Thus, population density was lower in the Americas than in Asia or Europe at the time of European contact.

Fig. 3. Reconstructed spatial extent of anthropogenic land use at 1500 CE and 1600 CE in P08 (Pongratz et al., 2008a), HYDE 3.1 (Klein Goldewijk et al., 2011) and KK10 (Kaplan et al., 2011). The datasets vary in their underlying assumptions about land use (see text for details), leading to visibly different extents (Table 1).

3. Extent of anthropogenic land use in 1492

It is widely accepted that the indigenous population of the Americas in 1492 CE did not live in a pristine landscape (e.g. Barlow et al., 2012; Denevan, 2001; Heckenberger et al., 2003; Levis et al., 2017; Whitmore and Turner, 1992). The extent and complexity of pre-Columbian agricultural systems is well-constrained on a local scale for many areas based on a variety of evidence including historical documents, aerial photography and newer remote sensing techniques as well as field based archaeological studies (Canuto et al., 2018; De Souza et al., 2018; Denevan, 2001; Hunter and Sluyter, 2015; Liebmann et al., 2016; Loughlin et al., 2018; Maezumi et al., 2018a; Palace et al., 2017; Whitmore and Turner, 1992). Given the diversity of pre-Columbian land-use practices in the Americas, we first discuss the common land-use practices within our seven regions, then synthesize published land use per capita estimates for the regions, and the extent of pre-Columbian land use in each region to calculate a best estimate of anthropogenic land use in 1492 CE in the Americas.

Fig. 4. Probability Density Functions showing frequencies of the indigenous population, in millions, of the Americas in 1500 CE (A); anthropogenic land use area in ha for 1500 CE (B); indigenous population after the Great Dying, in 1600–1650 CE (C); area of anthropogenic land abandoned due to the Great Dying, in Mha (D); carbon uptake in vegetation and soils over 100 years following anthropogenic land abandonment after the Great Dying, in Pg C (E); reduction in atmospheric carbon dioxide over 100 years following anthropogenic land abandonment after the Great Dying, in ppm (F); each derived from combining published regional estimates (see text for method). Bold lines indicate 25–75% interquartile range (IQR), narrow, dashed perpendicular line indicates the median, i.e. where cumulative frequencies pass 50%.

3.1. Land use practices in the pre-Columbian Americas

Caribbean. On most islands, semi-permanent or permanent anthropogenic land use was widespread by the time Europeans arrived (Ortiz Aguilu et al., 1991). Farming consisted mainly of mixed sweet potato (Ipomoea batatas), cassava (Manihot esculenta) and maize (Zea mays) as staple crops eaten together with protein sourced from fishing (Mickleburgh and Pag an-Jimenez, 2012  ). Stone-built terraces with drainage systems on the hillsides of Puerto Rico (Ortiz Aguilu et al., 1991), raised fields (Castilla-Beltran et al., 2018) and semi-permanent agriculture on Hispaniola (Richards, 2006) suggest that Caribbean land use was comparable to that in Central America and populous regions of Amazonia. Fire was periodically used for slash and burn systems and small clearings that increased after European settlement (Castilla-Beltran et al., 2018; Hooghiemstra et al., 2018).

Mexico. Large areas of land were under cultivation of maize fields, cacao and fruit orchards and house gardens (Whitmore and Turner, 1992). The remains of elaborate canal systems have been found in the wetlands of the Gulf Coast and the Basin of Mexico (Whitmore and Turner, 1992). The 700 km2 Valle del Mezquital in central Mexico (Fig. 2) was extensively managed with fields, some with trees as agroforestry-type systems, irrigated through complex systems of dams and weirs connecting waterbodies in the basin and terraces along the hillsides (Hunter and Sluyter, 2015; Whitmore and Turner, 1992). A distinctive form of Mayan slash and burn agriculture, milpa, is still practiced in Yucatan. It is a combination of clearing, burning, planting, cropping and short fallow cycles (~4 years) with small plots of about 4 ha per family (Turner and Brush, 1987).

Central America. Along with Mexico, Central America was one of the most intensely used landscapes in the western Hemisphere with agricultural systems similar to central Mexico (Whitmore and Turner, 1992). In Costa Rica, pollen records from Lake Zoncho (Fig. 2) indicate that the local population moved from slash and burn agriculture towards permanent agriculture in the later stages of the Mayan Postclassic period (1000e1539 CE) most likely to avoid laborious repetitive clearing of land, due to the epidemicsdriven decline in workforce (Clement and Horn, 2001). The abundant occurrence of charcoal in records suggests that fire was widely used as agricultural tool in Central America (Dull et al., 2010).

Inca Territory. Andean landscapes have been modified for over 8000 years (Dillehay et al., 2007), with evidence for the terracing of steep slopes to grow maize, quinoa (Chenopodium quinoa) and other staples appearing between 550 CE and 1000 CE (Kendall, 2013). While the actual extent of terracing fields in the Andean highlands is unknown, it was extensive, spanning steep and gentle slopes, and included elaborate canal systems for irrigation and runoff control (Denevan, 2001). Another indicator that the Inca significantly shaped their environment is evidence of managed afforestation to stabilize slopes (Chepstow-Lusty and Jonsson, 2000; Denevan, 2001). The use of fire for anthropogenic land use was widespread (Chepstow-Lusty et al., 1996).

Amazonia. A large body of research shows that locally intense landscape modifications occurred before Europeans arrived (e.g. Carson et al., 2014; Denevan, 2001; Erickson, 2010; Heckenberger et al., 2003; Mann, 2005; McMichael et al., 2012; Piperno et al., 2015). The alteration of forest by pre-Columbian groups in Amazonia ranged from the removal of non-useful plants and the protection of useful plants over simple seed dispersal and the relocation of seedlings to active fire management, the construction of earthen structures such as raised fields, ADE soils, homegardens and polyculture agroforestry (Heckenberger et al., 2003; Heckenberger and Neves, 2009; Levis et al., 2017, 2018; Maezumi et al., 2018a; Stahl, 2015; Woods et al., 2009). Modern forest composition in southwestern and eastern Amazonia (Levis et al., 2017; Ter Steege et al., 2013) and along floodplain habitats (De Souza et al., 2018; McMichael et al., 2017) are thought to be strongly influenced by pre-Columbian plant cultivation as well as incipient tree domestication, alongside crop domestication, including rice (Oryza sp.), cassava, peanut (Arachis hypogaea) and chili pepper (Capsicum sp.) (Hilbert et al., 2017; Watling et al., 2018)

Anthropogenic structures such as ADEs, earthen mounds (raised fields) and large (5 km) ring ditches, roads (Heckenberger et al., 2003) and fish weirs (Blatrix et al., 2018; Erickson, 2000) point to large (>1000 people) pre-Columbian settlements all over Amazonia including coastal French Guiana, the seasonal dry eastern Amazonia (Marajo Island and Santarem Region), central Amazonia’s Acre  province, the upper Xingu river in Brazil, the Llanos de Moxos, and the Beni Province of Bolivia (Carson et al., 2015; De Souza et al., 2018; Heckenberger et al., 2003; Piperno et al., 2015; Roosevelt et al., 1991; Stenborg et al., 2018; Whitney et al., 2014, Fig. 2). The dominant construction period for raised fields and ring ditches was between 500 CE and 1450 CE (Roosevelt, 2013). The continuity of the population history of these sites and their population density however remain uncertain and vary between sites (Denevan, 2014; Meggers, 2001). The workforce required to construct these structures remains debated as the method of construction is uncertain, ranging from frequent large scale burning (Carson et al., 2015; Whitney et al., 2014) to episodically clearing with stone axes and low severity burning for agroforestry-type systems (Clement et al., 2015; Erickson, 2010; Stahl, 2015). Alternatively, it has been suggested that a small workforce may has gradually constructed them very slowly over several thousand years on the margins of forest clearings (Watling et al., 2017).

North America. Pre-Columbian land use in the US and Canada varied from largely nomadic hunter-gatherers in the north towards semi-permanent and permanent agriculture further south (Abrams and Nowacki, 2008; Stinchcomb et al., 2011). Mississippian people grew maize, squash (Cucurbita pepo) and beans in alluvial areas and constructed large earthen mounds (Delcourt and Delcourt, 2004). Pre-Columbian societies in the southwestern North America built terraces and practised floodwater farming (Doolittle, 1992). In eastern North America people practised rainfed agriculture with sunflower (Helianthus annuus) and squash as staple crops (Smith, 2006), whilst altered forest composition through selective planting. Fire management was practised from the Midwest to the east coast (Abrams and Nowacki, 2008).

This gives an overview on how and to what degree indigenous people in the Americas used land, but it should be emphasized that not all sites with recorded human impact were populated at the time Europeans arrived. For example, in Mexico and Central America some parts of previously used cities and land were abandoned, or used at lower intensity, following widespread societal transformation in the Maya civilization (Beach et al., 2015; Guderjan et al., 2009; Tainter, 1988). Furthermore, sampling bias may be important particularly for Amazonia as ease of access to sites limits the sampling distribution across the region (McMichael et al., 2017).

3.2. New pre-Columbian per capita land use estimate

To calculate the area of pre-Columbian land use we need to know the typical land area used to provide sustenance for a person living in the Americas prior to 1492 CE, the per capita land use. Here we review per capita land use estimates for each of our seven regions (Table 2).

Caribbean. No direct per capita land use estimates are available for the Caribbean. Considering the similarity in agricultural systems, a similar per capita land use can be assumed to Mexico and Central America. The spatially explicit LUC dataset from Kaplan et al. (2011), KK10, applies a land use area per capita estimate of ~1.5 ha to Mexico, Central America and the Caribbean. The HYDE 3.1 LUC dataset instead provides a spatially explicit population dataset (Klein Goldewijk et al., 2011), and by dividing the land use area by the population, per capita land use was 1.2 ha in this region in 1492 CE. The third dataset, P08 (Pongratz et al., 2008a), assumes zero land use per capita in the Caribbean.

Mexico. The most common land use form at the time, the milpa system with its four-year cropping cycle, translates into 0.57e1 ha per capita (Borah and Cook, 1963; Bullock et al., 2013; Morfín and Storey, 2016; Turner and Brush, 1987). Estimates however vary substantially between sites, likely due to differing environmental conditions and family sizes (e.g. 1.5 ha per capita, Drucker and Heizer, 1960). KK10 assumes a land use per capita of ~1.5 ha per capita, about three times larger than that of HYDE 3.1, at 0.39 ha per capita. Unlike the other two LUC datasets P08 separates land use into lowland and highland land use. It applies a cropland per capita for highlands of 0.2 ha, based on Collins (1983), and for lowland land use just 0.05 ha per capita, which is one-tenth of the lowest milpa estimate.

Central America. Pre-Columbian land use systems in Central America are very similar to Mexico (Whitmore and Turner, 1992) and many studies do not distinguish between the two regions. For KK10 this region is assumed to have the same per capita land use as Mexico, 1.5 ha per capita, while HYDE 3.1 is slightly higher than for Mexico, at 0.43 ha per capita, while P08 assumes zero preColumbian land use anywhere in Central America (Fig. 3).

Inca Territory. Observations of the Aymara people have shown that contemporaneous indigenous cultivation in the Andes requires 0.2 ha per capita (Collins, 1983), although this is likely influenced by colonial and modern day practices and thus has probably a higher productivity than its pre-Columbian form (Boserup, 1965; Ruddiman and Ellis, 2009). Again, HYDE 3.1 has larger values, 1.4 ha per capita. KK10 has not published an estimate for this region nor provides a spatially explicit population dataset to establish a land use per capita estimate. P08 has low values, at 0.05 ha per capita in the lowlands (the same as P08 in Mexico) and every more efficient agriculture, at 0.2 ha per capita, in the highlands. The Inca Territory is the only region in the Americas where P08 assigns a third value, for pasture, at 0.57 ha per capita (for llamas and alpacas, the only domesticated grazers). The other two LUC datasets do not assign pasture to the pre-Columbian Americas, assuming no widespread use of domesticated grazers.

Amazonia. Highly productive ADEs, as found along riverine settings, make intensive semi-/permanent cultivation possible, with 0.2e0.3 ha per capita used (Woods et al., 2009), but such soils cover only ~3% of lowland Amazonia (McMichael et al., 2014). The larger part consists of less productive soils where extensive cultivation with long fallow (~15 years average, up to 30 years) was mostly employed to farm cassava. These long fallow systems require between 1.25 ha per capita (Clark and Uhl, 1987) and 2.2 ha per capita (Dull et al., 2010). Where shorter fallows are possible (e.g. Guiana, 4e6 years) estimates are lower: 1e1.3 ha per capita (McKey et al., 2010). Dull et al. (2010) report an average 1.2 ha per capita for lowland tropical forest regions of the Americas. Surprisingly, HYDE 3.1 gives a calculated rate of 0.9 ha per capita for Amazonia, a lower value than seen in either the Inca Territory or Mexico. Per capita values for KK10 are not derivable for Amazonia, while P08 give zero anthropogenic land use anywhere in Amazonia or its contiguous forested area before 1492 CE. North America. Published per capita land use estimates are not available for North America, but likely ranged from ~1 to 1.5 ha per capita for intensive slash and burn with short fallows (Doolittle, 1992) to several hectares per capita for large-scale fire-based forest clearing to support indigenous hunting strategies (Abrams and Nowacki, 2008; Anderson, 2006), based on similar farming systems in other regions in the Americas. For HYDE 3.1, we calculate 0.5 ha per capita for the US and Canada, more efficient than anywhere in Central or South America. P08 assume zero land use across the US and Canada, and the land use per capita estimates for KK10 are not derivable.

Rest of Americas. There are no direct estimates of per capita land use from Venezuela, Uruguay, Paraguay and Argentina. For HYDE 3.1 we calculate 1.9 ha per capita for Venezuela as well as Uruguay, 0.2 ha per capita for Paraguay and 0.4 ha per capita for Argentina, while again, P08 assumes zero land use anywhere in these four countries and land use per capita estimates for KK10 are not derivable.

Large: Table 2. Published land use per capita estimates (ha per capita) used for the different regions to calculate the pre-Columbian land use extent. Hemisphere is combined North and South America. Pongratz et al. (2008a) assumes different land uses for Highlands (H), Lowlands (L) and Pasture (P).

3.3. New pre-Columbian anthropogenic land use extent estimate

We group per capita estimates from 16 studies into the same regions as for the population estimate (section 2.2), with hemisphere estimates (Lewis and Maslin, 2015; Nevle and Bird, 2008; Ruddiman, 2007; Ruddiman and Ellis, 2009) applied to all regions (Table 2). This results in 63 regional estimates of per capita land use (Caribbean, n = 9; Mexico, n = 11; Central America, n = 9; Inca Territory, n = 9; Amazonia, n = 14; North America, n = 5; Rest of Americas, n = 6; Table 2), which we cross-multiply with the population estimates of the corresponding region to obtain every possible combination of pre-contact land use area estimates per region and sum them using the same methods described in section 2.2. Our median land use is 1.04 ha per capita (IQR 0.98-1.11 ha per capita), and the total extent of anthropogenic land use in 1492 CE is therefore 61.9 million ha (IQR 43.3e87.1 million ha; Table 1; Fig. 4B).

Our new per capita land use estimates are closest to Nevle and Bird (2008), 1 ha per capita, and the 1.3 ha estimated by Lewis and Maslin (2015). Our estimates are six times larger than the P08 value (0.18 ha per capita) and nearly double the HYDE 3.1 value (0.68 ha per capita). On the other extreme, Ruddiman’s (2007) figure is twice what we find (2 ha per capita), and the KK10 value is six times as large as our estimate, at 6.25 ha per capita. Our total anthropogenic land use area of 61.9 million ha in 1492 CE is very different from all three of the land-use database values utilized in modelling studies, P08, HYDE 3.1, and KK10. It is nine times as large as P08 (7.5 million ha) and a little over double the HYDE 3.1 value (27.7 million ha), but is six times smaller than the KK10 of 375 million ha (Table 1). Note that none of these land cover product values fall within the IQR derived from the empirical data, and diverge from each other on per capita land use (which varies thirty- five-fold) more than on the 1492 population (which vary by 50%). Only KK10 has agricultural activity that spans the same region as the literature shows, but employs a too large per capita land use (Fig. 3). We highlight that per capita land use estimates vary substantially on a regional scale and are not well constrained compared to population estimates. For example, per capita land use estimates in Mexico vary more than threefold despite being a well studied region. Our results show that pre-Columbian populations were large enough, and their per capita land use comparable to other regions of the world at that time (Europe, 0.5e1 ha per capita; China, ~0.5 ha per capita, Klein Goldewijk et al., 2017), to practice agriculture at a scale that transformed their environment to an extent that population changes after European arrival would have probably had a substantial impact on the environment.

4. The Great Dying

Accounts from contemporary eyewitnesses and documentary evidence from census data report a widespread collapse of the indigenous population over the decades after European contact (e.g. Cook, 1981, 1998; Diamond, 1997; Dobyns, 1993; Mann, 2005; Lewis and Maslin, 2015, 2018). Such a large population decline meant that the indigenous population would not have been able to manage their existing agricultural systems over such large areas. While other factors such as warfare, the enslavement of indigenous people, and hunger following social disintegration resulting from the loss of such a large fraction of societies meant even larger population losses (e.g. Cook, 1998; Mann, 2005; Dobyns, 1993; Lovell and Lutz, 1996), we focus on the epidemics as the main driver behind the majority of the deaths in the Great Dying.

4.1. Unprecedented mortality rates after European arrival

Existing evidence suggests that the indigenous population collapse was primarily caused by the introduction of pathogens unknown to the American continent (“virgin soil epidemics”) together with warfare and slavery (Black, 1992; Crosby, 1976; Dobyns, 1993; Joralemon, 1982; Mann, 2005; McNeill, 1977). Part of a wider Columbian Exchange of once-separate continental fauna and flora, these epidemics were introduced by European settlers and African slaves and were passed on to an indigenous population that had not been previously exposed to these pathogens and therefore did not initially possess suitable antibodies (Dobyns, 1993; Noymer, 2011; Walker et al., 2015). Such diseases included smallpox, measles, influenza, the bubonic plague, and later malaria, diphtheria, typhus and cholera. Most of these diseases originated from domesticated farm animals from Europe to which Native Americans had no prior exposure (Lewis and Maslin, 2018). The relative absence of American diseases arriving in Europe can therefore be explained by the low number of domesticated animals in the pre-contact Americas (Diamond, 1997; Lewis and Maslin, 2018; McNeill, 1977). Thus, influenza, smallpox, bubonic plague and other diseases ravaged the Americas, and not vice versa. Such diseases typically individually killed ~30% or more of the initial population. Hence a series of epidemics in rapid succession could have led to the loss of whole societies.

Overall, hemisphere wide post-epidemics population estimates range between 4.5 million and 14.4 million for 1600e1700 CE (Denevan, 1992a; Dobyns, 1966; McEvedy and Jones, 1977, Table 3). These studies are less clear on their assumptions than their reports of pre-contact population estimates. Furthermore, the rate of loss is strongly influenced by the chosen dates used to calculate it (Dobyns, 1966). Loss rates ranging from 40% (McEvedy and Jones, 1977) to 95% (Dobyns, 1966), with Denevan’s (1992a,b) more recent update producing a robust initial population estimate of 54 million (similar to our 60.5 million) that results in a 90% decline to 5.6 million in 1600 CE.

We review the published regional depopulation percentages and post-epidemics population estimates for each of our seven regions. Fig. 2 illustrates the spread of these diseases based the documented locations of disease outbreaks during the first century after European arrival.

Caribbean. Following the first disease outbreak in 1497 CE on the island of Hispaniola, the first census for the whole island was established in 1508 CE and arrived at 60,000 indigenous people (Moya Pons, 1992), a reduction of 85% within 16 years if Bartholomew Las Casas’ initial 4 million estimate was to be correct. A review of nadir populations, the estimated low point of a population, gives a population of 22,000 for the whole of the Caribbean by 1570 CE, corresponding to a 99% reduction from the 4 million estimate (Dobyns, 1966).

Mexcio. The first documented old World pathogens first reached the continent in southern Mexico (Cook, 1998). Here the most devastating outbreak occurred in 1520 CE when a single smallpox epidemic killed an estimated 30e50% of the indigenous population of Mexico (Cook and Borah, 1960; Dobyns, 1993; McCaa, 1995). Following the first comprehensive census in 1568 CE the population of central Mexico had already declined to 2.7 million (Sanders et al., 1979), which corresponds to a maximum decline of 87% within the first 50 years of European arrival when using a precontact population of 20 million (Whitmore and Turner, 1992). The nadir population of 1.5 million, equalling a total population decline of 93%, is reached in 1650 CE (Dobyns, 1966). Thus, a rapid population collapse of up 90% is plausible for the most populous parts of the country.

Central America. Combining the initial population estimates (315,000 to 2 million) and the 1550 CE population of 121,000e430,000, results in a 62e79% loss of Mayan population in Guatemala within 30 years of contact (Lovell and Lutz, 1995; Zamora, 1983). For all of Central America, Dobyns (1966) review on nadir populations gives a low point of 540,000 by 1650 CE, which corresponds to a depopulation of 90% based on pre-contact population of 6 million (Denevan, 1992a).

Inca Territory. Between 30% and 50% of indigenous populations in the Inca Territory reportedly died in the late 1520s as smallpox epidemics spread into the Inca kingdom (Cook, 1998). Some regions such as coastal Peru experienced only little initial decline (~10%, Shea, 1976), which supports the idea of waves of localized epidemics rather than subcontinent wide pandemics (Storey, 2012). By 1620 CE about 670,000 indigenous people remained in the Inca heartland, which is equivalent to a depopulation of 93% based on an initial population of 9 million (Cook, 1981).

Amazonia. Robust post-epidemics population estimates and depopulation percentages are sparse, with most studies applying a 90% population decline over a poorly defined post-contact period (Dobyns, 1966; Dull et al., 2010). A review on Amazonian archaeology concludes that in the Upper Xingu region, Brazil, an indigenous group (Xinguano) that occupied thousands of towns had been reduced to around 500 people by the time of the first comprehensive survey in the 1950s, corroborating a high depopulation rate estimate (Heckenberger and Neves, 2009).

North America. The first documented contact between indigenous groups and Europeans took place at the North American east coast in the 1530s when explorers reported dense Iroquois settlements. These settlements disappeared by the end of the 1600s (Sauer, 1980), 20 years before the first recorded smallpox epidemic that would kill 90% of the native population of New England (Davies, II, 2012). In the Jemez province, southwestern United States, an archaeological study found that a depopulation of 87% occurred within 100 years after contact (1541e1640 CE, Liebmann et al., 2016). By 1800 CE an estimated one million indigenous people were left in North America (Denevan, 1992a). Dobyns (1966) suggests the low point of 490,000 was reached in the 1930s, although it can be assumed that the initial waves of epidemics during the 1600s and 1700s were responsible for most of the depopulation.

Large: Table 3. Calculated depopulation percentages from published hemisphere-wide pre-contact population estimates and their corresponding posts-epidemics population estimates (2nd column) as well as relative to a standardized median 1600 CE population of 6 million people (3rd column).

Rest of the Americas. We found no depopulation estimates for Venezuela, Argentina, Uruguay and Paraguay

Only seven regional post-epidemic population estimates, or mortality rate estimates, have been published for 1600e1650 CE (Caribbean, n = 1; Mexico, n = 1; Central America, n = 1; Inca Territory, n = 2; Amazonia, n = 0; North America, n = 2; remainder, n = 0). It is therefore not possible to apply the same approach to estimating the population reduction as we used to calculate the pre-contact population and land use. Instead, we use the median of the estimates for the population of the Americas at 1600 CE, six million people (Denevan, 1992a; Dobyns, 1966; McEvedy and Jones, 1977; Klein Goldewijk et al., 2010; Krumhardt, 2010), and our calculated pre-contact population of 60.5 million (IQR 44.8-78.2 million) to estimate mortality. Based on these estimates 54.5 million people died (IQR 39.0-72.4 million; Fig. 4C), or 90% of the pre-contact population (IQR 87-92%). This is at the higher end of estimated losses (Table \tab:depopulationrate), but is consistent with those regions that are better studied (e.g. Central Mexico, Inca Territory) and when multiple waves of epidemics affecting the indigenous population, as well as warfare, slavery, and famine, are taken into account.

A mortality of up to 95% is higher than other virgin soil epidemics throughout history (Table 4, Crosby, 1972). While most of the other epidemics in history however were confined to a single pathogen and typically lasted for less than a decade, the Americas differed in that multiple pathogens caused multiple waves of virgin soil epidemics over more than a century. Those who survived influenza, may later have succumbed to smallpox, while those who survived both, may then have caught a later wave of measles. Hence, there were documented disease outbreaks in the Americas that killed 30% of the remaining indigenous population over 50 years after initial contact, i.e. between 1568 CE and 1605 CE (Black, 1992; Thornton, 1987). Indeed, high mortality rates from nonendemic diseases within indigenous communities were still common in the mid-1800s (Boyd, 1992). With at least eight major diseases documented (smallpox, influenza, measles, typhus, pneumonia, scarlet fever, malaria and yellow fever), a 30% mortality per new arrival e not uncommon for virgin soil epidemics e would result in a depopulation of 95% (Table 4). The depopulation is therefore much higher than might initially be expected from comparison with for example, the Black Death epidemic in Europe in the fourteenth century which was caused by a single pathogen (see Table 4).

For indigenous people that survived, immunity to most of the European diseases would be acquired during childhood of future generations (Black, 1992; Crosby, 1976; Dobyns, 1993; Jones, 2003; Noymer, 2011). Two main hypothesis, not mutually exclusive, have been proposed to explain why the depopulation continued until centuries after initial contact. The first is that the low genetic diversity between the indigenous hosts has facilitated the spread of potent pathogens (Black, 1992, 1994; Lewis, 2010; Lindenau et al., 2013; Salzano and Callegari-Jacques, 1988; Wang et al., 2007). The ancestors of the modern Native Americans migrated most likely from East Asia into North America (Goebel et al., 2008; MorenoMayar et al., 2018; Raghavan et al., 2015). Due to their small initial group size, the newly established population of the Americas had a lower level of genetic diversity compared to the original Asian population (“founder effect”) which would have lowered resistance against diseases from certain pathogens (Black, 1994; Lindenau et al., 2013). However, despite this there is no evidence for a causal relationship between genetic differences and the increased severity of the impacts of diseases (Walker et al., 2015).

A second “environmental hypothesis” (Walker et al., 2015) suggests that the effects of the European conquest and the consequences of the epidemics amplified the already devastating deadliness of the diseases (Crosby, 1976). Direct effects of the European conquest included warfare and the exploitation of indigenous populations (Livi-Bacci, 2008), including their forced relocation of people and slavery (Nevle and Bird, 2008). The consequences of famine from crop failure due to epidemic-driven labor shortages (Jones, 2003; Marr and Kiracofe, 2000) and lack of nursing care (Joralemon, 1982) led to a decline in birthrates following social stress and subsequently further population decline (Cook, 1998). Additionally, cultural factors such as no prior knowledge of quarantine procedures (Jones, 2003) and people moving along major trade routes and to areas more remote from direct European contact are likely to have contributed to the spread of epidemics. Thus, societal breakdown following ongoing war, slavery and famine due to shortage of healthy labor force likely increased the potency of epidemics in the Americas (Jones, 2003) with 90% of the population lost by 1600 CE (IQR 87-92%).

5. Extent of land use cessation and secondary vegetation succession

The decimation of Indigenous American populations by reoccurring waves of epidemics combined with societal disruption likely drove large reductions in agricultural activity. The absence of anthropogenic land use would mean that these areas would have undergone secondary succession. Using indigenous population decline estimates and their per capita land use we can estimate the extent to which vegetation succession in the Americas occurred following the Great Dying in the first century after European arrival, and how this may have impacted resulting carbon sequestration.

5.1. Secondary succession dynamics

The pathways of secondary succession will differ depending on the vegetation at abandonment, the proximity of sources of seeds, type of past anthropogenic land use, the soil, the climate, and exogenous disturbances such as fires (e.g. Chazdon et al., 2016; Griscom et al., 2017; Marín-Spiotta and Sharma, 2013; Poorter et al., 2016). However, in almost all cases an increase in plant biomass occurs over time, increasing the carbon stored on the land. For example, in an area where previously tropical forest has been cleared to plant annual crops microclimatic conditions are usually hotter and drier than suitable for most tropical forest species. When anthropogenic land use ceases, light-demanding, “pioneer” tree species that are more heat and moisture tolerant establish first, changing the microclimate and facilitating growth of typical tropical forest species (Guariguata and Ostertag, 2001; Rozendaal et al., 2017). Light availability in the understorey decreases as succession progresses, allowing the establishment of trees that will replace the pioneers, further increasing carbon stocks. However, in wetter forests, with their denser canopy than in drier forests, the transition from light-demanding pioneer species to later-successional forest species occurs at a faster rate compared to dry forests. Thus, while vegetation success increases carbon stocks and typically increases plant species diversity, rates depend on several factors, but typically, carbon stocks rise rapidly over the initial decades and then the rates of sequestration slows, with forest structure typically becoming similar to mature undisturbed systems within about 100 years (Martin et al., 2013), as does plant biodiversity (Chepstow Lusty et al., 2009; Platt Bradbury, 2000). However, structure and composition may never return to the prior vegetation state, only something somewhat structurally and functionally similar. One method of investigating successional changes following the Great Dying is to assess the paleo-ecological evidence.

Large: Table 4. Single virgin soil epidemics through history and their mortality.

5.2. Paleo-ecological evidence for widespread vegetation change

The widespread cessation of indigenous land use implies extensive secondary successional vegetation regrowth. This rapid change of land use is expected to cause fire use as a land management tool to decline, seen in sediment records as a reduction in charcoal. Secondary succession and implied increasing carbon stocks would be reflected in pollen records. We review the evidence for both in turn.

Charcoal. The use of fire as an agricultural tool for land clearance before the arrival of Europeans is well documented (Bush et al., 2008; Mayle and Power, 2008; Nevle et al., 2011), although there are exceptions such as fire-free agriculture documented French Guiana (Iriarte et al., 2012). Given the near-ubiquity of fire for agricultural production, any collapse in agricultural land use could be expected to be mirrored by a decreasing abundance of macrocharcoal in sediment cores near agricultural areas, record resolution permitting (Clark, 1988; Dull et al., 2010; Nevle and Bird, 2008; Nevle et al., 2011). Although it is possible that the prior land use could lead to forests with a differing structure and species composition from the pre-anthropogenic land use, meaning the new forests may have a differing natural fire regime than that prior to human impacts (e.g. Bush et al., 2007; Maezumi et al., 2018b). Furthermore, climatic changes associated with the drop in atmospheric CO2 and air temperatures may also affect rainfall and hence the probability of fires. Paleo-ecological reconstructions of biomass burning indicate a period of reduced fire activity in the Americas during the 1500s, although different studies attribute the reduction to different causes (Dull et al., 2010; Marlon et al., 2008, 2013; Nevle and Bird, 2008; Nevle et al., 2011, Fig. S1).

One group of studies (Dull et al., 2010; Nevle and Bird, 2008), examining charcoal records from sites across tropical Meso- and South America, attribute the decline in biomass burning to the indigenous population collapse. Another group of studies (Marlon et al., 2008, 2013; Power et al., 2013) argue that the driver was a change in the climate. The latter group suggest biomass burning was already in decline on a hemisphere-wide scale before Europeans arrived, following a shift towards wetter conditions that suppressed biomass burning globally. Of course, these two hypotheses are not mutually exclusive, the decline in charcoal deposits and biomass burning may be due to both climate and the human depopulation. Hence, it is possible that regional climate is driving a decline in biomass burning in the 1500s, but at a subset of sites where the human population also collapsed, this decline in biomass burning is being misattributed solely to climate. Alternatively, locations with an apparent loss of human population and concomitant decline in charcoal may be being solely attributed to people, rather than a mixture of direct human and climate forcing impacts. This attempt at defining one of two causal mechanisms is further complicated by the non-random location of study sites and that a decline in anthropogenic burning and increased carbon uptake from regrowing vegetation may amplify any climate-driven decline in burning via a feedback effect. The decline in biomass burning in the Americas is broadly consistent with the abandonment of agriculture following the Great Dying, although how much of this reduction is directly attributable the human population reduction is debated.

To analyse local fire histories to assess whether climate or humans are the cause for changes in the charcoal record we use our seven regions and select all records from the Global Charcoal Database that cover the period 1500e1600 CE at a resolution of <50 years, which is sufficient to capture potential population changes following the European contact epidemics (version 4.0.3, paleofire R-package, Blarquez et al., 2014). The sample sizes are: Mexico, n = 2; Central America, n = 4; Caribbean, n = 4; the Inca Territory, n = 4; Amazonia, n = 2; North America, n = 8 (Fig. 2). To compare the records we first normalize the biomass burning reconstruction to z-scores via Box-Cox transformation, where zi ¼ ðxi x01750Þ=sxð01750Þ, with x01750 being the charcoal value over the base period from year 0 until 1750 (i.e. pre-industrial) and sxð01750Þ the standard deviation over this period (Power et al., 2013). Of the selected 24 records, half (12) exhibit a persistent decline or low point in biomass burning during the 1500s (dating uncertainty 35e60 years), with most (10) located in Meso- and South America, and two in North America (see Fig. S1). Of the 12 records, biomass burning declines have been attributed to the population collapse following European epidemics in three case studies representing five records (Berrío et al., 2002; Chepstow-Lusty et al., 1998; Dull et al., 2010), one study notes that the record has an anthropogenic fire signature but does not attribute the decline in biomass burning to a cause (Bush et al., 2007), one study attributes changes in two records (of the same location) to climatic fluctuations (Kennedy et al., 2006), while for the other four no explanation on the cause of the decline was given (Carcaillet et al., 2006; HigueraGundy et al., 1999; Muller et al., 2003). All six records with attributed anthropogenic fire signatures are in regions with high preColumbian populations: Mexico (Dull et al., 2010), Central America, specifically Nicaragua (Dull et al., 2010); and the Andean highlands in Colombia (Berrío et al., 2002) and Peru (Bush et al., 2007; Chepstow-Lusty et al., 1998). Absence of human activity in two records from Hispaniola (Dominican Republic) led the authors conclude that climatic changes as the likely cause for the observed changes (Kennedy et al., 2006). The studies with records showing a decrease or low point in biomass burning during the 1500s tend to focus on the more distant past (Caribbean, North America, HigueraGundy et al., 1999; Muller et al., 2003) or on a different aspect of fire and landscape (North America, Carcaillet et al., 2006). Overall, at a series of sites with large pre-Contact populations charcoal records show a decline in biomass burning that authors attribute to rapid population declines after European contact. Our analysis shows that declining biomass burning in six locations throughout the Americas e 25% of records e was caused by the indigenous population decline, while at a further 25% of records showed a decline in burning not attributed to a local population decline (climate or unknown) with the remaining 50% of records show no discernible trend.

Pollen. From sedimentary pollen records, used to reconstruct the composition of local plant communities, it is possible to infer nearby land use and settlement history. For example, a change from presence to absence of maize or other cultivated plants, followed by an increase in successional plant species could indicate LUC due to the site being abandoned, either because of deaths or movement of the local population. A recent compilation of pollen records classifies 14 sites in Meso- and South America as influenced by human activity that cover the period 1500e1700 CE (Flantua et al., 2016, for locations of sites, see Fig. 2; for corresponding studies, see Table S4). At seven of these sites vegetation changes are attributed to depopulation following European contact, at four sites land use increased after European contact (Brugger et al., 2016; Bush et al., 2015; Goman and Byrne, 1998; Park et al., 2010), while at three sites abandonment occurred before 1492 CE (Whitney et al., 2013, 2014). In Central America, at Laguna Atezca site abandonment coincides with European arrival in 1519 CE (Conserva and Byrne, 2002), while in P atzcuaro basin the record indicates widespread deforestation until the Spanish arrived at 1541 CE (Platt Bradbury, 2000). In the Inca territory, a site in the southern Cauca Valley in Colombia shows reforestation occurred at around 1500 CE (± 25 years), which the authors link to depopulation at the time of Spanish arrival (Berrío et al., 2002). Further south near the city of Loja, land use ceased at two sites between 1600 CE (± 43 years) and 1660 CE (Niemann and Behling, 2010; Niemann et al., 2013). Near Cusco, in the Inca heartland, the important and widely planted tree Alnus, declined after 1500 CE (± 40 years) while successional species increase. After European arrival around 1530 CE these land cover changes became even more pronounced, with rapidly colonizing herbs and shrubs taking over previously maintained pastures, likely due to a reduction in anthropogenic land use (Chepstow-Lusty et al., 1998, 2009). At this site charcoal concentrations do not decline, possibly implying that the remaining population attempted to clear overgrown land for cultivation, but largely failing to keep the successional species at bay (ChepstowLusty et al., 2003). In northeast Bolivia, climate-driven expansion of dense canopy forest (1556 CE ±36 years) follows a period of extensive ring-ditch construction that previously kept the forest from encroaching into the savanna (Carson et al., 2014). Overall, 50% of sites (7 of 14) showed a change in vegetation that was attributed by the authors of the studies to a decline in land use following European contact after 1500 CE, while three sites (20%) show a decline in anthropogenic land use before 1492 CE and four sites (30%) show an increase in anthropogenic land use after European arrival. Thus, the pollen data tends to agree with the independent historical evidence that population declines led to vegetation succession at several locations throughout the Americas.

5.3. New estimate for the extent of abandoned agricultural area after the Great Dying

By multiplying the population that died in each region (regional population depopulation rate), cross-multiplied with the corresponding regional estimates of per capita land use we estimate 54.5 million people died, and a median area of 55.8 million ha of anthropogenic land use was abandoned and recovered to forest or another natural vegetation (IQR 39.0e78.4 million ha; Fig. 4D). This represents about 1% of the total landmass of the Americas being abandoned following European arrive and the spread of waves of pandemic diseases.

Our best estimate is closest to that of Lewis and Maslin (2015) at 65 million ha recovering and Nevle and Bird (2008) at 35e90 million ha. All other estimates lie outside our calculated interquartile range. Well above our estimated range is KK10 that gives a 250 million ha recovering, while Ruddiman (2007) and Faust et al. (2006) suggest 110 and 100 million ha respectively. Below our IQR are the other two LUC datasets: HYDE 3.1, assumes only 7 million ha recovering; while P08 reports the opposite pattern, with an increase in anthropogenic land use after European arrival, of 4 million ha (Table 1). The wide differences in published overall values are explained by the choice of low or high initial hemisphere population estimates combined with low or high mortality rates and low or high per capita land area requirements. At the low extreme, the P08 anthropogenic land area in 1600 CE is positive because an increase in pasture offsets the small recovery from a decline in cropland. In turn the cropland area decline is small because P08 assume a small initial indigenous population that is combined with a very low mortality rate and extremely low per capita land area estimates (Table 1). Indeed, in this model, there is an expansion in land use in northern Mexico, Central America and the southeastern coast of Brazil; Fig. 3). Similarly, HYDE 3.1 assumes a small initial population, a population decline smaller than historic documents suggest, and a relatively low per capita land use. At the high extreme, KK10 is driven by an extremely large per capita land use, some six times our central estimate, but is combined with initial population estimates and mortality rates that are similar to our central estimate. The strength of our approach is the absence of such covariance as all available combinations are evaluated, and suggests substantial LUC e the recovery of 56 million ha of land e followed the arrival of the Europeans in the Americas after 1492 CE.

6. Post-epidemic carbon uptake from secondary vegetation succession

The abandonment of 56 million ha of land in the 1500s following the deaths of 55 million people would have resulted in widespread vegetation succession coupled with a decrease in human driven fire activity. The new vegetation is likely to store substantially more carbon (thereby removing it from the atmosphere) over a relatively short timescale. For example, the above ground biomass (AGB) of tropical forest in Amazonia is 161 Mg C ha1 , fourteen times greater than maize cropland at 11 Mg C ha1 (Malhi et al., 2006; Yan et al., 2017). The recovery time of above and below ground biomass and carbon stocks is fast, typically 55e95 years depending on the system (Hughes et al., 1999; Martin et al., 2013; Poorter et al., 2016; Read and Lawrence, 2003).

6.1. Secondary vegetation succession carbon dynamics

While moisture, temperature, soil and disturbance events in- fluence the rate of biomass accumulation, increases are typically rapid in the first two decades after abandonment and slowing thereafter. Rates are typically lower in drier forests compared to wetter forests (Rozendaal et al., 2017), meaning mature tropical dry forests generally have a lower biomass than tropical wet forests (Poorter et al., 2016). Thus, continental-scale average biomass recovery rates may be improved by using biome or region-specific estimates.

A recent synthesis of Neotropical biomass data from forest inventory plots in secondary growth forests shows that fast AGB recovery is due to a high average net carbon uptake into live biomass, of 3.1 Mg C ha1 yr1 , with 75% of the expected final biomass accumulated in the first 20 years after abandonment, and a 90% recovery within 40e100 years (Poorter et al., 2016). Similarly, a synthesis of Mexican secondary tropical forest data finds that AGB carbon accumulates fastest during the first 25 years after abandonment (4.7±0.5 Mg C ha1 yr1 ) and slows down thereafter (to 3.0±0.7 Mg C ha1 yr1 , Orihuela-Belmonte et al., 2013). In dry forests carbon accumulation rates are around four times lower compared to wet forests (Poorter et al., 2016). Accumulation rates of lowland and montane tropical forests have been found to be similar (Fehse et al., 2002; Spracklen and Righelato, 2016). Recovering forests outside the tropics sequester substantially less carbon, just 0.7 Mg C ha1 yr1 for the U.S. (Williams et al., 2012). Differences in the rate of AGB accumulation appear to be mainly influenced by differing environmental conditions rather than the specific type of previous land use (Orihuela-Belmonte et al., 2013; Poorter et al., 2016).

Soil carbon changes follow a similar pattern: carbon stocks increase at a rate of 1.0 ±0.1 Mg C ha1 yr1 over the first 20e30 years before slowing down (Don et al., 2011; Marín-Spiotta and Sharma, 2013), although some studies report substantial increases in soil carbon until 70 years after abandonment (Pregitzer and Euskirchen, 2004). A review of soil carbon stores for low-to mid-latitudes found that reverting to forest from cropland increases carbon in the upper 20 cm of soil by 0.5 Mg C ha1 year1 (Shi et al., 2013). In areas where the natural vegetation is not forest, studies of the conversion of cropland to grassland provide a similar scenario to the abandonment of fields. In these studies, soil carbon increased by 1.2 Mg C ha1 yr1 (Don et al., 2011). Taken together, both vegetation recovery and soil carbon uptake following agricultural abandonment imply substantial carbon uptake in the 1500s across the estimated 56 million ha abandoned farmland. 6.2.

6.2. New estimate of carbon uptake from vegetation succession

To calculate the total carbon uptake following anthropogenic land use abandonment and resulting natural vegetation reestablishment we multiplied the increase in carbon stocks (per unit area per unit time) and the regrowth area estimate. We first compiled above ground (n ¼ 112) and soil carbon uptake (n ¼ 38) rates following abandonment within each of our seven regions (Table 5). From these studies that predominantly report short-term (<20 years) uptake rates, we then assume that these apply for the first 20 years, and for longer-term uptake, we assume this is at 8% of the short-term value over the next 80 years, following OrihuelaBelmonte et al. (2013). The 100-year uptake rate is then calculated as uptake rate 20 years þ uptake rate 80 years 0:08. We assign the carbon uptake rates to our seven regions and crosssum the aboveground carbon uptake rates (Caribbean, n = 18; Mexico, n = 23; Central America, n = 23; Inca Territory, n = 6; Amazonia, n = 18; North America, n = 1, Remainder regions, n = 23) and soil carbon uptake rates (Caribbean, n = 6; Mexico, n = 7; Central America, n = 7; Inca Territory, n = 1; Amazonia, n = 6; North America, n = 4; remainder, n = 7; Table 5) to calculate all possible combinations of net total carbon uptake rates per region. We cross-multiply these with the previously calculated abandoned areas for each region to calculate a frequency distribution of the total carbon sink (Fig. 4E).

The median carbon uptake is 7.4 Pg C (IQR 4.9-10.8 Pg C, 1 Pg = 1 1015 g) from the abandonment of 55.8 million ha of land in the Americas over the 100 years after 1517 CE (date of first documented mainland epidemic). A major source of uncertainty is the land that would have been under various stages of fallow states meaning their continued succession may have led to a more modest uptake of carbon than we estimate here. Furthermore, the large uncertainty in the per capita land use estimates also means that relatively small changes in this parameter result in substantial changes in abandoned land area and terrestrial carbon uptake. Our new estimate is in the lower end of the range of published estimates using area and median carbon uptake values, which range from 7 Pg C to over 20 Pg C a century after the arrival of Europeans and their diseases in the Americas (Table 6). Despite this, our estimate confirms other studies suggesting a globally significant uptake of carbon into the land surface at this time (Faust et al., 2006; Lewis and Maslin, 2015; Nevle and Bird, 2008; Ruddiman, 2007). Our results are most closely in line with two published studies: Nevle and Bird (2008) and Lewis and Maslin (2015), whose range both include 7.4 Pg C uptake (Table 6).

By contrast, the modelling results based on KK10, HYDE 3.1 and P08 land use change estimates are all outside the IQR uptake we estimate (4.9e10.8 Pg C). KK10 reports carbon uptake of 40 Pg C, well above our value, while both HYDE 3.1 and P08, give near zero uptake, and so are well below our central estimate and the lower limit of our IQR (Table 6). This is because only in KK10 is LUC of large enough magnitude, 248 million ha, due to the large land use per capita values (Kaplan et al., 2011, Table 1). Using HYDE 3.1 with the same model, Kaplan et al. (2011) did not find an observable perturbation of carbon stock due to the small area of abandoned land, at 6.9 million ha, ultimately due to a much lower depopulation estimate and lower per capita land-use rates than we find (Table 1). The zero carbon uptake in vegetation in Pongratz et al. (2011) is substantially smaller than our estimate. This is firstly because the low depopulation estimate and very small per capita land use combine to result in only to a small abandoned area of cropland. The carbon uptake from the LUC is then almost completely offset by an assumed legacy carbon emissions from the soil. This modelled result is the opposite of that found in metaanalyses of field studies, which report an increase in soil carbon almost immediately after agricultural abandonment in tropical forest regions (Don et al., 2011; Marín-Spiotta and Sharma, 2013), and is not seen in other model studies (Kaplan et al., 2011). Furthermore, if soil carbon release were to fully offset carbon stocks, forest reforestation would not work as a climate mitigation strategy. Measurements, monitoring and modelling shows that it does (e.g. Arora and Montenegro, 2011; Griscom et al., 2017; Smith et al., 2014).

While we show that there is substantial carbon uptake following the abandonment of agriculture in the Americas following European arrival, some 7.4 Pg C over the 100 years from 1517 CE, the two modelling studies covering this period do not replicate the results we find from our data synthesis. This is due to combinations of three different primary reasons (the size of the population loss, per capita land use, and legacy soil carbon fluxes), although KK10, if it were parametrized with lower per capita land use assumptions, would likely to be the closest to our assessment of the published literature.

7. Post-epidemic land use change impacts and the 16th century carbon cycle

Our estimate of a 7.4 Pg C increase in the carbon stored on land following the Great Dying, i.e. in the 100 years after 1517 CE, corresponds to a decline in atmospheric CO2 of approximately 3.5 ppm (Fig. 4E; IQR 2.3e5.1), where 2.12 Pg C equals 1 ppm CO2 over decadal timescales (Trudinger et al., 2002). The magnitude and timing of this increase in carbon stocks on the land surface is consistent with an abrupt 7e10 ppm decrease in atmospheric CO2 (14.8e21.2 Pg C) beginning in the 1500s and ending in the early 1600s, shown in two high resolution ice cores from Antarctica, the Law Dome and the West Antarctic Ice Sheet (WAIS) ice cores (Fig. 1, Ahn et al., 2012; MacFarling Meure et al., 2006). Thus a simple calculation implies that 35e50% of the 14.8e21.2 Pg C reduction can be attributed to the Great Dying, if all other radiative forcing agents and carbon sinks and sources are assumed to have remained constant.

Large: Table 5. Published above ground (A) and soil carbon (B) uptake rates for the different regions in our analysis and their means.
Large: Table 6. Published estimates of carbon uptake following conquest driven reforestation, including estimates from this study with interquartile range (IQR).

However, without corroborating evidence such analyses could merely be quantifying an unusual coincidence. Furthermore, other radiative forcing agents may have changed over this time period, and carbon sinks and sources will not remain constant given the land-use change perturbation to the carbon cycle following the abandonment of 56 million ha of anthropogenic land. Ideally, fully coupled Earth System Model simulations that alter the known forcing agents in factorial-type experiments, including (i) the atmospheric CO2 decline, (ii) changes in total solar irradiation, (ii) the impact of volcanic eruptions, and (iv) anthropogenic land use change, in the Americas and the rest of the world, are required to precisely attribute the changes to the Earth system and fully understand the anthropogenic and non-anthropogenic carbon cycle impacts over the 1500s. Such analyses are not available, therefore this section presents the available evidence, firstly to show that the observed atmospheric CO2 decline is largely caused by a large terrestrial sink, coinciding with the enhanced carbon uptake from LUC in the Americas and secondly, to show that the increased terrestrial sink is not primarily caused by other changes at the time, including total solar irradiation (TSI) or volcanic eruptions. Finally, we estimate the major carbon fluxes in the 1500s to balance the global carbon budget at this time to the assess the impact of the Great Dying on the Earth system. Such an investigation has assumed more importance recently, with the impact of the Columbian Exchange, including the Great Dying impacts on atmospheric CO2 concentrations, being proposed as the beginning of a new geological epoch dominated by human actions, termed the Anthropocene (Lewis and Maslin, 2015, 2018).

7.1. Land carbon uptake and the 1610 decline in atmospheric CO2

The Law Dome record shows a steep decrease in atmospheric CO2 of 10 ppm, beginning moderately at around 1530 CE, and intensifying from ~1570 CE with a clear low-point of 272 ppm at 1610 CE, followed by a small increase until 1650 CE, when consistently low levels are maintained until 1750 CE, the beginning of the Industrial Revolution (MacFarling Meure et al., 2006). The WAIS record also shows a rapid decrease in CO2 beginning around 1570 CE, followed after 1600 CE by a slower decrease to a minimum at 1650 CE, with low levels remaining until the beginning of the Industrial Revolution around 1750 CE (Ahn et al., 2012). While both show a similar 16th-century decline, more clearly seen when the records are plotted as anomalies, as in Fig. 5A, there are differences in magnitude of the two absolute minima, ~3 ppm, and the timing of that absolute minima, 1610 CE or 1650 CE. These differences are thought to be due to dating uncertainty, local CO2 variation and imprecise experimental uncertainty estimation (Ahn et al., 2012). Both records show an unusual decline and the lowest CO2 concentrations in the last two millennia (Ahn et al., 2012)

The isotopic signature of the CO2 in the ice cores shows this reduction was driven by an increase in terrestrial carbon uptake, indicated by a preferential uptake of 12C over 13C through photosynthesis, leading to 13C atmospheric CO2 enrichment (Bauska et al., 2015). Specifically, the increase in d13CeCO2, beginning after ~1520 CE, shows a near-constant rate of terrestrial carbon uptake until around 1610 CE and a cessation of the sink thereafter, most consistent with the pattern of CO2 reduction seen in the Law Dome record (Fig. 1). Additionally, d13C data from tropical shallow water sclerosponges also indicates that the decline in atmospheric CO2 is caused by an increase in the terrestrial carbon uptake (Bohm € et al., 2002), as do changes in carbonyl sulfide (COS) captured in Antarctic ice cores that are directly related to changes in photosynthesis (Rubino et al., 2016). Furthermore, a series of modelling studies also point to a carbon uptake on land at this time (Bauska et al., 2015; Joos et al., 1999; Rubino et al., 2016; Stocker et al., 2011; Trudinger et al., 2002). The mean of these model studies is an terrestrial carbon uptake of 19.9 Pg C (Table S5) Thus, the overall temporal pattern and magnitude of both the decline in CO2, changes in d13CeCO2, and data-driven models of the carbon cycle are consistent with uptake into terrestrial vegetation and soils. This unusual event of a rapid and large increase in terrestrial carbon stocks is consistent with a role for secondary succession following epidemics in the Americas after the arrival of Europeans, although there are other potential interpretations in the literature which we discuss in the next section.

7.2. Explaining the 1610 CE CO2 decline and land carbon uptake

To more fully understand the changes to the Earth system in the 1500s and 1600s we compile data from 1000 CE to 1800 CE for (A) the atmospheric CO2 concentration, (B) global surface air temperatures, and three potential drivers of these changes, (C) Total Solar Irradiation, (D) volcanic eruptions, and land carbon fluxes, as (E) the carbon flux to the land, (F) LUC in the Americas, and (G) LUC in the rest of the world (Fig. 5). Beyond visual comparisons that imply a strong role for LUC in the Americas, we can also assess the contribution of each of the potential drivers of changes in atmospheric CO2 and global surface air temperature by calculating their radiative forcing impact on air temperatures.

The lower atmospheric CO2 coincides with lower surface air temperatures (Fig. 5A). A global synthesis of over 500 paleoclimate records shows that 1577e1694 CE is the only period of significant global cooling within the past two millennia, and the only part of the Little Ice Age (LIA) that was global in extent (Fig. 5B from, Neukom et al., 2014). This analysis shows that the cooling over this period was 0.15+C. How does this compoare with the radiative forcing impact of the 3.5 ppm decline in atmospheric CO2 at this time? This would be expected to reduce radiative forcing by 0.07 W m2 (i.e. 5.35 W m2 lnð274:5 ppm CO/278 ppm CO2); 5.35 W m2 is the radiative forcing constant of CO2, Myhre et al., 1998). Scaling this decline in radiative forcing by a climate sensitivity parameter of 0.3e1.2 K [W m2 ] 1 (Myhre et al., 2013), suggests that LUC in the Americas over the 1500s and early 1600s is directly responsible for reducing global surface air temperature of 0.03e0.08+C. A second scaling approach gives a similar answer: using a fully coupled Earth System Model simulation reforesting the tropics in the 100 years after 2000 CE (Arora and Montenegro, 2011), but scaled linearly to our 56 million ha reforestation area, indicates a slightly greater reduction in temperature of 0.04-0.13+C. This suggest that 20-87%, with a mid-point of about 50%, of the 0.15+C reduction in temperatures can be attributed to the radiative forcing from the reduction in atmospheric CO2 resulting from the Great Dying.

The changes in TSI, shown in Fig. 5C appear unlikely to be the key driver of the 1610 CE CO2 decline and subsequent temperature reduction, despite such suggestions in the literature (Joos et al., 1999; Pongratz et al., 2008a; Trudinger et al., 2002). For TSI the Sporer Minimum (1416 € e1537 CE) and the Maunder Minimum (1645e1715 CE) are commonly suggested to explain the cool periods during LIA conditions seen in Europe, but these events leave a gap between the mid-1500s to mid-1600s, and so are not the cause of the increased land carbon stocks or the decline in temperature during this period (Schurer et al., 2013; Owens et al., 2017). Furthermore, TSI decreased substantially over the 1300s (Vieira et al., 2011), when no major change in land carbon uptake occurred (Bauska et al., 2015), suggesting there two variables are not well coupled (Fig. 5C cf 5E). While there was no long-term increase in TSI during the period of increasing terrestrial carbon stocks, interannual changes in TSI may have contributed to the CO2 and temperature decline, which has the potential to impact photosynthesis and terrestrial carbon uptake. For the period 1580e1610 TSI was 0.11e0.31 W m2 above the long-term 1000e1800 CE mean, thus of opposite sign to that necessary to explain lower temperatures, while at 1610 CE itself, TSI was 0.01e0.12 W m2 below 1000e1800 CE mean, suggesting some contribution to decreased temperatures at this time. Overall, the documented changes in TSI cannot explain the anomalously low CO2 levels over the late 1500s seen in the ice core records, nor lower global temperatures in the 1600s as seen in the paleoclimate temperature reconstructions, nor the larger increase in carbon stored on land (Schurer et al., 2013; Owens et al., 2017).

The impacts of aerosols from large volcanic eruptions, seen in Fig. 5D, can also be ruled out as the dominant cause of carbon uptake on land from 1520 CE to 1610 CE, despite such suggestions in the literature (Joos et al., 1999; Pongratz et al., 2011; Trudinger et al., 2002). Of the forty largest eruptions in the past 2500 years, only one occurred in the 1520e1610 CE period (Huaynaputina, Peru, in 1601 CE), whereas two occurred in the century prior (1453 CE, 1458 CE, both Kuwae, Vanuatu) and two in the century after it (1641 CE, Mt Parker, Philippines; 1695 CE, unknown location, Sigl et al., 2015). For documented tropical eruptions, which have the potential to impact the Earth system for longer, there were five during 1520e1610 CE compared to four during the following century; while only five extratropical eruptions occurred during 1520e1610 CE compared to 14 over the following century (Sigl et al., 2015). Probably the most similar 100 years in terms of total volcanic eruptions to 1510e1610 CE was 1610e1710 CE (Fig. 5D), which did not show a large increase terrestrial carbon uptake (Fig. 5E). Volcanic eruptions typically cause global temperatures to drop for a minimum of two to a maximum of ten years (Frolicher € et al., 2011; Sigl et al., 2015). Such short-term change of the expected radiative forcing does not reflect the temporal pattern of land carbon stock increase inferred from the isotopic signature of Antarctic ice-core CO2 records (Fig. 1) or double-deconvolution method (Fig. 5E). Overall, modelling studies investigating the impacts of volcanic eruptions suggest that they decreased surface air temperatures modestly over the period of the long LIA that affected the northern hemisphere between 1440 CE and 1920 CE, but are unlikely to be the primary cause of the increase in land carbon uptake over the 1500s (Schurer et al., 2013; Owens et al., 2017). While the case for the importance of LUC in the Americas is increasingly recognized there has been a suggestion that LUC outside of the Americas over 1500 CE to 1600 CE could mask the uptake of LUC in the Americas. This would lead to no detectable impact of LUC globally on the Earth System, and hence no impact on atmospheric CO2 levels or surface air temperature (Pongratz et al., 2011). Such a scenario, while possible, would not be compatible with the diverse evidence that there was net terrestrial carbon uptake at that time somewhere on Earth.

LUC emissions in the rest of the world would, however, need to be lower than the LUC uptake in the Americas for there to be a net sink, assuming all other factors remain constant. Extracting the LUC in Europe and Asia (the most populous and well-studied regions) from the three global LUC datasets, HYDE 3.1, KK10 and P08, each show that LUC emissions continued after 1500 CE but were considerably smaller (1.6±1.7 Pg C) than the carbon uptake in the Americas (7.4 Pg C, Fig. 5G). This suggests an overall terrestrial carbon sink from 1500 CE to 1650 CE, consistent with the diverse evidence of a land sink at this time. Overall, of the potential drivers of the decline in atmospheric CO2 in the 1500s and early 1600s, and the resulting decrease in global surface air temperature, the temporal pattern of LUC in the Americas fits best, rather than changes in the radiative forcing from TSI or volcanic eruptions, or the carbon cycle impact of LUC outside of the Americas (Fig. 5).

Fig. 5. Impact of changes in radiative forcing on atmospheric CO2 and temperature 1000–1800 CE. A) CO2 concentrations recorded in two Antarctic ice cores: Law Dome (grey, MacFarling Meure et al., 2006) and West Antarctic Ice Sheet (WAIS) Divide (blue, Ahn et al., 2012); B) Global mean temperature reconstruction from northern and southern hemisphere proxies, anomaly compared to the means from 1000 to 2000 CE (Neukom et al., 2014, grey, smoothed in black); C) Total Solar Irradiation anomaly compared to the means from 1000 to 2000 CE from Steinhilber et al. (2009, grey) and Vieira et al. (2011, green); D) volcanic, radiative forcing (grey, individual eruptions; black, smoothed, Sigl et al., 2015); E) Land-atmosphere flux based on double deconvolution of WAIS CO2 data (Bauska et al., 2015); F) Land use change (LUC) in the Americas from the LUC reconstructions KK10 (purple, Kaplan et al., 2011), HYDE 3.1 (blue, Klein Goldewijk et al., 2011), and P08 (grey, Pongratz et al., 2008a). The 1600 CE value has been replaced with the LUC estimate from this review (solid lines) and the original LUC reconstruction is plotted for comparison (dashed lines); G) LUC in the two other regions with considerable agrarian societies at the time, Asia and Europe based on the three LUC reconstructions. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

An alternative to assessing the radiative forcing impact of the LUC following the Great Dying, suggested in Ruddiman (2013), is to consider the surface air temperature decline over the 1500s and early 1600s, ~0.15+C and estimate the expected change in atmospheric CO2 from such a temperature decline. Therefore, any additional CO2 reduction that is not associated directly with the reduction in temperature, is assumed to be associated with LUC uptake of carbon. The coupled Climate-Carbon Cycle Model Intercomparison Project (C4MIP) models show that the temperature land surface feedback, termed gTL, is 40 to 70 Pg C K1 , meaning as temperatures warm, carbon is released from the land to the atmosphere (Cox et al., 2013). Mechanistically, this is because both gross primary productivity (GPP), i.e. photosynthesis, and autotrophic (plant) respiration (Ra), tend to decrease in cooler conditions, as does a similar set of processes in soils. But the response of Ra to cooling is less relative to the response of GPP, meaning overall net primary productivity increases and thereby carbon stocks on land. Hence a 0.15 K reduction in temperature equals an increased terrestrial uptake of 8.3 (±2.3) Pg C (mean from 40 0:15 to 0:70 0:15), or ~4 ppm CO2 over the 1500s and early 1600s. Thus, all else being equal, of the 7e10 ppm reduction in atmospheric CO2 seen in the ice core records, we would expect that the vegetation regeneration following the Great Dying accounts for the remainder, some 3e6 ppm CO2. This is consistent with the 3.5 ppm CO2 uptake we find in this study.

However, temperature changes are not the only change at this time. In addition to gTL changes in atmospheric CO2 also affect the carbon stored on land. The CO2 land surface feedback, termed bTL, is estimated to be 0.5e1 Pg C ppm1 CO2, meaning as CO2 increases so does land carbon storage (Cox et al., 2013). Mechanistically, this is because CO2 is a substrate for photosynthesis, thus lower CO2 means less photosynthesis, lower NPP, and lower carbon storage on land. Thus, this feedback reduces land uptake over the 1500s by 3.5e10 Pg C.

Combining both feedbacks (gTL þ bTL) from the decrease in surface temperature and CO2 over the 1500s, adds ~8.3 Pg C to the land due to lower surface temperatures, that is offset by ~6.4 Pg C (mean bTL 0.75 Pg C ppm1 a mid-value of 8.5 ppm decline) due to lower CO2, leaving a net uptake into the land surface of 1.9 Pg C or about 1 ppm CO2 reduction solely from lower temperatures and CO2 levels (Fig. 6). This level of terrestrial carbon uptake from the temperature and CO2 changes accounts only for 9e13% of the 7e10 ppm reduction in atmospheric CO2 and implies a dominant role for LUC in balancing the global carbon budget for the 1500s. Furthermore, by adding the net change of 1.9 Pg C from the temperature and CO2 feedbacks to our new estimate of LUC carbon uptake in the Americas, 7.4 Pg C, gives a total of 9.3 Pg C additional uptake into land in the 1500s, or 4.4 ppm CO2. Thus, using this method, some 45e65% of the 7e10 ppm decline in atmospheric CO2 can be explained by LUC in the Americas, after accounting for the Earth system feedbacks from the change in temperature and CO2 on the land surface at this time.

Through multiple routes we arrive at the conclusion that LUC in the Americas played an important role in driving lower atmospheric CO2 in the late 1500s and early 1600s. This is counter to conclusions reached by previous modelling studies that have not highlighted LUC as a dominant driver of terrestrial carbon uptake in the 1500s. Some of these are older publications (Joos et al., 1999; Trudinger et al., 2002), that may have not highlighted this factor because depopulation of the Americas has only more recently become well-known and incorporated into LUC datasets (Kaplan et al., 2011; Klein Goldewijk et al., 2011; Pongratz et al., 2008a,b). Other more recent studies that include LUC after 1500 CE ague that it was not sufficient to drive changes in atmospheric CO2 (Pongratz et al., 2011). This is partly a result of using old 1960s depopulation data combined with the assumption of extremely low per capita estimates of land use and partly due to legacy carbon fluxes that are not supported by observational studies (see section 6.2).

Fig. 6. Carbon fluxes from 1500 to 1600 CE relative to those from 1400 to 1500 CE. The first 3 bars, atmosphere, land-atmosphere and ocean-atmosphere fluxes are the mean from double deconvolution approaches (Bauska et al., 2015; Joos et al., 1999; Rubino et al., 2016; Trudinger et al., 2002) and a GCM (Pongratz et al., 2011), the fourth bar is the land-atmosphere flux estimated in this study, which is the sum of the four last bars, LUC (Americas), terrestrial carbon fluxes due to temperature (γTL), CO2 feedbacks (βTL), and LUC elsewhere. LUC (Americas) is calculated in this study; γTL and βTL are based on feedback strength from C4MIP models reported in Cox et al. (2013); LUC (elsewhere) is the average of KK10 (Kaplan et al., 2011), HYDE 3.1 (Klein Goldewijk et al., 2011), and P08 (Pongratz et al., 2008a) datasets Europe and Asia, the most populous and well-studied regions outside the Americas. Error bars are the min-max

Most recently, Rubino et al. (2016) argued that lower surface air temperatures drove terrestrial carbon uptake, rather than an increase in secondary forest succession. They analysed carbonyl sul- fide (COS) concentrations from ice-cores, which track GPP, demonstrating that synchronous GPP and Ra decline lead to an increase in land carbon uptake between 1520 CE and 1610 CE. Despite their conclusion, this result corroborates with our interpretation of the temperature response of the global carbon cycle at this time. Indeed, Rubino et al. (2016) explicitly estimate the land surface temperature feedback, gTL, as between 10 and 90 Pg C K1 , based on the decline in GPP, which is comparable to the C4MIP range used in this study (40 to 70 Pg C K1 ).

However, Rubino et al. (2016) then conclude that vegetation recovery from LUC is incompatible with the declining GPP, and so cannot be the cause of the change in land carbon stocks. Yet observations of terrestrial carbon stocks do not support such conclusion. The carbon stored in vegetation is determined by both the carbon input into a system e the net primary productivity, NPP, i.e. GPP minus Ra e and how long a unit of fixed carbon stays in the system, the carbon residence time (CRT). It is the CRT that lengthens during vegetation succession, from just 1e6 years in agricultural crops to ~50 years in many forests (Galbraith et al., 2013). Critically, GPP and NPP do not change radically over time as vegetation regenerates following cropland abandonment, yet carbon stocks rapidly increase (Anderson-Teixeira et al., 2012). Small changes in CRT are well known to have large impacts on vegetation carbon stocks, despite not typically being included in models (Ahlstrom et al., 2015 € ; Friend et al., 2014; Galbraith et al., 2013). Thus the assumption made by Rubino et al. (2016) that vegetation succession and declining GPP are incompatible is incorrect (succession continues even if temperature declines and GPP and NPP are lower). Furthermore, as COS has no known or theorized relationship to CRT, COS cannot be used in a straightforward manner to evaluate LUC impacts on land carbon stocks. Thus, lower GPP due to lower temperatures, and an increase in carbon stocks following LUC in the Americas, are not the mutually exclusive mechanisms described by Rubino et al. (2016).

Overall, the change in a forcing agent (a large abrupt LUC) driving a reduction of atmospheric CO2 (widespread vegetation regrowth), its location (carbon uptake on land), its timing (1520e1610 CE), and its magnitude (7e10 ppm CO2) each fit with the evidence for the impact of the arrival of Europeans in the Americas directly. By summing the impacts about half of the reduction in atmospheric CO2 can be attributed to LUC following the Great Dying (4.4 ppm of 7e10 ppm CO2 reduction).

7.3. The 16th century carbon budget

Given that the atmosphere, oceans and land surfaces are linked, a reconstruction of the global carbon budget of the century between 1510 CE and 1610 CE can provide a more comprehensive understanding of the role of LUC in the Americas on the Earth System at this time. Indeed, assessing the changes in atmospheric and oceanic carbon fluxes provides an independent assessment of the necessary uptake on land required to balance the global carbon budget. For the centuries prior to the 1500s, there was no change in atmospheric CO2 concentration (Fig. 5A), thus the carbon released from global LUC was offset by enhanced carbon uptake somewhere else in the Earth system, most likely in peatlands (Bauska et al., 2015; Ruddiman, 2007; Stocker et al., 2017). We therefore need to assess the change in the global carbon budget in the 1500s compared to the 1400s, which is detailed in Fig. 6.

The decline in atmospheric carbon, 14.8-22.2 Pg C, is well known (Fig. 6). Whether the oceans acted as a source or sink of carbon to the atmosphere relative to the 1400s is less clear (see Table S5). A number of double deconvolution studies and one GCM study have estimated the ocean response to the carbon uptake on land, giving a mean carbon source to the atmosphere of 5.5 Pg C relative to the 1400s (8:8 to 26.5 Pg C range, Bauska et al., 2015; Joos et al., 1999; Pongratz et al., 2011; Rubino et al., 2016; Trudinger et al., 1999, 2002, Fig. 6). While it should be noted that the results range from a sink to a large source, a strong ocean carbon sink over that period is not supported by the increase in d13CeCO2 (see section 7.1), thus a carbon source from oceans to atmosphere appears much more likely (Fig. 6). Such a carbon source from ocean to atmosphere requires an even larger carbon uptake on land to balance the carbon budget than the atmospheric CO2 decline alone suggests. Given the known decline in atmospheric carbon, 14.8e21.2 Pg C, and the addition of carbon from the oceans to the atmosphere (5.5 Pg C), the necessary mean land sink to balance the global carbon budget is approximately 23.5 Pg C over the 1500s (mid-range atmospheric carbon decline, 18 Pg C plus the ocean source).

A mean land sink of 23.5 Pg C being needed to balance the carbon budget of the 1500s is considerably larger than the carbon uptake from LUC in the Americas we compute above, at either 7.4 Pg C, or 9.3 Pg C uptake on land including the temperature and CO2 feedbacks. However, these values are an underestimate. From the 1400s to the 1500s LUC in the Americas shifted from a source of 2.5 ± 2.3 Pg C (mean of mid-ranges of KK10, HYDE3.1, and P08 land use of 19.2nbsp;± 17.7 million ha in the 1400s, assuming the 0.13 Pg C per 1 million ha derived in this study), to a sink of 7.4 Pg C following the Great Dying. Hence, a net carbon uptake from LUC in the Americas of 9.9 (5.1e15.6) Pg C in the 1500s compared to the 1400s. However, carbon was also lost from LUC elsewhere in the world. This added a net 1.6 (±1.7) Pg C to the atmosphere from the 1400s to the 1500s (mid-range KK10, HYDE 3.1, and P08 LUC in the most populous and well-studied regions, Europe and Asia, an increase from 64nbsp;± 7.1 million ha anthropogenic land in the 1400s to 76.1nbsp;± 20.5 million ha in the 1500s, a difference of 12.1nbsp;± 13.4 million ha, assuming of the same 0.13 Pg C per 1 million ha).

Thus the total LUC carbon storage on land in the 1500s compared to the 1400s is (i) the net uptake from LUC in the Americas of 9.9 Pg C (IQR 5.1e15.6 Pg C), (ii) a net release from LUC from the rest of the world of 1.6nbsp;± 1.7 Pg C, (iii) the land temperature feedback response increasing carbon storage on land by 8.3nbsp;± 2.3 Pg C, (iv) the land CO2 feedback leading to a release of 6.8nbsp;± 3.3 Pg C (Fig. 6). Combined, these fluxes lead to a net increase in carbon stocks on the land surface in the 1500s, relative to the 1400s, of 9.8 Pg C (ranging from a small source to large sink, 2:2 to 22.8 Pg C; Fig. 6). The lower error boundary, a terrestrial carbon source, is not supported from the d13CeCO2 evidence in ice cores and marine records (see section 7.1). Only those figures in the upper part of the range of estimates close the global carbon budget in the 1500s (requirement is a lower bound of 14.8 Pg C less atmospheric carbon plus the ocean source). This implies our estimate of the increase in carbon stocks following the Great Dying may be low.

Balancing the global carbon cycle in the 1500s implies substantial land sink, as even if the ocean was neutral and the atmospheric decline was at the lower boundary of 14.8 Pg C, the 9.8 Pg C land sink we calculate is smaller than the required land sink to close the carbon budget. At the other extreme, if the oceanic source was at the upper end of the range, 26.5 Pg C, and the atmospheric decline was also at the top of the range, 22.1 Pg C, then the land absorbed 48.6 Pg C in the 1500s. Of course, every term in the reconstructed 1500s carbon budget has a large uncertainty, as the land sink may be up to 22.8 Pg C, so it is possible to close the budget by combining the terms within these ranges. However, the likelihood, given an oceanic carbon source at this time, is that our central value for the land sink is underestimated. This may plausibly be because (i) the land temperature feedback is too weak, (ii) the CO2 feedback is too strong, or (iii) the per capita land use used in the LUC in the Americas calculation is too small. The latter may be particularly important, as the much larger per capita land use employed in the KK10 dataset may be closer to reality, which would substantially increase carbon uptake on land in the Americas (Kaplan et al., 2011). Indeed, if per capita land use was double our median estimate, at 2.08 ha per person, still a modest value, total uptake would be 18.8 Pg C. By contrast, the temperature and CO2 feedback parameters are unlikely to radically change with updated work. Overall, while the carbon budget in the 1500s is uncertain, it appears impossible to balance the global carbon cycle at this time without large terrestrial uptake including substantial land use change in the Americas as a result of the Great Dying.

8. Conclusion

We estimate that 55 million indigenous people died following the European conquest of the Americas beginning in 1492. This led to the abandonment and secondary succession of 56 million hectares of land. We calculate that this led to an additional 7.4 Pg C being removed from the atmosphere and stored on the land surface in the 1500s. This was a change from the 1400s of 9.9 Pg C (5 ppm CO2). Including feedback processes this contributed between 47% and 67% of the 15e22 Pg C (7e10 ppm CO2) decline in atmospheric CO2 between 1520 CE and 1610 CE seen in Antarctic ice core records. These changes show that the Great Dying of the Indigenous Peoples of the Americas is necessary for a parsimonious explanation of the anomalous decrease in atmospheric CO2 at that time and the resulting decline in global surface air temperatures. These changes show that human actions had global impacts on the Earth system in the centuries prior to the Industrial Revolution. Our results also show that this aspect of the Columbian Exchange e the globalisation of diseases e had global impacts on the Earth system, key evidence in the calls for the drop in atmospheric CO2 at 1610 CE to mark the onset of the Anthropocene epoch (Lewis and Maslin, 2015, 2018). We conclude that the Great Dying of the Indigenous Peoples of the Americas led to the abandonment of enough cleared land in the Americas that the resulting terrestrial carbon uptake had a detectable impact on both atmospheric CO2 and global surface air temperatures in the two centuries prior to the Industrial Revolution.


We thank Andrew Sluyter for discussions on the extent of preColumbian land use in Northern Mexico and its representation in LUC datasets, Jed Kaplan for providing the KK10 dataset, William A. Huber for statistical advice, and Joyce Chaplin and Matt Liebmann for discussions on an adequate name for the depopulation event. This work has been supported by the UK Natural Environmental Research Council Doctoral Training Partnership initiative (NE/ L002485/1).

Appendix A.

Supplementary data Supplementary data to this article can be found online at


Abrams, M.D., Nowacki, G.J., 2008. Native americans as active and passive promoters of mast and fruit trees in the eastern USA. Holocene 18, 1123e1137.

Ahlstrom, A., Xia, J., Arneth, A., Luo, Y., Smith, B., 2015. Importance of vegetation €
dynamics for future terrestrial carbon cycling. Environ. Res. Lett. 10, 054019.

Ahn, J., Brook, E.J., Mitchell, L., Rosen, J., McConnell, J.R., Taylor, K., Etheridge, D.,
Rubino, M., 2012. Atmospheric co2 over the last 1000 years: a high-resolution
record from the west antarctic ice sheet (wais) divide ice core. Glob. Biogeochem.
Cycles 26, 1e11.

Amiama, M., 1959. La poblacion de santo domingo. Clio 115, 118e133.

Amick, D.S., 2017. Evolving views on the pleistocene colonization of north America.
Quat. Int. 431, 125e151.
Anderson, R.C., 2006. Evolution and origin of the central grassland of north America
: climate, fire, and mammalian grazers. J. Torrey Bot. Soc. 133, 626e647.

Anderson-Teixeira, K.J., Snyder, P.K., Twine, T.E., Cuadra, S.V., Costa, M.H.,
DeLucia, E.H., 2012. Climate-regulation services of natural and agricultural
ecoregions of the americas. Nat. Clim. Change 2, 177e181.

Arora, V.K., Montenegro, A., 2011. Small temperature benefits provided by realistic
afforestation efforts. Nat. Geosci. 4, 514e518.

Barlow, J., Gardner, T.A., Lees, A.C., Parry, L., Peres, C.A., 2012. How pristine are
tropical forests? an ecological perspective on the pre-columbian human footprint
in amazonia and implications for contemporary conservation. Biol. Conserv.
151, 45e49.

Bauska, T.K., Joos, F., Mix, A.C., Roth, R., Ahn, J., Brook, E.J., 2015. Links between
atmospheric carbon dioxide, the land carbon reservoir and climate over the
past millennium. Nat. Geosci. 8, 383e387.

Beach, T., Luzzadder-Beach, S., Cook, D., Dunning, N., Kennett, D.J., Krause, S.,
Terry, R., Trein, D., Valdez, F., 2015. Ancient maya impacts on the earth’s surface:
an early anthropocene analog? Quat. Sci. Rev. 124, 1e30.

Beckermann, S., 1987. In: Turner II, B.L., Brush, S. (Eds.), Comparative Farming
Systems. Guildford Press, New York, pp. 55e94.

Benedictow, O.J., 2004. The Black Death 1346e1353: the Complete History. Boydell
Press, Rochester.

Berlinguer, G., 1993. The interchange of disease and health between the old and
new worlds. Int. J. Health Serv. 23, 703e715.

Berrío, J.C., Hooghiemstra, H., Marchant, R., Rangel, O., 2002. Late-glacial and holocene
history of the dry forest area in the south. J. Quat. Sci. 17, 667e682.

Biraben, J., 1979. Essai sur l’evolution du nombre des hommes. Popul. (French Ed. 
34, 13e25.

Black, F.L., 1992. Why did they die? Science 258, 1739e1741.

Black, F.L., 1994. An explanation of high death rates among new world peoples
when in contact with old world diseases. Perspect. Biol. Med. 37, 292e307.

Blarquez, O., Vanniere, B., Marlon, J.R., Daniau, A.L., Power, M.J., Brewer, S., Bartlein, P.J., 2014. Paleofire: an r package to analyse sedimentary charcoal records
from the global charcoal database to reconstruct past biomass burning.
Comput. Geosci. 72, 255e261.

Blatrix, R., Roux, B., Bearez, P., Prestes-Carneiro, G., Amaya, M., Aramayo, J.L., 
Rodrigues, L., Lombardo, U., Iriarte, J., De Souza, J.G., Robinson, M., Bernard, C.,
Pouilly, M., Durecu, M., Huchzermeyer, C.F., Kalebe, M., Ovando, A., McKey, D., 
2018. The unique functioning of a pre-columbian amazonian floodplain fishery.
Sci. Rep. 8, 1e16.

Bohm, F., Joachimski, M.M., Lehnert, H., Reitner, J., 2002. Evidence for preindustrial €
variations in the marine surface water carbonate system from coralline
sponges. Geochem. Geophys. Geosyst. 3

Bolt, J., van Zanden, J.L., 2014. The maddison project: collaborative research on
historical national accounts. Econ. Hist. Rev. 67, 627e651.

Borah, W., Cook, S.F., 1960. The Population of Central Mexico in 1548: An Analysis of
the Suma de Visitas de Pueblos, 1 ed. University of California Press, Berkeley.

Borah, W., Cook, S.F., 1963. The Aboriginal Population of Central Mexico on the Eve
of the Spanish Conquest. Nusmber no. 45 in. University of California publications
Ibero-Americana, University of Califronia Press, Berkeley.

Borah, W., Cook, S.F., 1969. Conquest and population: a demographic approach to
mexican history. Am. Philos. Soc. 113, 177e183.

Boserup, E., 1965. The conditions of agricultural growth. Popul. Stud. (NY) 20,

Boyd, R., 1992. Population decline from two epidemics on the northwest coast. In:
Verano, J.W., Ubelaker, D.H. (Eds.), Dis. Demogr. Am. Smithsonian Institution
Press, Washington, D.C., pp. 249e255

Brown, S., Lugo, A., 1992. Aboveground biomass estimates for tropical moist forests
of the brazilian amazon. Interciencia 17, 8e18.

Brugger, S.O., Gobet, E., van Leeuwen, J.F., Ledru, M.P., Colombaroli, D., van der
Knaap, W., Lombardo, U., Escobar-Torrez, K., Finsinger, W., Rodrigues, L.,
Giesche, A., Zarate, M., Veit, H., Tinner, W., 2016. Long-term man-environment
interactions in the bolivian amazon: 8000 years of vegetation dynamics. Quat.
Sci. Rev. 132, 114e128.

Bullock, M., Marquez, L., Hernandez, P., Ruíz, F., 2013. Paleodemographic age-atdeath
distributions of two mexican skeletal collections: a comparison of transition
analysis and traditional aging methods. Am. J. Phys. Anthropol. 152,

Bush, M.B., Silman, M.R., Listopad, C.M.C.S., 2007. A regional study of holocene
climate change and human occupation in peruvian amazonia. J. Biogeogr. 34,

Bush, M., Silman, M., McMichael, C., Saatchi, S., 2008. Fire, climate change and
biodiversity in amazonia: a late-holocene perspective. Philos. Trans. R. Soc. B
Biol. Sci. 363, 1795e1802.

Bush, M.B., Mosblech, N.a.S., Church, W., 2015. Climate change and the agricultural
history of a mid-elevation andean montane forest. Holocene 1.

Canuto, M.A., Estrada-Belli, F., Garrison, T.G., Houston, S.D., Acuna, M.J., Kov ~ 
ac, M., Marken, D., Nonded eo, P., Auld-Thomas, L., Castanet, C., Chatelain, D., 
Chiriboga, C.R., Drapela, T., Lieskovský, T., Tokovinine, A., Velasquez, A., 
Fernandez-Díaz, J.C., Shrestha, R., 2018. Ancient lowland maya complexity as
revealed by airborne laser scanning of northern Guatemala. Science 80, 361. eaau0137.

Carcaillet, C., Richard, P.J., Asnong, H., Capece, L., Bergeron, Y., 2006. Fire and soil
erosion history in east canadian boreal and temperate forests. Quat. Sci. Rev. 25,

Carson, J.F., Whitney, B.S., Mayle, F.E., Iriarte, J., Prumers, H., Soto, J.D., Watling, J.,
2014. Environmental impact of geometric earthwork construction in precolumbian
amazonia. Proc. Natl. Acad. Sci. 111, 10497e10502.

Carson, J.F., Watling, J., Mayle, F.E., Whitney, B.S., Iriarte, J., Prumers, H., Soto, J.D.,
2015. Pre-columbian land use in the ring-ditch region of the bolivian amazon.

Castilla-Beltran, A., Hooghiemstra, H., Hoogland, M.L., Pagan-Jimenez, J., van Geel, B., Field, M.H., Prins, M., Donders, T., Herrera Malatesta, E., Ulloa Hung, J.,
McMichael, C.H., Gosling, W.D., Hofman, C.L., 2018. Columbus’ footprint in
hispaniola: a paleoenvironmental record of indigenous and colonial impacts on
the landscape of the central cibao valley, northern Dominican Republic.
Anthropocene 22, 66e80.

Chazdon, R.L., Broadbent, E.N., Rozendaal, D.M.A., Bongers, F., Zambrano, A.M.A.,
Aide, T.M., Balvanera, P., Becknell, J.M., Boukili, V., Brancalion, P.H.S., Craven, D.,
Almeida-Cortez, J.S., Cabral, G.A.L., de Jong, B., Denslow, J.S., Dent, D.H.,
DeWalt, S.J., Dupuy, J.M., Duran, S.M., Espírito-Santo, M.M., Fandino, M.C., 
Cesar, R.G., Hall, J.S., Hern  andez-Stefanoni, J.L., Jakovac, C.C., Junqueira, A.B., 
Kennard, D., Letcher, S.G., Lohbeck, M., Martínez-Ramos, M., Massoca, P.,
Meave, J.A., Mesquita, R., Mora, F., Munoz, R., Muscarella, R., Nunes, Y.R.F., ~
Ochoa-Gaona, S., Orihuela-Belmonte, E., Pena-Claros, M., P ~ erez-García, E.A., 
Piotto, D., Powers, J.S., Rodríguez-Velazquez, J., Romero-Perez, I.E., Ruíz, J., 
Saldarriaga, J.G., Sanchez-Azofeifa, A., Schwartz, N.B., Steininger, M.K.,
Swenson, N.G., Uriarte, M., van Breugel, M., van der Wal, H., Veloso, M.D.M.,
Vester, H., Vieira, I.C.G., Bentos, T.V., Williamson, G.B., Poorter, L., 2016. Carbon
sequestration potential of second-growth forest regeneration in the Latin
american tropics. Sci. Adv 2.

Chepstow-Lusty, A., Jonsson, P., 2000. Inca agroforestry: lessons from the past.
AMBIO A J. Hum. Environ. 29, 322e328.

Chepstow-Lusty, A.J., Bennett, K.D., Switsur, V.R., Kendall, A., 1996. 4000 years of
human impact and vegetation change in the central peruvian andes – with
events parallelling the maya record? Antiquity 70, 824e833.

Chepstow-Lusty, A.J., Bennett, K.D., Fjeldså, J., Kendall, A., Galiano, W., Herrera, A.T.,
1998. Tracing 4,000 years of environmental history in the cuzco area, Peru, from
the pollen record. Mt. Res. Dev. 18, 159e172.

Chepstow-Lusty, A., Frogley, M.R., Bauer, B.S., Bush, M.B., Herrera, A.T., 2003. A late
holocene record of arid events from the cuzco region, Peru. J. Quat. Sci. 18,

Chepstow-Lusty, A.J., Frogley, M.R., Bauer, B.S., Leng, M.J., Boessenkool, K.P.,
Carcaillet, C., Ali, A.A., Gioda, A., 2009. Putting the rise of the inca empire within
a climatic and land management context. Clim. Past 5, 375e388.

Clark, C., 1967. Population Growth and Land Use. London.

Clark, J.S., 1988. Particle motion and the theory of charcoal analysis: source area,
transport, deposition, and sampling. Quat. Res. 30, 67e80.

Clark, K., Uhl, C., 1987. Farming, fishing, and fire in the history of the upper rio negro
region of Venezuela. Hum. Ecol. 15, 1e26.

Clement, R.M., Horn, S.P., 2001. Pre-columbian land-use history in Costa Rica: a
3000-year record of forest clearance, agriculture and fires from laguna zoncho.
Holocene 11, 419e426.

Clement, C.R., Denevan, W.M., Heckenberger, M.J., Junqueira, A.B., Neves, E.G.,
Teixeira, W.G., Woods, W.I., 2015. The domestication of amazonia before european
conquest. Proc. Biol. Sci. 282, 20150813.

Collins, J.L., 1983. Fertility determinants in a high andes community. Popul. Dev. Rev.
9, 61e75.

Colmenares, G., 2007. In: La formacion de la economia colonial (1500-1740),
pp. 21e60. Bogota: Planeta.

Conant, R.T., Cerri, C.E.P., Osborne, B.B., Paustian, K., 2017. Grassland management
impacts on soil carbon stocks: a new synthesis. Ecol. Appl. 27, 662e668. https://

Conserva, M.E., Byrne, R., 2002. Late holocene vegetation change in the sierra madre
oriental of central Mexico. Quat. Res. 58, 122e129.

Cook, D.N., 1981. Demographic Collapse: Indian Peru, 1 ed. Cambridge University
Press, Cambridge, pp. 1520e1620.

Cook, N.D., 1998. Born to Die. Disease and New World Conquest, 1 ed. Cambridge
University Press, Cambridge, pp. 1492e1650.

Cook, S.F., Borah, W., 1960. The Indian Population of Central Mexico, 1 ed. University
of California Press, Berkeley, pp. 1531e1610.

Cook, S.F., Borah, W., 1971. Essays in Population History, Vol. III: Mexico and California,
1 ed. University of California Press, Berkeley.

Cook, S.F., Simpson, L.B., 1948. The Populanon of Central mexico in the Sixteenth
Century. Ibero-Americana Ser. Los Angeles Berkeley (Univ. Cal. Press, p. 31.

Cordova, E., 1968. La encomienda y la desaparici  on de los indios en las antillas 
mayores. Caribb. Stud. 8, 23e49.

Cox, P.M., Pearson, D., Booth, B.B., Friedlingstein, P., Huntingford, C., Jones, C.D.,
Luke, C.M., 2013. Sensitivity of tropical carbon to climate change constrained by
carbon dioxide variability. Nature 494, 341e344.

Crosby, A.W., 1972. The Columbian Exchange: Biological and Cultural Consequences
of 1492, 1 ed. Greenwood Publishing Group, Santa Barbara.

Crosby, A.W., 1976. Virgin soil epidemics as a factor in the aboriginal depopulation
in America. William Mary Q. 33, 289e299.
Davies II, E.J., 2012. The americas, 1450e2000. In: Bentley, J.H. (Ed.), Oxford Handb.
World Hist, vol. 1. Oxford University Press, pp. 1e25.

De Souza, J.G., Schaan, D.P., Robinson, M., Barbosa, A.D., Aragao, L.E., Marimon, B.H., ~
Marimon, B.S., Da Silva, I.B., Khan, S.S., Nakahara, F.R., Iriarte, J., 2018. Precolumbian
earth-builders settled along the entire southern rim of the amazon.
Nat. Commun. 9

Delcourt, P., Delcourt, H., 2004. Prehistoric Native Americans and Ecological Change
Human Ecosystems in Eastern North America since the Pleistocene, 1 ed.
Cambridge University Press, Cambridge.

Denevan, W.M., 1970. The aboriginal population of western amazonia in relation to
habitat and subsistence. Rev. Geogr. 72, 61e86.

Denevan, W.M., 1992a. The Native Population of the Americas in 1492, 2 ed. University
of Wisconsin Press, Madison.

Denevan, W.M., 1992b. The pristine myth: the landscape of the americas in 1492.
Ann. Assoc. Am. Geogr. 82, 369e385.

Denevan, W., 2001. Cultivated Landscapes of Native Amazonia and the Andes, 1 ed.
Oxford Universtiy Press, New York.

Denevan, W., 2003. The native population of amazonia in 1492 reconsidered. Rev.
Indias R. I LXIII, 175e188.

Denevan, W.M., 2012. Rewriting the late pre-european history of amazonia. J. Lat.
Am. Geogr. 11, 9e24.

Denevan, W.M., 2014. Estimating amazonian indian numbers in 1492. J. Lat. Am.
Geogr. 13, 207e221.

Diamond, J., 1997. Guns, Germs, and Steel: the Fates of Human Societies. W.W.
Norton & Company., London.

Dillehay, T.D., Rossen, J., Andres, T.C., Williams, D.E., 2007. Preceramic adoption of
peanut, squash, and cotton in northern Peru. Science 80 (316), 1890e1893.

Dobyns, H.F., 1966. An Appraisal of Techniques with a New Hemispheric Estimate.

Dobyns, H.F., 1983. Their Number Become Thinned: Native American Population
Dynamics in Eastern North America. University of Tennessee Press, Tennessee.

Dobyns, H.F., 1993. Disease transfer at contact. Annu. Rev. Anthropol. 22, 273e291.

Don, A., Schumacher, J., Freibauer, A., 2011. Impact of tropical land-use change on
soil organic carbon stocks – a meta-analysis. Glob. Chang. Biol. 17, 1658e1670.

Doolittle, W.E., 1992. Agriculture in north America on the eve of contact: a reassessment.
Ann. Assoc. Am. Geogr. 82, 386e401.

Drucker, P., Heizer, R.F., 1960. A study of the milpa system of la venta island and its
archaeological implications. SW. J. Anthropol. 16, 36e45.

Dull, R.A., Nevle, R.J., Woods, W.I., Bird, D.K., Avnery, S., Denevan, W.M., 2010. The
columbian encounter and the little ice age: abrupt land use change, fire, and
greenhouse forcing. Ann. Assoc. Am. Geogr. 100, 755e771.

Durand, J.D., 1977. Historical estimates of world population : an evaluation. Popul.
Dev. Rev. 3, 253e296.

Erickson, C.L., 2000. An artificial landscape-scale fishery in the bolivian amazon.
Nature 408, 190e193.

Erickson, C.L., 2010. The transformation of environment into landscape: the historical
ecology of monumental earthwork construction in the bolivian amazon.
Diversity 2, 618e652.

Etter, A., 2000. Patterns of landscape transformation in Colombia , with emphasis in
the andean region patterns of landscape transformation in Colombia , with
emphasis in the andean region. Ambio 29, 432e441.

Faust, F.X., Gnecco, C., Mannstein, H., Stamm, J., 2006. Evidence for the postconquest
demographic collapse of the americas in historical co2 levels. Earth Interact. 10

Fehse, J., Hofstede, R., Aguirre, N., Paladines, C., Kooijman, A., Sevink, J., 2002. High
altitude tropical secondary forests: a competitive carbon sink? For. Ecol.
Manage. 163, 9e25.

Flantua, S.G.A., Hooghiemstra, H., Vuille, M., Behling, H., Carson, J.F., Gosling, W.D., Hoyos, I., Ledru, M.P., Montoya, E., Mayle, F., Maldonado, A., Rull, V., Tonello, M.S., Whitney, B.S., Gonz
alez-Arango, C., 2016. Climate variability and human impact in south America during the last 2000 years: synthesis and perspectives from pollen records. Clim. Past 12, 483e523.

Follett, R.F., Kimble, J., Lal, R., 2001. The Potential of US Grazing Lands to Sequester
Carbon and Mitigate the Greenhouse Effect. CRC Press, Chelsea, Michigan.

Fowler, W.R., 1988. La poblacion nativa de El Salvador al momento de la conquista 
espanola. Mesoamerica 15, 79 ~ e116.

Francey, R.J., Allison, C.E., Etheridge, D.M., Trudinger, C.M., Enting, I.G.,
Leuenberger, M., Langenfelds, R.L., Michel, E., Steele, L.P., 1999. A 1000-year high
precision record of d13c in atmospheric co2. Tellus Ser. B Chem. Phys. Meteorol.
51, 170e193.

Friend, A.D., Lucht, W., Rademacher, T.T., Keribin, R., Betts, R., Cadule, P., Ciais, P.,
Clark, D.B., Dankers, R., Falloon, P.D., Ito, A., Kahana, R., Kleidon, A., Lomas, M.R.,
Nishina, K., Ostberg, S., Pavlick, R., Peylin, P., Schaphoff, S., Vuichard, N.,
Warszawski, L., Wiltshire, A., Woodward, F.I., 2014. Carbon residence time
dominates uncertainty in terrestrial vegetation responses to future climate and
atmospheric co 2. Proc. Natl. Acad. Sci. 111, 3280e3285.

Frolicher, T.L., Joos, F., Raible, C.C., 2011. Sensitivity of atmospheric co2 and climate €
to explosive volcanic eruptions. Biogeosciences 8, 2317e2339.

Galbraith, D., Malhi, Y., Affum-Baffoe, K., Castanho, A.D., Doughty, C.E., Fisher, R.A.,
Lewis, S.L., Peh, K.S., Phillips, O.L., Quesada, C.A., Sonke, B., Lloyd, J., 2013. 
Residence times of woody biomass in tropical forests. Plant Ecol. Divers. 6,

Gilroy, J.J., Woodcock, P., Edwards, F.A., Wheeler, C., Baptiste, B.L.G., Medina
Uribe, C.A., Haugaasen, T., Edwards, D.P., 2014. Cheap carbon and biodiversity
co-benefits from forest regeneration in a hotspot of endemism. Nat. Clim.
Change 4, 503e507.

Goebel, T., Waters, M.R., O’Rourke, D.H., 2008. The late pleistocene dispersal of
modern humans in the americas. Science 319, 1497e1502.

Goldberg, A., Mychajliw, A.M., Hadly, E.A., 2016. Post-invasion demography of
prehistoric humans in south America. Nature 532, 232e235.

Goman, M., Byrne, R., 1998. A 5000-year record of agriculture and tropical forest
clearance in the tuxtlas, veracruz, Mexico. Holocene 8, 83e89.

Griscom, B.W., Adams, J., Ellis, P.W., Houghton, R.A., Lomax, G., Miteva, D.A.,
Schlesinger, W.H., Shoch, D., Siikamaki, J.V., Smith, P., Woodbury, P., Zganjar, C., €
Blackman, A., Campari, J., Conant, R.T., Delgado, C., Elias, P., Gopalakrishna, T.,
Hamsik, M.R., Herrero, M., Kiesecker, J., Landis, E., Laestadius, L., Leavitt, S.M.,
Minnemeyer, S., Polasky, S., Potapov, P., Putz, F.E., Sandermanc, J., Silvius, M.,
Wollenberg, E., Fargione, J., 2017. Natural climate solutions symposium. Proc.
Natl. Acad. Sci. 114, 1e6.

Guariguata, M.R., Ostertag, R., 2001. Neotropical secondary forest succession:
changes in structural and functional characteristics. For. Ecol. Manage. 148,

Guderjan, T., Beach, T., Luzzadder-Beach, S., Bozarth, S., 2009. Understanding the
causes of abandonment in the maya lowlands. Archaeol. Rev. Camb. 24, 99e121.
Hays, J., 2005. Epidemics and Pandemics: Their Impacts on Human History. ABCCLIO,
Santa Barbara.

Heckenberger, M., Neves, E.G., 2009. Amazonian archaeology. Annu. Rev. Anthropol.
38, 251e266.

Heckenberger, M.J., Petersen, J.B., Neves, E.G., 1999. Village size and permanence in
amazonia: two archaeological examples from Brazil. Lat. Am. Antiq. 10,

Heckenberger, M.J., Kuikuro, A., Kuikuro, U.T., Russell, J.C., Schmidt, M., Fausto, C.,
Franchetto, B., 2003. Amazonia 1492: pristine forest or cultural parkland? Science
80 (301), 1710e1714.

Hemming, J., 1995. Red Gold: Conquest of the Brazilian Indians, 2 ed. Papermac,

Henige, D., 1978. On the contact population of hispaniola: history as higher
mathematics. Hisp. Am. Hist. Rev. 58, 217e237.

Henige, D., 1998. Numbers from Nowhere: the American Indian Contact Population
Debate, 1 ed. University of Oklahoma Press.

Higuera-Gundy, A., Brenner, M., Hodell, D.A., Curtis, J.H., Leyden, B.W.,
Binford, M.W., 1999. A 10,300 14c yr record of climate and vegetation change
from Haiti. Quat. Res. 52, 159e170.

Hilbert, L., Neves, E.G., Pugliese, F., Whitney, B.S., Shock, M., Veasey, E., Zimpel, C.A.,
Iriarte, J., 2017. Evidence for mid-holocene rice domestication in the americas.
Nat. Ecol. Evol. 1, 1693e1698.

Hooghiemstra, H., Olijhoek, T., Hoogland, M., Prins, M., van Geel, B., Donders, T., Gosling, W., Hofman, C., 2018. Columbus’ environmental impact in the new
world: land use change in the yaque river valley, Dominican Republic. Holocene.

Hughes, R.F., Kauffman, J.B., Jaramillo, V.J., 1999. Biomass, carbon, and nutrient
dynamics of secondary forests in a humid tropical region of Mexico. Ecology 80,

Hunter, R., Sluyter, A., 2015. Sixteenth-century soil carbon sequestration rates based
on mexican land-grant documents. Holocene.

Iriarte, J., Power, M.J., Rostain, S., Mayle, F.E., Jones, H., Watling, J., Whitney, B.S., McKey, D.B., 2012. Fire-free land use in pre-1492 amazonian savannas. Proc.
Natl. Acad. Sci. 109, 6473e6478.

Jones, D.S., 2003. Virgin soils revisited. William Mary Q. 60, 703e742.

Joos, F., Meyer, R., Bruno, M., Leuenberger, M., 1999. The variability in the carbon
sinks as reconstructed for the last 1000 years. Geophys. Res. Lett. 26, 1437.

Joralemon, D., 1982. New world depopulation and the case of disease. J. Anthropol.
Res. 38, 108e127.

Kaplan, J.O., Krumhardt, K.M., Ellis, E.C., Ruddiman, W.F., Lemmen, C.,
Goldewijk, K.K., 2011. Holocene carbon emissions as a result of anthropogenic
land cover change. Holocene 21, 775e791.

Kendall, A., 2013. Applied archaeology in the andes: the contribution of prehispanic
agricultural terracing to environmental and rural development strategies.
In: IJM, D., FN, M. (Eds.), Humans Environ. New Archaeol. Perspect.
Twenty-First Century. Oxford University Press, Oxford, pp. 153e169.

Kennedy, L.M., Horn, S.P., Orvis, K.H., 2006. A 4000-year record of fire and forest
history from valle de bao, cordillera central, Dominican Republic. Palaeogeogr.
Palaeoclimatol. Palaeoecol. 231, 279e290.

Kimble, J., Birdsey, R., Heath, L., Lal, R., 2002. The Potential of US Forest Soils to
Sequester Carbon and Mitigate the Greenhouse Effect. CRC Press, Boca-Raton.

Klein Goldewijk, K., Beusen, A., Janssen, P., 2010. Long-term dynamic modeling of
global population and built-up area in a spatially explicit way: hyde 3.1. Holocene
20, 565e573.

Klein Goldewijk, K., Beusen, A., Van Drecht, G., De Vos, M., 2011. The hyde 3.1
spatially explicit database of human-induced global land-use change over the
past 12,000 years. Glob. Ecol. Biogeogr. 20, 73e86.

Klein Goldewijk, K., Dekker, S.C., van Zanden, J.L., 2017. Per-capita estimations of
long-term historical land use and the consequences for global change research.
J. Land Use Sci. 12, 313e337.

Kroeber, A., 1939. Cultural and natural areas of native north America. Univ. Calif.
Publ. Am. Archaeol. Ethnol. 38, 242.

Krumhardt, K.M., 2010. Methodology for Worldwide Population Estimates: 1000 BC
to 1850. Technical Report.

Lange, F., 1971. Una reevaluacion de la poblaci  on del norte de yucat  an en el tiempo 
del contacto espanol: 1528. Am. Indig. 31, 117 ~ e139.

Levis, C., Costa, F.R.C., Bongers, F., Pena-Claros, M., Clement, C.R., Junqueira, A.B., ~
Neves, E.G., Tamanaha, E.K., Figueiredo, F.O.G., Salomao, R.P., Castilho, C.V., Magnusson, W.E., Phillips, O.L., Guevara, J.E., Sabatier, D., Molino, J.F., Lopez, D.C., Mendoza, A.M., Pitman, N.C.A., Duque, A., Vargas, P.N., Zartman, C.E.,  Vasquez, R., Andrade, A., Camargo, J.L., Feldpausch, T.R., Laurance, S.G.W., Laurance, W.F., Killeen, T.J., Nascimento, H.E.M., Montero, J.C., Mostacedo, B., Amaral, I.L., Guimaraes Vieira, I.C., Brienen, R., Castellanos, H., Terborgh, J., Carim, M.d.J.V., Guimaraes, J.R.d.S., Coelho, L.D.S., Matos, F.D.D.A., Wittmann, F., Mogollon, H.F., Damasco, G., D avila, N., García-Villacorta, R., Coronado, E.N.H.,
Emilio, T., Filho, D.d.A.L., Schietti, J., Souza, P., Targhetta, N., Comiskey, J.A.,
Marimon, B.S., Marimon, B.H., Neill, D., Alonso, A., Arroyo, L., Carvalho, F.A., de
Souza, F.C., Dallmeier, F., Pansonato, M.P., Duivenvoorden, J.F., Fine, P.V.A.,
Stevenson, P.R., Araujo-Murakami, A., Aymard, C.,G.A., Baraloto, C., do
Amaral, D.D., Engel, J., Henkel, T.W., Maas, P., Petronelli, P., Revilla, J.D.C.,
Stropp, J., Daly, D., Gribel, R., Paredes, M.R., Silveira, M., Thomas-Caesar, R.,
Baker, T.R., da Silva, N.F., Ferreira, L.V., Peres, C.A., Silman, M.R., Ceron, C., 
Valverde, F.C., Di Fiore, A., Jimenez, E.M., Mora, M.C.P., Toledo, M., Barbosa, E.M.,
Bonates, L.C.d.M., Arboleda, N.C., Farias, E.d.S., Fuentes, A., Guillaumet, J.L.,
Jørgensen, P.M., Malhi, Y., de Andrade Miranda, I.P., Phillips, J.F., Prieto, A.,
Rudas, A., Ruschel, A.R., Silva, N., von Hildebrand, P., Vos, V.A., Zent, E.L., Zent, S.,
Cintra, B.B.L., Nascimento, M.T., Oliveira, A.A., Ramirez-Angulo, H., Ramos, J.F.,
Rivas, G., Schongart, J., Sierra, R., Tirado, M., van der Heijden, G., Torre, E.V., €
Wang, O., Young, K.R., Baider, C., Cano, A., Farfan-Rios, W., Ferreira, C.,
Hoffman, B., Mendoza, C., Mesones, I., Torres-Lezama, A., Medina, M.N.U., van
Andel, T.R., Villarroel, D., Zagt, R., Alexiades, M.N., Balslev, H., Garcia-Cabrera, K.,
Gonzales, T., Hernandez, L., Huamantupa-Chuquimaco, I., Manzatto, A.G.,
Milliken, W., Cuenca, W.P., Pansini, S., Pauletto, D., Arevalo, F.R., Reis, N.F.C.,
Sampaio, A.F., Giraldo, L.E.U., Sandoval, E.H.V., Gamarra, L.V., Vela, C.I.A., ter
Steege, H., 2017. Persistent effects of pre-columbian plant domestication on
amazonian forest composition. Science 80 (355), 925e931.

Levis, C., Flores, B.M., Moreira, P.A., Luize, B.G., Alves, R.P., Franco-Moraes, J., Lins, J.,
Konings, E., Pena-Claros, M., Bongers, F., Costa, F.R.C., Clement, C.R., 2018. How ~
people domesticated amazonian forests. Front. Ecol. Evol. 5

Lewis, C.M., 2010. Hierarchical modeling of genome-wide short tandem repeat (str)
markers infers native american prehistory. Am. J. Phys. Anthropol. 141,

Lewis, S.L., Maslin, M.A., 2015. Defining the anthropocene. Nature 519, 171e180.

Lewis, S.L., Maslin, M.A., 2018. The Human Planet: How We Created the Anthropocene,
1 ed. Pengiun, London.

Liebmann, M.J., Farella, J., Roos, C.I., Stack, A., Martini, S., Swetnam, T.W., 2016.
Native american depopulation, reforestation, and fire regimes in the southwest
United States, 1492e1900 ce. Proc. Natl. Acad. Sci. 113, E696eE704. https://

Lindenau, J.D., Salzano, F.M., Guimaraes, L.S.P., Callegari-Jacques, S.M., Hurtado, A.M., Hill, K.R., Petzl-Erler, M.L., Tsuneto, L.T., Hutz, M.H., 2013. Distribution
patterns of variability for 18 immune system genes in amerindians –
relationship with history and epidemiology. Tissue Antigens 82, 177e185.

Littman, R.J., 2009. The plague of athens: epidemiology and paleopathology. Mt.
Sinai J. Med. A J. Transl. Pers. Med. 76, 456e467.

Livi-Bacci, M., 2008. Conquest: the Destruction of the American Indios. Polity Press,

Loughlin, N.J., Gosling, W.D., Montoya, E., 2018. Identifying environmental drivers of
fungal non-pollen palynomorphs in the montane forest of the eastern andean
flank, Ecuador. Quat. Res. (United States) 89, 119e133.

Lovell, W.G., Lutz, C.H., 1995. Demography and Empire: A Guide to the Population
History of Spanish Central America, 1 ed. Westview Press, Boulder,
pp. 1500e1821.

Lovell, W.G., Lutz, C.H., 1996. A dark obverse: maya survival in Guatemala: 1520-
1994. Lat. Am. Geogr. 86, 398e407.

Lovell, W.G., Swezey, W., 1982. The population of southern Guatemala at Spanish
contact. Can. J. Anthropol. 3, 71e84.

MacFarling Meure, C., Etheridge, D., Trudinger, C., Steele, P., Langenfelds, R., Van
Ommen, T., Smith, A., Elkins, J., 2006. Law dome co2, ch4 and n2o ice core
records extended to 2000 years bp. Geophys. Res. Lett. 33, 2000e2003. https://

Maddison, A., 2001. The World Economy: A Millenial Perspective. OECD, Paris,
p. 384.

Maezumi, S.Y., Alves, D., Robinson, M., de Souza, J.G., Levis, C., Barnett, R.L., Almeida
de Oliveira, E., Urrego, D., Schaan, D., Iriarte, J., 2018a. The legacy of 4,500 years
of polyculture agroforestry in the eastern amazon. Native Plants 4, 540e547.

Maezumi, S.Y., Robinson, M., de Souza, J., Urrego, D.H., Schaan, D., Alves, D., Iriarte, J.,
2018b. New insights from pre-columbian land use and fire management in
amazonian dark earth forests. Front. Ecol. Evol. 6, 1e23.

Malhi, Y., Wood, D., Baker, T.R., Wright, J., Phillips, O.L., Cochrane, T., Meir, P.,
Chave, J., Almeida, S., Arroyo, L., Higuchi, N., Killeen, T.J., Laurance, S.G.,
Laurance, W.F., Lewis, S.L., Monteagudo, A., Neill, D.A., Vargas, P.N.,
Pitman, N.C.A., Quesada, C.A., Salomao, R., Silva, J.N.M., Lezama, A.T., Terborgh, J., ~
Martínez, R.V., Vinceti, B., 2006. The regional variation of aboveground live
biomass in old-growth amazonian forests. Glob. Chang. Biol. 12, 1107e1138.

Mann, C., 2005. New Revelations of the Americas before Columbus, 2 ed. Knopf,
New York, N.Y, p. 1491.

Marín-Spiotta, E., Sharma, S., 2013. Carbon storage in successional and plantation
forest soils: a tropical analysis. Glob. Ecol. Biogeogr. 22, 105e117.

Marlon, J.R., Bartlein, P.J., Carcaillet, C., Gavin, D.G., Harrison, S.P., Higuera, P.E.,
Joos, F., Power, M.J., Prentice, I.C., 2008. Climate and human influences on global
biomass burning over the past two millennia. Nat. Geosci. 1, 697e702. https://

Marlon, J.R., Bartlein, P.J., Daniau, A.l., Harrison, S.P., Maezumi, S.Y., Power, M.J.,
Tinner, W., Vanniere, B., 2013. Global biomass burning: a synthesis and review 
of holocene paleo fire records and their controls. Quat. Sci. Rev. 65, 5e25.

Marr, J.S., Kiracofe, J.B., 2000. Was the huey cocoliztli a haemorrhagic fever? Med.
Hist. 44, 341e362.

Martin, P.A., Newton, A.C., Bullock, J.M., 2013. Carbon pools recover more quickly
than plant biodiversity in tropical secondary forests. Proc. R. Soc. B Biol. Sci. 280, 20132236e20132236.

Mayle, F.E., Power, M.J., 2008. Impact of a drier early-mid-holocene climate upon
amazonian forests. Philos. Trans. R. Soc. Lond. B Biol. Sci. 363, 1829e1838.

McCaa, R., 1995. Spanish and Nahuatl views on smallpox and demographic catastrophe
in Mexico. J. Interdiscip. Hist. 25, 397e431.

McEvedy, C., Jones, R., 1977. Atlas of World Population History. Pengiun.

McKey, D., Rostain, S., Iriarte, J., Glaser, B., Birk, J.J., Holst, I., Renard, D., 2010. Precolumbian
agricultural landscapes, ecosystem engineers, and self-organized
patchiness in amazonia. Proc. Natl. Acad. Sci. U.S.A. 107, 7823e7828. https://

McMichael, C.H., Piperno, D.R., Bush, M.B., Silman, M.R., Zimmerman, A.R.,
Raczka, M.F., Lobato, L.C., 2012. Sparse pre-columbian human habitation in
western amazonia. Science 80 (336), 1429e1431.

McMichael, C.H., Palace, M.W., Bush, M.B., Braswell, B., Hagen, S., Neves, E.G.,
Silman, M.R., Tamanaha, E.K., Czarnecki, C., 2014. Predicting pre-columbian
anthropogenic soils in amazonia. Proc. R. Soc. B Biol. Sci. 281
10.1098/rspb.2013.2475, 20132475e20132475.

McMichael, C.N.H., Matthews-Bird, F., Farfan-Rios, W., Feeley, K.J., 2017. Ancient
human disturbances may be skewing our understanding of amazonian forests.
Proc. Natl. Acad. Sci. 114, 201614577.

McNeill, W.H., 1977. Plagues and People. Anchor, New York, N.Y.

Meggers, B., 1992. Disease and Demography in the Americas, 1 ed. Smithsonian
Institution Press, Washington, D.C.

Meggers, B., 2001. The continuing quest for el dorado: round two. Lat. Am. Antiq. 12,

Mickleburgh, H.L., Pagan-Jim  enez, J.R., 2012. New insights into the consumption of 
maize and other food plants in the pre-columbian caribbean from starch grains
trapped in human dental calculus. J. Archaeol. Sci. 39, 2468e2478. https://

Milner, G.R., Chaplin, G., 2010. Eastern north american population at ca. ad 1500.
Am. Antiq. 75, 707e726.

Montenegro, V., Mostajo de Muente, A., Mostajo de Muente, P., 1990. La Situacion
Poblacional Peruana: Balance Y Perspectivas. Technical Report. Instituto Andino
de Estudios en Poblacion y Desarrollo, Lima.

Moreno-Mayar, J.V., Potter, B.A., Vinner, L., Steinrücken, M., Rasmussen, S.,
Terhorst, J., Kamm, J.A., Albrechtsen, A., Malaspinas, A.S., Sikora, M.,
Reuther, J.D., Irish, J.D., Malhi, R.S., Orlando, L., Song, Y.S., Nielsen, R., Meltzer, D.J., Willerslev, E., 2018. Terminal pleistocene alaskan genome reveals first founding population of native americans. Nature 553, 203e207.

Morfín, L.M., Storey, R., 2016. Population history in precolumbian and colonial
times. In: Nichols, D.L., Rodríguez-Alegría, E. (Eds.), Oxford Handb. Aztecs, vol. 1.
Oxford University Press, Oxford, pp. 1e16.

Morison, S., 1948. Admiral of the Ocean Sea (Boston).

Morley, S., 1968. The Ancient Maya, 3 ed. Stanford University Press, Stanford.
Moya Pons, F., 1979. Datos para el estudio de la demografía aborigen en santo
domingo. Jahrb. Gesch. Lat. Am. 16, 1e12.

Moya Pons, F., 1992. The politics of forced indian labour in la espanola 1493 ~ AS ¸ 1520.
Antiquity 66, 130e139.

Muller, S.D., Richard, P.J., Guiot, J., De Beaulieu, J.L., Fortin, D., 2003. Postglacial
climate in the st. lawrence lowlands, southern quebec: pollen and lake-level 
evidence. Palaeogeogr. Palaeoclimatol. Palaeoecol. 193, 51e72.

Myhre, G., Highwood, E.J., Shine, K.P., Stordal, F., 1998. New estimates of radiative
forcing due to wel mixed greenhouse gases. Geophys. Res. Lett. 25, 2715e2718.

Myhre, G., Shindell, D., Breon, F.M., Collins, W., Fuglestvedt, J.S., Huang, J., Koch, D., 
Lamarque, J.F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G.,
Takemura, T., Zhang, H., 2013. Anthropogenic and natural radiative forcing. In:
Stocker, T.F., Qin, D., Plattner, G.K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A.,
Xia, Y., Bex, V., Midgley, P. (Eds.), Clim. Chang. 2013 – Phys. Sci. Basis, vol.
9781107057. Cambridge University Press, Cambridge, pp. 659e740. https:// (chapter 8).

Nepstad, D., Carvalho, G., Barros, A.C., Alencar, A., Capobianco, J.P., Bishop, J.,
Moutinho, P., Lefebvre, P., Silva, U.L., Prins, E., 2001. Road paving, fire regime
feedbacks, and the future of amazon forests. For. Ecol. Manage. 154, 395e407.

Neukom, R., Gergis, J., Karoly, D.J., Wanner, H., Curran, M., Elbert, J., Gonzalezrouco,
F., Linsley, B.K., Moy, A.D., Mundo, I., Raible, C.C., Steig, E.J., Ommen, T.V.,
Vance, T., Villalba, R., Zinke, J., Frank, D., 2014. Inter-hemispheric temperature
variability over the past millennium. Nat. Clim. Change 5, 362e367. https://

Nevle, R.J., Bird, D.K., 2008. Effects of syn-pandemic fire reduction and reforestation
in the tropical americas on atmospheric co2 during european conquest.
Palaeogeogr. Palaeoclimatol. Palaeoecol. 264, 25e38.

Nevle, R.J., Bird, D.K., Ruddiman, W.F., Dull, R.A., 2011. Neotropical human-landscape
interactions, fire, and atmospheric co2 during european conquest. Holocene 21,

Newson, L.A., 1996. The population of the amazon basin in 1492: a view from the
ecuadorian headwaters. Trans. Inst. Br. Geogr. 21, 5e26.

Niemann, H., Behling, H., 2010. Late holocene environmental change and human
impact inferred from three soil monoliths and the laguna zurita multi-proxi
record in the southeastern ecuadorian andes. Veg. Hist. Archaeobotany 19,

Niemann, H., Matthias, I., Michalzik, B., Behling, H., 2013. Late holocene human
impact and environmental change inferred from a multi-proxy lake sediment
record in the loja region, southeastern Ecuador. Quat. Int. 308e309, 253e264.

Noymer, A., 2011. Population decline in post-conquest America: the role of disease.
Popul. Dev. Rev. 37, 178e183.
Nunn, N., Qian, N., 2010. The columbian exchange: a history of disease, food, and
ideas. J. Econ. Perspect. 24, 163e188.

Ocampo, J.S., 1997. Historia economica de Colombia. Technical Report. Presidencia
de la Republica, Bogota.

Orihuela-Belmonte, D.E., De Jong, B.H.J., Mendoza-Vega, J., Van der Wal, J., PazPellat,
F., Soto-Pinto, L., Flamenco-Sandoval, A., 2013. Carbon stocks and accumulation
rates in tropical secondary forests at the scale of community, landscape
and forest type. Agric. Ecosyst. Environ. 171, 72e84.

Ortiz Aguilu, J.J., Rivera Melendez, J., Principe Jacome, A., Melendez Maiz, M., Lavergne
Colberg, M., 1991. Intensive Agriculture in Pre-columbian West Indies:
the Case for Terraces.

Owens, M.J., Lockwood, M., Hawkins, E., Usoskin, I., Jones, G.S., Barnard, L.,
Schurer, A., Fasullo, J., 2017. The maunder minimum and the little ice age: an
update from recent reconstructions and climate simulations. J. Sp. Weather Sp.
Clim. 7, A33.

Palace, M.W., McMichael, C.N., Braswell, B.H., Hagen, S.C., Bush, M.B., Neves, E., Tamanaha, E., Herrick, C., Frolking, S., 2017. Ancient amazonian populations left
lasting impacts on forest structure. Ecosphere 8.

Park, J., Byrne, R., Bohnel, H., Garza, R.M., Conserva, M., 2010. Holocene climate €
change and human impact, central Mexico: a record based on maar lake pollen
and sediment chemistry. Quat. Sci. Rev. 29, 618e632.

Pedersen, M.W., Ruter, A., Schweger, C., Friebe, H., Staff, R.A., Kjeldsen, K.K.,
Mendoza, M.L.Z., Beaudoin, A.B., Zutter, C., Larsen, N.K., Potter, B.A., Nielsen, R.,
Rainville, R.A., Orlando, L., Meltzer, D.J., Kjær, K.H., Willerslev, E., 2016. Postglacial
viability and colonization in north America’s ice-free corridor. Nature
537, 1e15.

Piperno, D.R., Dillehay, T.D., 2008. Starch grains on human teeth reveal early broad
crop diet in northern Peru. Proc. Natl. Acad. Sci. 105.

Piperno, D.R., Ranere, A.J., Holst, I., Iriarte, J., Dickau, R., 2009. Starch grain and
phytolith evidence for early ninth millennium b.p. maize from the central
balsas river valley, Mexico. Proc. Natl. Acad. Sci. 106, 5019e5024. https://

Piperno, D.R., McMichael, C., Bush, M.B., 2015. Amazonia and the anthropocene:
what was the spatial extent and intensity of human landscape modification in
the amazon basin at the end of prehistory? Holocene 25, 1588e1597. https://

Platt Bradbury, J., 2000. Limnologic history of lago de patzcuaro, michoacan, Mexico
for the past 48,000 years: impacts of climate and man. Palaeogeogr. Palaeoclimatol.
Palaeoecol. 163, 69e95.

Pongratz, J., Reick, C., Raddatz, T., Claussen, M., 2008a. A reconstruction of global
agricultural areas and land cover for the last millennium. Glob. Biogeochem.
Cycles 22. n/ad-n/a.

Pongratz, J., Reick, C.H., Raddatz, T., Claussen, M., 2008b. A Global Land Cover
Reconstruction AD 800 to 1992 – Technical Description. Technical Report. Max
Planck Institute for Meteorology, Hamburg.

Pongratz, J., Caldeira, K., Reick, C.H., Claussen, M., 2011. Coupled climate-carbon
simulations indicate minor global effects of wars and epidemics on atmospheric
co2 between ad 800 and 1850. Holocene 21, 843e851.

Poorter, L., Bongers, F., Aide, T.M., Almeyda Zambrano, A.M., Balvanera, P.,
Becknell, J.M., Boukili, V., Brancalion, P.H.S., Broadbent, E.N., Chazdon, R.L.,
Craven, D., de Almeida-Cortez, J.S., Cabral, G.A.L., de Jong, B.H.J., Denslow, J.S.,
Dent, D.H., DeWalt, S.J., Dupuy, J.M., Duran, S.M., Espírito-Santo, M.M., 
Fandino, M.C., Cesar, R.G., Hall, J.S., Hernandez-Stefanoni, J.L., Jakovac, C.C., 
Junqueira, A.B., Kennard, D., Letcher, S.G., Licona, J.c., Lohbeck, M., MarínSpiotta,
E., Martínez-Ramos, M., Massoca, P., Meave, J.A., Mesquita, R., Mora, F.,
Munoz, R., Muscarella, R., Nunes, Y.R.F., Ochoa-Gaona, S., de Oliveira, A.A., Ori- ~
huela-Belmonte, E., Pena-Claros, M., P ~ erez-García, E.A., Piotto, D., Powers, J.S., 
Rodríguez-Velazquez, J., Romero-P  erez, I.E., Ruíz, J., Saldarriaga, J.G., Sanchez- 
Azofeifa, A., Schwartz, N.B., Steininger, M.K., Swenson, N.G., Toledo, M.,
Uriarte, M., van Breugel, M., van der Wal, H., Veloso, M.D.M., Vester, H.F.M.,
Vicentini, A., Vieira, I.C.G., Bentos, T.V., Williamson, G.B., Rozendaal, D.M.A.,
2016. Biomass resilience of neotropical secondary forests. Nature 1e15. https://

Power, M., Mayle, F., Bartlein, P., Marlon, J., Anderson, R., Behling, H., Brown, K.,
Carcaillet, C., Colombaroli, D., Gavin, D., Hallett, D., Horn, S., Kennedy, L., Lane, C.,
Long, C., Moreno, P., Paitre, C., Robinson, G., Taylor, Z., Walsh, M., 2013. Climatic
control of the biomass-burning decline in the americas after ad 1500. Holocene
23, 3e13.

Pregitzer, K.S., Euskirchen, E.S., 2004. Carbon cycling and storage in world forests:
biome patterns related to forest age. Glob. Chang. Biol. 10, 2052e2077. https://

Pyle, J., 1976. A reexamination of aboriginal population estimates for Argentina. In:
Denevan, W.M. (Ed.), Nativ. Popul. Am. 1492, 1 ed. University of Wisconsin Press,
Madison, Wis, pp. 181e204 (chapter 6).

Raghavan, M., Steinrucken, M., Harris, K., Schiffels, S., Rasmussen, S., DeGiorgio, M.,
Albrechtsen, A., Valdiosera, C., Avila-Arcos, M.C., Malaspinas, A.S., Eriksson, A.,
Moltke, I., Metspalu, M., Homburger, J.R., Wall, J., Cornejo, O.E., MorenoMayar,
J.V., Korneliussen, T.S., Pierre, T., Rasmussen, M., Campos, P.F.,
Damgaard, P.d.B., Allentoft, M.E., Lindo, J., Metspalu, E., Rodriguez-Varela, R.,
Mansilla, J., Henrickson, C., Seguin-Orlando, A., Malmstrom, H., Stafford, T.,
Shringarpure, S.S., Moreno-Estrada, A., Karmin, M., Tambets, K., Bergstrom, A.,
Xue, Y., Warmuth, V., Friend, A.D., Singarayer, J., Valdes, P., Balloux, F.,
Leboreiro, I., Vera, J.L., Rangel-Villalobos, H., Pettener, D., Luiselli, D., Davis, L.G.,
Heyer, E., Zollikofer, C.P.E., Ponce de Leon, M.S., Smith, C.I., Grimes, V., Pike, K.A.,
Deal, M., Fuller, B.T., Arriaza, B., Standen, V., Luz, M.F., Ricaut, F., Guidon, N.,
Osipova, L., Voevoda, M.I., Posukh, O.L., Balanovsky, O., Lavryashina, M.,
Bogunov, Y., Khusnutdinova, E., Gubina, M., Balanovska, E., Fedorova, S.,
Litvinov, S., Malyarchuk, B., Derenko, M., Mosher, M.J., Archer, D., Cybulski, J.,
Petzelt, B., Mitchell, J., Worl, R., Norman, P.J., Parham, P., Kemp, B.M., Kivisild, T.,
Tyler-Smith, C., Sandhu, M.S., Crawford, M., Villems, R., Smith, D.G.,
Waters, M.R., Goebel, T., Johnson, J.R., Malhi, R.S., Jakobsson, M., Meltzer, D.J.,
Manica, A., Durbin, R., Bustamante, C.D., Song, Y.S., Nielsen, R., Willerslev, E.,
2015. Genomic evidence for the pleistocene and recent population history of
native americans. Science 80 (349).

Read, L., Lawrence, D., 2003. Recovery of biomass following shifting cultivation in
dry tropical forests of the yucatan. Ecol. Appl. 13, 85e97.

Rhoades, C.C., Eckert, G.E., Coleman, D.C., 2009. Soil carbon differences among
forest, agriculture, and secondary vegetation in lower montane Ecuador. Ecol.
Appl. 10, 497e505.

Richards, J.F., 2006. The Unending Frontier. An Environmental History of the Early
Modern World, 1 ed. University of California Press, Oakland, CA.

Roosevelt, A.C., 2013. The amazon and the anthropocene: 13,000 years of human
influence in a tropical rainforest. Anthropocene 4, 69e87.

Roosevelt, A.C., Housley, R.A., Da Silveira, M.I., Maranca, S., Johnson, R., 1991. Eighth
millennium pottery from a prehistoric shell midden in the brazilian amazon.
Science 80 (254), 1621e1624.

Rosenblat, A., 1954. La poblacion indigena y el mestizaje en America, 1 ed. Editorial
Nova, Buenas Aires.

Rosenblat, A., 1976. The population of hispaniola at the time of columbus. In:
Denevan, W. (Ed.), Nativ. Popul. Am. 1492, 1 ed. University of Wisconsin Press,
Madison, Wis, pp. 43e66 (chapter 2).

Rostain, S., 2008. The archaeology of the guianas: an overview. Handb. South Am.
Archaeol. 279e302.

Rozendaal, D.M., Chazdon, R.L., Arreola-Villa, F., Balvanera, P., Bentos, T.V.,
Dupuy, J.M., Hernandez-Stefanoni, J.L., Jakovac, C.C., Lebrija-Trejos, E.E.,
Lohbeck, M., Martínez-Ramos, M., Massoca, P.E., Meave, J.A., Mesquita, R.C.,
Mora, F., Perez-García, E.A., Romero-P  erez, I.E., Saenz-Pedroza, I., van 
Breugel, M., Williamson, G.B., Bongers, F., 2017. Demographic drivers of
aboveground biomass dynamics during secondary succession in neotropical dry
and wet forests. Ecosystems 20, 340e353.

Rubino, M., Etheridge, D.M., Trudinger, C.M., Allison, C.E., Rayner, P.J., Enting, I.,
Mulvaney, R., Steele, L.P., Langenfelds, R.L., Sturges, W.T., Curran, M.A.J., Smith, A.M., 2016. Low atmospheric co2 levels during the little ice age due to
cooling-induced terrestrial uptake. Nat. Geosci. 9, 691e694.

Ruddiman, W.F., 2005. Climate and human history. In: Plows, Plagues, Pet. Princeton
University Press, p. 9.

Ruddiman, W.F., 2007. The Early Anthropogenic Hypothesis: Challenges and Responses.

Ruddiman, W.F., 2013. Earth Transformed, 2 ed. W. H. Freeman.

Ruddiman, W.F., Ellis, E.C., 2009. Effect of per-capita land use changes on holocene
forest clearance and co2 emissions. Quat. Sci. Rev. 28, 3011e3015. https://

Salzano, F., Callegari-Jacques, S., 1988. South American Indians: a Case Study in
Evolution. Clarendon Press, Oxford.

Sanders, W.T., 1972. Population, agricultural his- tory, and societal evolution in
mesoamerica. In: Spooner, B. (Ed.), Popul. Growth Anthropol. Implic, 8 ed. MIT
Press, Cambridge, pp. 101e153.

Sanders, W.T., 1976. The population of the central mexican symbiotic region, the
basin of Mexico, and the teotihuacan valley in the sixteenth century. In:
Denevan, W. (Ed.), Nativ. Popul. Am. 1492. University of Wisconsin Press,
Madison, Wis, pp. 85e151 (chapter 4),.

Sanders, W.T., Parsons, J.R., Santley, R.S., 1979. The Basin of Mexico: the Ecological
Processes in the Evolution of a Civilization. New York Academic Press, New
York, N.Y.

Sapper, K., 1924. Die zahl und die volksdichte der indianischen bevolkerung in der
zeit vor der conquista und in der gegenwart. In: 21st Int. Congr. Am., the Hague,
pp. 95e104.

Sauer, C., 1980. Seventeenth Century North America. University of California Press,

Schaan, D.P., 2008. The nonagricultural chiefdoms of marajo island. In: 

Silverman, H., Isbell, W. (Eds.), Handb. South Am. Archaeol. Springer New York,
New York, NY, pp. 339e357 (chapter 19).

Schacht, R.M., 1984. The contemporaneity problem. Am. Antiq. 49, 678e695.

Schroth, G., D’Angelo, S.A., Teixeira, W.G., Haag, D., Lieberei, R., 2002. Conversion of
secondary forest into agroforestry and monoculture plantations in amazonia:
consequences for biomass, litter and soil carbon stocks after 7 years. For. Ecol.
Manage. 163, 131e150.

Schurer, A.P., Hegerl, G.C., Mann, M.E., Tett, S.F.B., Phipps, S.J., 2013. Separating
forced from chaotic climate variability over the past millennium. J. Clim. 26,

Shea, D., 1976. A defense of small population estimates for the central andes in
1520. In: Dene (Ed.), Nativ. Popul. Am. 1492, 2 ed. University of Wisconsin Press,
Madison, Wis, pp. 157e180 (chapter 4).

Shi, S., Zhang, W., Zhang, P., Yu, Y., Ding, F., 2013. A synthesis of change in deep soil
organic carbon stores with afforestation of agricultural soils. For. Ecol. Manage.
296, 53e63.

Sigl, M., Winstrup, M., McConnell, J.R., Welten, K.C., Plunkett, G., Ludlow, F.,
Büntgen, U., Caffee, M., Chellman, N., Dahl-Jensen, D., Fischer, H., Kipfstuhl, S.,
Kostick, C., Maselli, O.J., Mekhaldi, F., Mulvaney, R., Muscheler, R., Pasteris, D.R.,
Pilcher, J.R., Salzer, M., Schüpbach, S., Steffensen, J.P., Vinther, B.M.,
Woodruff, T.E., 2015. Timing and climate forcing of volcanic eruptions for the
past 2,500 years. Nature 523, 543e549.
Silver, W.L., Ostertag, R., Lugo, A.E., 2000. The potential for carbon sequestration
through reforestation of abandoned tropical agricultural and pasture lands.
Restor. Ecol. 8, 394e407.

Smith, N.J., 1999. The Amazon River Forest: A Natural History of Plants, Animals and
People, 1 ed. Oxford University Press, New York, N.Y.

Smith, B.D., 2006. Eastern north America as an independent center of plant domestication. Proc. Natl. Acad. Sci. 103, 12223e12228.

Smith, P., Bustamante, M., Ahammad, H., Clark, H., Dong, H., Elsiddig, E., Haberl, H.,
Harper, R., House, J., Jafari, M., Masera, O., Mbow, C., Ravindranath, N., Rice, C.W.,
Robledo Abad, C., Romanovskaya, A., Sperling, F., Tubiello, F., 2014. Agriculture,
forestry and other land use (afolu). In: Edenhofer, O., Pichs-Madruga, R.,

Sokona, Y., Farahani, E., Kadner, S., Seyboth, K., Adler, A., Baum, I., Brunner, S.,
Eickemeier, P., Kriemann, B., Savolainen, J., Schlomer, S., von Stechow, C., €
Zwickel, T., Minx, J. (Eds.), Clim. Chang. 2014 Mitig. Clim. Chang. Contrib. Work.
Gr. III to Fifth Assess. Rep. Intergov. Panel Clim. Chang, vol. 1. Cambridge University
Press, Cambridge & New York, NY, pp. 7340e7349.
10.1104/pp.900074 (chapter 11).

Spracklen, D.V., Righelato, R., 2016. Carbon storage and sequestration of re-growing
montane forests in southern Ecuador. For. Ecol. Manage. 364, 139e144. https://

Stahl, P.W., 2015. Interpreting interfluvial landscape transformations in the precolumbian
amazon. Holocene.

Steinhilber, F., Beer, J., Frohlich, C., 2009. Total solar irradiance during the holocene. €
Geophys. Res. Lett. 36, L19704.

Stenborg, P., Schaan, D.P., Figueiredo, C.G., 2018. Contours of the past: lidar data
expands the limits of late pre-columbian human settlement in the santarem
region, lower amazon. J. Field Archaeol. 43, 44e57.

Stinchcomb, G.E., Messner, T.C., Driese, S.G., Nordt, L.C., Stewart, R.M., 2011. Precolonial
(a.d. 1100-1600) sedimentation related to prehistoric maize agriculture
and climate change in eastern north America. Geology 39, 363e366. https://

Stocker, B.D., Strassmann, K., Joos, F., 2011. Sensitivity of holocene atmospheric co2
and the modern carbon budget to early human land use: analyses with a
process-based model. Biogeosciences 8, 69e88.

Stocker, B.D., Yu, Z., Massa, C., Joos, F., 2017. Holocene peatland and ice-core data
constraints on the timing and magnitude of co2 emissions from past land use.
Proc. Natl. Acad. Sci. 114, 1492e1497.

Storey, R., 2012. Population decline during and after conquest. In: Nichols, D.L. (Ed.),
Oxford Handb. Mesoamerican Archaeol. Oxford University Press, Oxford,
pp. 908e915.
Tainter, J.A., 1988. Collapse of Complex Societies, 1 ed. Cambridge University Press,
Cambridge & New York, NY.

Ter Steege, H., Pitman, N.C., Sabatier, D., Baraloto, C., Salomao, R.P., Guevara, J.E.,
Phillips, O.L., Castilho, C.V., Magnusson, W.E., Molino, J.F., Monteagudo, A.,
Vargas, P.N., Montero, J.C., Feldpausch, T.R., Coronado, E.N., Killeen, T.J.,
Mostacedo, B., Vasquez, R., Assis, R.L., Terborgh, J., Wittmann, F., Andrade, A.,
Laurance, W.F., Laurance, S.G., Marimon, B.S., Marimon, B.H., Vieira, I.C.G.,
Amaral, I.L., Brienen, R., Castellanos, H., Lopez, D.C., Duivenvoorden, J.F., 
Mogollon, H.F., Matos, F.D.D.A., D  avila, N., García-Villacorta, R., Diaz, P.R.S., 
Costa, F., Emilio, T., Levis, C., Schietti, J., Souza, P., Alonso, A., Dallmeier, F.,
Montoya, A.J.D., Piedade, M.T.F., Araujo-Murakami, A., Arroyo, L., Gribel, R.,
Fine, P.V., Peres, C.A., Toledo, M., Aymard, C.,G.A., Baker, T.R., Ceron, C., Engel, J., 
Henkel, T.W., Maas, P., Petronelli, P., Stropp, J., Zartman, C.E., Daly, D., Neill, D.,
Silveira, M., Paredes, M.R., Chave, J., Lima Filho, D.D.A., Jørgensen, P.M.,
Fuentes, A., Schongart, J., Valverde, F.C., Di Fiore, A., Jimenez, E.M., Mora, M.C.P., €
Phillips, J.F., Rivas, G., Van Andel, T.R., Von Hildebrand, P., Hoffman, B., Zent, E.L.,
Malhi, Y., Prieto, A., Rudas, A., Ruschell, A.R., Silva, N., Vos, V., Zent, S.,
Oliveira, A.A., Schutz, A.C., Gonzales, T., Nascimento, M.T., Ramirez-Angulo, H.,
Sierra, R., Tirado, M., Medina, M.N.U., Van Der Heijden, G., Vela, C.I., Torre, E.V.,
Vriesendorp, C., Wang, O., Young, K.R., Baider, C., Balslev, H., Ferreira, C.,
Mesones, I., Torres-Lezama, A., Giraldo, L.E.U., Zagt, R., Alexiades, M.N.,
Hernandez, L., Huamantupa-Chuquimaco, I., Milliken, W., Cuenca, W.P.,
Pauletto, D., Sandoval, E.V., Gamarra, L.V., Dexter, K.G., Feeley, K., LopezGonzalez,
G., Silman, M.R., 2013. Hyperdominance in the amazonian tree flora.
Science 80 (342).

Thornton, R., 1987. American Indian Holocaust and Survival: A Population History
since 1492, 1 ed. University of Oklahoma Press, Norman, OK.

Trudinger, C.M., Enting, I.G., Francey, R.J., Etheridge, D.M., Rayner, P.J., 1999. Longterm
variability in the global carbon cycle inferred from a high-precision co2
and d13c ice-core record. Tellus 51B, 233e248.

Trudinger, C.M., Enting, I.G., Rayner, P.J., Francey, R.J., 2002. Kalman filter analysis of
ice core data 2. double deconvolution of co2 and ??13c measurements.
J. Geophys. Res. Atmos. 107, 1e24.

Turner, B.L., Brush, S., 1987. Comparative Farming Systems. New York. Guilford
Press, New York.

Ubelaker, D.H., 1976. The sources and methodolgy for mooney’s estimates of north
american indian populations. In: Denevan, W.M. (Ed.), Nativ. Popul. Am. 1492, 2
ed. University of Wisconsin Press, Madison, Wis, pp. 243e288 (chapter 8).

Ubelaker, D.H., 2006. Environment, origins, and population. In: Ubelaker, D.H. (Ed.),
Handb. North Am. Indians 3. Smithsonian Institution Press, Washington D.C.,
pp. 694e701

Urquiza-Haas, T., Dolman, P.M., Peres, C.A., 2007. Regional scale variation in forest
structure and biomass in the yucatan peninsula, Mexico: effects of forest
disturbance. For. Ecol. Manage. 247, 80e90.

Verlinden, C., 1973. La population de l’amerique precolumbienne: Une question de
method. Methodol. l’histoire des Sci. Hum. Melanges en Honn. Fernand Braud pp. 453e462.

Vieira, L.E.A., Solanki, S.K., Krivova, N.A., Usoskin, I., 2011. Evolution of the solar
irradiance during the holocene. Astron. Astrophys. 531, A6.

Villamarín, J., Villamarín, J., 2000. Chiefdoms: the prevalence and persistence of
”senorios naturales ~ ” 1400 to european conquest. In: Salomon, F., Schwartz, S.
(Eds.), Cambridge Hist. Nativ. Peoples Am. Vol. III South Am. Part 1. Cambridge
University Press, Cambridge, pp. 577e667 (chapter 7).

Walker, R.S., Sattenspiel, L., Hill, K.R., 2015. Mortality from contact-related epidemics
among indigenous populations in greater amazonia. Sci. Rep. 5, 14032.

Wang, S., Lewis, C.M., Jakobsson, M., Ramachandran, S., Ray, N., Bedoya, G.,
Rojas, W., Parra, M.V., Molina, J.A., Gallo, C., Mazzotti, G., Poletti, G., Hill, K.,
Hurtado, A.M., Labuda, D., Klitz, W., Barrantes, R., Bortolini, M.C., Salzano, F.M.,
Petzl-Erler, M.L., Tsuneto, L.T., Llop, E., Rothhammer, F., Excoffier, L.,
Feldman, M.W., Rosenberg, N.A., Ruiz-Linares, A., 2007. Genetic variation and
population structure in native americans. PLoS Genet. 3, 2049e2067. https://

Watling, J., Iriarte, J., Mayle, F.E., Schaan, D., Pessenda, L.C.R., Loader, N.J., StreetPerrott,
F.A., Dickau, R.E., Damasceno, A., Ranzi, A., 2017. Impact of precolumbian builders on amazonian forests. Proc. Natl. Acad.  Sci. 114, 1868e1873.

Watling, J., Shock, M.P., Mongelo, G.Z., Almeida, F.O., Kater, T., De Oliveira, P.E., 
Neves, E.G., 2018. Direct archaeological evidence for southwestern amazonia as
an early plant domestication and food production centre. PLoS One 13,

Watson, R., Noble, I., Bolin, B., Ravindranath, N.H., Verardo, J., Dokken, D., 2000.
Land use, land use change, and forestry. In: A Spec. Rep. IPCC. Cambridge
University Press, Cambridge, p. 30.

Whitmore, T.M., 1991. A simulation of the 16th century population collapse in the
basin of Mexico. Ann. Assoc. Am. Geogr. 81, 464e487.

Whitmore, T.M., Turner, B.L., 1992. Landscapes of cultivation in mesoamerica on the
eve of the conquest. Ann. Assoc. Am. Geogr. 82, 402e425.

Whitney, B.S., Dickau, R., Mayle, F.E., Soto, J.D., Iriarte, J., 2013. Pre-columbian
landscape impact and agriculture in the monumental mound region of the
llanos de moxos, lowland Bolivia. Quat. Res. (United States) 80, 207e217.

Whitney, B.S., Dickau, R., Mayle, F.E., Walker, J.H., Soto, J.D., Iriarte, J., 2014. Precolumbian
raised-field agriculture and land use in the bolivian amazon. Holocene
24, 231e241.

Williams, E., 1970. From Columbus to Castro: the History of the Caribbean 1492-
1969. Vintage, London.

Williams, C.A., Collatz, G.J., Masek, J., Goward, S.N., 2012. Carbon consequences of
forest disturbance and recovery across the conterminous United States. Glob.
Biogeochem. Cycles 26, 1e13.

Wolf, E.R., 1982. Europe and the People without History. University of California
Press, Berkeley.

Woods, W.I., Teixeira, W.G., Lehmann, J., Steiner, C., WinklerPrins, A., Rebellato, L.,
2009. Amazonian Dark Earths: Wim Sombroek’s Vision, Amazonian Dark
Earths: Wim Sombroek’s Vision. Springer Netherlands, Dordrecht. https://

Woods, W.I., Denevan, W.M., Rebellato, L., 2013. Population estimates for anthropogenically
enriched soils (amazonian dark earths). In: Soils, Clim. Soc. University
Press of Colorado, Boulder. Archaeological Investigations in Ancient
America, pp. 1e20 (chapter 1).

Woomer, P., Palm, C., Alegre, J., Castilla, C., Cordeiro, D., Hairiah, K., Kotto-Same, J.,
Moukam, A., Reise, A., Rodrigues, V., van Noordwijk, M., 1999. Slash-and-burn
effects on carbon stocks in the humid tropics. In: Lal, R., Kimble, J., Stewart, B.
(Eds.), Glob. Clim. Chang. Trop. Ecosyst. CRC Press, Boca-Raton, pp. 99e115.

Yan, P., Pan, J., Zhang, W., Shi, J., Chen, X., Cui, Z., 2017. A high plant density reduces
the ability of maize to use soil nitrogen. PLoS One 12, 1e12.

Zambardino, R.A., 1980. Mexico’s population in the sixteenth century: demographic
anomaly or mathematical illusion? J. Interdiscip. Hist. 11, 1e27.

Zamora, E., 1983. Conquista y crisis demografica: la poblacion indígena del occi- 
dente de Guatemala en el siglo xvi. Mesoamerica 6, 291e338.