Below is the uncorrected machine-read text of this chapter, intended to provide our own search engines and external engines with highly rich, chapter-representative searchable text of each book. Because it is UNCORRECTED material, please consider the following text as a useful but insufficient proxy for the authoritative book pages.
10 Extinction and the Spatial Dynamics of Biodiversity DAVID JABLONSKI The fossil record amply shows that the spatial fabric of extinc- tion has profoundly shaped the biosphere; this spatial dimension provides a powerful context for integration of paleontological and neontological approaches. Mass extinctions evidently alter extinc- tion selectivity, with many factors losing effectiveness except for a positive relation between survivorship and geographic range at the clade level (confirmed in reanalyses of end-Cretaceous extinction data). This relation probably also holds during âânor- malââ times, but changes both slope and intercept with increasing extinction. The strong geographical component to clade dynamics can obscure causation in the extinction of a feature or a clade, owing to hitchhiking effects on geographic range, so that multifac- torial analyses are needed. Some extinctions are spatially com- plex, and regional extinctions might either reset a diversity ceiling or create a diversification debt open to further diversification or invasion. Evolutionary recoveries also exhibit spatial dynamics, including regional differences in invasibilty, and expansion of clades from the tropics fuels at least some recoveries, as well as biodiversity dynamics during normal times. Incumbency effects apparently correlate more closely with extinction intensities than with standing diversities, so that regions with higher local and global extinctions are more subject to invasion; the latest Ceno- zoic temperate zones evidently received more invaders than the tropics or poles, but this dynamic could shift dramatically if tropi- Department of Geophysical Sciences, University of Chicago, 5734 South Ellis Avenue, Chi- cago, IL 60637. 189
190â /â David Jablonski cal diversity is strongly depleted. The fossil record can provide valuable insights, and their application to present-day issues will be enhanced by partitioning past and present-day extinctions by driving mechanism rather than emphasizing intensity. T he inventory of life on Earth has always been determined at the most basic level by the difference between origination and extinc- tion. This fundamental macroevolutionary equation, richness = origination â extinction, has been formally applied in many ways and with many elaborations, but takes on special consequence when attempting to evaluate the processes shaping present-day biodiversity, where neither term in the right side of the equation can be observed directly. Some prog- ress has been made in modeling these parameters, but most approaches involve strong assumptions, require very large datasets, and carry large uncertainties [e.g., Paradis (2004), Ricklefs (2007a,b)]. The spatially explicit form of this equation, where richness is a local or regional pool of species or higher taxa, and immigration and emigration terms enter the right side of the equation (Jablonski et al., 2006), is important in many situations, from biotic responses to Pleistocene climate cycles and ongoing climate changes to recoveries from mass extinctions. However, this form is even more difficult to apply rigorously without historical data, and my empha- sis here will be on the fossil record. Few would argue against the idea that the spatial fabric of extinction has shaped, and will continue to shape, the biosphere in profound ways, but spatial effects have been neglected rela- tive to temporal patterns (partly because documentation is so challenging). I will argue that the insights beginning to emerge from spatially explicit approaches to ancient extinctions have significant implications for the dynamics of diversity of the past and in the future. This discussion will encompass a range of extinction intensities and focus on marine systems, where the fossil record is richest: application of these generalizations to terrestrial realms requires more study. Opinions are divided on whether the handful of mass extinctions of the geologic past are a separate class of intensities from the ââbackgroundââ extinction that constitutes the great bulk of geologic time [and the bulk of total extinction; Raup (1994)], but this is a secondary issue that can probably be resolved by factoring out the well-known secular decline in background extinction rates (Bambach et al., 2004; Jablonski, 2005; Stanley, 2007). As discussed below, extinction selectivities evidently shift between episodes of low and high extinction rates, and this selectivity is the more important issue for understanding the role of extinction in shaping past and future biotas. I will corroborate previous evidence for a strong spatial component to survivorship during major extinction events, present a multifactorial analysis of the end-Cretaceous (K-T) mass extinction in which geographic
Extinction and the Spatial Dynamics of Biodiversityâ /â 191 range emerges as the best predictor of survivorship in marine bivalves, and argue that such indirect effects are probably more important than gen- erally appreciated. I will also discuss regional variations in the balance of invasions and local origination in the aftermath of the K-T event, which are somewhat unexpected given that the extinction itself tended to increase biotic homogenization on a global scale by preferentially removing the more localized taxa. Invasions and extinctions are also important during times of âânormalââ extinction intensities, as I will illustrate with reference to the dynamics of the latitudinal diversity gradient. I will conclude with some implications for integrating insights for past and present-day extinc- tions and suggest that a powerful approach might involve comparative dissection in extinction patterns according to likely drivers. Throughout I will note gaps in our understanding that would benefit from combined study of modern and ancient systems. In this chapter, I will focus mainly on marine bivalves such as mus- sels, scallops, and cockles. Bivalves are becoming a model system for the analysis of large-scale biogeographic and evolutionary patterns (Crame, 2000, 2002; Jablonski et al., 2003b, 2006; Kidwell, 2005; Valentine et al., 2006; Krug et al., 2007) for several reasons. They are taxonomically rich but not unmanageable (â3,000 living and fossil genera), and their systematics are increasingly understood, so that taxonomic standardization and phylo- genetic treatment of heterogeneous data are feasible. They have diverse life habits, from filter-feeding to photosymbiosis and chemosymbiosis to carnivory. They occur at all depths from the intertidal zone to deep-sea trenches and from the tropics to the poles. They are abundant and often well preserved as fossils [although not all habitats and clades are equally represented; Valentine et al. (2006)], and they have diverse shell mineralo- gies and microstructures, which allows analyses to control statistically for, and thus factor out, some, although not all, of the biases in the fossil record (Kidwell, 2005; Valentine et al., 2006). These favorable attributes do not mean that the bivalve fossil record is perfect, and preservation and sampling biases must always be considered in large-scale analyses [see, for example, the variety of approaches in Alroy (2000), Foote (2003), Bush and Bambach (2004), Bush et al. (2004), Jablonski et al. (2006), and Smith (2007)]. However, our growing knowledge of living and fossil bivalves, including the taxonomic, preservational, and geographic factors that can distort their fossil record, makes this group an excellent vehicle for inte- grating present-day and paleontological diversity dynamics. EXTINCTION SELECTIVITY CHANGES AT THE MOST EXTREME EVENTS A broad array of organismic and clade-level traits enter into extinc- tion risk for present-day species. For example, in evaluating extinction
192â /â David Jablonski risk in present-day terrestrial vertebrates, Purvis and colleagues (Purvis et al., 2000b, 2005a; Davies et al., Chapter 14, this volume) found mixed, but significant, effects for body size, a consistent inverse relation between extinction risk and both abundance and geographic range, and either a positive relation or no effect for habitat specialization. Similar patterns are seen in the fossil record. For example, the geographic range is a significant determinant of Cretaceous and Cenozoic molluscan species duration or survivorship [Hunt et al. (2005) and Jablonski and Hunt (2006) and refer- ences therein], and Paleozoic crinoids show a significant positive relation between habitat breadth and species duration (Kammer et al., 1998). Pre- dictable interactions among factors can also be seen, although this aspect needs much more work. Molluscan genera containing many widespread species tend to be more extinction-resistant, with a median duration of 130 million years (Myr), than genera having just a few, localized species, which show a median duration of 32 Myr, and the genera with the other combinations give intermediate values (Jablonski, 2005). These are not theoretically surprising results, but it is encouraging that the paleontologi- cal outcomes so clearly match expectations. Extinction selectivity appears to change significantly at the most severe mass extinctions, however. The rules of survivorship changed during the K-T extinction, such that species-richness and species-level range failed to predict genus survivorship, singly or in concert [Jablonski (1986a, 2005) and see Kiessling and Baron-Szabo (2004) for comparable results for K-T corals]. In fact, survivorship of marine invertebrates in the K-T mass extinction is unrelated to a number of factors that have been shown or hypothesized to be important during more normal times. Besides the two already mentioned, these factors include local abundance, mode of larval development (which is in turn related to fecundity and species-level dispersal capability), estimated generation time, living position relative to the sedimentâwater interface, and trophic strategy (Jablonski, 2005). Despite this loss in effectiveness of a variety of organismic, species- and even clade-level traits, survivorship at mass extinction boundaries is not random. Every event seems to show some degree of selectivity, but one factor that seems to have promoted survival for most major groups and most mass extinctions is broad geographic distribution at the clade level (i.e., genera), regardless of species-level geographic ranges. This effect, which has been recorded for many groups and all of the major mass extinc- tions [see Jablonski (2005) for a tabulation], is again further corroborated in an extensively revised version of Jablonski and Raupâs (1995) data on K-T bivalves (Fig. 10.1A and B). This is more than a simple binary effect: bivalve genus extinction is inversely related to geographic range, with strong concordance between the new and old data (Fig. 10.1C). The 70% extinction suffered by the genera found in just one or two biogeographic
Extinction and the Spatial Dynamics of Biodiversityâ /â 193 A 70 60 Victims B 30 Survivors 25 50 n=172 n=117 Frequency Frequency 20 40 15 30 20 10 10 5 0 0 1 3 5 7 9 11 13 15 1 3 5 7 9 11 13 15 Number of Provinces Number of Provinces C 0.8 D 35 Rudists 0.7 2007 Extinxtion Intensity 30 0.6 25 (=pachyodont hinge) 0.5 Frequency 20 0.4 15 0.3 1995 10 0.2 0.1 5 0 0 1 2 3 4 5 6 7 8 1 3 5 7 9 11 13 15 Number of Provinces Number of Provinces FIGURE 10.1â Spatial effects in the end-Cretaceous (K-T) mass extinction for ma- rine bivalve genera. (A, B, and D) Victims of the K-T extinction (A) tended to be significantly less widespread than surviving bivalve genera (B), as measured by the number of biogeographic provinces they occurred in during the Maastrichtian zpq9990837740001.g.eps stage just before the event [MannâWhitney U test, P = 0.00001; new analysis based on an extensive, in-progress revision and update of Jablonski and Raup (1995), omitting rudist bivalves (D) as before, note that their inclusion as narrow-ranging victims would strengthen this result]. Adding provinces to fill gaps in observed geographic ranges strengthens the separation between victims and survivors (75% of victims are unchanged in range size and their median range is unchanged at two provinces; 60% of survivors are unchanged and their median increases from four to five provinces). Some caution is needed, because the proportion of survi- vors is likely to increase with phylogenetic analysis and further taxonomic stan- dardization of early Cenozoic bivalves, but the major pattern is unlikely to change. (C) Significant inverse relation between extinction intensity and the number of biogeographic provinces occupied by bivalve genera during the K-T extinction (Spearman rank test, P < 0.01). Solid line indicates analysis of revised dataset (n = 289 genera). Broken line indicates analysis of previous version of dataset [Jablonski and Raup (1995); n = 297 genera; 28 genera were added and 36 genera were removed in the revision]. (D) Loss of a major adaptation (the pachyodont hinge) by hitchhiking on geographic distribution. The unique pachyodont hinge structure disappeared with the extinction of these genera at the K-T boundary, signaling the termination of the rudist bivalves (Order Hippuritoida).
194â /â David Jablonski provinces is significantly higher than the 20% losses seen for genera found in eight or more provinces [of a global total of 16; following Jablonski and Raup, 1995)]. That said, even 20% represents a major, and highly unusual, drawdown of diversity in this most extinction-resistant part of the biota, equivalent to losing 20% of the most widespread genera in the sea today, such as the mussels (Mytilus, Modiolus) and the scallops (Pecten, Chlamys). [Although not ideal in some respects, analyses were conducted at the provincial scale rather than based on occurrences at individual localities, because clades are distributed not along simple linear coastlines, thereby undermining the use of linear distances or simple latitude/longitude extremes. Binning by province also damps some aspects of sampling and taxonomic uncertainty at the genus level, the range-endpoints of present- day molluscan genera tend to cluster at province boundaries (Campbell and Valentine, 1977; Roy et al., 1996), and the results are robust to different approaches to quantifying province-based range sizes.] Multifactorial analyses corroborate the importance of clade-level dis- tribution in determining survival during mass extinctions and show the value of testing for interaction among factors. For example, if variables are treated independently in the updated K-T dataset, geographic range remains the most important factor in clade survivorship, but species rich- ness also appears to play a significant role (and body size is insignificant as a survivorship predictor). However, multiple logistic regression models taking the three variables simultaneously into account, using Akaikeâs Information Criterion (AIC) as a basis for model selection (Burnham and Anderson, 2002), shows species richness to covary with range such that when range is factored out, richness has an insignificant effect on survi- vorship (P = 0.85, as opposed to P = 0.002 for clade range, in the multiple- factor model). Body size also enters into the multiple-factor model as a weak, but significant, variable, but the multiple-factor model does not have significantly more explanatory power than the geographic range model alone, according both to the similar AIC weights (Table 10.1) and a likelihood ratio test [P = 0.09; see Hosmer and Lemeshow (2000)]. Multi- variate approaches will help clarify patterns of extinction selectivity, even if, as here, they show that survivorship virtually collapses to the single variable of geographic range for K-T bivalves. The overlapping variation in range size among the victims and survivors suggests, however, that additional factors, or strong stochastic elements, enter into the fates of individual clades. Widespread clades are probably always extinction-resistant compared with narrow-ranging relatives. However, during times of low extinction intensity, range is evidently just one significant feature among many, becoming increasingly important as the crowd of factors influencing taxon duration falls away as intensity mounts. How the relation between range
Extinction and the Spatial Dynamics of Biodiversityâ /â 195 TABLE 10.1â Testing Models for Bivalve Genus Survivorship During the K-T Mass Extinction No. of Models Parameters AIC Weight P G** + R + B* 4 347.9 0.58 0.002/0.85/0.03 Geographic range*** 2 348.6 0.40 e-6 Species richness*** 2 356.2 0.02 0.0001 Body size 2 373.5 e-6 0.94 NOTES: When geographic range (G), species richness (R), and body size (B) are analyzed as independent factors, G is the most important factor, but R is also significant. When the three are analyzed together, R is not a significant factor. Note that the combined model is not significantly better than geographic range alone according to the AIC (for model selection, which essentially weighs the adding of parameters against the improved explanatory power of each model). *P < 0.05; **P < 0.01; ***P < 0.001. and extinction risk varies with extinction intensity is not known and is difficult to assess. The few available data suggest that the relation between extinction probability has a steep slope during times of low extinction intensity, with the most widespread genera suffering negligible extinction at those times (Payne and Finnegan, 2007; Powell, 2007a) (Fig. 10.2A). The simplest view would be that perturbations generally operate at too small a spatial scale to affect these most widespread elements of the global biota. In the major mass extinctions, the y-intercept increases, so that a greater fraction of taxa are lost from all range classes, and the slope probably decreases (Fig. 10.2B). This configuration is much more demanding of the data, so that sparse or noisy data may fail to capture that shallower A âBackgroundâ Extinction B Mass Extinction Extinction Intensity Extinction Intensity Geographic Range Geographic Range FIGURE 10.2â The inverse relation between geographic range and extinction risk appears to vary with severity of extinction. Conceptual model for this variation, zpq9990837740002.g.eps such that both slope and intercept may change between times of background (A) and mass extinction (B).
196â /â David Jablonski slope. We know little about whether the slope and intercept change con- tinuously or shift abruptly at thresholds. For obvious reasons, including some very practical ones relating to the present-day biota, this would be good to know. We also know relatively little about the determinants of geographic range size at the clade level. Organismic traits such as dispersal ability and ecological strategy must play a role, but interactions with biogeographic context, clade history, and many other factors, including the way that clades extend their ranges by speciation across barriers, serve to decouple geographic ranges at the species and clade levels. For example, the geo- graphic ranges of the 213 marine bivalve genera present today at shelf depths on the eastern Pacific margin from Point Barrow, Alaska to Cape Horn, Chile are not significantly related to the median or maximum ranges of their respective constituent species (Jablonski, 2005). Genera can attain broad ranges via a few widespread species, a mosaic of nonoverlapping but narrow-ranging species, or any combination thereof, each apparently equivalent in a mass extinction event (although this equivalence deserves further study). Genus ranges are not simply species attributes writ large, but involve a dynamic that is set at another hierarchical level, by the com- plexities of speciation, species extinction, and range expansion. This strong spatial component to extinction selectivity suggests that survival can be determined by features that are not tightly linked to the organismic and species-level traits that are favored, indeed shaped, during times of lower extinction intensities. Even well-established clades or adap- tations could be lost simply because they are not associated with those few features that enhance survivorship during unusual, and geologically brief, high-intensity events. As discussed below, the removal of incumbents and the subsequent diversification of formerly marginal taxa are essential ele- ments of the evolutionary dynamic fueled by major extinctions [see also Erwin (2001), Jablonski (2001, 2005), and Erwin, Chapter 9, this volume]. These results also suggest that hitchhiking effects may be mistaken for direct selectivity more often than generally appreciated. Biological traits tend to covary, even across hierarchical levels, and so selection on one feature can drag others along with it, hampering efforts to pinpoint cause and byproduct. Such hitchhiking was detected for bivalve species rich- ness in Table 10.1. Whenever widespread or restricted taxa tend to occupy nonrandom regions of phenotype space, for example in body sizes, trophic habits, or metabolic rates within a major group, hitchhiking becomes a real possibility, an interesting interaction across hierarchical levels where the extinction probabilities of organism-level characters are conditioned by clade-level properties (Jablonski, 2007). For example, the rudist bivalves of the Cretaceous seas (Order Hippuritoida) represent an extreme repat-
Extinction and the Spatial Dynamics of Biodiversityâ /â 197 terning of the bivalve body plan, with highly modified conical shells and a unique, pachyodont hinge structure (Skelton, 1985; Seilacher, 1998). This clade and its bizarre growth form, including the pachyodont hinge, disappeared in the K-T mass extinction (Steuber et al., 2002), but this loss was probably related less to any disadvantage inherent in the remarkable hinge apparatus than to the restricted ranges of rudist clades (Fig. 10.1D), and perhaps to their reliance in at least some instances on photosymbionts [Seilacher (1998), but see Steuber (2000)], as seen in modern reef-build- ing corals. If the range-frequency distribution of rudists played a role in their demise, with correlated morphologies carried along, then we would expect a similar pattern for the other bivalve orders. This is in fact the case: the five bivalve orders with median genus ranges of one or two provinces suffered significantly more severe K-T bottlenecks (median = 93% genus extinction) than the four orders with median genus ranges of three or more provinces (median = 32% genus extinction; Spearmanâs rank correlation of median genus range and extinction intensity for orders = 0.74, P = 0.02), as predicted by the hitchhiking argument for rudists. More detailed analyses must await a morphometric or discrete-character study combined with a well-resolved phylogeny of bivalve genera, and these results suggest that such studies would be worthwhile. The hitchhiking of such striking adaptations on the less flamboyant features that actually determine extinction resistance is probably perva- sive, both during background times [hence the large literature on com- parative methods and phylogenetic autocorrelation; e.g., Freckleton et al. (2002) and Paradis (2005)] and during mass extinctions. For example, marine bryozoan genera with complex colonies suffer more severely dur- ing mass extinctions than simple genera, but colony complexity is also inversely related to genus-level geographic range (Anstey, 1978, 1986), which may well be the ultimate basis for differential survival during the end-Ordovician mass extinction. The end-Ordovician extinction also preferentially removed snails with broad selenizones providing access to the mantle cavity, and planktonic graptolites with multiple stipes creating complex pendant colonies; the K-T extinction also preferentially removed bivalves with schizodont hinges (trigonioids), echinoids with elongate rostra (a clade of holasteroids), cephalopods with complex sutures (ammo- nites), and a major clade of birds with foot bones that fused from the ankles to the toes (Enantiornithes). All of these losses or severe bottlenecks are more likely to represent correlations, not necessarily with geographic range, but with some other organismic or higher-level factor, rather than direct selectivity on the most striking morphology or functional trait. These extinctions nonetheless truncated or rechanneled evolutionary tra- jectories through morphospace, and additional examples are plentiful.
198â /â David Jablonski SOME EXTINCTIONS ARE SPATIALLY COMPLEX The K-T extinction is remarkably homogenous on a global scale, except perhaps for greater intensity in tropical carbonate settings (Raup and Jablonski, 1993; Jablonski, 2005). However, other extinction events, particularly those that are less severe on a global scale, tend to show more spatial structure. For example, the mid-Cretaceous (end-Cenomanian) marine extinction appears to have been concentrated in northern Europe and the Western Interior Seaway of North America. The smaller events in the geologic record must be interpreted critically, because at least some of them may represent, or at least be heavily overprinted by, sampling variations [e.g., Foote (2003) and Smith and McGowan (2007)]. However, a seemingly genuine extinction pulse or regional series of pulses occurs in the oceans near the start of the Pleistocene. These regional extinctions are generally taken to represent a culling of taxa unable to cope with the onset of rapid, high-amplitude climate swings and oceanographic shifts that typify the Pleistocene. They vary in intensity and occur at slightly different times among regions [e.g., Kitamura et al. (2000), Monegatti and Raffi (2001), Todd et al. (2002), Smith and Roy (2006), Rivadeneira and Marquet (2007)], perhaps owing to regional variations in the tim- ing of oceanographic transitions toward a glacial state [e.g., Ravelo et al. (2004)]; the many subsequent glacial-interglacial cycles evidently drove few extinctions in marine or terrestrial settings [e.g., Huntley (2005)]. This spatial structuring of Plio-Pleistocene extinctions is interesting from many perspectives, but perhaps the most urgent need is to under- stand the dynamical consequences of these extinctions, which bear directly on the path of future biodiversity. Do these events reflect the setting of a new regional diversity level, such that taxa capable of weathering the volatile Pleistocene climate regime are more generalized and thus structure a biota capable of accommodating fewer species and clades? [See Valentine (1983) and Valentine et al. (2008) for a view of diversity-dependent factors that would favor this explanation.] Or do some regions incur a ââdiversi- fication debt,ââ a more positive analog to the extinction debt sometimes inferred for modern biotas squeezed into refugia too small to accommodate their present richness? The rapid recovery of diversity in the Caribbean, which evidently suffered more severely than the tropical eastern Pacific just on the other side of the Panama Isthmus (Todd et al., 2002), suggests that at least some Plio-Pleistocene extinctions involve diversification debts rather than diversity resettings. [By this argument, the anomalously low diversity of the Southeast Pacific molluscan fauna (e.g., Rivadeneira and Marquet, 2007) is a transient effect rather than a permanent biogeographic feature, attributable perhaps to the lower rates of diversity accumulation in extratropical regions; cf. Jablonski et al. (2006).] However, the spatially explicit form of the fundamental macroevolutionary equation shows that
Extinction and the Spatial Dynamics of Biodiversityâ /â 199 regional diversity can accrue either by in situ origination or invasion (immigration). The high present-day rate of anthropogenic introductions in marine systems will likely outstrip regional evolutionary recovery by an order of magnitude or more, exacerbated of course by other anthropogenic stresses [e.g., Ruiz et al. (2000), Jackson et al. (2001), LÃ¶tze et al. (2006)]. This paleontological perspective on regional marine biodiversity adds another element to the urgency of slowing anthropogenic homogenization of marine biotas, if the diversification debt of some regions indeed makes them more susceptible to invasions (and see below). RECOVERIES ARE ALSO SPATIALLY COMPLEX Most research on spatial dynamics has focused on extinctions, but evi- dence is accumulating for a spatial component to recoveries as well. The raw evolutionary material that survives the mass extinction filter is crucial in shaping the postextinction world. However, the evolutionary novelties and the ecological restructuring that emerge in the postextinction interval, including the little-appreciated process of sorting survivors into winners and losers (Jablonski, 2002), may be as important as the extinction filter in determining the long-term trajectory of individual clades. Returning to the spatially explicit form of the fundamental macroevo- lutionary equation, we can partition a regional biota after, e.g., the K-T mass extinction, into survivors of the event, locally evolved new taxa, and invaders. The four regions with the best marine molluscan records in the 10 Myr after the K-T event differ significantly in their recovery dynamics (Jablonski, 1998). For example, the North American biota was much more subject to invasions during the early stages of the post-Cretaceous recov- ery, and contains significantly fewer novel taxa, than the other regions. However, the survivor components (the regional extinction intensities) are indistinguishable among regions (Fig. 10.3). Spatial heterogeneity has recently being detected for other recoveries as well. For example, after the end-Ordovician extinction â445 Myr ago, Baltica, a continental plate centered â30Â° south, had the lowest extinction intensity and the lowest invasion intensity among marine invertebrates, whereas Laurentia, which was straddling the equator, showed a tighter bunching of the three faunal components (Krug and Patzkowsky, 2007). Spatial heterogeneity in recovery from the huge end-Permian extinction has been reported for some groups [e.g., brachiopods and bivalves (Chen et al., 2005; Bonuso and Bottjer, 2008)] but not others [e.g., ammonoids (McGowan, 2005), but see Brayard et al. (2006)]. The K-T example is striking in its failure to show the inverse rela- tion between local survivorship and invasion that is generally expected [e.g., Fridley et al. (2007)] and is observed for the end-Ordovician. One
200â /â David Jablonski Survival Origination Invasion Percent of Regional Diversity 70 60 50 40 30 20 10 0 N America N Europe N Africa Pakistan FIGURE 10.3â Regional variation in molluscan recovery dynamics after the K-T extinction. Four postextinction biotas are partitioned into their surviving, newly zpq9990837740003.g.eps originating, or invading components. The Gulf and Atlantic Coastal Plain was subject to significantly more intense invasion than the other three regions (95% binomial confidence intervals). See Jablonski (1998) for details. can speculate that this difference is attributable to the greater proxim- ity of North America to the K-T impact site in Yucatan, but how does that proximity translate to great invasions without imposing exceptional extinction? Qualitative, rather than quantitative, losses might account for the greater invasibilty of North America, but this hypothesis has not been tested. Alternatively, the correlation between extinction intensity and subsequent invasion may begin to break down above some threshold extinction level. One clue may come from another spatial difference: the short-lived evolutionary burst of a few taxa in North America [ââbloom taxaââ of Hansen (1988)] that is evidently absent in northern Europe, North Africa, or Pakistan and India (Jablonski, 1998). North Americaâs bloom taxa and the invasion pulse are almost certainly linked, presumably indi- cating a more profound ecological and evolutionary disturbance in North America than elsewhere, but this needs to be examined more closely. Whatever the ultimate cause, the fossil record pinpoints a theoretically interesting but pragmatically disquieting gap in our understanding of the extinctionâinvasion relationship. Given the current acceleration of both processes, and the pressure to establish reserves for remaining biodiver- sity, this relationship deserves more attention.
Extinction and the Spatial Dynamics of Biodiversityâ /â 201 EXTINCTION INFLUENCES SPATIAL DYNAMICS ACROSS LATITUDE Extinctions apparently promote not only invasion but evolutionary diversification in the fossil record, the classic case being the impressive radiation of mammals after the demise of the (nonavian) dinosaurs and other marine and terrestrial vertebrates at the end of the Cretaceous [e.g., Alroy (1999, 2000) and Cifelli and Gordon (2007); for general discussions, see Erwin (2001) and Jablonski (2001, 2005, 2008)]. These macroevolu- tionary observations are often seen as two sides of the same coin, as they intersect nicely with ecological work on the potential for incumbency or priority effects to resist extinction or damp diversification (Jablonski, 2007, 2008). The three most prevalent explanations for both invasions and diversifications today and in the geologic past are (i) extinction or at least suppression of incumbents, already mentioned, (ii) superior competitive ability of the invaders, not least owing to their escape from their own competitors, predators, and pathogens when they enter a new area (Sax et al., 2007) (although this may be a transient effect and therefore less likely to play a macroevolutionary role), and (iii) changing climatic and other environmental conditions, such as those that drove the extensive invasions and reshuffling of Pleistocene biotas. One of the most pervasive spatial patterns of invasions, seemingly independent of mass extinction events, underlies the marine latitudinal diversity gradient, wherein morphologies, species, and higher taxa are richest in the tropics and decline toward the poles. Although the gradi- ent has been known for a long time and is documented for many groups and regions, the processes underlying this pervasive biodiversity pattern remain poorly understood (Hillebrand, 2004; Mittelbach et al., 2007). The ââout of the tropicsââ model for the marine gradient has taxa preferentially originating in the tropics, and then expanding their latitudinal ranges over time without actually abandoning their tropical cradle (Jablonski et al., 2006). The tropics are thus a diversity source, containing both young and old taxa, which accumulate to high richness. The poles are a diversity sink, mainly containing older taxa that have moved in from lower latitudes, and the temperate zones have intermediate richness and taxon ages, at least for the marine invertebrates where direct fossil evidence is available. If this model is generally true, then invasion is a basic factor in the latitudinal deployment of life on Earth. The out-of-the-tropics model is strongly supported in the marine bivalve fossil record. For each of three time slices (Pleistocene, Pliocene, and late Miocene), roughly twice as many bivalve genera first occur in the tropics than in higher latitudes (Jablonski et al., 2006). Because the extra- tropical fossil record is far better sampled than that of the tropics (Allison and Briggs, 1993; Jackson and Johnson, 2001; Bush and Bambach, 2004;
202â /â David Jablonski Jablonski et al., 2006; Valentine et al., 2006), the tropical values must be underestimates of their true origination rates, and the extratropical values must be overestimates. Further, most of the genera that first appeared in the tropics over the past 11 Myr have since spread to higher latitudes. This dynamic accounts for the striking inverse relation between diversity in a latitudinal bin and the median age of the genera in that bin: most of the geologically old genera at high latitudes also occur in the tropics, but the young genera are concentrated at low latitudes, decreasing the low- latitude median value significantly. The number of taxa that expand out of the tropics is impressive, par- ticularly given that these clades are invading new climate zones, travers- ing a gradient of increasing physical challenges for most taxa. Further, these extratropical expansions occurred in the face of progressive global refrigeration, culminating in the full-blown glacial cycles of the Pleisto- cene. However, while clades regularly left the tropics, few, if any, of the analyzed cohort have expanded above â50Â° north or south latitude. If the relative invasibility of the temperate and polar zones over the past 11 Myr was underlain by regional variation in background extinction intensity (E), we would expect, not the usual two-bin model, tropical E < extratropical E, or low-latitude E < polar E (Goldberg et al., 2005; Jablonski et al., 2006; Roy and Goldberg, 2007), or the monotonic latitudinal trend in extinction rates assumed by many others, but a hump-shaped pattern with an extinction maximum at midlatitudes. A preliminary test of these alternatives did find a humped extinction pattern with latitude for Northern Hemisphere bivalve genera in the latest Cenozoic, with global plus regional extinction totaling â9% in the tropics, â20% in the temperate zone, and â12% in the Arctic (Valentine et al., 2008). This result suggests that the temperate zones are invasible on geologi- cal timescales because they suffer the highest extinction rates, at least in global climate states approaching our own. Thus, even if climate does not directly set standing diversity, its fluctuations, which are greatest both in temperate latitudes today and during Pleistocene climate swings [e.g., Jansson and Dynesius (2002), Ravelo et al. (2004), and Lyle et al. (2008)], may set the pattern of extinction intensities. The data are not yet sufficient to study these dynamics in detail, but the relation between midlatitude thermal variability (which coincides with fluctuations in many additional factors) and extinction patterns clearly deserve further scrutiny. The poles are doubtless demanding places to live, but taxa evolve to cope with the challenges; Valentine and colleagues (Valentine, 1983; Valentine et al., 2008) suggest they do this by becoming highly generalized trophically and argue that these broad niche dimensions are what tend to block invasions and allow them to weather glacial episodes subtidally, as they avoid seasonal extremes today. In any case, invasion resistance is apparently not a func-
Extinction and the Spatial Dynamics of Biodiversityâ /â 203 tion of standing diversity alone, but of regional extinction rates, suggest- ing a significant role for incumbency [see also Vermeij (2005)]. The polar mollusks may ultimately have come out of the tropics as well, but if so this must have occurred before the 11-Myr window presently available (which would be consistent with the much older genus ages seen at the poles). The interplay of extinction, origination, and immigration is complex, and of course it need not be in a steady state. Terrestrial animals may well follow a different dynamic (Hawkins et al., 2007b; Weir and Schluter, 2007; Wiens, 2007a). Marine organisms can move down the continental shelf when ice forms at the surface, but terres- trial animals, plants, and fungi do not have that luxury when confronted with a kilometer-thick ice sheet. High-latitude extinction and recoloniza- tion are thus almost certainly more important factors on land. Whatever the spatial dynamic near the poles, however, the tropics appear to be a crucial reservoir for biodiversity, with a subset of low-latitude clades expanding out of the tropics over geological timescales. Although this pat- tern is most readily detected in the shallowest part of the geologic record, thus falling entirely within times away from the major mass extinctions, some evidence suggests that postextinction recoveries are also fueled by the tropics, on land (Kerp et al., 2006) and in the oceans (Jablonski, 1998; Brayard et al., 2006; Krug and Patzkowsky, 2007). The tropics thus appear to be key to the generation and maintenance of global biodiversity across a wide range of boundary conditions. INTEGRATING PALEONTOLOGICAL AND NEONTOLOGICAL PERSPECTIVES I have touched on four spatial aspects of ancient extinctions that should be integrated with theoretical and applied approaches to the present-day biota. The fossil record amply demonstrates that the spatial fabric of extinction has profoundly shaped the biosphere. First, broad geographic range probably always buffers clades from extinction, but it becomes most important and clear-cut as the suite of other factors that enhance species and genus survival during normal times become ineffec- tive. It is not yet clear whether the selectivity regime changes steadily with increasing extinction intensity or as a step function (Jablonski, 2005). More intense extinctions may tend to be less selective, which might explain the failure of intrinsic factors to predict extinction risk in the some of the most heavily stressed elements of the modern biota, such as freshwater fishes, amphibians, and Australian marsupials (Duncan and Lockwood, 2001; Fisher et al., 2003). This shift to a strong spatial component in survivorship during major extinction events greatly increases the likelihood of hitchhiking effects.
204â /â David Jablonski Organismic traits can rise or fall according to the strength of their link- age to broad geographic range or other factors promoting survivorship through those bottlenecks, lending a highly stochastic element to the expansion or demise of individual adaptations or clades. Thinking about the present day, these linkages are unlikely to promote factors beneficial to, or even desirable for, humans or the ecosystems they hope to conserve (see Jackson, Chapter 1, this volume). Given that narrow-ranging genera cannot have wide-ranging species, the net effect must be to deplete spe- cialists in favor of weedy generalists, but this pattern can be ameliorated by survival of clades whose broad ranges arise from the far-flung deploy- ment of individually localized species. Second, the fossil record is rich in regional extinction events of inter- mediate intensities, and these can provide insights into present-day bio- diversity issues. For example, the Cenozoic history of todayâs biodiversity hotspots and coldspots (relative to expectations for their latitudes, for example) may help to predict the potential of these regions to accom- modate further diversification, or alternatively to be subject to biotic invasions. Third, the rules of successful recovery are poorly known, but are important for our understanding of both the larger outlines of the his- tory of life and the future of modern diversity. The inordinate production of evolutionary novelties during recoveries suggests that postextinction dynamics do not simply involve an immediate return to business as usual. At the same time the spatial heterogeneity of recoveries, with significant invasions driving some of the regional patterns, requires a more careful look at the dynamics if we want to avoid biotic homogenization even after the reduction of the pressures on the modern biota. This could be another highly fruitful area at the intersection of paleontology and conservation biology. Fourth, invasion has always been an evolutionary fact of life (Vermeij, 2005), even across biogeographic barriers and against climate gradients. The out-of-the-tropics model suggests an evolutionary approach to model- ing biotic responses to future climate changes and attests to the evolution- ary consequences of the stresses on tropical biotas today. If the tropics are the engine of global biodiversity, then driving tropical populations into extinction will have a global effect, by cutting off the primary source of new taxa for all latitudes. Further, if invasibility is more closely tied to extinction than to diversity per se, then there is the possibility of a reversal of the diversity flow, increasing the influx of invaders from higher lati- tudes. A tropical diversity crisis, now or in the geologic past, has profound long-term evolutionary consequences at a truly global scale. Simply comparing the magnitude of the extinction occurring today, which is undoubtedly severe, with ancient intensities is not the most fruit-
Extinction and the Spatial Dynamics of Biodiversityâ /â 205 ful way to draw on the insights of the fossil record, or catalyze integrative research. A better approach might be to recognize that present-day extinc- tions have many drivers, and then to test for common patterns of selectiv- ity on that basis: partitioning present-day extinction mechanisms should permit a clearer application of insights from the fossil record. Extinction selectivity probably does vary with driving mechanism to some extent. In birds, for example, habitat loss preferentially removes specialized and small-bodied taxa (but does not select on generation time), whereas exploi- tation and introduced predators preferentially remove large-bodied and long-lived taxa (Owens and Bennett, 2000). Such systematic variations in selectivity help explain the apparent contradictions in and among analyses of present-day extinction risk [e.g., Purvis et al. (2000b) and FrÃ©ville et al. (2007)], and similar arguments can be made for differences among ancient extinctions as well. As several authors have noted, extinction drivers have probably compounded over human history, with exploitation perhaps the most important in early phases, species invasions rising in frequency with the era of European exploration, and finally habitat alteration on a global scale accelerating with increased human population pressure, pollution, and climate change [e.g., Purvis et al. (2000b) and Jackson et al. (2001)]. Many paleontological perturbations are probably most analogous to present-day habitat loss and could be explored in comparative fashion on that basis. Others will more closely correspond to the introduction of enemies, as when provinces collide or novel predation mechanisms evolve (see Barnosky, Chapter 12, this volume). The particular combina- tion of pressures seen today may be unique, just as they may have been for the K-T or end-Permian events. For example, today the ordinary biotic response to climate change (range translocation) is disallowed or at least severely curtailed over much of the planet owing to occupation or con- version of suitable habitat or migration corridors by humans and their artifacts. The unique combination of forces behind each major extinction puts a premium on focusing on first principles rather than extinction- specific patterns, underscoring the need for integrative research. It also underscores the need to take a hard look at the roles of incumbency and hitchhiking effects, to separate large-scale artifacts or byproducts from the underlying drivers. Extinction thresholds presumably exist for todayâs biota, beyond which whole systems collapse and most selectivity factors drop out, as seen for major events of the geologic past. Identifying such thresholds among environments, clades, and regions using fossil data, as another basis for avoiding them in the future, would be a valuable undertaking. More generally, spatially explicit approaches to the fossil record have great potential for new insights into diversity dynamics, not just in the geologic past, but in the present day as well. The integration of paleon-
206â /â David Jablonski tological and neontological insights takes on special urgency with the acceleration of extinction rates in the modern world, and the incorporation of the spatial dimension offers a powerful vehicle for that integration. ACKNOWLEDGMENTS I thank John Avise, Francisco Ayala, and Stephen Hubbell for inviting me to participate in this colloquium; J. Alroy, J. C. Avise, P. G. Harnik, S. M. Kidwell, M. J. Novacek, and J. W. Valentine for manuscript reviews; K. Roy and J. W. Valentine for many fruitful discussions and collaborations; A. Z. Krug and P. G. Harnik for discussions; and P. G. Harnik for assis- tance with multiple logistic regression. The National Science Foundation, the National Aeronautics and Space Administration, and the John Simon Guggenheim Foundation supported this research and synthesis.