SEARCH SEARCH

Article Search

266 tocWhat are the best modern analogs for ancient South American mammal communities? Evidence from ecological diversity analysis (EDA)

Angeline M. Catena and Darin A. Croft

Article number: 23(1):a03
https://doi.org/10.26879/962
Copyright Society for Vertebrate Paleontology, February 2020

Author biographies
Plain-language and multi-lingual abstracts
PDF version

Submission: 27 January 2019. Acceptance: 14 January 2020.

ABSTRACT

Ecological diversity analysis (EDA) is a technique that uses ecological attributes of mammals to reconstruct the community structure and habitat of a fossil locality. EDAs of South American paleofaunas have generally relied on modern comparative datasets from that continent, but modern faunas from other continents may be more appropriate models considering the high-level taxonomic differences that exist between modern and fossil South American mammal communities. To test this hypothesis, we selected five, well-sampled fossil localities for which independent paleoenvironmental data (e.g., paleosols, ichnofossils) have been published: four from South America (La Venta, Colombia; Quebrada Honda, Bolivia; Santa Cruz, Argentina; Tinguiririca, Chile) and one from Europe (Rümikon, Switzerland). We coded the extinct species from these sites, as well as ca. 2,450 modern mammal species, for three ecological attributes: diet (eight categories), locomotor habit (six categories), and body mass (six categories). Percentages of species in each attribute category were used to compare the five paleofaunas to 179 modern faunas from six continents using correspondence analysis, hierarchical clustering, similarity percentage, and classification trees. The four South American paleofaunas were found to be most similar to Afrotropical, Indo-Malayan, and Palearctic modern faunas and, similarly, the European paleofauna most resembled faunas from a biome not currently present in the Palearctic. Our study highlights important differences in community structure between ancient and modern South American mammal faunas and suggests that modern mammalian communities from other continents are better analogues for ancient South American communities than modern South American faunas.

Angeline M. Catena. Diablo Valley College, 321 Golf Club Road, Pleasant Hill, California, 94523-1529, USA. Angeline.Catena@case.edu
Darin A. Croft. Department of Anatomy, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio, 44196-4930, USA. Darin.Croft@case.edu

Keywords: Cenozoic; ecology; community structure; South America; paleoenvironment; neotropics

Final citation: Catena, Angeline M. and Croft, Darin A. 2020. What are the best modern analogs for ancient South American mammal communities? Evidence from ecological diversity analysis (EDA). Palaeontologia Electronica, 23(1):a03. https://doi.org/10.26879/962
palaeo-electronica.org/content/2020/2909-south-america-mammal-eda

Copyright: February 2020 Society of Vertebrate Paleontology.
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
creativecommons.org/licenses/by/4.0/creativecommons.org/licenses/by-nc-sa/4.0/

INTRODUCTION

Paleontologists use a variety of approaches to reconstruct terrestrial habitats, climates, and ecological communities (Dark, 2008; Piperno, 2008; Shuman, 2009; Wright, 2017; Croft et al., 2018b). For example, sedimentological analyses can be used to determine the sources of sediments and characteristics of a depositional environment (Boggs, 2006; Nichols, 2009, Wright, 2017). Paleopedology can be used to deduce environmental factors such as mean annual precipitation and paleotopographic position (Caudill et al., 1996; Cerling et al., 1989; Retallack, 1994; Kraus, 1999; Retallack, 2001; Driese and Ober, 2005; Retallack, 2005; Kraus and Riggins, 2007), and ichnology can be used to determine soil moisture regimes, primary productivity, and shifts in ecosystems related to changes in environmental conditions (Rhoads, 1975; Retallack, 1984; Hasiotis, 2002, 2007; Kraus and Riggins, 2007; Smith et al., 2008; Hembree 2018). Plant fossils, including macrofossils (wood, leaves), pollen, and phytoliths, can provide taxonomic and ecomorphological information pertaining to floral communities, vegetative cover, and climate (Hooghiemstra, 1989; Wing and Greenwood, 1994; Gregory-Wodzicki et al., 1998; Dunn et al., 2015; Wright 2017; Kovarovic et al., 2018; Mander and Punyasena, 2018; Peppe et al., 2018; Strömberg et al., 2018). Ecological attributes and/or environmental constraints of vertebrates can be used to infer properties of the environments in which they lived (Fleming, 1973; Andrews et al., 1979; Eisenberg, 1981; Brown, 1995; Blain et al., 2008; Finlayson et al., 2011; Louchart et al., 2001; Cruz et al., 2016; Vermillion et al., 2018).

Under ideal circumstances, all of the aforementioned techniques can be used to create a multifaceted paleoenvironmental reconstruction of a fossil locality (Su and Croft, 2018). Unfortunately, this is seldom possible. Sedimentology has only existed as a distinct branch of geology for several decades (Nichols, 2009); consequently, the sedimentology (and paleopedology) of many historically older fossil localities remain unknown or are described in insufficient detail for stand-alone facies and environmental analyses (Bronger and Catt, 1989; Kraus and Bown, 1986; Retallack, 2001; Stanley, 2005; Matheos and Raigemborn, 2012). Continental ichnology is also a relatively new field, and ichnofossils are often overlooked (Voorhies, 1975; Howard, 1977; Buatois and Mángano, 1995; Hembree et al., 2012). Plant macrofossils and pollen are not typically preserved in the same sediments that preserve terrestrial vertebrates, and although this is not true of phytoliths, their abundance can vary considerably depending on sedimentological conditions (Solomon and Silkworth, 1986; Hjelmroos and Franzén, 1994; Janis, 2000; Wright, 2017). In this context, ecological diversity analysis (EDA) can be an effective alternative for reconstructing terrestrial paleoenvironments, particularly for faunas that are already well sampled and described. EDA has the added benefit of providing insights into the ecological structure of past communities.

Ecological diversity analysis is a technique based on the concept that the ecological adaptations of an animal are reflected in its morphology (Losos, 1990; Saunders and Barclay, 1992; Wainwright and Reilly, 1994; Werdelin, 1996). Studies of modern mammal communities have shown that proportions of mammals in different ecomorphological categories vary with habitat (Fleming, 1973; Andrews et al., 1979; Eisenberg, 1981; Brown, 1995). Although the biological attributes of a fossil species cannot be observed directly, such attributes can be inferred with some confidence by comparing craniodental and postcranial anatomy to that of modern species (Kappelmann, 1988; Scott et al., 1990; DeGusta and Vrba, 2003, 2005; Kovarovic et al., 2018). By repeating this process for all species in a fossil assemblage and comparing the resulting community structure to that of modern faunas, the paleohabitat of a fossil site can be reconstructed (Andrews and Van Couvering, 1975; Nesbit Evans et al., 1979; Van Couvering 1980). This type of ecomorphological approach is considered to be taxon-independent because the proportions of species occupying broad ecological niches within a particular habitat type tend to be similar regardless of the taxonomic affinities of the species involved (Fleming, 1973; Andrews and Nesbet Evans, 1979; Andrews et al., 1979 Damuth, 1992). In theory, an EDA of any fauna for which the ecological characteristics of its species can be determined should faithfully reflect the habitat in which those species were living.

Ecological diversity analysis was originally developed and applied to late Cenozoic hominin-bearing sites in Africa (e.g., Andrews and Lord, 1979; Andrews et al., 1979; Van Couvering 1980), but the technique has also been used to reconstruct South American paleohabitats (e.g., Kay and Madden, 1997a, b; Flynn et al., 2003; Croft et al., 2008; Kay et al., 2012; Lintulaakso and Kovarovic, 2016). The study of these geologically older South American mammal faunas comes with a formidable set of challenges due to the continent’s biogeographic history. For the majority of the Cenozoic, South America was geographically and therefore biotically isolated from other landmasses and developed a highly endemic terrestrial mammal fauna (Patterson and Pascual, 1968; Fooden, 1972; Simpson, 1980; Marshall et al., 1983; Flynn and Wyss, 1998; Croft, 2016). During the Paleocene and most of the Eocene, South American mammal faunas consisted mainly of metatherians (opossums and many extinct relatives), several groups of endemic “ungulates” (notoungulates, astrapotheres, pyrotheres, and litopterns), and xenarthrans (relatives of sloths, armadillos, and anteaters) (Patterson and Pascual, 1968; Simpson, 1980; Marshall et al., 1983; Flynn and Wyss, 1998). Caviomorph rodents and platyrrhine primates reached South America during the middle to late Eocene via waif dispersal from Africa (Antoine et al., 2012; Bond et al., 2015) and subsequently became major components of most ecosystems (Flynn and Wyss, 1998; Poux et al., 2006; Vuectich et al., 2010). Until the late Cenozoic, when North American mammals began to appear in South America (Webb, 1976; Woodburne, 2010; Cione et al., 2015), nearly all terrestrial, nonvolant mammals in South America pertained exclusively to these groups. During the Great American Biotic Interchange (GABI; Stehli and Webb, 1985) a protracted interval of faunal exchange between North and South America (Stehli and Webb, 1985; Webb, 2006; Woodburne et al., 2006; Woodburne, 2010), six orders and 18 families of North American mammals emigrated to South America including camelids (camels), cervids (deer), felids (cats), canids (dogs), and cricetid rodents (New World mice and rats) (for details, see Webb, 2006, Woodburne, 2010, Cione et al., 2015). A Pleistocene megafaunal extinction eliminated ~83% of South American megafaunal (> 44 kg) genera (Koch and Barnoksy, 2006), including the last representatives of two orders of endemic ungulates (Notoungulata, Litopterna) (Cifelli, 1985; Marshall and Cifelli, 1990; Croft, 2000; Bond et al., 2001; Koch and Barnosky, 2006; Fiedel, 2009; Barnosky and Lindsey, 2010; Villavivencio et al., 2015). These two events substantially changed the high-level taxonomic composition of South American mammal communities relative to modern ones, with the result that modern South American mammal faunas bear very little resemblance to those that predominated for much of the Cenozoic (Croft, 2012).

Given the aforementioned differences between fossil and modern South American mammal faunas, one might reasonably question whether modern South American mammal faunas are the best ecological analogues for South American paleofaunas (cf. Croft and Townsend, 2005). Even though general resemblances in community structure may be sufficient for reconstructing paleohabitat (e.g., Kay and Madden, 1997a, b; Flynn et al., 2003; Croft et al., 2008), limiting the modern comparative dataset to the same biogeographic realm may preclude a more nuanced interpretation of community structure (Lintulaakso and Kovarovic, 2016).

The goal of this present study is to test if modern faunas from outside of South America are better analogues for South American fossil mammal communities than modern South American faunas. This proposition is evaluated via EDA using a global comparative dataset and five well-sampled paleofaunas, four from South America and one from Europe. These paleofaunas were chosen because they represent different geographic areas and depositional environments and have been the subject of paleoenvironmental reconstruction based on other types of evidence. We make three predictions: Prediction 1: our paleoenvironmental reconstructions based this global comparative dataset will be congruent with reconstructions based on other types of data; this is based on the previously-established efficacy of EDAs for reconstructing paleoenvironments. Prediction 2: the community structures of the South American paleofaunas will most closely resemble those of modern faunas from continents other than South America; this is based on observed taxonomic and body size differences between fossil and modern South American faunas. Prediction 3: any sampling biases (due to preservation or collection methods) shared among all five paleofaunas will not override underlying ecological signals (i.e., all five faunas will not plot together in multivariate analyses); this is based on the assumption that sampling biases should not be significant in these faunas, which were selected a priori to minimize such factors.

FOSSIL LOCALITIES

Localities Rationale

The paleofaunas chosen for this EDA (described below) derive from a small number of South American localities that have relatively precise stratigraphic control and an extensive (and up-to-date) faunal list. Within this limited subset of well-sampled localities, they were chosen because they represent distinct geographic areas (nearly the entire latitudinal range of South America), temporal intervals (early Oligocene to middle Miocene), paleohabitats (relatively arid grassland/shrubland to moist broadleaf forest), and depositional environments (floodplains, rivers, lakes, lahar deposits) (Croft, 2016). Thus, any similarities among these four sites can only reasonably be due to two factors: (1) taxonomic characteristics shared by South American paleofaunas; or (2) systematic sampling biases common to fossil mammal faunas in general. The European locality of Rümikon was included to distinguish between these two possibilities. Since the mammals of Rümikon are taxonomically distinct from those of the South American localities (no families and only two orders of mammals are shared between them), any features common to all five fossil localities must reflect some other sampling bias, such as time-averaging or under-sampling species with particular ecological characteristics.

La Venta, Colombia

s figure1La Venta (ca. 3° 20’ N, 75° 20’ W) is located in the Magdalena Valley of central Colombia (Figure 1) (Guerrero, 1997); it is composed of sedimentary units within the Honda Group, and the fossil bearing units date to the late middle Miocene (13.5-11.5 Ma, Serravallian Age; Flynn et al., 1997). The Monkey Beds is an approximately 15 m-thick, extremely fossiliferous interval within the Honda Group that has been constrained to an approximately 15,000-year interval within Polarity Chron C5AAn (Kay and Madden, 1997a), which spans 13.18-13.03 Ma based on Ogg (2012). La Venta is the richest pre-Pleistocene fossil locality within the tropical latitudes of South America (Croft, 2016); 66 species of non-volant mammals have been identified from La Venta in addition to a diversity of non-mammal vertebrates, and 50 mammal species are preserved within the Monkey Beds (Appendix 1). Based on paleobotany, paleopedology, and a previous EDA, the paleoenvironment of La Venta has been reconstructed as a river-associated, multi-storied tropical forest environment with a discontinuous canopy (Kay and Madden, 1997a). La Venta’s mean annual precipitation (MAP) is conservatively estimated at 1,500-2,000 mm (Kay and Madden, 1997a, b; Spradley et al., 2019).

Quebrada Honda, Bolivia

Quebrada Honda (ca. 21° 57’ S, 65° 25’ W) is located in southern Bolivia, approximately 65 km southwest of the city of Tarija (Figure 1), and its fossils derive from sediments of the Honda Group of Bolivia (unrelated to the Honda Group of Colombia) (MacFadden and Wolff, 1981). Quebrada Honda is coeval with La Venta; the entire sedimentary section of Quebrada Honda has been dated to the late middle Miocene (13.0-11.9 Ma; Serravallian Age), and the lower, most fossiliferous zones also fall within Polarity Chron C5AAn (MacFadden et al., 1990). At present, 33 species of non-volant mammals have been identified from Quebrada Honda (Appendix 1). Based on paleopedology and ichnology, the paleoenvironment of the most fossiliferous intervals of Quebrada Honda have recently been reconstructed as representing a seasonal, sub-humid to semi-arid wooded grassland or savanna-like environment with a MAP of ca. 1,000 mm (Catena et al., 2017).

Santa Cruz, Argentina

The localities of the coastal portion of the Santa Cruz Formation (ca. 50° 17' W, 51° 36' S) (Figure 1) are located in the southern Argentinian province of the same name and sample the late early Miocene (~17.9-16.2 Ma; Burdigalian Age; Perkins et al., 2012; Vizcaíno et al., 2012). Fossiliferous Levels 1-7 are fossil-rich stratigraphic intervals that crop out in four discrete but closely-spaced localities (Vizcaíno et al., 2012) and have been dated to a narrow time interval that spans 17.5-17.4 Ma (Fleagle et al., 2012; Perkins et al., 2012). Santa Cruz is one of the richest fossil localities in South America; over 100 mammal species have been reported (Croft, 2016; but see also Croft, 2013), and 48 species of non-volant mammals have been identified from Fossiliferous Levels 1-7 (Appendix 1). Based on sedimentology, ichnology, paleobotany, paleopedology, and an EDA, the paleoenvironment of Santa Cruz has been reconstructed as a mosaic of open temperate humid to semi-arid forest environments with pronounced seasonality and a MAP of approximately 1,000-1,500 mm (Kay et al., 2012; Raigeborna et al., 2019; Spradley et al., 2019).

Tinguiririca, Chile

Tinguiririca is located in central Chile near the town of Termas del Flaco (34° 59’ S, 70° 26’ W; Wyss et al., 1990; Flynn et al., 2003) (Figure 1). The fossiliferous horizons are composed of volcaniclastic sediments of the Abanico Formation (Charrier et al., 1996; Flynn et al., 2003) and have been dated to the early Oligocene (~31.5 Ma; Rupelian Age; Flynn et al., 2003, 2012). Twenty-six species of mammals have been identified from the site (Appendix 1). Faunal-based techniques, such as cenograms and EDAs, are presently the only options for paleohabitat reconstruction at Tinguiririca, as other biotic and sedimentary features, like phytoliths and paleosols, are not preserved (Flynn et al., 2003; Su and Croft, 2018). Tinguiririca has been interpreted as an open-habitat environment (Wyss et al., 1994; Flynn et al., 2003; Croft et al., 2008) with a MAP of approximately 600 mm (Flynn et al., 2003).

Rümikon, Switzerland

Rümikon is located in northern Switzerland near its border with Austria, approximately 25 km northeast of Zurich (47° 29’ N, 8° 47’ E; NOW Database) (Figure 1). The sediments of Rümikon are part of the Upper Freshwater Molasse (Kälin and Kempf, 2009), which is the uppermost formation of the Molasse Group sediments within the North Alpine Foreland Basin of Switzerland (Matter et al., 1980; Labhart, 2005). The sediments of Rümikon are placed in an early middle Miocene interval from 14.6-14.2 Ma (Langhian Age) based faunal correlations with Swiss and Bavarian localities within the Upper Freshwater Molasse that have been chronostratigraphically dated (Kälin and Engesser, 2001; Kälin and Kempf, 2009). Thirty-seven species of non-volant mammals have been identified from Rümikon (Appendix 1; Kälin and Kempf, 2009; NOW Database, 2017). Based on paleobotanical data and the paleobiology of fossil rodents, the habitats of MN 6 localities in the Upper Freshwater Molasse have been reconstructed as humid to sub-humid mixed and forested environments (Hantke, 1984; Kälin and Kempf, 2009). A general MAP range of approximately 1,300-1,500 mm has been estimated for the fossil localities of the Upper Freshwater Molasse (Hantke, 1984; Kälin and Kempf, 2009).

METHODOLOGY

Faunal Lists

The faunal lists for La Venta, Quebrada Honda and Tinguiririca are from Croft (2016), updated based on Engelman et al. (2017, 2018), Brandoni et al. (2018), McGrath et al. (2018), and Wyss et al. (2018). The faunal list for Santa Cruz is from Kay et al. (2012), and the faunal list for Rümikon is from the NOW database (http://www.helsinki.fi/science/now/), which was downloaded on November 01, 2017.

Modern faunal lists from 179 ecoregions are taken directly from the World Wild Life Fund’s WildFinder Database (WFDB) (World Wildlife Fund 2006; full database can be downloaded at http:// www.worldwildlife.org/science/data.cfm). The WFDB contains presence/absence data for terrestrial vertebrates compiled from field guides and scientific works; the data are based on the ranges of extant species and do not include species that are present as human commensals or were introduced by humans (World Wildlife Fund, 2006). In order to reflect the historic, pre-disturbed, or potential vegetation of each ecoregion, the WFDB uses historic ranges of species with approximate distributions at 1,500 A.D. instead of current distributions, when possible (World Wildlife Fund, 2006).

Biogeographic realms are the largest biogeographic division of Earth; they share and are broadly characterized by the evolutionary history of the flora and fauna contained within them (Reid et al., 2005). The WFDB divides Earth into eight terrestrial biogeographic realms based on Olson et al. (2001), six of which were used in this study: Australasia (Appendix 2); Afrotropic (Appendix 3); Indo-Malaya (Appendix 4); Nearctic (Appendix 5); Neotropic (Appendix 6), which includes Central America and the tropical area of North America in addition to all of South America; and the Palearctic (Appendix 7, Figure 1). Oceana and Antarctic biogeographic realms were excluded due to their lack of terrestrial (non-marine) mammals. The WFDB (i.e., Olson et al., 2001) also recognizes 14 terrestrial biomes, which are floral and faunal communities characterized by adaptations to a particular environment (Campbell, 1996). Each biogeographic realm includes one or more of these biomes, and 11 of these biomes were included in this study. The polar boreal forest/taiga and tundra biomes (Biomes 6 and 11, respectively) were excluded, as the paleofaunas selected for this analysis are from tropical and temperate, rather than polar regions; the mangroves biome (Biome 14) was also excluded, as those ecoregions are poorly mapped and lack adequate species data (World Wildlife Fund, 2006). Each biome in a biogeographic realm is subdivided into one or more ecoregions, which are relatively large areas (mean = 150,000 km2, median = 56,300 km2) that represent distinct biotas (Dasmann, 1973, 1974; Udvardy, 1975; Olson and Dinerman, 1998; Olson et al., 2001) and accurately describe terrestrial biodiversity patterns of vertebrates and plants (Smith et al., 2018). Each ecoregion is unique and identified by a six-digit alphanumeric code. As an example of this hierarchical structure, AA0105 represents the Central Range Montane Rain Forests ecoregion, which is a particular area of Biome 1 (tropical and subtropical moist broadleaf forests; middle pair of digits) within the Australasia biogeographic realm (prefix AA).

The ecoregions used in this study were selected to encompass maximum habitat variation while maintaining relatively equal representation among biomes and biogeographic realms. This was achieved by selecting up to five ecoregions for each biome within each biogeographic realm. For example, for the Indo-Malaya biogeographic realm, five ecoregions were selected within both Biome 1 and Biome 2. All ecoregions were included if five or fewer were represented in a particular biome in a particular biogeographic realm; if more than five ecoregions were available, five were chosen at random using a random number generator. Since precipitation is a major determinant of vegetational structure (Nicholson et al., 1990; Kay and Madden, 1997b), a dataset containing a 30-year average of mean annual precipitation (MAP) with spatial resolution of approximately 1 km2 was obtained from WorldClim Version 2 for ArcGIS (Fick and Hijmans, 2017) for the 179 ecoregions used in this study (Table 1; Appendix 2Appendix 3Appendix 4Appendix 5Appendix 6, Appendix 7). A linear regression with these data was used to estimate MAP for fossil localities (described below).

To avoid using taxonomically under-sampled ecoregions with inaccurate faunal lists, only ecoregions with at least 20 species of non-volant mammals were included. This is an arbitrary threshold based on the observation that most terrestrial mammal faunas include at least 30 species (incudling a variable proportion of volant mammals; Ceballos and Ehrlich, 2006; see also Mares and Willig, 1994), and that most well-sampled South American fossil localities include at least 20 species (Croft, 2006). As with other ecological diversity analyses (e.g., Andrews et al., 1979; Reed, 1998; Kay and Madden 1997a, b; Kovarovic et al., 2002; Mendoza et al., 2004; Louys et al., 2011; Lintulaakso and Kovarovic, 2016), volant mammals (i.e., bats, order Chiroptera) and fully aquatic mammals (i.e., manatees/dugongs and whales, order Sirenia and infraoder Cetacea, respectively) were excluded from our analyses.

Ecological Variables

All extant and extinct mammal species were coded for three ecological attributes - body size, primary locomotor style, and diet - based primarily on the published literature (Kay and Madden, 1997a, b; Kay et al., 1997; Nowak, 1999; Flynn et al., 2003; MacDonald, 2006; Vizcaíno et al., 2006, 2012; Croft et al., 2008; Kay et al., 2012; Toledo et al., 2014; Bradham et al., 2015; Gainsbury et al., 2018; Appendix 2Appendix 3Appendix 4Appendix 5Appendix 6, Appendix 7, and Appendix 8; Table 2). Data from congeners were used in cases were such ecological data were unavailable.

Species were assigned to one of seven logarithmic body mass categories ranging from 0-10 g to 1,000-10,000+ kg (Table 2). The distribution of mammalian body sizes is right skewed, with smaller taxa (≤ 10 kg) naturally occurring at higher frequencies than larger taxa (Caughley 1987; Gardezi and da Silva, 1998; Kozlowski and Gawelczyk, 2002). A logarithmic scale was used to account for this greater frequency of small taxa and to more accurately characterize body mass ranges of the fauna.

The primary locomotor style ecological attribute included six categories: Arboreal, Scansorial, Terrestrial, Semi-fossorial, Fossorial, and Semi-aquatic (Table 2). These categories parallel those used by Kay and Madden (1997a, b) and Kay et al. (2012), except that we did not discriminate among terrestrial, cursorial, and ambulatory species (all were classified as Terrestrial).

The diet ecological attribute included eight categories (Table 2). Mammals display a wide range of dietary preferences that can vary seasonally and with resource availability (Jarman, 1971; Fedriani et al., 2001; Sinclair and Zeppelin, 2002). To more accurately characterize this dietary variation, we used a new method to code for diet in which each species was coded for two equally-weighted diet categories. Species with a relatively homogeneous diet (e.g., giant anteater, Myrmecophaga tridactyla), were coded twice for the same diet category (Diet 1 = In; Diet 2 = In), whereas species with a heterogeneous diet (e.g., monito del monte, Dromiciops gliroides), were coded for two different diet categories (Diet 1 = Fr; Diet 2 = In). Species with a diet spanning three or more diet categories were double-coded as omnivorous (e.g., common opossum, Didelphis marsupialis, Diet 1 = Om; Diet 2 = Om), and species with a single primary diet but a secondary diet spanning three or more categories (e.g., nine-banded armadillo, Dasypus novemcinctus), were coded as omnivorous for one of the diet categories (Diet 1 = In; Diet 2 = Om). Because these categories were weighted equally, designating the primary diet as “Diet 1” or “Diet 2” had no effect on analyses. Since each species was coded twice, the denominator for computing the percentage of species in each dietary category was twice the number of species in each ecoregion (2N). The resulting percentages of each diet category were used in the subsequent analyses.

Data Analyses

The variables used in the analyses were percentages of species in each ecological attribute category within each fauna (Appendix 8, Appendix 9). These 21 variables (one for each ecological attribute category) were analyzed in PAST v. 3.1 (Hammer et al., 2001) using correspondence analysis (CA), linear regression, hierarchical cluster analysis (HCA), and similarity percentage (SIMPER) analysis. These variables were also analyzed via classification trees (CTs) in JMP Pro 13 (SAS Institute, Cary, NC). All analyses were performed with and without the Size 0 variable to determine whether a lack of very small species in a fossil fauna (as commonly occurs in the absence of screen-washing) might bias EDA-based paleohabitat interpretations. The analyses used in this study do not assume multivariate normal distributions of the data (Hammer and Harper, 2006). Nonetheless, all analyses were also performed with and without an arcsine data transformation. There were no discernable differences in the results between the transformed and non-transformed data; the results of the analyses using the non-transformed data are presented here.

Correspondence analysis is an exploratory technique used to visualize trends and groupings in a multivariate dataset (Hammer and Harper, 2006; LeGendre and LeGendre, 2012). In our study, the CA plotted faunas and ecological variables in the same space and maintained correspondence between the two, permitting ecological variables most highly associated with particular faunas to be determined. The CA results were also used with linear regression analysis to infer paleoprecipitation; MAP data for modern faunas was regressed on first axis score (CA1 scores), and paleoprecipitation was estimated based on the CA1 values of the fossil faunas. We used these results to test all three of our predictions. The CA and the associated linear regression were used to broadly compare the results of this study to those of the previous studies and were also used to compare the paleofaunas to the modern faunas grouped by biome and biogeographic realm. The positions of the paleofaunas on the CA axes relative to each other were used to test Prediction 3.

Hierarchical cluster analysis is an exploratory technique used to identify trends and groupings in a multivariate data set based on a similarity index (Hammer and Harper, 2006; Legendre and Legendre, 2012). HCA determines which faunas are most similar to one another in terms of their ecological variables (Hammer and Harper, 2006). Unweighted pair group method with arithmetic mean (UPGMA) clustering was used in the HCAs. This is an agglomerative clustering method in which clusters are based on the average distance between all members of the group; it is commonly used in data sets with ecological variables (Hammer and Harper, 2006; Legendre and Legendre, 2012). The chord similarity measure was used in the HCA; it is used to compare faunal assemblages by species abundances (Legendre and Legendre, 2012) and has greater accuracy in high-dimensional datasets than other similarity indexes, such as Euclidian distance and Average distance (Shirkhorshidi et al., 2015). Six HCAs, one including all five paleofaunas and five with one paleofauna each, were performed in order to visualize the relationships of the fossil and modern faunas as a group and individually. The HCAs with individual paleofaunas were used to test which biomes and specific ecoregions the fossil faunas most closely resembled; they were also used to further assess Prediction 2.

A SIMPER analysis quantifies and ranks the average percentage contribution of different variables to the dissimilarity between samples (Clarke, 1993). SIMPER analyses were used to assess which ecological variables were primarily responsible for similarities and differences between the individual paleofaunas and the groups of modern faunas within each of the 11 biomes. The chord similarity measure was also used in the SIMPER analyses. Dissimilarities between the paleofaunas and the modern faunas (grouped by biomes) were used to test Prediction 3.

Classification trees belong to the more general category of statistical analyses known as classification and regression trees (Breiman et al., 1984), which are used to construct predictive models for complex multivariate systems (De’ath and Fabricus, 2000). We used three CTs to create models and further explore the relationships between the ecological variables and the habitats represented by the modern faunas. The first CT model (CT1) was built to classify modern faunas by vegetative cover; ecoregions of Biomes 1, 2, 3, 4, and 5 were coded as “closed”, those of Biomes 7, 8, 9, 10, and 12 were coded as “mixed”, and those of Biome 13 were coded as “open”. The second CT model (CT2) was built to classify modern faunas by biogeographic realm (six categories), and the third CT model (CT3) was built to classify the modern faunas by biome (11 categories). The CT models were then used to predict the vegetative cover, biogeographic realm, and biome of the five paleofaunas. As five is the maximum number of selected ecoregions from within the same biome and biogeographic realm, terminal leaves with fewer than five faunas were pruned in order to avoid overfitting the CT models. The predictive models as well as the relationships between specific ecological variables and the modern faunas at the biome and biogeographic level were used to further test all three of our predictions.

RESULTS

Correspondence Analysis

s figure2Four axes of the CA collectively account for two-thirds of the total variation within the dataset (Table 3) and are discussed herein. Axis one (CA1) has an eigenvalue of 0.065 (27.4% of the total variation) and generally represents a continuum from closed to open habitats (Table 3; Figure 2A). The Arboreal, Frugivorous, Size 3, Scansorial, and Size 2 variables have the greatest negative correlation on CA1 (Figure 2B); higher proportions of species in these ecological attribute categories are found in in mammalian communities from closed habitats (Kay et al., 1997b). The Fossorial, Semi-fossorial, Size 0, Graminivorous, and Granivorous variables have greatest positive correlation on CA1 (Figure 2B); higher proportions of species in four of these ecological attribute categories (Fossorial, Semi-fossorial, Graminivorous, and Granivorous) are found in open habitats (Guo et al., 1997; Kinlaw, 1999; Gibson, 2009; Davidson et al., 2012). The Australasia and Afrotropic biogeographic realms have the broadest ranges on CA1, whereas the Nearctic and Palearctic biogeographic realms have the narrowest ranges; all of the Nearctic faunas and all but three of the Palearctic faunas have positive CA1 scores. Faunas from Biomes 1, 2, 3, and 4 mostly have negative scores on CA1; those from Biomes 5, 7, 9, and 10 largely plot near the origin on CA1; and those from Biomes 8, 12, and 13 have negative scores on CA1. Among paleofaunas, Quebrada Honda has a positive score on CA1, while La Venta, Santa Cruz, Tinguiririca, and Rümikon have negative scores on CA1.

Axis two (CA2) has an eigenvalue of 0.042, accounts for 17.5% of the total variation, and generally represents a continuum of faunas with many to relatively few large species (Table 3; Figure 2A). The Size 6, Size 5, Size 4, Semi-aquatic, and Carnivorous variables have the greatest negative correlation on CA2, and the Size 0, Semi-fossorial, Arboreal, Size 1, and Fossorial variables have the greatest positive correlation on CA2 (Figure 2B). The Afrotropic and Indo-Malaya biogeographic realms have the broadest ranges on CA2, whereas the Neotropic and Nearctic biogeographic realms have the narrowest ranges on CA2. All Australasian faunas and all but one Nearctic fauna have positive scores on CA2; all but six faunas from the Afrotropics and Indo-Malaya have negative scores on CA2. Faunas from Biomes 1 and 3 largely have positive scores on CA2; those from Biomes 2, 4, 7, 8, 12, and 13 largely plot toward the origin of CA2; and those from Biome 9 have negative scores on CA2. The four South American paleofaunas have negative scores on CA2, while Rümikon has a positive score on CA2.

When both CA1 and CA2 are considered together, the fossil faunas plot as follows: La Venta falls outside the range of all modern biogeographic realms and biomes; Quebrada Honda and Santa Cruz plot within the range of the Afrotropic biogeographic realm and Biome 9 (flooded grasslands and savannas); Tinguiririca plots within the ranges of the Afrotropic and Indo-Malaya biogeographic realms and Biome 7 (tropical and subtropical grasslands, savannas, and shrublands); Rümikon plots within the ranges of Biomes 1, 2, 3, 4, 7, and 10 of the Australasia, Afrotropic, Neotropic, and Palearctic biogeographic realms (Figure 2A).

s figure3Axis three (CA3) has an eigenvalue of 0.029 and accounts for 12.0% of the total variation (Table 3, Figure 3A). The Size 6, Insectivorous, Graminivorous, Size 1, and Ramivorous variables have the greatest correlation on the negative side of CA3; the Granivorous, Semi-fossorial, Carnivorous, Omnivorous and Lignivorous variables have the greatest correlation on positive side of CA3 (Figure 3B). CA3 generally discriminates among biogeographic realms. All of the Nearctic and all but one of the Indo-Malayan faunas plot on the positive side of CA3; all of the Afrotropical and all but three of the Australasian faunas plot on the negative side of CA3. The Australasia and Palearctic biogeographic realms have the broadest ranges on CA3, and the Afrotropic, Neotropic, and Indo-Malaya biogeographic realms have the narrowest ranges on CA3. La Venta, Quebrada Honda, Rümikon, and Tinguiririca have positive scores of CA3, and Santa Cruz has a negative score on CA3.

Axis four (CA4) has an eigenvalue of 0.023 and accounts for 9.6% of the total variation (Table 3). No major trends are discernable among the biomes and biogeographic realms on CA4. The Ramivorous, Semi-fossorial, Graminivorous, Size 2, and Size 4 variables have the greatest correlation on the negative side of CA4; the Size 0, Semi-aquatic, Scansorial, Lignivorous, and Size 6 variables have the greatest correlation on the positive side of CA4; (Figure 3B). La Venta, Quebrada Honda, Rümikon, and Tinguiririca have positive scores on CA4, while Santa Cruz has a negative score. The Palearctic and Neotropic biogeographic realms have the broadest ranges on CA4, and Afrotropic and Nearctic biogeographic realms have the narrowest range on CA4.

On CA3 and CA4, Quebrada Honda, Santa Cruz, and Tinguiririca do not plot within the range of any biogeographic realm or biome. La Venta plots within the ranges of the Australasia biogeographic realm and Biome 10 (montane grasslands and shrublands), while Rümikon plots within the ranges of the Afrotropic, Neotropic, and Palearctic biogeographic realms and all biomes (Figure 3A).

s figure4Linear Regression

An ordinary least squares regression of MAP on CA1 score results in a predictive equation with a negative slope and an r2 value of 0.61 (Figure 4). Based on the CA1 scores of the paleofaunas and the 95% confidence interval of the regression equation, MAP is estimated at 1,917-2,035 mm for La Venta, 1,210-1,286 mm for Santa Cruz, 1,115-1,214 mm for Rümikon, 1,079-1,186 mm for Tinguiririca, and 877-1,032 mm for Quebrada Honda. Regressions of MAP on CA2, CA3, and CA4 resulted in r2 values of 0.01, 0.01, and 0.04, respectively.

Hierarchical Cluster Analysis

s figure5When all five paleofaunas are included in the HCA, the four South American paleofaunas form an exclusive cluster that is the second most basal cluster of the dendrogram (Figure 5); Rümikon is the basal branch within a large cluster of 51 ecoregions from 10 biomes (mostly Biomes 10 and 12, represented by nine ecoregions each), from the Afrotropic (n=20), Palearctic (n=19), Indo-Malaya (n=11), and Nearctic (n=1) biogeographic realms. When analyzed individually using HCA (i.e., without other fossil localities), the position of Rümikon is the same. In contrast, the positions of the four South American paleofaunas are distinct. La Venta plots among modern faunas, as the basalmost branch in a cluster of seven Australasian ecoregions and three Indo-Malayan ecoregions; five of these ecoregions represent Biome 1, one represents Biome 3, and one represents Biome 10. Quebrada Honda plots as the basalmost branch among a small cluster of four Afrotropical ecoregions from Biome 9. Santa Cruz does not cluster with any modern faunas due to high proportions of species in the Size 4, Scansorial, and Size 5 ecological variables. Tinguiririca also does not cluster with any modern faunas, the result of high proportions of species in the Terrestrial, Ramivorous, Graminivorous, and Size 3 variables.

Similarity Percentage

La Venta is most similar to Biomes 2, 1, and 7 (dissimilarity = 0.26, 0.29, 0.29, respectively; Appendix 9) and most dissimilar to Biomes 13, 12, and 8 (dissimilarity = 0.51, 0.48, 0.45, respectively). Collectively, the Size 1, Ramivorous, Arboreal, and Semi-fossorial variables account for most of the dissimilarity between La Venta and these latter biomes.

The Quebrada Honda fauna is most similar to Biomes 9, 7, and 2 (dissimilarity = 0.30, 0.31, 0.31, respectively; Appendix 9). The Quebrada Honda fauna is most dissimilar to Biomes 1, 12, and 13 (dissimilarity = 0.47, 0.46, 0.40, respectively). Collectively, the Arboreal, Size 1, Size 4, and Ramivorous variables account for most of the dissimilarity between Quebrada Honda and these latter biomes.

The Santa Cruz fauna is most similar to Biomes 2, 9, and 7 (dissimilarity = 0.36, 0.39, 0.41, respectively; Appendix 9). The Santa Cruz fauna is most dissimilar to Biomes 12, 13, and 1 (dissimilarity = 0.53, 0.52, 0.49, respectively). Collectively, the Ramivorous, Size 1, Size 4, and Arboreal variables account for the most dissimilarity between Santa Cruz and these latter biomes.

The Tinguiririca fauna is most similar to Biomes 7, 2, and 10 (dissimilarity = 0.36, 0.37, 0.38, respectively; Appendix 9). The Tinguiririca fauna is most dissimilar to Biomes 12, 1, and 13 (dissimilarity = 0.52, 0.51, 0.48, respectively). The Ramivorous, Arboreal, and Size 1, and Insectivorous variables account for the most dissimilarity between Tinguiririca and these latter biomes.

The Rümikon fauna is most similar to Biomes 7, 5, and 4 (dissimilarity = 0.15, 0.16, 0.17, respectively; Appendix 9). The Rümikon fauna is most dissimilar to Biomes 1, 9, and 2 (dissimilarity = 0.24, 0.23, 0.21, respectively). Collectively, the Size 3, Arboreal, Carnivorous, and Size 2 variables account for the most dissimilarity between Rümikon and these latter faunas.

Classification Trees

s figure6The CT for vegetative cover (CT1, Figure 6A) has an r2 value of 0.73 and resulted in misclassifications rates of 12%, 17%, and 33% for the ecoregions classified as closed, open, and mixed, respectively (Table 4.1). The Arboreal variable defines the first split in the tree; it separates 17 of the closed, 46 of the mixed, and all 24 of the open ecoregions to the left side of the tree; 57 of the closed and 35 of mixed ecoregions are separated to the right side of the tree. On the left side of the tree, the Semi-aquatic variable defines the second split; it separates four closed, three mixed, and 16 open ecoregions to the left branch and 13 closed, 43 mixed, and eight open ecoregions to the right branch. On the right side of the tree, the Size 2 variable defines the second split; it separates nine closed and 20 mixed ecoregions to the left branch and 48 closed and 15 mixed ecoregions to the right branch. In this model, La Venta is predicted to represent a closed habitat, Quebrada Honda, Santa Cruz, and Rümikon are predicted to represent mixed habitats, and Tinguiririca is predicted to represent an open habitat (Figure 6A).

The CT used to predict biogeographic realms (CT2, Figure 6B) has an r  value of 0.88 (Table 4.2). The Size 6 variable defines the first split, separating 12 of the Indo-Malayan, one of the Palearctic, and all 29 of the Afrotropical ecoregions to the left side of the tree; 13 Indo-Malayan, 25 Palearctic and all Australasian, Nearctic, and Neotropical ecoregions (N = 29, 25, and 35, respectively) are separated to the right side of the tree. On the left side of the tree, the Size 3 variable defines the second split, with 23 of the Afrotropical ecoregions separated to the left branch; 12 Indo-Malayan, six Afrotropical, and one of the Palearctic ecoregions are separated to the right branch. On the right side of the tree, the Size 5 variable defines the second split; 28 Palearctic, 13 Indo-Malayan, and four Nearctic are separated to the left branch, whereas 35 Neotropical, 29 Australasian, 21 Nearctic, and seven Palearctic ecoregions are separated to the right branch. In this model, La Venta, Quebrada Honda, and Tinguiririca are predicted to be within the Indo-Malaya biogeographic realm, and Rümikon and Santa Cruz are predicted to be within the Palearctic biogeographic realm (Figure 6B).

s figure7The CT used to predict biomes (CT3, Figure 7) has an r2 value of 0.64 (Table 4.3). The Fruit variable defines the first split; all ecoregions from Biomes 8, 12, and 13, along with 1-17 ecoregions from all other biomes are separated to the left side of the tree, while 21 ecoregions from Biome 1 and 1-10 ecoregions from Biomes 2, 3, 4, 5, 7, 9, and 10 are separated to the right side of the tree. On the left side of the tree, the Semi-aquatic variable defines the second split; 20 of the remaining ecoregions from Biome 13 and 1-7 ecoregions from all Biomes except Biome 4 are separated to the left branch, while all 17 remaining ecoregions from Biome 4 and 1-12 ecoregions from all Biomes except Biome 1 are separated to the right branch. On the right side of the tree, the Size 0 variable defines the second split. Nine ecoregions from Biomes 1 and 7 and 1-7 ecoregions from Biomes 2, 3, 4, 7, 9 and 10 are separated to the left branch while 12 ecoregions from Biome 1 and 1-6 ecoregions from Biomes 3, 4, 5, 7, and 10 are separated to the right branch. In this model, La Venta is predicted to be within Biome 1, Santa Cruz is predicted to be within Biome 2, Quebrada Honda is predicted to be within Biome 9, Tinguiririca is predicted to be within Biome 10, and Rümikon is predicted to be within Biome 7. When the Size 0 variable is removed from the analysis, Rümikon is predicted to be within Biome 10; all other predictions are unchanged (Figure 7).

PALEOHABITAT RECONSTRUCTIONS

The five types of analyses described above result in paleohabitat interpretations that are broadly consistent with those of previous investigations (Hantke, 1984; Kay and Madden, 1997a, b; Flynn et al., 2003; Croft et al., 2008; Kälin and Kempf, 2009; Kay et al., 2012; Catena et al., 2017). Thus, they confirm our first prediction that these EDA-based paleoenvironmental reconstructions will be congruent with reconstructions based on other types of data. Nevertheless, detailed examination of the mammalian community structures of these sites reveals that they are not analogous to modern fauns from their corresponding biogeographic realm (i.e., Neotropic and Palearctic). Rather, they more closely resemble modern sites from other biogeographic realms. This confirms our second prediction and illustrates the importance of using a global comparative dataset when investigating ancient mammal community paleoecology, especially for fossil sites in biogeographic realms that have experienced significant changes in climate and/or faunal composition between the time of deposition and the Recent.

La Venta, Colombia

The Hierarchical Cluster Analysis (HCA) and Classification Tree (Figure 5, Figure 6B) indicate that the ecological structure of La Venta most closely resembled those of modern mammal communities from the Indo-Malaya biogeographic realm. Even though La Venta plots within the range of Australasian faunas on Axis 3 and 4 of the Correspondence Analysis (Figure 3A), this is primarily due to its high proportion of insectivorous species and low proportion of predatory species, a pattern also noted by Croft (2006) and attributed to a lack of placental carnivores (Croft, 2006; see also Croft et al., 2018a). However, La Venta differs from all Australasian faunas in CT2 (Figure 6B) in the proportion of species in mass categories Size 5 and Size 6; such species are absent in Australasian ecoregions but present in Indo-Malayan ecoregions. Therefore, we consider the paleofauna of La Venta to more closely resemble faunas of Indo-Malaya than those of Australasia in terms of overall ecological structure.

The results of the SIMPER analysis and CT3 (Appendix 9; Figure 7) indicate that the paleofauna of La Venta resembles a modern fauna from Biome 1, moist tropical to subtropical broadleaf forests. This biome is characterized by low temperature variability, high rainfall (MAP generally > 2,000 mm), and semi-evergreen to evergreen deciduous vegetation (Holdridge et al., 1971). Biome 1 ecoregions are mainly present in the modern Indo-Malaya, Neotropic, and Australasia biogeographic realms. Our MAP estimate of 1,917-2,035 mm for La Venta is consistent with those of Biome 1 ecoregions as well as Kay and Madden’s (1997a, b) estimate of 1,500-2,000 mm for the site. The classification tree for vegetative cover (CT1; Figure 6A) indicates that the paleofauna of La Venta resembles that of a modern closed vegetation habitat. Overall, the paleohabitat of La Venta is reconstructed as a humid, closed tropical forest, also consistent with the findings of Kay and Madden (1997a, b).

Quebrada Honda, Bolivia

The collective results of the CA (Figure 2A, Figure 3A) and HCA (Figure 5) indicate that the paleofauna of Quebrada Honda had an ecological structure most similar to that of a modern mammal community from the Afrotropics. While CT2 (Figure 6B) predicts Quebrada Honda to more closely resemble faunas from Indo-Malaya, this is largely due to similar proportions of Size 3 species. Quebrada Honda is distinguished from Indo-Malayan faunas on CA2 (Figure 2A) by its high proportions of Size 4, Terrestrial, and Graminivorous species. In this analysis, faunal communities with similar proportions of species within those ecological categories are only present in modern Afrotropical faunas. We therefore consider the fauna of Quebrada Honda to have an ecological structure most similar to Afrotropical mammal communities.

The results of the HCA, SIMPER, and CT3 (Appendix 9; Figure 5, Figure 7) indicate that the paleofauna of Quebrada Honda had a community structure similar to a modern fauna from Biome 9, flooded grasslands and savannas. This biome is characterized by tropical and subtropical climates, seasonal rainfall generally ranging from 500-1,300 mm/year, vegetation with a non-continuous canopy, and seasonal to year-round flooding (Holdridge et al., 1971; González-Jiménez, 1979). Ecoregions from Biome 9 mainly occur in the modern Afrotropic, Neotropic, and Palearctic biogeographic realms.

The estimated MAP of 877-1,032 mm for Quebrada Honda based on linear regression is within the MAP range of Biome 9 ecoregions used in this study and is also consistent with the previous MAP estimate of ~1,000 mm by Catena et al. (2017). The presence of alluvial systems and periodic of flooding events at Quebrada Honda has previously been indicated by sedimentological, paleosol, and ichnofossil analyses (MacFadden and Wolff, 1981; Croft et al., 2011; Catena et al., 2017) and is compatible with the presence of a freshwater turtle at the site (Acanthochelys; Cadena et al., 2015). However, other aquatic animals such as crocodilians have not yet been discovered at Quebrada Honda. The Classification Tree for vegetative cover (CT1, Figure 6A) indicates that the mammal community of Quebrada Honda had an ecological structure similar to that of a modern habitat with a mixed vegetative cover. Overall, the paleohabitat of Quebrada Honda is reconstructed as a subtropical, river-associated savanna with a relatively open, mixed vegetative cover. This is also consistent with the paleoenvironmental findings of Catena et al. (2017).

Santa Cruz, Argentina

The results of the CT2 (Figure 6B) indicate that the paleofauna of Santa Cruz had an ecological structure similar to a modern mammal community from the Palearctic. While Santa Cruz is within the range of Afrotropical faunas in the CA (Figure 2A, Figure 3A), it differs from Afrotropical ecoregions in its low proportion of species within the Size 6 category. The results of the SIMPER and CT3 (Appendix 9, Figure 7) analyses indicate the fauna of Santa Cruz was ecologically similar to a modern fauna from Biome 2, dry tropical to subtropical broadleaf forests. This biome is characterized by low temperature variability, predominantly deciduous trees, and strongly seasonal rainfall generally ranging from 500-2,000 mm/yr but sometimes greater (Holdridge et al., 1971). Ecoregions from Biome 2 occur mainly in the modern Neotropic and Indo-Malaya biogeographic realms and are absent in the modern Palearctic biogeographic realm.

The estimated MAP of 1,210-1,286 mm for Santa Cruz based on linear regression is within the range of MAPs from ecoregions of Biome 2 used in this study and is consistent with the previous MAP estimate of 1,000-1,500 mm of Kay et al. (2012). CT1 (Figure 6A) indicates that, based on the structure of its fauna, the paleohabitat of Santa Cruz resembled a modern habitat with a mixed vegetative cover. Overall the paleohabitat of Santa Cruz is reconstructed as a subtropical, mixed forest environment, which is consistent with the previous paleoenvironmental findings of Kay et al. (2012) and Raigemborna et al. (2019).

Tinguiririca, Chile

The results of the CA and CT2 (Figure 2A, Figure 3A, Figure 6B) indicate that the paleofauna of Tinguiririca had an ecological structure most similar to that of modern mammal communities from Indo-Malaya. The results of the SIMPER and CT3 (Appendix 9; Figure 7) indicate that the paleofauna had a community structure similar to a modern fauna from Biome 10, the montane grasslands and shrublands. This biome is characterized by montane elevation generally ranging from 1,500-3,500 m, temperate, humid to semi-arid climates, and open, shrubby vegetation (Holdridge et al., 1971; North et al., 2016). Ecoregions from Biome 10 occur mainly in the modern Afrotropic, Neotropic, and Palearctic biogeographic realms; a single ecoregion from Biome 10 is present in the Indo-Malaya biogeographic realm. The estimated MAP of 1,079-1,186 mm for Tinguiririca based on linear regression is within the range of MAPs from ecoregions of Biome 2 used in this study is higher than the previous MAP estimate of ~600 mm by Flynn et al. (2003). The discrepancy between these MAP estimates could be due, in part, to the higher proportion of taxa coded for Ramivory (i.e., Browse) in this analysis compared to Flynn et al. (2003). In our analysis, taxa with high-crowned teeth (notohippids, archaeohyracids, and the interatheriine Santiagorothia) were coded both Ramivorous and Graminivorous (i.e., as both browsers and grazers) instead of just Graminivours (grazers) as in Flynn et al. (2003); this was done to lessen the chance of misinterpreting diet based on tooth crown height (see Croft et al., 2008). Hypsodonty (having high-crowned teeth) in ungulates is considered an adaptation to abrasive diets, such as those with high amounts silica and/or abiotic material that is inadvertently consumed in the form of dust and grit (Janis, 1988; Pérez Barbería and Bordon, 2001; Williams and Kay, 2001; Reguero et al., 2010). In modern faunas, higher relative abundances of browsing taxa (i.e., ungulates) are associated with higher annual rainfall, generally around 1,000 mm annually (Bell, 1971; Cumming, 1982; McNaughton and Georgiadis, 1986). In our analysis, ecoregions with faunas that have the highest proportions of Ramivorous species (≥ 25%) have an average MAP of 1,417 and a median MAP of 1,017 mm. Mean annual precipitation estimates for Patagonian sites equivalent in age to Tinguiririca based on phytolith ecomorphology and stable carbon isotopes are quite low (ca. 200-700 mm; Dunn et al., 2015; Kohn et al., 2015). The relatively arid conditions indicated by these studies imply that high-crowned herbivores from Tinguiririca may have been open habitat specialists (i.e., pure grazers) rather than mixed feeders, as coded here.

The Classification Tree for vegetative cover (Figure 6A) indicates that the paleofauna of Tinguiririca resembles that of a modern open vegetation habitat. Overall the paleoenvironment of Tinguiririca is reconstructed as a temperate, open shrubland; despite the differences in estimated MAP, this is consistent with the previous paleoenvironmental findings of Flynn et al. (2003) and Croft et al. (2008).

Rümikon, Switzerland

The results of the CA, HCA (Figure 2A, Figure 3A, Figure 6), and CT2 (Figure 6B) indicate the paleofauna of Rümikon had an ecological structure most similar to that of a modern mammal community from the Palearctic. The results of the SIMPER and CT3 (Appendix 9, Figure 7) indicate that the paleofauna of Rümikon had a community structure similar to a modern fauna from Biome 7, tropical and subtropical grasslands, savannas, and shrublands. This biome is characterized by semi-arid to sub-humid climates, seasonal rainfall generally ranging from 500-1,300 mm/year, and an open to mixed vegetative cover dominated by shrubs and scattered trees (González-Jiménez, 1979). Biome 7 ecoregions mainly occur in the modern Afrotropic and Neotropic biogeographic realms and are absent from the modern Palearctic biogeographic realm. The estimated MAP of 1,115-1,214 mm for Rümikon based on the linear regression is lower than the general estimate of 1,300-1,500 mm for all fossil localities within the Swiss Upper Freshwater Molasse of Hantke (1984) but within the range of the MAPs from ecoregions of Biome 7 used in this study.

Classification Tree 1 (Figure 6A) indicates that, based on the structure of its fauna, the paleohabitat of Rümikon resembles a modern habitat with a mixed vegetative cover. Overall the paleohabitat of Rumikon is reconstructed as a temperate, mixed woodland, which is consistent with the previous interpretations of Kälin and Kempf (2009). Paleoenvironmental interpretations for units within the Upper Freshwater Molasse, including Rümikon, vary from humid, closed forests to drier, more open habitats (Kälin and Kempf, 2009). Rümikon and other fossil localities within the MN6 units of the Upper Freshwater Molasse contain some genera that are interpreted to indicate humid, forested environments, such as the flying squirrel Albanesia; however, these units contain a greater abundance of specimens from genera that are interpreted to be indicators of drier, more open conditions, such as the ground squirrel Heteroxerus and the lagomorph Lagopsis (Kälin and Kempf, 2009). This has been attributed to a local environmental and climatic change toward more open, less humid habitats (Kälin and Kempf, 2009).

DISCUSSION

Biotic History, Contingency, and Non-Analog Faunas

Perhaps the most noteworthy result of this study is that all four of the South American paleofaunas analyzed most closely resemble modern faunas from outside of (rather than within) the Neotropic biogeographic realm in all analyses. These results confirm our second prediction and demonstrate that the ecological structure of South American paleofaunas was quite different from that of modern mammal communities of that continent. This difference is largely driven by proportions of species in two diet categories, Carnivorous and Ramivorous, and three body mass categories, Size 1, 5, and 6. With the exception of Santa Cruz, South American paleofaunas have fewer carnivorous species than any modern fauna from the Neotropic biogeographic realm; only faunas from Australasia typically have similarly low proportions of carnivorous species (cf. Croft, 2006). Conversely, the South American paleofaunas we analyzed have higher proportions of ramivorous (browsing) species than any modern Neotropical fauna used in this study; these partly or wholly ramivorous taxa are primarily large (size categories 4-6) xenarthrans (ground sloths) and ungulates (notoungulates, litopterns, astrapotheres). On average, these taxa represent of 41% of the species in these paleocommunities. In modern faunas, ramivorous species are primarily represented by extant ungulates such as deer (cervids) and peccaries (tayassuids) and several groups of rodents (e.g., cricetids and echimyids). These taxa account for only 18% of species in modern South American mammal communities, on average. The decrease in the proportion of ramivorous species in modern South American mammal faunas could be related to changes in global climate over the course of the Cenozoic. Janis et al. (2000) suggested that declining temperatures and associated lower primary productivity could account for fewer browsing species in modern North American mammal faunas compared to those of the Miocene. A similar shift could have taken place in South America. The extinction of large (> 44 kg), folivorous (browsing and/or grazing) species during the late Cenozoic (Koch and Barnoksy, 2006; Cione et al., 2009) is another possible, non-exclusive explanation for this difference, as most of the browsers in these paleofaunas were large-bodied taxa.

Another difference between the South American paleofaunas used in this study and modern South American faunas is the lower proportions of Size 1 species. While this could partially reflect a preservation or collection bias (see below for further discussion), this lower proportion could also reflect differences in community structure between pre-, and post-GABI faunas. In modern South American faunas, the species in the Size 1 mass category are predominantly cricetid rodents, a group that dispersed to South America by the late Miocene and today accounts for well over 400 South American species (Barbiere et al., 2016). Prior to the rise of cricetids, Size 1 ecological niches in South America were mainly occupied by metatherians (e.g., argyrolagoids; Zimicz, 2011), which were far less diverse. Most native South American rodents (caviomorphs) were larger or much larger than most modern cricetids (cf. Walton, 1997). It is certainly the case that many small species of metatherians and rodents are unrepresented in the South American fossil record, but it is unlikely that such undiscovered species rival the taxonomic diversity represented by modern cricetids. Perhaps a paucity of very small species in pre-GABI faunas is one of the reasons cricetids experienced such an explosive radiation in South America (Parada et al., 2013; Maestri et al., 2017).

On the other end of the body mass scale, with the exception of Quebrada Honda, South American paleofaunas also have higher proportions of Size 5 species than modern faunas from the Neotropic biogeographic realm. Additionally, no modern fauna within this realm has any species in the Size 6 category, while the species in this category account for 8% and 6% of the species of La Venta and Quebrada Honda, respectively (Appendix 10). The absence of such very large species in modern Neotropical faunas can largely be attributed to the Pleistocene megafaunal extinction (Koch and Barnosky, 2006; Cione et al., 2009; Barnosky and Lindsey, 2010; Villavivencio et al., 2015). The only faunas with similar proportions of species within the Size 5 and 6 variables are from the Afrotropics and Indo-Malaya, the two biogeographic realms that did not experience a significant megafaunal extinction event at the end of Pleistocene (~100 ka) (Barnosky et al., 2004; Sandom et al., 2014; Smith et al., 2018).

Taken together, the ecological characteristics of South American paleofaunas suggest that they are not closely analogous to modern Neotropical faunas. Although the idea of non-analog faunal associations was originally developed in the context of Pleistocene mammal associations (Graham, 2005; Soligo and Andrews, 2005; Semken et al., 2010), it has subsequently been used more generally to refer to faunas that are not closely similar to any modern fauna. The distinctive nature of South American paleofaunas is clearly evident in Figure 5, in which they form a cluster distinct from all modern faunas. This cannot reasonably be interpreted as a sampling bias common to fossil faunas in general, given that: (1) thorough sampling was the primary criterion used to select the fossil faunas analyzed; and (2) the other (non-South American) paleofauna included in the analysis plots separately, among modern faunas. Rather, this cluster of South American paleofaunas is distinguished by its unique ecological characteristics (mainly, a low proportion of Size 1 species and a high proportion of ramivorous species). In this sense, South American paleocommunities can be accurately characterized as broadly non-analog compared to those present today. This conclusion, which is based on the entire mammal community, parallels similar conclusions arising from study of the early Miocene carnivore/predator guild of South America (Croft et al., 2018a). Our findings are also concordant with the idea that non-analog vegetational communities were present in southern South America during the Eocene and Oligocene (Dunn et al., 2015; Kohn et al., 2015). As more and better data ecological data are gathered about South American paleofaunas and localities, we suspect that many other non-analog aspects of South American ecological communities will come to light.

Changes in Global Climate

Changes in climate over the course of the Cenozoic have changed the distributions of biomes and thus faunal communities within biogeographic realms (Kutzbach et al., 1997; Gates, 2003; Hardy, 2003; Peñuelas and Boada, 2003; Ortiz-Jaureguizar and Cladera, 2006; Rowan et al., 2019). The faunal resemblance of Rümikon to Biome 7 (tropical and subtropical grasslands, savannas, and shrublands), a biome absent from the modern Palearctic biogeographic realm, is almost certainly due in large part to the underlying climatic history of Central Europe. In this region, global cooling immediately following the Mid-Miocene Climatic Optimum led to decreased precipitation and more open habitats (Hanke, 1984; Van der Meulen and Daams, 1992; Kälin and Kempf, 2009). Modern central Europe is generally humid and forested, primarily composed of Biomes 4 and 5 (Kottek et al., 2006; Peel et al., 2007). The modern humid conditions in central Europe are largely a result of the relatively recent (~5-3 Ma) closure of the Central American Seaway and the subsequent strengthening of the Atlantic Gulf Stream (Raymo et al., 1989; Murdock et al., 1997; Haug and Tiedemann, 1998; Cederbom et al., 2004). Geographic shifts in biomes have occurred during the Cenozoic in all biogeographic realms, and no single modern biogeographic realm includes all biomes; one to three biomes are missing from each. Thus, restricting an EDA of a paleofauna to a modern comparative dataset from the same biogeographic realm could result in an inaccurate reconstruction due to a priori limits on the breadth of biomes to which it is being compared.

Fidelity of the Fossil Record

The paleohabitat reconstructed for Rümikon based on this EDA differs from those inferred for the South American paleofaunas. This confirms our third prediction of no systematic, sampling bias among paleofaunas that results in spurious similarities between all five paleofaunas. Rümikon differs from the South American sites partly in its higher proportion of species within the Size 1 category (Appendix 10). As noted previously, South American faunas are characterized by relatively few species of this size (3-10%), lower than that of any modern faunas used in this analysis (Appendix 10). By contrast, 35% of species at Rümikon are in the Size 1 category, a proportion within 4% of the mean for most other biogeographic realms (all except Indo-Malaya, at 20%). The relatively low proportion of Size 1 species in the South American paleofaunas could be the result of a preservation bias against small species, a common phenomenon in the fossil record (e.g., Wolff, 1975; Damuth, 1982; Kidwell and Flessa, 1995). However, it seems unlikely that such a bias would affect all four South American paleofaunas relatively uniformly but not Rümikon, particularly considering the geographic and temporal breadth of the South American faunas and the diversity of depositional environments they represent. As noted above, the low proportion of Size 1 species in middle Cenozoic South American paleofaunas is more likely due the continent’s peculiar biogeographical history than a systematic bias in preservation or collection methods.

In a similar vein, the congruence between the analytical results obtained by including and excluding the Size 0 variable indicates that - even if a preservation or collection bias against very small (≤ 10 g) species is a common feature of Cenozoic mammal paleofaunas - it is unlikely to have a major effect on reconstructing community structure and paleohabitat. For the four South American paleofaunas, reconstructed habitats were the same with and without the Size 0 variable in all analyses. Nearly the same was true for Rümikon; excluding the Size 0 variable did not affect the overall habitat reconstruction, though the paleofauna was grouped with Biome 10 rather than Biome 7 in one analysis (the Classification Tree for biome prediction, CT3). However, this had no effect on the overall paleohabitat reconstruction, as Rümikon lacks Size 0 species, and Rümikon was also found to have a similar faunal composition to Biome 7 in the SIMPER analysis (Appendix 9).

Scales of Comparison

When performing a paleoecological analysis such as this, an important factor to consider is the geographic area represented by both the modern and fossil faunas (Kovarovic et al., 2018). Ideally, the areas of the two should be similar, though this can be challenging to achieve in practice due to a variety of factors. For example, modern species lists typically derive from national parks and other protected areas, which vary substantially in size; these most often represent gamma diversity (the diversity of a landscape or geographic area) rather than alpha diversity (diversity within a single community, typically < 10 km2; Whittaker, 1977; Stoms and Estes, 1993) and can differ from one another in area by several orders of magnitude (cf. Kay and Madden, 1997). The area represented by a fossil list can also vary depending on several factors; limited outcrop exposure can lead to underestimates of geographic area, while time-averaging can have the opposite effect, particularly in the context of habitats that shifted spatially through time. As discussed extensively by Croft (2013), it has not yet been determined how these different spatial scales can affect ecological diversity analyses of fossil faunas. Similarly, it is not possible to know the degree to which our findings are the result of using regional faunas as opposed to more restricted ones (e.g., Kay et al., 2012; Spradley et al., 2019). These approaches are likely complimentary to one another; studies based on more restricted geographic areas may be more suitable for incorporating finer-scale habitat and climate data, while those based on larger areas may be more appropriate for revealing general or large-scale patterns.

CONCLUSIONS AND FUTURE DIRECTIONS

Accurate paleohabitat reconstructions are essential for understanding the evolution and ecologies of paleofaunas. Our study is the first EDA of South American paleofaunas using a comparative dataset of modern temperate and tropical faunas from multiple biogeographic realms. The paleohabitat reconstructions for the South American paleofaunas indicate that modern mammalian faunas from other continents represent better ecological analogues for these mammalian paleocommunities than those from within the Neotropic biogeographic realm. The similarities between Rümikon, a European fossil fauna, and faunas from a biome not currently represented in the Palearctic biogeographic realm further underscores the importance of using global comparative datasets for accurate paleohabitat reconstructions, even in areas that have not experienced faunal changes of the same magnitude as in South America.

We envision several ways in which the accuracy and precision of ecological diversity analyses could be improved in the future, besides using a global comparative dataset. The incorporation of smaller-scale modern faunal communities from within the predicted biomes and biogeographic realms, such as those in nature reserves and natural parks, could be used to increase the precision of the paleohabitat reconstructions through the use of faunal lists based on alpha (locality), rather than gamma (regional) data and species richness. The integration of modern temperature and elevation data from the habitats of the smaller-scale faunal communities could also further refine the habitat reconstructions to a local level and be used to test hypotheses regarding the effects of local environmental differences, the spatial dynamics, diversification, and provinciality of paleofaunas. This method and this global comparative dataset can also provide essential habitat context when other paleoenvironmental indicators are not present and can be used in conjunction with other methods, such as paleopedology and ichnology, to create robust, multidisciplinary paleoenvironmental interpretations.

ACKNOWLEDGEMENTS

We thank K. Kovarovic and one anonymous reviewer for comments and suggestions that significantly improved and expanded upon the scope of this paper. We thank M. Sameh and B. Townsend for inspiring some of the ideas and techniques used in this study, D. Su for input and insight, and R. Patel for support and for assistance with programming R. This research was supported by Case Western Reserve University and the National Science Foundation (EAR 0958733 and EAR 1423058 to DAC).

REFERENCES

Andrews, P. and Evans, E.N. 1979. The environment of Ramapithecus in Africa. Paleobiology, 5:22-30. https://doi.org/10.1017/s0094837300006266

Andrews, P., Lord, J.M., and Evans E.M.N. 1979. Patterns of ecological diversity in fossil and modern mammalian faunas. Biological Journal of the Linnean Society, 11:177-205. https://doi.org/10.1111/j.1095-8312.1979.tb00034.x

Andrews, P. and Van Couvering, J.H. 1975. Palaeoenvironments in the East African Miocene. Approaches to Primate Paleobiology, 5:62-103.

Antoine, P.O., Marivaux, L., Croft, D.A., Billet, G., Ganerød, M., Jarmillo, C., Martin, T., Orliac, M.v.J., Tejada, J., Altamirano, A.J., Dunanthon, F., Fanjat, G.g., Rousse, S., and Gismondi, R.S. 2012. Middle Eocene rodents from Peruvian Amazonia reveal the pattern and timing of caviomorph origins and biogeography. Proceedings of the Royal Society B: Biological Sciences, 279:1319-1326. https://doi.org/10.1098/rspb.2011.1732

Barbiere, F., Ortiz, P.E., and Pardiñas, U.F.J. 2016. Evolución de los roedores Cricetidae en América del Sur: una visión desde el registro paleontológico. Contribuciones del Museo Argentino de Ciencias Naturales, 6:335-347.

Barnosky, A.D., Koch, P.L., Feranec, R.S., Wing, S.L., and Shabel, A.B. 2004. Assessing the causes of Late Pleistocene extinctions on the continents. Science, 306:70-75. https://doi.org/10.3410/f.1021999.250627

Barnosky, A.D. and Lindsey, E.L. 2010. Timing of Quaternary megafaunal extinction in South America in relation to human arrival and climate change. Quaternary International, 217:10-29. https://doi.org/10.1016/j.quaint.2009.11.017

Bell, R.H.V. 1971. A grazing ecosystem in the Serengeti. Scientific American, 225:86-93.

Bertrand, O.C., Flynn, J.J., Croft, D.A., and Wyss, A.R. 2012. Two new taxa (Caviomorpha, Rodentia) from the early Oligocene Tinguiririca fauna (Chile). American Museum Novitates, 3750:1-6. https://doi.org/10.1206/3750.2

Blain, H-A., Bailon, S., Cuenca-Bescós, G. 2008, The Early-Middle Pleistocene palaeoenvironmental change based on the squamate reptile and amphibian proxies at the Gran Dolina site, Atapuerca, Spain. Palaeogeography, Palaeoclimatology, Palaeoecology, 261:177-192. https://doi.org/10.1016/j.palaeo.2008.01.015

Boggs Jr., S. 2006. Principles of Sedimentology and Stratigraphy (fourth edition). Pearson Education, California.

Bond, M., Perea, D., Ubilla, M., and Tauber, A. 2001. Neolicaphrium recens Frenguelli, 1921, the only surviving Proterotheriidae (Litopterna, Mammalia) into the South American Pleistocene. Palaeovertebrata, 30:37-50.

Bond, M., Tejedor, M.F., Campbell Jr., K. E., Chornogubsky, L., Novo, N., and Goin, F. 2015. Eocene primates of South America and the African origins of New World monkeys. Nature, 520:538. 

Bradham, J., Flynn, J.J., Croft, D.A., and Wyss, A.R. 2015. New notoungulates (Notostylopidae and basal toxodontians) from the early Oligocene Tinguiririca fauna of the Andean Main Range, central Chile. American Museum Novitates, 3841:1-24. https://doi.org/10.1206/3841.1

Brandoni, D., Carlini, A.A., Anaya, F., Gans, P., and Croft, D.A. 2018. New remains of Megathericulus patagonicus Ameghino (Xenarthra, Tardigrada) from the Serravallian (middle Miocene) of Bolivia; chronologic and biogeographic implications. Journal of Mammalian Evolution, 25:327-337. https://doi.org/10.1007/s10914-017-9384-y

Breiman, L., Friedman, J.H., Olshen, R.A., and Stone, C.G. 1984. Classification and Regression Trees. Wadsworth International Group, Belmont, California, USA.

Bronger, A. and Catt. J.A. 1989. Paleosols: problems of definition, recognition and interpretation, p. 1-7. In Bronger, A. and Catt, J.A. (eds.), Paleopedoloqy, Nature and Application of Paleosols, Catena Supplement. Cremlingen, Germany.

Brown, J.H. 1995. Macroecology (first edition). University of Chicago Press, Chicago, United States.

Buatois, L.A. and Mángano, M.G. 1995. The paleoenvironmental and paleoecological significance of the lacustrine Mermia ichnofacies: an archetypical subaqueous nonmarine trace fossil assemblage. Ichnos, 4:151-161. https://doi.org/10.1080/10420949509380122

Cadena, E.A., Anaya, F., and Croft, D.A. 2015. Giant fossil tortoise and freshwater chelid turtle remains from the middle Miocene, Quebrada Honda, Bolivia: Evidence for lower paleoelevations for the southern Altiplano. Journal of South American Earth Science, 64:190-198. https://doi.org/10.1016/j.jsames.2015.10.013

Campbell, N.A. 1996. Biology (fourth edition). The Benjamin/Cummings Publishing Company, Inc., Menlo Park, California.

Catena, A.M., Hembree, D.I., Saylor, B.Z., Anaya, F., and Croft, D.A. 2017. Paleosol and ichnofossil evidence for significant Neotropical habitat variation during the late middle Miocene (Serravallian). Palaeogeography, Palaeoclimatology, Palaeoecology, 459:423-439. https://doi.org/10.1016/j.palaeo.2017.09.024

Caudill, M.R., Driese, S.G., and Mora, C.I. 1996. Preservation of a paleo-vertisol and an estimate of Late Mississippian paleoprecipitation. Journal of Sedimentary Research, 66:58-70. https://doi.org/10.1306/d42682b1-2b26-11d7-8648000102c1865d

Caughley, G. 1987. The distribution of eutherian body weights. Oecologia, 74:319-320. https://doi.org/10.1007/bf00379377

Ceballos, G. and Ehrlich, P.R. 2006. Global mammal distributions, biodiversity hotspots, and conservation. Proceedings of the National Academy of Sciences, 103:19374-19379. https://doi.org/10.1073/pnas.0609334103

Cederbom, C.E., Sinclair, H.D., Schlunegger, F., and Rahn, M.K. 2004. Climate-induced rebound and exhumation of the European Alps. Geology, 32:709-712. https://doi.org/10.1130/g20491.1

Cerling, T.E., Quade, J., Wang, Y., and Bowman, J.R. 1989. Carbon isotopes in soils and palaeosols as ecology and palaeoecology indicators. Nature, 341:138-139. https://doi.org/10.1038/341138a0

Charrier, R., Wyss, A.R., Flynn, J.J. Swisher III, C.C., Norell, M.A., Zapatta, F., McKenna, M.C., and Novacek, M.J. 1996. New evidence for late Mesozoic-early Cenozoic evolution of the Chilean Andes in the upper Tinguiririca valley (35° S), central Chile. Journal of South American Earth Science, 9:393-422. https://doi.org/10.1016/s0895-9811(96)00035-1

Cifelli, R.L. 1985. South American ungulate evolution and extinction, p. 249-266. In Stehli, F.G. and Webb, S.D. (eds.), The Great American Biotic Interchange. Springer, Boston. https://doi.org/10.1007/978-1-4684-9181-4_9

Cione, A.L., Gasparini, G.M., Soibelzon, E., Soibelzon, L.H., and Tonni, E.P. 2015. The GABI in southern South America, p. 71-96. In Cione, A.L., Gasparini, G.M., Soibelzon, E., Soibelzon, L.H., and Tonni, E.P. (eds.), The Great American Biotic Interchange: A South American Perpective. Springer, Dordechet. https://doi.org/10.1007/978-94-017-9792-4_3

Cione, A.L., Tonii, E.P., and Soibelzon, L. 2009. Did humans cause the late Pleistocene-early Holocene mammalian extinctions in South America in a context of shrinking open areas? p. 125-144. In Haynes, G. (ed.), American Megafaunal Extinctions at the end of the Pleistocene. Springer, New York. https://doi.org/10.1007/978-1-4020-8793-6_7

Clarke, K.R. 1993. Non‐parametric multivariate analyses of changes in community structure. Austral Ecology, 18:117-143. https://doi.org/10.1111/j.1442-9993.1993.tb00438.x

Croft, D.A. 2000. Archaeohyracidae (Mammalia: Notoungulata) from the Tinguiririca Fauna, Central Chile, and the Evolution and Paleoecology of South American Mammalian Herbivores. Unpublished PhD Thesis, University of Chicago, Department of Organismal Biology and Anatomy.

Croft, D. A. 2006. Do marsupials make good predators? Insights from predator-prey diversity ratios. Evolutionary Ecology Research, 8:1193-1214.

Croft, D.A. 2007. The middle Miocene (Laventan) Quebrada Honda Fauna, southern Bolivia, and a description of its notoungulates. Palaeontology, 50:277-303. https://doi.org/10.1111/j.1475-4983.2006.00610.x

Croft, D.A. 2012. A synthesis of Cenozoic Neotropical mammal evolution in South America: Biogeography and influences from higher latitudes. Journal of Vertebrate Paleontology, SVP Program and Abstracts Book, 2013:236.

Croft, D.A. 2013. What constitutes a fossil mammal community in the early Miocene Santa Cruz Formation? Journal of Vertebrate Paleontology, 33:401-409. https://doi.org/10.1080/02724634.2013.722154

Croft, D.A. 2016. Horned Armadillos and Rafting Monkeys: The Fascinating Fossil Mammals of South America. Indiana University Press, Bloomington, Indiana.

Croft, D.A. and Anaya, F. 2006. A new middle Miocene hegetotheriid (Notoungulata: Typotheria) and a phylogeny of the Hegetotheriidae. Journal of Vertebrate Paleontology, 26:387-399. https://doi.org/10.1671/0272-4634(2006)26[387:anmmhn]2.0.co;2

Croft, D.A., Bond, M., Flynn, J.J., Reguero, M., and Wyss, A.R. 2003. Large archaeohyracids (Typotheria, Notoungulata) from central Chile and Patagonia including a revision of Archaeotypotherium. Fieldiana: Geology (N.S.), 49:1-38. https://doi.org/10.5962/bhl.title.5264

Croft, D.A., Chick, J., and Anaya, F. 2011. New middle Miocene caviomorph rodents from Quebrada Honda, Bolivia. Journal of Mammal Evolution, 18:245-268. https://doi.org/10.1007/s10914-011-9164-z

Croft, D.A., Engelman, R.K., Dolgushina, T., and Wesley, G. 2018a. Diversity and disparity of sparassodonts (Metatheria) reveal non-analogue nature of ancient South American mammalian carnivore guilds. Proceedings of the Royal Society B: Biological Sciences, 285:20172012. https://doi.org/10.1098/rspb.2017.2012

Croft, D.A., Flynn, J.J., and Wyss, A.R. 2008. The Tinguiririca Fauna of Chile and the early stages of "modernization" of South American mammal faunas. Arquivos do Museu Nacional, Rio de Janeiro, 66:191-211.

Croft, D.A., Su, D.F., and Simpson, S.W. 2018b. Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Springer (Vertebrate Paleobiology and Paleoanthropology Series), Cham, Switzerland. https://doi.org/10.1007/978-3-319-94265-0

Croft, D.A. and Townsend, K.E. 2005. Inferring habitat for the late early Miocene Santa Cruz Fauna (Santa Cruz Province, Argentina) using ecological diversity analysis. Journal of Vertebrate Paleontology, 25:48A.

Cruz, J.A., Arroyo-Cabrales, J., and Reynoso, V.H. 2016. Reconstructing the paleoenvironment of Loltún Cave, Yucatán, Mexico, with Pleistocene amphibians and reptiles and their paleobiogeographic implications. Revista Mexicana de Ciencias Geológicas, 33:342-354.

Cumming, D.H.M. 1982. The influence of large herbivores on savanna structure in Africa, p. 217-245. In Huntley, B.J. and Walker, B.H. (eds.), Ecology of Tropical Savannas. Springer, Berlin. https://doi.org/10.1007/978-3-642-68786-0_11

Damuth, J. 1982. Analysis of the preservation of community structure in assemblages of fossil mammals. Paleobiology, 8:434-446. https://doi.org/10.1017/s009483730000717x

Damuth, J.D. 1992. Taxon-free characterization of animal communities, p. 183-203. In Behrensmeye, A.K., Damuth, J.D., DiMichele, W.A., Potts, R., Sues, H.D., and Wing, S.L. (eds.), Terrestrial Ecosystems Through Time. University of Chicago, Chicago.

Dark, P. 2008. Paleoenvironmental reconstruction, methods, p. 1787-1790. In Pearsall, D.M. (ed.), Encyclopedia of Archaeology. Oxford, UK. https://doi.org/10.1016/b978-012373962-9.00226-0

Dasmann, R.F. 1973. A system for defining and classifying natural regions for purposes of conservation. A progress report. IUCN Occasional Paper 7, International Union for Conservation of Nature and Natural Resources, Morges, Switzerland.

Dasmann, R.F. 1974. Biotic provinces of the world: Further development of a system for defining and classifying natural regions for purposes of conservation. IUCN Occasional Paper 9, International Union for Conservation of Nature and Natural Resources, Morges, Switzerland.

Davidson, A.D., Detling, J.K., and Brown, J.H. 2012. Ecological roles and conservation challenges of social, burrowing, herbivorous mammals in the world's grasslands. Frontiers in Ecology and the Environment, 10:477-486. https://doi.org/10.1890/110054

De'ath, G. and Fabricius, K.E. 2000. Classification and regression trees: a powerful yet simple technique for ecological data analysis. Ecology, 81:3178-3192. https://doi.org/10.2307/177409

DeGusta, D. and Vrba, E. 2003. A method for inferring paleohabitats from the functional morphology of bovid astragali. Journal of Archaeoogical Science, 30:1009-1022. https://doi.org/10.1016/s0305-4403(02)00286-8

DeGusta, D. and Vrba, E. 2005. Methods for inferring paleohabitats from the functional morphology of bovid phalanges. Journal of Archaeoogical Science, 32:1099-1113. https://doi.org/10.1016/j.jas.2005.02.010

Driese, S.G. and Ober, E.G. 2005. Paleopedologic and paleohydrologic records of precipitation seasonality from Early Pennsylvanian "underclay" paleosols, USA. Journal of Sedimentary Research, 75:997-1010. https://doi.org/10.2110/jsr.2005.075

Dunn, R.E., Strömberg, C.A., Madden, R.H., Kohn, M.J., and Carlini, A.A. 2015. Linked canopy, climate, and faunal change in the Cenozoic of Patagonia. Science, 347:258-261. https://doi.org/10.1126/science.1260947

Eisenberg, J.F. 1981. The Mammalian Radiations. University of Chicago Press, Chicago.

Engelman, R.K., Anaya, F., and Croft, D.A. 2015. New specimens of Acyon myctoderos (Metatheria, Sparassodonta) from Quebrada Honda, Bolivia. Ameghiniana, 52:204-225. https://doi.org/10.5710/amgh.19.11.2014.2803

Engelman, R.K., Anaya, F., and Croft, D.A. 2017. New palaeothentid marsupials (Paucituberculata) from the middle Miocene of Quebrada Honda, Bolivia, and their implications for the palaeoecology, decline and extinction of the Palaeothentoidea. Journal of Systematic Palaeontology, 15:787-820. https://doi.org/10.1080/14772019.2016.1240112

Engelman, R.K., Anaya, F., and Croft, D.A. 2018. Australogale leptognathus gen. et. sp. nov., a second species of small sparassodont (Mammalia: Metatheria) from the middle Miocene locality of Quebrada Honda, Bolivia. Journal of Mammalian Evolution, 2018:1-8. https://doi.org/10.1007/s10914-018-9443-z

Engelman. R.K. and Croft, D.A. 2014. A new small-bodied sparassodont (Mammalia: Metatheria) from the middle Miocene locality of Quebrada Honda, Bolivia. Journal of Vertebrate Paleontology, 34:672-688. https://doi.org/10.1080/02724634.2013.827118

Fedriani, J.M., Fuller, T.K., and Sauvajot, R.M. 2001. Does availability of anthropogenic food enhance densities of omnivorous mammals? An example with coyotes in southern California. Ecography, 24:325-331. https://doi.org/10.1111/j.1600-0587.2001.tb00205.x

Fick, S.E. and Hijmans, R.J. 2017. WorldClim 2: new 1‐km spatial resolution climate surfaces for global land areas. International Journal of Climatology, 37:4302-4315. https://doi.org/10.1002/joc.5086

Fiedel, S. 2009. Sudden deaths: the chronology of terminal Pleistocene megafaunal extinction, p. 21-37. In Haynes, G. (ed.), American Megafaunal Extinctions at the end of the Pleistocene. Springer, New York. https://doi.org/10.1007/978-1-4020-8793-6_2

Finlayson, C., Carrión, J., Brown, K., Finlayson, G., Sánchez-Marco, A., Fa, D., Rodríguez-Vidal, J., Fernández, S., Fierro, E., Bernal-Gómez, M., and Giles-Pacheco, F. 2011. The Homo habitat niche: using the avian fossil record to depict ecological characteristics of Palaeolithic Eurasian hominins. Quaternary Science Reviews, 30:1525-1532. https://doi.org/10.1016/j.quascirev.2011.01.010

Fleagle, J.G., Perkins, M.E., Heizler, M.T., Nash, B., Bown, T.M., Tauber, A.A., Dozo, M.T., and Tejedor, M.F. 2012. Absolute and relative ages of fossil localities in the Santa Cruz and Pinturas Formations, p. 41-58. In Vizcaíno, S.F., Kay, R.F, and Bargo, M.S. (eds.), Early Miocene Paleobiology in Patagonia: High-Latitude Paleocommunities of the Santa Cruz Formation. Cambridge University Press. https://doi.org/10.1017/cbo9780511667381.004

Fleming, T.H. 1973. Numbers of mammal species in North and Central American forest communities. Ecology, 54:555-563. https://doi.org/10.2307/1935340

Flynn, J.J., Charrier, R., Croft, D.A., and Wyss, A.R. 2012. Cenozoic Andean faunas: Shedding new light on South American mammal evolution, biogeography, environments, and tectonics, p. 51-75. In Patterson, B.D. and Costa, L.P. (B.D. Patterson and L.P. Costa (eds.), Bones, Clones, and Biomes: The History and Geography of Recent Neotropical Mammals. University of Chicago Press, Chicago.

Flynn, J., Guerrero, J., and Swisher III, C.C. 1997. Geochronology of the Honda group, p. 4-59. In Kay, R.F., Madden, R.H., Cifelli, R.H., and Flynn, J.J. (eds.), Vertebrate Paleontology in the Neotropics: The Miocene Fauna of La Venta, Colombia. Smithsonian Institution Press, Washington, DC.

Flynn, J.J., Wyss, A.R., Croft, D.A., and Charrier, R. 2003. The Tinguiririca Fauna, Chile: biochronology, paleoecology, biogeography, and a new earliest Oligocene South American Land Mammal "Age". Palaeogeography, Palaeoclimatology, Palaeoecology, 195:229-259. https://doi.org/10.1016/s0031-0182(03)00360-2

Fooden, J. 1972. Breakup of Pangaea and isolation of relict mammals in Australia, South America, and Madagascar. Science, 175:894-898. https://doi.org/10.1126/science.175.4024.894

Forman, R.T.T. and Godron, M. 1981. Patches and structural components for a landscape ecology. BioScience, 31:733-740. https://doi.org/10.2307/1308780

Gainsbury, A.M., Tallowin, O.J., and Meiri, S. 2018. An updated global data set for diet preferences in terrestrial mammals: testing the validity of extrapolation. Mammal Review, 48:160-167. https://doi.org/10.1111/mam.12119

Gardezi, T. and da Silva, J. 1999. Diversity in relation to body size in mammals: a comparative study. American Naturalist, 153:110-123. https://doi.org/10.2307/2463899

Gates, D.M. 1993. Climate Change and Its Biological Consequences. Sinauer, Massachusetts.

Gibson, D.J. 2009. Grasses and Grassland Ecology. Oxford University Press, Oxford, UK.

González-Jiménez, E. 1979. Primary and secondary productivity in flooded savannas, p. 620-625. In Spredding, C.R.W. (ed.), Tropical Grazing Land Ecosystem: a ‘State-of-Knowledge’ report prepared by UNESCO. Paris. 

Graham R.W. 2005. Quaternary mammal communities: Relevance of the individualistic response and non-analogue faunas. The Paleontological Society Papers, 11:141-58. https://doi.org/10.1017/s1089332600001297

Gregory-Wodzicki, K.M., McIntosh, W.C., and Velasquez, K. 1998. Climate and tectonic implications of the late Miocene Jakokkota flora, Bolivian Altiplano. Journal of South American Earth Sciences, 11:533-560. https://doi.org/10.1016/s0895-9811(98)00031-5

Guerrero, J. 1997. Stratigraphy, sedimentary environments, and the Miocene uplift of the Colombian Andes, p. 13-43. In Kay, R.F., Madden, R.H., Cifelli, R.H., and Flynn, J.J. (eds.), Vertebrate Paleontology in the Neotropics: The Miocene Fauna of La Venta, Colombia. Smithsonian Institution Press, Washington, DC.

Guo, Q., Rundel, P.W., and Goodall, D.W. 1998. Horizontal and vertical distribution of desert seed banks: patterns, causes, and implications. Journal of Arid Environments, 38:465-478. https://doi.org/10.1006/jare.1997.0353

Hammer, Ø. and Harper, D.A.T. 2006. Paleontological Data Analysis, Blackwell Publishing, Malden, MA, USA. https://doi.org/10.1002/9780470750711

Hammer, Ø., Harper, D.A.T., and Ryan, P.D. 2001. PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica, 4.1.4:1-9. http://palaeo-electronica.org/2001_1/past/issue1_01.htm

Hantke, R. 1984. Floreninhalt, biostratigraphische Gliederung und Paläoklima der mittelmiozänen Oberen Süßwassermolasse (OSM) der Schweiz und ihrer nördlichen Nachbargebiete. Heimatliche Schriftenreihe des Landkreises Günzburg, 2:47-53.

Hardy, J.T. 2003. Climate Change: Causes, Effects, and Solutions. John Wiley & Sons, Chichester, UK.

Hasiotis, S.T. 2002. Continental Trace Fossils. SEPM Short Course Notes 51, Tulsa, Oklahoma.

Hasiotis, S.T. 2007. Continental ichnology: fundamental processes and controls on trace fossil distribution, p. 268-284. In Miller III, W. (ed.), Trace Fossils: Concepts, Problems, Prospects. Elsevier, Amsterdam. https://doi.org/10.1016/b978-044452949-7/50142-x

Haug, G.H. and Tiedemann, R. 1998. Effect of the formation of the Isthmus of Panama on Atlantic Ocean thermohaline circulation. Nature, 393:673-676. https://doi.org/10.1038/31447

Hembree, D. 2018. The role of continental trace fossils in Cenozoic paleoenvironmental reconstructions, p. 185-214. In Croft, D.A., Simpson, S.W., and Su, D.F (eds.), Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Springer (Vertebrate Paleobiology and Paleoanthropology Series), Cham, Switzerland. https://doi.org/10.1007/978-3-319-94265-0_10

Hembree, D.I., Johnson, L.M., and Tenwalde, R.W. 2012. Neoichnology of the desert scorpion Hadrurus arizonensis: burrows to biogenic cross lamination. Palaeontologia Electronic a, 15.1.10A:1-34. https://doi.org/10.26879/292
palaeo-electronica.org/content/2012-issue-1-articles/192-neoichnology-of-scorpions

Hjelmroos, M. and Franzen, L.G. 1994. Implications of recent long-distance pollen transport events for the interpretation of fossil pollen records in Fennoscandia. Review of Palaeobotany and Palynology, 82:175-189. https://doi.org/10.1016/0034-6667(94)90029-9

Holdridge, L.R., Grenke, W.C., Hatheway, W.H., Liang, T., and Tosi, J.A. 1971. Forest Environments in Tropical Life Zones: A Pilot Study. Oxford Pergamon Press, Oxford.

Hooghiemstra, H. 1989. Quaternary and Upper-Pliocene glaciations and forest development in the tropical Andes: evidence from a long high-resolution pollen record from the sedimentary basin of Bogotá, Columbia. Palaeogeography, Palaeoclimatology, Palaeoecology, 72:11-26. https://doi.org/10.1016/0031-0182(89)90129-6

Howard, J.D. 1975. The sedimentological significance of trace fossils, p. 131-146. In Frey, R.W. (ed.), The Study of Trace Fossils. Springer-Verlag, New York. https://doi.org/10.1007/978-3-642-65923-2_8

Janis, C.M. 1988. Estimation of tooth volume and hypsodonty indices in ungulate mammals, and the correlation of these factors with dietary preference, p 367-387. In Russell, E., Santoro, J.P., and Signogneau-Russell, D. (eds.), Teeth Revisited: Proceedings VII International Symposium on Dental Morphology, Memoires du Museum National d’Histoire Naturelle, Serie C.

Janis, C.M., Damuth, J., and Theodor, J.M. 2000. Miocene ungulates and terrestrial primary productivity: where have all the browsers gone? Proceedings of the National Academy of Sciences, 97:7899-7904. https://doi.org/10.1073/pnas.97.14.7899

Jarman, P.J. 1971. Diets of large mammals in the woodlands around Lake Kariba, Rhodesia. Oecologia, 8:157-178. https://doi.org/10.1007/bf00345811

JMP®, Version 13.0.0. 1989-2007. SAS Inc., Cary, NC.

Kälin, D. and Engesser, B. 2001. Die jungmiozäne Säugetierfauna vom Nebelbergweg bei Nunningen (Kanton Solothurn, Schweiz). Schweizerische Paläontologische Abhandlungen, 121:1-61.

Kälin, D. and Kempf, O. 2009. High-resolution stratigraphy from the continental record of the Middle Miocene Northern Alpine Foreland Basin of Switzerland. Neues Jahrbuch für Geologie und Paläontologie-Abhandlungen, 254:177-235. https://doi.org/10.1127/0077-7749/2009/0010

Kappelman, J. 1988. Morphology and locomotor adaptations of the bovid femur in relation to habitat. Journal of Morphology, 198:19-130. https://doi.org/10.1002/jmor.1051980111

Kay, R.F., Madden R.H., Cifelli R.L., and Flynn J.H. 1997. Vertebrate Paleontology in the Neotropics: The Miocene Fauna of La Venta, Colombia. Smithsonian Institution Press, Washington, D.C.

Kay, R.F. and Madden, R.H. 1997a. Paleogeography and paleoecology, p. 520-550. In Kay, R.F., Madden, R.H., Cifelli, R.H., and Flynn, J.J. (eds.), Vertebrate Paleontology in the Neotropics: The Miocene Fauna of La Venta, Colombia, Smithsonian Institution Press, Washington, DC.

Kay, R.F. and Madden, R.H. 1997b. Mammals and rainfall: paleoecology of the middle Miocene at La Venta (Colombia, South America). Journal of Human Evolution, 32:161-199. https://doi.org/10.1006/jhev.1996.0104

Kay, R.F., Vizcaíno, S.F., and Bargo, M.S. 2012. A review of the paleoenvironment and paleoecology of the Miocene Santa Cruz Formation, p. 331-365. In Vizcaíno, S.F., Kay, R.F., and Bargo, M.S. (eds.), Early Miocene Paleobiology in Patagonia: High-Latitude Paleocommunities of the Santa Cruz Formation. Cambridge University Press, Cambridge. https://doi.org/10.1017/cbo9780511667381.018

Kidwell, S. M. and Flessa, K.W. 1995. The quality of the fossil record: populations, species, and communities. Annual Review of Ecology and Systematics, 26:269-299. https://doi.org/10.1146/annurev.ecolsys.26.1.269

Kinlaw, A.L. 1999. A review of burrowing by semi-fossorial vertebrates in arid environments. Journal of Arid Environments, 41:127-145. https://doi.org/10.1006/jare.1998.0476

Koch, P.L. and Barnosky, A.D. 2006. Late Quaternary extinctions: state of the debate. Annual Review of Ecology and Systematics, 37:215-260. https://doi.org/10.1146/annurev.ecolsys.34.011802.132415

Kohn, M.J., Strömberg, C.A., Madden, R.H., Dunn, R.E., Evans, S., Palacios, A., and Carlini, A.A. 2015. Quasi-static Eocene-Oligocene climate in Patagonia promotes slow faunal evolution and mid-Cenozoic global cooling. Palaeogeography, Palaeoclimatology, Palaeoecology, 435:24-37. https://doi.org/10.1016/j.palaeo.2015.05.028

Kottek, M., Grieser, J., Beck, C., Rudolf, B., and Rubel, F. 2006. World map of the Köppen-Geiger climate classification updated. Meteorologische Zeitschrift, 15:259-263. https://doi.org/10.1127/0941-2948/2006/0130

Kovarovic, K., Andrews P., and Aiello L. 2002. The palaeoecology of the Upper Ndolanya Beds at Laetoli, Tanzania. Journal of Human Evolution, 43:395-418.

Kovarovic, K., Su, D.F., and Lintulaakso, K. 2018. Mammal community structure analysis p. 351-372. In Croft, D.A., Simpson, S.W., and Su, D.F (eds.), Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Springer (Vertebrate Paleobiology and Paleoanthropology Series), Cham, Switzerland.

Kozłowski, J. and Gawelczyk, A. 2002. Why are species’ body size distributions usually skewed to the right? Functional Ecology, 16:419-432. https://doi.org/10.1046/j.1365-2435.2002.00646.x

Kraus, M.J. 1999. Paleosols in clastic sedimentary rocks: their geologic applications. Earth Science Review, 47:41-70. https://doi.org/10.1016/s0012-8252(99)00026-4

Kraus, M.J. and Bown, T.M. 1986. Paleosols and time resolution in alluvial stratigraphy, p. 180-207. In Wright, V.P. (ed.), Paleosols: Their Recognition and Interpretation, Blackwell Scientific, Oxford.

Kraus, M.J. and Riggins, S. 2007. Transient drying during the Paleocene-Eocene Thermal Maximum (PETM): analysis of paleosols in the Bighorn Basin, Wyoming. Palaeogeography, Palaeoclimatology, Palaeoecology, 245:444-461. https://doi.org/10.1016/j.palaeo.2006.09.011

Kutzbach, J.E., Ruddiman, W.F., and Prell, W.L. 1997. Possible effects of Cenozoic uplift and CO2 lowering on global and regional hydrology, p 149-170. In Ruddiman, W.F. (ed.), Tectonic Uplift and Climate Change. Springer, Boston. https://doi.org/10.1007/978-1-4615-5935-1_7

Labhart, T.P. 2005. Geologie der Schweiz, seventh edition. Ott Verlag, Bern, Switzerland.

Legendre, P. and Legendre, L. 2012. Numerical Ecology (third edition). Elsevier, Amsterdam.

Lintulaakso, K. and Kovarovic, K. 2016. Diet and locomotion, but not body size, differentiate mammal communities in worldwide tropical ecosystems. Palaeogeography, Palaeoclimatology, Palaeoecology, 454:20-29. https://doi.org/10.1016/j.palaeo.2016.04.012

Losos, J.B. 1990. Ecomorphology, performance capability, and scaling of West Indian Anolis lizards: an evolutionary analysis. Ecological Monographs, 60:369-388. https://doi.org/10.2307/1943062

Louchart, A., Wesselman, H., Blumenschine, R.J., Hlusko, L.J., Njau, J.K., Black, M.T., Asnake, M., and White, T.D. 2009. Taphonomic, avian and small-vertebrate indicators of Ardipithecus ramidus habitat. Science, 326:66-66e4. https://doi.org/10.1126/science.1175823

Louys, J., Meloro, C., Elton, S., Ditchfield, P., and Bishop, L.C. 2011. Mammal community structure correlates with arboreal heterogeneity in faunally and geographically diverse habitats: implications for community convergence. Global Ecology and Biogeography, 20:717-729. https://doi.org/10.1111/j.1466-8238.2010.00643.x

Macdonald, D. 2006. The Encyclopedia of Mammals. Oxford University Press, Oxford.

MacFadden, B.J., Anaya, F., Perez, H., Naeser, C.W., Zeitler, P.K., and Campbell Jr., K.E. 1990. Late Cenozoic paleomagnetism and chronology of Andean basins of Bolivia: evidence for possible oroclinal bending. Journal of Geology, 98:541-555. https://doi.org/10.1086/629423

MacFadden, B.J. and Wolff, R.G. 1981. Geological investigations of Late Cenozoic vertebrate-bearing deposits in southern Bolivia. Anais do II Congresso Latino-Americano de Paleontología, 2:765-778.

Maestri, R., Monteiro, L.R., Fornel, R., Upham, N.S., Patterson, B.D., and de Freitas, T.R.O. 2017. The ecology of a continental evolutionary radiation: is the radiation of sigmodontine rodents adaptive? Evolution, 71:610-632. https://doi.org/10.1111/evo.13155

Mander, L. and Punyasena, S.W. 2018. Fossil pollen and spores in paleoecology, p. 215-234. In Croft, D.A., Simpson, S.W., and Su, D.F (eds.), Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Springer (Vertebrate Paleobiology and Paleoanthropology Series), Cham, Switzerland. https://doi.org/10.1007/978-3-319-94265-0_11

Mares, M.A. and Willig, M.R. 1994. Inferring biome associations of Recent mammals from samples of temperate and tropical faunas: paleoecological considerations. Historical Biology, 8:31-48.

Marshall, L.G. and Cifelli, R.L. 1990. Analysis of changing diversity patterns in Cenozoic land mammal age faunas, South America. Palaeovertebrata, 19:169-210.

Marshall, L.G., Hoffster, R., and Pascual, R. 1983. Mammals and stratigraphy: geochronology of the mammal-bearing Tertiary of South America. Palaeovertebrata Mémoir Extraordinaire, 1-93.

Matheos, S.D. Raigemborn, M.S. 2012. Sedimentology and paleoenvironment of the Santa Cruz Formation, p. 59-82. In Vizcaíno, S.F., Kay, R.F, and Bargo, M.S. (eds.), Early Miocene Paleobiology in Patagonia: High-Latitude Paleocommunities of the Santa Cruz Formation, Cambridge University Press, Cambridge. https://doi.org/10.1017/cbo9780511667381.005

Matter, A., Homewood, P.W., Caron, C., Van Stuijvenberg, J., Weidmann, M., and Winkler, W. 1980. Flysch and molasse of central and western Switzerland, p. 261-293. In Trümpy, R. (ed.), Geology of Switzerland, A Guide Book. Schweiz Geologica Kommissione, Basel, New York.

McGrath, A., Anaya, F., and Croft, D.A. 2018. Two new macraucheniids (Mammalia: Litopterna) from the late middle Miocene (Laventan SALMA) of Quebrada Honda, Bolivia. Journal of Vertebrate Paleontology, 38:e1461632. https://doi.org/10.1080/02724634.2018.1461632

McNaughton, S.J. and Georgiadis, N.J. 1986. Ecology of African grazing and browsing mammals. Annual Review of Ecology and Systematics, 17:39-65. https://doi.org/10.1146/annurev.ecolsys.17.1.39

Mendoza, M., Goodwin, B., and Criado, C. 2004. Emergence of community structure in terrestrial mammal-dominated ecosystems. Journal of Theoretical Biology, 230:203-214. https://doi.org/10.1016/j.jtbi.2004.05.002

Murdock, T.Q., Weaver, A J., and Fanning, A.F. 1997. Paleoclimatic response of the closing of the Isthmus of Panama in a coupled ocean‐atmosphere model. Geophysical Research Letters, 24:253-256. https://doi.org/10.1029/96gl03950

Nesbit Evans, E.M., Van Couvering, J.A.H., and Andrews, P. 1981. Paleoecology of Miocene sites in western Kenya. Journal of Human Evolution, 10:99-116.

Nichols, G. 2009. Sedimentology and Stratigraphy. John Wiley and Sons, Chichester, UK.

Nicholson, S.E., Davenport, M.L., and Malo, A.R. 1990. A comparison of the vegetation response to rainfall in the Sahel and East Africa, using Normalized Difference Vegetation Index from NOAA AVHRR. Climate Change, 17:209-241. https://doi.org/10.1007/bf00138369

North, M., Collins, B.M., Safford, H., and Stephenson, N.L. 2016. Montane forests, p. 553-577. In Mooney, H. and Zavaleta, E. (eds.), Ecosystems of California, University of California Press, Berkeley.

(The) NOW Community. 2017. New and Old Worlds Database of Fossil Mammals (NOW). Retrieved November 01, 2017 from http://www.helsinki.fi/science/now/.

Nowak, R.M. 1999. Walker’s Mammals of the World, Volume 2 (sixth edition). John Hopkins University Press, London, United Kingdom.

Ogg, J.G. 2012. Geomagnetic Polarity Time Scale, p. 85-113. In Gradstein, F.M., Ogg, J.G., Schmitz, M., and Ogg, G. (eds.), The Geologic Time Scale 2012. Elsevier, Amsterdam.

Olson, D.M. and Dinerstein, E. 1998. The Global 200: a representation approach to conserving the Earth’s most biologically valuable ecoregions. Conservation Biology, 12:502-515. https://doi.org/10.1046/j.1523-1739.1998.012003502.x

Olson, D.M., Dinerstein, E., Wikramanayake, E.D., Burgess, N.D., Powell, G.V., Underwood, E.C., D’amico, J.A., Itoua, I., Strand, H.E., Morrison, J.C., and Loucks, C.J. 2001. Terrestrial ecoregions of the world: a new map of life on Earth: a new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity. BioScience, 51:933-938.

Ortiz-Jaureguizar, E. and Cladera, G.A. 2006. Paleoenvironmental evolution of southern South America during the Cenozoic. Journal of Arid Environments, 66:498-532. https://doi.org/10.1016/j.jaridenv.2006.01.007

Parada, A., Pardiñas, U.F., Salazar-Bravo, J., D’Elía, G., and Palma, R.E. 2013. Dating an impressive Neotropical radiation: molecular time estimates for the Sigmodontinae (Rodentia) provide insights into its historical biogeography. Molecular Phylogenetics and Evolution, 66:960-968. https://doi.org/10.1016/j.ympev.2012.12.001

Patterson, B. and Pascual, R. 1968. The fossil mammal fauna of South America. Quarterly Review of Biology, 43:409-451. https://doi.org/10.1086/405916

Peel, M.C., Finlayson, B.L., and McMahon, T.A. 2007. Updated world map of the Köppen-Geiger climate classification. Hydrology and Earth System Sciences Discussions, 4:439-473. https://doi.org/10.5194/hessd-4-439-2007

Peppe, D.J., Baumgartner, A., Flynn, A., and Blonder, B., 2018. Reconstructing paleoclimate and paleoecology using fossil leaves, p. 289-317. In Croft, D.A., Simpson, S.W., and Su, D.F (eds.), Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Springer (Vertebrate Paleobiology and Paleoanthropology Series), Cham, Switzerland. https://doi.org/10.31233/osf.io/stzuc

Peñuelas, J. and Boada, M. 2003. A global change‐induced biome shift in the Montseny mountains (NE Spain). Global Change Biology, 9:131- 140. https://doi.org/10.1046/j.1365-2486.2003.00566.x

Pérez-Barbería, F.J., Gordon, I.J., and Nores, C., 2001. Evolutionary transitions among feeding styles and habitats in ungulates. Evolutionary Ecology Research, 3:221-230.

Perkins, M.E., Fleagle, J.G., Heizler, M.T., Nash, B., Bown, T.M., Tauber, A., and Dozo, M.T. 2012. Tephrochronology of the Miocene Santa Cruz and Pinturas Formations, Argentina, p. 23-40. In Vizcaíno, S.F., Kay, R.F, and Bargo, M.S. (eds.), Early Miocene Paleobiology in Patagonia: High-Latitude Paleocommunities of the Santa Cruz Formation, Cambridge University Press, Cambridge. https://doi.org/10.1017/cbo9780511667381.003

Piperno, D.R. 2008. Paleoenvironmental reconstruction in the lowland Neotropics, p. 1772-1778. In Pearsall, D.M. (ed.), Encyclopedia of Archaeology. Elsevier, Oxford. https://doi.org/10.1016/b978-012373962-9.00227-2

Poux, C., Chevret, P., Huchon, D., De Jong, W.W., and Douzery, E.J. 2006. Arrival and diversification of caviomorph rodents and platyrrhine primates in South America. Systematic Biology, 55: 228-244. https://doi.org/10.1080/10635150500481390

Raigemborna, M.S., Krapovickas, V., Zucol, A.F., Zapata, L., Beilinson, E., Toledo, N., Perry, J., Lizzoli, S., Martegani, L., Tineo, D.E., and Passeggi, E. 2019. Paleosols and related soil-biota of the early Miocene Santa Cruz Formation (Austral-Magallanes Basin, Argentina): a multidisciplinary approach to reconstructing ancient terrestrial landscapes. Latin American Journal of Sedimentology and Basin Analysis, 25:117-148.

Raymo, M.E., Ruddiman, W.F., Backman, J., Clement, B.M., and Martinson, D.G. 1989. Late Pliocene variation in Northern Hemisphere ice sheets and North Atlantic deep water circulation. Paleoceanography, 4:413-446. https://doi.org/10.1029/pa004i004p00413

Reed, K.E. 1998. Using large mammal communities to examine ecological and taxonomic structure and predict vegetation in extant and extinct assemblages. Paleobiology, 24:384-408.

Reid, W.V., Mooney, H.A., Cropper, A., Capistrano, D., Carpenter, S.R., Chopra, K., Dasgupta, P., Dietz, T., Duraiappah, A.K., Hassan, R., Kasperson, R., Leemans, R., May, R.M., McMichael, A.J., Pingali, P., Samper, C., Scholes, R., Watson, R.T., Zakri, A.H., Shidong, Z., Ash, N.J., Bennett, E., Kumar, P., Lee, M.J., Raudsepp-Hearne, C., Simons, H., Thonell, J., and Zurek, M.B. 2005. Ecosystems and Human Well-being - Synthesis: A Report of the Millennium Ecosystem Assessment. Island Press, Washington, D.C.

Reguero, M.A., Candela, A.M., and Cassini, G.H. 2010. Hypsodonty and body size in rodent-like notoungulates, p. 362-374. In Madden, R.H., Carlini, A.A., Vucetich, M.G., and Kay. R.F. (eds.), The Paleontology of Gran Barranca: Evolution and Environmental Change through the Middle Cenozoic of Patagonia. Cambridge University Press, New York.

Reguero, M., Croft, D.A., Flynn, J.J., and Wyss, A.R. 2003. Small archaeohyracids (Notoungulata, Typotheria) from Chubut, Argentina and central Chile: trans-Andean temporal correlation. Fieldiana: Geology (N.S.), 48:1-17. https://doi.org/10.5962/bhl.title.5341

Retallack, G.J. 1984. Trace fossils of burrowing beetles and bees in an Oligocene paleosol, Badlands National Park, South Dakota. Journal of Paleontology, 58:571-592. 

Retallack, G.J. 1994. The environmental factor approach to the interpretation of paleosols. SSSA Special Publication, 33:31-31.

Retallack, G.J. 2001. Soils of the Past: An Introduction to Paleopedology, (second edition). Blackwell Science Ltd., Oxford. https://doi.org/10.1002/9780470698716

Retallack, G.J. 2005. Pedogenic carbonate proxies for amount and seasonality of precipitation in paleosols. Geology, 33:333-336. https://doi.org/10.1130/g21263.1

Rhoads, D.C. 1975. The paleoecological and environmental significance of trace fossils, p. 147-160. In Frey, R.W. (ed.), The Study of Trace Fossils. Springer-Verlag, New York. https://doi.org/10.1007/978-3-642-65923-2_9

Rowan, J., Beaudrot, L., Franklin, J., Reed, K.E., Smail, I.E., Zamora, A., and Kamilar, J.M. 2019. Geographically divergent evolutionary and ecological legacies shape mammal biodiversity in the global tropics and subtropics. Proceedings of the National Academy of Sciences, 201910489. https://doi.org/10.1073/pnas.1910489116

Sandom, C., Faurby, S., Sandel, B., and Svenning, J.C. 2014. Global late Quaternary megafauna extinctions linked to humans, not climate change. Proceedings of the Royal Society B: Biological Sciences, 281:20133254. https://doi.org/10.1098/rspb.2013.3254

Saunders, M.B. and Barclay, M.R. 1992. Ecomorphology of insectivorous bats: a test of predictions using two morphologically similar species. Ecology, 73:1335-1345. https://doi.org/10.2307/1940680

Scott, K.M. 1990. Postcranial dimensions of ungulates as predictors of body mass, p. 301-335. In Damuth, J. and MacFadden, B.J. (eds.), Body Size in Mammalian Paleobiology. Cambridge University Press, Cambridge.

Semken H.A., Jr., Graham, R.W., and Stafford, T.W., Jr. 2010. AMS 14C analysis of Late Pleistocene non-analog faunal components from 21 cave deposits in southeastern North America. Quaternary International, 217:240-55. https://doi.org/10.1016/j.quaint.2009.11.031

Shirkhorshidi, A.S., Aghabozorgi, S., and Wah, T.Y. 2015. A comparison study on similarity and dissimilarity measures in clustering continuous data. PloS One, 10:e0144059. https://doi.org/10.1371/journal.pone.0144059

Shuman, B.N. 2009. Approaches to paleoclimate reconstruction, p. 179-184. In Elias, S.A. and Mock, C.J. (eds.) Encyclopedia of Quaternary Science (second edition). Elsevier, Amsterdam. https://doi.org/10.1016/b0-44-452747-8/00008-9

Shmida, A. and Wilson, M.V. 1985. Biological determinants of species diversity. Journal of Biogeography, 12:1-20. https://doi.org/10.2307/2845026

Simpson, G.G. 1980. Splendid Isolation: The Curious History of South American Mammals. Yale University Press, New Haven, Connecticut.

Sinclair, E.H. and Zeppelin, T.K. 2002. Seasonal and spatial differences in diet in the western stock of Steller sea lions (Eumetopias jubatus). Journal of Mammology, 83:973-990. https://doi.org/10.1644/1545-1542(2002)083<0973:sasdid>2.0.co;2

Smith, F.A., Smith, R.E.E., Lyons, S.K., and Payne, J.L. 2018. Body size downgrading of mammals over the late Quaternary. Science, 360: 310-313. https://doi.org/10.1126/science.aao5987

Smith, J.J. and Hasiotis, S.T. 2008. Traces and burrowing behaviors of the cicada nymph Cicadetta calliope: neoichnology and paleoecological significance of extant soil-dwelling insects. Palaios, 23:503-513. https://doi.org/10.2110/palo.2007.p07-063r

Smith, J.R., Letten, A.D., Ke, Po-Ju., Anderson, C.B., Hendershot, J.N., Dhami, M.K., Dloot, G.A., Grainger, T.N., Howard, M.E., Morrison, B.M.L., Routh, D., San Juan, P.A., Mooney, H.A., Mordecai, E.A., Crowther, T.W., and Daily, G.C. 2018. A global test of ecoretions. Nature Ecology and Evolution, 2:1889-1896. https://doi.org/10.1038/s41559-018-0709-x

Soligo, C. and Andrews P. 2005. Taphonomic bias, taxonomic bias and historical non-equivalence of faunal structure in early hominin localities. Journal of Human Evolution, 49:2006-229. https://doi.org/10.1016/j.jhevol.2005.03.006

Solomon, A.M. and Silkworth, A.B. 1986. Spatial patterns of atmospheric pollen transport in a montane region. Quaternary Research, 25:150-162. https://doi.org/10.1016/0033-5894(86)90053-0

Spradley, J.P., Glazer, B.J., and Kay, R.F. 2019. Mammalian faunas, ecological indices, and machine-learning regression for the purpose of paleoenvironmental reconstruction in the Miocene of South American. Palaeogeography, Palaeoclimatology, Palaeoecology, 518:155-171. https://doi.org/10.1016/j.palaeo.2019.01.014

Stanley, S.M. 2005. Earth System History. W. H. Freeman and Company, New York.

Stehli F.G. and Webb, S.D. 1985. The Great American Biotic Interchange. Topics in Geobiology, Plenum Press, New York

Stoms, D.M. 1994. Scale dependence of species richness maps. Professional Geographer, 46:346-358. https://doi.org/10.1111/j.0033-0124.1994.00346.x

Stoms, D.M. and Estes, J.E. 1992. A remote sensing research agenda for mapping and monitoring biodiversity. International Journal of Remote Sensing, 10:1839-1860. https://doi.org/10.1080/01431169308954007

Strömberg, C.A., Dunn, R.E., Crifò, C., and Harris, E.B. 2018. Phytoliths in paleoecology: analytical considerations, current use, and future directions, p. 235-287. In Croft, D.A., Simpson, S.W., and Su, D.F (eds.), Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Springer (Vertebrate Paleobiology and Paleoanthropology Series), Cham, Switzerland. https://doi.org/10.1007/978-3-319-94265-0_12

Su, D.F. and Croft, D.A. 2018. Making sense of the evidence: synthesizing paleoecological data, p. 395-404. In, Croft, D.A., Su, D.F., and Simpson, S.W. (eds.), Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Cham, Switzerland.  https://doi.org/10.1007/978-3-319-94265-0_18

Toledo, N., Cassini, G. H., Vizcaíno, S.F., and Bargo, M.S. 2014. Mass estimation of Santacrucian sloths from the early Miocene Santa Cruz formation of Patagonia, Argentina. Acta Palaeontologica Polonica, 59:267-280.

Udvardy M.D.F. 1975. A classification of the biogeographical provinces of the world. Morges (Switzerland): International Union of Conservation of Nature and Natural Resources. IUCN Occasional Paper, 18. https://doi.org/10.1017/s0376892900019226

Van Couvering, J.A.H. 1980. Community evolution in East Africa during the late Cenozoic p. 272-298. In Behrensmeyer, A.K. and Hill, A.P. (eds.), Fossils in the Making: Vertebrate Taphonomy and Paleoecology. The University of Chicago Press, Chicago.

Van der Meulen, A.J. and Daams, R. 1992. Evolution of early-middle Miocene rodent faunas in relation to long-term palaeoenvironmental changes. Palaeogeography, Palaeoclimatology, Palaeoecology, 93:227-253.

Vermillion, W.A., Polly, P.D., Head, J.J., Eronen, J.T., and Lawing, A.W. 2018. Ecometrics: a trait-based approach to paleoclimate and paleoenvironmental reconstruction, p. 373-394. In Croft, D.A., Simpson, S.W., and Su, D.F (eds.), Methods in Paleoecology: Reconstructing Cenozoic Terrestrial Environments and Ecological Communities. Springer (Vertebrate Paleobiology and Paleoanthropology Series), Cham, Switzerland. https://doi.org/10.1007/978-3-319-94265-0_17

Villavicencio, N.A., Lindsey, E.L., Martin, F.M., Borrero, L.A., Moreno, P.I., Marshall, C.R., and Barnosky, A.D. 2016. Combination of humans, climate, and vegetation change triggered Late Quaternary megafauna extinction in the Última Esperanza region, southern Patagonia, Chile. Ecography, 39:125-140. https://doi.org/10.1111/ecog.01606

Vizcaíno, S.F., Bargo, M.S., and Cassini, G.H. 2006. Dental occlusal surface area in relation to body mass, food habits and other biological features in fossil xenarthrans. Ameghiniana, 43:11-26.

Vizcaíno, S.F., Kay, R.F., and Bargo, M.S. 2012. Early Miocene Paleobiology in Patagonia: High-Latitude Paleocommunities of the Santa Cruz Formation. Cambridge University Press, Cambridge. https://doi.org/10.1017/cbo9780511667381

Voorhies, M.R. 1975. Vertebrate burrows, p. 325-350. In Frey, R.W. (ed.), The Study of Trace Fossils. Springer-Verlag, New York. https://doi.org/10.1007/978-3-642-65923-2_15

Vucetich, M.G., Vieytes, E.C., Pérez, M.E., and Carlini, A.A. 2010. The rodents from La Cantera and the early evolutionof caviomorphs in South America, p. 189-201. In Madden, R.H., Carlini, A.A., Vucetich, M.G., and Kay. R.F. (eds.), The Paleontology of Gran Barranca: Evolution and Environmental Change through the Middle Cenozoic of Patagonia. Cambridge University Press, New York.

Wainwright, P.C. and Reillly, S.M. 1994. Ecological Morphology: Integrative Organismal Biology. University of Chicago Press, Chicago.

Walton, A.H. 1997. Rodents, p. 392-409. In Kay, R.F., Madden, R.H., Cifelli, R.H., and Flynn, J.J. (eds.), Vertebrate Paleontology in the Neotropics: The Miocene Fauna of La Venta, Colombia, Smithsonian Institution Press, Washington, D.C.

Webb, S.D. 1976. Mammalian faunal dynamics of the great American interchange. Paleobiology, 4:27-40. https://doi.org/10.1017/s0094837300004802

Webb, S.D. 1985. Faunal interchange between North and South America. Acta Zoologica Fennica, 170:177-178.

Webb, S.D. 2006. The great American biotic interchange: patterns and processes. Annals of the Missouri Botanical Garden, 93:245-257.

Werdelin, L. 1996. Carnivoran ecomorphology: a phylogenetic perspective, p. 582-624. In Gittleman, J.L. (ed.), Carnivore Behavior, Ecology and Evolution, Volume 2. Cornell University Press, Ithaca.

Whittaker, R.H. 1971. Evolution and measurement of species diversity. Taxon, 21:215-251. https://doi.org/10.2307/1218190

Whittaker, R.H. 1977. Evolution of species diversity in land communities. Evolutionary Biology, 10:1-67. 

Williams, H.S.H., and Kay, R. 2001. A comparative test of adaptive explanations for hypsodonty in ungulates and rodents. Journal of Mammalian Evolution, 8:207-229.

Wing, S.L. and Greenwood, D.R. 1994. Fossils and fossil climate: the case for equable continental interiors in the Eocene. p. 243-252. In Allen, J.R.L., Hoskins, B.K., Sellwood, B.W., and Spicer, R.A. (eds.), Paleoclimates and Their Modelling with Special Reference to the Mesozoic Era. Springer, Dordrecht. https://doi.org/10.1007/978-94-011-1254-3_5

Wolff, R.G. 1975. Sampling and sample size in ecological analyses of fossil mammals. Paleobiology, 1:195-204. https://doi.org/10.1017/s0094837300002384

Woodburne, M.O. 2010. The Great American Biotic Interchange: dispersals, tectonics, climate, sea level and holding pens. Journal of Mammal Evolution, 17:245-264. https://doi.org/10.1007/s10914-010-9144-8

Woodburne, M.O., Cione, A.L., and Tonni, E.P. 2006. Central American provincialism and the great American biotic interchange, p. 73-101. In Carranza-Castañeda, O. and Lindsay E.H. (eds.), Advances in Late Tertiary Vertebrate Paleontology in Mexico and the Great American Biotic Interchange. Publicación Especial del Instituto de Geología y Centro de Geociencias de la Universidad Nacional Autónoma de México, Mexico.

World Wildlife Fund. 2006. WildFinder: Online Database of Species Distributions, accessed January 2006. https://www.worldwildlife.org/publications/wildfinder-database

Wright, D.K. 2017. Paleoenvironmental science: methods. Oxford Research Encyclopedia of African History. Downloaded 14 March 2019. https://doi.org/10.1093/acrefore/9780190277734.013.211

Wyss, A.R., Flynn, J.J., and Croft, D.A. 2018. New Paleogene notohippids and leontiniids (Toxodontia: Notoungulata; Mammalia) from the early Oligocene Tinguiririca Fauna of the Andean Main Range, Central Chile. American Museum Novitates, 3903:1-42. https://doi.org/10.1206/3903.1

Wyss, A.R., Flynn, J.J., Norell, M., Swisher, C.C., Novacek, M.J., McKenna, M.C., and Charrier, R., 1994. Paleogene mammals from the Andes of central Chile: a preliminary taxonomic, biostratigraphic, and geochronologic assessment. American Museum Novitates, 3098:1-31.

Wyss, A.R., Norell, M.A., Flynn, J.J., Novacek, M.J., Charrier, R., McKenna, M.C., Swisher, C.C., Frasinett, D., and Jin, M. 1990. A new early Tertiary mammal fauna from central Chile: implications for Andean stratigraphy and tectonics. Journal of Vertebrate Paleontology, 10:518-522. https://doi.org/10.1080/02724634.1990.10011835

Zimicz, N. 2011. Patrones de desgaste y oclusión en el sistema masticatorio de los extintos Argyrolagoidea (Marsupialia, Polydolopimorphia, Bonapartheriiformes). Ameghiniana, 48:358-379. https://doi.org/10.5710/amgh.v48i3(472)

 

 

 

Related Articles