ORIGINAL RESEARCH article
Long-Term Vegetation Dynamics in a Megadiverse Hotspot: The Ice-Age Record of a Pre-montane Forest of Central Ecuador
- 1School of Environment, Earth and Ecosystem Sciences, The Open University, Milton Keynes, United Kingdom
- 2Instituto de Ciencias de la Tierra Jaume Almera, Consejo Superior de Investigaciones Científicas, Barcelona, Spain
- 3Estación Biológica de Pindo-Mirador, Universidad Tecnológica Equinoccial, Quito, Ecuador
- 4Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, Amsterdam, Netherlands
Tropical ecosystems play a key role in many aspects of Earth system dynamics currently of global concern, including carbon sequestration and biodiversity. To accurately understand complex tropical systems it is necessary to parameterise key ecological aspects, such as rates of change (RoC), species turnover, dynamism, resilience, or stability. To obtain a long-term (>50 years) perspective on these ecological aspects we must turn to the fossil record. However, compared to temperate zones, collecting continuous sedimentary archives in the lowland tropics is often difficult due to the active landscape processes, with potentially frequent volcanic, tectonic, and/or fluvial events confounding sediment deposition, preservation, and recovery. Consequently, the nature, and drivers, of vegetation dynamics during the last glacial are barely known from many non-montane tropical landscapes. One of the first lowland Amazonian locations from which palaeoecological data were obtained was an outcrop near Mera (Ecuador). Mera was discovered, and analysed, by Paul Colinvaux in the 1980s, but his interpretation of the data as indicative of a forested glacial period were criticised based on the ecology and age control. Here we present new palaeoecological data from a lake located less than 10 km away from Mera. Sediment cores raised from Laguna Pindo (1250 masl; 1°27′S, 78°05′W) have been shown to span the late last glacial period [50–13 cal kyr BP (calibrated kiloyears before present)]. The palaeoecological information obtained from Laguna Pindo indicate that the region was characterised by a relatively stable plant community, formed by taxa nowadays common at both mid and high elevations. Miconia was the dominant taxon until around 30 cal kyr BP, when it was replaced by Hedyosmum, Asteraceae and Ilex among other taxa. Heat intolerant taxa including Podocarpus, Alnus, and Myrica peaked around the onset of the Last Glacial Maximum (c. 21 cal kyr BP). The results obtained from Laguna Pindo support Colinvaux’s hypothesis that glacial cooling resulted in a reshuffling of taxa in the region but did not lead to a loss of the forest structure. Wide tolerances of the plant species occurring to glacial temperature range and cloud formation have been suggested to explain Pindo forest stability. This scenario is radically different than the present situation, so vulnerability of the tropical pre-montane forest is highlighted to be increased in the next decades.
The degree to which the structure and composition of vegetation in tropical South America has been altered in response to high magnitude past global climate change has been long debated (Haffer, 1969; Liu and Colinvaux, 1985; Bush et al., 1990; Absy et al., 1991; Heine, 1994). Revealing the sensitivity of tropical forests to past climate change is the only way in which empirical data can be obtained into how this complex biodiverse region is likely to respond to projected future climate change (Cox et al., 2000; Myers et al., 2000; Malhi and Wright, 2004; IPCC, 2013). Furthermore, it is only by exploring the fossil record that we can parameterise the speed of change that the vegetation has experienced with in the past and consequently gain an idea of the rate at which in may be able to change in the future. Palaeoecology contains powerful tools, such as fossil pollen analysis, with which the dynamics of the vegetation communities through time can be unravelled (Von Post, 1916). The global Last Glacial Maximum (LGM) period (26.5–19 kyr BP; Clark et al., 2009) saw temperatures in the South American tropics of between 4 and 5°C and up to 8°C cooler than modern in the Andes (Bush et al., 2007), and 4–7°C cooler than the modern Amazon (Liu and Colinvaux, 1985; Bush et al., 1990). The LGM-to-modern warming to which tropical South America was subject to over the past 20 ka is, therefore, equivalent to the projected magnitude of change for the next century (IPCC, 2013). In addition, during the last glacial period precipitation (Mosblech et al., 2012) and landscape processes (Loughlin et al., 2018) are likely to have contributed to vegetation change. However, due to a paucity of study sites little is known about the structure and composition of tropical South American glacial vegetation and how it changed during the last glacial period (Colinvaux et al., 1996; Flantua et al., 2015).
The first evidence of glacial vegetation obtained through fossil pollen analysis comes from an outcrop on the Ecuadorian eastern Andean flank near the town of Mera (Liu and Colinvaux, 1985). The glacial sediments from the Mera section were interpreted as containing fossils from a mixture of lowland vegetation and other taxa that live nowadays at higher elevations, and were used to infer a temperature decrease of around 4.5°C compared to present-day. These data were received with scepticism by some researchers, especially regarding: (1) the chronology of the section (Heine, 1994), and (2) the tolerance range of ecological conditions of some of the taxa identified in the record (Gentry, 1993). In fact, some of these debates are still ongoing (Cárdenas et al., 2011a,b; Puyasena et al., 2011). The use of outcrops represents the primary source of sedimentary archives in a very geomorphologically active region (Hall et al., 2008; Lombardo, 2014, 2016). However, outcrops often only represent a short time window and do not contain sediments extending up to the present so consequently interpretation can sometimes be challenging. In addition, outcrops on the eastern Andean flank normally contain interbedded layers of organic (pollen-rich) and inorganic (tephra-like) layers (Cárdenas et al., 2014; Loughlin et al., 2018). This mixture of processes that lead to the sediments’ deposits can easily compromise the continuity of the record by containing numerous sedimentary gaps or hiatuses between the different layers, preventing the study of the vegetation changes in a continuous, dynamic fashion.
Here we present for the first time a continuous vegetation dynamics record from a lacustrine sequence of a mid-elevation Ecuadorian forest located within the diversity hotspot of the eastern Andean flank, in western Amazonia. The glacial dynamics will be explored based on pollen analysis, and supported by charcoal and stable isotope analyses. The aim is to reconstruct the full to late glacial vegetation dynamics, from around 50 to 13 cal kyr BP (calibrated kiloyears before present), as well as to derive potential palaeoclimatic and palaeoecological inferences. The study focuses on the potential changes that might have occurred around the LGM, and it was prompted by the lack of glacial lacustrine palaeoecological studies in the region and the unknown responses of its unique and endangered ecosystem to potential environmental changes. Emphasis will be placed on diversity and stability dynamics, as well as to identify the drivers that have triggered such dynamism. Final details about the sensitivity and resilience of the glacial forest that preceded the current plant community as well as the nature of the forest components will be discussed.
Materials and Methods
Laguna Pindo (1°27′S–78°05′W) is a small shallow lake (c. 1.2 m depth), roughly circular shaped (c. 40 m diameter), located in the Pastaza province near the town of Mera at an elevation of 1248 masl (Figure 1). Mean annual temperature is about 20.8°C with little seasonal variation, annual precipitation can reach up to 4800 mm per year (Ferdon, 1950; Hijmans et al., 2005). Currently the lake is not directly fed by an in-flow and has no visible out-flow; the lake receives water from surface run-off and via direct precipitation, with a rough estimation of a small catchment of around 2–3 km (Figure 1). The study site is positioned in the Andean foothills on a steep slope dropping down to the Pastaza river basin (Figure 1), there are no obvious geomorphological causes for the escarpment of the lake and we hypothesise it is tectonic in origin (Matthews-Bird et al., 2017). The underlying geology of the eastern Andean flank is composed primarily of metamorphic rocks of Palaeozoic to Jurassic age (Aspden and Litherland, 1992; de Berc et al., 2005). The rocks were metamorphosed during the late Cretaceous and the Paleocene, and subsequently were overlain by volcanic/volcaniclastic formations of late Miocene to Quaternary origin (Bernal et al., 2011, 2012). Specifically, the inorganic sediment recovered at the bottom of the sequence has been classified as basalt (Matthews-Bird et al., 2017). Laguna Pindo is a mid-elevation site at the transition between the high elevation páramo vegetation and lowland Amazonia rain forest. The site lies within an area classified as lower montane rain forest (Harling, 1979), or pre-montane forest (Sierra, 1999), just below the lower limit of the cloud forest. The lake is in an advanced stage of filling in, with abundant aquatic plants (Eleocharis maculosa, Cyperaceae) growing within the lake. The lake is completely surrounded by a closed belt of vegetation to the water’s edge. The closed forest surrounding the site has a canopy of 15–25 m high; the dominant species belong to families of Melastomataceae, Araceae, Cecropiaceae, Euphorbiaceae, Myrtaceae, Rubiaceae, Myristicaceae, Asteraceae, and Mimosaceae. Lianas, epiphytes (Bromeliaceae, Orchidaceae), and tree ferns are also common (Jørgensen and León-Yánez, 1999). A survey of the vegetation belt at the lake’s edge is presented in Table 1. The lake is remote and currently beyond the influence of direct human activity such as agriculture and urbanisation, except for the presence of a biological station nearby.
FIGURE 1. Study area. (A) Map of Ecuador, showing the location, main Quaternary active volcanoes near the study site (circles), the speleothem Santiago record (triangle) and other sedimentary archives mentioned in the text (squares): (A) Vinillos (Loughlin et al., 2018), and (B) Mera and Laguna Pindo (the present study). (B) Aerial photograph of the specific study area showing the detailed location and proximity of Laguna Pindo (star) and the closest sequence Mera (square). Main geographical features including river (italics) and towns (bold) are also shown. Year of the photograph: 1990; scale 1:60,000 (bar = 1 km). Photograph obtained from the Ecuadorian Geographic Military Institute (IGM). (C) Picture of the lake where it can be observed the closed vegetation surrounding the lake’s edge (February 2013).
TABLE 1. List of main vegetation taxa currently present surrounding Laguna Pindo based on rough field survey by C. X. Luzuriaga in 2013 and Luzuriaga (2007).
A sediment core was extracted from the deepest point of the lake in January 2013 using a cam-modified Livingstone piston corer (Livingstone, 1955; Colinvaux et al., 1999). The sediment core recovered had a total length of 924 cm. This study presents a multi-proxy investigation (sediment characteristics, pollen, charcoal, and stable isotope) of the lower section of the sediment core (514–924 cm), which has been dated to the last glacial period. Eighteen samples were selected through the entire sequence and sent to the NERC Radiocarbon Facility, SUERC, East Kilbride, Scotland for radiocarbon analysis by accelerator mass spectrometery (Table 2). An age-depth model was constructed using the statistical package “clam” in R (Blaauw, 2010) using the calibration curve SHCal.13.14c (Hogg et al., 2013).
TABLE 2. Conventional (yr BP) and calibrated (cal kyr BP) radiocarbon data used in construction of chronologies for Laguna Pindo.
Samples for pollen analysis (1 cc of wet sediment and 1 cm thickness) were processed using standard methods including KOH, HCl, and HF digestions, acetolysis and mounting/storing in glycerin jelly (Faegri and Iversen, 1989). Lycopodium tablets (University of Lund batch n° 124961; 12,542 spores/tablet) were added before chemical processing (Stockmarr, 1971). Counting was conducted until a minimum of 300 pollen and spores and the saturation of diversity (Rull, 1987). The pollen sum included all pollen types with the exception of aquatic plants (Cyperaceae, Myriophyllum, Sagittaria, Utricularia). Identification was based on the reference collection held at The Open University (United Kingdom), and regional floras and atlases (e.g., Roubik and Moreno, 1991; Colinvaux et al., 1999; Bush and Weng, 2006). Given the diversity of the study area, non-identified morphotypes were coded with the acronym UPP (Unidentified Pollen grain from Pindo). Charcoal particles were identified and counted in the same palynological slides, only particles >5 μm were considered and two different classes were established based on size: (1) small particles (>5–100 μm), indicative of regional fires due to easy dispersion by wind, and (2) big particles (>100 μm), indicative of fires occurred more in a local scale (Whitlock and Larsen, 2001).
Analysis of the stable isotopes δ13C and δ15N was performed at a 4–10 cm sampling interval. Samples for δ13C and δ15N were obtained from ∼0.6 g sample aliquots that were homogenised and treated sequentially with 0.1 and 1 M HCl for 24 h, before being rinsed to neutrality with Milli-Q water (18.2 MΩ cm-1). Each step, involving a change of reagent or water, was preceded by centrifugation to prevent the loss of fine material in suspension. The isotopic composition (δ13C and δ15N) of the dried re-homogenised residues was then determined using a Thermo Flash HT elemental analyser equipped with a Thermo zero-blank device coupled to a Thermo MAT 253 mass spectrometer (EA-MS). Results are expressed following the guidelines for the reporting of stable isotope measurement results (Coplen, 2011).
Pollen diagram, diversity measures, and cluster analyses were performed in R version 3.82 using the packages “vegan” 2.3-5 (Oksanen et al., 2013) and “rioja” (Juggins, 2017). For the cluster analysis the dataset used was the percentage data after square root transformation and the cluster method used was “average” (after calculating the dissimilarity). Zonation was obtained by CONISS using the broken stick method to determine the significant zones (Bennett, 1996). RoC as defined by Urrego et al. (2009) were calculated using the R package PaleoMas (Correa-Metrio et al., 2010). Diversity measures include N0 (richness or species number, also called S), N1 and N2, calculated following Hill (1973). For calculating the indices, the dataset used was the raw data without downweight of rare taxa in order to capture the total diversity values.
Sediment Description and Chronology
The sediment recovered from Laguna Pindo consists mostly of peat and clay with different levels of organic content, and frequent wood remains interbedded in the sediment. Based on the differences found, five sedimentary units were defined (Figure 2 and Table 3). The glacial interval was found in the sedimentary record from 514 cm downwards, corresponding to the Units 4 and 5, and the bottom section of Unit 3, delimited by the presence of a hiatus (Figure 2).
FIGURE 2. Simplified sediment description, age-depth model and accumulation rates calculated. Age-depth model was calculated only considering the organic sediment where proxies were preserved. Dash line in the age-depth model marks the position of the hiatus, and the green band represents an Holocene tephra.
An age-depth model for Laguna Pindo was constructed in Clam.R based on eighteen radiocarbon dates samples (Table 2). The best fit was obtained with a linear interpolation (Figure 2), allowing the calculation of the RoC following Urrego et al. (2009) (Table 4). The sedimentation rate was found to be highly variable, ranging from 0.0009 to 0.25 cm yr-1, with an average of 0.074 cm yr-1 (Figure 2). For the glacial interval, the sedimentation rate found was at least an order of magnitude slower compared to the Holocene section, ranging between 0.02 and 0.003 cm yr-1 for most of the interval, and showing the lowest value of 0.0009 cm yr-1 for the upper part of the glacial section (late Glacial).
TABLE 4. Ecological metrics of Laguna Pindo glacial vegetation based on pollen data: rates of change (RoC) were calculated following Urrego et al. (2009), and diversity indices (N0, N1, N2, and the ratio N2/N0) following Hill (1973).
Pollen grains were found in Laguna Pindo until 822 cm depth, coinciding with the beginning of the sedimentary Unit 5 (Figure 2), which was barren for any biological remain. The pollen diagram of Laguna Pindo during the glacial interval delimits four significant zones based on differences in the most abundant taxa, only taxa occurring at percentages higher than 10% are represented (Figure 3). Charcoal particles have been calculated both as concentration and influx values. The two orders of magnitude difference in the sedimentation rate between the top and bottom of the sequence is generating an artefact in the charcoal influx curve, which is masking the values attained at the top of the sequence (Figure 3). In this sense, results will be described following the concentration curve, although the interpretation of both curves will be provided in the next section.
FIGURE 3. Pollen diagram of the glacial sequence of Laguna Pindo expressed in percentages, including microcharcoal particles results expressed both in concentration (particles cm-3) and influx (particles cm-2 yr-1) values. All the taxa percentages are shown with the same scale (equal width means equal % values) and outline values represent ×10 exaggeration. The grey shade shows the LGM core interval (22–18 cal kyr BP). Taxa have been grouped and coloured following Figure 4.
Pollen Zone PIG-1: From 822 to 756 cm; >50–48.8 cal kyr BP
The oldest section of Laguna Pindo is marked by a decrease to the top of the zone of Hedyosmum and Asteraceae, and in a minor extent Alnus and some unidentified morphotypes such as UPP36, UPP42, and UPP25. At the same time, Miconia shows the opposite trend. Total pollen concentration ranges from 257,000 to 1,000,000 pollen grains cc-1. The charcoal record during this interval is low (Figure 3).
Pollen Zone PIG-2: From 756 to 611 cm; 48.8–37.3 cal kyr BP
PIG-2 shows the inverse relationship between the dominant taxa of the previous zone Hedyosmum and Miconia, with low and high values, respectively (Figure 3). Weinmannia appears during this zone and Alnus, Myrica, Myrsine, Podocarpus, and UPP36 increase the abundance attained in the previous zone. On the contrary, taxa including Sapindaceae, Urticales, and unidentified UPP42 and UPP25 greatly decrease or disappears from the record. Begonia, Ilex, and Melastomataceae are also abundant at the beginning of the zone, but start decreasing towards the upper section. This zone is characterised in the upper half by the maximum values of charcoal, both small and big size (indicative of regional and local fires, respectively). Regarding pollen concentration, the values in this zone are also the highest of the record ranging from 436,000 to 1,279,000 pollen grains cc-1.
Pollen Zone PIG-3: From 611 to 532 cm; 37.3–19.6 cal kyr BP
The most dramatic change of the vegetation surrounding Pindo during glacial time corresponds to the decrease observed in Miconia during this zone to values below 20%, until the middle of the section. Following the previous zone, Hedyosmum shows an opposite character to Miconia, evident in this zone with a subtle but solid increase. Coeval to the minimum value of Miconia, a peak in Ilex is observed, followed by a steadier increase towards the top of the zone. UPP42 and UPP45 reappear during this zone, whereas Papilionidae, Weinmannia, and Myrtaceae decrease. The charcoal curve follows the high values attained at the end of the previous zone until approximately 32 cal kyr BP when they suddenly decrease to half of the particles abundance. Total pollen concentration in PIG-3 ranges from 277,000 to 905,000 pollen grains cc-1.
Pollen Zone PIG-4: From 532 to 514 cm; 19.6–13.5 cal kyr BP
The post-LGM and late Glacial interval of Laguna Pindo is characterised by low values of both Hedyosmum and Miconia compared to the rest of the sequence. UPP25 peaks at the beginning of the zone and disappears again. Asteraceae and Ilex attain during this zone their highest values, and Myrtaceae shows its lowest abundance. Charcoal particles present in this zone are among the minimum values of the entire sequence, and the same occurs with the pollen concentration, ranging from 211,000 to 544,000 pollen grains cc-1.
Additional Metrics for Plant Dynamics
Several tests were also run to get a better idea of the palynological dynamics. A cluster analysis of the represented taxa (percentages above 10%) was performed to see the grouping formed by the most abundant pollen morphotypes with similar distributions along the sequence (Figure 4). The plot shows three different groups well defined, with only Miconia belonging to the first group; Hedyosmum, Asteraceae, and Ilex forming the second group; and the rest of morphotypes (19) included in Group 3. Stable isotopes of C and N were plotted stratigraphically against age (Figure 5), providing values between -30 and -25‰ for δ13C, indicative of mostly C3 land plants, and between 0 and 5‰ for δ15N, representing a mixed primary production source formed by aquatic and terrestrial plants (Meyers and Teranes, 2001). δ13C sees an increasing trend starting around 25 cal kyr BP. δ15N shown an earlier increasing trend around 43 cal kyr BP including a brief drop around 19 cal kyr BP. Finally, some diversity metrics were calculated following Hill (1973), as previous works have highlighted their suitability for pollen data (Finsinger et al., 2017; Gosling et al., 2017). Based on a total of 134 recognisable pollen morphotypes for the glacial palynological assemblage of Laguna Pindo, samples N0 vary between 40 and 57 morphotypes, N1 and N2 values rank from 17 up to 24 and from 7 to 18, respectively, and the ratio N2/N0 ranges between 0.13 and 0.39 (Table 4).
FIGURE 4. Cluster analysis through the “average” method of the most abundant pollen taxa found in Laguna Pindo during the glacial interval. Colours define the groups resulted: Group 1 (blue), Group 2 (orange), and Group 3 (green).
FIGURE 5. Laguna Pindo metrics and additional data framed in a regional palaeoclimatic context (raw data downloaded from NOAA): (a) GISP2 temperature reconstruction of Greenland (Alley, 2000), in which are shown the Heinrich (H) and Dansgaard–Oeschger (D/O) events that occurred in the time interval under study; (b) precipitation reconstruction based on Santiago speleothem record (Mosblech et al., 2012); (c) percentages of the groups resulted in the cluster analysis of Laguna Pindo most abundant palynological taxa (Figure 4): with dotted, dashed, and continuous lines showing the different groups (Group 1, 2, and 3 of Figure 4, respectively); (d) local fires curve of Laguna Pindo based on big (>100 μm) charcoal particles, expressed in concentration and influx values (continuous and discontinuous lines, respectively); (e) δ13C of Laguna Pindo, expressed in ‰; and (f) δ15N of Laguna Pindo, expressed in ‰. The dark grey shade represents the hiatus of Laguna Pindo, calculated as the calibrated ages of the dated samples bracketing the sedimentary gap (at 512 and 514 cm depth), and the light grey bands mark the pollen zones obtained. Vertical dashed lines show the LGM core interval following Figure 3.
Palaeoecological Interpretation and Discussion
In order to understand the ecological dynamics of the forests recorded in Laguna Pindo, a palaeoclimatic background is needed. Here, independent archives for temperature (Greenland ice core record) and precipitation (Santiago speleothem, Ecuador) reconstructions will be used and placed in a regional context (Figure 5). A palaeotemperature record from Greenland is shown instead of the closer record of the Cariaco Basin (offshore Venezuela) as both are equally representatives of the North Atlantic Ocean palaeotemperature (Haug et al., 2001; Deplazes et al., 2013), and the Greenland ice record is expressed directly in temperature degrees (Alley, 2000). Vegetation dynamics will be compared with nearby records when possible and in addition, as Laguna Pindo is located in NW Amazonia at the boundary between Amazon and Andean forests just below the narrow band of cloud forests, comparison with long records retrieved in ecotonal areas from SW and NE Amazonia will be also contemplated.
Glacial Vegetation at Laguna Pindo
The small catchment size of Laguna Pindo and the persistent abundance of woody taxa (based on pollen and isotopic signals) throughout the glacial period suggest that the Mera region was continuously covered by forest during this period (Figures 3, 5). The large number of pollen taxa found in the glacial record from Laguna Pindo (>130 terrestrial morphotypes) reflects the high biodiversity that characterises western Amazonia and the eastern Andean flank. In this sense, the palynological assemblage of the glacial Pindo forest shows similar values of diversity (N1) to modern pollen traps of tropical locations in Bolivia and Ghana (Gosling et al., 2017). It is noteworthy to highlight the high values of diversity obtained in such a present-day closed canopy to the shoreline and small catchment of the lake, an environmental setting prone to collect only evidence from a very local spatial scale (Figure 1; Jacobson and Bradshaw, 1981). In addition, it is interesting to observe that despite the homogeneous nature of the richness (N0) values along the sequence, the maximum values for all the indices calculated were attained in a single sample, with an estimated age just at the end of the LGM (Table 4).
The forest taxa occurring had varied through time based on different climatic and ecological requirements, as forests have been doing since previous glacial-interglacial periods in the tropics and elsewhere (Bush et al., 2004a; Cárdenas et al., 2011a). In this sense, full to late glacial forest in Laguna Pindo was characterised by a mix of taxa that live nowadays in both mid and high elevations (>1000 masl). The continuous occurrence of taxa such as Alnus, Hedyosmum, Myrica, Podocarpus, and Weinmannia not present in today’s forest (Table 1) suggest a colder climate than nowadays, whereas the presence of warm indicators like Arecaceae point to the existence of a no-analogue plant community (Williams and Jackson, 2007). Besides the mixed composition of the glacial forest, Laguna Pindo palynological record is in agreement with previous sequences in another important aspect, the community stability and tolerance during the glacial interval. RoC were calculated and provided very low results of magnitude around ×10-4 (Table 4). During full glacial conditions, the structure of the forest did not change compared to those observed nowadays, a pre-montane forest. The vegetation shift in composition is due to the coexistence during the ice age of taxa with different climatic requirements in the present-day.
Dynamics of Vegetation at Laguna Pindo
Laguna Pindo recorded a remarkably continuous lower montane forest cover during the last glaciation, despite several disturbance events or climatic episodes. Disturbance events likely occurred even in the absence of a clear imprint in the sedimentary archive. For instance, a high peak of local fires is located around 43–33 cal kyr BP, or slightly before (with peaks around 52 and 45 cal kyr BP) when expressed in influx values (Figures 3, 5). The occurrence of fires in the wet western Amazonia/eastern Andean flank prior to human arrival were probably caused by volcanic activity (Cárdenas et al., 2011a; Loughlin et al., 2018). However, during this interval no tephras have been observed in Laguna Pindo sequence coeval to the charcoal peak, although volcanic activity is known to have occurred in the region at this time (Keen, 2015; Loughlin et al., 2018). Consequently, a climatic origin for the fires recorded in Laguna Pindo cannot not be entirely ruled out because precipitation levels are thought to have fluctuated in the region throughout the last glacial and both periods with high charcoal in Laguna Pindo record are coeval with drier intervals in Santiago (Mosblech et al., 2012; Figure 5b). Regardless of the origin of the fires, the most conspicuous change observed in the vegetation occurred once these stopped, around 27 cal kyr BP also coeval to a drier interval in Santiago (Figures 5b,d), which occurred during a time period where the speleothem record was decoupled with insolation (Mosblech et al., 2012).
Based on the trends observed in the vegetation dynamics of the upper section of the sequence, late glacial interval was characterised in Laguna Pindo by a possible gradual change towards the Holocene (Figures 3, 5). In this sense, it can be observed that the shifts in taxa abundance were mostly recorded prior or during LGM, and that the dynamics during the late Glacial were minor, until the record stopped at the beginning of the younger Dryas cold reversal (YD; 12.9–11.7 cal kyr BP). The smoothness observed in the diagram curves could be due to either insensitivity or time-lags of the taxa occurring along the lake shore during this time, or because the climatic change itself was mild at Pindo (i.e., not enough for crossing the tolerance thresholds of the occurring species). The nearby Santiago speleothem record (Figures 1, 5b) suggested a wet late glacial interval (Mosblech et al., 2012) whereas the palaeotemperature reconstruction of the Greenland ice record showed a more stable interval (Figure 5a), so it is suggested that both precipitation and temperature could have played a role in the late glacial dynamics of Pindo vegetation. In this sense it can be argued that wet conditions and a stable temperature trend without abrupt extreme values facilitated the late Glacial stability of the vegetation around Pindo, a location that receives >4000 mm yr-1 and has a stable annual temperature around 20°C nowadays. Regarding temperature, comparing the palynological groups and the Greenland temperature reconstruction trends (Figures 5a,c), it can be observed, that, some changes did occur before the LGM. Thus, during the interval around 30–20 cal kyr BP and especially between Heinrich events H3 and H2, Group 1 (Miconia) and 2 (Hedyosmum, Ilex, and Asteraceae) attained their minimum and maximum values, respectively (Figure 5c). Miconia is a genus with more than 200 species with different ecological and climatic tolerances in Ecuador and is very abundant in the surroundings of Laguna Pindo in the present-day (Table 1; Jørgensen and León-Yánez, 1999). However, some heat-intolerant taxa such as Alnus, Podocarpus, Myrica, and Hedyosmum peaked around the LGM core interval (22–18 cal kyr BP; Figure 3). Some of the changes in the abundances of Miconia and the heat intolerant taxa could be temperature-driven. If so, these shifts in abundances would be in agreement with Lago Consuelo record, a cloud forest location in the southern hemisphere (Peru/Bolivia), which also reported a gradual transition from glacial towards Holocene forests, and was preliminary interpreted in terms of temperature rather than precipitation as the major driver of vegetation changes (Bush et al., 2004b).
With respect to the different taxa sensitivity, the key factor to consider here is the proximity to the occurrence of environmental conditions’ thresholds for the given species. Such ecological proximity could be modified through time due to ecosystem interactions, climatic shifts or feedbacks. Longer records have shown forest stability in more than one glaciation (Bush et al., 2004a). Considering the time scales of glacial versus interglacial duration, it is logical to think that the glacial vegetation was more adapted to cold conditions than to the warm characteristic of the Holocene. Moreover, Bush (2002) hypothesised that as the climate warms, the elevation of cloud formation on the Andean flank would have increased. In Lago Consuelo, the authors suggested a higher stability of glacial forest to dry events as a result of cloud cover formation (Urrego et al., 2010). Therefore, it is proposed here that the forest composition of Pindo during the glaciation was the result of complex environmental interactions, including at least these two factors: (i) the forest taxa had climatic tolerances that include the temperature range (and other parameters like CO2 levels) that occurred during the glaciation, and (ii) these lower temperatures of the glacial period facilitated the occurrence of clouds in the wet Laguna Pindo, buffering the effects of a decreased available moisture during the drier intervals. This way, the ecological or climatic threshold was not crossed even during the occurrence of abrupt events resulting in forest stability. The assumption of glacial tolerance range would imply a higher sensitivity to climate change of the forest during interglacials such as the present one. Holocene palaeoecological studies have manifested a high dynamism of plant communities through the entire Amazon basin (e.g., Flantua et al., 2016). Given the current projections on climate change, it is expected that cloud coverage of the eastern Andean flank will continue moving upwards and narrowing as the temperature rises (Bush, 2002; IPCC, 2013). The cloud migration would progressively intensify the vulnerability of large areas of pre-montane and cloud forests, as the buffering effect of the cloud cover moves to higher elevations. The lack of buffer, together with an increased human impact (i.e., increment in land use involving deforestation), would have dramatic consequences for the already threatened biodiversity hotspot of the Ecuadorian eastern Andean flank.
Dynamics of Glacial Vegetation in the Eastern Andean Flank
The first evidence of non-Andean glacial age forests in tropical South America was from a sedimentary sequence obtained by Paul Colinvaux near the town of Mera, Ecuador (Liu and Colinvaux, 1985) just around 10 km from Laguna Pindo (Figure 1). The Mera record was an outcrop exposed by road cutting where a temperature drop of about 4.5°C relative to modern was estimated based on the pollen assemblage found (Liu and Colinvaux, 1985). The estimate of cooling was based on the occurrence of heat-intolerant taxa such as Alnus, Hedyosmum, and Podocarpus. Based on the Mera palynological record, it was suggested that during glacial times tropical forests did not disappear, but were reconfigured and included taxa today only found at much higher elevations. The taxa mixture implied that forests species behaved individually, and not grouped by associations and/or belts, a view those days still not widely accepted outside North America (Whittaker, 1951; van der Hammen, 1974; Davis et al., 1986). Moreover, the glacial forest persistence hypothesis resulting from Mera directly confronted the biogeographic hypothesis of glacial refugia proposed by Haffer (1969) (Colinvaux, 1998; Colinvaux et al., 2000).
The controversial nature of the fossil pollen record from Mera by Liu and Colinvaux (1985) about the structure and composition of the glacial tropical vegetation resulted in several challenges to the interpretation. First of all, the outcrop section contained a large amount of unidentified taxa, as Neotropical palynology was still in its infancy (Colinvaux et al., 1999). Second, the volcanic nature of the inorganic sediment could allow transportation of material from elsewhere, including wood remains of the coniferous tree Podocarpus that were found interbedded in Mera’s sediments (Gentry, 1993). Third, this potentially transported material could be much older than the age reported by radiocarbon dating due to contamination and then the record could not be glacial (Heine, 1994). And finally, the uniqueness of a single-site results located within one of the defined refuge areas prevented a general statement of glacial cooling instead of aridity as the main driver of vegetation changes during the ice-age in the tropical forests of South America (Haffer, 1969). Most of these concerns were addressed by the finding of a second glacial record from another outcrop located at San Juan Bosco, just 160 km to the south from Mera, and the setup of modern pollen traps systematic sampling in different locations of Ecuador and Brazil (Bush et al., 1990, 2001; Bush, 1991; Bush and Weng, 2006). San Juan Bosco record showed a mixed assemblage of pollen grains from taxa nowadays living in mid and high elevations as Mera did (Bush et al., 1990). However, the outcrop nature of the record still compromised the potential occurrence of reworked material from volcanic events transportation and hence contamination with older sediments (Heine, 1994).
Subsequent palaeoecological records recovered from lowland tropical sites in South America have supported Liu and Colinvaux’s assertion that taxa currently found at higher, colder elevations coexisted with modern “warm” lowland taxa in the lowlands during the last glacial period (Haberle and Maslin, 1999; van’t Veer et al., 2000; Bush et al., 2004b; Whitney et al., 2011). The closest of the other glacial records to the Mera and San Juan Bosco sites is Vinillos, which is 100 km north of Pindo (Figure 1), and is also characterised by a very similar pollen glacial forest assemblage (Loughlin et al., 2018). Finally, the glacial palynological assemblage is not restricted to the last glaciation, as it has been found in previous glacial stadials as well in the region (Cárdenas et al., 2011a). However, all these evidence share the outcrop nature of the record, with the dubious provenance of the sediments, the dating limitations and the potential contamination problems. Laguna Pindo is the first lacustrine record on the westernmost edge of the equatorial Amazonia that registers the last glacial dynamics. By being lacustrine in origin and with such a small catchment, scepticism about the continuity of the sediment, the spatial scale information of the data provided and the volcanic transport contamination can be finally ruled out. Given the volume of the glacial records obtained so far for the Ecuadorian Andean flank, there is a general agreement on the heat intolerant taxa found mixed with current heat tolerant were common and abundant during the last ice age in the eastern Andean flank in western Amazonia.
Laguna Pindo sequence contains the story of a pre-montane forest in the biodiversity hotspot of the eastern Andean flank of Ecuador during the last glaciation, for the first time obtained from a continuous lacustrine record. The glacial forest of Laguna Pindo has been described as a mix of taxa living nowadays in mid and high elevations. Heat-intolerant taxa including Podocarpus, Alnus, or Myrica showed maximum values around the start of the LGM, replacing the previous dominant taxon, Miconia. However, the forest was characterised by stability, in contrast to the Holocene dynamism of Amazon plant communities (Flantua et al., 2016). This stability can be understood in terms of forest structure, being a pre-montane forest during full glacial conditions as it is nowadays, instead of compositional, based on the different climatic and ecological requirements of the occurring species. Wide tolerance ranges to glacial conditions and cloud formation have been proposed as key drivers to maintain Pindo stability. Given the ongoing climate change and the potential scenarios projected, it is more than likely that eventually the clouds that act as buffer of unfavourable climatic events such as droughts will disappear and move upwards, increasing the vulnerability of these tropical forests. Hence, it is essential to adopt management strategies and conservation measures in order to keep the currently desired ecosystem services of this biodiversity hotspot. Finally, the mature or “old-grown” forest concept in the eastern Andean flank or western Amazonia should be revised in order to establish accurate future projections and resilience estimates.
EM and WG lead the field work to recover the sediments. HK analysed the charcoal record and created the GIS figure in Figure 1. CL provided a vegetation inventory survey. All authors participate in the discussion during the writing of the manuscript, which was lead by EM.
This study has been developed under the auspices of project FORSENS, funded by the Natural Environment Research Council of United Kingdom (grant NE/J018562/1 to EM). Radiocarbon dating was supported by the NERC Radiocarbon Facility NRCF010001 (allocation number 1682.1112) awarded to FORSENS project by EM and performed by Pauline Gulliver (NRCF East Kilbride).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Field sampling was performed thanks to the required permits for research investigations provided by the Ecuadorian Ministry of Environment (ref. 14-2012-IC-FLO-DPAP-MA), with the support of Susana León-Yánez. Frazer Matthews-Bird helped with core sampling and the core recovery was possible at Laguna Pindo thanks to the Universidad Tecnológica Equinoccial de Ecuador (UTE), the local Council of the Province Pastaza, and the Estación Biológica de Pindo Mirador staff, with special emphasis in Doña Gloria Quichimbo. Special thanks to Emily Sear, Mabs Gilmour, and Simona Nicoara for the technical support, and Stephen Brooks, Mark Bush, and Iain Gilmour for their participation in the project. The work and suggestions of three anonymous referees and the editor Urs Feller greatly improved the quality of the manuscript. This paper is in memorial of Paul Colinvaux, who could not find the Ecuadorian lake of his dreams in the search of the Ice-Age in the American Equator.
Absy, M., Clief, A., Fournier, M., Martin, L., Servant, M., Sifeddine, A., et al. (1991). Mise en evidence de quatre phases d’ouverture de la fôret dense dans le sud-est de L’Amazonie au cours des 60,000 dernières annés: premiere comparaison avec d’autres regions tropicals. C. R. Acad. Sci. Paris 312, 673–678.
Bernal, C., Christophoul, F., Darrozes, J., Soula, J., Baby, P., and Burgos, J. (2011). Late glacial and Holocene avulsions of the Río Pastaza megafan (Ecuador-Peru): frequency and controlling factors. Int. J. Earth Sci. 100, 1759–1782. doi: 10.1007/s00531-010-0555-9
Bernal, C., Christophoul, F., Soula, J., Darrozes, J., Bourrel, L., Laraque, A., et al. (2012). Gradual diversions of the Río Pastaza in the Ecuadorian piedmont of the Andes from 1906 to 2008: role of tectonics, alluvial fan aggradation, and ENSO events. Int. J. Earth Sci. 101, 1913–1928. doi: 10.1007/s00531-012-0752-9
Bush, M. B. (1991). Modern pollen-rain data from South and Central America: a test of the feasibility of fine-resolution lowland tropical palynology. Holocene 1, 162–167. doi: 10.1177/095968369100100209
Bush, M. B., Colinvaux, P. A., Wiemann, M. C., Piperno, D. R., and Liu, K.-B. (1990). Late Pleistocene temperature depression and vegetation change in Ecuadorian Amazonia. Quat. Res. 34, 330–345. doi: 10.1016/0033-5894(90)90045-M
Bush, M. B., De Oliveira, P. E., Colinvaux, P. A., Miller, M. C., and Moreno, J. E. (2004a). Amazonian paleoecological histories: one hill, three watersheds. Palaeogeogr. Palaeoclimatol. Palaeoecol. 214, 359–393. doi: 10.1016/S0031-0182(04)00401-8
Bush, M. B., Hanselman, J. A., and Hooghiemstra, H. (2007). “Andean montane forests and climate change,” in Tropical Rainforest Responses to Climatic Change, 2nd Edn, eds M. B. Bush, J. R. Flenley, and W. D. Gsoling (Chichester: Springer-Praxis), 35–60. doi: 10.1007/978-3-540-48842-2
Bush, M. B., Moreno, E., de Oliveira, P. E., Asanza, E., and Colinvaux, P. A. (2001). The influence of biogeographic and ecological heterogeneity of Amazonian pollen spectra. J. Trop. Ecol. 17, 729–743. doi: 10.1017/S0266467401001547
Cárdenas, M. L., Gosling, W. D., Pennington, R. T., Poole, I., Sherlock, S. C., and Mothes, P. (2014). Forests of the tropical eastern Andean flank during the middle Pleistocene. Palaeogeogr. Palaeoclimatol. Palaeoecol. 393, 76–89. doi: 10.1016/j.palaeo.2013.10.009
Cárdenas, M. L., Gosling, W. D., Sherlock, S. C., Poole, I., Pennington, R. T., and Mothes, P. (2011b). Response to comment on “The response of vegetation on the Andean flank in western Amazonia to Pleistocene climate change”. Science 333:1825. doi: 10.1126/science.1207525
Cárdenas, M. L., Gosling, W. D., Sherlock, S. C., Poole, I., Pennington, R. T., and Mothes, P. (2011a). The response of vegetation on the Andean flank in western Amazonia to Pleistocene climate change. Science 331, 1055–1058. doi: 10.1126/science.1197947
Colinvaux, P. A., De Oliveira, P. E., and Bush, M. B. (2000). Amazonian and neotropical plant communities on glacial time-scales: the failure of the aridity and refuge hypothesis. Quat. Sci. Rev. 19, 141–169. doi: 10.1016/S0277-3791(99)00059-1
Colinvaux, P. A., De Oliveira, P. E., Moreno, J. E., Miller, M. C., and Bush, M. B. (1996). A long pollen record from lowland Amazonia: forest and cooling in glacial times. Science 274, 85–88. doi: 10.1126/science.274.5284.85
Correa-Metrio, A., Cabrera, K. R., and Bush, M. B. (2010). Quantifying ecological change through discriminant analysis: a paleoecological example from the Peruvian Amazon. J. Veg. Sci. 21, 695–704. doi: 10.1111/j.1654-1103.2010.01178.x
Cox, P. M., Betts, R. A., Jones, C. D., Spall, S. A., and Totterdell, I. J. (2000). Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model. Nature 408, 184–187. doi: 10.1038/35041539
de Berc, S. B., Soula, J., Baby, P., Souris, M., Christophoul, F., et al. (2005). Geomorphic evidence of active deformation and uplift in a modern continental wedge-top-foredeep transition: example of the eastern Ecuadorian Andes. Tectonophysics 399, 351–380. doi: 10.1016/j.tecto.2004.12.030
Deplazes, G., Lückge, A., Peterson, L. C., Timmermann, A., Hamann, Y., Hughen, K. A., et al. (2013). Links between tropical rainfall and North Atlantic climate during the last glacial period. Nat. Geosci. 6, 213–217. doi: 10.1038/ngeo1712
Finsinger, W., Morales-Molino, C., Gałka, M., Valsecchi, V., Bojovic, S., and Tinner, W. (2017). Holocene vegetation and fire dynamics at Crveni Potok, a small mire in the Dinaric Alps (Tara National Park, Serbia). Quat. Sci. Rev. 167, 63–77. doi: 10.1016/j.quascirev.2017.04.032
Flantua, S. G. A., Hooghiemstra, H., Grimm, E. C., Behling, H., Bush, M. B., González-Arango, C., et al. (2015). Updated site compilation of the Latin American Pollen Database. Rev. Palaeobot. Palynol. 223, 104–115. doi: 10.1016/j.revpalbo.2015.09.008
Flantua, S. G. A., Hooghiemstra, H., Vuille, M., Behling, H., Carson, J., Gosling, W., et al. (2016). Climate variability and human impact on the environment in South America during the last 2000 years: synthesis and perspectives. Clim. Past 12, 483–523. doi: 10.5194/cp-12-483-2016
Gentry, A. H. (1993). A Field Guide to the Families and Genera of Woody Plant of Northwest South America (Colombia, Ecuador, Peru) with Supplementary Notes on Herbaceous Taxa. Washington, DC: Conservation International.
Gosling, W. D., Julier, A. C. M., Adu-Bredu, S., Djagbletey, G. D., Fraser, W. T., Jardine, P. E., et al. (2017). Pollen-vegetation richness and diversity relationships in the tropics. Veg. Hist. Archaeobot. doi: 10.1007/s00334-017-0642-y
Haberle, S. G., and Maslin, M. A. (1999). Late Quaternary vegetation and climate change in the Amazon basin based on a 50,000-year pollen record from the Amazon Fan, ODP Site 932. Quat. Res. 51, 27–38. doi: 10.1006/qres.1998.2020
Hall, M. L., Samaniego, P., Le Pennec, J. L., and Johnson, J. B. (2008). Ecuadorian Andes volcanism: a review of Late Pliocene to represent activity. J. Volcanol. Geotherm. Res. 176, 1–6. doi: 10.1016/j.jvolgeores.2008.06.012
Haug, G., Hughen, K. A., Sigman, D. M., Peterson, L. C., and Röhl, U. (2001). Southward migration of the Intertropical convergence zone through the holocene. Science 293, 1304–1308. doi: 10.1126/science.1059725
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., and Jarvis, A. (2005). Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978. doi: 10.1002/joc.1276
Hogg, A. G., Hua, Q., Blackwell, P. G., Niu, M., Buck, C. E., Guilderson, T. P., et al. (2013). SHCal13 Southern Hemisphere calibration, 0–50,000 cal yr BP. Radiocarbon 55, 1889–1903. doi: 10.2458/azu_js_rc.55.16783
IPCC (2013). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds T. F. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S. K. Allen, J. Boschung, et al. Cambridge: Cambridge University Press, 1535.
Juggins, S. (2017). rioja: Analysis of Quaternary Science Data, R Package Version (0.9-15). Available at: http://cran.r-project.org/package=rioja
Loughlin, N. J. D., Gosling, W. D., Coe, A. L., Gulliver, P., Mothes, P., and Montoya, E. (2018). Landscape-scale drivers of glacial ecosystem change in the montane forests of the eastern Andean flank. Palaeogeogr. Palaeoclimatol. Palaeoecol. 489, 198–208. doi: 10.1016/j.palaeo.2017.10.011
Matthews-Bird, F., Brooks, S. J., Gosling, W. D., Gulliver, P., Mothes, P., and Montoya, E. (2017). Aquatic community response to volcanic eruptions on the Ecuadorian Andean flank: evidence from the palaeoecological record. J. Paleolimnol. 58, 437–453. doi: 10.1007/s10933-017-0001-0
Meyers, P. A., and Teranes, J. L. (2001). “Sediment organic matter,” in Tracking Environmental Change Using Lake Sediments. Physical and Geochemical Methods, Vol. 2, eds J. P. Smol, H. J. B. Birks, and W. M. Last (Dordrecht: Kluwer), 239–270. doi: 10.1007/0-306-47670-3_9
Mosblech, N. A. S., Bush, M. B., Gosling, W. D., Hodell, D., Thomas, L., van Calsteren, P., et al. (2012). North Atlantic forcing of Amazonian precipitation during the last ice age. Nat. Geosci. 5, 817–820. doi: 10.1038/ngeo1588
Oksanen, J., Guillaume Blanchet, F. R. K., Legendre, P., Minchin, P. R., O’Hara, B. R., Simpson, G. L., et al. (2013). Vegan: Community Ecology Package. R Package Version 2.0-10. Available at: https://cran.r-project.org/package=vegan
Puyasena, S. W., Dalling, J. W., Jaramillo, C., and Turner, B. L. (2011). Comment on “The response of vegetation on the Andean flank in Western Amazonia to Pleistocene climate change”. Science 333:1825. doi: 10.1126/science.1207525
Urrego, D. H., Bush, M. B., Silman, M. R., Correa-Metrio, A., Ledru, M. P., Mayle, F. E., et al. (2009). “Millennial-scale ecological changes in tropical South America since the Last Glacial Maximum,” in Past Climate Variability from the Last Glacial Maximum to the Holocene in South America and Surrounding Regions, eds F. Vimieux, F. Sylvestre, and M. Khodri (Paris: Springer).
van’t Veer, R., Islebe, G. A., and Hooghiemstra, H. (2000). Climate change during the Younger Dryas chron in northern South America: a test of the evidence. Quat. Sci. Rev. 19, 1821–1835. doi: 10.1016/S0277-3791(00)00093-7
Whitlock, C., and Larsen, C. (2001). “Charcoal as afire proxy,” in Tracking Environmental Change Using Lake Sediments. Terrestrial, Algal, and Siliceous Indicators, Vol. 3, eds J. P. Smol, H. J. B. Birks, and W. M. Last (Dordrecht: Kluwer), 75–98. doi: 10.1007/0-306-47668-1_5
Whitney, B. S., Mayle, F. E., Punyasena, S. W., Fitzpatrick, K. A., Burn, M. J., Guillen, R., et al. (2011). A 45 kyr palaeoclimate record from the lowland interior of tropical South America. Palaeogeogr. Palaeoclimatol. Palaeoecol. 307, 177–192. doi: 10.1016/j.palaeo.2011.05.012
Keywords: diversity dynamics, eastern Andean flank, Last Glacial Maximum, neotropics, palaeoecology, stability, vulnerability, western equatorial Amazonia
Citation: Montoya E, Keen HF, Luzuriaga CX and Gosling WD (2018) Long-Term Vegetation Dynamics in a Megadiverse Hotspot: The Ice-Age Record of a Pre-montane Forest of Central Ecuador. Front. Plant Sci. 9:196. doi: 10.3389/fpls.2018.00196
Received: 16 November 2017; Accepted: 01 February 2018;
Published: 20 February 2018.
Edited by:Urs Feller, University of Bern, Switzerland
Reviewed by:Irene Tunno, Lawrence Livermore National Laboratory (DOE), United States
Walter Finsinger, Centre National de la Recherche Scientifique (CNRS), France
Juanma Rubiales, Universidad Politécnica de Madrid (UPM), Spain
Copyright © 2018 Montoya, Keen, Luzuriaga and Gosling. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.