1. Introduction
The Yalahau region of northern Quintana Roo, Mexico (Figure 1), is defined by a network of 174 freshwater wetlands beginning at the north coast and extending south for approximately 50 km (Weidie, 1985; Fedick and Taube, 1995). This extensive freshwater wetland system represents a water-rich, biologically diverse ecosystem that stands in sharp contrast to the generally arid northern Yucatán Peninsula. The Yalahau Regional Human Ecology Project (YRHEP) has conducted archaeological and ecological research in the region since 1993 (Fedick and Taube, 1995). Research activities of the YRHEP include: a regional archaeological survey; detailed mapping of selected archaeological sites; excavation of ancient residential trash deposits; analysis of floral, faunal, and geological resources utilized; soil studies; ancient and modern agricultural practices; economic activities of the historical period; and environmental reconstruction (e.g., Mathews, 1998; Fedick et al., 2000; Morrison, 2000; Andersen, 2001; Rissolo, 2001; Rissolo et al., 2005; Fedick and Mathews, 2005; Glover, 2006; Rissolo and Mathews, 2006; Wollwage, 2008; Sedov et al., 2007, 2008; Sorensen, 2010; Solleiro-Rebolledo et al., 2011; Glover, 2012; Wollwage et al., 2012; Chmilar, 2013; Leonard, 2013; Gust, 2016).
The Yalahau region is interesting for several reasons. First is the settlement history. Settlement studies in the Yalahau region indicate an atypical occupation history that saw initial settlement in the Middle Preclassic period (700-200 BC), widespread use of the region in the Late Preclassic/Early Classic period (ca. 100 BC - AD 350/450), near abandonment of the region for the subsequent 800 years, and a smaller-scale reoccupation during the Postclassic period (AD 1200-1521; Fedick and Mathews, 2005; Rissolo et al.,2005; Glover and Stanton, 2010; Glover, 2012). The long Classic period abandonment of the Yalahau region is opposite the trend in neighboring areas, which saw Classic period political and demographic expansion of northern lowland polities such as Cobá to the south (Folan et al., 1983, 2009; Robles Castellanos, 1990), Chikinchel to the west (Kepecs, 1998), and Ek Balám to the southwest (Bey et al., 1998).
Second is exploitation of wetlands. Initial research by Fedick et al. (2000), and more recent work by Leonard (2013) and Chmilar (2013), documented anthropogenic landscapes within the wetlands characterized by constructed rock alignments of various configurations in diverse physical settings within the wetland swales (Figures 2 and 3). Twenty-six of approximately 174 wetlands have been surveyed to date, and 20 of the surveyed wetlands contain 516 rock alignments (for full details of the wetland surveys see Fedick, 1998a; Fedick et al., 2000; Leonard, 2013). Rock alignments are typically comprised of a single row of uncut limestone boulders and slabs; they range in length from 2 to 700 m (average 27 m) and average 0.46 m wide by 0.25 m tall. The alignments were originally suggested to have functioned in ancient times to modify soil and hydrologic conditions, in order to increase the production of useful wetland resources or for agricultural purposes (Fedick, 1998b; Fedick et al., 2000; Andersen, 2001). Evidence for ancient Maya wetland agriculture is well known from northern Belize (Turner and Harrison, 1983; Lambert et al., 1984; Hammond et al., 1987; Pohl, 1990; Jacob, 1995; Berry and McAnany, 2007; Beach et al., 2009), southern Quintana Roo (Gliessman et al.,1983), and Campeche (Siemens and Puleston, 1972) but usually takes the form of canals and raised fields. The use of rock alignments in the Yalahau region therefore is a different type of wetland agro-engineering, and along with rock alignments/walkways recorded by Hixon (2011), is the only documented evidence for wetland manipulation in the northern Maya Lowlands.
Third is the issue of environmental change. There is abundant evidence from lake core studies and geoarchaeological excavations that the Maya were confronted with short and long term environmental change, as well as their own impacts on the landscape (Beach et al., 2015; Dunning et al., 2015). For example, dry periods at the end of the Preclassic period are said to contribute to sociopolitical reorganization of this time (Carrillo-Bastos et al., 2010; Kennett et al., 2012; Douglas et al., 2015;). Likewise, a series of droughts in the Terminal Classic coincides with abandonment of numerous large urban/civic centers (Hodell et al., 1995, 2001, 2005, 2007; Haug et al., 2003; Hodell et al., 2005, 2007; Webster et al., 2007; Medina-Elizalde et al., 2010;). In low lying coastal plain areas where the water table floats on a sea water intrusion (Perry et al., 2003), rising Holocene sea level (Fairbanks, 1989; Lambeck and Chappell, 2001; Toscano and Macintyre, 2003) pushed up groundwater levels and caused increased flooding and soil aggradation in wetlands (Alcala-Herrera et al., 1994; Pope et al., 1996; Beach et al., 2009). In at least one case (Albion Island, Belize), continuing water table rise was cited as the cause for abandonment of farming in the Late Preclassic period (Pohl et al.,1990). “Maya clays” that accumulated in basins and lakes largely in the Preclassic period have been interpreted as evidence that a relatively small number of early farmers were responsible for a considerable amount of upland erosion due to deforestation, and the creation of anthropogenic soils within wetlands (Curtis et al., 1998; Hansen et al., 2002; Brenner et al., 2003; Beach et al., 2006; Dunning et al., 2006; Anselmetti et al., 2007).
In the Yalahau region, Fedick et al. (2000) suggested rising water levels may have impacted use of wetlands, and this in turn may have been a contributing factor in the Early Classic abandonment of the region (see also Wollwage, 2008; Wollwage et al., 2012; Fedick, 2014). This is supported by geoarchaeological investigations at the El Edén wetland (Figure 1), where one profile near the wetland edge shows modern wetland soil formation atop earlier dryland soils, indicating recent lateral expansion of this wetland (Solleiro-Rebolledo et al., 2011). Near the Yalahau coastal archaeological site of Vista Alegre, shallow water sediment cores indicated four sea level high stands that may have drowned freshwater springs along the coast and would have increased flood levels at inland wetlands (Beddows et al., 2016).
The El Edén reconstruction was based on a model of wetland soil development, initially described by Hodell et al. (2007: Figure 12), in which alternating bands of micritic muds and organic matter reflect periods of greater and lesser flooding (Solleiro-Rebolledo et al., 2011). The soil horizons dominated by carbonates formed over periods (decade to century scale), probably during the Holocene (although no instrumental dating is available), when wetlands experienced regular annual flooding as they do today - flooded with 1-2 m of water for ~7 months out of the year (Figure 4). These micritic muds are a product of algae (periphyton), which form extensive mats when the wetland is flooded, and which precipitate micrite out of the carbonate-enriched floodwaters (Sedov et al.,2007, 2008). Organic-dominated horizons, alternatively, formed during periods of less flooding (Hodell et al., 2007; Solleiro-Rebolledo et al., 2011). In these environmental conditions, algal activity is suppressed and the soil formation process is controlled more by the growth, death, and decomposition of rooting plants (Hodell et al., 2007).
It is in this context of complex interplay of human- and natural-induced environmental impacts that this study considers the use and abandonment of wetland farming in the Yalahau region. Of particular interest is how changing flood patterns in the wetlands relate to the Early Classic abandonment of the Yalahau region. The aim of this work is to use the El Edén soil sequence, and the model used to explain it (Hodell et al., 2007: Figure 12; Solleiro-Rebolledo et al., 2011; see Figure 4), as the basis for interpreting the timing and nature of changing hydrologic conditions at other wetlands in the region. This in turn will provide a context for understanding Maya use of wetlands and cycles of settlement and abandonment of the region. Five wetland soil-sediment profiles from across the region are studied - with analysis focused on soil micromorphology.
2. Materials and Methods
2.1 Environmental setting
The Yalahau region is located in northern Quintana Roo, Mexico (Figure 1), and covers approximately 3000 km2. The region is defined by a large network of freshwater wetlands that form within long, north-south trending depressions beginning at the north coast and continuing for at least 50 km south (Weidie, 1985).
The limestone platform of the Yucatan Peninsula consists of Early Cretaceous deposits, Middle Cretaceous deposits in the central and northern portions, and Cenozoic limestones up to 1200 meters thick (Weidie, 1985). The heavily karstified surface geology of the Yalahau region, and of the northern and eastern margin of the Yucatan Peninsula in general, consists of Miocene/Pliocene deposits (upper Tertiary) and coastal Holocene deposits of the Carillo Puerto Formation (Lauderdale et al., 1979; Rebolledo-Vieyra et al.,2000; Purdy et al., 2003; SGM, 2007; Gondwe, 2010). The Carillo Puerto Formation averages about 13 meters in thickness, and consists of poorly consolidated coquina, interbedded with calcarenite, capped with a caliche crust (Rebolledo-Vieyra et al., 2000).
The wetlands that characterize the Yalahau region are surface dissolution features reflecting the structure of the underlying Holbox Fracture Zone (Tulaczyk, 1993). This fracture zone is of undefined origin, but the faults are thought to create exceptionally permeable features that act as channels for groundwater movement (Weidie, 1985; Perry et al., 2003;). The approximately 174 wetlands range in size from about 0.005 km2 to 8 km2 (three wetlands exceed this, the largest, an extremely long wetland, being 31 km2) that cover an area of about 134 km2 (Fedick and Mathews, 2005). The region also contains over 100 cenotes, limestone sinkholes that typically reach the water table (Bell, 1998; Glover, 2006). Water found in wetlands andcenotesrepresents exposure of the water table (Tulaczyk, 1993). The aquifers of the northern Yucatan Peninsula have developed in the highly permeable rock dominated by Tertiary limestones and dolostones (Perry et al., 2003). The freshwater aquifers of the northern peninsula are underlain by a seawater intrusion, and as a result, the surface elevation of the water table is primarily controlled by sea level (Perry et al., 2003).
Soils of the Yalahau region include Rendzinc Leptosols in upland areas and Leptic Calcisols in wetlands (Sedov et al., 2008; Solleiro-Rebolledo et al., 2011). The wetland Calcisols are shallow and include substantial amounts of secondary carbonates, contrary to the usual image of wetlands containing thick deposits of organic peats and gleys. The relief is flat with altitudes <20 m asl (Lugo-Hubp et al., 1992), although there is an absolute difference in elevation of 1.93 m (Charvet, 2009). Average annual rainfall is 1295.8 mm, about 75% of which falls during a rainy season from May to October (Comisión Nacional del Agua, 2008). During the dry season from November to April, wetlands appear as open grasslands - soil may be moist but water is generally restricted to small ponds and lagoons within wetlands. During the rainy season, wetlands typically flood with 1-2 m of water. The upland forests of the region are characterized as dry tropical medium-statured, with notable species includingchico zapote(Manilkara zapota), black chechem (Metopium brownei), gumbo limbo (Bursera simaruba), and guano palm (Thrinax radiata). Wetlands include both swamp forests dominated palo tinto (Haematoxylum campechianum),and large open areas dominated by sedges such as sawgrass (Cladium jamaicense) and spikerushes (Eleocharissp.) (Schultz, 2003).
2.2 Fieldwork
Over the course of fieldwork between 2010 and 2011, 30 soil profiles were sampled across 12 wetlands. From these 30 profiles, 5 were selected that represent the following criteria: the range of thicknesses of the soil-sediment layer, the range of types of surface soil horizons and buried soil horizons, and a variety of landscape positions within wetlands (some of which are spatially related to rock alignments). The 5 profiles included in this study come from four wetlands: Zanja (2 profiles), Esperanza, Media Luna, and Pichol (Figure 5). In the case of Zanja, Esperanza, and Media Luna, sediment cores were extracted using 7.62 cm PVC pipe manually pushed into the ground. At Pichol, a soil pit was excavated and stratigraphy was described and sampled in situ.
2.3 Profile locations
The Zanja SP-1 core (87°17’50”W, 21°24’24”N) was retrieved from a relatively low elevation position near the center of the wetland dominated by sawgrass vegetation (Figure 1). Water was encountered at 0.10 m below surface and soil depth to bedrock was greater than 6 m (measured using a soil probe). The core was inserted to a depth of 1.6 m below surface and the compacted core length was 0.98 m. Periphyton was not observed on the surface here. Only one rock alignment was identified at the Zanja wetland, and it is located nearly 20 km south of the SP-1 core location.
The Zanja SP-3 core (87°17’55”W, 21°23’58”N) was retrieved from a relatively high elevation position at the west edge of the wetland with calabash (Crescentiasp.) and tasiste (Acoelorraphe wrightii) vegetation (Figure 1). Water was encountered at 0.18 m below surface and soil depth to bedrock was 0.38 m. The core was inserted to bedrock with no soil compaction. Periphyton was not observed on the surface here. Only one rock alignment was identified at the Zanja wetland, and it is located nearly 20 km south of the SP-3 core location.
The Esperanza SP-2 core (87°25’05”W, 20°58’44”N) was retrieved from a relatively low elevation position in the wetland dominated byScirpussp. vegetation (Figure 1). Water was encountered at 0.55 m below surface and soil depth to bedrock was 2.1 m. The core was inserted to a depth of 1.2 m and the compacted core length was 1.1 m. Periphyton was present on the surface. Three rock alignments were recorded at Esperanza. Alignment RA-1 is located 30 m north of the SP-2 core.
The Media Luna SP-1 core (87°25’05”W, 21°14’54”N) was retrieved from a relatively low-intermediate elevation position dominated by spikerush (Eleocharissp.) vegetation (Figure 1). Water was encountered at 0.15 m below surface and soil depth to bedrock was 1.21 m. The core was inserted to a depth of 1 m and the compacted core length was 0.67 m. Periphyton was present on the surface. Sixty-two rock alignments were recorded at Media Luna. The closest alignment to the SP-1 core is 25 m away, with many other rock alignments in the vicinity (Figure 6). The general distribution of rock alignments in this wetland shows alignments located near the margins of the open wetland, frequently connecting groups of topographic rises.
The Pichol SP-2 profile pit (87°24’00”W, 21°09’56”N) was excavated in a relatively intermediate elevation position in the wetland with sedge (Scirpussp.) vegetation (Figure 1). Water was not encountered and soil depth to bedrock was 0.52 m. The soil pit was excavated down to bedrock. Periphyton was present on the surface. Fifty-four rock alignments were recorded at Pichol. SP-2 was excavated 5 m from rock alignment RA-32, and numerous other alignments are in the vicinity. Similar to Media Luna, rock alignments at Pichol are found near the margins of the open wetland, connecting small bedrock islands.
2.4 Laboratory analyses
The cores and field samples were sent to the labs of Dr. Sedov and Dr. Solleiro-Rebolledo, Institute of Geology at the Universidad Nacional Autónoma de México, where they were opened, photographed, described, and sampled. Immediately upon opening the cores, it was clear the soils represented a highly heterogeneous system. The horizons were separated and designated according to the IUSS Working Group (2015).
Morphological and micromorphological analyses were applied to all cores. Thin sections were prepared from undisturbed samples taken at the time cores were opened, then observed and described with an Olympus petrographic microscope. Plane polarized light, crossed polarized light, and occasionally reflected light, was used to observe micromorphological features in the thin sections. Three cores (Zanja SP-1, Esperanza SP-2, and Media Luna SP-1) were selected for instrumental dating. Ten radiocarbon dates were run on bulk organic sediment from these 3 cores (Figure 7,Table 1). The radiocarbon samples were processed at the University of Arizona Accelerator Mass Spectrometry Laboratory and the University of California Irvine W. M. Keck Carbon Cycle Accelerator Mass Spectrometry (AMS) Laboratory. Radiocarbon date calibration was computed using the IntCal09 dataset and CABLIB 6.0 software package (Stuiver et al., 2013).
3. Results
3.1 Macromorphology
Macromorphologically, all profiles can be divided into an upper and lower unit (Figures 8 and 9). The upper units are fairly similar to one another, while the lower units are quite variable. The upper units of Zanja SP-1 and Esperanza SP-2 are characterized by alternating layers of whitish carbonate material and dark material enriched in organic matter. Shells are common in these layers and even the organic layers incorporate carbonates. The upper units of Media Luna SP-1, Pichol SP-2, and Zanja SP-3 display both humified and poorly decomposed organic matter mixed with carbonates and some shells.
In Media Luna, Pichol, and Zanja SP-3, lower unit soils exhibit clayey texture, little to no carbonates, and well-decomposed organic matter (proto-humus in some cases). The Zanja SP-1 lower unit has a slight reaction to HCl, some shells, and organic matter is mostly proto-humus. The structure is massive, breaking into angular blocks. The Esperanza SP-2 lower unit is white in color, has shells and a strong reaction to HCl, and also has a massive structure that breaks into angular blocks. Macroscopically, therefore, the lower and upper units of Zanja SP-1 and Esperanza SP-2 show similar morphology.
According to the IUSS Working Group WRB (2015), a horizon designation was established in each unit (Figure 5).
The upper unit of Media Luna SP-1 profile has O1 (0-12 cm), O2 (12-21 cm), and A (21-27 cm) horizons; while the lower unit shows the following horizons: 2A (27-39 cm), 2Bg (39-51 cm), 2Btg (51-66 cm).
Pichol SP-2 profile is constituted by A (0-10 cm, the upper unit), 2B1 (10-20 cm), 2B2 (20- 32 cm), and 2B3 (32-55 cm) horizons.
The Zanja SP-3 upper unit has BA and Bk horizons (0-4 cm, 4-11 cm, respectively). The lower unit exhibits 2A (11-28 cm) and 2C (>28 cm) horizons.
Esperanza SP-2 has in the upper unit A (0-19 cm), Bk (19-22 cm), 2A (22-30 cm), 2Bk (30-61 cm), a transitional layer (61-71 cm), 3AB (71-100 cm), and 4A (100-105 cm) horizons. The lower unit is only constituted by the 4Bk horizon (>105 cm).
The Zanja SP 1 upper unit has the following horizons: Bk (0-18 cm), 2A (18-34 cm), 2Bk (34-46), 3A (46-48 cm), 3Bk (48-58 cm), 4A (58-61 cm), 4Bk (61-66 cm), 5A (66-68 cm), 5Bk (68-79 cm), 6A (79-84 cm), 6Bk (84-85 cm). The lower unit of Zanja SP-1 is represented by the 7A and 7Bk horizons at a depth of 89-96 cm below surface.
3.2 Micromorphology of the upper units
The main micromorphological components of the upper units are biogenic micrite, sparite, freshwater mollusk shells, and organic matter in the form of partly decomposed fragments of plant tissues (colloidal organic material) (Figure 10a, b). They are found in the upper units of all cores, however, proportions vary. In the organic horizons, for instance in the O1 and O2 horizons of the Media Luna SP-1 profile, plant residues with intact cellular structure are observable; some pores of these residues are filled with micrite (Figures 10c and 10d). In the calcic horizons (e.g., 3Bk horizon of the Zanja SP-1) complete shells are found in a micritic groundmass locally containing dark organic pigment (Figures 10e and 10f). These components reflect the modern seasonal flood regime.
3.3.1 Media Luna SP-1
The 2A horizon is characterized by organic matter and clay, with no observable carbonates (Figure 11a). The groundmass is very fine, with some relatively large charcoal particles. Frequent infillings of tabular gypsum crystals were observed in cracks around the groundmass (Figures 11b and 11c). Localization in the pores (while unobserved in the groundmass inside the peds), as well as perfect crystalline shape, point to gypsum neoformation.
Development of gypsic pore pedofeatures is frequently associated with evaporation of gypsum-saturated groundwater (Poch et al., 2010). The key point is that gypsum is 10 times more soluble than carbonates. Thus, in these environmental conditions it is very rare for gypsum to remain undissolved in crystal form. The presence of gypsum therefore may indicate drier past conditions. However, infillings of neoformed gypsum are also known to be the oxidation product of pyrite (Miedema et al., 1974). In this case it could be related to the aeration of swampy soils where pyrite was formed earlier under anoxic conditions due to fluctuation of the water table. This latter scenario does not necessarily imply strong aridization. The discovery of framboidal pyrite in the lower unit of Esperanza SP-2 supports this version.
The 2Bg and 2Btg horizons are characterized by a clay groundmass with a patchy color pattern; pale clayey areas alternating with rusty brown ferruginous mottles points to the redistribution of iron, due to gleyzation processes (Figure 11d).
3.3.2 Zanja SP-1
The 7A horizon shows a high concentration of organic matter, a large part of which has been converted to “protohumus” - dark collomorphic material showing no or weak remains of cellular structure (although some pieces have tissue morphology). The groundmass includes areas impregnated with organic matter with frequent small cracks and pores, as well as a few large cracks. Other areas of the groundmass with more carbonates show medium-compact blocks separated by fractures. Unlike in the upper units, the calcite is mostly sparitic, and shell content is high. The interpretation of the lower unit of Zanja SP-1 is that organic matter had time for considerable decomposition - this implies at least some periods of an oxidative soil environment and limited flooding.
3.3.3 Zanja SP-3
Although the 2A horizon has less carbonates compared with the overlying Bk horizon (see upper unit description), there is much more organic matter here, and it is highly decomposed. In some places, there are sparitic carbonate nodules or larger crystals as well as micritic nodules. The most important soil feature in the lower unit here, however, is the area with large anhedral calcite crystals forming compact infillings (Figures 12a and 12b). Solleiro-Rebolledo et al. (2011) published similar results from peripheral polygenic profiles at El Edén where carbonates were found filling pores and cementing aggregates enriched with humus. Such coarse-grained cement of interlocking calcite crystals is typical for the phreatic environment (Durand et al., 2010). These carbonates are very different from those formed by wetland floods with algae - instead, they grow slowly under permanent saturation by groundwater.
3.3.4 Esperanza SP-2
In the 4Bk horizon micrite was not observed; carbonate is mostly in the form of sparitic calcite crystals (sand-size) that are well spread out and evenly distributed. The crystals have a well-defined shape, indicating they were formedin situ. Some crystals display etched pits - a dissolution feature. This indicates a highly dynamic horizon where calcite is neoformed but also dissolved (after being neoformed/precipitated). Mollusk shells are present but most are fragmented. Some organic matter is present, appearing as small dark dots (humus punctuations). A few loose infillings of black globular particles that are orange and gold under reflected light represent neoformed pyrite (framboidal crystals), pointing to an anoxic, permanently water-saturated environment (Figures 12c and 12d).
3.3.5 Pichol SP-2
The 2B1 horizon is characterized by nearly pure clay material, pale yellow in color - depleted of ferruginous pigment - due to gleying. The groundmass is compact in general, but some pores and channels occur, as do a few excrements of mesofauna (species undetermined). Some humus pigmentation is mixed with clay, and some decomposed organic particles (nearly humus) are incorporated into the groundmass. In some places, there are large ferruginous nodules, indicating redistribution of iron. A large round iron oxide nodule with clear concentric growth rings indicates that it is growing in situ (Figures 12e and 12f). Rounded ferruginous concretions are a common feature of the lower unit horizons of Pichol SP-2.
Horizon 2B2 is compact and clayey, with some brown ferruginous pigment - fine iron oxides mixed with clay result in this color. There are also signs of gleyzation - some iron nodules and areas with paler color. However, it was not saturated permanently since the horizon is only partially gleyed. No carbonates (including shells) or gypsum were observed in the 2B2 horizon.
Lab ID | Depth (cm) | Material | δ13C | 14 C Age BP | 14 C Age BP | 2σ Calibrated Age Ranges |
---|---|---|---|---|---|---|
Zanja Core 1 | ||||||
AA 92181 | 24 | org-rich soil | -27.2 | 38 ± 39 | 51 ± 82 | AD 1692-1728 (p = .231) AD 1811-1920 (p = .725) AD 1952- 1955 (p = .045) |
UCIAMS 104165 | 47 | org-rich soil | * | 1550 ± 15** | 1457 ± 42** | AD 432-495 (p = .605) AD 503-560 (p = .395) |
AA 92182 | 59 | org-rich soil | -27.8 | 1137 ± 35 | 1050 ± 53 | AD 781-791 (p = .026) AD 807-986 (p = .974) |
UCIAMS 104167 | 67 | org-rich soil | * | 1375 ± 15 | 1375 ± 15 | AD 643-667 (p = 1) |
UCIAMS 104168 | 79 | org-rich soil | * | 1695 ± 15 | 1595 ± 23 | AD 260-283 (p = .129) AD 323-407 (p = .871) |
AA 92183 | 88 | org-rich soil | -27.6 | 2098 ± 39 | 2073 ± 51 | 344-323 BC (p = .025) 205-37 BC (p = .955) 30-21 BC (p = .009) 11-2 BC (p = .011) |
Esperanza Core 2 | ||||||
UCIAMS 104164 | 30 | org-rich soil | * | 1175 ± 15 | 1107 ± 33 | AD 779-794 (p = .106) AD 798-894 (p = .886) AD 929-932 (p = .007) |
UCIAMS 104166 | 74 | org-rich soil | * | 2630 ± 15 | 2757 ± 5 | 815-793 BC (p = 1) |
AA 92179 | 99 | org-rich soil | -26.9 | 3115 ± 40 | 3332 ± 48 | 1493-1476 BC (p = .025) 1460-1292 BC (p = .968) 1277-1272 BC (p = .007) |
Media Luna Core 1 | ||||||
AA 92180 | 37 | org-rich soil | -28.7 | 1706 ± 37 | 1627 ± 54 | AD 248-413 (p = 1) |
*Dates are corrected for isotopic fractionation but δ13C values are not reported to clients.
**Date represents an age reversal and is excluded from analysis.
Horizon 2B3 is similar to the above, but the gleyzation is stronger and ferruginous nodules and concretions are more abundant. This horizon also lacks micritic carbonates and shells, similar to those horizons observed in the well-drained soils in the uplands of the adjacent forest, as reported by Sedov et al. (2007, 2008). Therefore, this soil layer formed under semi-hydromorphic, semi-swampy conditions (while the soil in Esperanza SP-2 [4Bk horizon] with pyrite was permanently saturated).
3.4 Chronological frame of the studied profiles
Ten radiocarbon dates were collected from 3 different cores that span the past 3332 years cal BP (Table 1; see Figure 7). Due to the absence of identifiable charcoal or other reliably ancient organic materials, all dates were run on bulk organic sediment. The Zanja core is the most well-stratified sequence, with 6 dates ranging from 2073 ± 51 cal BP ([AA 92183] 205-37/344-323/11-2/30-21 BC) at the base of the core to 51 ± 82 cal BP ([AA 92181] AD 1811-1920/1692-1728/1952-1955) near the top. Aside from the date of 1457 ± 42 cal BP at a depth of 47 cm ([UCIAMS 104165] AD 432-495/503-560), the radiocarbon dates from Zanja SP-1 are in chronological order. Three dates from Esperanza SP-2 range from 3332 ± 48 cal BP ([AA 92179] 1460-1292/1493-1476/1277-1272 BC) to 1107 ± 33 cal BP ([UCIAMS 104164] AD 798-894/779-794/929-932) and are in chronological order. A single determination [AA 92180] from Media Luna SP-1 dates 1627 ± 54 cal BP (AD 248-413). We acknowledge that relying on a single date from Media Luna is problematic, and recognize there is a need for further dating of the site.
4. Discussion
4.1 Upper vs. lower units
The macroscopic and microscopic results clearly show the separation of the cores into two major units (see Figures 8 and 9). Similar to the soils studied at El Edén (Solleiro-Rebolledo et al., 2011), the upper units show an alteration of fine carbonates and organic materials, as well as shell. Under the microscope, the fine carbonates are mostly micritic, and the organic matter is mainly detritus (intact plant fragments) with some humified parts or protohumus substance. As described by Hodell et al. (2007) and Solleiro-Rebolledo et al. (2011), the micritic muds identify flooding events because algae responsible for micrite precipitation only develop on the water surface. Organic materials, on the other hand, mostly come from higher/rooting plants, which may tolerate some flooding, but mostly develop under conditions of less flooding.
The lower units demonstrate a great variety of features and formation processes (despite appearing similar macromorphologically). Micromorphological examination of soil thin sections identified gley properties, highly organic materials with coarse-grained calcite, gypsum neoformation, and in one case, in the 4Bk horizon of Esperanza SP-2, precipitation of pyrite ( Figure 12c and 12d). None of the lower units have micritic muds, and it is therefore concluded that they were not flooded for any significant periods of time. All of the lower units, however, show signs of water saturation and swampy pedogenesis. There is evidence of gleyzation with the accumulation of organic matter (mostly well decomposed) that is partly mixed with clay. Also present are neoformed minerals that also show a hydrogenic origin. Pyrite can be formed only in a very anoxic environment, thus in water-saturated (i.e. reduced) contexts. If exposed to oxygen, pyrite forms iron oxide, and if it interacts with sulfuric acid, it can form gypsum (Miedema et al.,1974; Poch et al., 2010).
Coarse-grained carbonate crystals are typical of hydrogenic carbonates precipitated from groundwater. They are not biogenic, and grow slowly under permanent saturation by groundwater. Conditions in the wetland during formation of the lower unit therefore appear to have been swampy, but with limited flooding. This formation process occurs similarly in the gleyic horizons (2Bg and 2Btg horizons) seen in the lower unit at Media Luna.
The presence of gypsum is particularly interesting. Gypsum is typically associated with saline groundwater, but in this case, since pyrite was identified elsewhere, it could be that gypsum appeared due to pyrite oxidation and further interaction of sulfuric acid with carbonates. The deep pyrite in Esperanza SP-2 (4Bk horizon) was permanently saturated, and therefore did not oxidize. In shallower horizons (such as Media Luna 2A horizon), pyrite could have been present, oxidized, and all that remains visible now is gypsum.
4.2 Initial interpretation correlating Yalahau settlement history, wetland modification, and paleoenvironmental sequences
As discussed above, the lower soil units reflect a variety of ancient soils, which have been buried by the more recent, and relatively uniform, carbonate-enriched flood-water derived upper units. It should be noted that although only 5 profiles were studied in detail, 25 others show an even greater variety of pedogenic features. Many of these 25 profiles have organic strata beneath carbonate layers, indicating that numerous other wetlands were swampy but unflooded in the past. However, as indicated by the radiocarbon dates (see Figure 7, Table 1), the transitions between the lower and upper units did not occur at the same time. In Esperanza, the transition occurred in the Early Preclassic period (prior to 1460 to 1292 BC) while in Zanja SP-1 it occurred at the end of the Late Preclassic period (between 205-37 BC and 323-407 AD), and in Media Luna SP-1 it occurred early in the Classic period (after 248-413 AD).
Regarding the presence of rock alignments in the studied wetlands and association with the studied cores, Media Luna and Pichol contain far more rock alignments than Zanja and Esperanza (1 and 3 respectively). While no stratigraphic relationship has been established and alignments have not been dated, we suspect that the concentration of alignments at Media Luna and Pichol relates to a later onset of flooding than in wetlands such as Esperanza, where flooding began in the Early Preclassic. In Media Luna and Pichol, cores were extracted in close proximity to rock alignments, and it is supposed that alignments are stratigraphically associated with lower unit soil cover when at least the margins of the wetland rarely flooded, and alignments in lower areas of these wetlands may have been buried by upper unit soils, as flooding increased through the late Holocene and Calcisols developed. In other wetlands not studied here, rock alignments are common at wetlands with lower unit organic strata.
In the context of Maya settlement history and agricultural exploitation of wetlands, it is suggested, as a most likely scenario, that the absence of surface flooding in certain wetlands during the Preclassic period was advantageous to the Yalahau Maya. In the Preclassic period, many wetlands such as Media Luna and Pichol were more like broad, soil-filled depressions that were excellent for farming during the dry season. It is suggested that the Maya constructed the majority of rock alignments during this time, primarily in swampy but unflooded wetlands, to function as low-slope terraces to trap soil (as earlier hypothesized), or to protect young crops from rapid surface runoff into the depressions, near the end of the rainy season (for more on potential cultivation methods, see Fedick et al. 2000:145). These particular wetlands generally did not flood in the Late Preclassic because the water table was perhaps ~1m lower than it is today (Fedick et al., 2000:143). Swampy, but unflooded, conditions prevailed as evidenced by the lower units observed in the cores.
Subsequently, sometime in the Early Classic, wetlands began to flood more regularly, with algae communities developing and starting to build up calcareous muds. Gradually the number of wetlands suitable for dry season farming dwindled and impacted the regional subsistence system, and gradually the region was abandoned sometime around AD 400. A similar explanation is given for the Late Preclassic period abandonment of wetland farming at Albion Island in Belize (Pohl et al., 1990).
When Maya people returned to the Yalahau region in the Late Postclassic (AD 1200-1542), there is good evidence from the ancient settlement of T’isil (Fedick et al., 2012) that people took advantage of the more regularly flooded wetland and associated resources. For example, at T’isil (Figure 1), fish bone, turtle shell, and apple snails were recovered from Postclassic deposits, as were fish-net weights made on Postclassic ceramic sherds (Goossens, 2005).
In near-coastal areas around the Yucatan peninsula, numerous studies have documented late Holocene sea rise (Fairbanks, 1989; Lambeck and Chapell, 2001; Toscano and Macintyre, 2003) and consequent water table rise (Perry et al., 2003). At wetlands in northern Belize, soil-sediment sequences record this water table rise, as well as Maya responses such as canal-digging and raised field construction (Turner and Harrison, 1983; Alcala-Herrera et al., 1994; Pope et al., 1996; Beach et al., 2009). Elsewhere around the Peninsula, sediment cores extracted from lakes have been interpreted to reflect climate changes (Leyden et al., 1996; Brenner et al., 2003; Hodell et al., 2007). For example, at Punta Laguna in the northeastern Yucatan Peninsula (just south of the Yalahau region), alternating bands of organic- and carbonate-enriched strata were interpreted to reflect lower and higher lake stands, suggested to be a product of lesser and greater rainfall (Hodell et al., 2007). Periods of lesser rainfall (droughts) in turn have been linked to Maya sociopolitical reorganization during the Late Preclassic (Carrillo-Bastos et al., 2010; Kennett et al., 2012; Douglas et al., 2015) as well as to the Terminal Classic collapse (Hodell et al., 1995, 2001; Haug et al., 2003; Hodell et al., 2005, 2007; Webster et al.,2007; Medina-Elizalde et al., 2010).
In the Yalahau region, the soil sequences considered in this study provide a record of changing hydrologic conditions and, potentially, a record of anthropogenic manipulation of wetland soilscapes. According to previous soil studies in the region (Solleiro-Rebolledo et al., 2011; Wollwage et al., 2012), wetlands were impacted both by a gradual rise in water table over the last several thousand years, as well as changes in climate including droughts and moist periods. Water table rise may have occurred in starts and stops in accordance with sea rise, as documented by Beddows et al. (2016), independent of rainfall patterns such as droughts and moist periods such as those recorded by Hodell et al. (1995, 2001, 2005, 2007) and many others. Nevertheless, the overall trend has been one of rising sea level, increased flooding and horizontal expansion of wetland areas, and aggradation of wetland soilscapes. In the profiles studied, this is recorded by the transition from lower units to upper units. Unique soil features identified in the lower units of the cores reflect swampy but generally unflooded conditions. Subsequently, increases in flooding and water table level lead to the formation of strata of primarily micritic groundmass - the upper units.
It is in the context of these changing hydrologic conditions and shifts from swampy to more regularly flooded conditions that we must consider the archaeology: specifically, the fact that the Maya constructed hundreds of rock alignments in 20 of 26 wetlands thus far surveyed. It is important to consider that the transition from lower to upper units occurred at different times in different wetlands, probably dependent on individual wetland elevation and surrounding topography. This is evident from the radiocarbon dates. For example, in the Esperanza wetland, the transition from swampy conditions to regular annual flooding appears to have occurred very early on (prior to 1460 to 1292 BC), long before we have evidence of settlement in the region. This does not mean that it was not used by the Maya as a resource of the subsistence economy. However, the paucity of rock alignments (only 3) suggests it was not suitable for manipulation with rock alignments, or that rock alignments were not needed. In Zanja SP-1 the transition occurred at the end of the Late Preclassic period (between 205-37 BC and 323-407 AD), and in Media Luna SP-1 it occurred early in the Classic period (after 248-413 AD).
5. Conclusions
The research presented in this study helps to visualize wetland landscapes of the Yalahau region and how these landscapes changed over three millennia of Maya occupation. The findings demonstrate that the wetland ecosystems of today; soils, biota, and hydrology, are very different than those encountered by the earliest Maya settlers. The findings, while preliminary in nature, also suggest that a significant regional change in wetland flooding regimes occurred, with lower soil units representing a lower water table and variable degrees of saturation or flooding, and upper soil units representing a higher water table and more uniform flooding.
The micromorphological observations of the study cores of the Yalahau region were used to identify the various types of pedogenetic features and to distinguish two units: an upper unit dominated by biogenic micrite, and a lower unit characterized by a variety of unique soil features such as sparitic calcite crystals, pyrite, gypsum, and gleyic properties. Based on these observations, it is suggested that modern soil forming processes are fairly uniform across the region and controlled by regular annual flooding (~7 months of >1 m deep flooding each year), while in the past (Preclassic period), flooding was less regular and wetland environments were more swampy and varied.
It is suggested that the Preclassic depressions encountered by the Maya were a prime agricultural resource that allowed the region to flourish during the Late Preclassic period when the regional population level peaked (Fedick and Mathews, 2005; Rissolo et al.,2005; Glover and Stanton, 2010; Glover, 2012). During this time, as shown by the radiocarbon dated soil sequences studied here, it is evident that at least some wetlands such as Media Luna, were swampy but unflooded. Media Luna and Pichol (with a similar profile to Media Luna) are wetlands with some of the highest concentrations of rock alignments. It is suggested that the concentration of alignments at Media Luna and Pichol relates to a later onset of flooding than in wetlands such as Esperanza, where flooding began in the Early Preclassic. In Media Luna and Pichol, cores were extracted in close proximity to rock alignments, and it is further suspected that alignments are stratigraphically associated with lower unit soil cover when at least the margins of the wetland rarely flooded. In other wetlands that were surveyed and cored, but not studied here, rock alignments are common at wetlands with lower unit organic strata (similar to Media Luna and Pichol).
Future soils research should continue investigating the diversity of soil resources that exist throughout the region, specifically focusing on the development of well-dated paleoenvironmental reconstructions of changing wetland flood conditions in many more wetlands. In particular, additional dates are needed from Media Luna to better anchor the chronological framework of that site. Future investigations should also examine how variable hydrological histories of individual wetlands may be influenced by localized topography and drainage patterns, in addition to regional water table fluctuations resulting from sea level change. It may also be possible to better define the distribution of rock alignments in relation to wetlands with varying hydrological histories, and to identify and date the stratigraphic position of rock alignments.