Hostname: page-component-586b7cd67f-2brh9 Total loading time: 0 Render date: 2024-11-23T16:22:16.101Z Has data issue: false hasContentIssue false

Stratigraphic paleobiology

Published online by Cambridge University Press:  07 March 2024

Steven M. Holland*
Affiliation:
Department of Geology, University of Georgia, Athens, Georgia 30602-2501, U.S.A.
Mark E. Patzkowsky
Affiliation:
Department of Geosciences, Pennsylvania State University, University Park, Pennsylvania 16802-2714, U.S.A.
Katharine M. Loughney
Affiliation:
Department of Earth and Environmental Sciences and Watershed Studies Institute, Murray State University, Murray, Kentucky 42071, U.S.A.
*
Corresponding author: Steven M. Holland; Email: [email protected]

Abstract

Stratigraphic paleobiology uses a modern understanding of the construction of the stratigraphic record—from beds to depositional sequences to sedimentary basins—to interpret patterns and guide sampling strategies in the fossil record. Over the past 25 years, its principles have been established primarily through forward numerical modeling, originally in shallow-marine systems and more recently in nonmarine systems. Predictions of these models have been tested through outcrop-scale and basin-scale field studies, which have also revealed new insights. At multi-basin and global scales, understanding the joint development of the biotic and sedimentary records has come largely from macrostratigraphy, the analysis of gap-bound packages of sedimentary rock. Here, we present recent advances in six major areas of stratigraphic paleobiology, including critical tests in the Po Plain of Italy, mass extinctions and recoveries, contrasts of shallow-marine and nonmarine systems, the interrelationships of habitats and stratigraphic architecture, large-scale stratigraphic architecture, and the assembly of regional ecosystems. We highlight the potential for future research that applies stratigraphic paleobiological concepts to studies of climate change, geochemistry, phylogenetics, and the large-scale structure of the fossil record. We conclude with the need for more stratigraphic thinking in paleobiology.

Type
Invited Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press on behalf of Paleontological Society

Non-technical Summary

Stratigraphic paleobiology uses a modern understanding of the construction of the stratigraphic record—from beds to depositional sequences to sedimentary basins—to interpret patterns and guide sampling strategies in the fossil record. Over the past 25 years, its principles have been established primarily through forward numerical modeling, originally in shallow-marine systems and more recently in nonmarine systems.

Introduction

Stratigraphic paleobiology is the application of modern concepts of stratigraphic accumulation—such as event deposition and the formation of beds, sequence-stratigraphic architecture, and sedimentary basin analysis—to the interpretation of the fossil record (Patzkowsky and Holland Reference Patzkowsky and Holland2012; Holland and Loughney Reference Holland and Loughney2021). It stems from the awareness that the fossil record is not simply the history of life: it is also the history of preservation, which is governed largely by the processes of stratigraphic accumulation. As a result, any interpretation of the fossil record must be grounded in an understanding of the sedimentary rock record. Without such a grounding, misinterpretations of the history of life are likely.

The principles of stratigraphic paleobiology for marine systems are well established from a large and growing series of modeling and field studies (summarized in Patzkowsky and Holland Reference Patzkowsky and Holland2012; Holland Reference Holland2016, Reference Holland2020, Reference Holland2023). Similar studies for nonmarine systems are still in their early stages (summarized in Holland and Loughney Reference Holland and Loughney2021). Stratigraphic paleobiology has delivered important implications for community paleoecology, patterns of morphological evolution, patterns in diversity, and biotic events, such as mass extinctions and biotic invasions. Stratigraphic paleobiology provides a means for interpreting all of these topics. In addition, it enables sampling strategies that allow one to distinguish patterns in the fossil record arising from biological change from those created by the processes of stratigraphic accumulation. Given the body of work and previous reviews, the goal of this paper is not a comprehensive coverage of the subject. Instead, we describe several key concepts, highlight recent advances in six important areas, outline four promising future research directions, and end by arguing for more stratigraphic thinking in paleobiology.

Key Concepts

Several core aspects of stratigraphic paleobiology center on species ecology and stratigraphic architecture. The intersection of these two components produces characteristic patterns in the stratigraphic occurrence of fossils that have broad applications to how the fossil record is interpreted.

Ecological Gradients

Although numerous factors control the ecological distribution of species, the application of gradient analysis (Whittaker Reference Whittaker1967) considerably simplifies the coordinated effects of these factors. Ordination methods such as nonmetric multidimensional scaling (NMS) and detrended correspondence analysis (DCA; Kruskal Reference Kruskal1964; Hill and Gauch Reference Hill and Gauch1980; Jongman et al. Reference Jongman, Ter Braak and Van Tongeren1995; Legendre and Legendre Reference Legendre and Legendre1998; Borcard et al. Reference Borcard, Gillet and Legendre2018) also provide a means for identifying the principal sources of variation in the ecological composition of communities (Fig. 1). For example, in shallow-marine (<200 m) settings, numerous physical and chemical parameters that affect the distribution of species vary with water depth, such as grain size, bed shear stress, oxygen, sunlight, salinity, temperature, and nutrients. As a result, water depth is the primary ecological gradient in marine benthic communities and is a proxy for a complex chemical and physical gradient that also includes species interactions (see reviews in Patzkowsky and Holland Reference Patzkowsky and Holland2012; Holland Reference Holland2023). In terrestrial settings, the variation in temperature, precipitation, and soil moisture with elevation cause it to be the principal ecological gradient (Holland and Loughney Reference Holland and Loughney2021). These gradients are ubiquitous in modern and ancient marine and nonmarine systems.

Figure 1. Coding ordination scores of samples and taxa can reveal the origin of the axes, shown with an example from the C2 sequence of the Upper Ordovician of the Cincinnati Arch, USA (data from Holland and Patzkowsky Reference Holland and Patzkowsky2007). A, Sample scores coded by lithofacies reveal that detrended correspondence analysis (DCA) axis 1 is correlated with water depth. B, Taxon scores coded by life habit demonstrate that DCA axis 2 is correlated with substrate characteristics, with burrowers associated with soft muds at low axis 2 scores and attached species associated with shelly gravels at high axis 2 scores.

Species distributions along these ecological gradients can be described with species response curves, which describe species’ abundance or probability of occurrence as a function of the position along an ecological gradient (e.g., water depth, elevation; Gauch and Whittaker Reference Gauch and Whittaker1976; Austin Reference Austin1987; Jongman et al. Reference Jongman, Ter Braak and Van Tongeren1995). Most commonly, the shape of these response curves is Gaussian or a skewed Gaussian; bimodal, multimodal, irregular, and flat response curves are far less common (Minchin Reference Minchin1989). The importance of these response curves is that the abundance and probability of occurrence will hit a peak (called peak abundance) at a particular position along a gradient (the preferred environment), and they will decrease away from this position depending on the environmental tolerance of a species. In shallow-marine systems, these response curves are sometimes surprisingly tight. For example, changes in water depth of even a few meters are accompanied by substantial changes in species abundance (Horton et al. Reference Horton, Edwards, Lloyd, Shennan and Andrews2000; Scarponi and Kowalewski Reference Scarponi and Kowalewski2004; Brown and Larina Reference Brown and Larina2019). Species in deeper-water systems would be expected to have broader tolerances, although species tolerances in offshore settings can sometimes be remarkably narrow (e.g., Holland et al. Reference Holland, Miller, Meyer and Dattilo2001).

Stratigraphic Architecture

That water depth and elevation are the principal ecological gradients takes on special significance, because water depth and elevation change systematically in the stratigraphic record through changes in the rates of accommodation and sedimentation. Accommodation, the space in which sediment can accumulate, is generated by tectonic subsidence, plus eustasy in marine systems and lake level in lacustrine systems. The combination of these drivers is what creates sequence-stratigraphic architecture (Van Wagoner et al. Reference Van Wagoner, Mitchum, Campion and Rahmanian1990; Catuneanu Reference Catuneanu2006; Catuneanu et al. Reference Catuneanu, Abreu, Bhattacharya, Blum, Dalrymple, Eriksson and Fielding2009; Neal and Abreu Reference Neal and Abreu2009), specifically the formation of subaerial unconformities and other hiatal surfaces, the lateral distribution of sedimentary facies, transgressions and regressions, and variations in sedimentation rates. In marine and coastal systems, these drivers produce systems tracts characterized by their positions, internal stacking patterns of strata, and shoreline trajectories. In nonmarine systems, systems tracts are characterized by position, channel stacking patterns, types of paleosols, and various other facies characteristics (Catuneanu Reference Catuneanu2006; Holland and Loughney Reference Holland and Loughney2021).

Although there are several approaches to defining these systems tracts (Catuneanu et al. Reference Catuneanu, Abreu, Bhattacharya, Blum, Dalrymple, Eriksson and Fielding2009), the most common model includes a lowstand systems tract (LST), transgressive systems tract (TST), highstand systems tract (HST), and falling-stage systems tract (FSST; Fig. 2A). In nonmarine systems, a low-aggradation systems tract (LAST) and high-aggradation systems tract (HAST) are the most widely used (Martinsen et al. Reference Martinsen, Ryseth, Helland-Hansen, Flesche, Torkildsen and Idil1999; Catuneanu Reference Catuneanu2006; Holland and Loughney Reference Holland and Loughney2021; Loughney and Holland Reference Loughney and Holland2023; Rogers et al. Reference Rogers, Eberth and Ramezani2023; Fig. 2B). Systems tracts record distinct patterns of change in water depth and elevation, including progressive trends as well as surfaces of abrupt change. For example, the HST displays increasingly rapid upward shallowing punctuated by relatively minor surfaces of abrupt deepening (flooding surfaces); the LST is similar, but the upward trend slows rather than accelerates. The TST displays net upward deepening, composed in detail of major flooding surfaces of abrupt deepening that separate intervals of minor shallowing. The FSST displays net upward shallowing punctuated by surfaces of abrupt shallowing (surfaces of forced regression).

Figure 2. Stratigraphic cross sections showing systems tracts and surfaces in marine and nonmarine settings. A, Cross section is model output from the basin simulation model sedflux (Hutton and Syvitski Reference Hutton and Syvitski2008). FSST, falling-stage systems tract; HST, highstand systems tract; LST, lowstand systems tract; TST, transgressive systems tract. B, Cross section is model output from the basin simulation model strataR (Holland Reference Holland2022a). Black lines are evenly spaced timelines, which can be used to infer aggradation rates. HAST, high-accommodation systems tract; LAST, low-accommodation systems tract.

Stratigraphic architecture varies markedly depending on the spatial pattern and temporal changes in subsidence rate, the rates and timescales of eustatic sea-level change, and the nature of how sediment is produced and transported, itself reflecting climate, topography, source area, and biogenic production (Catuneanu Reference Catuneanu2006). Understanding these controls has been a primary and ongoing focus of research in sedimentary geology for decades. For example, stratigraphic architecture can differ on opposite sides of a sedimentary basin depending on subsidence rates and the relative contributions of siliciclastic and carbonate sediment (Tomašových et al. 2022).

This is significant for the fossil record, because the stratigraphic record in which it is housed is highly structured and nonuniform, although the structure varies for all the reasons mentioned earlier. Sedimentation rates vary widely and predictably. Water depth and elevation change continuously, progressively in some cases, and abruptly in others. Subaerial unconformities occur in predictable places laterally and vertically. As a result, time is not preserved uniformly in strata, and preserved habitats are always in flux. Moreover, the stratigraphic record varies markedly and predictably along depositional dip (Fig. 3). This view of stratigraphy radically changes how we approach the fossil record, because constancy in sedimentation rates and ecological settings can rarely be assumed. The probability of preservation of fossil species is therefore rarely uniform but is instead in constant flux. Any interpretation of the fossil record must be based on this reality. Moreover, these changes are not random but are instead highly structured.

Figure 3. Depositional strike and dip, illustrated with the coast of Egypt at the Nile Delta. Map from Google Earth Pro. Depositionally updip areas are toward the sediment source (typically landward), and depositionally downdip areas are distal to the sediment source and therefore more commonly sediment starved.

Recent Advances

Tests of Key Concepts in the Po Plain of Italy

Paleobiology gets increasingly more challenging into deep time for two reasons. The first is decreasing geochronological precision in older rocks. For example, typical uncertainties in Pleistocene strata using 14C and 230Th are less than ±2 kyr, and lower than ±200 yr for dates younger than 15 ka (e.g., Leigh Reference Leigh, Greenberg and Collins2016; Railsback et al. Reference Railsback, Brook, Liang, Voarintsoa, Cheng and Edwards2018). In contrast, U–Pb age uncertainties in the Late Ordovician range from ±170 kyr to ±340 kyr (Ling et al. Reference Ling, Zhan, Wang, Wang, Amelin, Tang and Liu2019), up to 20 times the duration of the Holocene. The second issue is that the ecology of species is increasingly poorly known into deep time. In the younger parts of the fossil record, fossil species may still be extant, and their ecology today is commonly a reasonable proxy for their ecology in the fossil record. In deeper time, proxies for extinct species are at higher taxonomic levels, and niche evolution (see Blois et al. 2024) makes inferences even more tenuous.

An extraordinary series of studies in late Pleistocene to Holocene deposits of the Po River delta in northern Italy avoids both complications through its geologically young age (Scarponi and Kowalewski Reference Scarponi and Kowalewski2004, Reference Scarponi and Kowalewski2007; Scarponi et al. Reference Scarponi, Kaufman, Amorosi and Kowalewski2013, Reference Scarponi, Azzarone, Kusnerik, Amorosi, Bohacs, Drexler and Kowalewski2017; Wittmer et al. Reference Wittmer, Dexter, Scarponi, Amorosi and Kowalewski2014; Huntley and Scarponi Reference Huntley and Scarponi2015; Kowalewski et al. Reference Kowalewski, Wittmer, Dexter, Amorosi and Scarponi2015; Nawrot et al. Reference Nawrot, Scarponi, Azzarone, Dexter, Kusnerik, Wittmer, Amorosi and Kowalewski2018). These studies provide critical tests of many aspects of stratigraphic paleobiology.

The studies are built on a well-delineated sequence-stratigraphic framework assembled from cores arrayed across depositional dip and depositional strike (Amorosi et al. Reference Amorosi, Centineo, Colalongo, Paasini, Sarti and Vaiani2003, Reference Amorosi, Colalongo, Fiorini, Fusco, Pasini, Vaiani and Sarti2004). Within these cores, mollusks were counted at multiple horizons (Scarponi and Kowalewski Reference Scarponi and Kowalewski2004). Ultimately, the dataset included 16 cores and 131,780 specimens in 234 species (Wittmer et al. Reference Wittmer, Dexter, Scarponi, Amorosi and Kowalewski2014; Kowalewski et al. Reference Kowalewski, Wittmer, Dexter, Amorosi and Scarponi2015), although the original dataset was less than a fifth of that. Counts of foraminifera, ostracods, and pollen were also analyzed in some studies (Amorosi et al. Reference Amorosi, Colalongo, Fiorini, Fusco, Pasini, Vaiani and Sarti2004, Reference Amorosi, Rossi, Scarponi, Vaiani and Ghosh2014). 14C-age dates were obtained throughout, supplemented by dated pollen correlations to other European sites. Selected shells were dated individually with amino acid racemization and calibrated with 14C (Scarponi et al. Reference Scarponi, Kaufman, Amorosi and Kowalewski2013) with precisions of ±100–300 yr. Because most of these species are extant in the Mediterranean, their modern ecological distributions relative to substrate and water depth are well understood (Pérès and Picard Reference Pérès and Picard1964; Pérès Reference Pérès1967; Dominici and Scarponi Reference Dominici and Scarponi2020) and accessible through a database operated by the New Technologies, Energy, and Environmental Agency (ENEA; Wittmer et al. Reference Wittmer, Dexter, Scarponi, Amorosi and Kowalewski2014).

One of the original goals of these studies was to test interpretations made from ordinations of fossil assemblages, specifically that axis 1 of some ordination methods is commonly correlated with water depth (e.g., Holland et al. Reference Holland, Miller, Meyer and Dattilo2001; Miller et al. Reference Miller, Holland, Meyer and Dattilo2001). The most common of these methods are DCA and NMS. The significance of axis 1 is that it reflects the greatest source of variation in the composition of fossil assemblages.

Salient Results

The Po Plain studies reach several significant conclusions, and they are encouraging news for stratigraphic paleobiology in deep time. First, indirect ordination techniques like DCA perform just as well as direct ordination techniques in their ability to recover ecological gradients (Wittmer et al. Reference Wittmer, Dexter, Scarponi, Amorosi and Kowalewski2014). This is important in deep time, where independent estimates of water-depth preferences of taxa are unavailable, rendering direct ordination impossible. This also confirms interpretations of ancient fossil assemblages, where a water-depth gradient is inferred based on lithofacies (e.g., Cisne and Rabe Reference Cisne and Rabe1978) or sequence-stratigraphic architecture (e.g., Miller et al. Reference Miller, Holland, Meyer and Dattilo2001; Holland and Patzkowsky Reference Holland and Patzkowsky2007). Indirect ordination can provide depth estimates with uncertainties as small as ±3 m, suggesting that these methods are generally useful (Scarponi and Kowalewski Reference Scarponi and Kowalewski2004). These estimates allow one to track changes in water depth through time and lateral changes in water depth (Scarponi and Kowalewski Reference Scarponi and Kowalewski2004).

Second, these studies validate interpretations based on sequence-stratigraphic architecture alone and predictions of time resolution in the fossil record. For example, the Po Plain studies demonstrate that the TST indeed records upward deepening and that the HST records upward shallowing (Scarponi and Kowalewski Reference Scarponi and Kowalewski2004). Similarly, these methods demonstrate that upward within a TST, time averaging increases, the frequency of depositional events decreases, and net accumulation rates decrease; HSTs show the opposite trends (Scarponi et al. Reference Scarponi, Kaufman, Amorosi and Kowalewski2013). Notably, patterns of time averaging differ on the opposite side of the Po Basin, where mixed carbonate–siliciclastic deposition dominates, underscoring the differences in carbonate and siliciclastic systems (Tomašových et al. 2022; Belanger and Bapst Reference Belanger and Bapst2023). Moreover, ordination-based depth proxies can allow sequence-stratigraphic surfaces and vertical water-depth trends to be recognized in lithologically uniform facies (Holland et al. Reference Holland, Miller, Meyer and Dattilo2001; Amorosi et al. Reference Amorosi, Rossi, Scarponi, Vaiani and Ghosh2014). The Po Plain studies also demonstrate that the residence time of shells near the sediment surface increases offshore, with taphonomic degradation taking place over centennial scales proximally and millennial scales distally (Scarponi et al. Reference Scarponi, Azzarone, Kusnerik, Amorosi, Bohacs, Drexler and Kowalewski2017).

Third, these studies demonstrate the need for sampling along depositional dip (Fig. 3). In the original study, which analyzed only three cores along depositional strike, the authors found that late TST samples were far more uniform in composition than those of the early TST or the HST (Scarponi and Kowalewski Reference Scarponi and Kowalewski2004). When cores along depositional dip were added, this pattern disappeared, because a greater range of late TST environments could be captured, increasing the apparent variability of the communities (Wittmer et al. Reference Wittmer, Dexter, Scarponi, Amorosi and Kowalewski2014). Without this sampling along depositional dip, one might misinterpret this pattern as increasing community variability through time. Other studies have also demonstrated the extent to which time averaging varies along depositional dip (Ritter et al. Reference Ritter, Erthal, Kosnik, Kowalewski, Coimbra, Caron and Kaufman2023) and that along-strike variation can be as great as along-dip variation (Tomašových et al. Reference Tomašových, Kidwell and Barber2016).

Fourth, these studies demonstrate the substantial resilience of communities to high-amplitude sea-level changes. Mollusk communities are dominated by the same species through two glacioeustatic cycles, with similar water-depth gradients in alpha diversity, dominance, species composition, and specimen abundance (Kowalewski et al. Reference Kowalewski, Wittmer, Dexter, Amorosi and Scarponi2015). Even so, communities have their limits, and a comparison of modern-day assemblages to Pleistocene and early Holocene assemblages reveals stark differences, underscoring the intensity of anthropogenic impacts on shallow-marine ecosystems (Kowalewski et al. Reference Kowalewski, Wittmer, Dexter, Amorosi and Scarponi2015). Hypoxia in particular can have a substantial impact on the structure of ecological gradients (Tomašových et al. Reference Tomašových, Albano, Fuksi, Gallmetzer, Haselmair, Kowalewski, Nawrot, Nerlović, Scarponi and Zuschin2020), which may have particular importance in the application of stratigraphic paleobiology to mass extinctions.

Finally, reconstructing the tempo of extinction from patterns of last occurrences has been an important implication of stratigraphic paleobiology (Holland and Patzkowsky Reference Holland and Patzkowsky2015; Holland Reference Holland2020; Zimmt et al. Reference Zimmt, Holland, Finnegan and Marshall2021). By leveraging the Po Plain dataset, Nawrot et al. (Reference Nawrot, Scarponi, Azzarone, Dexter, Kusnerik, Wittmer, Amorosi and Kowalewski2018) show that if a mass extinction occurred today, the fossil record of the Po Plain would produce a highly misleading picture of the extinction (Nawrot et al. Reference Nawrot, Scarponi, Azzarone, Dexter, Kusnerik, Wittmer, Amorosi and Kowalewski2018). Specifically, stratigraphically produced clusters of last occurrences arise that do not coincide with times of elevated extinction (Fig. 4). Their study has sobering implications for interpreting nearly all ancient extinction intervals.

Figure 4. Patterns of fossil occurrence through two transgressive–regressive sequences. The two sequences each contain 10 high-frequency parasequences, upward-shallowing cycles bounded by flooding surfaces across which water depth increases rapidly. Fifty species having the same peak abundance and depth tolerance but varying in preferred depth are shown sorted by their preferred depths. Dots depict fossil occurrences; solid lines show fossil ranges; and thin dashed gray lines show when species are extant (although possibly not preserved). Fossil occurrences systematically track changes in water depth, and major flooding surfaces (blue lines) have clusters of first and last occurrences. Note that the interval of elevated extinction does not have an obviously greater number of last occurrences. Most species that go extinct during the period of elevated extinction have their last occurrences at two major flooding surfaces immediately below the window of elevated extinction. Other last occurrences are spread through the elevated extinction interval and are not distinguishable from last occurrences caused by changing water depths. The full magnitude of the extinction would not be apparent until the next transgressive–regressive sequence (not shown). The timing of the extinction can be known only by studying correlative columns along the depositional dip.

Mass Extinctions and Recoveries

Determining the drivers of mass extinction and the timing of recoveries hinges on a careful and nuanced reading of the stratigraphic and fossil records to infer cause and effect. A careful reading of the records is crucial, because as the geologic timescale is refined and calibrated to a precision of a few tens of thousands of years, it is possible to estimate rates of processes at timescales thought to be impossible not that long ago. For these reasons, many studies collect and interpret data at the outcrop scale, commonly from single stratigraphic columns. Some have even advocated that mass extinctions are best studied in a single stratigraphic column (Lucas Reference Lucas2017). Yet it is precisely at the outcrop scale where interpreting the timing, tempo, and causes of mass extinction and recovery are potentially most fraught. Because stratigraphic paleobiological modeling studies of mass extinctions are so far limited to marine settings, our discussion will focus on the marine record.

Models of the Fossil Record

Owing to the distribution of species along environmental gradients and their characteristic unimodal response curves, coupled with the systematic stratigraphic changes in water depth and accumulation rate, species occurrences are markedly structured and non-uniform. Notably, the first and last occurrences of fossils will be clustered at certain surfaces, such as subaerial unconformities, marine flooding surfaces, and surfaces of forced regression (Fig. 4). These clusters of first and last occurrences will form even if extinction rate is not elevated above background levels (Holland and Patzkowsky Reference Holland and Patzkowsky2015; Nawrot et al. Reference Nawrot, Scarponi, Azzarone, Dexter, Kusnerik, Wittmer, Amorosi and Kowalewski2018; Holland Reference Holland2020). This nonrandom distribution of fossil occurrences presents a particular challenge for paleontologists, because it means that the fossil record cannot be taken at face value. Models of the fossil record have been essential for guiding the interpretation of patterns of first and last occurrences during intervals of mass extinction and recovery.

A fossil record of species occurrences can be simulated by combining the output from three different models (Holland Reference Holland1995, Reference Holland2020, Reference Holland2023; Holland and Patzkowsky Reference Holland and Patzkowsky1999, Reference Holland and Patzkowsky2002, Reference Holland and Patzkowsky2015). First, a random-branching model of evolution with specified origination and extinction rates produces a set of species with known times of origination and extinction. These rates may be uniform through time or may be allowed to vary to simulate a mass extinction. Second, a Gaussian species response model assigns random values of preferred water depth, depth tolerance, and peak abundance to species. These three parameters determine the shape of the species response curve, and the curve's height specifies the probability of the species occurring at each point along its distribution. Third, a basin-scale sediment-accumulation model, including sea-level fluctuations, tectonic subsidence, and sediment dispersal and accumulation, is used to simulate a stratigraphic record. Several models produce realistic stratigraphic geometries, including marine systems tracts, as well as important stratigraphic surfaces such as subaerial unconformities, surfaces of forced regression, transgressive surfaces, and marine flooding surfaces (e.g., Hutton and Syvitski Reference Hutton and Syvitski2008). By combining these three models, a stratigraphic column can be sampled at any point in the basin, and for each time step, species occurrences are determined by their water depth–dependent probability of collection.

A primary finding of these modeling studies is that clusters of first and last occurrences are common at particular stratigraphic surfaces, even when origination and extinction rates are constant. If the stratigraphic context of these clusters is not recognized, they are likely to be misinterpreted as pulses of extinction and origination. For example, a cluster of last occurrences of shallow-water taxa is likely to form at major flooding surfaces within the TST, partly from the associated water-depth change and partly from the low sediment-accumulation rate. Similar clusters of last occurrences of deeper-water taxa can occur at surfaces of forced regression within the FSST. Clusters of last occurrences also typically occur at subaerial unconformities.

When mass extinctions are added to these simulations, clusters of first and last occurrences form as expected at sequence-stratigraphic surfaces, although their magnitude is greater than in simulations lacking elevated extinction rates (Holland and Patzkowsky Reference Holland and Patzkowsky2015; Holland Reference Holland2020). These clusters may lie within the interval of elevated extinction rates but may also precede it. Multiple clusters are typically formed and are likely to be misinterpreted as pulsed mass extinction, particularly if the stratigraphic context is not considered. Within the stratigraphic interval of elevated extinction rates, the numbers of last occurrences are generally elevated above background levels. To recognize such elevated numbers, one must compare this interval with another stratigraphic interval not associated with a mass extinction.

Mass Extinction Records Are Consistent with Models

Extensive reviews of the literature on the stratigraphic paleobiology of mass extinctions (see Holland and Patzkowsky Reference Holland and Patzkowsky2015; Holland Reference Holland2020) indicate that nearly all mass extinctions, with one notable exception, are associated with clusters of last occurrences at predictable stratigraphic surfaces, making it difficult to interpret the causes, timing, and duration of the mass extinctions. In many cases, such clusters of last occurrences are associated with major flooding surfaces, and multiple clusters can form at the major flooding surfaces that constitute TSTs (Fig. 4; e.g., Palmer Reference Palmer1984; Farabegoli et al. Reference Farabegoli, Perri and Posenato2007). For some extinctions, particularly those that have been studied in depositionally downdip areas, clusters are associated with surfaces of forced regression in the FSST (e.g., Harper et al. Reference Harper, Hammarlund and Rasmussen2014). Clusters have also been reported from surfaces that are subaerial unconformities (e.g., Finney et al. Reference Finney, Berry, Cooper, Ripperdan, Sweet, Jacobson, Soufiane, Achab and Noble1999). We are generally skeptical that clusters associated with sequence-stratigraphic surfaces record elevated extinction rates, particularly where they are documented from a few geographically concentrated columns. Demonstrating that they reflect elevated extinction rates would require showing that the clusters exist in depositionally updip or downdip areas where that surface is not developed. Without such a demonstration, the fossil record is consistent with a protracted period of extinction and turnover, with pulses of last occurrences being stratigraphically generated.

The only exception to these patterns is the Cretaceous/Paleogene (K/Pg) mass extinction, where the principal cluster of last occurrences does not obviously coincide with a sequence-stratigraphic surface, except in depositionally updip areas. Although a purported earlier pulse of extinction has been reported from Antarctica (Tobin Reference Tobin2017), there is strong evidence that this cluster of last occurrences lies at a flooding surface. The strong agreement between the stratigraphic expression of most mass extinctions with the modeling results suggests that most mass extinctions represent long intervals of increased extinction rate of tens to hundreds of thousands of years or more, rather than a rapid pulse extinction as represented by the K/Pg event.

Strategies for Interpretation

Knowing how stratigraphic architecture controls fossil occurrences provides a strategy for interpreting clusters of first and last occurrences in terms of mass extinction and recovery. Because clusters of last occurrences will occur just below abrupt changes in facies at major stratigraphic surfaces, the surest approach is to study additional columns depositionally updip and downdip in search of the correct facies above the stratigraphic surface where the species could occur. If the species cannot be found above the stratigraphic surface, then it supports the interpretation of extinction. If some of the species are found above the stratigraphic surface, then it suggests that stratigraphic architecture is guiding the pattern, or that the magnitude and tempo of extinction may need to be reinterpreted. For example, trilobite species in New York, USA, are found above the Ptychaspid biomere boundary where they were thought from other localities to have gone extinct (Landing et al. Reference Landing, Westrop, Kröger and English2011). This demonstrates that the last occurrence in those other localities is likely controlled by sequence-stratigraphic architecture, indicating that the biomere extinction event was more protracted than a face-value reading of the fossil record would suggest. Similarly, in a study of multiple sections along a depositional transect in New York and Pennsylvania, USA, brachiopod species thought to have gone extinct in the Late Devonian (Frasnian–Fammenian) extinction were found above the lower Kellwasser (Pier et al. Reference Pier, Brisson, Beard, Hren and Bush2021). We would argue that the extinction was more protracted than a simple reading of last occurrences would suggest and that sampling along depositional dip demonstrates its protracted nature. In the Jurassic (lower Toarcian) of eastern Spain, a cluster of last occurrences of some species is at the transgressive surface downdip but at the maximum flooding surface updip (Danise et al. Reference Danise, Clémence, Price, Murphy, Gómez and Twitchett2019). Again, this demonstrates the need to study multiple columns across depositional dip to determine the tempo of the extinction and to help constrain its cause.

In some cases, it is impossible to track the proper facies for all species up and down depositional dip (Fig. 3), because that facies is removed by erosion or buried by overlying strata. Simulations suggest that it is still possible to infer the history of extinction from clusters of last occurrences produced by stratigraphic architecture by filtering out the species whose last occurrences coincide with the basin-wide loss of their preferred facies (Zimmt et al. Reference Zimmt, Holland, Finnegan and Marshall2021). More studies that specifically evaluate stratigraphic controls on the record of extinction are needed. For example, an analysis of the Silurian Mulde event in England found that depositional environment was the strongest control on changes in conodont assemblages (Jarochowska et al. Reference Jarochowska, Ray, Röstel, Worton and Munnecke2018).

Similarly, interpretations of the tempo and timing of recoveries are also subject to the control of stratigraphic architecture on first occurrences. Because most mass extinctions are expressed as clusters of last occurrences in the late HST, FSST, LST, and earliest TST (Holland Reference Holland2020), the recoveries must begin in the latter part of the TST or the earliest HST. Alternatively, the mass extinction may record a period of faunal turnover in which species origination is elevated as the extinction progresses, a prediction supported at least by ecological neutral theory models (Hubbell Reference Hubbell2001). Moreover, if diversity declines offshore (e.g., Sepkoski Reference Sepkoski1988), the TST will preserve progressively lower-diversity habitats, and the progradation accompanying the following HST would return to shallower-water, higher-diversity depositional environments. In this situation, there would be an initial cluster of first occurrences of deep-water species followed by a long gap in the lower-diversity deeper-water facies, overlain by a string of first occurrences as shallower species first appear in the HST. If read naively as a direct record of recovery rather than as the expected stratigraphic pattern of fossil occurrences in the TST and overlying HST, the pattern would be misinterpreted as indicating a long recovery period. Many descriptions of recovery fit this pattern (summarized in Holland Reference Holland2020) and could explain the common view of long recoveries following mass extinctions. As with mass extinctions, understanding the duration of recovery intervals requires studying multiple columns along depositional dip. For example, recovery from the early Toarcian extinction seems to have occurred earlier in shallow-water sections depositionally updip compared with deeper-water settings depositionally downdip (Danise et al. Reference Danise, Clémence, Price, Murphy, Gómez and Twitchett2019). Study of additional columns may narrow the duration of the recovery interval (Atkinson et al. Reference Atkinson, Little and Dunhill2023). Some studies indicate that local ecosystem recovery can occur considerably faster than global diversity, because recovery of local diversity can be accommodated through species migration, whereas rebound of global diversity requires speciation (Johnson and Ellis Reference Johnson and Ellis2002; Christie et al. Reference Christie, Holland and Bush2013; Lyson et al. Reference Lyson, Miller, Bercovici, Weissenburger, Fuentes, Clyde and Hagadorn2019; Wilf et al. Reference Wilf, Carvalho and Stiles2023). More studies on the stratigraphic paleobiology of recoveries are needed to know what the fossil record can reveal about the resiliency of ecosystems and their rates of change following mass extinctions. Essential to future analyses are studies of stratigraphic architecture that establish a baseline of ecological conditions leading up to extinction (Clement and Tackett Reference Clement and Tackett2021).

Marine versus Nonmarine Systems

Until recently, the stratigraphic paleobiology of marine systems has been far better understood than that of nonmarine systems. This may partly reflect the substantial differences in fossil abundance (Regan et al. Reference Regan, Rogers and Holland2022) and, as a result, the generally greater ease with which large ecological datasets can be created with marine fossils. In addition, numerical models for generating hypotheses for shallow-marine systems have been well established for nearly 30 years (Holland Reference Holland1995, Reference Holland, Erwin and Wing2000, Reference Holland2020; Holland and Patzkowsky Reference Holland and Patzkowsky1999, Reference Holland and Patzkowsky2002, Reference Holland and Patzkowsky2015), whereas comparable models for fluvial systems have only recently been created (Holland and Loughney Reference Holland and Loughney2021; Holland Reference Holland2022b, Reference Holland2023). Comparable models for lacustrine, eolian, and other nonmarine systems have not yet been developed.

Bioclastic Accumulations

One exception to this disparity in how these systems have been studied is the origin of bioclastic accumulations. In both systems, fossil concentrations are controlled by the rates of bioclast input (and preservation) and the rates of sedimentation (Kidwell Reference Kidwell1986; Tomašových et al. Reference Tomašových, Fürsich and Olszewski2006a,Reference Tomašových, Fürsich and Wilmsenb). As a result, bioclastic concentrations are predicted for stratigraphic settings where rates of sedimentation are low, such as surfaces of downlap, onlap, truncation, and backlap (Kidwell Reference Kidwell, Einsele, Ricken and Seilacher1991a,Reference Kidwell, Allison and Briggsb). In marine systems, shell and bone concentrations occur in all four settings (Kidwell Reference Kidwell, Einsele, Ricken and Seilacher1991a,Reference Kidwell, Allison and Briggsb; Rogers and Kidwell Reference Rogers, Kidwell, Rogers, Eberth and Fiorillo2007). Such associations of bioclastic horizons with stratigraphically significant surfaces are much weaker in nonmarine systems, as is their relationship to the scale or duration of these surfaces (Rogers and Kidwell Reference Rogers and Kidwell2000, Reference Rogers, Kidwell, Rogers, Eberth and Fiorillo2007). Moreover, bioclastic concentrations in nonmarine rocks are most likely to occur where an erosional surface intersects underlying rock containing bioclasts (Rogers and Kidwell Reference Rogers and Kidwell2000; Rogers and Brady Reference Rogers and Brady2010).

Comparative Architecture

For shallow-marine and fluvial settings, numerical stratigraphic paleobiological models are important sources of testable hypotheses. Although shallow-marine and fluvial systems can be simulated with a similar random-branching model, they differ substantially in the simulation of species along ecological gradients and the simulation of stratigraphic architecture.

In both settings, species are distributed along ecological gradients that are expressed stratigraphically: water depth in marine settings and elevation in nonmarine settings. This distribution is a key component of stratigraphic paleobiology, and probably the most underemphasized and poorly recognized. Neither elevation nor water depth directly controls where species occur, but both gradients are highly correlated with the physical and chemical factors that affect species. In marine settings, species occurrence also varies markedly with substrate consistency and oxygen concentration, which commonly do not correlate with water depth (Holland and Patzkowsky Reference Holland and Patzkowsky2007; Scarponi et al. Reference Scarponi, Azzarone, Kusnerik, Amorosi, Bohacs, Drexler and Kowalewski2017). In nonmarine settings, species preservation varies markedly between channel and floodplain facies (Behrensmeyer and Hook Reference Behrensmeyer, Hook, Behrensmeyer, Damuth, DiMichele, Potts, Sues and Wing1992).

Architecturally, shallow-marine and nonmarine settings differ in several important ways. The shallow-marine stratigraphic record is highly cyclic, with repeated changes in water depth and sediment-accumulation rate. In many cases, this is manifested as parasequences, generally shallowing-upward cycles bounded by flooding surfaces, which are sharp contacts characterized by deeper-water facies overlying shallower-water facies (Van Wagoner et al. Reference Van Wagoner, Mitchum, Campion and Rahmanian1990; but see Catuneanu and Zecchin [Reference Catuneanu and Zecchin2020] for a critique of the parasequence concept). These are in turn arranged into parasequence sets that display net upward shallowing (progradational), deepening (retrogradational), no net depth change (aggradational), or that step seaward and downward into the basin (degradational; Van Wagoner et al. Reference Van Wagoner, Mitchum, Campion and Rahmanian1990; Neal and Abreu Reference Neal and Abreu2009). Degradationally stacked units also tend to have surfaces of forced regression, where relatively deeper-water facies are abruptly overlain by shoreface facies. Based on their stacking patterns and positions relative to one another, parasequence sets form four systems tracts within unconformity-bounded depositional sequences (Van Wagoner et al. Reference Van Wagoner, Mitchum, Campion and Rahmanian1990; Catuneanu Reference Catuneanu2006; Catuneanu et al. Reference Catuneanu, Abreu, Bhattacharya, Blum, Dalrymple, Eriksson and Fielding2009). Moreover, depositional sequences and systems tracts can comprise smaller-scale depositional sequences instead of parasequences.

The fluvial nonmarine stratigraphic record is structured differently. First, fluvial systems are divided into LASTs and HASTs; Martinsen et al. Reference Martinsen, Ryseth, Helland-Hansen, Flesche, Torkildsen and Idil1999), which divide a spectrum of possible nonmarine architectures (Holland and Loughney Reference Holland and Loughney2021; Loughney and Holland Reference Loughney and Holland2023). As the rate of accommodation increases, fluvial channel bodies switch from multistory to single-story, and they become increasingly isolated within floodplain deposits. Incised valleys are increasingly less favored, abandoned channels become more common, paleosols transition from well drained to hydromorphic, and ponds and lakes are increasingly favored (Holland and Loughney Reference Holland and Loughney2021: fig. 19). Although fluvial systems have subaerial unconformities, as do shallow-marine systems, there are no nonmarine counterparts to the flooding surfaces and surfaces of forced regression of marine systems.

Patterns of Fossil Occurrences

The combinations of stratigraphic architecture and ecological gradients unique to shallow-marine and fluvial systems produce distinct patterns in fossil occurrences. In marine systems, several types of surfaces display abrupt changes in species occurrences, including flooding surfaces, surfaces of forced regression, condensed sections, and subaerial unconformities (Holland Reference Holland1995, Reference Holland, Erwin and Wing2000). In many cases, these surfaces have distinct ecological signatures, such as the abrupt appearance of a suite of deeper-water species just above a surface characterized by the similarly abrupt disappearance of shallower-water species at flooding surfaces. This creates not only clusters of first and last occurrences in marine settings, but also abrupt changes in community composition and attributes of communities (such as richness, evenness), as well as abrupt changes in the degree of time averaging and taphonomic modification, as those are also correlated with marine ecological gradients. Such clusters can also be modified over centimeters to decimeters by sediment mixing through bioturbation, as can fossil occurrences in general (Tomašových et al. Reference Tomašových, Kidwell and Dai2023b). In addition, marine deposits contain intervals of progressive change in fossil occurrences that reflect longer-term changes in water depth that are entirely predictable from sequence-stratigraphic architecture. For example, HSTs are characterized by increasingly progradational stacking patterns, which in turn drive accelerating rates of shallowing reflected in patterns of fossil occurrences.

Fluvial systems differ in several important ways. Fluvial systems record progressive changes in paleoelevation predicted to be reflected in fossil occurrences (Holland and Loughney Reference Holland and Loughney2021; Holland Reference Holland2022b, Reference Holland2023). For example, transgression of the shore drives a progressive decrease in paleoelevation, producing a stratigraphic shift toward coastal versus inland species. Transgression also accelerates fluvial aggradation rates, favoring a shift to single-story channels encased in greater amounts of floodplain deposits. As a result, fossil occurrences will be increasingly dominated by species more likely to be preserved in floodplain settings, and the taphonomic mode will also shift toward one characteristic of floodplains. The opposite patterns will occur where regression is driven by sediment supply (i.e., normal regression, as opposed to forced regression, which is driven by a relative fall in sea level). In fluvial systems, the subaerial unconformity is the only surface type likely to display abrupt changes in fossil occurrences.

Large-Scale Structure

At the coarsest scale, shallow-marine and nonmarine systems differ in other important ways. For nonmarine systems, the 84% of land area today that is eroding (Nyberg and Howell Reference Nyberg and Howell2015) does not lie within a sedimentary basin and cannot leave a fossil record. Moreover, modern sedimentary basins overwhelmingly record elevations of 0–600 m (Holland et al. Reference Holland, Loughney and Cone2023). Of basin types that are likely to survive into deep time, only extensional basins preserve elevations up to 2000–3000 m above sea level. Piggyback and intrarc basins also span elevations up to 3000–4000 m, but because they are embedded within mountainous areas undergoing long-term uplift, they are not likely to survive into deep time. As a result, the nonmarine fossil record will systematically preserve only a particular portion of Earth's biodiversity. Large areas will leave no fossil record, including biodiversity hotspots, commonly at high elevations (Orme et al. Reference Orme, Davies, Burgess, Eigenbrod, Pickup, Olson and Webster2005). One implication of this is that many nonmarine clades may have long ghost lineages and that fossil-based divergence times might commonly be substantially shorter than those based on molecular evidence.

The marine fossil record faces a different type of selectivity at these coarsest scales. Of modern coastlines, only 28% are depositional, the others are either rocky coasts or coasts with patchy sediment unlikely to leave a long-term stratigraphic record (Nyberg and Howell Reference Nyberg and Howell2016). Many of the erosional coasts are in high northern latitudes and are likely non-depositional owing to glacial rebound. As such, the proportion of depositional coasts may have been substantially higher in non-glacial times. Regardless, some substantial proportion of coasts throughout geologic time was likely rocky with a biota adapted to hard substrates and unlikely to leave a fossil record (Johnson Reference Johnson1988; but see Johnson et al. Reference Johnson, Skinner and MacLeod1988). Limited hard-substrate biotas are known from hardgrounds and reefs in carbonate settings, from boulders and gravel at ravinement surfaces, and from shells and wood (Taylor and Wilson Reference Taylor and Wilson2003; Dominici and Forli Reference Dominici and Forli2021). Marine areas also selectively preserve a range of water depths, specifically depths of up to at most a few hundred meters, the limit to which continental lithosphere can be depressed below sea level owing to its buoyancy. Although deposition certainly occurs in the deep ocean and is readily accessed through deep-sea cores, most deep-ocean deposits are destroyed through subduction. Except for comparatively rare ophiolites and other isolated fragments, nearly all oceanic lithosphere is Early Jurassic or younger (Granot Reference Granot2016; Rowley Reference Rowley2018). As a result, the macrofossil record of the deep ocean is much sparser than that of the shallow ocean. Taxa originating in the deep sea might also have long ghost lineages and spuriously short fossil-based divergence times, especially before the Jurassic.

Architecture and Fossil Preservation

A key component of stratigraphic paleobiology is the recognition that patterns of fossil occurrences track the presence of particular sedimentary facies through stratigraphic successions. Organisms live in their preferred environments, and fossil occurrences in the rock record partly reflect these preferences (Patzkowsky and Holland Reference Patzkowsky and Holland2012).

However, only a limited portion of sedimentary environments and their inhabitants have the potential for preservation in the rock record. The taphonomic factors that contribute to the preservation and concentration of fossils are an important control on fossil occurrences. Some environments are unlikely to have a fossil record because they are unattractive habitats, such as sandy deserts for most plants and large vertebrates. Others are unlikely to have a fossil record because environmental conditions are not conducive to preservation of organic remains, such as acidic tropical forest soils (Behrensmeyer et al. Reference Behrensmeyer, Kidwell, Gastaldo, Erwin and Wing2000). The physical setting of sedimentary environments is also an important factor in their long-term preservation potential and greatly constrains the types of environments and biotas known in the rock record. For example, the record of high-elevation areas, caves, and deep-sea environments becomes progressively poorer into deep time.

Facies may lack a fossil record because they have low potential for fossil preservation. Preservation potential varies by environment and by organism, such that some environments have the potential to preserve a wide range of organisms (e.g., offshore transition, abandoned channels, stratified lakes), whereas others have low potential to preserve organic remains (e.g., rocky alluvial fans; Behrensmeyer and Hook Reference Behrensmeyer, Hook, Behrensmeyer, Damuth, DiMichele, Potts, Sues and Wing1992; Behrensmeyer et al. Reference Behrensmeyer, Kidwell, Gastaldo, Erwin and Wing2000; Holland and Loughney Reference Holland and Loughney2021). Even species that are abundant in a particular habitat may be underrepresented or absent in the fossil record because their remains are physically or chemically prone to alteration or decomposition in that environment (Behrensmeyer et al. Reference Behrensmeyer, Kidwell, Gastaldo, Erwin and Wing2000).

In marine settings, many facies are likely to preserve macroinvertebrates. Their hard parts and the fact that their remains contribute to the rock that preserves them result in fossil-rich accumulations that are tied to relative sea-level fluctuations (Kidwell Reference Kidwell, Allison and Briggs1991b). Marine communities also modify their environments through feedbacks that can be related to stratigraphic drivers. Models of hard-part accumulation in relation to sedimentation show how fossil concentrations form in response to stratigraphic controls and also through the biological response of fauna and the rates of shell production (Kidwell Reference Kidwell1985, Reference Kidwell1986; Tomašových et al. Reference Tomašových, Fürsich and Olszewski2006a, Reference Tomašových, Dominici, Nawrot, Zuschin, Nawrot, Dominici, Tomašových and Zuschin2023a). When sedimentation rate decreases relative to the input of skeletal material, burrowing infauna decrease and encrusting epifauna increase in abundance with the changing character of the microenvironment (Kidwell and Jablonski Reference Kidwell, Jablonski, Tevesz and McCall1983). When sedimentation rate increases relative to the input of skeletal material, the opposite trend occurs (Kidwell Reference Kidwell1985, Reference Kidwell1986). These feedbacks alter the environment by making it more attractive to encrusters with increasing skeletal input and more attractive to burrowers with decreasing skeletal input. In this way, the fauna of a given bathymetric zone can vary with the rates of terrigenous sediment input and shell production. Sediment mixing can also be an important control on the pattern of fossil occurrences at scales of centimeters to meters, depending on the depth of bioturbation (Tomašových et al. Reference Tomašových, Kidwell, Alexander and Kaufman2019).

In nonmarine settings, where fossil occurrences and concentrations are rare compared with marine systems, patterns of fossil occurrences tied to particular facies have long been noted (Behrensmeyer and Hook Reference Behrensmeyer, Hook, Behrensmeyer, Damuth, DiMichele, Potts, Sues and Wing1992; Holland and Loughney Reference Holland and Loughney2021). The occurrence of such facies can be stratigraphically controlled, although the larger geologic setting of different successions can yield differing patterns of fossil occurrences within facies (Behrensmeyer et al. Reference Behrensmeyer, Willis and Quade1995). For example, the preservation of fossil plants is strongly controlled by moisture, as well-drained sediments promote rapid aerobic respiration of organic matter, and plant fossils are typically known from a narrow range of facies representing lacustrine, swampy, or low-lying and wet environments (Behrensmeyer and Hook Reference Behrensmeyer, Hook, Behrensmeyer, Damuth, DiMichele, Potts, Sues and Wing1992; Gastaldo and Demko Reference Gastaldo, Demko, Allison and Bottjer2011). Vertebrate fossil preservation can be variable in facies that typically form in fluvial environments, and whether fossils occur in channel or floodplain facies may vary among sedimentary basins. Wet floodplain and abandoned-channel deposits are common settings for vertebrate preservation and have the potential to preserve rare taxa compared with better-drained floodplain deposits (Behrensmeyer Reference Behrensmeyer, Ethridge, Flores and Harvey1987; Therrien and Fastovsky Reference Therrien and Fastovsky2000; Loughney et al. Reference Loughney, Fastovsky and Parker2011). In some successions, channel deposits are the dominant setting for vertebrate assemblages, such as in the extensive Miocene Siwalik sequence of Pakistan (Badgley and Behrensmeyer Reference Badgley and Behrensmeyer1980; Behrensmeyer Reference Behrensmeyer, Ethridge, Flores and Harvey1987). In other nonmarine successions, bone concentrations are found predominantly in floodplain paleosols, such as in the Eocene Willwood Formation of Wyoming, USA (Bown and Kraus Reference Bown and Kraus1981), although it has been questioned whether this material accumulated on a soil surface or in sediment that was subsequently pedogenically modified (Rogers and Brady Reference Rogers and Brady2010).

The environment of deposition is important to fossil preservation, but subsequent landscape processes can determine the long-term preservation or architectural context of fossil occurrences. For example, plant fossil preservation is especially subject to subsequent landscape evolution and the position of the water table after deposition and burial (Gastaldo and Demko Reference Gastaldo, Demko, Allison and Bottjer2011). Plants that might have been initially buried in moist sediments that favor fossil preservation can be subsequently oxidized by a later fall in the position of the water table, which can prevent fossilization or result in lower-quality moldic preservation. The highest preservational potential occurs when accommodation and sedimentation increase following an episode of low or even negative rates of accommodation, and fossils are predicted to occur in such intervals, called expansion zones (sensu Martinsen et al. Reference Martinsen, Ryseth, Helland-Hansen, Flesche, Torkildsen and Idil1999). In the Cretaceous Two Medicine and Judith River Formations of Montana, fossiliferous concentrations occur along discontinuity surfaces that represent the reworking of previously deposited material from surrounding coastal plain facies (Rogers and Kidwell Reference Rogers and Kidwell2000). Interestingly, the probabilities of occurrence of invertebrate, plant, and vertebrate fossils do not differ among floodplain and channel facies of the Judith River Formation, although plant fossils are overall more abundant in deposits of the HAST than the LAST (Regan et al. Reference Regan, Rogers and Holland2022).

Where the constraints on fossil preservation are systematic, they can greatly affect not only the reconstructions of ancient environments and biota but also patterns of turnover and extinction. The fossil record of the Karoo Basin in South Africa is doubly affected by sparse facies bearing macroflora and their occurrence within condensed successions rife with lacunae (Gastaldo and Bamford Reference Gastaldo and Bamford2023). As a result, hiatuses and the greatly differing contexts of vertebrate and plant preservation create patterns of fossil occurrences that are stratigraphically governed yet misinterpreted as reflecting the history of the end-Permian extinction.

Large-Scale Architecture

The hierarchical nature of sequence-stratigraphic processes means that many of the concepts of stratigraphic paleobiology apply at a range of scales, from the outcrop to the continent. The major drivers of stratigraphic architecture (i.e., tectonic subsidence, climate and its control on sediment supply, eustatic sea level) can be detected at the basin scale but operate over wider areas and very long temporal ranges. As online resources have expanded, paleobiologists have been able to use extensive databases to evaluate patterns of stratigraphic architecture and fossil occurrences over greater distances and longer spans of time, as well as the long-term covariation of the rock and fossil records.

Sea-level cycles of varying scales and temporal durations are a well-known feature of the Phanerozoic rock record from which depositional sequences at all scales have been described (Sloss Reference Sloss1963). Sloss (Reference Sloss1976) considered the large-scale controls on the deposition and preservation of cratonic deposits through time as linked to a global mechanism unrelated to eustasy, presaging the concept of dynamic topography by mantle convection (Lovell Reference Lovell2010).

Raup (Reference Raup1976) was among the first to consider patterns of marine diversity in the fossil record as a function of the preservation of the Phanerozoic rock record. Although he considered sampling to account for the correlation of sedimentary rock area and volume and fossil diversity, subsequent workers have investigated the extent to which the fates of the rock and fossil records are linked and the large-scale processes that control them. For example, short-term patterns in diversity appear to be highly shaped by rock availability (Peters and Foote Reference Peters and Foote2001). Macrostratigraphy has greatly enabled the analysis of these questions by quantifying temporal changes in gap-bound rock packages, allowing them to be analyzed similarly to species’ ranges (Foote Reference Foote, Erwin and Wing2000; Peters Reference Peters2005, Reference Peters2006; Peters et al. Reference Peters, Quinn, Husson and Gaines2022). The Macrostrat database (Peters et al. Reference Peters, Husson and Czaplewski2018) has since evolved into an essential tool for quantitatively analyzing changes in the rock record over regional to continental scales, and its integration with the Paleobiology Database (PBDB; www.paleobiodb.org) enables fossil occurrence and diversity patterns to be analyzed within this framework.

Peters (Reference Peters2005) used these methods to evaluate the origination and extinction of marine sedimentary rock packages and marine animal genera over the Phanerozoic, showing them to be highly correlated. This strong correlation suggests that, while the fossil record is inherently linked to the rock record and the preservation of habitable facies, both records are shaped by similar drivers over large temporal and spatial scales; this is known as the common cause hypothesis (Peters Reference Peters2005; Peters et al. Reference Peters, Quinn, Husson and Gaines2022). Moreover, extinction rates of marine packages and genera in the Phanerozoic of North America are strongly correlated, whereas origination rates are not (Heim and Peters Reference Heim and Peters2011). This finding indicates that large-scale relationships between stratigraphic architecture and genus first and last occurrences differ from those observed and predicted at smaller scales (Holland Reference Holland1995; Peters and Heim Reference Peters and Heim2011), and common drivers may shape package truncations and biotic extinctions. Similarly, diversity and extinction of planktonic foraminifera in the Atlantic Ocean are highly correlated with the number of sedimentary packages, whereas origination is not (Peters et al. Reference Peters, Kelly and Fraass2013). This discrepancy suggests that a biological signal is being preserved, because a dominant stratigraphic control would be expected to affect origination, extinction, and diversity equally.

The correlation between package truncation and extinction rates highlights the importance of evaluating taxon last occurrences at a variety of scales and within a stratigraphic-paleobiological context (Heim and Peters Reference Heim and Peters2011; Holland and Patzkowsky Reference Holland and Patzkowsky2015). At the global scale, this approach has also been used to examine diversity in relation to supercycles of supercontinent formation and breakup (Zaffos et al. Reference Zaffos, Finnegan and Peters2017). The diversity of marine genera is correlated with increasing continental fragmentation through the Phanerozoic, particularly following the breakup of Rodinia and Pangaea. Intervals of low diversity correspond with coalescing continental movements, a relationship that may have contributed to the persistently low marine diversity through the early Mesozoic following the end-Permian extinction.

A macrostratigraphic approach also illuminates the occurrence of stromatolites (Peters et al. Reference Peters, Husson and Wilcots2017). Normalizing for rock volume shows stromatolites increasing in the Paleoproterozoic (ca. 2500 Ma), dominating through most of the Proterozoic (2500–800 Ma, and declining through the late Proterozoic (700–541 Ma). Moreover, despite the conventional wisdom that stromatolites increase in abundance following mass extinctions, a macrostratigraphic approach shows that stromatolite occurrence is not well correlated with times following mass extinctions. Lithologically, stromatolite occurrence is most strongly associated with dolomite, which may poorly preserve other fossils. Macrostratigraphy has also demonstrated the strong tie between the Cambrian explosion and the profound changes in geochemical cycling as the Great Unconformity was transgressed and buried (Peters and Gaines Reference Peters and Gaines2012; Peters et al. Reference Peters, Quinn, Husson and Gaines2022).

In the nonmarine record, macrostratigraphy can be used to understand diversity and landscape evolution in relation to tectonics and climate. For example, Loughney et al. (Reference Loughney, Badgley, Bahadori, Holt and Rasbury2021) evaluated changes in Neogene mammalian species richness in relation to sedimentation history during the tectonic evolution of the Basin and Range province in western North America. The history of sediment accumulation tracks landscape evolution over the formation of the Basin and Range and partly determines the timing of deposition of fossil-bearing sedimentary units. In contrast, spikes in species richness occurred during intervals of global warming, indicating that species’ responses to climate may be enhanced in topographically complex regions and are decoupled from preservation.

Assembly of Regional Ecosystems

At any time in Earth's history, global diversity is an amalgam of biogeographic provinces, latitudinal and elevational diversity gradients, and biodiversity hotspots and coolspots. Changes in global diversity over time must therefore reflect changes in this biogeographic structure, and the processes that control the spatial and temporal structure of diversity in local and regional ecosystems, such as origination, extinction, dispersal, and extirpation, are a major focus of research (Jablonski et al. Reference Jablonski, Roy and Valentine2006; Patzkowsky Reference Patzkowsky2017; Close et al. Reference Close, Benson, Saupe, Clapham and Butler2020; Saupe Reference Saupe2023).

The repeated deposition of sedimentary packages containing depth or elevation gradients provides a time–environment framework at subsequence scale necessary to study the main patterns, processes, and rates of change in regional diversity in deep time. With a time–environment framework, it is possible to evaluate the completeness of the environmental gradient, so that the absence of a taxon is interpreted as either a true absence or as missing due to a failure to preserve the preferred gradient position. Taxon occurrences outside the depositional basin indicate whether last occurrences represent true extinctions or extirpations, or likewise, whether first occurrences represent originations or immigrations. If the time–environment framework spans an environmental shift, the effect of the environmental change on the regional ecosystem can be evaluated, because the ecological gradient can be quantified before, during, and after the environmental shift (Patzkowsky and Holland Reference Patzkowsky and Holland2007; Dominici and Danise Reference Dominici, Danise, Nawrot, Dominici, Tomašových and Zuschin2022). In addition, it is possible to track individual taxa through time to test for niche conservatism and niche evolution (Holland and Zaffos Reference Holland and Zaffos2011).

Pulsed Taxonomic Turnover

Taxonomic turnover in global and local marine ecosystems has not been continuous through time, but instead was pulsed with long intervals (several million years) of persistence in marine communities arrayed along a depth gradient, interrupted by short intervals (a few hundred thousand years) of rapid, ecosystem-wide turnover (Boucot Reference Boucot1978, Reference Boucot1996; Brett and Baird Reference Brett, Baird, Erwin and Anstey1995; Brett et al. Reference Brett, Ivany, Bartholomew, DeSantis, Baird and Konigshof2009). Notably, these patterns were thought to be widespread throughout the fossil record, and turnover pulses were interpreted as pulses in both origination and extinction.

More recent work on turnover in regional marine ecosystems using a time–environment framework found that these patterns of turnover were much more nuanced and that they have implications for the processes that structured regional ecosystems in the past (Patzkowsky and Holland Reference Patzkowsky and Holland1993, Reference Patzkowsky and Holland1997, Reference Patzkowsky and Holland1999; Holland Reference Holland and Brett1997). In some cases, turnover was dominated by pulses in extinction or origination, but not by pulses in both. Based on knowledge of stratigraphic occurrence of taxa outside the basin, some extinctions were shown to be extirpations and some originations to be immigrations. Turnover was usually associated with shifts in regional environmental conditions, but extinctions and extirpations and originations and immigrations sometimes occurred over more than one depositional sequence. These patterns indicate that the duration of the environmental change was also spread out over multiple sequences and that taxa had individual responses to environmental change, rather than all taxa responding at once.

Changes in Gradient Structure

The shift in thought to view fossil assemblages as positions along environmental gradients (e.g., Fig. 1) rather than as parts of individual communities led to new ways to quantify ecological gradients and how they change through time (Cisne and Rabe Reference Cisne and Rabe1978; Patzkowsky Reference Patzkowsky1995; Ricklefs Reference Ricklefs2008). The predominate theme in these studies is that ecological gradients persist with slow turnover for hundreds of thousands to millions of years, often in the face of significant environmental change, although the rate of turnover can vary between regions and time periods and among taxa (Brett and Baird Reference Brett, Baird, Erwin and Anstey1995; Patzkowsky and Holland Reference Patzkowsky and Holland1997, Reference Patzkowsky and Holland1999; Olszewski and Patzkowsky Reference Olszewski and Patzkowsky2001; Ivany et al. Reference Ivany, Brett, Wall, Wall and Handley2009; Zuschin et al. Reference Zuschin, Harhauser, Hengst, Mandic and Roetzel2014). These prolonged intervals of slower turnover suggest that taxa track their preferred environments as base level rises and falls, and it raises questions about the kind and intensity of environmental change that is required for the gradient structure to break down.

This framework can be used to address many more questions about patterns and rates of ecosystem persistence. Intervals of pulsed turnover have also been shown to vary along an onshore-to-offshore gradient, with offshore communities showing higher levels of extinction than onshore communities (Jarochowska et al. Reference Jarochowska, Munnecke, Frisch, Ray and Castagner2016; Danise and Holland Reference Danise and Holland2017). This pattern is contrary to other claims of higher extinction in onshore Paleozoic and Mesozoic communities (Sepkoski Reference Sepkoski1987; Kiessling and Aberhan Reference Kiessling and Aberhan2007), pointing to a disconnect in the direction of onshore–offshore gradients in extinction rate between short-term pulsed extinction compared with long-term background extinction. Similarly, Tomašových et al. (Reference Tomašových, Dominici, Zuschin and Merle2014) used a time–environment framework to compare marine depth gradients from two time intervals, the Eocene and Plio–Pleistocene, to show that the onshore–offshore gradient in taxonomic turnover (higher onshore, lower offshore) is not manifested at shorter timescales (less than 5 Myr), but it does show up over longer timescales (ca. 50 Myr). Bonelli and Patzkowsky (Reference Bonelli and Patzkowsky2008, Reference Bonelli and Patzkowsky2011) found little extinction associated with the onset of the late Paleozoic Ice Age (LPIA) in the Illinois Basin, despite a major change in gradient structure, with a well-defined gradient on the pre-LPIA carbonate platform followed by a more poorly defined gradient dominated by eurytopes on a more steeply dipping carbonate ramp after the onset of the LPIA.

Niche Conservatism

The persistence of ecological gradients implies that marine benthic taxa conserve their niches over millions of years. By using multivariate methods to identify biotic gradients, it is possible to characterize taxon response curves (preferred environment, peak abundance, environmental breadth) and to compare these parameters through time to test for niche conservatism. Studies that have done this (Holland and Patzkowsky Reference Holland and Patzkowsky2004; Holland and Zaffos Reference Holland and Zaffos2011; Tomašových et al. Reference Tomašových, Dominici, Zuschin and Merle2014; Dominici and Danise Reference Dominici, Danise, Nawrot, Dominici, Tomašových and Zuschin2022) have found that niches are conserved for several million years over multiple depositional sequences and even through mass extinction events (Brisson et al. Reference Brisson, Pier, Beard, Fernandes and Bush2023). These studies are consistent with others, suggesting that niche conservatism is common throughout the Phanerozoic among marine invertebrate groups (Hopkins et al. Reference Hopkins, Simpson and Kiessling2014; Saupe et al. Reference Saupe, Hendricks, Portell, Dowsett, Haywood, Hunter and Lieberman2014; Antell et al Reference Antell, Fenton, Valdes and Saupe2021). High-resolution stratigraphic correlation and environmental niche modeling in the Upper Ordovician of Ohio, USA, indicates that brief intervals of niche evolution have occurred during biotic invasion followed by niche conservatism (Forsythe and Stigall Reference Forsythe and Stigall2023). These studies indicate the great potential for resolving ecological and evolutionary processes at the sequence and subsequence scale.

Community Saturation and Diversity Partitioning

Although recognizing continuous ecological change along environmental gradients is a hallmark of stratigraphic paleobiology and has led to an improved understanding of processes, lithofacies can be used to identify depositional environments (e.g., shallow subtidal, deep subtidal) along the environmental gradient. Defining discrete environments can be used to advantage to ask specific questions about limits to diversity within environments and the partitioning of diversity within and between environments (Bambach Reference Bambach1977; Sepkoski Reference Sepkoski1988; Patzkowsky and Holland Reference Patzkowsky and Holland2007; Holland Reference Holland2010; Dominici and Danise Reference Dominici, Danise, Nawrot, Dominici, Tomašových and Zuschin2022). This approach was used to investigate whether richness in Middle and Late Ordovician marine communities was saturated with species by tracking shallow and deep subtidal environments across multiple sequences. A positive relationship between local and regional diversity over 10 depositional sequences suggested that Late Ordovician marine communities were not saturated with species (Patzkowsky and Holland Reference Patzkowsky and Holland2003). Diversity partitioning among shallow and deep subtidal communities showed that most of the regional diversity increase caused by the Richmondian invasion was accommodated by increases in beta diversity, rather than by increases in local diversity (Patzkowsky and Holland Reference Patzkowsky and Holland2007), indicating plenty of ecospace for new taxa to inhabit the environmental gradient. In a further comparison of faunal provinces across Laurentia, it was shown that during a major biotic invasion, niche parameters of invading taxa were highly conserved between the donor area and the recipient area, suggesting the local communities were open to invasion and the invading taxa maintained their environmental niches and could easily enter the communities they invaded (Patzkowsky and Holland Reference Patzkowsky and Holland2016).

The water-depth gradient in marine environments has been explored extensively, but this is a complex gradient with many factors that covary with depth, and there are other environmental axes that could be specifically explored, such as lithology, grain size, latitude, elevation, substrate consistency, salinity, oxygenation, and nutrients, to name a few (Holland and Patzkowsky Reference Holland and Patzkowsky2007; Belanger and Garcia Reference Belanger and Garcia2014; Danise et al. Reference Danise, Twitchett and Little2015; Zuschin et al. Reference Zuschin, Nawrot, Harzhauser, Mandic and Tomašových2017; Brown and Larina Reference Brown and Larina2019; Dominici and Danise Reference Dominici, Danise, Nawrot, Dominici, Tomašových and Zuschin2022; Bryant and Belanger Reference Bryant and Belanger2023). There is a very rich array of questions to explore on the origin and assembly of regional ecosystems and how they change over time. These questions can be best addressed by using a time–environment framework built on depositional sequences, by quantifying ecological gradients, and by collecting and analyzing environmental information to characterize the multiple sources of variation in community structure.

Future Prospects

Climate Change

Numerous studies document the effects of modern and ancient climate change on biotic community composition, evolution, and mass extinction. Climate is one important control on species distributions and stratigraphic architecture, although most studies investigating these relationships do not explicitly consider them from the perspective of stratigraphic paleobiology. Even so, there is great potential in studying climate-related patterns of biotic change using the principles of stratigraphic paleobiology.

Studies investigating the relationship between climate change and species diversity patterns yield variable results, underscoring that biodiversification events are commonly shaped by a “perfect storm” of contributing factors (Jablonski and Edie Reference Jablonski and Edie2023). For example, although the magnitude and rate of Phanerozoic temperature changes are positively correlated with marine animal extinction rates (Song et al. Reference Song, Kemp, Tian, Chu, Song and Dai2021), the extinction sensitivity of Phanerozoic marine invertebrates lacks a strong relationship with climate (Reddin et al. Reference Reddin, Kocsis and Kiessling2019). Global and regional analyses of mammalian diversity (Alroy et al. Reference Alroy, Koch, Zachos, Erwin and Wing2000; Barnosky and Carrasco Reference Barnosky and Carrasco2002) likewise show no strong relationship to climate change. Even so, mammalian faunal composition does track global climate trends (Figueirido et al. Reference Figueirido, Janis, Perez-Claros, de Renzi and Palmqvist2012). It is perhaps not surprising that something as extensive as the climate system should have unexpected or differing effects on biological systems, and there likely are aspects of the climate–biosphere connection that we either do not understand or cannot measure (Reddin et al. Reference Reddin, Kocsis and Kiessling2019). Moreover, recent climate changes demonstrate that local climate patterns need not track global patterns: for example, although global temperatures are currently increasing, the temperature of some regions is in net decline (Hansen et al. Reference Hansen, Ruedy, Sato and Lo2010). As such, analyzing regional biotic changes in light of regional rather than global climate proxies is essential (Raffi et al. Reference Raffi, Stanley and Marasti1985; Garilli Reference Garilli2011; Dominici et al. Reference Dominici, Benvenuti, Garilli, Uchman, Pollina and David2020; Borghi and Garilli Reference Borghi and Garilli2022). Similarly, marine-derived paleoclimate proxies do not necessarily reflect nonmarine climate, which may explain some discrepancies among climate and nonmarine diversity patterns (Alroy et al. Reference Alroy, Koch, Zachos, Erwin and Wing2000).

The availability of global climate-proxy datasets allows large-scale analyses with diversity data from the PBDB. For most analyses involving Cenozoic climate, Zachos et al. (Reference Zachos, Pagani, Sloan, Thomas and Billups2001) has been the favored source, with Westerhold et al. (Reference Westerhold, Marwan, Drury, Liebrand, Agnini, Anagnostou and Barnet2020) ascending with updated records; Veizer and Prokoph (Reference Veizer and Prokoph2015) and Scotese et al. (Reference Scotese, Song, Mills and van der Meer2021) provide reconstructions spanning the Phanerozoic. Although these are useful and powerful resources for investigating change through time, the patterns cannot necessarily be taken at face value. A stratigraphic paleobiological perspective that considers the spatial distribution of sedimentary basins and their intersections with biomes and provinces is crucial for understanding how the record is assembled and how it can be interpreted (Holland Reference Holland2016). Similarly, studies should reflect the implications of the large-scale selectivity of the fossil record, such as the dominance of low-elevation habitats in the nonmarine fossil record and the preferential preservation of particular shallow-marine habitats (Nyberg and Howell Reference Nyberg and Howell2015, Reference Nyberg and Howell2016; Holland et al. Reference Holland, Loughney and Cone2023).

For large-scale analyses of climate and biodiversity, macrostratigraphy provides an opportunity to investigate the coupled effects of changing climate on diversity, sea level, and rock availability (Peters et al. Reference Peters, Quinn, Husson and Gaines2022). For example, Late Ordovician extinction and sedimentary packages appear to be jointly driven by climate and eustasy (Finnegan et al. Reference Finnegan, Heim, Peters and Fischer2012). This approach may also be useful in establishing other biological correlates of climate change, such as ecological and functional traits.

Geochemistry and Stratigraphic Paleobiology

Global geochemical proxies of environmental conditions are often matched against fossil occurrences to infer environmental causes of extinction (e.g., Saltzman et al. Reference Saltzman, Edwards, Adrain and Westrop2015), often with little consideration of how stratigraphic architecture affects fossil occurrences. Even so, geochemical proxies are also affected by variation in sedimentation rates, facies variability, and hiatuses, just as are fossil occurrences (Railsback et al. Reference Railsback, Holland, Hunter, Jordan, Diaz and Crowe2003; Trampush and Hajek Reference Trampush and Hajek2017; Hay et al. Reference Hay, Creveling, Hagen, Maloof and Huybers2019; Pulsipher et al. Reference Pulsipher, Schiffbauer, Jeffrey, Huntley, Fike and Shelton2021; Prave et al. Reference Prave, Kirsimäe, Lepland, Fallick, Kreitsmann, Deines and Romashkin2022), so they cannot be interpreted as literal indicators of the timing and intensity of environmental change. Subaerial exposure and diagenesis are additional sources of variability in proxy records. Studies of modern carbonate environments indicate that the variability among facies in δ13C values of carbonate sediment grains equals that seen in the deep-time record (Geyman and Maloof Reference Geyman and Maloof2021). Simulated stratigraphic sections using these modern values show abrupt jumps in δ13C values across facies shifts at flooding surfaces. In the deep-time record, δ13C shifts at flooding surfaces would be expected to co-occur with clusters of first or last occurrences, and this could easily be misinterpreted as a cause-and-effect relationship rather than just the expectation across abrupt shifts in facies. On the other hand, fossil occurrences and geochemical proxies may not respond similarly to environmental gradients, so that even if a specific environmental shift caused an extinction, it may not coincide stratigraphically with fossil occurrences in a way that permits an interpretation as simple cause and effect.

Integrating local environmental proxies with fossil occurrences in measured stratigraphic sections holds promise for making direct inferences about environmental causes of extinction. For example, marine anoxia likely contributed to many regional extinctions in the geologic past, and its effects on benthic organisms should be observable. However, the expansion of ancient marine anoxic zones is often marked by the appearance of black shales at a flooding surface. As a result, the coincident extinction of shallow-water taxa cannot be confirmed until those taxa can be shown to be absent in shallow-water facies above the flooding surface. Geochemical proxies of local seafloor anoxia can help identify oxygen fluctuations where facies do not change substantially (Pruss et al. Reference Pruss, Jones, Fike, Tosca and Wignall2019). Pairing geochemical proxies for local marine anoxia with fossil occurrences in a time–environment framework built on depositional sequences could reveal the direct effects of seafloor oxygen fluctuations on benthic communities. Similar studies integrating stratigraphic paleobiology and geochemistry are needed to test for specific oceanographic or climatic variables as causes of extinction.

Phylogenetics and Stratigraphic Paleobiology

Phylogenetic methods are central to many of the basic questions in evolutionary paleoecology surrounding the origin, diversification, and environmental distribution of clades. Phylogenetic trees scaled against geologic time can be used to estimate rates of speciation, extinction, and morphological evolution. Adding information about abundance and habitat preference allows an evaluation of how rates of diversification and morphological evolution are determined by environment (Ricklefs Reference Ricklefs2006; Lamsdell et al. Reference Lamsdell, Congreve, Hopkins, Krug and Patzkowsky2017). Stratigraphic paleobiological principles can support and enhance phylogenetic studies in two important ways. First, estimating the time of origin and rates of speciation and extinction in diversifying clades is based on the record of species’ first and last occurrences. Early development of phylogenetic models to make these estimates were based on simplistic assumptions of fossil recovery, specifically uniform preservation rates (Barido-Sottani et al. Reference Barido-Sottani, Pett, O'Reilly and Warnock2019). Fossil recovery is demonstrably nonuniform on both short timescales (<10 Myr), where fossil recovery is controlled by stratigraphic architecture and environmental distributions of species (Patzkowsky and Holland Reference Patzkowsky and Holland2012), and longer timescales (>10 Myr), where fossil recovery is controlled by tectonic processes that construct and destroy sedimentary basins (Holland Reference Holland2016). Models incorporating more realistic fossil recovery rates based on stratigraphic architecture are needed (e.g., Barido-Sottani et al. Reference Barido-Sottani, Pett, O'Reilly and Warnock2019). Case studies of specific clades should utilize empirical estimates of fossil preservation using databases of fossil and lithologic occurrences, such as the PBDB and Macrostrat (Wagner and Marcot Reference Wagner and Marcot2013; Congreve et al. Reference Congreve, Patzkowsky and Wagner2021).

Second, phylogenies of diversifying clades can be mapped onto the time–environment framework provided by a stratigraphic-paleobiological approach. Doing this would help explain, for example, the distribution of niche conservatism and niche expansion among subclades in maintaining environmental gradients in diversity (Sepkoski Reference Sepkoski1991; Jablonski et al. Reference Jablonski, Roy and Valentine2006; Ricklefs Reference Ricklefs2006; Tomašových et al. Reference Tomašových, Dominici, Zuschin and Merle2014). Merging phylogenetics with stratigraphic paleobiology promises to accelerate progress in understanding the ecological and evolutionary processes of clade origin and diversification and the assembly of regional biotas in an environmental context.

Large-Scale Structure of the Fossil Record

The PBDB has permitted extraordinary insight into long-term and large-scale patterns in the fossil record for more than 20 years. A principal challenge to interpreting the biological meaning of this archive is coming to terms with the temporal and geographic inhomogeneities in the data. Despite substantial efforts by an ever-growing number of paleontologists, large portions of the published fossil record have not been entered, and the PBDB remains an uneven sampling of the known fossil record (Ye and Peters Reference Ye and Peters2023). At large scales (global, Phanerozoic), the hope is that the sampling is even enough that patterns are biologically meaningful, but at smaller geographic and shorter temporal scales, uneven sampling of the known record becomes more of a concern.

Even if these effects were overcome by an exhaustive catalogue of the known fossil record, the inference of biological meaning to patterns in the PBDB must increasingly hinge on a stratigraphic-paleobiological outlook, specifically that even over large spatial (basin to globe) and temporal (tens to hundreds of millions of years) scales, the fossil record is systematically structured by how the stratigraphic record is constructed. Macrostratigraphy offers the most promise at these scales (Peters et al. Reference Peters, Quinn, Husson and Gaines2022).

One such challenge is the global temporal correlation in the available stratigraphic record. For example, Ager (Reference Ager1973) pointed out the lateral persistence of facies—that certain facies are extraordinarily widespread at particular times in Earth's history, such as Upper Cretaceous chalk, Tithonian limestone, Pennsylvanian coal, and Mississippian limestone, and many others. Although this partly reflects global climate and the configuration of continents, it must be remembered that the presence of some facies implies that others must have existed. For example, Upper Devonian black shales are widespread across North America, yet the presence of these offshore shales implies that coastal marine environments must also have existed. Owing to their subsequent erosion or inaccessibility, it may be impossible to sample them (Smith et al. Reference Smith, Gale and Monks2001), and as a result, the regional Late Devonian drop in diversity may be largely a reflection of the facies available for study. This can be addressed only by considering the nature of the available (and sampled) stratigraphic record. Similarly, the formation of sedimentary basins—and the types of basins, which differ greatly in their extent and life spans—is likely to be highly correlated through geologic time, owing to global cycles in supercontinent formation and breakup (Holland Reference Holland2016). This not only affects global and continental diversity (Holland Reference Holland2018), but it also affects what life can be preserved as fossils and where that can happen (Holland Reference Holland2016). Because of these global correlations, it can be difficult to distinguish whether patterns in the fossil record represent what has been preserved in the stratigraphic record or patterns in the history of life (Miller Reference Miller, Erwin and Wing2000). A promising area of investigation is understanding how the history of life—particularly biogeography—is driven by plate tectonics (e.g., Holland Reference Holland2018; Jablonski and Edie Reference Jablonski and Edie2023), but also how the formation and cessation of sedimentary basins selectively preserve this history of life (Holland Reference Holland2016; Holland et al. Reference Holland, Loughney and Cone2023). A stratigraphic paleobiological perspective at these wide spatial and long temporal scales offers the best chance for extracting the history of life from the fossil record.

On the Need for Stratigraphic Thinking in Paleobiology

One of our growing concerns, also raised by two reviewers, is the apparent decline of field-based research in paleobiology and the limited application of sedimentary geology in interpreting the fossil record. Similar concerns have been raised previously, notably in Alan Shaw's Butterfingered Handmaiden (Shaw Reference Shaw1971) and Reuben Ross's editorial in the first issue of Palaios (Ross Reference Ross1986). Rather than a decline per se, this may reflect the shifting balance of intellectual emphasis in an interdisciplinary field like paleobiology (Newell and Colbert Reference Newell and Colbert1948). Even so, it is concerning.

Our concern is not a criticism of other approaches to our science. Paleobiology has flourished because it has embraced diverse approaches: field-based, museum-based, numerical simulations, and database studies; indeed, we have actively participated in these approaches. We view all these as good for the discipline; they have greatly contributed to the relevancy and richness of paleobiology, and this 50th anniversary issue testifies to that. Even so, over the 30 years we have published on stratigraphic paleobiology, we have noticed that relatively few research groups have pursued active field-based studies. It is unclear why. Perhaps it reflects the difficulty and expense of fieldwork. Maybe it is because gathering data in the field can be slower and more complicated than working with specimens in a museum or records in a database. It could be linked to the downscaling or rebranding of many university geology departments. That it might be a lack of training has been a prime motivation behind the field conference in stratigraphic paleobiology that we have developed for graduate students. Regardless, we feel that paleobiology would be on a stronger footing if fieldwork and the development of new data were pursued more commonly.

More broadly, we have been concerned by the stratigraphic simplicity of much paleobiological research, regardless of whether that research is based in the field. One of the goals of stratigraphic paleobiology is to build a broad recognition that the fossil record is not simply the history of life; the fossil record is at least equally the history of sedimentation. The fossil record is not simply a time series. Patterns in the fossil record should always be evaluated for how they have been generated, from the interaction of ecology and stratigraphic architecture to the global intersection of biogeographic realms with the development of sedimentary basins. Too often, paleobiologists assume that a literal reading of the fossil record is permissible, that the probability of preservation is uniform, that hiatuses and facies changes are absent or unimportant, and that the preserved geographic distribution of species faithfully reflects their original distribution. Many studies make these statements explicitly, but our goal is not to single out individual examples for rebuke but to say that these assumptions about the fossil record are usually incorrect to an extent that substantially compromises interpretations of the fossil record. Any study that asserts them but does not demonstrate them should be viewed with this point in mind. We encourage paleobiological studies to incorporate a modern understanding of the stratigraphic record: how it is assembled and how it shapes patterns throughout the fossil record.

Acknowledgments

We appreciate the helpful and constructive comments of Associate Editor A. Tomašových, reviewers S. Dominici and R. Rogers, and an anonymous reviewer. We also thank many colleagues for their conversations about these topics: C. Badgley, A. Behrensmeyer, S. Danise, M. Deckman, M. Foote, B. Gastaldo, K. Irizarry, S. Khatri, S. Kidwell, A. Miller, P. Monarrez, S. Peters, A. Regan, R. Rogers, and R. Warnock.

Competing Interests

The authors declare no competing interests.

References

Literature Cited

Ager, D. V. 1973. The nature of the stratigraphical record. Macmillan, London.Google Scholar
Alroy, J., Koch, P. L., and Zachos, J. C.. 2000. Global climate change and North American mammalian evolution. Pp. 259288 in Erwin, D. H. and Wing, S. L., eds. Deep time: Paleobiology's perspective. Paleontological Society, Lawrence, Kans.Google Scholar
Amorosi, A., Centineo, M. C., Colalongo, M. L., Paasini, G., Sarti, G., and Vaiani, S. C.. 2003. Facies architecture and latest Pleistocene–Holocene depositional history of the Po Delta (Comacchio area), Italy. Journal of Geology 111:3956.CrossRefGoogle Scholar
Amorosi, A., Colalongo, M. L., Fiorini, F., Fusco, F., Pasini, G., Vaiani, S. C., and Sarti, G.. 2004. Palaeogeographic and palaeoclimatic evolution of the Po Plain from 150-ky core records. Global and Planetary Change 40:5578.CrossRefGoogle Scholar
Amorosi, A., Rossi, V., Scarponi, D., Vaiani, S. C., and Ghosh, A.. 2014. Biosedimentary record of postglacial coastal dynamics: high-resolution sequence stratigraphy from the northern Tuscan coast (Italy). Boreas 43:939954.CrossRefGoogle Scholar
Antell, G. S., Fenton, I. S., Valdes, P. J., and Saupe, E. E.. 2021. Thermal niches of planktonic foraminifera are static throughout glacial–interglacial climate change. Proceedings of the National Academy of Sciences USA 118:e2017105118.CrossRefGoogle ScholarPubMed
Atkinson, J. W., Little, C. T. S., and Dunhill, A. M.. 2023. Long duration of benthic ecological recovery from the early Toarcian (Early Jurassic) mass extinction event in the Cleveland Basin, UK. Journal of the Geological Society 180:jgs2022–126.CrossRefGoogle Scholar
Austin, M. P. 1987. Models for the analysis of species’ response to environmental gradients. Vegetatio 69:3545.CrossRefGoogle Scholar
Badgley, C., and Behrensmeyer, A. K.. 1980. Paleoecology of Middle Siwalik sediments and faunas, northern Pakistan. Palaeogeography, Palaeoclimatology, Palaeoecology 30:133155.CrossRefGoogle Scholar
Bambach, R. K. 1977. Species richness in marine benthic habitats through the Phanerozoic. Paleobiology 3:152167.CrossRefGoogle Scholar
Barido-Sottani, J., Pett, W., O'Reilly, J. E., and Warnock, R. C. M.. 2019. FossilSim: an R package for simulating fossil occurrence data under mechanistic models of preservation and recovery. Methods in Ecology and Evolution 10:835840.CrossRefGoogle Scholar
Barnosky, A. D., and Carrasco, M. A.. 2002. Effects of Oligo–Miocene global climate changes on mammalian species richness in the northwestern quarter of the USA. Evolutionary Ecology Research 4:811841.Google Scholar
Behrensmeyer, A. K. 1987. Miocene fluvial facies and vertebrate taphonomy in northern Pakistan. In Ethridge, F. G., Flores, R. M., and Harvey, M. D., eds. Recent developments in fluvial sedimentology. Society of Economic Paleontologists and Mineralogists Special Publication 39:169176.CrossRefGoogle Scholar
Behrensmeyer, A. K., and Hook, R. W.. 1992. Paleoenvironmental contexts and taphonomic modes. Pp. 15136 in Behrensmeyer, A. K., Damuth, J. D., DiMichele, W. A., Potts, R., Sues, H.D., and Wing, S. L., eds. Terrestrial ecosystems through time. University of Chicago Press, Chicago.Google Scholar
Behrensmeyer, A. K., Willis, B. J., and Quade, J.. 1995. Floodplains and paleosols of Pakistan Neogene and Wyoming Paleogene deposits: a comparative study. Palaeogeography, Palaeoclimatology, Palaeoecology 115(1–4):3760.CrossRefGoogle Scholar
Behrensmeyer, A. K., Kidwell, S. M., and Gastaldo, R. A.. 2000. Taphonomy and paleobiology. Pp. 103147 in Erwin, D. H. and Wing, S. L., eds. Deep time: Paleobiology's perspective. Paleontological Society, Lawrence, Kansas.Google Scholar
Belanger, C. L., and Bapst, D. W.. 2023. Simulating our ability to accurately detect abrupt changes in assemblage-based paleoenvironmental proxies. Palaeontologia Electronica 26:a24.Google Scholar
Belanger, C. L., and Garcia, M. V.. 2014. Differential drivers of benthic foraminiferal and molluscan community composition from a multivariate record of early Miocene environmental change. Paleobiology 40:398416.CrossRefGoogle Scholar
Bonelli, J. R., and Patzkowsky, M. E.. 2008. How are global patterns of faunal turnover expressed at regional scales? Evidence from the Upper Mississippian (Chesterian Series), Illinois Basin, USA. Palaios 23:760772.CrossRefGoogle Scholar
Bonelli, J. R., and Patzkowsky, M. E.. 2011. Taxonomic and ecologic persistence across the onset of the Late Paleozoic ice age: evidence from the Upper Mississippian (Chesterian Series), Illinois Basin, United States. Palaios 26:517.CrossRefGoogle Scholar
Borcard, D., Gillet, F., and Legendre, P.. 2018. Numerical ecology with R. Springer, New York.CrossRefGoogle Scholar
Borghi, E., and Garilli, V.. 2022. Climate-driven diversity changes of Mediterranean echinoids over the last 6 Ma. Acta Palaeontologica Polonica 67:781805.CrossRefGoogle Scholar
Boucot, A. J. 1978. Community evolution and rates of cladogenesis. Evolutionary Biology 11:545655.Google Scholar
Boucot, A. J. 1996. Epilogue. Palaeogeography, Palaeoclimatology, Palaeoecology 127:339359.CrossRefGoogle Scholar
Bown, T. M., and Kraus, M. J.. 1981. Vertebrate fossil-bearing paleosol units (Willwood Formation, Lower Eocene, Northwest Wyoming, U.S.A.): implications for taphonomy, biostratigraphy, and assemblage analysis. Palaeogeography, Palaeoclimatology, Palaeoecology 34:3156.CrossRefGoogle Scholar
Brett, C. E., and Baird, G. C.. 1995. Coordinated stasis and evolutionary ecology of Silurian to Middle Devonian faunas in the Appalachian Basin. Pp. 285315 in Erwin, D. H. and Anstey, R. L., eds. New approaches to speciation in the fossil record. Columbia University Press, New York.Google Scholar
Brett, C. E., Ivany, L. C., Bartholomew, A. J., DeSantis, M. K., and Baird, G. C.. 2009. Devonian ecological–evolutionary subunits in the Appalachian Basin: a revision and a test of persistence and discreteness. In Konigshof, P., ed. Devonian change: case studies in palaeogeography and palaeoecology. Geological Society Special Publications 314:736.Google Scholar
Brisson, S. K., Pier, J. Q., Beard, J. A., Fernandes, A. M., and Bush, A. M.. 2023. Niche conservatism and ecological change during the Late Devonian mass extinction. Proceedings of the Royal Society of London B 290:20222524.Google ScholarPubMed
Brown, G. M., and Larina, E.. 2019. Environmental controls on shallow subtidal molluscan death assemblages on San Salvador Island, The Bahamas. Palaeogeography, Palaeoclimatology, Palaeoecology 527:1424.CrossRefGoogle Scholar
Bryant, R., and Belanger, C. L.. 2023. Spatial heterogeneity in benthic foraminiferal assemblages tracks regional impacts of paleoenvironmental change across Cretaceous OAE2. Paleobiology 49:431453.CrossRefGoogle Scholar
Catuneanu, O. 2006. Principles of sequence stratigraphy. Elsevier, New York.Google Scholar
Catuneanu, O., and Zecchin, M.. 2020. Parasequences: allostratigraphic misfit in sequence stratigraphy. Earth-Science Reviews 208:103289.CrossRefGoogle Scholar
Catuneanu, O., Abreu, V., Bhattacharya, J. P., Blum, M. D., Dalrymple, R. W., Eriksson, P. G., Fielding, C. R., et al. 2009. Towards the standardization of sequence stratigraphy. Earth-Science Reviews 92:133.CrossRefGoogle Scholar
Christie, M., Holland, S. M., and Bush, A. M.. 2013. Contrasting the ecological and taxonomic consequences of extinction. Paleobiology 39:538559.CrossRefGoogle Scholar
Cisne, J. L., and Rabe, B. D.. 1978. Coenocorrelation: gradient analysis of fossil communities and its applications in stratigraphy. Lethaia 11:341364.CrossRefGoogle Scholar
Clement, A. M., and Tackett, L. S.. 2021. Facies stacking and distribution in the Gabbs Formation (Late Triassic, west-central Nevada, U.S.A.): an environmental baseline to the end-Triassic carbonate crisis. Sedimentary Geology 425:106021.CrossRefGoogle Scholar
Close, R. A., Benson, R. B. J., Saupe, E. E., Clapham, M. E., and Butler, R. J.. 2020. The spatial structure of Phanerozoic marine animal diversity. Science 368:420424.CrossRefGoogle ScholarPubMed
Congreve, C. R., Patzkowsky, M. E., and Wagner, P. J.. 2021. An early burst in brachiopod evolution corresponding with significant climatic shifts during the Great Ordovician Biodiversification Event. Proceedings of the Royal Society of London B 288:20211450.Google ScholarPubMed
Danise, S., and Holland, S. M.. 2017. Faunal response to sea-level and climate change in a short-lived seaway: Jurassic of the Western Interior, USA. Palaeontology 60:213232.CrossRefGoogle Scholar
Danise, S., Twitchett, R. J., and Little, C. T. S.. 2015. Environmental controls on Jurassic marine ecosystems during global warming. Geology 43:263266.CrossRefGoogle Scholar
Danise, S., Clémence, M.-E., Price, G. D., Murphy, D. P., Gómez, J. J., and Twitchett, R. J.. 2019. Stratigraphic and environmental control on marine benthic community change through the early Toarcian extinction event (Iberian Range, Spain). Palaeogeography, Palaeoclimatology, Palaeoecology 524:183200.CrossRefGoogle Scholar
Dominici, S., and Danise, S.. 2022. Mediterranean onshore–offshore gradient in the composition and temporal turnover of benthic molluscs across the middle Piacenzian Warm Period. In Nawrot, R., Dominici, S., Tomašových, A., and Zuschin, M., eds. Conservation paleobiology of marine ecosystems. Geological Society of London Special Publication 529:365394.Google Scholar
Dominici, S., and Forli, M.. 2021. Lower Pliocene molluscs from southern Tuscany (Italy). Bollettino della Società Paleontologica Italiana 60:6998.Google Scholar
Dominici, S., and Scarponi, D.. 2020. Paleobiology in the Mediterranean. Bollettino della Società Paleontologica Italiana 59:17.Google Scholar
Dominici, S., Benvenuti, M., Garilli, V., Uchman, A., Pollina, F., and David, A.. 2020. Pliocene–Pleistocene stratigraphic paleobiology at Altavilla Milicia (Palermo, Sicily): tectonic, climatic and eustatic forcing. Bollettino della Società Paleontologica Italiana 59:5783.Google Scholar
Farabegoli, E., Perri, M. C., and Posenato, R.. 2007. Environmental and biotic changes across the Permian–Triassic boundary in western Tethys: the Bulla parastratotype, Italy. Global and Planetary Change 55:109135.CrossRefGoogle Scholar
Figueirido, B., Janis, C. M., Perez-Claros, J. A., de Renzi, M., and Palmqvist, P.. 2012. Cenozoic climate change influences mammalian evolutionary dynamics. Proceedings of the National Academy of Sciences USA 109:722727.CrossRefGoogle ScholarPubMed
Finnegan, S., Heim, N. A., Peters, S. E., and Fischer, W. W.. 2012. Climate change and the selective signature of the Late Ordovician mass extinction. Proceedings of the National Academy of Sciences USA 109:68296834.CrossRefGoogle ScholarPubMed
Finney, S. C., Berry, W. B. N., Cooper, J. D., Ripperdan, R., Sweet, W. C., Jacobson, S. R., Soufiane, A., Achab, A., and Noble, P. J.. 1999. Late Ordovician mass extinction: a new perspective from stratigraphic sections in central Nevada. Geology 27:215218.2.3.CO;2>CrossRefGoogle Scholar
Foote, M. 2000. Origination and extinction components of taxonomic diversity: general problems. Pp. 74102 in Erwin, D. H. and Wing, S. L., eds. Deep time: Paleobiology's perspective. Paleontological Society, Lawrence, Kans.Google Scholar
Forsythe, I. J., and Stigall, A. L.. 2023. Insights for modern invasion ecology from biotic changes of the Clarksville Phase of the Richmondian Invasion (Ordovician, Katian). Paleobiology 49:493508.CrossRefGoogle Scholar
Garilli, V. 2011. Mediterranean Quaternary interglacial molluscan assemblages: palaeobiogeograpical and palaeoceanographical responses to climate change. Palaeogeography, Palaeoclimatology, Palaeoecology 312:98114.CrossRefGoogle Scholar
Gastaldo, R. A., and Bamford, M. K.. 2023. The influence of taphonomy and time on the paleobotanical record of the Permian–Triassic transition of the Karoo Basin (and elsewhere). Journal of African Earth Sciences 204:104960.CrossRefGoogle Scholar
Gastaldo, R. A., and Demko, T. M.. 2011. The relationship between continental landscape evolution and the plant-fossil record: long term hydrologic controls on preservation. Pp. 248285 in Allison, P. A. and Bottjer, D. J., eds. Taphonomy: process and bias through time. Springer, Dordrecht.Google Scholar
Gauch, H. G. Jr., and Whittaker, R. H.. 1976. Simulation of community patterns. Vegetatio 33:1316.CrossRefGoogle Scholar
Geyman, E. C., and Maloof, A. C.. 2021. Facies control on carbonate δ13C on the Great Bahama Bank. Geology 49:10491054.CrossRefGoogle Scholar
Granot, R. 2016. Palaeozoic oceanic crust preserved beneath the eastern Mediterranean. Nature Geoscience 9:701705.CrossRefGoogle Scholar
Hansen, J., Ruedy, R., Sato, M., and Lo, K.. 2010. Global surface temperature change. Reviews of Geophysics 48:RG4004.CrossRefGoogle Scholar
Harper, D. A. T., Hammarlund, E. U., and Rasmussen, C. M. Ø.. 2014. End Ordovician extinctions: a coincidence of causes. Gondwana Research 25:12941307.CrossRefGoogle Scholar
Hay, C. C., Creveling, J. R., Hagen, C. J., Maloof, A. C., and Huybers, P.. 2019. A library of early Cambrian chemostratigraphic correlations from a reproducible algorithm. Geology 47:457460.CrossRefGoogle Scholar
Heim, N. A., and Peters, S. E.. 2011. Covariation in macrostratigraphic and macroevolutionary patterns in the marine record of North America. Geological Society of America Bulletin 123:630640.CrossRefGoogle Scholar
Hill, M. O., and Gauch, H. G. Jr. 1980. Detrended correspondence analysis: an improved ordination technique. Vegetatio 42:4758.CrossRefGoogle Scholar
Holland, S. M. 1995. The stratigraphic distribution of fossils. Paleobiology 21:92109.CrossRefGoogle Scholar
Holland, S. M. 1997. Using time–environment analysis to investigate faunal events in the Upper Ordovician of the Cincinnati Arch. Pp. 309334 in Brett, C. E., ed. Paleontological event horizons: ecological and evolutionary implications. Columbia University Press, New York.Google Scholar
Holland, S. M. 2000. The quality of the fossil record: a sequence stratigraphic perspective. Pp. 148168 in Erwin, D. H. and Wing, S. L., eds. Deep time: Paleobiology's perspective. Paleontological Society, Lawrence, Kans.Google Scholar
Holland, S. M. 2010. Additive diversity partitioning in palaeobiology: revisiting Sepkoski's question. Palaeontology 53:12371254.CrossRefGoogle Scholar
Holland, S. M. 2016. The non-uniformity of fossil preservation. Philosophical Transactions of the Royal Society of London B 371:20150130.CrossRefGoogle ScholarPubMed
Holland, S. M. 2018. Diversity and tectonics: predictions from neutral theory. Paleobiology 44:219236.CrossRefGoogle Scholar
Holland, S. M. 2020. The stratigraphy of mass extinctions and recoveries. Annual Review of Earth and Planetary Sciences 48:7597.CrossRefGoogle Scholar
Holland, S. M. 2022a. strataR: an R package for simulating the nonmarine stratigraphic and fossil record. https://github.com/hms-github/strataR.Google Scholar
Holland, S. M. 2022b. The structure of the nonmarine fossil record: predictions from a coupled stratigraphic–paleoecological model of a coastal basin. Paleobiology 48:372396.CrossRefGoogle Scholar
Holland, S. M. 2023. The contrasting controls on the occurrence of fossils in marine and nonmarine systems. Bollettino della Società Paleontologica Italiana 62:125.Google Scholar
Holland, S. M., and Loughney, K. M.. 2021. The stratigraphic paleobiology of nonmarine systems. Elements of paleontology. Cambridge University Press, Cambridge.CrossRefGoogle Scholar
Holland, S. M., and Patzkowsky, M. E.. 1999. Models for simulating the fossil record. Geology 27:491494.2.3.CO;2>CrossRefGoogle Scholar
Holland, S. M., and Patzkowsky, M. E.. 2002. Stratigraphic variation in the timing of first and last occurrences. Palaios 17:134146.2.0.CO;2>CrossRefGoogle Scholar
Holland, S. M., and Patzkowsky, M. E.. 2004. Ecosystem structure and stability: middle Upper Ordovician of central Kentucky, USA. Palaios 19:316331.2.0.CO;2>CrossRefGoogle Scholar
Holland, S. M., and Patzkowsky, M. E.. 2007. Gradient ecology of a biotic invasion: biofacies of the type Cincinnatian Series (Upper Ordovician), Cincinnati, Ohio region, USA. Palaios 22:392407.CrossRefGoogle Scholar
Holland, S. M., and Patzkowsky, M. E.. 2015. The stratigraphy of mass extinction. Palaeontology 58:903924.CrossRefGoogle Scholar
Holland, S. M., and Zaffos, A.. 2011. Niche conservatism along an onshore–offshore gradient. Paleobiology 37:270286.CrossRefGoogle Scholar
Holland, S. M., Miller, A. I., Meyer, D. L., and Dattilo, B. F.. 2001. The detection and importance of subtle biofacies within a single lithofacies: the Upper Ordovician Kope Formation of the Cincinnati, Ohio region. Palaios 16:205217.2.0.CO;2>CrossRefGoogle Scholar
Holland, S. M., Loughney, K. M., and Cone, M.. 2023. Preferential preservation of low-elevation biotas in the nonmarine fossil record. Geology 51:111114.CrossRefGoogle Scholar
Hopkins, M. J., Simpson, C., and Kiessling, W.. 2014. Differential niche dynamics among major marine invertebrate clades. Ecology Letters 17:314323.CrossRefGoogle ScholarPubMed
Horton, B. P., Edwards, R. J., and Lloyd, J. M.. 2000. Implications of a microfossil-based transfer function in Holocene sea-level studies. In Shennan, I. and Andrews, J., eds. Holocene land–ocean interaction and environmental change around the North Sea. Geological Society Special Publications 166:4154.Google Scholar
Hubbell, S. P. 2001. The unified neutral theory of biodiversity and biogeography. Princeton University Press, Princeton, N.J.Google Scholar
Huntley, J. W., and Scarponi, D.. 2015. Geographic variation of parasitic and predatory traces on mollusks in the northern Adriatic Sea, Italy: implications for the stratigraphic paleobiology of biotic interactions. Paleobiology 41:134153.CrossRefGoogle Scholar
Hutton, E. W. H., and Syvitski, J. P. M.. 2008. Sedflux 2.0: an advanced process–response model that generates three-dimensional stratigraphy. Computers and Geosciences 34:13191337.CrossRefGoogle Scholar
Ivany, L. C., Brett, C. E., Wall, H. L. B., Wall, P. D., and Handley, J. C.. 2009. Relative taxonomic and ecologic stability in Devonian marine faunas of New York State: a test of coordinated stasis. Paleobiology 35:499524.CrossRefGoogle Scholar
Jablonski, D. J., and Edie, S. M.. 2023. Perfect storms shape biodiversity in time and space. Evolutionary Journal of the Linnean Society 2:19.CrossRefGoogle Scholar
Jablonski, D. J., Roy, R. K., and Valentine, J. W.. 2006. Out of the tropics: evolutionary dynamics of the latitudinal diversity gradient. Science 314:102106.CrossRefGoogle ScholarPubMed
Jarochowska, E., Munnecke, A., Frisch, K., Ray, D. C., and Castagner, A.. 2016. Faunal and facies changes through the mid Homerian (late Wenlock, Silurian) positive carbon isotope excursion in Podolia, western Ukraine. Lethaia 49:170198.Google Scholar
Jarochowska, E., Ray, D. C., Röstel, P., Worton, G., and Munnecke, A.. 2018. Harnessing stratigraphic bias at the section scale: conodont diversity in the Homerian (Silurian) of the Midland Platform, England. Palaeontology 61:5776.CrossRefGoogle Scholar
Johnson, K. R., and Ellis, B.. 2002. A tropical rainforest in Colorado 1.4 million years after the Cretaceous–Tertiary boundary. Science 296:23792383.CrossRefGoogle ScholarPubMed
Johnson, M. E. 1988. Why are ancient rocky shores so uncommon? Journal of Geology 96:469480.CrossRefGoogle Scholar
Johnson, M. E., Skinner, D. F., and MacLeod, K. G.. 1988. Ecological zonation during the carbonate transgression of a Late Ordovician rocky shore (northeastern Manitoba, Hudson Bay, Canada). Palaeogeography, Palaeoclimatology, Palaeoecology 65:93114.CrossRefGoogle Scholar
Jongman, R. H. G., Ter Braak, C. J. F., and Van Tongeren, O. F. R., eds. 1995. Data analysis in community and landscape ecology. Cambridge University Press, Cambridge.CrossRefGoogle Scholar
Kidwell, S. M. 1985. Palaeobiological and sedimentological implications of fossil concentrations. Nature 318:457460.CrossRefGoogle Scholar
Kidwell, S. M. 1986. Models for fossil concentrations: paleobiologic implications. Paleobiology 12:624.CrossRefGoogle Scholar
Kidwell, S. M. 1991a. Condensed deposits in siliciclastic sequences: expected and observed features. Pp. 682695 in Einsele, G., Ricken, W., and Seilacher, A., eds. Cycles and events in stratigraphy. Springer-Verlag, Berlin.Google Scholar
Kidwell, S. M. 1991b. The stratigraphy of shell concentrations. Pp. 211290 in Allison, P. A. and Briggs, D. E. G., eds. Taphonomy: releasing the data locked in the fossil record. Plenum, New York.CrossRefGoogle Scholar
Kidwell, S. M., and Jablonski, D.. 1983. Taphonomic feedback: ecological consequences of shell accumulation. Pp. 195248 in Tevesz, M. J. S. and McCall, P. L., eds. Biotic interactions in recent and fossil benthic communities. Plenum, New York.CrossRefGoogle Scholar
Kiessling, W., and Aberhan, M.. 2007. Environmental determinants of marine benthic biodiversity dynamics through Triassic–Jurassic time. Paleobiology 33:414434.CrossRefGoogle Scholar
Kowalewski, M., Wittmer, J. M., Dexter, T. A., Amorosi, A., and Scarponi, D.. 2015. Differential responses of marine communities to natural and anthropogenic changes. Proceedings of the Royal Society of London B 282:20142990.Google ScholarPubMed
Kruskal, J. B. 1964. Multidimensional scaling by optimizing goodness of fit to a nonmetric hypothesis. Psychometrika 29:127.CrossRefGoogle Scholar
Lamsdell, J. C., Congreve, C. R., Hopkins, M. J., Krug, A. Z., and Patzkowsky, M. E.. 2017. Phylogenetic paleoecology: tree-thinking and ecology in deep time. Trends in Ecology and Evolution 32:452463.CrossRefGoogle ScholarPubMed
Landing, E., Westrop, S. R., Kröger, B., and English, A. M.. 2011. Left behind—delayed extinction and a relict trilobite fauna in the Cambrian–Ordovician boundary succession (east Laurentian platform, New York). Geological Magazine 148:529557.CrossRefGoogle Scholar
Legendre, P., and Legendre, L.. 1998. Numerical ecology, 2nd English ed. Elsevier, Amsterdam.Google Scholar
Leigh, D. S. 2016. Multi-millennial record of erosion and fires in the Southern Blue Ridge Mountains, USA. Pp. 167202 in Greenberg, C. H. and Collins, B. S., eds. Natural disturbances and historic range of variation. Springer International, Switzerland.CrossRefGoogle Scholar
Ling, M., Zhan, R., Wang, G., Wang, Y., Amelin, Y., Tang, P., Liu, J., et al. 2019. An extremely brief end Ordovician mass extinction linked to abrupt onset of glaciation. Solid Earth Sciences 4:190198.CrossRefGoogle Scholar
Loughney, K. M., and Holland, S. M.. 2023. The eight architectural zones of nonmarine basins. Geosphere 19:16671689.CrossRefGoogle Scholar
Loughney, K. M., Fastovsky, D. E., and Parker, W. G.. 2011. Vertebrate fossil preservation in blue paleosols from the Petrified Forest National Park, Arizona, with implications for vertebrate biostratigraphy in the Chinle Formation. Palaios 26:700719.CrossRefGoogle Scholar
Loughney, K. M., Badgley, C., Bahadori, A., Holt, W. E., and Rasbury, E. T.. 2021. Tectonic influence on Cenozoic mammal richness and sedimentation history of the Basin and Range, western North America. Science Advances 7:113.CrossRefGoogle ScholarPubMed
Lovell, B. 2010. A pulse in the planet: regional control of high-frequency changes in relative sea level by mantle convection. Journal of the Geological Society 167:637648.CrossRefGoogle Scholar
Lucas, S. G. 2017. The best sections method of studying mass extinctions. Lethaia 50:465466.CrossRefGoogle Scholar
Lyson, T. R., Miller, I. M., Bercovici, A. D., Weissenburger, K., Fuentes, A. J., Clyde, W. C., Hagadorn, J. W., et al. 2019. Exceptional continental record of biotic recovery after the Cretaceous–Paleogene mass extinction. Science 366:977983.CrossRefGoogle ScholarPubMed
Martinsen, O. J., Ryseth, A., Helland-Hansen, W., Flesche, H., Torkildsen, G., and Idil, S.. 1999. Stratigraphic base level and fluvial architecture: Ericson Sandstone (Campanian), Rock Springs Uplift, SW Wyoming, USA. Sedimentology 46:235259.CrossRefGoogle Scholar
Miller, A. I. 2000. Conversations about Phanerozoic global diversity. Pp. 5373 in Erwin, D. H. and Wing, S. L., eds. Deep time: Paleobiology's perspective. Paleontological Society, Lawrence, Kans.Google Scholar
Miller, A. I., Holland, S. M., Meyer, D. L., and Dattilo, B. F.. 2001. The use of faunal gradient analysis for intraregional correlation and assessment of changes in sea-floor topography in the type Cincinnatian. Journal of Geology 109:603613.CrossRefGoogle Scholar
Minchin, P. R. 1989. Montane vegetation of the Mt. Field massif, Tasmania: a test of some hypotheses about properties of community patterns. Vegetatio 83:97110.CrossRefGoogle Scholar
Nawrot, R., Scarponi, D., Azzarone, M., Dexter, T. A., Kusnerik, K. M., Wittmer, J. M., Amorosi, A., and Kowalewski, M.. 2018. Stratigraphic signatures of mass extinctions: ecological and sedimentary determinants. Proceedings of the Royal Society of London B 285:20181191.Google ScholarPubMed
Neal, J., and Abreu, V.. 2009. Sequence stratigraphy hierarchy and the accommodation succession method. Geology 37:779782.CrossRefGoogle Scholar
Newell, N. D., and Colbert, E. H.. 1948. Paleontologist—biologist or geologist? Journal of Paleontology 22:264267.Google Scholar
Nyberg, B., and Howell, J. A.. 2015. Is the present the key to the past? A global characterization of modern sedimentary basins. Geology 43:643646.CrossRefGoogle Scholar
Nyberg, B., and Howell, J. A.. 2016. Global distribution of modern shallow marine shorelines: implications for exploration and reservoir analogue studies. Marine and Petroleum Geology 71:83104.CrossRefGoogle Scholar
Olszewski, T. D., and Patzkowsky, M. E.. 2001. Measuring recurrence of marine biotic gradients: a case study from the Pennsylvanian–Permian Midcontinent. Palaios 16:444460.2.0.CO;2>CrossRefGoogle Scholar
Orme, C. D. L., Davies, R. G., Burgess, M., Eigenbrod, F., Pickup, N., Olson, V. A., Webster, A. J., et al. 2005. Global hotspots of species richness are not congruent with endemism or threat. Nature 436:10161019.CrossRefGoogle ScholarPubMed
Palmer, A. R. 1984. The biomere problem: evolution of an idea. Journal of Paleontology 58:599601.Google Scholar
Patzkowsky, M. E. 1995. Gradient analysis of Middle Ordovician brachiopod biofacies: biostratigraphic, biogeographic, and macroevolutionary implications. Palaios 10:154179.CrossRefGoogle Scholar
Patzkowsky, M. E. 2017. Origin and evolution of regional biotas: a deep time perspective. Annual Review of Earth and Planetary Sciences 45:471495.CrossRefGoogle Scholar
Patzkowsky, M. E., and Holland, S. M.. 1993. Biotic response to a middle Ordovician paleoceanographic event in eastern North America. Geology 21:619622.2.3.CO;2>CrossRefGoogle Scholar
Patzkowsky, M. E., and Holland, S. M.. 1997. Patterns of turnover in Middle and Upper Ordovician brachiopods of the eastern United States: a test of coordinated stasis. Paleobiology 23:420443.CrossRefGoogle Scholar
Patzkowsky, M. E., and Holland, S. M.. 1999. Biofacies replacement in a sequence stratigraphic framework: Middle and Upper Ordovician of the Nashville Dome, Tennessee, USA. Palaios 14:301323.CrossRefGoogle Scholar
Patzkowsky, M. E., and Holland, S. M.. 2003. Lack of community saturation at the beginning of the Paleozoic plateau: the dominance of regional over local processes. Paleobiology 29:545560.2.0.CO;2>CrossRefGoogle Scholar
Patzkowsky, M. E., and Holland, S. M.. 2007. Diversity partitioning of a Late Ordovician marine biotic invasion: controls on diversity in regional ecosystems. Paleobiology 33:295309.CrossRefGoogle Scholar
Patzkowsky, M. E., and Holland, S. M.. 2012. Stratigraphic paleobiology: understanding the distribution of fossil taxa in time and space. University of Chicago Press, Chicago.CrossRefGoogle Scholar
Patzkowsky, M. E., and Holland, S. M.. 2016. Biotic invasion, niche stability, and the assembly of regional biotas in deep time: comparison between faunal provinces. Paleobiology 42:359379.CrossRefGoogle Scholar
Pérès, J. 1967. The Mediterranean benthos. Oceanography and Marine Biology Annual Review 5:449533.Google Scholar
Pérès, J., and Picard, J.. 1964. Nouveau manuel de bionomie benthique de la mer Méditerranée, Edition revue et augmentée. Recueil des Travaux de la Station Marine d'Endoume 31:1137.Google Scholar
Peters, S. E. 2005. Geologic constraints on the macroevolutionary history of marine animals. Proceedings of the National Academy of Sciences USA 102:1232612331.CrossRefGoogle ScholarPubMed
Peters, S. E. 2006. Macrostratigraphy of North America. Journal of Geology 114:391412.CrossRefGoogle Scholar
Peters, S. E., and Foote, M.. 2001. Biodiversity in the Phanerozoic: a reinterpretation. Paleobiology 27:583601.2.0.CO;2>CrossRefGoogle Scholar
Peters, S. E., and Gaines, R. R.. 2012. Formation of the “Great Unconformity” as a trigger for the Cambrian Explosion. Nature 484:363366.CrossRefGoogle ScholarPubMed
Peters, S. E., and Heim, N. A.. 2011. Stratigraphic distribution of marine fossils in North America. Geology 39:259262.CrossRefGoogle Scholar
Peters, S. E., Kelly, D. C., and Fraass, A. J.. 2013. Oceanographic controls on the diversity and extinction of planktonic foraminifera. Nature 493:398403.CrossRefGoogle ScholarPubMed
Peters, S. E., Husson, J. M., and Wilcots, J.. 2017. The rise and fall of stromatolites in shallow marine environments. Geology 45:487490.CrossRefGoogle Scholar
Peters, S. E., Husson, J. M., and Czaplewski, J.. 2018. Macrostrat: a platform for geological data integration and deep-time earth crust research. Geochemistry, Geophysics, Geosystems 19:13931409.Google Scholar
Peters, S. E., Quinn, D. P., Husson, J. M., and Gaines, R. R.. 2022. Macrostratigraphy: insights into cyclic and secular evolution of the Earth–life system. Annual Review of Earth and Planetary Sciences 50:419449.CrossRefGoogle Scholar
Pier, J. Q., Brisson, S. K., Beard, J. A., Hren, M. T., and Bush, A. M.. 2021. Accelerated mass extinction in an isolated biota during Late Devonian climate changes. Scientific Reports 11:24366.CrossRefGoogle Scholar
Prave, A. R., Kirsimäe, K., Lepland, A., Fallick, A. E., Kreitsmann, T., Deines, Yu. E., Romashkin, A. E., et al. 2022. The grandest of them all: the Lomagundi–Jatuli Event and Earth's oxygenation. Journal of the Geological Society 179:jgs2021–036.CrossRefGoogle Scholar
Pruss, S. B., Jones, D. S., Fike, D. A., Tosca, N. J., and Wignall, P. B.. 2019. Marine anoxia and sedimentary mercury enrichments during the Late Cambrian SPICE event in northern Scotland. Geology 47:475478.CrossRefGoogle Scholar
Pulsipher, M. A., Schiffbauer, J. D., Jeffrey, M. J., Huntley, J. W., Fike, D. A., and Shelton, K. L.. 2021. A meta-analysis of the Steptoean Positive Carbon Isotope Excursion: the SPICEraq database. Earth-Science Reviews 212:103442.CrossRefGoogle Scholar
Raffi, S., Stanley, S. M., and Marasti, R.. 1985. Biogeographic patterns and Plio-Pleistocene extinction of Bivalvia in the Mediterranean and southern North Sea. Paleobiology 11:368388.CrossRefGoogle Scholar
Railsback, L. B., Holland, S. M., Hunter, D. M., Jordan, E. M., Diaz, J. R., and Crowe, D. E.. 2003. Controls on geochemical expression of subaerial exposure in Ordovician limestones from the Nashville Dome, Tennessee, U.S.A. Journal of Sedimentary Research 73:790805CrossRefGoogle Scholar
Railsback, L. B., Brook, G. A., Liang, F., Voarintsoa, N. R. G., Cheng, H., and Edwards, R. L.. 2018. A multi-proxy climate record from a northwestern Botswana stalagmite T suggesting wetness late in the Little Ice Age (1810–1820 CE) and drying thereafter in response to changing migration of the tropical rain belt or ITCZ. Palaeogeography, Palaeoclimatology, Palaeoecology 506:139153.CrossRefGoogle Scholar
Raup, D. M. 1976. Species diversity in the Phanerozoic: an interpretation. Paleobiology 2:289297.CrossRefGoogle Scholar
Reddin, C. J., Kocsis, Á. T., and Kiessling, W.. 2019. Climate change and the latitudinal selectivity of ancient marine extinctions. Paleobiology 45:7084.CrossRefGoogle Scholar
Regan, A. K., Rogers, R. R., and Holland, S. M.. 2022. Quantifying controls on the occurrence of nonmarine fossils. Geology 50:12871290.CrossRefGoogle Scholar
Ricklefs, R. E. 2006. Evolutionary diversification and the origin of the diversity–environment relationship. Ecology 87:S3S13.CrossRefGoogle ScholarPubMed
Ricklefs, R. E. 2008. Disintegration of the ecological community. American Naturalist 172:741750.CrossRefGoogle ScholarPubMed
Ritter, M. D. N., Erthal, F., Kosnik, M. A., Kowalewski, M., Coimbra, J. C., Caron, F., and Kaufman, D. S.. 2023. Onshore–offshore trends in the temporal resolution of molluscan death assemblages: how age–frequency distributions reveal Quaternary sea-level history. Palaios 38:148157.CrossRefGoogle Scholar
Rogers, R. R., and Brady, M. E.. 2010. Origins of microfossil bonebeds: insights from the Upper Cretaceous Judith River Formation of north–central Montana. Paleobiology 36:80112.CrossRefGoogle Scholar
Rogers, R. R., and Kidwell, S. M.. 2000. Associations of vertebrate skeletal concentrations and discontinuity surfaces in continental and shallow marine records: a test in the Cretaceous of Montana. Journal of Geology 108:131154.CrossRefGoogle Scholar
Rogers, R. R., and Kidwell, S. M.. 2007. A conceptual framework for the genesis and analysis of vertebrate skeletal concentrations. Pp. 163 in Rogers, R. R., Eberth, D. A., and Fiorillo, A. R., eds. Bonebeds: genesis, analysis, and paleobiological significance. University of Chicago Press, Chicago.CrossRefGoogle Scholar
Rogers, R. R., Eberth, D. A., and Ramezani, J.. 2023. The “Judith River–Belly River problem” revisited (Montana-Alberta-Saskatchewan): new perspectives on the correlation of Campanian dinosaur-bearing strata based on a revised stratigraphic model updated with CA-ID-TIMS U-Pb geochronology. GSA Bulletin. https://doi.org/10.1130/B36999.1.CrossRefGoogle Scholar
Ross, R. J. Jr. 1986. Who's holding the watch? Palaios 1:1.CrossRefGoogle Scholar
Rowley, D. B. 2018. Oceanic axial depth and age–depth distribution of oceanic lithosphere: comparison of magnetic anomaly picks versus age–grid models. Geosphere 11:2143.Google Scholar
Saltzman, M. R., Edwards, C. T., Adrain, J. M., and Westrop, S. R.. 2015. Persistent oceanic anoxia and elevated extinction rates separate the Cambrian and Ordovician radiations. Geology 43:807810.CrossRefGoogle Scholar
Saupe, E. E. 2023. Explanations for latitudinal diversity gradients must invoke rate variation. Proceedings of the National Academy of Sciences USA 120:e2306220120.CrossRefGoogle ScholarPubMed
Saupe, E. E., Hendricks, J. R., Portell, R. W., Dowsett, H. J., Haywood, A., Hunter, S. J., and Lieberman, B. S.. 2014. Macroevolutionary consequences of profound climate change on niche evolution in marine molluscs over the past three million years. Proceedings of the Royal Society of London B 281:20141995.Google ScholarPubMed
Scarponi, D., and Kowalewski, M.. 2004. Stratigraphic paleoecology: bathymetric signatures and sequence overprint of mollusk associations from upper Quaternary sequences of the Po Plain, Italy. Geology 32:989992.CrossRefGoogle Scholar
Scarponi, D., and Kowalewski, M.. 2007. Sequence stratigraphic anatomy of diversity patterns: Late Quaternary benthic mollusks of the Po Plain, Italy. Palaios 22:296305.CrossRefGoogle Scholar
Scarponi, D., Kaufman, D., Amorosi, A., and Kowalewski, M.. 2013. Sequence stratigraphy and the resolution of the fossil record. Geology 41:239242.CrossRefGoogle Scholar
Scarponi, D., Azzarone, M., Kusnerik, K. M., Amorosi, A., Bohacs, K., Drexler, T. M., and Kowalewski, M.. 2017. Systematic vertical and lateral changes in quality and time resolution of the macrofossil record: insights from Holocene transgressive deposits, Po coastal plain, Italy. Marine and Petroleum Geology 87:128136.CrossRefGoogle Scholar
Scotese, C. R., Song, H., Mills, B. J. W., van der Meer, D. G.. 2021. Phanerozoic paleotemperatures: the Earth's changing climate during the last 540 million years. Earth-Science Reviews 215:103503.CrossRefGoogle Scholar
Sepkoski, J. J. Jr. 1987. Environmental trends in extinction during the Paleozoic. Science 235:6466.CrossRefGoogle ScholarPubMed
Sepkoski, J. J. Jr. 1988. Alpha, beta, or gamma: where does all the diversity go? Paleobiology 14:221234.CrossRefGoogle ScholarPubMed
Sepkoski, J. J. Jr. 1991. A model of onshore–offshore change in faunal diversity. Paleobiology 17:5877.CrossRefGoogle Scholar
Shaw, A.B. 1971. The butterfingered handmaiden. Journal of Paleontology 45:15.Google Scholar
Sloss, L. L. 1963. Sequences in the cratonic interior of North America. Geological Society of America Bulletin 74:93114.CrossRefGoogle Scholar
Sloss, L. L. 1976. Areas and volumes of cratonic sediments, western North America and eastern Europe. Geology 4:272276.2.0.CO;2>CrossRefGoogle Scholar
Smith, A. B., Gale, A. S., and Monks, N. E. A.. 2001. Sea-level change and rock-record bias in the Cretaceous: a problem for extinction and biodiversity studies. Paleobiology 27:241253.2.0.CO;2>CrossRefGoogle Scholar
Song, H., Kemp, D. B., Tian, L., Chu, D., Song, H., and Dai, X.. 2021. Thresholds of temperature change for mass extinctions. Nature Communications 12:4694.CrossRefGoogle ScholarPubMed
Taylor, P. D., and Wilson, M. A.. 2003. Palaeoecology and evolution of marine hard substrate communities. Earth-Science Reviews 62:1103.CrossRefGoogle Scholar
Therrien, F., and Fastovsky, D. E.. 2000. Paleoenvironments of early theropods, Chinle Formation (Late Triassic), Petrified Forest National Park, Arizona. Palaios 15:194211.2.0.CO;2>CrossRefGoogle Scholar
Tobin, T. S. 2017. Recognition of a likely two phased extinction at the K–Pg boundary in Antarctica. Scientific Reports 7:16317.CrossRefGoogle ScholarPubMed
Tomašových, A., Fürsich, F. T., and Olszewski, T. D.. 2006a. Modeling shelliness and alteration in shell beds: variation in hardpart input and burial ratees leads to opposing predictions. Paleobiology 32:278298.CrossRefGoogle Scholar
Tomašových, A., Fürsich, F. T., and Wilmsen, M.. 2006b. Preservation of autochthonous shell beds by positive feedback between increased hardpart-input rates and increased sedimentation rates. Journal of Geology 114:287312.CrossRefGoogle Scholar
Tomašových, A., Dominici, S., Zuschin, M., and Merle, D.. 2014. Onshore–offshore gradient in metacommunity turnover emerges only over macroevolutionary time-scales. Proceedings of the Royal Society of London B 281:20141533.Google ScholarPubMed
Tomašových, A., Kidwell, S. M., and Barber, R. F.. 2016. Inferring skeletal production from time-averaged assemblages: skeletal loss pulls the timing of production pulses towards the modern period. Paleobiology 42:5476.CrossRefGoogle Scholar
Tomašových, A., Kidwell, S. M., Alexander, C.R., and Kaufman, D. S.. 2019. Millennial-scale age offsets within fossil assemblages: result of bioturbation below the taphonomic active zone and out-of-phase production. Paleoceanography and Paleoclimatology 34:954977.CrossRefGoogle Scholar
Tomašových, A., Albano, P. G., Fuksi, T., Gallmetzer, I., Haselmair, A., Kowalewski, M., Nawrot, R., Nerlović, V., Scarponi, D., and Zuschin, M.. 2020. Ecological regime shift preserved in the Anthropocene stratigraphic record. Proceedings of the Royal Society of London B 287:20200695.Google ScholarPubMed
Tomašových, A. I. Gallmetzer, A. Haselmair, and M. Zuschin. 2022. Inferring time averaging and hiatus durations in the stratigraphic record of high-frequency depositional sequences. Sedimentology 69:10831118.CrossRefGoogle Scholar
Tomašových, A., Dominici, S., Nawrot, R., and Zuschin, M.. 2023a. Temporal scales, sampling designs and age distributions in marine conservation paleobiology. In Nawrot, R., Dominici, S., Tomašových, A., and Zuschin, M., eds. Conservation paleobiology of marine ecosystems. Geological Society of London Special Publication 529:139.Google Scholar
Tomašových, A., Kidwell, S. M., and Dai, R.. 2023b. A downcore increase in time averaging is the null expectation from the transit of death assemblages through a mixed layer. Paleobiology 49:527562.CrossRefGoogle Scholar
Trampush, S. M., and Hajek, E. A.. 2017. Preserving proxy records in dynamic landscapes: modeling and examples from the Paleocene–Eocene Thermal Maximum. Geology 45:967970.CrossRefGoogle Scholar
Van Wagoner, J. C., Mitchum, R. M., Campion, K. M., and Rahmanian, V. D.. 1990. Siliciclastic sequence stratigraphy in well logs, cores, and outcrops. AAPG Methods in Exploration Series, No. 7. American Association of Petroleum Geologists, Tulsa, Okla.CrossRefGoogle Scholar
Veizer, J., and Prokoph, A.. 2015. Temperatures and oxygen isotopic composition of Phanerozoic oceans. Earth-Science Reviews 146:92104.CrossRefGoogle Scholar
Wagner, P. J., and Marcot, J. D.. 2013. Modelling distributions of fossil sampling rates over time, space and taxa: assessment and implications for macroevolutionary studies. Methods in Ecology and Evolution 4:703713.CrossRefGoogle Scholar
Westerhold, T., Marwan, N., Drury, A. J., Liebrand, D., Agnini, C., Anagnostou, E., Barnet, J. S. K., et al. 2020. An astronomically dated record of Earth's climate and its predictability over the last 66 million years. Science 369:13831387.CrossRefGoogle ScholarPubMed
Whittaker, R. H. 1967. Gradient analysis of vegetation. Biological Reviews 49:207264.CrossRefGoogle Scholar
Wilf, P., Carvalho, M. R., and Stiles, E.. 2023. The end-Cretaceous plant extinction: heterogeneity, ecosystem transformation, and insights for the future. Cambridge Prisms: Extinction 1:e14.Google Scholar
Wittmer, J. M., Dexter, T. A., Scarponi, D., Amorosi, A., and Kowalewski, M.. 2014. Quantitative bathymetric models for Late Quaternary transgressive–regressive cycles of the Po Plain, Italy. Journal of Geology 122:649670.CrossRefGoogle Scholar
Ye, S., and Peters, S. E.. 2023. Bedrock geological map predictions for Phanerozoic fossil occurrences. Paleobiology 49:394413.CrossRefGoogle Scholar
Zachos, J., Pagani, M., Sloan, L., Thomas, E., and Billups, K.. 2001. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292:686693.CrossRefGoogle ScholarPubMed
Zaffos, A., Finnegan, S., and Peters, S. E.. 2017. Plate tectonic regulation of global marine animal diversity. Proceedings of the National Academy of Sciences USA 114:56535658.CrossRefGoogle ScholarPubMed
Zimmt, J. B., Holland, S. M., Finnegan, S., and Marshall, C. R.. 2021. Recognizing pulses of extinction from clusters of last occurrences. Palaeontology 64:120.CrossRefGoogle Scholar
Zuschin, M., Harhauser, M., Hengst, B., Mandic, O., and Roetzel, R.. 2014. Long-term ecosystem stability in an Early Miocene estuary. Geology 42:710.CrossRefGoogle Scholar
Zuschin, M., Nawrot, R., Harzhauser, M., Mandic, O., and Tomašových, A.. 2017. Taxonomic and numerical sufficiency in depth- and salinity-controlled marine paleocommunities. Paleobiology 43:463478.CrossRefGoogle Scholar
Figure 0

Figure 1. Coding ordination scores of samples and taxa can reveal the origin of the axes, shown with an example from the C2 sequence of the Upper Ordovician of the Cincinnati Arch, USA (data from Holland and Patzkowsky 2007). A, Sample scores coded by lithofacies reveal that detrended correspondence analysis (DCA) axis 1 is correlated with water depth. B, Taxon scores coded by life habit demonstrate that DCA axis 2 is correlated with substrate characteristics, with burrowers associated with soft muds at low axis 2 scores and attached species associated with shelly gravels at high axis 2 scores.

Figure 1

Figure 2. Stratigraphic cross sections showing systems tracts and surfaces in marine and nonmarine settings. A, Cross section is model output from the basin simulation model sedflux (Hutton and Syvitski 2008). FSST, falling-stage systems tract; HST, highstand systems tract; LST, lowstand systems tract; TST, transgressive systems tract. B, Cross section is model output from the basin simulation model strataR (Holland 2022a). Black lines are evenly spaced timelines, which can be used to infer aggradation rates. HAST, high-accommodation systems tract; LAST, low-accommodation systems tract.

Figure 2

Figure 3. Depositional strike and dip, illustrated with the coast of Egypt at the Nile Delta. Map from Google Earth Pro. Depositionally updip areas are toward the sediment source (typically landward), and depositionally downdip areas are distal to the sediment source and therefore more commonly sediment starved.

Figure 3

Figure 4. Patterns of fossil occurrence through two transgressive–regressive sequences. The two sequences each contain 10 high-frequency parasequences, upward-shallowing cycles bounded by flooding surfaces across which water depth increases rapidly. Fifty species having the same peak abundance and depth tolerance but varying in preferred depth are shown sorted by their preferred depths. Dots depict fossil occurrences; solid lines show fossil ranges; and thin dashed gray lines show when species are extant (although possibly not preserved). Fossil occurrences systematically track changes in water depth, and major flooding surfaces (blue lines) have clusters of first and last occurrences. Note that the interval of elevated extinction does not have an obviously greater number of last occurrences. Most species that go extinct during the period of elevated extinction have their last occurrences at two major flooding surfaces immediately below the window of elevated extinction. Other last occurrences are spread through the elevated extinction interval and are not distinguishable from last occurrences caused by changing water depths. The full magnitude of the extinction would not be apparent until the next transgressive–regressive sequence (not shown). The timing of the extinction can be known only by studying correlative columns along the depositional dip.