Open Access
Translator Disclaimer
5 March 2021 The effects of geographic range size and abundance on extinction during a time of “sluggish”' evolution
Michelle M. Casey, Erin E. Saupe, Bruce S. Lieberman
Author Affiliations +

Geographic range size and abundance are important determinants of extinction risk in fossil and extant taxa. However, the relationship between these variables and extinction risk has not been tested extensively during evolutionarily “quiescent” times of low extinction and speciation in the fossil record. Here we examine the influence of geographic range size and abundance on extinction risk during the late Paleozoic (Mississippian–Permian), a time of “sluggish” evolution when global rates of origination and extinction were roughly half those of other Paleozoic intervals. Analyses used spatiotemporal occurrences for 164 brachiopod species from the North American midcontinent. We found abundance to be a better predictor of extinction risk than measures of geographic range size. Moreover, species exhibited reductions in abundance before their extinction but did not display contractions in geographic range size. The weak relationship between geographic range size and extinction in this time and place may reflect the relative preponderance of larger-ranged taxa combined with the physiographic conditions of the region that allowed for easy habitat tracking that dampened both extinction and speciation. These conditions led to a prolonged period (19–25 Myr) during which standard macroevolutionary rules did not apply.


Determining the correlates of extinction is crucial to understanding macroevolutionary processes operating on geologic timescales (Jablonski 1986; McKinney 1997; Kiessling and Aberhan 2007; Payne and Finnegan 2007; Meseguer et al. 2015; Saupe et al. 2015) and for identifying taxa at potential risk of extinction today (Lee and Jetz 2011; Finnegan et al. 2015; Kiessling et al. 2019; Smits and Finnegan 2019). The International Union for Conservation of Nature Red List, an indicator of the global conservation status of biological species, uses both abundance and geographic range size as determinants of extinction risk for modern species (IUCN 2001). Abundance and geographic range size are often positively correlated (Gaston 1994; Gaston et al. 1997; Holt et al. 1997; Harnik 2011) but are not always equally important determinants of extinction risk (e.g., Kiessling and Aberhan 2007; Payne et al. 2011; Harnik et al. 2012). For example, examining Neogene bivalves in the Pacific, Stanley (1986) found that abundance, and not geographic range size, was a strong predictor of extinction selectivity. Similarly, Payne et al. (2011) showed an inverse association between abundance and extinction risk that was largely independent of geographic range size. Others, however, have found geographic range size to be a strong predictor of extinction risk, with this relationship not explained by the positive correlation between abundance and geographic range size (e.g., Kiessling and Aberhan 2007; Harnik et al. 2012).

Univariate analyses cannot glean the relative importance of abundance and geographic range size as predictors of extinction risk, because these variables are measured in different units and are rarely evaluated together in a multivariate framework that adequately establishes their relative importance (Harnik et al. 2012). Harnik et al. (2012) is an exception, in which the authors analyzed the association between extinction risk and traits related to taxon rarity (geographic range size, habitat breadth, and local abundance) using a multivariate framework and a global database of Phanerozoic marine fossil genera. The authors found that geographic range size was a strong predictor of extinction risk, whereas abundance had little effect. However, the authors were unable to calculate the relationship between geographic range size and extinction risk during the late Carboniferous and early Permian due to low numbers of generic extinctions during this interval. Consequently, here we assemble spatiotemporal fossil occurrence data to fill this knowledge gap and evaluate the relative strength of abundance and geographic range size as correlates of extinction during the late Paleozoic ice age (LPIA) using brachiopods from the midcontinent of the United States.

The late Paleozoic is interesting from a macroevolutionary perspective, because it was a time of “sluggish” evolution when global rates of origination and extinction were low (Stanley and Powell 2003; Powell 2005; Segessenman and Kammer 2018; Balseiro and Halpern 2019; Kolis and Lieberman 2019) in spite of dramatic, cyclical changes in climate and the environment (Parrish 1993; Olszewski and Patzkowsky 2001a,b; Raymond and Metz 2004; Horton et al. 2012; Balseiro 2016). The LPIA was the longest-lasting glacial period of the Phanerozoic, which lasted from 320 to 260 Ma. During the LPIA, Southern Hemisphere glacial–interglacial cycles were governed by Milankovitch orbital cycles (Raymond and Metz 2004). The resulting glacioeustatic sea-level changes produced cyclothems (Montañez and Poulsen 2013) characterized by vacillation of deep-marine, shallow-marine, and terrestrial environments (Heckel et al. 1994; Heckel 2008).

Late Paleozoic brachiopods within the midcontinent of the United States show little taxonomic turnover (Olszewski and Patzkowsky 2001a,b), making them excellent targets for a study of macroevolutionary dynamics during this “quiescent” interval. Brachiopods are also extremely abundant and typically well preserved (Foote and Sepkoski 1999). Recent museum digitization efforts and ongoing initiatives, such as those funded by the Advancing the Digitization of Biodiversity Collections (ADBC) directorate of the U.S. National Science Foundation (NSF), allow for assembly of brachiopod datasets from which geographic range and abundance data can be derived (Page et al. 2015). Moreover, extinction and origination rates for brachiopods at a global scale match broadly those observed in a range of marine shelled invertebrate taxa (Powell 2005).

The midcontinent of North America is an excellent target region to study the late Paleozoic, because it contains the best-studied record of the greenhouse/icehouse transitions during the Carboniferous and Permian (Heckel 1977, 1986; Heckel et al. 1994; Lupia and Armitage 2013). The stratigraphy of the region is well constrained and, in many cases, correlated across state lines (Heckel 1986; Olszewski and Patzkowsky 2003). The North American midcontinent was located near the equator for much of the Carboniferous, where changes associated with glaciation and global cooling were likely pronounced (Algeo and Heckel 2008). Previous work has shown that cooling temperatures and increased seasonality, which marked the onset of glaciation during the Mississippian ca. 323 Ma (Parrish 1993), led to preferential loss of tropical brachiopod species with narrow latitudinal ranges (Powell 2005). According to Powell (2005), brachiopod survivors had broader latitudinal distributions, potentially reflecting broader thermal tolerances, longer taxonomic durations, and larger populations. The absence of narrowly distributed and quickly evolving genera from the tropics led to a change in macroevolutionary dynamics, in which broadly adapted genera with long stratigraphic durations were common globally, not only in extratropical latitudes (Brezinski 1988; Stanley and Powell 2003; Powell 2005; Segassenman and Kammer 2018).

Powell (2005) proposed that the low extinction and origination rates during the Late Mississippian to middle Permian were due, in part, to the higher proportion of brachiopod taxa with larger latitudinal ranges, potentially reflecting broader ecological tolerances. The same characteristics that typically buffer against extinction (large geographic range size; broad ecological tolerance; and abundant, stable populations) also inhibit isolation of populations that would lead to speciation, potentially impeding generation of new taxa, themselves usually characterized by small geographic distributions at higher risk of extinction (Vrba 1980; Jablonski 1986; Stanley 1986; Foote 2007; Antell et al. 2020). Thus, reduced variation and skew toward larger range sizes may have rendered this variable an ineffective predictor of extinction risk during the Late Mississippian to middle Permian. Under these circumstances, abundance, rather than geographic range size, may be a better predictor of extinction risk and macroevolutionary dynamics during the LPIA. Although geographic range size and abundance frequently covary, they do not exhibit a one-to-one relationship and have been shown to predict extinction risk to varying degrees (e.g., Stanley 1986; Kiessling and Aberhan 2007; Payne et al. 2011; Harnik et al. 2012). We hypothesize that abundance is a better predictor of brachiopod extinction risk in the midcontinent and test this hypothesis using generalized linear models and 32,766 spatiotemporal occurrences over nine stages from the Chesterian to Leonardian (encompassing the LPIA).

Materials and Methods

Spatiotemporal Occurrence Data.—Individual, global specimen occurrences for brachiopod species present in the Carboniferous–Permian of the midcontinent of the United States were obtained from multiple, spatially explicit databases, including the Division of Invertebrate Paleontology, Biodiversity Institute (KUMIP), the Yale University Peabody Museum of Natural History (YPM), and the Paleobiology Database (PBDB). A total of 32,766 specimen occurrence records were obtained from existing databases and digitization efforts, comprising 4998 records from the PBDB and 27,768 records from the KUMIP and YPM. Only midcontinental species were targeted. However, the entire geographic range of these midcontinental species was reconstructed globally to capture true extinctions rather than local extirpations. Specimen records from KUMIP and YPM were chosen because these institutions have large numbers of brachiopod specimens from the midcontinent, with a high degree of stratigraphic and geographic control. We retained only occurrences with geographic uncertainty radii under 50 km. Records that lacked species-level identifications (those including sp., cf., aff., or ? in the species designation) or those that did not have stratigraphic information to the level of formation were removed. Taxonomy was standardized and updated using Muir-Wood and Copper (1960), Hoare (1961), Moore (1964), Williams et al. (1965), and Carter and Carter (1970), which represent key references on brachiopods from this region and time interval.

To ensure that distributional data were derived from geologic units of similar ages, a stratigraphic database for the midcontinent was generated from an extensive survey of the primary literature (see Supplementary Material). Informal North American stages (e.g., Kinderhookian–Guadalupian; Heckel and Clayton 2006; Menning et al. 2006) were followed to allow for the highest temporal resolution while maintaining the greatest sample size. Occurrences that could not be classified confidently to a North American stage were removed from analysis. Origination and extinction rates in each LPIA North American stage (Chesterian–Leonardian) were calculated using the modified gap-filler method of Alroy (2015) in R v. 3.4.1 (R Core Team 2017) using the divDyn v. 0.8.0 package (Kocsis et al. 2019). The modified gap-filler method of Alroy (2015) was used to minimize potential bias (e.g., those associated with the Signor-Lipps effect) while maintaining precision and accuracy. Rates from stages with <100 occurrences were not recorded but were used in calculations. To facilitate comparison with published species-level Paleozoic rates (e.g., Stigall 2010; Kolis and Lieberman 2019), we additionally calculated extinction and speciation rates as per-million-year rates using the equations of Foote (2000).

We assigned each species a chronostratigraphic duration based on the longest combination of (1) the chronostratigraphic duration reported in the PBDB; (2) the lithostratigraphic units the species occurred in according to Carter and Carter (1970); and (3) the lithostratigraphic units containing the species within museum databases. The chronostratigraphic ages of the lithostratigraphic units listed in these sources were determined using the compiled stratigraphic literature (see Supplementary Material) and the National Geologic Map Database (or Geolex; U.S. Geological Survey. n.d.). These chronostratigraphic durations were consulted when calculating per-million-year extinction and speciation rates, such that species with gaps in their records were never spuriously counted as extinctions. We tested the correlation between extinction and speciation rates using a generalized linear model (GLM) with a Gaussian distribution in R v. 3.4.1. To test for the effect of differential preservation on macroevolutionary rates, speciation and extinction rates were correlated with sampling intensity (estimated as the total number of brachiopod occurrences per stage; Powell 2005) using a GLM with a Gaussian distribution in R v. 3.4.1. A positive correlation between evolutionary rates and sampling intensity could indicate that macroevolutionary rates were driven by differences in sampling intensity rather than biological patterns (Powell 2005).

We ensured that paleogeographic analyses were conducted on spatially unique occurrences by culling each species to a single occurrence within an occupied grid cell of 0.025° × 0.025° resolution for each stage (equivalent to ∼3 km at the equator). This procedure removed artificial inflation of spatially unique occurrences caused by differences in georeferencing protocols among institutions and individuals; for example, two museum workflows yielding different decimal degree latitude and longitude estimates for the same collection locality. We removed singletons from analyses (i.e., any species by time bin with n = 1 spatially unique occurrence) to limit the effect of poorly sampled taxa. After removal of singletons, there remained 29,720 individual occurrences (4915 of which were spatially unique) belonging to 164 species in 83 genera.

The resulting species by time bin datasets were imported into ArcGIS v. 10.5.1, and the present-day latitude/longitude coordinates were rotated to their paleo-position using the University of Texas Institute for Geophysics (UTIG) plate model in PaleoWeb v. 1.0 (Rothwell Group Inc). Paleo-coordinate reconstructions were performed based on the start of the appropriate stage in millions of years before present (Fig. 1). Due to the equatorial position of Laurasia in the late Paleozoic, we applied the South American Albers Equal Area Conic map projection to reconstructed paleolatitude and paleolongitude occurrences before range-size metrics were calculated. The South American projection was chosen given the preponderance of coordinates located in the region occupied by present-day Brazil once they were rotated to their paleolatitude and paleolongitude.

Geographic Range-Size Metrics.—We quantified geographic range size and abundance for each unique species by stage. Geographic range was quantified using two different metrics: minimum convex hull and latitudinal range (for an example, see Fig. 1). A convex hull is a two-dimensional metric that estimates the amount of area occupied by a taxon and was calculated as the median of a series of convex hulls produced by exhaustive jackknifing of individual occurrence points (e.g., Stigall and Lieberman 2006; Hendricks et al. 2008; Myers and Lieberman 2011; Darroch and Wagner 2015; Saupe et al. 2015; Darroch et al. 2020). This method has been shown to be especially efficacious for quantifying ranges of fossil taxa (Darroch and Saupe 2018; Darroch et al. 2020). If a species was characterized by only two spatially unique occurrences, a 10 km buffer was applied to each occurrence, and the area of the resulting line substituted for the convex hull (Hendricks et al. 2008; Myers and Lieberman 2011). Each convex-hull geographic range estimate was log transformed for normality.

Latitudinal range (maximum observed paleolatitude minus minimum observed paleolatitude) is also used commonly to characterize geographic range sizes (Powell 2005, 2007; Foote and Miller 2013; Finnegan et al. 2016; Balseiro and Halpern 2019; Darroch et al. 2020). Unlike a convex hull, latitudinal range is a linear metric and may reflect breadth of thermal tolerance (Jackson 1974; Stanley and Powell 2003; Powell 2007; Sunday et al. 2012), with larger latitudinal ranges potentially indicating greater thermal tolerances. Latitudinal range estimates were square-root transformed for normality.


Paleogeographic reconstruction of spatially unique occurrence locations created using the PALEOMAP Paleo Atlas for GPlates v. 3 (Scotese 2016) for illustration purposes only. We used the University of Texas Institute for Geophysics (UTIG) plate model in PaleoWeb 1.0 for our analyses. Geographic range size was measured as a latitudinal range (yellow line with brackets) and as a convex hull (yellow polygon). Analyses used the median convex-hull area of all possible convex hulls created by jackknifing occurrences, which reduced the impact of geographic outliers on geographic range size. A, Desmoinesia muricatina (n = 118 spatially unique points) from the Desmoinesian. B, Neochonetes transversalis (n = 343 spatially unique points) from the Virgilian.


Abundance was quantified as the number of occurrences for a species within a stage. Calculating population size is not easy, even for modern organisms (He and Gaston 2000). However, occurrence data have been shown to provide adequate estimates of average abundance for both fossil and modern organisms (Buzas et al. 1982; Kunin 1998; Alroy 2000; He and Gaston 2003). Abundance estimates were square-root transformed for normality.

All three variables (convex hull, latitudinal range, and abundance) were calculated in R v. 3.4.1 (R Core Team 2017) using the packages sp v. 1.3-2 (Pebesma and Bivand 2005; Bivand et al. 2013) and PBSmapping v. 2.72.1 (Schnute et al. 2013).

Mixed-Effects Models.—We explored the nature of the relationship between species' extinction and abundance/range size metrics using generalized linear mixed-effects models (LMM) using both transformed (log or square root) unstandardized variables and standardized variables. For the latter, variables were standardized at stage level by dividing the value for each species by the maximum value of that variable in the stage (Foote et al. 2008). For example, the log geographic range size for Wellerella truncata during the Virgilian was divided by the maximum log-transformed geographic range size of any species from the Virgilian. Therefore, standardization scaled all three variables to vary between 0 and 1 in each stage. The results for the standardized variable analysis are included in the Supplementary Material (Supplementary Tables S3 and S4). For both unstandardized and standardized analyses, stages were limited to those associated with the LPIA (e.g., Chesterian–Leonardian). The time interval and the Linnaean family rank associated with each record were included as random effects in the LMM. We included these random effects because species from the same time bin and family were expected to pseudoreplicate one another due to shared environmental conditions/spatial sampling structure and evolutionary history, respectively. These random effects obviated the need to correct for multiple testing, because a single coefficient was calculated for each fixed effect over the entire time series.

We estimated the influence of the following three fixed effects in every combination: latitudinal range, abundance, and convex-hull area. This resulted in a total of eight possible models, including a model with no fixed effects added (i.e., a horizontal line). We implemented mixed-effects models with functions in the package lme4 v. 1.1-21 (Venables and Ripley 2002; Bates et al. 2015).

Akaike's information criterion (AIC) was estimated instead of corrected AIC, given the large number of species in the dataset, to determine the subset of models that received nontrivial weight (ΔAIC < 10; for derivation and discussion of weights, see Burnham and Anderson 2002). Models are listed in order of relative performance. For fixed effects of models, we recorded coefficient estimates on the logit scale and the 95% confidence interval bounds from the likelihood profile. Models were checked for overdispersion and linearity between continuous predictor variables and the logit of the outcome.

Temporal Analysis of Trends in Range Size and Abundance.—Foote (2007) suggested that genera tend to exhibit declining geographic range sizes preceding their extinction. Consequently, we evaluated temporal trends in geographic range size, latitudinal range, and abundance throughout a species' lifetime, including occurrence records from the Carboniferous to Permian (Kinderhookian–Lopingian). Unlike the generalized linear mixed effects framework, which accounted for shared environmental conditions and spatial sampling structure by specifying stage as a random effect, this analysis compared geographic range size and abundance estimates for a single species between stages characterized by differing area of rock outcrop. To account for effects of rock availability on our three metrics (geographic range size, latitudinal range, and abundance), variables were standardized at the stage level as described earlier (Foote et al. 2008). Temporal trends were evaluated using both standardized and unstandardized variables. Results for unstandardized variables are reported in the Supplementary Material (Supplementary Table S5).

Species present in three or more consecutive stages were determined to decline before extinction if they conformed to the following conditions: (1) the terminal value was lower than the value in the immediately preceding stage; and (2) the terminal value was lower than the mean for the species (Fig. 2). The number of species in decline was tabulated for each metric (convex hull, latitudinal range, and abundance). A one-tailed binomial test on the proportion of species in decline versus the number of stable and increasing species determined whether this proportion was significantly higher than 50%. That is, the test examined whether more species experienced declines in abundance or geographic range size before extinction than species that increased or remained unchanged; effect size was evaluated using Cohen's g.


Speciation and Extinction Rates.—Rates of speciation and extinction were low for most of the stages associated with the LPIA (Supplementary Table S1). For the Chesterian to Missourian stages, proportional extinction ranged between 0.014 and 0.154. Highest proportional extinction was in the Wolfcampian at 0.884, preceded by the second-highest proportional extinction of 0.363 in the Virgilian. Proportional extinction was lowest for the Missourian at 0.014. The median proportional extinction rate for the LPIA stages was 0.091. Extinction rates measured as per-million-year rates were lowest in the Morrowan (0.004) and highest in the Virgilian (0.047). As with proportional extinction, the highest per-million-year extinction rates were found in the Virgilian and Wolfcampian, but their relative ranks were reversed. For context, background extinction rates derived from species-level brachiopod data for the Middle Devonian range from 0.2 to 0.6 per million years (Stigall 2010), while late Paleozoic extinction rates derived from species-level cephalopod data range from 0.07 to 0.34 per million years (Kolis and Lieberman 2019). Thus, per-million-year extinction rates for all stages are below the range of Middle Devonian per-million-year background rates and comparable to or lower than other species-level late Paleozoic studies. Note that this dataset could not capture the Serupkhovian mass extinction (Chesterian), because we focused on LPIA taxa and omitted any taxa entirely absent from the Pennsylvanian or Permian.


Standardized geographic range size and abundance through time for the brachiopod species Cancrinella boonensis. The dashed horizontal line indicates the mean value for the species over this interval. A, Standardized geographic range size, measured as convex-hull area (km2) over maximum convex-hull area (km2) for the stage. Cancrinella boonensis specimen IP.008072 Yale Peabody Museum of Natural History; photo by J. Utrup, 2011. B, Standardized geographic range size, measured as species latitudinal range over maximum latitudinal range for the stage. C, Standardized abundance, measured as species abundance over maximum abundance for the stage. Abundance is the only metric that meets the criteria for decline before extinction—i.e., the terminal value is lower than the preceding value and below the species mean. Stage abbreviations: O, Osagean; Me, Meramecian; C, Chesterian; Mo, Morrowan; A, Atokan; D, Desmoinesian; M, Missourian; V, Virgilian; Wol, Wolfcampian.


Proportional speciation ranged from 0.154 to 0.486 and was always higher than proportional extinction, except during the Wolfcampian. The highest proportional speciation, 0.486, occurred during the Virgilian. Median proportional speciation rate for the LPIA was 0.173. Proportional speciation and extinction rates are not significantly correlated (GLM t = 0.66, p = 0.53). Per-million-year speciation rates ranged between 0.008 in the Wolfcampian and 0.054 in the Morrowan and are low relative to per-million-year extinction rates derived from species-level data on late Paleozoic cephalopods (0.03 to 0.41; Kolis and Lieberman 2019). Per-million-year speciation and extinction rates are not significantly correlated (GLM t = 0.11, p = 0.92).

Table 1.

Coefficient estimates and confidence intervals (CI) for mixed-effects models using unstandardized variables. Coefficient estimates (logit scale) and confidence intervals for each predictor in the eight models, with extinction status as response. Models were built on 377 range size/abundance records for 164 species from nine Phanerozoic stages.

Abundance and latitudinal range were square-root transformed; convex hull was log transformed for normality. The cumulative weight of the best-supported models was 99.7%. Confidence intervals, however, are wide, because there were relatively few extinctions even when pooled across intervals (88 extinctions compared with 289 survivals across the nine stages). Confidence intervals inclusive of zero are in bold.


To consider whether variations in extinction and speciation rate might be related to changes in sampling intensity, we examined the relationship between sampling intensity and speciation and extinction rates using GLMs (Powell 2005). No significant correlations were found between sampling intensity (estimated as the total number of brachiopod occurrences in each stage) and either extinction rate (GLM t = 0.41, p = 0.69) or speciation rate (GLM t = 1.23, p = 0.25) (Supplementary Table S2), indicating that patterns in origination and extinction likely reflect real biological fluctuations, not simply changes in preservation or sample availability. However, it is worth noting that edge effects are likely to inflate origination and extinction rates at both the start and end of the examined time interval (Foote 2000); extending the study interval to account for edge effects would likely reduce the rates for both the Chesterian and Wolfcampian reported herein.

Mixed-Effects Models.—The three assessed predictors of extinction risk (convex hull, latitudinal range, and abundance) were not highly collinear, based on the variance inflation factor (VIF). The VIF represents the proportion of variance in one predictor explained by all the other predictors (<2.5, with 1 being no collinearity; Zorro et al. 2010) and was 2.74, 2.18, and 1.41 for convex hull, latitudinal range, and abundance, respectively.

Mixed-effects models indicate that abundance is the strongest predictor of extinction risk. The three highest ranking mixed-effects models all included abundance as a predictor (Table 1). Moreover, abundance was the only predictor in the top three models to have confidence interval estimates exclusive of zero (Table 1) and was the sole predictor included in the best-supported model (Table 2). By contrast, range size–only models received little support (Table 1). Both latitudinal range and convex-hull area performed best when included as predictors with abundance (Tables 1 and 2, models ranked 2 and 3), but confidence intervals for these predictors were almost always inclusive of zero (Table 1). Of the best-supported models, only one model included measures of geographic range with confidence intervals that were exclusive of zero (model 5; convex hull as the only predictor). Results remain virtually unchanged when conducted using standardized variables (Supplementary Tables S3 and S4).

Table 2.

Relative performance of all mixed-effects models using unstandardized variables. Models predicted species extinction status based on measures of abundance and range size (latitudinal range and convex hull). Weights are the ratio of ΔAIC from a given model to the sum of ΔAIC values across all candidate models and are interpreted as the probability that a given model is the “best” (minimizes the Kullback–Leibler discrepancy) of the candidate models, given the data.


Temporal Analysis of Trends in Range Size and Abundance.—No statistical support was found for declining geographic range size in species before their extinction. In particular, the proportion of species with range size declines before extinction was no higher than the proportion of species without such declines, regardless of whether latitudinal range or convex hulls were used to quantify geographic range (Table 3). By contrast, species were found to decline in abundance in the interval leading to extinction significantly more than half the time (proportion of species displaying declining abundance = 62%, p = 0.05; Table 3), although Cohen's g suggests the effect size of this difference can be classified as small (g = 0.12). Results remain unchanged when conducted using unstandardized variables (Supplementary Table S5).


Late Paleozoic brachiopods from the North American midcontinent provide a diverse and abundant record of marine life during a distinctive time in Earth history associated with profound climatic oscillations. This time interval is especially noteworthy for the low extinction and speciation rates displayed in marine invertebrates (Stanley and Powell 2003; Powell 2005; Segessenman and Kammer 2018; Balseiro and Halpern 2019; Kolis and Lieberman 2019). Our analysis of macroevolutionary rates (Supplementary Table S1) is consistent with results from previous studies indicating this was a time of sluggish macroevolution (Stanley and Powell 2003; Powell 2005; Segessenman and Kammer 2018; Balseiro and Halpern 2019; Kolis and Lieberman 2019). In particular, per-million-year extinction and speciation rates for LPIA intervals are low (Supplementary Table S1), consistent with previously published work (Stanley and Powell 2003; Powell 2005). Per-million-year speciation rates are low throughout the interval, comparable to per-million-year speciation rates from cephalopods during the same interval (Kolis and Lieberman 2019). Similarly, the slight increase in extinction during the Late Pennsylvanian (Virgilian) and early Permian (Wolfcampian) is consistent with patterns of brachiopod extinction found previously by Olzewski and Patzkowsky (2001a), who attributed the increase to gradual sea-level fall and a shift toward more arid conditions evidenced by the increased appearance of evaporites. The weak to absent correlation between extinction and speciation rates is somewhat unusual (Stanley et al. 1990) and adds further credence to the notion this time period is exceptional with respect to macroevolutionary dynamics.

Table 3.

One-tailed binomial test results for the probability that species decline in range size and abundance before extinction using standardized variables. Species present in three or more consecutive stages were determined to decline before extinction if they conformed to the following conditions: (1) the terminal value was lower than the value in the immediately preceding stage, and (2) the terminal value was lower than the mean for the species. Values include abundance and geographic range size measured as either convex hull or latitudinal range.



Geographic range size and abundance estimates for species that survive and go extinct at stage boundaries (n = 88 extinct, n = 289 survive). A total of 164 unique species were analyzed across nine stages (Chesterian–Leonardian). Analyses used the median convex-hull area (km2) of all possible convex hulls created by jackknifing occurrences and was log transformed. Latitudinal range and abundance were square-root transformed from original count and degree latitude measurements, respectively. Species trait data were pooled, such that multiple measurements for species were recorded if they occurred in more than one stage.


Species that went extinct during the LPIA had smaller geographic range sizes and lower abundances (Fig. 3) than survivors. Although geographic range size has received robust support as a determinant of extinction risk during many, if not most, geologic intervals (Jablonski 1986; McKinney 1997; Payne and Finnegan 2007; Harnik et al. 2012; Saupe et al. 2015), we found that abundance, rather than geographic range size, was a strong predictor of extinction risk for midcontinental brachiopod species during the late Paleozoic of the North American midcontinent. These patterns receive further support from our analyses, which suggest that species' geographic range sizes did not decline before extinction but their abundance did (Table 3, Fig. 2). Our results reinforce previous work on LPIA invertebrates from the midcontinent. For example, Kolis and Lieberman (2019) found that geographic range sizes for cephalopod species did not correlate with extinction rates. Using a multiple linear regression, Powell (2007) demonstrated that abundance and geographic range size contributed equally to genus duration in midcontinental brachiopods. Although Powell's (2007) results for the importance of abundance are similar to those reported here, our results may differ for geographic range because of differences in statistical methodology, taxonomic level, and/or the use of summary metrics across the entire duration of the genus to characterize correlates of duration, rather than extinction risk.

The abundance–survivorship relationship reported herein could be confounded by ecological or physiological differences among species not controlled for in this analysis. Other potential correlates of extinction risk could include differences in basal metabolic rate (Strotz et al. 2018), trophic ecology (Norris 1992), factors affecting population structure (Norris 1992), variations in dispersal ability (Powell 2007; Birand et al. 2012), or local abundance (i.e., average relative abundance across samples or localities; Stanley 1986). Many of these factors are difficult to ascertain for brachiopods. For example, brachiopods utilize a range of larval development strategies (Thayer 1981; James et al. 1992; Peck and Robinson 1994) that cannot be inferred directly for extinct species in most cases (Valentine and Jablonski 1983; Rowell 1986), which makes dispersal ability difficult to estimate. It is also possible that implicit or explicit biases in the way fossils are collected and/or the way museum collections and databases are built fundamentally impugn the use of species abundance in these types of analyses.

Stanley and Powell (2003) and Powell (2005) suggested the low rates of speciation and extinction in the late Paleozoic represented a “new macroevolutionary state” relative to the rest of the Paleozoic. This dampened species turnover has been attributed to the loss of extinction-prone taxa earlier in the Paleozoic, which resulted in a higher proportion of widely distributed and long-duration brachiopod genera in the paleo-tropics, not only at higher latitudes (Brett and Baird 1995; Powell 2005). The physiographic conditions of the late Paleozoic sea in the North American midcontinent—a gently-sloping ramp with few geographic barriers to latitudinal movement—likely aided dispersal by marine organisms and decreased opportunities for both extinction and speciation (Heckel 1986; Olszewski and Patzkowsky 2003). Thus, enhanced abiotic potential for broad ranges in the late Paleozoic sea of the North American midcontinent, combined with the biotic characteristics of most brachiopods present, may have yielded a fauna that lacked a prominent geographic component of extinction susceptibility. Instead, the abundance–extinction relationship observed during the LPIA could reflect the vulnerability of populations with low abundance to the proximate causes of extinction, namely demographic stochasticity, environmental stochasticity, and genetic deterioration (Goodman 1987; Lande 1993; Pimm et al. 1993; Wissel et al. 1994; Henle et al. 2004a,b). Abundance may play an important role in regulating extinction risk during biodiversity crises or extinctions that, like the late Paleozoic, are characterized by low rates of origination (e.g., the Late Devonian biodiversity crisis or the end-Triassic mass extinction; Bambach et al. 2004). Such claims, however, require additional research.

Extinction and speciation rates are frequently correlated with geographic range size (Vrba 1980; Jablonski 1986; Stanley 1986, 1990b; Eldredge 1989; Jablonski and Roy 2003), and small-ranged taxa are often culled preferentially during extinction events of different scales (Jablonski 1986; Payne and Finnegan 2007; Powell 2007; Clapham and Payne 2011). Therefore, background and mass extinctions could lead to preferential loss of volatile taxa with intrinsically high rates of origination and extinction over time (Gilinksy 1994; Lieberman and Melott 2013). During “normal” evolutionary times, the extinction-prone, narrowly distributed taxa removed by extinction would be replaced after several million years, as new speciation events tend to generate small-ranged taxa (Foote 2007; Liow and Stenseth 2007; Antell et al. 2020). The eventual reappearance of these small-ranged taxa would have therefore reinstated geographic range size as an important determinant of extinction risk. However, the relative dearth of narrowly distributed taxa during the late Paleozoic, combined with the distinctive physiographic nature of the region studied, could have conspired to suppress both extinction and speciation rates and hinder emergence from this distinctive “macroevolutionary state”. This proposed mechanism shares many similarities with Stanley's (1990a) explanation for the apparent periodicity of mass extinctions, in which the suggested recurrence of extinction events is limited by the length of time necessary to recover extinction-prone or vulnerable taxa.

The abundance–extinction risk relationship uncovered here may be unique to our study system of the North American midcontinent during the LPIA, which focused on brachiopods, tropical latitudes, and epicontinental seaways. For example, study of brachiopod dynamics from midlatitudes and open ocean–facing systems found high regional turnover in response to the LPIA (Balseiro 2016) and regional extirpation selectivity associated with both genus range size and body size (Balseiro and Halpern 2019). These differing dynamics suggest that brachiopod faunas outside tropical, epicontinental seaways may have responded to the LPIA differently than those in the U.S. midcontinent. Indeed, macroevolutionary dynamics of epicontinental seas may differ from those operating along ocean-facing shelves (Miller and Foote 2009), such that extinction risk patterns could differ in open-shelf communities of the same age or in post-Paleozoic periods when epicontinental seas decrease in prevalence (Peters 2007). However, if the abundance–extinction risk relationship found herein is extendable to other regions and taxa during the late Paleozoic, this could add further nuance to Jablonksi's (1986) “alternation of macroevolutionary regimes,” which hypothesizes that mass extinctions are unique and not simply intensifications of background extinction dynamics. The possibility exists that macroevolutionary processes shift between intervals of background extinction and of mass extinction, but also between times of background extinction and intervals of sluggish turnover.


We thank M. Powell and M. Foote for thoughtful and constructive reviews. We thank museum staff for help with brachiopod collections, including S. Butts (YPM), J. Utrup (YPM), U. Farrell (KUMIP at the time, now N. Lopez Carranza). We thank T. West for assistance with ArcGIS and PaleoWeb. Financial support for this project was provided by NSF DEB 1256993, EF 1206757, and DBI 1602067 to B.S.L., and Leverhulme grant DGR01020 to E.E.S. This material is based upon work supported while working at the National Science Foundation.

Literature Cited


Algeo, T. J., and P. H. Heckel. 2008. The Late Pennsylvanian midcontinent sea of North America: a review. Palaeogeography, Palaeoclimatology, Palaeoecology 268:205–221. Google Scholar


Alroy, J. 2000. New methods for quantifying macroevolutionary patterns and processes. Paleobiology 26:707–733. Google Scholar


Alroy, J. 2015. A more precise speciation and extinction rate estimator. Paleobiology 41:633–639. Google Scholar


Antell, G. S., W. Kiessling, M. Aberhan, and E. E. Saupe. 2020. Marine biodiversity and geographic distributions are independent on large scales. Current Biology 30:115–121. Google Scholar


Balseiro, D. 2016. Compositional turnover and ecological changes related to the waxing and waning of glaciers during the late Paleozoic ice age in ice-proximal regions (Pennsylvanian, western Argentina). Paleobiology 42:335–357. Google Scholar


Balseiro, D., and K. Halpern. 2019. Immigration and extirpation selectivity patterns of brachiopods and bivalves across a Carboniferous glacial to non-glacial transition (Pennsylvanian, central western Argentina) and their influence in building the biotic bathymetric gradient. Palaeogeography, Palaeoclimatology, Palaeoecology 516:11–21. Google Scholar


Bambach, R. K., A. H. Knoll, and S. C. Wang. 2004. Origination, extinction, and mass depletions of marine diversity. Paleobiology 30:522–542. Google Scholar


Bates, D., M. Mächler, B. Bolker, and S. Walker. 2015. Fitting linear mixed-effects models using lme4. arXiv preprint arXiv:1406.5823. Google Scholar


Birand, A., A. Vose, and S. Gavrilets. 2012. Patterns of species ranges, speciation, and extinction. American Naturalist 179:1–21. Google Scholar


Bivand, R. S., E. J. Pebesma, and V. Gomez-Rubio. 2013. Applied spatial data analysis with R, 2nd ed. Springer, New York. Google Scholar


Brett, C. E., and Baird, G.C. 1995. Coordinated stasis and evolutionary ecology of Silurian to Middle Devonian faunas in the Appalachain Basin. Pp. 285–315 in D. H. Erwin and R. L. Anstey, eds. New approaches to speciaton in the fossil record. Columbia University Press, New York. Google Scholar


Brezinski, D. K. 1988. Appalachian Carboniferous trilobites. Journal of Paleontology 62:934–945. Google Scholar


Burnham, K. P., and D. R. Anderson. 2002. Model selection and multimodel inference: a practical information-theoretic approach, 2nd ed. Springer, New York. Google Scholar


Buzas, M. A., C. F. Koch, S. J. Culver, and N. F. Sohl. 1982. On the distribution of species occurrence. Paleobiology 8:143–150. Google Scholar


Carter, J. L., and R. C. Carter. 1970. Bibliography and index of North America Carboniferous brachiopods (1898–1968), Vol. 128. Geological Society of America, Boulder, Colo. Google Scholar


Clapham, M. E., and J. L. Payne. 2011. Acidification, anoxia, and extinction: a multiple logistic regression analysis of extinction selectivity during the Middle and Late Permian. Geology 39:1059–1062. Google Scholar


Darroch, S. A. F., and E. E. Saupe. 2018. Reconstructing geographic range-size dynamics from fossil data. Paleobiology 44:25–39. Google Scholar


Darroch, S. A. F., and P. J. Wagner. 2015. Response of beta diversity to pulses of Ordovician-Silurian mass extinction. Ecology 96:532–549. Google Scholar


Darroch, S. A. F., M. M. Casey, G. S. Antell, A. Sweeney, and E. E. Saupe. 2020. High preservation potential of paleogeographic range size distributions in deep time. American Naturalist. Scholar


Eldredge, N. 1989. Macroevolutionary dynamics. McGraw-Hill, New York. Google Scholar


Finnegan, S., S. C. Anderson, P. G. Harnik, C. Simpson, D. P. Tittensor, J. E. Byrnes, Z. V. Finkel, D. R. Lindberg, L. H. Liow, R. Lockwood, and H. K. Lotze. 2015. Paleontological baselines for evaluating extinction risk in the modern oceans. Science 348:567–570. Google Scholar


Finnegan, S., C. M. Ø. Rasmussen, and D. A. T. Harper. 2016. Biogeographic and bathymetric determinants of brachiopod extinction and survival during the Late Ordovician mass extinction. Proceedings of the Royal Society of London B 283:20160007. Google Scholar


Foote, M. 2000. Origination and extinction components of taxonomic diversity: Paleozoic and post-Paleozoic dynamics. Paleobiology 26:578–605. Google Scholar


Foote, M. 2007. Symmetric waxing and waning of marine invertebrate genera. Paleobiology 33:517–529. Google Scholar


Foote, M., and A. I. Miller. 2013. Determinants of early survival in marine animal genera. Paleobiology 39:171–192. Google Scholar


Foote, M., and J. J. Sepkoski. 1999. Absolute measures of the completeness of the fossil record. Nature 398:415–417. Google Scholar


Foote, M., J. S. Crampton, A. G. Beu, and R. A. Cooper. 2008. On the bidirectional relationship between geographic range and taxonomic duration. Paleobiology 34:421–433. Google Scholar


Gaston, K. J. 1994. Rarity. Chapman and Hall, London. Google Scholar


Gaston, K. J., Tim M. Blackburn, and John H. Lawton. 1997. Inter-specific abundance-range size relationships: an appraisal of mechanisms. Journal of Animal Ecology 66:579–601. Google Scholar


Gilinsky, N. L. 1994. Volatility and the Phanerozoic decline of background extinction intensity. Paleobiology 20: 445–458. Google Scholar


Goodman, D. 1987. The demography of chance extinction. Viable Populations for Conservation 11:34. Google Scholar


Harnik, P. G. 2011. Direct and indirect effects of biological factors on extinction risk in fossil bivalves. Proceedings of the National Academy of Sciences USA 108:13594–13599. Google Scholar


Harnik, P. G., C. Simpson, and J. L. Payne. 2012. Long-term differences in extinction risk among the seven forms of rarity. Proceedings of the Royal Society of London B 279:4969–4976. Google Scholar


He, F., and K. J. Gaston. 2000. Estimating species abundance from occurrence. American Naturalist 156:553–559. Google Scholar


He, F., and K. J. Gaston. 2003. Occupancy, spatial variance, and the abundance of species. American Naturalist 162:366–375. Google Scholar


Heckel, P. H. 1977. Origin of phosphatic black shale facies in Pennsylvanian cyclothems of mid-continent North America. AAPG Bulletin 61:1045–1068. Google Scholar


Heckel, P. H. 1986. Sea-level curve for Pennsylvanian eustatic marine transgressive-regressive depositional cycles along midcontinent outcrop belt, North America. Geology 14:330–334. Google Scholar


Heckel, P. H. 2008. Pennsylvanian cyclothems in Midcontinent North America as far-field effects of waxing and waning of Gondwana ice sheets. Pp. 275–290 In C. R. Fielding, T. D. Frank, and J. L. Isbell, eds. Resolving the Late Paleozoic ice age in time and space. Geological Society of America, Boulder, Colorado.. Google Scholar


Heckel, P. H., and G. Clayton. 2006. The Carboniferous System. Use of the new official names for the subsystems, series, and stages. Geologica Acta 4:403–407. Google Scholar


Heckel, P. H., J. M. Dennison, and F. R. Ettensohn. 1994. Evaluation of evidence for glacio-eustatic control over marine Pennsylvanian cyclothems in North America and consideration of possible tectonic effects. Tectonic and eustatic controls on sedimentary cycles. SEPM, Concepts in Sedimentology and Paleontology 4:65–87. Google Scholar


Hendricks, J. R., B. S. Lieberman, and A. L. Stigall. 2008. Using GIS to study palaeobiogeographic and macroevolutionary patterns in soft-bodied Cambrian arthropods. Palaeogeography, Palaeoclimatology, Palaeoecology 264:163–175. Google Scholar


Henle, K., K. F. Davies, M. Kleyer, C. Margules, and J. Settele. 2004a. Predictors of species sensitivity to fragmentation. Biodiversity and Conservation 13:207–251. Google Scholar


Henle, K., S. Sarre, and K. Wiegand. 2004b. The role of density regulation in extinction processes and population viability analysis. Biodiversity and Conservation 13:9–52. Google Scholar


Hoare, R. D. 1961. Desmoinesian Brachiopoda and Mollusca from southwestern Missouri. University of Missouri Studies 36. Google Scholar


Holt, R., J. Lawton, K. J. Gaston, and T. Blackburn. 1997. On the relationship between range size and local abundance: back to basics. Oikos 78:183–190. Google Scholar


Horton, D. E., C. J. Poulsen, I. P. Montañez, and W. A. DiMichele. 2012. Eccentricity-paced late Paleozoic climate change. Palaeogeography, Palaeoclimatology, Palaeoecology 331:150–161. Google Scholar


International Union for the Conservation of Nature. 2001. IUCN Red List Categories and Criteria, Version 3.1. IUCN, Gland, Switzerland. Google Scholar


Jablonski, D. 1986. Background and mass extinctions: the alternation of macroevolutionary regimes. Science 231.4734:129–133. Google Scholar


Jablonski, D., and K. Roy. 2003. Geographical range and speciation in fossil and living molluscs. Proceedings of the Royal Society of London B 270:401–406. Google Scholar


Jackson, J. B. C. 1974. Biogeographic consequences of eurytopy and stenotopy among marine bivalves and their evolutionary significance. American Naturalist 108:541–560. Google Scholar


James, M. A., A. D. Ansell, M. J. Collins, G. B. Curry, L. S. Peck, and M. C. Rhodes. 1992. Biology of living brachiopods. Advances in Marine Biology 28:175–387. Google Scholar


Kiessling, W., and M. Aberhan. 2007. Geographical distribution and extinction risk: lessons from Triassic–Jurassic marine benthic organisms. Journal of Biogeography 34:1473–1489. Google Scholar


Kiessling, W., N. B. Raja, V. J. Roden, S. T. Turvey, and E. E. Saupe. 2019. Addressing priority questions of conservation science with palaeontological data. Philosophical Transactions of the Royal Society of London B 374:20190222. Google Scholar


Kocsis, A. T., C. J. Reddin, J. Alroy, and W. Kiessling. 2019. The R package divDyn for quantifying diversity dynamics using fossil sampling data. Methods in Ecology and Evolution 10:735–743. Google Scholar


Kolis, K. M., and B. S. Lieberman. 2019. Using GIS to examine bio-geographic and macroevolutionary patterns in some late Paleozoic cephalopods from the North American Midcontinent Sea. PeerJ 7:e6910. Google Scholar


Kunin, W. E. 1998. Extrapolating species abundance across spatial scales. Science 281:1513–1515. Google Scholar


Lande, R. 1993. Risks of population extinction from demographic and environmental stochasticity and random catastrophes. American Naturalist 142:911–927. Google Scholar


Lee, T. M., and W. Jetz. 2011. Unravelling the structure of species extinction risk for predictive conservation science. Proceedings of the Royal Society of London B 278:1329–1338. Google Scholar


Lieberman, B. S., and A. L. Melott. 2013. Declining volatility, a general property of disparate systems: from fossils, to stocks, to the stars. Palaeontology 56:1297–1304. Google Scholar


Liow, L. H., and N. C. Stenseth. 2007. The rise and fall of species: implications for macroevolutionary and macroecological studies. Proceedings of the Royal Society of London B 274:2745–2752. Google Scholar


Lupia, R., and J. L. Armitage. 2013. Late Pennsylvanian–Early Permian vegetational transition in Oklahoma: palynological record. International Journal of Coal Geology 119:165–176. Google Scholar


McKinney, M. L. 1997. Extinction vulnerability and selectivity: combining ecological and paleontological views. Annual Review of Ecology and Systematics 28:495–516. Google Scholar


Menning, M., A. S. Alekseev, B. I. Chuvashov, V. I. Davydov, F. X. Devuyst, H. C. Forke, T. A. Grunt, L. Hance, P. H. Heckel, N. G. Izokh, Y. G. Jin, P. J. Jones, G. V. Kotlyar, H. W. Kozur, T. I. Nemyrovska, J. W. Schneider, X. D. Wang, K. Weddige, D. Weyer, and D. M. Work. 2006. Global time scale and regional stratigraphic reference scales of Central and West Europe, East Europe, Tethys, South China, and North America as used in the Devonian–Carboniferous–Permian Correlation Chart 2003 (DCP 2003). Palaeogeography, Palaeoclimatology and Palaeoecology 240:318–372. Google Scholar


Meseguer, A. S., J. M. Lobo, R. Ree, D. J. Beerling, and I. Sanmartín. 2015. Integrating fossils, phylogenies, and niche models into bio-geography to reveal ancient evolutionary history: the case of Hypericum (Hypericaceae). Systematic Biology 64:215–232. Google Scholar


Miller, A. I., and M. Foote. 2009. Epicontinental seas versus open-ocean settings: the kinetics of mass extinction and origination. Science 326:1106–1109. Google Scholar


Montañez, I. P., and C. J. Poulsen. 2013. The Late Paleozoic ice age: an evolving paradigm. Annual Review of Earth and Planetary Sciences 41:629–656. Google Scholar


Moore, R. C. 1964. Paleoecological aspects of Kansas Pennsylvanian and Permian cyclothems. Pp. 287–380 in Symposium on cyclic sedimentation, Vol. 169. Kansas Geological Survey Bulletin, Lawrence. Google Scholar


Muir-Wood, H. M., and G. A. Cooper. 1960. Morphology, classification and life habits of the Productoidea (Brachiopoda), Vol. 81. Geological Society of America, New York. Google Scholar


Myers, C. E., and B. S. Lieberman. 2011. Sharks that pass in the night: using geographical information systems to investigate competition in the Cretaceous Western Interior Seaway. Proceedings of the Royal Society of London B 278:681–689. Google Scholar


Norris, R. D. 1992. Extinction selectivity and ecology in planktonic foraminifera. Palaeogeography, Palaeoclimatology, Palaeoecology 95:1–17. Google Scholar


Olszewski, T. D., and M. E. Patzkowsky. 2001a. Evaluating taxonomic turnover: Pennsylvanian–Permian brachiopods and bivalves of the North American Midcontinent. Paleobiology 27:646–668. Google Scholar


Olszewski, T. D., and M. E. Patzkowsky. 2001b. Measuring recurrence of marine biotic gradients: a case study from the Pennsylvanian–Permian Midcontinent. Palaios 16:444–460. Google Scholar


Olszewski, T. D., and M. E. Patzkowsky. 2003. From cyclothems to sequences: the record of eustasy and climate on an icehouse epeiric platform (Pennsylvanian–Permian, North American Midcontinent). Journal of Sedimentary Research 73:15–30. Google Scholar


Page, L. M., B. J. MacFadden, J. A. Fortes, P. S. Soltis, and G. Riccardi. 2015. Digitization of biodiversity collections reveals biggest data on biodiversity. BioScience 65:841–842. Google Scholar


Parrish, J. T. 1993. Climate of the supercontinent Pangea. Journal of Geology 101:215–233. Google Scholar


Payne, J. L., and S. Finnegan. 2007. The effect of geographic range on extinction risk during background and mass extinction. Proceedings of the National Academy of Sciences USA 104:10506–10511. Google Scholar


Payne, J. L., S. Truebe, A. Nützel, and E. T. Chang. 2011. Local and global abundance associated with extinction risk in late Paleozoic and early Mesozoic gastropods. Paleobiology 37:616–632. Google Scholar


Pebesma, E., and R. S. Bivand. 2005. S classes and methods for spatial data: the sp package. R News 5:9–13. Google Scholar


Peck, L. S., and K. Robinson. 1994. Pelagic larval development in the brooding Antarctic brachiopod Liothyrella uva. Marine Biology 120:279–286. Google Scholar


Peters, S. E. 2007. The problem with the Paleozoic. Paleobiology 33:165–181. Google Scholar


Pimm, S. L., Jared Diamond, T. M. Reed, G. J. Russell, and J. Verner. 1993. Times to extinction for small populations of large birds. Proceedings of the National Academy of Sciences USA 90:10871–10875. Google Scholar


Powell, M. G. 2005. Climatic basis for sluggish macroevolution during the late Paleozoic ice age. Geology 33:381–384. Google Scholar


Powell, M. G. 2007. Geographic range and genus longevity of late Paleozoic brachiopods. Paleobiology 33:530–546. Google Scholar


Raymond, A., and C. Metz. 2004. Ice and its consequences: glaciation in the Late Ordovician, Late Devonian, Pennsylvanian-Permian, and Cenozoic compared. Journal of Geology 112.6:655–670. Google Scholar


R Core Team. 2017. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.  https://www.R-project.orgGoogle Scholar


Rowell, A. J. 1986. The distribution and inferred larval dispersion of Rhondellina dorei: a new Cambrian brachiopod (Acrotretida). Journal of Paleontology 60:1056–1065. Google Scholar


Saupe, E. E., H. Qiao, J. R. Hendricks, R. W. Portell, S. J. Hunter, J. Soberón, and B. S. Lieberman. 2015. Niche breadth and geographic range size as determinants of species survival on geological time scales. Global Ecology and Biogeography 24:1159–1169. Google Scholar


Schnute, J. T., N. Boers, R. Haigh, C. Grandin, A. Johnson, P. Wessel, and F. Antonio. 2013. PBSmapping: mapping fisheries data and spatial analysis tools, R package version 2. Accessed 15 December 2017. Google Scholar


Scotese, C. R. 2016. PALEOMAP PaleoAtlas for GPlates and the PaleoData Plotter Program, PALEOMAP Project. Accessed 16 June 2020. Google Scholar


Segessenman, D. C., and T. W. Kammer. 2018. Testing reduced evolutionary rates during the Late Palaeozoic Ice Age using the crinoid fossil record. Lethaia 51:330–343. Google Scholar


Smits, P., and S. Finnegan. 2019. How predictable is extinction? Forecasting species survival at million-year timescales. Philosophical Transactions of the Royal Society of London B 374:20190392. Google Scholar


Stanley, S. M. 1986. Population size, extinction, and speciation: the fission effect in Neogene Bivalvia. Paleobiology 12:89–110. Google Scholar


Stanley, S. M. 1990a. Delayed recovery and the spacing of major extinctions. Paleobiology 16:401–414. Google Scholar


Stanley, S. M. 1990b. The general correlation between rate of speciation and rate of extinction: fortuitous causal linkages. Pp. 103–127 in R. M. Ross and W. D. Allmon, eds. Causes of evolution: a paleontological perspective. University of Chicago Press, Chicago. Google Scholar


Stanley, S M., and M. G. Powell. 2003. Depressed rates of origination and extinction during the late Paleozoic ice age: a new state for the global marine ecosystem. Geology 31:877–880. Google Scholar


Stigall, A. L. 2010. Invasive species and biodiversity crises: testing the link in the Late Devonian. PLoS ONE 5:e15584. Google Scholar


Stigall, A. L., and B. S. Lieberman. 2006. Quantitative palaeobiogeography: GIS, phylogenetic biogeographical analysis, and conservation insights. Journal of Biogeography 33:2051–2060. Google Scholar


Strotz, L. C., E. E. Saupe, J. Kimmig, and B. S. Lieberman. 2018. Metabolic rates, climate and macroevolution: a case study using Neogene molluscs. Proceedings of the Royal Society of London B 285:20181292. Google Scholar


Sunday, J. M., A. E. Bates, and N. K. Dulvy. 2012. Thermal tolerance and the global redistribution of animals. Nature Climate Change 2:686–690. Google Scholar


Thayer, C. W. 1981. Ecology of living brachiopods. Studies in Geology, Notes for a Short Course 5:110–126. Google Scholar


U.S. Geological Survey. n.d. National Geologic Map Database: Geolex Search. Accessed 13 October 2017. Google Scholar


Valentine, J. W., and D. Jablonski. 1983. Larval adaptations and patterns of brachiopod diversity in space and time. Evolution 37:1052–1061. Google Scholar


Venables, W. N., and B. D. Ripley. 2002. Random and mixed effects. Pp. 271–300 in W. N. Venables and B. D. Ripley, eds. Modern applied statistics with S. Springer, New York. Google Scholar


Vrba, E. S. 1980. Evolution, species and fossils: how does life evolve. South African Journal of Science 76:61–84. Google Scholar


Williams, A., A. J. Rowell, H. M. Muir-Wood, C. W. Pitrat, H. Schmidt, F. G. Stehli, D. V. Ager, A. D. Wright, G. F. Elliott, T. W. Amsden, M. J. S. Rudwick, K. Hatai, G. Biernat, D. J. McLaren, A. J. Boucot, J. G. Johnson, R. D. Stanton, R. E. Grant, and H. M. Jope. 1965. Brachiopoda, Part H. R. C. Moore, ed. Treatise on invertebrate paleontology. Geological Society of America, New York. Google Scholar


Wissel, C., T. Stephan, and S. H. Zaschke. 1994. Modelling extinction and survival of small populations. Pp. 67–103 in Minimum animal populations. H. Remmert, ed. Springer, Berlin. Google Scholar


Zuur, A. F., E. N. Ieno, and C. S. Elphick. 2010. A protocol for data exploration to avoid common statistical problems. Methods in Ecology and Evolution 1:3–14. Google Scholar
© The Author(s), 2020. Published by Cambridge University Press on behalf of The Paleontological Society. This is an Open Access article, distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike licence (, which permits non-commercial re-use, distribution, and reproduction in any medium, provided the same Creative Commons licence is included and the original work is properly cited. The written permission of Cambridge University Press must be obtained for commercial re-use.
Michelle M. Casey, Erin E. Saupe, and Bruce S. Lieberman "The effects of geographic range size and abundance on extinction during a time of “sluggish”' evolution," Paleobiology 47(1), 54-67, (5 March 2021).
Accepted: 9 October 2020; Published: 5 March 2021

Get copyright permission
Back to Top