Next Article in Journal
The Asian Summer Monsoon: Teleconnections and Forcing Mechanisms—A Review from Chinese Speleothem δ18O Records
Previous Article in Journal
The Geology of the Middle Cam Valley, Cambridgeshire, UK
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Apparent Resilience of the Dry Tropical Forests of the Nicaraguan Region of the Central American Dry Corridor to Variations in Climate Over the Last C. 1200 Years

1
Long-Term Ecology Laboratory, Department of Zoology, University of Oxford, Oxford OX1 3SZ, UK
2
Department of Geology and Environmental Geosciences, Northern Illinois University, DeKalb, IL 60115, USA
3
School of Geography and Environmental Sciences, University of Southampton, Highfield, Southampton SO17 1BJ, UK
*
Author to whom correspondence should be addressed.
Submission received: 29 May 2019 / Revised: 7 July 2019 / Accepted: 11 July 2019 / Published: 18 July 2019

Abstract

:
The Central American Dry Corridor (CADC) is the most densely populated area of the Central American Isthmus and is subject to the greatest variability in precipitation between seasons. The vegetation of this region is composed of Dry Tropical Forests (DTF), which are suggested to be highly susceptible to variations in climate and anthropogenic development. This study examines the vulnerability of past DTF surrounding the Asese peninsula, Nicaragua to climatic and anthropogenic disturbances over the past c. 1200 years. Past vegetation, climate, burning, and animal abundance were reconstructed using proxy analysis of fossil pollen, diatoms, macroscopic charcoal, and Sporormiella. Results from this research suggest that DTF have been highly resilient to past climatic and anthropogenic perturbations. Changes in DTF structure and composition appear to be linked to the abundance and intensity of fire. Pre-Columbian anthropogenic impacts on DTF are not detected in the record; however, DTF taxa decline slightly after European contact (1522 C.E.). Overall the DTF for the Nicaraguan region of the CADC were found to be highly resilient to both climatic and anthropogenic disturbances, suggesting that this region will continue to be resilient in the face of future population expansion and climatic variation.

1. Introduction

Significant losses in food production across the Central American Isthmus (Isthmus of Tehuantepec south to the Isthmus of Darien) resulting from a deficit in precipitation at the beginning of the harvest in 2015 rendered an estimated 2.2 million people at risk of moderate or severe food insecurity [1]. This episode is indicative of the vulnerability of the flora, fauna and human population of the Central American Dry Corridor (CADC) to seasonal patterns of rainfall, impacting agrarian practices, food and water security [2,3]. Given this recent episode of apparent vulnerability of the CADC to hydroclimatic changes, the aim of this research was to identify past hydroclimatic changes and assess their impacts upon the dominant vegetation type, Dry Tropical Forests (DTF), for the Central Pacific lowlands of Nicaragua spanning the last c. 1200 years (Figure 1).

1.1. Climate

Previous research examining the spatial patterns and mechanisms driving precipitation change across the Central American Isthmus over the past 2000 years have focused upon the impacts and timings of hydrological extremes including the Medieval Climate Anomaly (MCA) 950–1250 C.E. and Little Ice Age (LIA) 1400–1850 C.E. [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20] and the transition between the two [21]. These changes in hydroclimate are reported to be primarily driven by the interactions of three climatic forces: (i) the position of the Inter Tropical Convergence Zone (ITCZ); (ii) the North Atlantic Oscillation; and (iii) the El Niño Southern Oscillation (ENSO) [22]. Warmer sea surface temperatures (SST) in the Caribbean and Gulf of Mexico are inferred for the early part of the MCA [23], and these warmer SST are associated with a more northerly position of the ITCZ as inferred from the Cariaco Basin titanium record [24]. Cooler SST in the Caribbean and Gulf of Mexico were associated with a southward shift in the ITCZ, and a weakening of the summer monsoon during the LIA [5].
Hydroclimate records from Central America spanning the MCA and LIA are reported to be spatially and temporally heterogeneous, with many records lacking sufficient chronological control or resolution to examine this period in detail [22]. Previous work reconstructing past precipitation across the CADC during this time period has been limited with only a single terrestrial δ18O record published from within this region spanning the past 1400 years [16]. This δ18O record indicates comparatively wetter conditions during the MCA and drier conditions during the LIA.

1.2. Vegetation

A range of vegetation types are found in the CADC which can be broadly categorised as: (i) dry tropical forests (DTF); (ii) moist tropical forests (MTF); (iii) coniferous tropical forests (CTF); and (iv) mangroves [25]. These formations are primarily distributed according to climate (temperature and precipitation), soil type (e.g., volcanic, alluvial), and elevation (highlands and lowlands). Past research reconstructing vegetation dynamics within the CADC has been somewhat limited, particularly for DTF and for the last 1200 years. Past anthropogenic influences are typically inferred from agricultural taxa such as Zea mays, anthropogenic burning, and evidence from nearby archaeological sites. In comparison, greater forest cover, reduced burning, and an absence of agricultural taxa, appear to coincide with intervals of evidence for reduced population [26]. In contrast, intervals with known archaeological evidence for more people are associated with greater forest clearance (using fire), an abundance of agricultural crops and the increase in abundance of economic or useful taxa [27,28,29].
In total three sites have been examined from the DTF to date indicating that there have been vegetation responses to known forcing over this interval in time. All three of these studies by Dull [27,28,29], looked at the vegetation structure and composition of the Ahuachapan region in the CADC of El Salvador, and found that anthropogenic burning between c. 500 B.C.E.–1500 C.E. resulted in a more open grassland environment at the expense of woody taxa. Following this disturbance, Urticales and Mimosoideae were found to have established and were attributed to the population collapse of the region. It is apparent that decreases in arboreal cover are commonly associated with anthropogenic activity, or drier conditions, and more fire; and vice versa for wetter conditions [27].
The resilience of the DTF to climatic and anthropogenic impacts is currently poorly understood [30]. DTF once covered more than 40% of the total area of all tropical forests [31]. They are one of the most threatened habitats globally and are argued to deserve a high priority for conservation [32,33,34]. It is suggested that either the DTF (i) will be sensitive to climatic change because they are already limited by water and are close to a habitat threshold; or (ii) they will be resilient to changes in rainfall because they are adapted to predictable seasonal drought. Similarly, DTF will either be (i) vulnerable to anthropogenic activities, such as agrarian practices, land conversion and population density because of land use change; or (ii) they will be resilient to anthropogenic factors due to composition adaptation over time [35]. Disturbance, change, and evolution through time can lead to the reorganisation and renewal of the ecosystem [36]. We define resilience by the capacity of a system (biotic or abiotic) to be perturbed while remaining in the same state.
In order to fill these knowledge gaps (detect the impacts of climate and people upon the biota residing within the CADC), this study reconstructed past hydroclimate, biomass burning, DTF dynamics and herbivorous mammal presence, over the last 1200 years (c. 840–2004 C.E). These results were then examined within the context of the archaeological record from the Central Pacific Nicaraguan region. The aim of this study was to (i) identify the past hydroclimatic changes within the Nicaraguan region of CADC; (ii) assess how the DTF responded to identified climatic shifts; and (iii) assess how anthropogenic activities have impacted the DTF of Central Pacific Nicaragua.

1.3. Study Area

The CADC is subject to severe and prolonged periods of drought. The dominant type of vegetation in the lowland (<1000 m) areas are DTF [37]. These occur from Chiapas (Mexico) and extend south across the central and lowland areas of Pacific Guatemala through El Salvador, Honduras, Nicaragua, and Costa Rica [38]. DTF occur: (i) in frost-free environments with mean annual temperatures greater than 17 °C; (ii) mean annual rainfall of 250–2000mm; and (iii) overall greater potential of evapotranspiration to precipitation [32]. DTF are subject to a long dry season (5–8 months) and a shorter wet season [2].
The most common taxa found to comprise DTF include: Bursera simaruba; Caesalpinia coraria; Calycophyllum candidissimum; Guazuma ulmifolia; Gyrocarpus americanus; Haematoxylon brasiletto; Luehea candida; Lysiloma kellermani; and Phyllostylon brasiliensis [39].
These formations are typically smaller in stature, lower in biomass and less complex floristically and structurally than tropical rain forests [31]. The canopy is usually closed with no substantial opening in the crown (<20%) reaching a height of approximately 15–20 m [40]. Most plants in the DTF cease their vegetative growth in the dry season and shed their leaves, while other produce fruits and spread their seed [40].
Nicaragua is currently covered by 2500 km2 of DTF, which represents around 2% of the total forest cover [41]. Prior to European contact, savannas were created within the dry tropical forests by inhabiting human populations, especially along the Pacific coast [42,43]. The DTF were more intensively anthropogenically managed compared to the MTF due to the abundance of more productive soils for crop rotation, also known as milpa [31,44]); however, management practices were inefficient in clearing land, allowing for immediate regeneration after several years of cultivation [44].
The Pre-Columbian population inhabiting Nicaragua (before 1522 C.E.) is unknown; however, several attempts to estimate total population have been made, with upper estimates of c. 800,000 [45] to 1,000,000 [46,47,48] and lower estimates of 100,000 [49]. With European contact and the Spanish conquest of Nicaragua, indigenous populations were reported to have been significantly reduced from 500,000 to 200,000 by 1535 C.E. [50]. Latter conquests and occupation of the Meseta Central of Costa Rica (in 1570 C.E.) reported further depopulation of up to 90% down to c. 61,000 though warfare, slavery and introduction of epidemic diseases [45,51,52]. By the end of the 16th Century there were only around 500 Spanish colonists living in Nicaragua [53]. Populations are suggested to have only started to increase again by c. 1800 C.E (c. 83,000) [51] and reach upper estimates of Pre-Columbian contact levels (c. 1,000,000) by c. 1950 [54]. The current population is around six times (c. 6,272,133) that of the maximum estimates for pre-contact levels [55].
Small-scale deforestation began in 1522 C.E. with the arrival of the Spanish [42,56,57,58], and continues through present-day. Pasture conversion for cattle was the main cause of deforestation in Nicaragua, which is thought to have impacted DTF more extensively than MTF [32,59,60]. After the Central American Act of Independence in 1821 C.E. landless peasants began to clear land in the DTF [58,60]. DTF were popular settlement areas because land was easier to clear, soil fertility was higher, diseases were rare, and high-quality timbers could be found within these forests [31,59]. In the early 1900s timber extraction, primarily from the DTF, increased for export [57,61]. Targeted species included: Swietenia humilis, Cedrela odorata, Pachira quinata, Dalbergia retusa, and Guaiacum sanctum [61]. The major causes of DTF loss in recent years are due to agricultural expansion for coffee plantations, crop fields and further ranches [62]. Only fragments of DTF remain today [41].

1.4. El Gancho

El Gancho (11.906°N, 85.918°W, 44 m above sea level) is a small, shallow (1.1 m), closed basin lake, situated on top of the Asese peninsula. The Asese peninsula protrudes into the northwest quadrant of Cocibolca (Lake Nicaragua) next to the city of Granada founded in 1524 C.E. (Figure 1). El Gancho is thought to have been formed sometime after 140–345 C.E. when the Asese peninsula was created as the result of a c. 57 km2 debris avalanche originating from the northeast flank of the Mombacho Volcano [16,63,64]. Sedimentological analysis suggests that sediments have continuously accumulated in El Gancho since its formation [16,64]. Total annual rainfall is 1298 mm and mostly falls between May–October [65]. Observations in the field between 2014–2017 C.E. capture variations in lake level derived from precipitation and evaporation: (i) before, (ii) during, and (iii) after a two-year ENSO event (Figure 2).
The vegetation assemblage currently on the Asese peninsula is comprised of taxa belonging to the DTF [25,66,67]. The vegetation is comprised of a higher canopy, primarily deciduous trees: e.g., Bursera simaruba, Gliricidia sepium, Diospyros acapulcensis, Cochlospermum vitifolium, and Myriocarpa bifurca; and an understory which is typically comprised of sparse and scattered evergreen or deciduous shrubs, and grass tufts, e.g., Amaranthus spinosus, Compositae, and Poaceae [33,39,66,67,68,69]. The understory is typically burnt by landowners to cleanse unwanted weeds and ticks and to prepare the land for agriculture or ranching [69].
The area surrounding the Asese peninsula, was once DTF but now predominantly agricultural comprising cattle ranching and the growing of corn (Zea mays), beans (Leguminosae) and rice (Oryza). Most of this agriculture is grown for export [70].
There are few wild animals left in this area due to extensive land management (e.g., cattle ranching and tourism); however, the mantled howler (Alouatta palliata), domesticated dogs (Canis), various species of iguana (Iguana) and rabbit (Sylvilagus) can be found on the Asese peninsula and surrounding islands (Las Isletas). The Asese peninsula is sparsely populated with only one or two dwellings; however, the city of Granada (less than 10 km away), is home to c. 124,000 people [71].

2. Methods

2.1. Field and Sampling Techniques

In June 2004 a bathymetric survey of El Gancho was conducted using a handheld sonar to identify the deepest part of the lake (1.1 m) from which a 277 cm composite core, with overlapping sections, was retrieved using a Livingstone piston corer [72]. The upper ~168 cm of the core contains massive, very dark brown to black, organic-rich sediments with no notable changes in the stratigraphy [16]. The top 34 cm of the core were extruded in the field at 0.5 cm intervals to ensure sediments did not mix. Forty-two subsamples (1 g wet weight) were extracted at approximately 4 cm intervals down through the upper 168 cm of this sequence for biological proxy analysis: macro-charcoal, pollen, and dung fungal spores (Sporormiella); and twenty-two subsamples (0.1 g wet weight) for diatom analysis at an interval of 8 cm (3.5–14.5 range). The biological proxies in the lower 109 cm of the core were found to be degraded and were therefore not analysed further.

2.2. Chronology

An age-depth model was constructed using five calibrated radiocarbon dates conducted on charcoal fragments published in Stansell et al. [16]. The radiocarbon dates were (i) recalibrated against the IntCal13 radiocarbon dataset [73]; and (ii) converted to calendar ages and modelled using OxCal 4.3 (Table 1). The Bayesian model applied used a Poisson distribution and made no assumptions of the priors. The Poisson parameter k0 was taken as 1 depth−1, the interpolation rate as 0.7 depth−1 and the log10(k/k0) as U(−2,2). The sedimentation rate per year was calculated from the results of the age-depth model and used in the calculation of pollen concentrations and charcoal influx.

2.3. Diatom Analysis

Fossil diatoms were used to infer past changes in lake level and area. They were cleaned and prepared using standard procedures as outlined in Berglund [74] and followed the Oxford Long-Term Ecology Laboratory diatom protocol [75]. Samples were mounted using Naphrax resin, chosen for its high refractive index, 1.73 [76]. Counting and identification of diatoms were conducted at 400x and 1000x magnification. A total of 300 frustules per level were counted and morphologically identified (see supplementary materials Table S1) using published keys [77,78,79,80,81,82,83]. Annual diatom concentrations were calculated using the evaporative tray method and results from the age-depth model [74]. The diatom samples were matched to the sampling resolution of the palynological data through linear interpolation for ease of comparative analysis.

2.4. Fossil Pollen and Dung-Spore Analysis

Fossil pollen was used to reconstruct the composition of past DTF vegetation. Fossil pollen extraction and preparation followed standard palynological procedures applying the Oxford Long-Term Ecology Laboratory protocol [84]. Silicone oil was used as the mounting agent to allow for rotation of grains, easing identification. Samples were spiked with known concentrations of Lycopodium spores (batch No. 212761 or 3862) to calculate pollen accumulation rates [85]. Counting and identification of pollen grains were conducted at 400x and 1000x magnification. For each level, a minimum of 300 terrestrial pollen grains were counted. Morphological identification (see supplementary materials Table S2) was achieved using (i) pollen databases [86,87,88]; (ii) published plates: [89,90]; and (iii) botanical reference materials from the OxLEL reference collections and specimens collected in the field (see supplementary materials Figure S1). Amaranthus spinosus was removed from the total palynological counts due to overrepresentation within the assemblage [91]. Total abundance of Amaranthus spinosus averaged to 73% of total pollen abundance with a maximum abundance of 90%. Overrepresentation has been attributed to growth within the catchment area and also because species in the Amaranthaceae family are known for their high production of pollen and seed [92]. The abundance of coprophilous dung fungal spores (Sporormiella) were used to indicate herbivorous mammal abundance. Sporormiella spores were counted and morphologically identified (from their sigmoid aperture) on the same slides [93,94].

2.5. Macroscopic Charcoal Analysis

Macroscopic fossil charcoal (>150 µm) was used as a proxy to indicate past occurrence of local fires, representing burning within a 10 km radius of the lake basin [95,96,97,98,99,100]. Fragments were separated from the sample using a sieve [101]. All charcoal over 150 µm was counted at 10x magnification.

2.6. Data Handling

Pollen and diatom counts were converted to percentages while Sporormiella and macroscopic charcoal were converted to influx using the exotic marker (Lycopodium) and the sedimentation rate [95,102,103,104,105]. To identify discrete zones in the palynological assemblage, constrained hierarchical clustering upon the palynological data set was conducted following the broken stick model [106]. Before performing all ordinations analyses a square-root transform was applied to the percentage data to normalise the distribution. This method of transformation was chosen because it can be applied to data sets containing zero values [107]. To check which ordination method was most appropriate for the palynological and diatom data sets, detrended correspondence analysis was conducted and the number of standard deviations away from the mean was assessed [107]. Principal component analysis (PCA) was selected and used to infer similarities between samples and taxa applying a singular value decomposition of the centred but not scaled data matrix. Canonical correspondence analysis (CCA) conducted to quantify the relationship between environmental variables and the palynological and diatom assemblage data. δ18O data from Stansell et al. [16] was matched through interpolation to the resolution of the palynological and diatom data sets for the CCA. Ellipses representing the discreet Zones were calculated using the parameterization (cos(theta + d/2), cos(theta − d/2)), where cos(d) is the correlation of the parameters was applied [108]. These were conducted at a confidence level of 95%. These analyses and the presentation of data was performed using R [109], applying packages Vegan [110] and Rioja [111].

3. Results

3.1. Chronology and Resolution

Re-calibration and modelling of the five radiocarbon dates published by Stansell et al. [16] indicated that the recovered composite core encompasses the past ~1700 years in continuous sedimentation (Figure 3). Distance between subsamples for palynological analysis represent an average resolution of c. 30 years. The age-depth model presents good overall agreement (index = 102.9) with an average age range of 190 years surrounding the median modelled age. The model agreement index varied between 99.1 and 110.6 (Table 1).

3.2. Palaeoecological Trends: Diatoms

The diatom assemblage is well represented by a dynamic assemblage of facultatively planktonic, planktonic, benthic and epiphytic taxa (Figure 4). Changes within the diatom assemblage broadly reflect the palynologically defined Zones which have been superimposed.
Zone 1 (168–117 cm, 5 samples, 850–1180 C.E.). Benthic (32.2–60.7%) and ephyphitic (18.5–2.3%) taxa dominate this first zone (Mastogloia smithii, Nitzschia frustulum, Navicula longicephala, Nitzschia amphibia, Encyonema silesiacum and Anomoeoneis vitrea), but decrease in abundance after 1030 C.E. (90.8–48.1%). Cyclotella meneghiniana rises from 1030–1150 C.E. (6.1–20.3%) as does Pseudostaurosira brevistriata (2.4–23.4%). Total frustule abundance is lowest in this zone (6.6–2.5 x106 frustules 1 g3 dry sediment).
Zone 2 (107.5–58 cm, 10 samples, 1180–1550 C.E.). Benthic and ephyphitic taxa diminish. Pseudostaurosira brevistriata becomes the dominant taxa between 1150–1230 C.E. (23.4–59.9%). After 1340 C.E. Aulacoseira ambigua (47.3%) succeeds Pseudostaurosira brevistriata (27.4%) in abundance, dominating the assemblage until the end of the zone (41.2%). Cyclotella meneghiniana rises between 1450–1530 C.E. (12.6–32.3%). Frustule counts gradually increase (2.4–4 x106 frustules 1 g3 dry sediment).
Zone 3 (55.5–3.5 cm, 7 samples, 1550–2004 C.E.). From 1600 C.E. Pseudostaurosira brevistriata (45.8%) sucedes Aulacoseira ambigua (29%). All other taxa diminish except for Nitzschia frustulum which increases from 1730–1970 C.E. (3.1–11.8%). Cyclotella meneghiniana remains in high abundance until 1600 C.E. (12.1%). Diatom frustules increase from 1570–1600 C.E. (8,863,231–45,552,185 frustules 1 g3 dry sediment) and 1810–1930 C.E. (37,353,262–95,983,203 frustules 1 g3 dry sediment), and decrease in number from 1700–1810 C.E. (63,663,242–37,353,262 frustules 1 g3 dry sediment).
Ordination of the diatom assemblage data denotes clear taxonomic associations with each of the superimposed palynological Zones (Figure 5). The PCA explains 40.2% of the assemblage variance along Axis 1 and 23.4% along Axis 2 (Figure 5a). Zone 1 is predominantly associated with benthic and ephyphitic taxa (e.g., Mastogloia smithii and Nitzschia frustulum, Navicula longicephala, Nitzschia amphibia, Encyonema silesiacum and Anomoeoneis vitrea); Zone 2 is most associated with Aulacoseira and Navicula taxa; while Zone 3 defined by Pseudostaurosira brevistriata and Pinnularia viridis. Results from the CCA show a clear significant association with δ18O along the first axis. More positive δ18O are most associated with Zones 2 and 3 (Figure 5b). The δ18O record presented in Stansell et al. (2013) was interpreted to represent drier conditions with more positive δ18O values and wetter conditions with more negative δ18O values.

3.3. Palaeoecological Trends: Pollen, Sporormiella and Macroscopic Charcoal

Zonation of the palynological sequence indicates that three numerically derived zones can be recognised between 840–2004 C.E. (Figure 6 and Figure 7). DTF taxa (e.g., Bursera simaruba, Brosimum, Protium and Terminalia) are well represented throughout the entire 1200 years analysed, comprising 48–95% of the total palynological assemblage; however, the DTF community composition varies through time (Figure 6 and Figure 7).
Zone 1 (168–117 cm, 9 samples, 850–1180 C.E.). Bursera simaruba (27.1–34.5%), is abundant in the first zone of this sequence, together with Leguminosae (11.5–15.5%), Protium (10–12.1%) and Mimosa (2.9–8.6%) from the DTF taxa. Taxa comprising the herbaceous understory (Poaceae, Compositeae, Cyperaceae and Spananthe paniculata) are abundant (33.3–17.3%) particularly at 890 C.E. (36.7%), 1000 C.E. (39.3%) and between 1090–1120 C.E. (36.1–38.2%). Non-DTF arboreal taxa are represented in low abundance by Pinus (1.4–4.4%), Alchornea (0–1.4%) and Roupala montana (0–0.9%). Sporormiella influx is high at 980 C.E. (2.6 spores cm−2yr−1) between 1030–1090 C.E. (4.3–3.3 spores cm−2yr−1) and at 1150 C.E. (4.5 spores cm−2yr−1). Macroscopic charcoal influx is high (10.3–38.1 particles cm−2yr−1), peaking between 920–980 C.E. (33.2–29.8 particles cm−2yr−1) and at 1090 C.E. (38.1 particles cm−2yr−1). Pollen influx is low (45,676–109,410 grains cm−2yr−1).
Zone 2 (107.5–58 cm, 18 samples, 1180–1550 C.E.). DTF taxa increase from 1180–1370 C.E. (75–95.3%) and persistently dominate the assemblage in abundances greater than 70% through to the end of this zone. High abundance DTF taxa include Bursera simaruba (6.5–45.6%), Brosimum (0–38.7%), Terminalia (4.1–16.4%), Protium (0–13%) and Celtis (0–16.7%). Celtis peaks at 1230 C.E. (16.7%). Coussapoa is persistent in low abundances between 1230–1530C.E (2.4–4.4%), at which point Anacardiaceae (3.3%) and Mimosa (3.3%) supersede. Herbaceous taxa decline and remain in low abundance until after 1440 C.E. (20.8–3.6%). From 1480 C.E., herbaceous taxa begin to increase through to the end of this zone (5.3–24.4%). Capparis (Caper) appears for the first time at c. 1500 C.E. accompanied by an increase in Poaceae. (2.6–6.7%). The abundance of non-DTF arboreal taxa remain in low abundance with the exception of a short-lived peak from Roupala montana at 1320 C.E. (23.7%). Sporormiella influx decreases from the beginning of this Zone (2.7 spores cm−2yr−1) and remains in low abundance between 1230–1480 C.E (0–0.7 spores cm−2yr−1). After 1480 C.E. Sporormiella influx increases substantially (0.7–2.7 spores cm−2yr−1). Macroscopic charcoal influx remains high until 1300 C.E. (30.6 particles cm−2yr−1) where it transitions down to a lower influx by 1400 C.E. (5.1 particles cm−2yr−1). Between 1400–1530 C.E. macroscopic charcoal influx remains comparatively low (0.1–7.9 particles cm−2yr−1). Pollen influx increases particularly after 1300 C.E. (266,775 grains cm−2yr−1) and is particularly high between 1440–1480 C.E. (1,413,110–1,333,889 grains cm−2yr−1).
Zone 3 (55.5–3.5 cm, 15 samples, 1550–2004 C.E.). DTF taxa decrease from 1560–1640 C.E. (83.6–48.4%) and increase from 1640–1750 C.E. (48.4–80%). The DTF assemblage is predominantly comprised of Bursera simaruba (6.4–36.7%), Brosimum (4.4–33.3%), Protium (0–9.6%), Leguminosae (2.2–7.1%), Anacardiaceae (1.4–10.1%), and Rhus (0–5.9%). Poaceae and Cyperaceae increase in abundance permanently after 1530 C.E. and 1600 C.E. Cyperaceae peaks at 1600 C.E. (15.8%) and 1640 C.E. (25.8%). DTF taxa decrease in abundance between 1750–1850C.E. reflected in the reduction of Bursera simaruba (36.7–19.5%), Brosimum (10–2.4%), and Mimosa (13.3–0%). Herbaceous taxa increase from 1750–1850 C.E. (16.7–46.3%). DTF abundance increases again from after 1850 C.E. (48.8%) through to 1970 C.E. (79.1%) with the increased abundance of Bursera simaruba (19.5–29.8%), Protium (4.9–10.5%), Leguminosae (2.4–11.9%), Rhus (2.4–10.5%), and Terminalia (0–6%). Sporormiella influx remains high between 1550–1570 C.E. (2–2.8 spores cm−2yr−1) and is high in 1650 C.E. (3.4 spores cm−2yr−1). After 1650 C.E. Sporormiella influx declines (3.4–0.8 spores cm−2yr−1). The abundance of macroscopic charcoal decreases towards the present from 1560–1970 C.E (3–0.03 particles cm−2yr−1). Pollen influx reduces overall (682,686–145,700 grains cm−2yr−1) but peaks at 1650 C.E. and (1,773,492 grains cm−2yr−1) and 1850 C.E. (1,391,104 grains cm−2yr−1).
PCA of the pollen assemblage displays taxa most associated with each Zone and explains 16.4% of variance on Axis 1 and 12.6% of variance on Axis 2 (Figure 7a). Zone 1 is most associated with Protium, Leguminosae, Compositeae and Spananthe paniculate; Zone 2 with Bursera simaruba, Brosimum, Terminalia, Coussapoa and Celtis; and Zone 3 with Rhus, Trema, Cordia, Capparis, Curtella americana and Mimosa.
CCA of the palynological assemblage delineates that macroscopic charcoal, Sporormiella and δ18O are all statistically significant drivers of variation within the palynological assemblage. Macroscopic charcoal is most associated with Zones 1 and 2 along the first axis, orthogonally placed to δ18O which is most associated with Zone 3 (Figure 7b). Sporormiella is associated in part with all Zones along the second axis.

4. Discussion

4.1. Hydroclimate Changes in the Nicaraguan Region of the Central American Dry Corridor during the Medieval Climate Anomaly and Little Ice Age

Interpretation of the diatom assemblage indicates that the Nicaraguan region of the CADC was wettest during the MCA (950–1250 C.E.) and drier during the transition to the LIA (1250–1400 C.E.) and throughout the LIA (1400–1800 C.E.). The transition between the MCA and LIA surrounding El Gancho is marked by a shift from wetter to drier conditions around 1300 C.E.
The high abundance of benthic and epiphytic diatom taxa between 840–1180 C.E. (Zone 1) indicates that wetter conditions were apparent in the region in the earliest part of this record during the MCA. The shallow basin that El Gancho currently occupies is poorly constrained and is surrounded by low lying topography (Figure 2). Subject to wetter conditions, El Gancho’s basin is suggested to have expanded laterally, increasing in area but not necessarily depth. El Gancho is assumed to have been at its maximum expanse and depth during this period time. The high abundance of benthic and epiphytic diatom taxa, as well as Cyperaceae from the palynological assemblage, therefore, probably reflect the permanent inundation of water for a larger area by occupying the available habitat niche (e.g., Cyperaceae and Potamogeton) [27,28,29].
The time surrounding the transition between the MCA and LIA (c. 1200–1400 C.E.) leading into the first half of the LIA (1400–1600 C.E.) is defined by the replacement of benthic and epiphytic diatoms by the increasing abundance of diatom taxa indicative of changing conditions (e.g., Pseudostaurosira brevistriata and Aulacoseira ambigua), inferred to indicate a trend towards drier conditions (Figure 4 and Figure 5). While the entire lake habitat is suitable for benthic and epiphytic diatoms to thrive, we suggest that these taxa were outcompeted by Pseudostaurosira brevistriata and Aulacoseira spp.
The rapid expansion and retraction of this water body between the wet season and the dry season is thought to have facilitated continuous biannual disturbance promoting a habitat more conducive for these two taxa to dominate the assemblage. This evidence further suggests that El Gancho retreated in area abandoning the large shallow expanse surrounding the central depocenter. After c. 1200 C.E. diatoms Aulacoseira spp. and Pseudostaurosira brevistriata comprise between 50–90% of the diatom assemblage. The succession of Aulacoseira ambigua at 1340–1600 C.E. is suggested to represent a rapid reduction in lake area, whereby the water occupying El Gancho retreated to the deepest part of its basin, abandoning the expansive shallow periphery. The succession back to facultatively planktonic taxa (Pseudostaurosira brevistriata) after 1600 C.E. is inferred to indicate that El Gancho continued to reduce in size, losing depth and available habitat for the planktonic diatom taxa to thrive (Figure 4). Pseudostaurosira brevistriata is a facultatively planktonic taxa and therefore we suggest that the available habitat niche provided conditions amiable to this taxa’s success over the other taxa present, such as Aulacoseira ambigua, in effect diluting the abundance of the other taxa. We previously explored the salinity and nutrient groupings for all diatom taxa and did not find conclusive trends associated with the most likely scenario presented in our manuscript.
The interpretation of the diatom record supports the δ18O precipitation reconstruction presented by Stansell et al. [16] which also indicates wetter conditions for the MCA, particularly between c. 950–1300 C.E. and drier conditions thereafter. The response of the diatom assemblage to drier conditions during the transition between the MCA and LIA appears to lag behind the δ18O presented in Stansell et al. [16] by c. 50 years; this could represent a biological lag response to hydroclimatic changes (Figure 4).

4.2. How Did the Dry Tropical Forests Respond to These Wet–Dry Climatic Shifts?

Three successional stages can be discerned from the last 1200 years of vegetation change in the DTF occupying the Asese peninsula (Figure 6). The first zone corresponds with wetter conditions from c. 800–1200 C.E., while the second (c. 1200–1550 C.E.) and third Zones (c. 1550–2004 C.E.) with relatively drier conditions (Figure 6).
The wettest period of the last 1200 years (c. 800–1250 C.E.) is primarily comprised of DTF taxa. The dominance of these taxa in the pollen diagram are taken to indicate a relatively closed canopy. Perhaps unusually with this forest type, there is evidence for intervals of greater burning. We suggest that this is in response to the wetter conditions from c. 800–1250 C.E. facilitating greater accumulation of vegetative biomass. This would have enabled larger fuel loads to develop and burn between seasons [96,100]. Fire in the CADC generally occurs in the transitions between the dry season and wet season as biomass dries out [112]. There is also an increased abundance of Poaceae, Compositeae, Cyperaceae and Spananthe paniculate during the MCA, which can represent burning facilitated by some openings in the canopy [112,113,114].
At the transition into relatively drier conditions (1250–1400 C.E.) there is a change in the composition of DTF taxa including the succession of Brosimum, Terminalia and Celtis (Figure 6). The highest abundance of these DTF taxa occurs under these drier conditions. In conjunction with this change in composition, reduced burning during this time occurred, probably associated with lower fuel loads. This reduction in burning (c. 1300 C.E.) appears to have resulted in an increasingly closed canopy, blocking light to the understory and the reduction in the presence of open ground taxa. Evidence of this closed canopy is also presented in the increase in the abundance of diatom flora Aulacoseira spp. which is suggested to be responsive to low-light levels, thriving in a well-developed forest with closed canopies [115].

4.3. Anthropogenic Activity in the Dry Tropical Forests: Population and Land Management

Pre-Columbian human occupation of the Asese peninsula has been documented from at least c. 600 C.E. through to 1250 C.E. with a hiatus during the Late Sapoa period c. 1100–1150 C.E. [116]. The nearby occupation of two municipalities within the Masaya region, as well as the occupation of Ometepe, was continuous from c. 200 B.C.E. through to the present day [52,53,117].
Despite this continuous Pre-Columbian occupation, there is no evidence of substantial deforestation in the pollen record of the Asese peninsula during this interval in time (c. 850–1522 C.E.); however, evidence of forest clearing is available from other nearby paleo-records. Deforestation on Ometepe Island, for example, from 775 C.E indicates that areas of Nicaragua were undergoing deforestation during Pre-Columbian times [118].
While there is no palynological evidence of Pre-Columbian agriculture (e.g., Zea mays), the abundance of fruiting trees and edible herbs are likely to have been utilised by local people for materials and sustenance. Amaranthus spinosus, Anacardiaceae, Bursera simaruba, Brosimum, Capparis, Celtis, Leguminosae and Protium are all examples of edible flora found in El Gancho, and were also eaten as part of the Pre-Columbian diet [119]. Amaranthus spinosus is a pioneering herbaceous plant with a similar nutritional value to that of spinach, and dominates the palynological record throughout [120].
The introduction of Capparis in the palynological record from c. 1500 C.E. possibly provides evidence for European contact. Capparis (commonly known as a caper) places its origins surrounding the Mediterranean in Europe, this taxon is not native to the Americas and would only have arrived with the Spanish conquistadors [121]. In addition, the reduction of select arboreal taxa, e.g., Bursera simaruba, and Protium, could reflect deforestation for building materials and the creation of pastures; however, this is also within the range of natural variability [122,123].
Evidence of herbivorous mammalian presence from c. 950–1250 C.E. aligns with the inferred wetter conditions brought about during the MCA. There is no current evidence for Pre-Columbian domestication of the native fauna, and therefore we suggest that greater vegetative growth during this period of time provided an increased food source for deer, peccary, and tapir allowing for their numbers to swell. Reductions in dung fungal spore abundance after the shift towards drier conditions, between the MCA and LIA, suggests that herbivorous mammals were largely absent from the Asese peninsula after 1250 C.E. until the arrival of the Spanish. Dung fungal spore evidence indicates that animal abundance on the Asese peninsula increased after the arrival of the Spanish through until 1650 C.E. This possibly represents the import of cattle brought into Granada (established 1524 C.E.) from Europe. These cattle would possibly have been kept on the Asese peninsula, a natural pen, until they could be sold at market [124]. The rise in population and clearance of land for agrarian practices leading up to 1800 C.E. may well therefore account for the apparent reduction of DTF taxa in and around the Asese peninsula from c. 1700–1850 C.E. During this time the abundance of DTF taxa reduce and economically important taxa such as Anacardiaceae increase.

4.4. Resilience of the Dry Tropical Forests Surrounding the Asese Peninsula

The reconstruction of the DTF surrounding the Asese peninsula, over the past 1200 years, show little evidence of vulnerability to climatic and anthropogenic impacts. While the composition of taxa comprising the DTF assemblage changes in response to the identified hydroclimatic and anthropogenic perturbations, the ecological structure of the DTF has remained remarkably resilient. High seasonal variation within the CADC has exposed the residing DTF to biannual disturbance since the modern configuration of the North American Monsoon c. 4000 B.C.E. [21]. Cole et al. [125] identifies high frequency disturbance as the most significant factor in determining resilience in tropical forests, which is clear from the reconstructions of DTF around the Asese peninsula. Biannual disturbance caused by seasonal variation over the last four millennia is the likely reason why these DTF are so resilient, allowing these DTF to adapt and increase their resilience over time [126].
Studies of satellite data [127] and meta-analyses [35] examining resilience of vegetation in this region in recent decades to climatic perturbations suggest that DTF are sensitive to changes in precipitation/drought intensity, frequency, and/or timing; however, the results presented from El Gancho suggest that the DTF around the Asese peninsula are not particularly sensitive to climatic and anthropogenic fluctuations, at least not over the past 1200 years. While there is some evidence for the reorganisation of taxa between the palynological Zones, likely responding to more frequent or severe drought or variability in rainfall [128,129,130], the palynological abundance of DTF taxa in the record from El Gancho never fall below 50%, and exceeds 90% from c. 1350–1420 C.E., when fire reduces and climate switches to overall drier conditions (Figure 6). Further, these findings suggest that the vegetation surrounding the Asese peninsula has not been subject to a tipping point leading to a switch to a different stable state and vegetation type during the time frame studied [30].
Specific DTF taxa such as Bursera simaruba and Brosimum, which are prominent in the palynological record from El Gancho, are known to have a high wood density comprised of greater cell wall material and narrow vessels. These limit their hydraulic efficiency but increase their resistance to drought-induced water scarcity [131,132]. High wood density has previously been demonstrated to be an important factor in enabling persistence of woody taxa during intervals of drought [130]. These species may also resist extended periods of drought by tapping into subsoil water reserves with their deep roots [35,133,134].

5. Conclusions

The results from this research suggest that DTF have been highly resilient to past climatic and anthropogenic perturbations over the past 1200 years (Table 2). Interpretations of changes in climate between the MCA and LIA support wetter conditions during the MCA (950–1250 C.E.) and drier conditions in the transition between the MCA and LIA (1250–1400 C.E.) and the LIA (1400–1850 C.E.). The proxy analysis of the diatom record independently supports this interpretation of the δ18O first presented by Stansell et al. [16].
Changes in DTF composition appear to be linked to the abundance and intensity of fire with a high amount of burning during the MCA and low burning thereafter c. 1300 C.E. Pre-Columbian anthropogenic impacts on DTF are not detected in the palynological record; however, DTF taxa are noted to decline slightly after European contact (1522 C.E.).
After European contact, and through until c. 1650 C.E., there is an increase in dung fungal spore abundance which possibly represents the increased import of cattle brought into the city of Granada.
It is apparent that over the past 1200 years there have been significant fluctuations in climatic and anthropogenic impacts; however, evidence from El Gancho suggests that changes in the structure of DTF surrounding the Asese peninsula have responded with little consequence. Drought resilient traits of the dominant DTF taxa (e.g., Bursera simaruba and Brosimum) including high wood density and long roots may have helped this community resist climatic and anthropogenic impacts over the past 1200 years.
Rapid agricultural growth from c. 1800 C.E. has allowed for larger populations to be supported; however, it has also increased the reliance on agriculture. Areas that were once covered by DTF that have since been converted to pastures and plantations do appear to be highly threatened by seasonal hydrological fluctuations. An agrarian based society in the CADC is therefore suggested to be far more vulnerable to hydroclimatic changes than a society that primarily relies upon hunting and gathering for sustenance.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2571-550X/2/3/25/s1, Table S1: Images of diatom taxa at 400x magnification identified from El Gancho, Table S2: Images of palynological taxa at 400x magnification identified from El Gancho, Figure S1: Amaranthus spinosus field specimen. Collected on 20 May 2015 from 11°54’22.4”N 85°55’08.4”W.

Author Contributions

W.J.H. conceptualized the presented idea with input from N.S., S.N., and K.J.W. W.J.H. conducted all palaeoecological lab. work, statistical analysis, and age-depth modelling. S.N. and K.J.W. verified the analytical methods. N.S., S.N., and K.J.W. helped supervise the project. W.J.H. drafted the final manuscript. All authors discussed the results contributing to the final manuscript.

Funding

This research was funded by the Natural Environment Research Council of the United Kingdom NE/L002612/1, US National Science Foundation EAR-1502989, the Geological Society of America, and St Edmund Hall, University of Oxford.

Acknowledgments

I would like to give thanks to Julio B. Gomez Martinez, Veronica Lacayo de Gomez and Julio B. Gomez Lacayo, and Maurice Hodgson for their support while conducting field work in Nicaragua, and Gillian Petrokofsky, Peter Long and Carol Adolf from the Oxford Long-Term Ecology Lab., and Ian Matthews from Royal Holloway.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wold Food Programme. Analysis of the Impact of the drought on food security in Guatemala, El Salvador, and Honduras. World Food Programme—Latin America and the Caribbean. Available online: https://documents.wforg/stellent/groups/public/documents/ena/wfp277948.pdf?_ga=2.194671341.1933972640.1559121603-515777643.1559121603 (accessed on 25 April 2018).
  2. Imbach, P.; Locatelli, B.; Zamora, J.C.; Fung, E.; Calderer, L.; Molina, L.; Ciais, P. Impacts of Climate Change on Ecosystem Hydrological Services of Central America: Water Availability. In Climate Change Impacts on Tropical Forests in Central America; Routledge: Abingdon, UK, 2015; pp. 81–106. [Google Scholar]
  3. Durán-Quesada, A.M.; Gimeno, L.; Amador, J. Role of moisture transport for Central American precipitation. Earth Syst. Dyn. 2017, 8, 147. [Google Scholar] [CrossRef]
  4. Hodell, D.A.; Curtis, J.H.; Brenner, M. Possible role of climate in the collapse of Classic Maya civilization. Nature 1995, 375, 391. [Google Scholar] [CrossRef]
  5. Hodell, D.A.; Brenner, M.; Curtis, J.H.; Medina-Gonzalez, R.; Can, E.I.C.; Albornaz-Pat, A.; Guilderson, T.P. Climate change on the Yucatan Peninsula during the little ice age. Quat. Res. 2005, 63, 109–121. [Google Scholar] [CrossRef]
  6. Hodell, D.A.; Brenner, M.; Curtis, J.H. Climate and cultural history of the northeastern Yucatan Peninsula, Quintana Roo, Mexico. Clim. Chang. 2007, 83, 215–240. [Google Scholar] [CrossRef]
  7. Curtis, J.H.; Hodell, D.A.; Brenner, M. Climate variability on the Yucatan Peninsula (Mexico) during the past 3500 years, and implications for Maya cultural evolution. Quat. Res. 1996, 46, 37–47. [Google Scholar] [CrossRef]
  8. Curtis, J.H.; Brenner, M.; Hodell, D.A.; Balser, R.A.; Islebe, G.A.; Hooghiemstra, H. A multi-proxy study of Holocene environmental change in the Maya lowlands of Peten, Guatemala. J. Paleolimnol. 1998, 19, 139–159. [Google Scholar] [CrossRef]
  9. Rosenmeier, M.F.; Hodell, D.A.; Brenner, M.; Curtis, J.H.; Martin, J.B.; Anselmetti, F.S.; Ariztegui, D.; Guilderson, T.P. Influence of vegetation change on watershed hydrology: Implications for paleoclimatic interpretation of lacustrine δ18O records. J. Paleolimnol. 2002, 27, 117–131. [Google Scholar] [CrossRef]
  10. Lachniet, M.S.; Burns, S.J.; Piperno, D.R.; Asmerom, Y.; Polyak, V.J.; Moy, C.M.; Christenson, K. A 1500-year El Niño/Southern Oscillation and rainfall history for the Isthmus of Panama from speleothem calcite. J. Geophys. Res. Atmos. 2004, 109. [Google Scholar] [CrossRef]
  11. Webster, J.W.; Brook, G.A.; Railsback, L.B.; Cheng, H.; Edwards, R.L.; Alexander, C.; Reeder, P.P. Stalagmite evidence from Belize indicating significant droughts at the time of Preclassic Abandonment, the Maya Hiatus, and the Classic Maya collapse. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2007, 250, 1–17. [Google Scholar] [CrossRef]
  12. Metcalfe, S.; Breen, A.; Murray, M.; Furley, P.; Fallick, A.; McKenzie, A. Environmental change in northern Belize since the latest Pleistocene. J. Quat. Sci. 2009, 24, 627–641. [Google Scholar] [CrossRef]
  13. Carrillo-Bastos, A.; Islebe, G.A.; Torrescano-Valle, N.; González, N.E. Holocene vegetation and climate history of central Quintana Roo, Yucatán Península, Mexico. Rev. Palaeobot. Palynol. 2010, 160, 189–196. [Google Scholar] [CrossRef]
  14. Medina-Elizalde, M.; Burns, S.J.; Lea, D.W.; Asmerom, Y.; von Gunten, L.; Polyak, V.; Vuille, M.; Karmalkar, A. High resolution stalagmite climate record from the Yucatán Peninsula spanning the Maya terminal classic period. Earth Planet. Sci. Lett. 2010, 298, 255–262. [Google Scholar] [CrossRef]
  15. Velez, M.I.; Curtis, J.H.; Brenner, M.; Escobar, J.; Leyden, B.W.; Popenoe de Hatch, M. Environmental and cultural changes in highland Guatemala inferred from Lake Amatitlán sediments. Geoarchaeology 2011, 26, 346–364. [Google Scholar] [CrossRef]
  16. Stansell, N.D.; Steinman, B.A.; Abbott, M.B.; Rubinov, M.; Roman-Lacayo, M. Lacustrine stable isotope record of precipitation changes in Nicaragua during the Little Ice Age and Medieval Climate Anomaly. Geology 2013, 41, 151–154. [Google Scholar] [CrossRef]
  17. Kennett, D.J.; Breitenbach, S.F.; Aquino, V.V.; Asmerom, Y.; Awe, J.; Baldini, J.U.; Bartlein, P.; Culleton, B.J.; Ebert, C.; Jazwa, C.; et al. Development and disintegration of Maya political systems in response to climate change. Science 2012, 338, 788–791. [Google Scholar] [CrossRef] [PubMed]
  18. Akers, D.; Brook, G.A.; Railsback, L.B.; Liang, F.; Iannone, G.; Webster, J.W.; Reeder, P.; Cheng, H.; Edwards, R.L. An extended and higher-resolution record of climate and land use from stalagmite MC01 from Macal Chasm, Belize, revealing connections between major dry events, overall climate variability, and Maya sociopolitical changes. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2016, 459, 268–288. [Google Scholar] [CrossRef]
  19. Correa-Metrio, A.; Vélez, M.I.; Escobar, J.; St-Jacques, J.M.; López-Pérez, M.; Curtis, J.; Cosford, J. Mid-elevation ecosystems of Panama: Future uncertainties in light of past global climatic variability. J. Quat. Sci. 2016, 31, 731–740. [Google Scholar] [CrossRef]
  20. Smyth, M.P.; Dunning, N.P.; Weaver, E.M.; van Beynen, P.; Zapata, D.O. The perfect storm: Climate change and ancient Maya response in the Puuc Hills region of Yucatán. Antiquity 2017, 91, 490–509. [Google Scholar] [CrossRef]
  21. Graham, N.E.; Hughes, M.K.; Ammann, C.M.; Cobb, K.M.; Hoerling, M.P.; Kennett, D.J.; Kennett, J.P.; Rein, B.; Stott, L.; Wigand, E.; et al. Tropical Pacific–mid-latitude teleconnections in medieval times. Clim. Chang. 2007, 83, 241–285. [Google Scholar] [CrossRef]
  22. Metcalfe, S.E.; Barron, J.A.; Davies, S.J. The Holocene history of the North American Monsoon: ‘Known knowns’ and ‘known unknowns’ in understanding its spatial and temporal complexity. Quat. Sci. Rev. 2015, 120, 1–27. [Google Scholar] [CrossRef]
  23. Richey, J.N.; Poore, R.Z.; Flower, B.P.; Quinn, T.M. 1400 yr multiproxy record of climate variability from the northern Gulf of Mexico. Geology 2007, 35, 423–426. [Google Scholar] [CrossRef] [Green Version]
  24. Haug, G.H.; Hughen, K.A.; Sigman, D.M.; Peterson, L.C.; Röhl, U. Southward migration of the intertropical convergence zone through the Holocene. Science 2001, 293, 1304. [Google Scholar] [CrossRef] [PubMed]
  25. Dinerstein, E.; Olson, D.; Joshi, A.; Vynne, C.; Burgess, N.D.; Wikramanayake, E.; Hahn, N.; Palminteri, S.; Hedao, P.; Noss, R.; et al. An ecoregion-based approach to protecting half the terrestrial realm. BioScience 2017, 67, 534–545. [Google Scholar] [CrossRef] [PubMed]
  26. Ford, A.; Nigh, R. Origins of the Maya forest garden: Maya resource management. J. Ethnobiol. 2009, 29, 213–236. [Google Scholar] [CrossRef]
  27. Dull, R.A. An 8000-year record of vegetation, climate, and human disturbance from the Sierra de Apaneca, El Salvador. Quat. Res. 2004, 61, 159–167. [Google Scholar] [CrossRef]
  28. Dull, R.A. A Holocene record of Neotropical savanna dynamics from El Salvador. J. Paleolimnol. 2004, 32, 219–231. [Google Scholar] [CrossRef]
  29. Dull, R.A. Evidence for Forest Clearance, Agriculture, and Human-Induced Erosion in Precolumbian El Salvador. Ann. Assoc. Am. Geogr. 2007, 97, 127–141. [Google Scholar] [CrossRef]
  30. Willis, K.J.; Jeffers, E.S.; Tovar, C. What makes a terrestrial ecosystem resilient? Science 2018, 359, 988–989. [Google Scholar] [CrossRef]
  31. Murphy, P.G.; Lugo, A.E. Ecology of tropical dry forest. Annual. Rev. Ecol. Syst. 1986, 17, 67–88. [Google Scholar] [CrossRef]
  32. Janzen, D. Management of habitat fragments in a tropical dry forest: Growth. Ann. Mo. Bot. Gard. 1988, 75, 105–116. [Google Scholar] [CrossRef]
  33. Gillespie, T.W.; Grijalva, A.; Farris, C.N. Diversity, composition, and structure of tropical dry forests in Central America. Plant Ecol. 2000, 147, 37–47. [Google Scholar] [CrossRef]
  34. Miles, L.; Newton, A.C.; DeFries, R.S.; Ravilious, C.; May, I.; Blyth, S.; Kapos, V.; Gordon, J.E. A global overview of the conservation status of tropical dry forests. J. Biogeogr. 2006, 33, 491–505. [Google Scholar] [CrossRef]
  35. Allen, K.; Dupuy, J.M.; Gei, M.G.; Hulshof, C.; Medvigy, D.; Pizano, C.; Salgado-Negret, B.; Smith, C.M.; Trierweiler, A.; Van Bloem, S.J.; et al. Will seasonally dry tropical forests be sensitive or resistant to future changes in rainfall regimes? Environ. Res. Lett. 2017, 12, 2. [Google Scholar] [CrossRef]
  36. Carpenter, S.; Walker, B.; Anderies, J.M.; Abel, N. From Metaphor to Measurement: Resilience of What to What? Ecosystems 2001, 4, 765–781. [Google Scholar] [CrossRef]
  37. Chiabai, A. (Ed.) Climate Change Impacts on Tropical Forests in Central America: An Ecosystem Service Perspective; Routledge: Abingdon, UK, 2015. [Google Scholar]
  38. Guevara-Murua, A.; Williams, C.A.; Hendy, E.J.; Imbach, P. 300 years of hydrological records and societal responses to droughts and floods on the Pacific coast of Central America. Clim. Past 2018, 14, 175–191. [Google Scholar] [CrossRef] [Green Version]
  39. Taylor, B.W. An outline of the vegetation of Nicaragua. J. Ecol. 1963, 51, 27–54. [Google Scholar] [CrossRef]
  40. González, R.B. Tree Species Diversity and Regeneration of Tropical Dry Forest in Nicaragua. Ph.D. Thesis, Swedish University of Agricultural Sciences, Umea, Sweden, 2005. [Google Scholar]
  41. Harcourt, C.S.; Sayer, J.A. The Conservation Atlas of Tropical Forest. The Americas; Simon & Schuster: New York, NY, USA, 1996; pp. 206–211. [Google Scholar]
  42. Denevan, W.M. The pristine myth: The landscape of the Americas in 1492. Annu. Assoc. Am. Geogr. 1992, 82, 369–385. [Google Scholar] [CrossRef]
  43. Cooke, R.; Ranere, A.J. Prehistoric human adaptations to the seasonally dry forests of Panama. World Archaeol. 1992, 24, 114–133. [Google Scholar] [CrossRef]
  44. Daniels, A.E.; Painter, K.; Southworth, J. Milpa imprint on the tropical dry forest landscape in Yucatan, Mexico: Remote sensing and field measurement of edge vegetation. Agric. Ecosyst. Environ. 2008, 123, 293–304. [Google Scholar] [CrossRef]
  45. Newson, L.A. Indian Survival in Colonial Nicaragua; University of Oklahoma Press: Norman, OK, USA, 1987. [Google Scholar]
  46. Denevan, W.M. The Native Population of the Americas in 1492; University of Wisconsin Press: Madison, WI, USA, 1987. [Google Scholar]
  47. Randell, D.H. The Indian Slave Trade and Population of Nicaragua during the Sixteenth Century. In The Native Population of the Americas in 1942; Denevan, W.M., Ed.; University of Wisconsin Press: Madison, WI, USA, 1976. [Google Scholar]
  48. Lovell, W.G.; Lutz, C.H. The Historical Demography of Colonial Central America. In Yearbook: Conference of Latin Americanist Geographers Conference of Latin Americanist Geographers; University of Texas Press: Austin, TX, USA, 1992; pp. 127–138. [Google Scholar]
  49. Kroeber, A.L. Cultural and Natural Areas of Native North America; University of California Publications in American Archaeology and Ethnology; University of California Press: Berkeley, CA, USA, 1939; Volume 38. [Google Scholar]
  50. MacLeod, M.J. Spanish Central America: A Socioeconomic History; University of California Press: Berkeley, CA, USA, 2010; pp. 1520–1720. [Google Scholar]
  51. Newson, L.A. The Cost of Conquest: Indian Decline in Honduras under Spanish Rule; Dellplain Latin American Studies; Westview Press: Boulder, CO, USA, 1986. [Google Scholar]
  52. Lange, F.W.; Haberland, W. The Archaeology of Pacific Nicaragua; University of New Mexico Press: Albuquerque, NM, USA, 1992. [Google Scholar]
  53. Healy, P.; Pohl, M. Archaeology of the Rivas Region, Nicaragua; Wilfrid Laurier University Press: Waterloo, ON, Canada, 1980. [Google Scholar]
  54. Vinogradov, A. Demography: The Population of the Countries of the World from Most Ancient Times to the Present Days; Statistical Tables, Russian Edition; CreateSpace Independent Publishing Platform: Scotts Valley, CA, USA; Volume 1, ISBN 1496106733.
  55. Worldometers. Nicaragua Population. 2014. Available online: http://www.worldometers.info/world-population/nicaragua-population/ (accessed on 25 April 2018).
  56. Janzen, D. Costa Rican Natural History; The University of Chicago Press: Chicago, IL, USA, 1983. [Google Scholar]
  57. Calvo-Alvarado, A.; McLennan, B.; Sanchez-Azofeifa, A.; Garvin, T. Deforestation and forest restoration in Guanacaste, Costa Rica: Putting conservation policies in context. J. For. Ecol. Manag. 2009, 258, 931–940. [Google Scholar] [CrossRef]
  58. Griscom, H.; Ashton, M.S. Restoration of dry tropical forests in Central America: A review of pattern and process. For. Ecol. Manag. 2011, 261, 1564–1579. [Google Scholar] [CrossRef]
  59. Murphy, G.; Lugo, A.E. Dry forests of Central America and the Caribbean. Seas. Dry Trop. For. 1995, 9–34. [Google Scholar] [CrossRef]
  60. Heckadon Moreno, S. Panama’s Expanding Cattle Front: The Santeno Campesinos and the Colonization of the Forests. Ph.D. Thesis, University of Essex, Colchester, UK, 1984. [Google Scholar]
  61. Sabogal, C. Regeneration of tropical dry forests in Central America, with examples from Nicaragua. J. Veg. Sci. 1992, 3, 407–414. [Google Scholar] [CrossRef]
  62. Roldan, H. Recursos Forestales y Cambio en el Uso de la Tierra; Republica de Nicaragua: Santiago, Chile, 2001; p. 73. [Google Scholar]
  63. Van Wyk de Vries, B.; Francis, W. Catastrophic collapse at stratovolcanoes induced by gradual volcano spreading. Nature 1997, 387, 387–390. [Google Scholar] [CrossRef]
  64. Stansell, N.D. Radiocarbon ages for the timing of debris avalanches at Mombacho Volcano, Nicaragua. Bull. Volcanol. 2013, 75, 1–4. [Google Scholar] [CrossRef]
  65. Montenegro-Guillén, S. Lake Cocibolca/Nicaragua: Lake Basin Management Initiative Experience and Lessons Learned Brief USAID; World Bank: Washington, DC, USA, 2003; p. 12. [Google Scholar]
  66. Beard, J.S. Montane vegetation in the Antilles. Caribb. For. 1942, 3, 61–74. [Google Scholar]
  67. Beard, J.S. Climax vegetation in tropical America. Ecology 1944, 25, 127–158. [Google Scholar] [CrossRef]
  68. Atwood, J.T. A floristic study of volcán Mombacho department of Granada, Nicaragua. Ann. MO Botanical Garden 1984, 71, 191–209. [Google Scholar] [CrossRef]
  69. Bellefontaine, R.; Petrucci, Y. Management of Natural Forests of Dry Tropical Zones: FAO Conservation Guide 32; Food and Agriculture Organization of the United Nations: Rome, Italy, 2000; Available online: http://www.fao.org/3/w4442e/w4442e00.htm (accessed on 16 July 2019).
  70. Biondi-Morra, B.N. Hungry Dreams: The Failure of Food Policy in Revolutionary Nicaragua, 1979–1990; Cornell University Press: Ithaca, NY, USA, 1993. [Google Scholar]
  71. UNIDA. Población Total 1993, Estimada al 30 de Junio del Año. Available online: http://www.inide.gob.ni/estadisticas/Cifras%20municipales%20a%C3%B1o%202012%20INIDE.pdf (accessed on 22 August 2018).
  72. Livingstone, D.A. A lightweight piston sampler for lake deposits. Ecology 1955, 36, 137–139. [Google Scholar] [CrossRef]
  73. Reimer, J.; Bard, E.; Bayliss, A.; Beck, J.W.; Blackwell, G.; Ramsey, C.B.; Brown, D.M.; Buck, C.E.; Edwards, R.L.; Friedrich, M.; et al. Selection and treatment of data for radiocarbon calibration: An update to the International Calibration (IntCal) criteria. Radiocarbon 2013, 55, 1923. [Google Scholar] [CrossRef]
  74. Berglund, B.E. Handbook of Holocene Palaeoecology and Palaeohydrology. J. Ecol. 1987, 75, 755. [Google Scholar]
  75. OxLEL, Diatom Preparation Procedure. Oxford Long-Term Ecology Lab. Protocols. Available online: https://oxlel.zoo.ox.ac.uk/wp-content/uploads/2014/04/OxLEL-diatom-protocol.pdf (accessed on 13 June 2018).
  76. Fleming, W.D. Naphrax: A synthetic mounting medium of high refractive Index new and improved methods of preparation. J. R. Microsc. Soc. 1954, 74, 42–44. [Google Scholar] [CrossRef] [PubMed]
  77. Gasse, F. East African Diatoms: Taxonomy, Ecological Distribution; Cramer: Berlin, Germany, 1986; Volume 11. [Google Scholar]
  78. Round, F.E.; Crawford, R.M.; Mann, D.G. Diatoms: Biology and Morphology of the Genera; Cambridge University Press: Cambridge, UK, 1990. [Google Scholar]
  79. Krammer, K.; Lange-Bertalot, H. Bacillariophyceae: Teil 3: Centrales, Fragilariaceae, Eunotiaceae (Suesswasserflora von Mitteleuropa); Spektrum Akademischer Verlag: Berlin, Germany, 1991. [Google Scholar]
  80. Krammer, K.; Lange-Bertalot, H. Bacillariophyceae: Teil 4: Achnanthaceae, Kritische Erganzungen Zu Achnanthes S.L., Navicula S.Str., Gomphonema (Suesswasserflora Von Mitteleuropa); Spektrum Akademischer Verlag: Berlin, Germany, 1991. [Google Scholar]
  81. Krammer, K.; Lange-Bertalot, H. Bacillariophyceae: 2. Teil: Bacillariaceae, Epithemiaceae, Surirellaceae (Suesswasserflora von Mitteleuropa); Spektrum Akademischer Verlag: Berlin, Germany, 1997. [Google Scholar]
  82. Lange-Bertalot, H.; Krammer, K. Süßwasserflora von Mitteleuropa, Bd. 02/1: Bacillariophyceae, 1. Teil: Naviculaceae, A: Text; B: Tafeln; Spektrum Akademischer Verlag: Berlin, Germany, 1997. [Google Scholar]
  83. Lange-Bertalot, H.; Krammer, K.; Bate, N. Süßwasserflora Von Mitteleuropa, Bd. 02/5: Bacillariophyceae: Teil 5: English And French Translation Of The Keys; Spektrum Akademischer Verlag: Berlin, Germany, 2000. [Google Scholar]
  84. OxLEL, Pollen Preparation Procedure. Oxford Long-Term Ecology Lab. Protocols. Available online: https://oxlel.zoo.ox.ac.uk/wp-content/uploads/2016/10/OxLEL-Fossil-pollen-extraction-protocol.pdf (accessed on 13 June 2018).
  85. Leyden, B.W. Pollen evidence for climatic variability and cultural disturbance in the Maya lowlands. Anc. Mesoam. 2002, 13, 85–101. [Google Scholar] [CrossRef]
  86. Bush, M.B.; Weng, C. Introducing a new (freeware) tool for palynology. J. Biogeogr. 2007, 34, 377–380. [Google Scholar] [CrossRef]
  87. APSA. The Australasian Pollen and Spore Atlas V1.0; Australian National University: Canberra, Australia, 2007; Available online: http://apsa.anu.edu.au/ (accessed on 13 June 2018).
  88. Martin, A.C.; Harvey, W.J. The Global Pollen Project: A new tool for pollen identification and the dissemination of physical reference collections. Methods Ecol. Evol. 2017, 8, 892–897. [Google Scholar] [CrossRef]
  89. Roubik, D.W.; Moreno, P. Pollen and spores of Barro Colorado Island (Panama). Pollen Spores Barro Colo. Isl. 1991, 4, 36. [Google Scholar]
  90. Willard, D.A.; Bernhardt, C.E.; Weimer, L.; Cooper, S.R.; Gamez, D.; Jensen, J. Atlas of pollen and spores of the Florida Everglades. Palynology 2004, 28, 175–227. [Google Scholar] [CrossRef]
  91. Bush, M.B. Neotropical plant reproductive strategies and fossil pollen representation. Am. Nat. 1995, 145, 594–609. [Google Scholar] [CrossRef]
  92. Bensch, C.N.; Horak, M.J.; Peterson, D. Interference of redroot pigweed (Amaranthus retroflexus L.), Palmer amaranth (A. palmeri), and common waterhemp (A. rudis) in soybean. Weed Sci. 2003, 51, 37–43. [Google Scholar] [CrossRef]
  93. Davis, O.K.; Shafer, D.S. Sporormiella fungal spores, a palynological means of detecting herbivore density. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2006, 237, 40–50. [Google Scholar] [CrossRef]
  94. Baker, A.G.; Cornelissen, P.; Bhagwat, S.A.; Vera, F.W.; Willis, K.J. Quantification of population sizes of large herbivores and their long-term functional role in ecosystems using dung fungal spores. Methods Ecol. Evol. 2016, 7, 1273. [Google Scholar] [CrossRef]
  95. Gavin, D.G.; Brubaker, L.B.; Lertzman, K.P. An 1800-year record of the spatial and temporal distribution of fire from the west coast of Vancouver Island, Canada. Can. J. For. Res. 2003, 33, 573–586. [Google Scholar] [CrossRef]
  96. Lynch, J.A.; Clark, J.S.; Stocks, B.J. Charcoal production, dispersal, and deposition from the Fort Providence experimental fire: Interpreting fire regimes from charcoal records in boreal forests. Can. J. For. Res. 2004, 34, 1642. [Google Scholar] [CrossRef]
  97. Peters, M.E.; Higuera, P.E. Quantifying the source area of macroscopic charcoal with a particle dispersal model. Quat. Res. 2007, 67, 304–310. [Google Scholar] [CrossRef]
  98. Higuera, E.; Peters, M.E.; Brubaker, L.B.; Gavin, D.G. Understanding the origin and analysis of sediment-charcoal records with a simulation model. Quat. Sci. Rev. 2007, 26, 1790. [Google Scholar] [CrossRef]
  99. Higuera, E.; Gavin, D.G.; Bartlein, J.; Hallett, D.J. Peak detection in sediment–charcoal records: Impacts of alternative data analysis methods on fire-history interpretations. Int. J. Wildland Fire 2011, 19, 996. [Google Scholar] [CrossRef]
  100. Anderson, L.; Wahl, D. Two Holocene paleofire records from Peten, Guatemala: Implications for natural fire regime and prehispanic Maya land use. Glob. Planet. Chang. 2016, 138, 82–92. [Google Scholar] [CrossRef]
  101. Whitlock, C.; Anderson, R.S. Fire History Reconstructions Based on Sediment Records from Lakes and Wetlands. In Fire and Climatic Change in Temperate Ecosystems of the Western AMERICAS; Springer: New York, NY, USA, 2003; pp. 3–31. [Google Scholar]
  102. Maher, L.J. Statistics for microfossil concentration measurements employing samples spiked with marker grains. Rev. Palaeobot. Palynol. 1981, 32, 153–191. [Google Scholar] [CrossRef]
  103. Bennett, J.P. Concepts of mathematical modeling of sediment yield. Water Resour. Res. 1974, 10, 485–492. [Google Scholar] [CrossRef]
  104. Bennett, K.D.; Willis, K.J. Pollen. In Tracking Environmental Change Using Lake Sediments: Terrestrial, Algal, and Siliceous Indicators; Springer: Dordrecht, The Netherlands, 2001; Volume 3. [Google Scholar]
  105. Whitlock, C.; Larsen, C. Charcoal as a Fire Proxy. In Tracking Environmental Change Using Lake Sediments; Springer: Dordrecht, The Netherlands, 2002; pp. 75–97. [Google Scholar]
  106. Bennett, K.D. Determination of the number of zones in a biostratigraphical sequence. New Phytol. 1996, 132, 155–170. [Google Scholar] [CrossRef]
  107. Ter Braak, C.J.; Prentice, I.C. A Theory of Gradient Analysis. Adv. Ecol. Res. 1988, 18, 271–317. [Google Scholar]
  108. Murdoch, D.J.; Chow, E.D. A graphical display of large correlation matrices. Am. Stat. 1996, 50, 178–180. [Google Scholar]
  109. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2012. [Google Scholar]
  110. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, R.; O’hara, R.B.; Simpson, G.L.; Solymos, P.; Stevens, M.H.H.; Wagner, H.; et al. Package ‘vegan’. Available online: https://cran.ism.ac.jp/web/packages/vegan/vegan.pdf (accessed on 25 April 2018).
  111. Juggins, S. Rioja: Analysis of Quaternary Science Data, R Package Version 0.9-5, The Comprehensive R Archive Network. 2015. Available online: https://cran.r-project.org/web/packages/rioja/index.html (accessed on 25 April 2018).
  112. Brown, K.J.; Hebda, R.J. Coastal rainforest connections disclosed through a Late Quaternary vegetation, climate, and fire history investigation from the Mountain Hemlock Zone on southern Vancouver Island, British Colombia, Canada. Rev. Palaeobot. Palynol. 2003, 123, 247–269. [Google Scholar] [CrossRef]
  113. Bush, M.B.; De Oliveira, E.; Colinvaux, A.; Miller, M.C.; Moreno, J.E. Amazonian paleoecological histories: One hill, three watersheds. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2004, 214, 359–393. [Google Scholar] [CrossRef]
  114. Figueroa-Rangel, B.L.; Willis, K.J.; Olvera-Vargas, M. Late-Holocene successional dynamics in a transitional forest of west-central Mexico. Holocene 2012, 22, 143–153. [Google Scholar] [CrossRef]
  115. Wang, L.; Lu, H.; Liu, J.; Gu, Z.; Mingram, J.; Chu, G.; Li, J.; Rioual, P.; Negendank, J.F.; Han, J.; et al. Diatom-based inference of variations in the strength of Asian winter monsoon winds between 17,500 and 6000 calendar years BP. J. Geophys. Res. Atmos. 2008, 113. [Google Scholar] [CrossRef]
  116. McCafferty, G.G.; McCafferty, S.D. Bling Things: Ornamentation and Identity in Pacific Nicaragua. In Identity Crisis: Archaeological Perspectives on Social Identity; Chacmool Archaeological Association: Calgary, AB, Canada, 2011; pp. 243–252. [Google Scholar]
  117. Roman-Lacayo, M. Social and Environmental Risk and the Development of Social Complexity in Precolumbian Masaya, Nicaragua. Ph.D. Thesis, University of Pittsburgh, Pittsburgh, PA, USA, 2013. [Google Scholar]
  118. Avnery, S.; Dull, R.A.; Keitt, T.H. Human versus climatic influences on late-Holocene fire regimes in southwestern Nicaragua. Holocene 2011, 21, 699–706. [Google Scholar] [CrossRef]
  119. White, C.D. Reconstructing Ancient Maya Diet; University of Utah Press: Salt Lake City, UT, USA, 1999. [Google Scholar]
  120. Schnetzler, K.A. Food Uses and Amaranth Product Research: A Comprehensive Review. In Amaranth Biology, Chemistry, and Technology; CRC Press: Boca Raton, FL, USA, 2018; pp. 155–184. [Google Scholar]
  121. Rivera, D.; Inocencio, C.; Obon, C.; Carreno, E.; Reales, A.; Alcaraz, F. Archaeobotany of capers (Capparis)(Capparaceae). Veg. Hist. Archaeobot. 2002, 11, 295–314. [Google Scholar] [CrossRef]
  122. Kaimowitz, D. Policies Affecting Deforestation for Cattle in Central America. In Sustainable Agriculture in Central America; Palgrave Macmillan: London, UK, 1997; pp. 51–66. [Google Scholar]
  123. McNeil, C.L. Deforestation, agroforestry, and sustainable land management practices among the Classic period Maya. Quat. Int. 2012, 249, 19–30. [Google Scholar] [CrossRef]
  124. Burns, E.B. Patriarch and Folk: The Emergence of Nicaragua, 1798–1858; Harvard University Press: Cambridge, MA, USA, 1991. [Google Scholar]
  125. Cole, L.E.; Bhagwat, S.A.; Willis, K.J. Recovery and resilience of tropical forests after disturbance. Nat. Commun. 2014, 5, 3906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  126. Li, D.; Wu, S.; Liu, L.; Zhang, Y.; Li, S. Vulnerability of the global terrestrial ecosystems to climate change. Glob. Chang. Biol. 2018, 24, 4095. [Google Scholar] [CrossRef] [PubMed]
  127. Seddon, A.W.; Macias-Fauria, M.; Long, R.; Benz, D.; Willis, K.J. Sensitivity of global terrestrial ecosystems to climate variability. Nature 2016, 531, 229. [Google Scholar] [CrossRef] [PubMed]
  128. Chadwick, R.; Good, P.; Martin, G.; Rowell, D.P. Large rainfall changes consistently projected over substantial areas of tropical land. Nat. Clim. Chang. 2015, 6, 177–181. [Google Scholar] [CrossRef]
  129. Feng, X.; Porporato, A.; Rodriguez-Iturbe, I. Changes in rainfall seasonality in the tropics. Nat. Clim. Chang. 2013, 3, 811. [Google Scholar] [CrossRef]
  130. Greve, P.; Orlowsky, B.; Mueller, B.; Sheffield, J.; Reichstein, M.; Seneviratne, S.I. Global assessment of trends in wetting and drying over land. Nat. Geosci. 2014, 7, 716. [Google Scholar] [CrossRef]
  131. Brodribb, T.J.; Holbrook, N.M.; Gutierrez, M.V. Hydraulic and photosynthetic co-ordination in seasonally dry tropical forest trees. Plant Cell Environ. 2002, 25, 1435. [Google Scholar] [CrossRef]
  132. Scholz, G.; Liebner, F.; Koch, G.; Bues, C.T.; Günther, B.; Bäucker, E. Chemical, anatomical and technological properties of Snakewood (Brosimum guianense (Aubl.) Huber). Wood Sci. Technol. 2007, 41, 673. [Google Scholar] [CrossRef]
  133. Borchert, R. Soil and stem water storage determine phenology and distribution of tropical dry forest trees. Ecology 1994, 75, 1437. [Google Scholar] [CrossRef]
  134. David, T.S.; Pinto, C.A.; Nadezhdina, N.; David, J.S. Water and forests in the Mediterranean hot climate zone: A review based on a hydraulic interpretation of tree functioning. For. Syst. 2016, 25, 1–14. [Google Scholar] [CrossRef]
Figure 1. Nicaraguan region of the Central American Dry Corridor (ellipsis) and location of Granada, the Asese peninsula and El Gancho.
Figure 1. Nicaraguan region of the Central American Dry Corridor (ellipsis) and location of Granada, the Asese peninsula and El Gancho.
Quaternary 02 00025 g001
Figure 2. El Gancho lake level and vegetation cover. (A) Pre-El Nino Southern Oscillation: March 2014; (B) El Nino Southern Oscillation: May 2015; (C) Post- El Nino Southern Oscillation: March2017.
Figure 2. El Gancho lake level and vegetation cover. (A) Pre-El Nino Southern Oscillation: March 2014; (B) El Nino Southern Oscillation: May 2015; (C) Post- El Nino Southern Oscillation: March2017.
Quaternary 02 00025 g002
Figure 3. Age-depth model for El Gancho.
Figure 3. Age-depth model for El Gancho.
Quaternary 02 00025 g003
Figure 4. Diatom assemblage reconstructed from El Gancho with δ18O data from Stansell et al. [16].
Figure 4. Diatom assemblage reconstructed from El Gancho with δ18O data from Stansell et al. [16].
Quaternary 02 00025 g004
Figure 5. Principal component analysis of the diatom assemblage (A) Canonical correspondence analysis of the diatom assemblage (B) Zones derived from the broken stick model conducted on the palynological assemblage represented by ellipses at a confidence of 95%. Zone 1 = green triangles; Zone 2 = yellow squares; Zone 3 = blue circles.
Figure 5. Principal component analysis of the diatom assemblage (A) Canonical correspondence analysis of the diatom assemblage (B) Zones derived from the broken stick model conducted on the palynological assemblage represented by ellipses at a confidence of 95%. Zone 1 = green triangles; Zone 2 = yellow squares; Zone 3 = blue circles.
Quaternary 02 00025 g005
Figure 6. Palynological assemblage, Sporormiella, and macroscopic charcoal reconstructed from El Gancho, Nicaragua.
Figure 6. Palynological assemblage, Sporormiella, and macroscopic charcoal reconstructed from El Gancho, Nicaragua.
Quaternary 02 00025 g006
Figure 7. Principal component analysis of the palynological assemblage (A) Canonical correspondence analysis of the palynological assemblage (B) Zones derived from the broken stick model represented by ellipses at a confidence of 95%. Zone 1 = green triangles; Zone 2 = yellow squares; Zone 3 = blue circles.
Figure 7. Principal component analysis of the palynological assemblage (A) Canonical correspondence analysis of the palynological assemblage (B) Zones derived from the broken stick model represented by ellipses at a confidence of 95%. Zone 1 = green triangles; Zone 2 = yellow squares; Zone 3 = blue circles.
Quaternary 02 00025 g007
Table 1. Measured, calibrated and modelled 14C ages from El Gancho.
Table 1. Measured, calibrated and modelled 14C ages from El Gancho.
Lab #Measured Age (14C cal B.P.)Depth (cm)2σ Calibrated Age Range (A.D.)2σ Median Calibrated Age2σ OxCal Modelled Age (A.D.)Model Agreement IndexMedian Modelled Age
UCI-19881630± 3582.251286140013431288140599.11346.5
UCI-19882860± 35117.25104612601153105312531051153
UCI-198831100± 30162.258871013950775101290.7893.5
UCI-227661640± 40212.75266538402333530110.6431.5
UCI-227671770± 30226.25123345234186385103285.5
Table 2. Chronology of palynological zones, archaeological periods [52,53,117] local populations [45,46,47,48,51,116,117] and climatic shifts [16].
Table 2. Chronology of palynological zones, archaeological periods [52,53,117] local populations [45,46,47,48,51,116,117] and climatic shifts [16].
Quaternary 02 00025 i001

Share and Cite

MDPI and ACS Style

Harvey, W.J.; Stansell, N.; Nogué, S.; Willis, K.J. The Apparent Resilience of the Dry Tropical Forests of the Nicaraguan Region of the Central American Dry Corridor to Variations in Climate Over the Last C. 1200 Years. Quaternary 2019, 2, 25. https://0-doi-org.brum.beds.ac.uk/10.3390/quat2030025

AMA Style

Harvey WJ, Stansell N, Nogué S, Willis KJ. The Apparent Resilience of the Dry Tropical Forests of the Nicaraguan Region of the Central American Dry Corridor to Variations in Climate Over the Last C. 1200 Years. Quaternary. 2019; 2(3):25. https://0-doi-org.brum.beds.ac.uk/10.3390/quat2030025

Chicago/Turabian Style

Harvey, William J., Nathan Stansell, Sandra Nogué, and Katherine J. Willis. 2019. "The Apparent Resilience of the Dry Tropical Forests of the Nicaraguan Region of the Central American Dry Corridor to Variations in Climate Over the Last C. 1200 Years" Quaternary 2, no. 3: 25. https://0-doi-org.brum.beds.ac.uk/10.3390/quat2030025

Article Metrics

Back to TopTop