Abstract
A significant gap in our understanding of global change involves the linkages between historical land-use and land-cover change (LULCC), Holocene terrestrial carbon (C) pool fluxes, and climate change. To address that research problem, this method uses land grants for sheep ranches awarded in early colonial (1521–1620 CE) Mexico to quantify the amount of land converted from cropland to pastureland. Soil is the largest terrestrial C pool, and converting cropland to pastureland significantly increases soil C sequestration rates, thereby reducing atmospheric C. The land grants and associated archival documents contain location-specific information about soils, vegetation, hydrology, and other variables that make it possible to map the ranches in a Geographic Information System (GIS) and quantify LULCCs and terrestrial C pools over time and space. The results demonstrate the utility of such archival documents to research on Holocene global change, indicate that LULCC during Mexico’s colonial period increased the region’s soil C sequestration rate, and confirm previous research that has suggested that transformations associated with European colonization of the Americas acted as an anthropogenic contributor to the period of moderate cooling known as the ‘Little Ice Age’ (1550–1850 CE).
Introduction
The influence of pre-industrial land-use and land-cover change (LULCC) on the global climate system is becoming more clearly understood (He et al., 2014; Ruddiman, 2003; Ruddiman et al., 2014). Many conclusions are based on global-scale models that calculate change in the extent of cropland and pasture using historically reconstructed population densities as the independent variable (Boyle et al., 2011; Kaplan et al., 2011; Klein Goldewijk et al., 2011; Pongratz et al., 2008). Wide disagreement among these models because of differing demographic assumptions and lack of agroecological specificity has resulted in more attention to empirical LULCC research (Ellis, 2011; Yan et al., 2013). One source of data has been sedimentary charcoal records from lake and wetland cores, but their interpretation has likewise led to disagreement. For example, some have argued that sedimentary charcoal records for the Neotropics demonstrate less forest clearance through burning as indigenous populations collapsed beginning in 1492, forest regeneration through the 19th century, and increased terrestrial carbon (C) sequestration that contributed significantly to the period of atmospheric cooling known as the ‘Little Ice Age’ (1550–1850 CE; Dull et al., 2010; Nevle and Bird, 2008; Nevle et al., 2011). Others, in contrast, have interpreted the sedimentary charcoal record as demonstrating that biomass burning and atmospheric C levels had already begun to decline before 1492, and cannot therefore be ascribed to European colonization of the Americas, and that the ‘Little Ice Age’ actually forced LULCC in the Americas (Power et al., 2012).
One approach to improving interpretation of both the models based on demographic data and empirical data such as sedimentary charcoal involves the use of LULCC data systematically collected by colonial officials over extensive areas and now preserved as land-grant documents in archives. This research uses land grants awarded in an area of central Mexico in the 16th and 17th centuries to reconstruct LULCC, to estimate associated soil C sequestration rates, and to relate them to climate change. The study area is the southeastern quarter of the Valle del Mezquital, a semi-arid tropical highland where New Spain’s viceroys awarded 47 sheep ranches between 1535 and 1610. Its colonial land grants and associated archival documents contain location-specific information about the sites of the grants that include eyewitness observations of soils, vegetation, hydrology, previous habitations, and other LULCC variables. Descriptions of landmarks and boundaries allow the land grants to be relocated and mapped in a Geographic Information System (GIS) to quantify LULCCs and terrestrial C pools over time and space.
Study area
The southeastern Valle del Mezquital is a semi-arid basin of ~700 km2 located in the highlands of central Mexico, just north of Mexico City (Figure 1; Hunter, 2010). Extreme land degradation, quite evident to even a casual observer, is a defining characteristic of the entire Valle del Mezquital (Figure 2). In some places, gullying and unchecked sheet erosion have turned once vegetated slopes into wastelands representative of major degradation of Earth’s 4.5 million km2 of tropical highlands (Hughes, 2001; Nyssen et al., 2009; Richards, 2003). The causes and timing of that soil erosion remain controversial, with some blaming precolonial agricultural practices (Cook, 1949), others overgrazing by sheep during the first century of the colonial period (1521–1620; Melville, 1994), others agricultural modernization during the late colonial and early postcolonial periods (Butzer, 1992, 1993, 1996), and yet others a more complex process related to the cessation of maintenance of extensive systems of precolonial terracing during the demographic collapse of the 16th century (Hunter, 2009).

Central Mexico and the Valle del Mezquital.

Soil pedestals, a deeply incised streambed, and arid-zone spiny vegetation near Tecozautla, Hidalgo, Mexico, July 2013.
By late precolonial times (1400–1521), the basin was densely settled, with irrigated fields along valley bottoms and semi-terraces stepping up hillsides (Gerhard, 1993; Hunter, 2009; Mastache and Cobean, 2001; Melville, 1994). Warfare, epidemic diseases, and colonization caused a demographic collapse similar to elsewhere in the Americas during the first century of colonization, with sudden population decline, abandonment of cropland, deterioration of canals and terraces, and awards of emptied lands to Spaniards and indigenous communities to use as sheep pasture (Hoekstra, 1992; Hunter, 2009, 2014; Prem, 1992). During the remainder of the colonial period (1621–1810) and the early postcolonial period (1810–1910), sheep ranching and agave cultivation became the principal land uses (De Mendizábal, 1947; Fournier-García and Mondragón, 2003; Hernández Mogica, 2000; Johnson, 1977; Konrad, 1980). The soil erosion that characterizes the Valle del Mezquital came to serve as an emblem of land degradation, rural poverty, and out-migration characteristic of so much of Mexico. Since early in the 20th century, untreated storm runoff and sewage water has flowed northward from Mexico City through a complex system of canals, tunnels, and reservoirs to irrigate roughly 90,000 ha of rehabilitated agricultural fields, albeit with health risks for the local population and urban markets that consume the produce (Luchano-Constantino et al., 2009; Ortega-Larrocea et al., 2001; Tortajada and Castelán, 2003).
Materials and methods
In 1521, Spanish Conquistadors captured the Aztec capital Tenochtitlán and established the colony of New Spain, roughly equivalent to present-day Mexico. In the 1530s, with demand growing for hides, tallow, meat, and wool, colonists and native communities alike increasingly sought titles for livestock ranches. Throughout much of the early colonial period, the viceroys of New Spain awarded land grants for various land uses: agriculture, lime kilns, and ranching. The colonial land-grant register is now preserved in the Archivo General de la Nación in Mexico City (Sluyter, 2002). The associated site descriptions allow many land grants to be relocated and mapped using GIS (Hunter, 2010). These maps then facilitate area calculations of LULCCs over time and space for an entire century.
In general, as detailed in the literature, the land grants contain the awardee’s name, the name of the jurisdiction where the ranch was to be sited, a description of the site intended to reduce ambiguity about its borders, and the calendar date the grant was notarized by the scribe in Mexico City (Hunter, 2009, 2014; Hunter and Sluyter, 2011; Sluyter, 1999, 2002). The land grants for sheep ranches also stipulate a north–south oriented square of 776 ha and a stocking rate of 2000 head as the minimum necessary to retain grant ownership (Hunter, 2010). The site descriptions vary greatly in detail and may reference conspicuous landscapes features such as large trees or waterholes, toponyms of nearby settlements or mountains, the names of adjacent landowners, roads and stock trails, streambeds and whether they carry water, and vegetation and soil characteristics.
Overall, mapping ranches distributed in land grants yields a locationally and temporally precise, spatially extensive, fairly contiguous set of eyewitness descriptions of an area’s landscape over the course of several decades. These maps temporally cover the critical transitional century from precolonial to colonial land uses (1520s–1620s) and the onset of the ‘Little Ice Age’. The coverage of the grants clearly must be biased in favor of areas supporting at least some grassland, that resource being what prospective ranchers sought. Yet, despite that caveat, the grants do cover extensive portions of Mexico and by definition cover the lands that would have been turned into pastureland. As another caveat, informal grazing might have impacted the vegetation of a particular site before a rancher obtained a grant to it, but the land inspector’s descriptions are contemporaneous with the major influx of livestock in any particular district within a region and thus describe its LULCC over the course of several decades of ranch expansion. As a general principle, then, if the earliest land grants in a district describe intact vegetation cover and land grants awarded several decades later that occupy the interstices among existing ranches in that same district continue to describe intact vegetation cover, the conclusion must be that introduced livestock did not degrade the land. If the later grants describe soil erosion and vegetation degradation in comparison with the earlier ones, the conclusion must be that livestock, perhaps in concert with climate changes and other factors, did contribute to that transformation. The role of climate versus humans and their livestock can be further distinguished by applying the method to a large enough area to discriminate between local-scale LULCC caused by ranching and regional-scale LULCC caused by climate change.
Through a combination of field and archival research, the ranches can be mapped in a GIS (Hunter, 2010, 2014). This approach permits the reconstruction of spatio-temporal patterns of LULCC that occurred in the study area during the first century of colonization as the native population collapsed and cropland was converted into pastureland, as also occurred in other parts of New Spain. The quantifiable LULCCs over time, together with the different rates of soil C sequestration between these types, permits the calculation of changes in the study area’s total soil C sequestration over nearly a century.
Soil is the largest terrestrial C pool, and numerous studies and meta-analyses confirm that converting cropland to pastureland significantly increases soil organic C (SOC) sequestration rates (Don et al., 2011; Guo and Gifford, 2002; Johnson, 1992; Post and Kwon, 2000). This increase occurs because after tillage ceases, SOC is not as readily oxidized and lost to the atmosphere through CO2 emissions. In Central Mexico’s tropical dry-forest biome, the total C pool (above- and below-ground biomass and soil) of a pasture is 31% larger than that of a cultivated field: 108.3 Megagrams (Mg) C ha−1 and 82.7 Mg C ha−1, respectively (Jaramillo et al., 2003; Ordóñez et al., 2008). This figure is only slightly outside the well-established, more general range of 25–30% SOC loss following cultivation (Houghton et al., 2012). Because precolonial farmers in Mexico did not use plows, such modern estimates might somewhat inflate the differences in SOC between cropland and rangeland during the 16th-century transition from precolonial to colonial agroecologies.
Results
In this initial test of that method, 44 land grants awarded between 1535 and 1610 were mapped in a GIS (Figure 3). These land grants established 47 sheep ranches (three of the grants awarded two ranches each). They were all awarded within six cabeceras, colonial administrative districts with identically named head towns: Apaxco, Hueypoxtla, Tequixquiac, Tezcatepec, Tlapanaloya, and Tuzantlalpa. In the GIS, each ranch is represented by a north–south oriented square of 776 ha, the legally stipulated shape and size for sheep ranches. Overlap occurs because of the rudimentary surveying of the time and because ranches were not fenced (Aguilar-Robledo, 2003; Hunter and Sluyter, 2011; see also Wong et al., 2012). Agricultural semi-terraces visible in orthophotographs were digitized as polygons in the GIS. These polygons were then draped over the ranch polygons to reveal the areas of LULCC. To better track and visualize LULCC over time, the 47 ranches are divided into five date ranges (gaps between ranges reflect that no known ranches were awarded during those intervals).

Screenshot of GIS showing overlap between sheep ranches and agricultural semi-terraces.
The study area’s reconstructed LULCC history reveals three zones of colonial rural land use. Agricultural production occupied a humid elevational zone along the footslopes. Ranching mainly occurred on steeper, higher slopes. An intermediate zone marks LULCC in which roughly 5950 ha of semi-terraces were converted to pastoralism through land granting of cropland putatively abandoned because of depopulation. Because a hectare of pasture sequesters 25.6 Mg more SOC than a hectare of cultivated land, this LULCC increased the study area’s total C pool by at least 152,300 Mg (Figure 4).

LULCC and increase in the study area’s soil C sequestration.
Discussion
The 152,300 Mg increase in the study area’s SOC sequestration may represent a considerable underestimate, as supported by three lines of evidence. First, at the time the Spaniards established their colony, the study area had a dispersed settlement pattern, which implies that farmers were living near their terraced fields to maintain them (Evans, 1985; Gerhard, 1993; Sanders et al., 1979). Today, hillsides with scrub-covered terrace embankments are interspersed with wide swaths of bare bedrock, the result of the colonial period’s demographic collapse, and erosion because of terrace abandonment, lack of maintenance, and degradation (Hunter, 2014). The semi-terraces visible on the orthophotos used to estimate the conversion of agricultural land to pasture are therefore but the remnants of a much more highly managed and extensive precolonial agricultural landscape. The C sequestered in the areas of eroded agricultural terraces would have been interred as the topsoil moved downslope and was buried under the mineral parent material that followed it soon after (Van Oost et al., 2007).
Second, just as there were more semi-terraces in the 16th century than can be mapped in the GIS, there might also have been more ranches. An unknown number of land grants might no longer be preserved in any archive, and some scholars believe that on average, at least 25% of New Spain’s land grants are missing (Hunter, 2010; Prem, 1992; Simpson, 1952); others disagree (Sluyter, 2002). If true for the study region, the occurrence of more land grants and more ranches could increase the area of overlap between pasture and agricultural land used to estimate the LULCC.
Finally, because ranch boundaries did not adhere to their legally prescribed shape and size, each ranch possessed a unique shape and size that generally exceeded the prescribed square of 776 ha (Hunter, 2010). Again, the greater area of pasture could increase the area of overlap between pasture and agricultural land used to estimate the LULCC.
Global-scale vegetation models that rely on historical population density estimates and assumed per capita land use have shown large discrepancies in their calculations of LULCC and its relation to Holocene C fluxes (Boyle et al., 2011; Kaplan et al., 2011; Klein Goldewijk et al., 2011; Olofsson and Hickler, 2008; Pongratz et al., 2008). This empirical method provides a more detailed and realistic analysis of C sequestration at the sub-regional scale that can be used to refine global-scale models. Although this study’s archival and cartographic method is labor and time intensive, a series of similar small-scale studies where appropriate land grant data exist would permit an informed estimate of LULCC and the historical C budgets for entire regions. Extending the finding from the study area to the entire Valle del Mezquital in simple linear fashion suggests that roughly 1.5 million Mg of terrestrial C was sequestered as a result of conversion of cropland to pastureland in the early colonial period. By using land grants to actually map LULCCs over the course of a century, this method also overcomes the problem of uncertain population density estimates and per capita land use, numbers that drive the global-scale models and explain their discrepancies (Yan et al., 2013).
Conclusion
This article contributes to a growing literature that seeks to demystify the social and environmental processes by which human societies of the past have contributed to global climate anomalies such as the ‘Little Ice Age’ through impacts on the C cycle. The method therefore potentially extends the use of such local, context-specific LULCC data further into the past, complementing sedimentary records and models based on demographic estimates. Up until this point, the major studies of this type have focused on the US Great Plains, have a time depth of ~140 years, and use detailed population and agricultural censuses to run biogeochemical simulation models of ecosystem change (Gutmann et al., 2005; Hartman et al., 2011; Parton et al., 2005, 2013). Although similar census data do not exist for colonial New Spain, this research project shows that insights into fluxes in terrestrial C stocks can nonetheless be reached using land grant data. The Spaniards awarded such land grants throughout their colonies, from California and Florida to Argentina (Sluyter, 2012). Additional LULCC reconstructions in those other areas would reduce uncertainties over pre-industrial, human-induced changes in the global C cycle (Klein Goldewijk and Verburg, 2013). These reconstructions would also clarify whether and how regional-scale LULCCs triggered global-scale climate teleconnections detectable by climate models (Pitman et al., 2009). Wider application of this paper’s method would therefore provide a fuller understanding of the ecological consequences of European colonization of the Americas over the past five centuries, its LULCC consequences, and its relationship to climate change.
Footnotes
Funding
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.
