White Pine County, Nevada

The upper Cambrian succession in eastern Nevada was deposited on a distally steepened carbonate ramp (Osleger and Read, Reference Osleger and Read1993) that lay on a passive margin formed by Neoproterozoic rifting (e.g., Merdith et al., Reference Merdith, Williams, Collins, Tetley and Mulder2021). The Windfall Formation is divided into three units, in ascending order, the Barton Canyon Limestone, Catlin, and Bullwhacker members (Adrain and Westrop, Reference Adrain and Westrop2004, fig. 2). Figure 1 summarizes the stratigraphy and correlation of the study interval.

Barton Canyon Limestone Member.—The Barton Canyon Limestone Member forms a conspicuous low cliff of light gray carbonates (Fig. 8.1, 8.2) between the more recessive slopes formed by the underlying Dunderberg Formation and the overlying Catlin Member. It is fully exposed only at section CHC-2, where it is 17 m thick. The Barton Canyon Limestone Member consists mostly of thick, stylo-bedded, bioturbated lime mudstone to wackestone (Fig. 8.6) with mm- to cm-thick horizons of bioclastic packstone (Fig. 8.5). Bioturbation is recorded by irregular, often dolomitic swirls (Fig. 8.5, 8.6). An interval with cm- to dm-thick echinoderm–trilobite grainstone is present at the base of the member (Fig. 8.7), and a 10–15 cm thick condensed bioclastic shell bed at the top (Figs. 8.3, 9) includes the Irvingella “major” Zone (6.5 cm thick at STR; 7.5 cm thick at CHC-1), which defines the base of the Sunwaptan Stage and marks the onset of faunal change (e.g., Westrop and Cuggy, Reference Westrop and Cuggy1999). The Barton Canyon Limestone Member is mostly light gray in color (Fig. 8.1, 8.2), but becomes dark gray below (30 cm below at CHC-1 and 2.2 m below at STR) the I. “major” Zone. Miller et al. (Reference Miller, Evans, Dattilo, Derby, Fritz, Longacre, Morgan and Sternbach2012, p. 802) reported a similar color change in the correlative interval of the Sneakover Limestone Member of the Orr Formation in west-central Utah.

Figure 8. Barton Canyon Limestone and Catlin members of the Windfall Formation. (1) Cliff formed by the Barton Canyon Limestone Member overlain by recessive slope of the Catlin Member, section STR, North Egan Range (Fig. 1.1); figure in circle for scale. (2) Cliff formed by Barton Canyon Limestone Member at section CHC-1, Barton Canyon, Cherry Creek Range (Fig. 1.2); backpack in white circle for scale. (3) Uppermost Barton Canyon Limestone Member and lower Catlin Member, section STR; white arrow shows position of the Irvingella “major” coquina; backpack in white circle for scale. (4) Thin-bedded, dark-gray to black, silty lime mudstone, lower Catlin Member, approximately 13.4–13.7 m above the base of section STR; hammer for scale. (5) Naturally polished surface of Barton Canyon Limestone Member, STR 5.5, showing bioturbated lime wackestone and thin seams of bioclastic packstone, with abundant trilobite sclerites; fingertip for scale. (6) Polished slab of bioturbated lime mudstone to wackestone with dolomitic burrow-mottles, CHC-2-210 m. (7) Dm-thick echinoderm grainstone interbedded with bioturbated lime mudstone to wackestone, lower Barton Canyon Limestone Member, section CHC-1.

Brady and Rowell (Reference Brady and Rowell1976) characterized the Barton Canyon Limestone Member as having formed in a widespread area of shallow carbonate deposition that blanketed the shelf. They showed that it passed landward (west-central Utah) into high-energy shoal deposits of cross-bedded grainstone. To the west, in central Nevada, the Barton Canyon Limestone Member carbonates are bounded by continental slope facies (Brady and Rowell, Reference Brady and Rowell1976) that include debris flows and soft-sediment slumping (Taylor and Cook, Reference Taylor and Cook1976). The large areal extent (more than 40,000 km², according to Brady and Rowell, Reference Brady and Rowell1976) of a rather monotonous, subtidal carbonate facies between west-central Utah and east-central Nevada implies a low gradient to the carbonate ramp.

As noted by Brady and Rowell (Reference Brady and Rowell1976), lime mudstone–wackestone facies represents an open-shelf environment below fair-weather wave base. Echinoderms and trilobites record normal marine conditions. Thin concentrations of bioclasts in packstone layers do, however, indicate that there was at least minor winnowing of the sea floor, as does spar-filled shelter porosity. As such, the Barton Canyon Limestone Member was probably deposited between fair-weather and storm wave base.

In a broader context, the appearance of clean carbonates above the mixed carbonate–fine siliciclastic succession of the Dunderberg Formation implies a transgression that sequestered siliciclastic sediment in terrigenous source areas due to a rise in base level (see also Evans et al., Reference Evans, Miller, Dattilo and Swanson2003). Miller et al. (Reference Miller, Evans, Dattilo, Derby, Fritz, Longacre, Morgan and Sternbach2012) identified this transgression as the Sauk III transgression, which was a “major” flooding event that can be recognized widely in Laurentian North America. A second, equally abrupt sea-level rise terminated subtidal carbonate deposition in east-central Nevada, with the appearance of the Catlin Member above the Barton Canyon Limestone Member (see also Evans et al. [2003, p. 29 and fig. 7], who interpreted the Irvingella “major” Zone in the House Range of Utah as recording a sea-level rise, and the top of the Barton Canyon Limestone Member in the North Egan Range as a drowning unconformity [Evans et al., Reference Evans, Miller, Dattilo and Swanson2003, fig. 8E]). The onset of each of these sea-level rises is marked by bioclastic accumulations.

Bioclastic accumulations.—Bioclastic accumulations occur at the bottom and top of the member. From their stratigraphic positions at the onset of sea-level rise, they represent onlap shell beds (Zecchin, Reference Zecchin2007; Zecchin and Catuneanu, Reference Zecchin and Catuneanu2013). Such beds form under conditions of sediment starvation and usually rest on ravinement surfaces in siliciclastic successions (Zecchin, Reference Zecchin2007). Sediment starvation will also accompany drowning in carbonate systems as the carbonate factory is shut down (Catuneanu, Reference Catuneanu2022).

The basal meter of the Barton Canyon Limestone Member is bioclastic packstone, grainstone, and rudstone with dm-thick intercalations of stylolitic lime mudstone–wackestone (Fig. 8.7). Bioclastic intervals lack primary sedimentary structures like cross-bedding but yield trilobites of the Elvinia Zone. The unit is interpreted as recording higher energy conditions in the initial phase of the Sauk III transgression. Above the basal meter, in the lower half of the member, cm-thick intervals of grainstone occur as minor interbeds in the wackestone–packstone facies.

The shell bed at the very top of the Barton Canyon Limestone Member is a 10–20 cm thick, condensed interval of trilobite-rich, bioclastic grainstone, packstone, and rudstone that includes the Irvingella “major” Zone. It is overlain by the deeper water facies of the Catlin Member (Fig. 8.3). In detail, this shell bed is divided internally into four (CHC; Fig. 9.1) or three (STR; Fig. 9.2) layers that are separated by planar to irregular hardgrounds, indicating that it is strongly condensed. At both localities, the basal layers (layer 1 in Fig. 9.1, 9.2) yield trilobites of the Elvinia Zone, with the overlying layers containing the I. “major” Zone. At CHC-1 (Fig. 9.1), the Elvinia Zone (layer 1) is mostly bioturbated wackestone–packstone, capped by black lime mudstone at an irregular contact. The base of the I. “major” Zone is a thin (cm) grainstone layer (layer 2) that is succeeded by an interval of trilobite packstone to rudstone (layer 3) that comprises most of the zone. It appears to include at least one undulating internal hardground. The uppermost part of the bed (layer 4) is another thin (cm) bioclastic grainstone to rudstone.

Figure 9. Condensed shell bed (including the Irvingella “major” coquina that marks the base of the Sunwaptan Stage and the onset of faunal change) at the top of the Barton Canyon Limestone Member. Short white and black arrows show internal hardgrounds. (1) Sample CHC-1-0. Each of the four layers was sampled separately for carbon isotopes and trilobites. The upper half of the slab is the Irvingella “major” Zone and is divisible into three layers (2–4). Layer 2 is bioclastic grainstone; layer 3 is bioclastic pack- and rudstone; layer 4 is bioclastic grain- to rudstone. Layer 1 yields a trilobite fauna of the Elvinia Zone and includes light gray lime mudstone to wackestone and black lime mudstone. Trilobite abundance data for layers 1, 3, and 4 are shown in Figure 14. (2) Sample STR 11.9–12.1. Three layers of bioclastic grain- to rudstone were sampled separately for trilobites, each of which is separated by a well-defined, irregular hardground (short black arrows); layer 1 contains the of the Elvina Zone; layers 2 and 3 comprise the Irvingella “major” coquina. Trilobite abundance data for each of the layers are shown in Figure 14.

Section STR is about 48 km south of section CHC, and there are differences in the microstratigraphy of the shell beds. The basal layer (layer 1) of the shell bed at STR (Fig. 9.2) is a bioclastic grainstone to rudstone that yields the Elvinia fauna, and is separated from two overlying layers (layers 2 and 3) with the I. “major” fauna by a well-defined hardground. The I. “major” Zone consists of cm-thick bioclastic rudstone with large trilobite sclerites (predominantly Irvingella cranidia) with spar-filled shelters beneath them. The two layers are separated by a gently undulating hardground. As such, the shell bed at STR is more strongly winnowed and lacks a packstone component.

Although welded onto the shallow-water carbonates at the top of the Barton Canyon Limestone Member, the shell beds are genetically the initial deposits of the drowning succession. This indicates that the extinctions began after sea-level rise. The color change from light gray to dark gray near the top of the carbonate succession may herald the deepening that led to the platform drowning.

Catlin Member.—The Catlin Member records a marked change in facies from the underlying Barton Canyon Limestone Member (Fig. 8.1, 8.3), which is accompanied by an equally sharp turnover in the trilobite faunas, beginning with the onlap shell bed. The dominant lithology of the Catlin Member consists of cm-thick dark gray to black, cherty lime mudstone and calcisiltite, both of which include laminae of terrigenous silt (Fig. 8.4). Minor bioturbation is present, but undisturbed plane lamination becomes increasingly prevalent up section. Centimeter-thick echinoderm–trilobite grainstone and packstone interbeds yield low-diversity trilobite assemblages with abundant olenid trilobites. The facies change at the base of the Catlin Member, with reduced bioturbation and preservation of primary sedimentary lamination, is consistent with deepening and a shift towards dysoxic conditions. At the same time, the presence of trilobites indicates that there was sufficient oxygen at least periodically for respiration (e.g., Landing and Westrop, Reference Landing and Westrop2015, p. 985, 986; see also Dahl et al., Reference Dahl, Siggaard-Andersen, Schovsbo, Persson, Husted, Hougård, Dickson, Kjær and Nielsen2019, for geochemical evidence of episodic oxygenation of generally dysoxic to anoxic shale facies in Baltica) and maintenance of strongly calcified exoskeletons that, unlike those of modern marine arthropods, were rebuilt entirely without any resorption as part of the molting process (Miller and Clarkson, Reference Miller and Clarkson1980; Brandt, Reference Brandt2002). At horizons that yield trilobites, levels of dissolved oxygen must have been sufficient to meet these metabolic constraints.

Sequence stratigraphy.—The base of the Barton Canyon Limestone Member clearly marks a shut-down in the siliciclastic supply and records a major transgression that can be recognized across Laurentia in the Elvinia Zone (e.g., Lochman-Balk, Reference Lochman-Balk1970). The basal interval of the member with echinoderm–trilobite grainstone layers may represent a small-scale transgressive systems tract, with a small-scale highstand succession recorded by the overlying lime mudstone–wackestone facies. This depositional sequence is terminated by further deepening that can be interpreted as a drowning unconformity or type 3 sequence boundary, as defined by Schlager (Reference Schlager, Crevello, Wilson, Sarg and Read1989, Reference Schlager2005) (see also Evans et al., Reference Evans, Miller, Dattilo and Swanson2003). Note that Catuneanu (Reference Catuneanu2022, p. 355) considered drowning unconformities to be a type of flooding surface rather than a sequence boundary. Rather than forming at a surface of wave ravinement typical of siliciclastic successions (e.g., Zecchin and Catuneanu, Reference Zecchin and Catuneanu2013), the condensed shell bed at the top of the Barton Canyon Limestone Member (Fig. 9) sits on the drowning surface and is the product of sediment starvation as the carbonate factory shut down. The implications of stratigraphic condensation on our perception of the extinctions are explored in more detail below.

Central Oklahoma

As noted above, sampling of the Honey Creek Formation in the Slick Hills (section KR1) failed to yield a reliable carbon isotope curve. Although we were able to generate a curve from archival collections from Stitt's (Reference Stitt1971b) Royer Ranch (RR) section in the Arbuckle Mountains, changes in land ownership prevented us from gaining access to this or any other of Stitt's localities. Stitt's (Reference Stitt1971b, p. 65) lithologic log of the relevant interval of RR described the Honey Creek Formation as a single unit that is 31.1 m (102 feet) thick. Stitt's general description of the lithologies as “glauconitic trilobite–pelmatozoan biosparite and biomicrite,” and the presence of “fine sand to silt-sized quartz” indicate that the succession is broadly comparable to section KR1, which will be used to augment the facies description of the Honey Creek Formation. In the Slick Hills, deposition of the Honey Creek Formation was influenced by islands of the Southern Oklahoma Archipelago (Donovan and Stephenson, Reference Donovan, Stephenson and Johnson1991) that were sources of siliciclastics, and there are rapid lateral facies changes from siliciclastic-rich to carbonate-dominated successions (e.g., see sections through the lower Honey Creek Formation in Westrop et al., Reference Westrop, Waskiewicz Poole and Adrain2010, fig. 1). Donovan and Bucheit (Reference Donovan, Bucheit and Johnson2000, fig. 2) placed their Ring Top Island at the north end of Ring Top Mountain on the Caddo–Comanche county line, roughly 0.5 km north of section KR1 (Westrop et al., Reference Westrop, Waskiewicz Poole and Adrain2010, fig. 1). However, this island was apparently submerged prior to deposition of the Honey Creek Formation (Donovan and Bucheit, Reference Donovan, Bucheit and Johnson2000), which is consistent with the carbonate-rich succession at KR1. Sections (KR2, KR3) along the flanks of Ring Top Mountain, no more than 1 km to the southwest of KR1, have abundant sandstone (Westrop et al., Reference Westrop, Waskiewicz Poole and Adrain2010, fig. 1), some of which includes grainstone to rudstone lenses that represent starved bioclastic ripples, which is laterally equivalent to carbonate intervals in KR1. There were clearly other islands that continued to be sources of siliciclastics and, indeed, there were at least two other islands in the vicinity of Ring Top Mountain (Donovan and Bucheit, Reference Donovan, Bucheit and Johnson2000, fig. 2).

As noted by Osleger and Read (Reference Osleger and Read1993), the initial drowning in central Oklahoma is marked by the Reagan Sandstone, which sits unconformably on the Carlton Rhyolite and passes gradationally upward into the sandy, glauconitic skeletal carbonates of the Honey Creek Formation (see also Donovan et al., Reference Donovan, Ayan, Bucheit and Johnson2000). This is the general Elvinia Zone sea-level rise that is recorded by the Barton Canyon Limestone Member and correlatives in Nevada and Utah. Progressive drowning of the archipelago persisted through the Taenicephalus Zone and was not completed until deposition of the Fort Sill Formation, which overlies the Honey Creek Formation (Donovan and Bucheit, Reference Donovan, Bucheit and Johnson2000).

Honey Creek Formation.—In contrast to the sharp environmental change in the interval of extinction in east-central Nevada, the Honey Creek Formation records a continuous facies succession, and the dominant pattern is lateral change influenced by proximity to islands that supplied siliciclastic sediment. The sequence stratigraphy, sedimentary facies, and biostratigraphy of the Honey Creek Formation in the Slick Hills will be documented in detail elsewhere. However, at section KR1, pre-extinction lithologies (Fig. 10.1; see Westrop et al., Reference Westrop, Waskiewicz Poole and Adrain2010, fig. 1, for a stratigraphic column of the Elvinia Zone interval) consist of rippled, variably sandy, glauconitic, trilobite–echinoderm grainstone to rudstone with thin siliciclastic drapes that are often accentuated by pressure solution. An underlying sandstone-dominated interval includes beds with starved bioclastic ripples. Bidirectional cross-lamination is common, which Donovan (Reference Donovan and Johnson2000) and Donovan and Bucheit (Reference Donovan, Bucheit and Johnson2000) interpreted in both the Honey Creek Formation and underlying siliciclastics of the Reagan Formation in the Slick Hills as tidally influenced.

Figure 10. Honey Creek Formation, Ring Top Mountain, Kimbell Ranch, Comanche County, Oklahoma, section KR1 (see Fig. 11 for stratigraphic column; see Westrop et al., Reference Westrop, Waskiewicz Poole and Adrain2010, fig. 1, for locality information and a stratigraphic column for the pre-extinction interval). (1) Rippled, glauconitic trilobite–echinoderm grainstone with thin siliciclastic drapes, Elvinia Zone, ~10.5 m above the base of the section, and ~9.5 m below the base of the Irvingella “major” Zone; pen and hammer for scale; note cross-lamination below pen. (2) Glauconitic trilobite–brachiopod–echinoderm grain- to rudstone, Irvingella “major” Zone, 22.85–24 m above the base of the section; rectangle shows field of view in (4); pencil for scale. (3) Rippled, glauconitic trilobite–echinoderm grainstone with thin siliciclastic drapes, Taenicephalus Zone; hammer is 33.4 m above the base of the section, and 11.4 m above the base of the Irvingella “major” Zone. (4) Close-up of area of rectangle in (2) showing orthid brachiopod valves with spar-filled shelters. (5) Irvingella shell bed, I. “major” Zone, Honey Creek Formation, Bally Mountain, Slick Hills, Oklahoma, collected ~1 km south along strike from the section described by Blackwell and Westrop (Reference Blackwell and Westrop2023, fig. 1); the surface is crowded with cranidia of Irvingella; scale bar represents 1 cm. (6) Eoorthis shell bed, Parabolinoides Subzone, collection RR 142.3, Royer Ranch section, Arbuckle Mountains, Oklahoma (Stitt, Reference Stitt1971b); surface is dominated by valves of the brachiopod Eoorthis indianola (Walcott, Reference Walcott1905) (see Freeman and Stitt, Reference Freeman and Stitt1996, p. 364, 365), with scattered cranidia of the trilobite Orygmaspis (Parabolinoides); scale bar represents 1 cm.

Trilobite shell beds with abundant sclerites are present in the Elvinia Zone and extend into the I. “major” (e.g., Fig. 10.5) and Taenicephalus zones, albeit with sharp changes in genus composition. Orthid brachiopods enter the succession in the extinction interval in both the Slick Hills and the Arbuckle Mountains, where they also form shell beds (Fig. 10.4, 10.6) that are interbedded with trilobite beds (Fig. 11). Although in a different environmental setting than the onlap shell beds above the drowning unconformity at the top of the Barton Canyon Limestone Member in Nevada, formation of shell accumulations in the Honey Creek Formation still implies some degree of condensation or sediment starvation (e.g., Datillo et al., Reference Dattilo, Brett, Tsujita and Fairhurst2008). In other respects, there is little change in the environment, with trilobite–echinoderm grainstone as the dominant lithology (Fig. 10.2) that persists above the extinction interval (Fig. 10.3). An exception at KR1 is marked by the appearance of dm-thick interbeds of bioturbated sandstone through a five-meter interval in the upper Taenicephalus Zone (Fig. 11). This is interpreted as a progradational package recording a minor relative sea-level fall.

Figure 11. Stratigraphic column for the upper Elvinia Zone–Taenicephalus Zone interval of the Honey Creek Formation at locality KR1, showing the distribution of trilobite genera, trilobite shell beds, and orthid brachiopod shell beds. Orthid shell beds have a limited stratigraphic distribution of ~7.5 m, appearing in the I. “major” Zone and disappearing in the upper part of the Taenicephalus Zone.

Biofacies, abundances, and diversity changes

Westrop and Cuggy (Reference Westrop and Cuggy1999) presented a quantitative analysis of biofacies changes at the three Cambrian trilobite extinction events (bases of the Steptoean, Sunwaptan, and Skullrockian stages). They showed that ecological differentiation of trilobite faunas declined in the extinction intervals, along with a steady fall in species diversity. We employed a similar methodology using mostly new field collections to compare faunal turnover in Oklahoma and Nevada. Because the species-level systematics of the trilobites and associated agnostid arthropods is still under study, the analysis is conducted at the genus level. We agree with Hendricks et al. (Reference Hendricks, Saupe, Myers, Hermsen and Allmon2014) that ecological and evolutionary inferences from genus-level patterns need to be made carefully. For all of our new collections and those derived from the literature, abundances of genera were calculated using the minimum number of individuals method (Gilinsky and Bennington, Reference Gilinsky and Bennington1994).

Cluster analysis of biofacies.— Figure 12 shows the results of cluster analysis of log-transformed abundance data using Euclidian distance as a dissimilarity metric and Ward's method for linkage (the untransformed abundance data are presented in Table 1). This combination is a robust, space-conserving method that minimizes chaining in the dendrogram (McCune and Grace, Reference McCune and Grace2002).

Figure 12. Two-way cluster analysis with collections in Q-mode order and genera in R-mode order. Analysis was performed in PAST (Hammer et al., Reference Hammer, Harper and Ryan2001), using Ward's method and log-transformed genus abundances; genus abundances are expressed as percentages in the cluster diagram by a graded series of black circles. Untransformed abundance data are presented in Table 1. See Figures 4–7 for illustrations of genera. Sample localities are indicated by the following abbreviations: STR = Steptoe Ranch; CHC-1 and CHC-2 = Barton Canyon; KR1 = Ring Top Mountain; RR = Royer Ranch. Six biofacies defined in the Q-mode clustering are Camaraspis (Ca), Anechocephalus (An), Irvingella–Comanchia (Ir–Co), Parabolinoides (Pa), Mendoparabolina (Me), and Taenicephalus (Ta).

Table 1. Trilobite abundance data for collections used in the cluster analysis (Fig. 12). Abundances were calculated using the minimum number of individuals method (Gilinsky and Bennington, Reference Gilinsky and Bennington1994), which in practice meant the maximum number of either cranidia or pygidia for each genus in each collection.

The pre-extinction faunas (Figs. 4, 7.1–7.3, 7.9, 7.11) fall into two distinct clusters that represent collections from Nevada (Anechocephalus Biofacies) and Oklahoma (Camaraspis Biofacies). Although there are some shared genera (e.g., Elvinia; Figs. 4.8, 7.11), the dominant genera in Nevada (Fig. 4.9, 4.10) and Oklahoma (Fig. 7.1) do not co-occur, so that the differentiation between these deep- and shallow-subtidal settings is robust at the species level. Agnostid arthropods (e.g., Biciragnostus, Fig. 4.7; Homagnostus, Fig. 5.5, 5.6; Pseudagnostus, Figs. 5.10, 5.11, 6.10, 6.11) are persistent to locally abundant in biofacies throughout the study interval in the outer shelf of Nevada, but are rare in interior sites.

Collections from the I. “major” Zone from Nevada (Fig. 5) and Oklahoma (Fig. 7.4, 7.5, 7.8) form a single cluster (Irvingella–Comanchia Biofacies). This is consistent with Westrop and Cuggy's (Reference Westrop and Cuggy1999) analysis, which showed that biofacies based on genera became more widely distributed in intervals of faunal turnover (see also Ludvigsen and Westrop, Reference Ludvigsen and Westrop1983). The occurrence of the Irvingella “major” Zone fauna over a broad swathe of Laurentia North America has been known for almost 75 years (Wilson and Frederickson, Reference Wilson and Frederickson1950). Irvingella is in need of revision, but it is represented by similar species in Nevada (Fig. 5.4) and Oklahoma (Fig. 7.4). Comanchia, which is common in almost all collections assigned to the Irvingella–Comanchia Biofacies occurs as distinct species in Nevada (Fig. 5.1, 5.2) and Oklahoma (Fig. 7.5), as does Bartonaspis, a genus that is confined to the I. “major” Zone (Westrop and Adrain, Reference Westrop and Adrain2007). The data indicate some degree of faunal differentiation at the species level, but this may reflect a more general pattern. There is evidence that geographically and environmentally arrayed species groups may be a feature of trilobite faunas both between (e.g., Adrain and Westrop, Reference Adrain and Westrop2005, p. 379) and during (e.g., Westrop and Adrain, Reference Westrop and Adrain2007) Cambrian extinctions (see also Westrop et al., Reference Westrop, Landing and Dengler2018, for commentary on pseudocryptic trilobite species). At minimum, these data indicate that processes responsible for speciation prior to the onset of extinction remained active during turnover of clades and changes in their distribution in the extinction interval. It seems unlikely that the extinction is a consequence of a reduced rate of speciation.

The overlying Parabolinoides Subzone and correlatives continue the pattern of increased uniformity of shelf biofacies distribution, which was also documented by Westrop and Cuggy (Reference Westrop and Cuggy1999, figs. 14, 15; Orygmaspis Biofacies); see also the Orygmaspis (Parabolinoides) Biofacies of Chatterton and Gibb, Reference Chatterton and Gibb2016, in outer-shelf, deep-subtidal facies of the McKay Group of British Columbia. Biofacies in both the inner and outer shelf are dominated by olenids that broadly conform to current diagnoses of Orygmaspis (Parabolinoides) (e.g., Westrop, Reference Westrop1986a; Chatterton and Gibb, Reference Chatterton and Gibb2016), although this group of species is in need of revision. As in the older Irvingella–Comanchia Biofacies, there is a contrast in distributions at clade and species levels, with distinct species of O. (Parabolinoides) in Nevada (Fig. 6.8, 6.9) and Oklahoma (Fig. 7.6), and a stratigraphic succession of species in Alberta (Westrop, Reference Westrop1996, fig. 5). Although some reports still view olenids as diagnostic of low-oxygen environments (e.g., LeRoy et al., Reference LeRoy, Gill, Sperling, McKenzie and Park2021), the environmental distribution of taxa such as O. (Parabolinoides) across broad shelves into environments that are clearly well oxygenated shows the limitations of this approach.

The remaining biofacies are younger than the main interval of extinction and faunal replacement. The Taenicephalus Biofacies is represented in our data set by only one collection from Oklahoma, but it is typical of an association of genera (Taenicephalus, Fig. 7.7; Orygmaspis, Fig. 7.10) that is widely distributed in the upper part of the Taenicephalus Zone (Westrop and Cuggy, Reference Westrop and Cuggy1999, figs. 14, 15). The correlative interval in Nevada is unfossiliferous, but an outer-shelf, deep-subtidal Taenicephalus–Kendallina Biofacies has been reported from western Canada by Chatterton and Gibb (Reference Chatterton and Gibb2016). The systematics of Taenicephalus needs revision. New data from Oklahoma (S.R. Blackwell and Westrop, unpublished) indicate that one widely reported species, T. shumardi (Hall, Reference Hall1863), may actually represent three distinct species. However, it would not be surprising if Taenicephalus also proves to represent an array of geographically and environmentally segregated species.

The youngest biofacies in Nevada, the Mendoparabolina Biofacies is from a stratigraphic interval (Drumaspis Subzone of Stitt, Reference Stitt1971b, and correlatives) that was not sampled at section RR in Oklahoma. The Mendoparabolina Biofacies is dominated by olenids (Fig. 6.3, 6.4), with Loganellus (Fig. 6.1, 6.2) and Drumaspis (Fig. 6.6, 6.7) also important. The latter is part of a more inclusive, mostly Laurentian shelf clade (Elviniidae) that includes such genera as Irvingella and Elvinia that disappear during the extinction interval. Loganellus is typical of shelf-margin sites (e.g., Ludvigsen et al., Reference Ludvigsen, Westrop and Kindle1989) and is currently assigned to Family Idahoidae (e.g., Westrop, Reference Westrop1995). Idahoiids are well represented in Sunwaptan strata (e.g., Westrop, Reference Westrop1986a; Ludvigsen et al., Reference Ludvigsen, Westrop and Kindle1989), and sclerites of Noelaspis? n. sp. from the Barton Canyon Limestone Member (Fig. 4.3, 4.4) pulls the group down into the pre-extinction interval of the Steptoean Stage.

Diversity.—Westrop and Cuggy (Reference Westrop and Cuggy1999, figs. 16, 17) showed that species diversity declined through a series of biofacies replacements in the Sunwaptan, reaching a minimum in the Taenicephalus Zone. Because species-level systematics is still in progress, we have analyzed diversity in Nevada and Oklahoma at the genus level, but the results conform to Westrop and Cuggy's (Reference Westrop and Cuggy1999) conclusions. In both regions, rarefaction (Fig. 13) shows that genus diversity drops progressively from the Elvinia Zone into the Taenicephalus Zone. In post-Taenicephalus Zone collections from Nevada, there is a modest rebound that brings genus diversity back to levels comparable to those of the onset of faunal change in the I. “major” Zone. However, diversity remains below pre-extinction levels.

Figure 13. Rarefaction curves (calculated in PAST; Hammer et al., Reference Hammer, Harper and Ryan2001) for collections used in the cluster analysis. Curves for Nevada and Oklahoma show the same general pattern of declining numbers of genera from the Elvinia Zone into the Irvingella “major” Zone, with lowest numbers in the succeeding Taenicephalus Zone. In Nevada, collections from the post-Taenicephalus interval show a modest rebound to levels comparable to the I. “major” Zone.

Irvingella shell beds.—One of the characteristic features of the onset of faunal turnover and replacement is a spike in the abundance of Irvingella, which forms shell beds (e.g., Fig. 10.5) in a variety of environmental settings (e.g., Wilson and Frederickson, Reference Wilson and Frederickson1950; Westrop and Cuggy, Reference Westrop and Cuggy1999). Figure 14 shows the profound change in abundance of Irvingella in both shallow- and deep-subtidal settings in Oklahoma and Nevada. At all three sections, Irvingella is a minor component of the fauna prior to the onset of faunal change (lower set of bar charts in Fig. 14). Dominant genera of the Elvinia Zone faunas are lost during the I. “major” Zone (two upper rows of bar charts in Fig. 14), with the appearance of immigrants such as Comanchia (see Westrop and Cuggy, Reference Westrop and Cuggy1999, p. 349, for discussion), and all samples show a dramatic increase in the abundance of Irvingella. Estimating absolute abundances from shell beds is no simple matter. Condensation certainly played a role in onlap shell bed formation at the top of the Barton Canyon Limestone Member. Sediment starvation is implicated in the generation of shell accumulations in other settings (Dattilo et al., Reference Dattilo, Brett, Tsujita and Fairhurst2008) and likely influenced shell bed formation in the Honey Creek Formation. In addition, taphonomic sorting can lead to major changes in relative abundances (for a case involving Irvingella, see Westrop, Reference Westrop1986b, figs. 8, 9). However, widespread shell accumulations also imply a supply of raw materials and hence at least periodically substantial population sizes in the local environment.

Figure 14. Bar charts showing genus abundances (%) in the upper Elvina Zone and Irvingella “major” Zone in Nevada and Oklahoma. Samples for STR and CHC-1 are taken from the condensed interval at the top of the Barton Canyon Limestone Member (Fig. 9); layer 1 in each case is from the top of the Elvinia Zone and the overlying layers are from the I. “major” Zone. Samples from the Honey Creek Formation are from section KR1 (Figs. 10, 11) at Ring Top Mountain, Comanche County, Oklahoma (Westrop et al., Reference Westrop, Waskiewicz Poole and Adrain2010, fig. 1). The base of the I. “major” Zone is sample KR1 22. Note that Comanchia is present in the I. “major” Zone at CHC but Parabolinoides is absent. The reverse is true for samples from STR.

Irvingella is not the only trilobite taxon to become widespread and abundant in low-diversity biofacies during and immediately after the extinctions. Orygmaspis (Parabolinoides) becomes a dominant taxon across much of the shelf following the disappearance of Irvingella at the base of the Parabolinoides Zone, as does Taenicephalus higher in the succession (e.g., the Orygmaspis and Taenicephalus Biofacies of Westrop and Cuggy, Reference Westrop and Cuggy1999, figs. 14, 15). Low diversity and numerical dominance of biofacies by a single genus is a persistent feature of the extinctions.

Orthid brachiopod shell beds.—Orthid brachiopods appear in the I. “major”–lower Taenicephalus zones in Oklahoma (see Freeman and Stitt, Reference Freeman and Stitt1996, for treatment of the systematics and biostratigraphy) and, like Irvingella, form shell beds (Fig. 10.2, 10.4, 10.6). Orthids are abundant only through a few meters of strata (e.g., Fig. 11), and then become rare. The brachiopod “blooms” are a shallow-water signature of the extinction interval, and orthids are entirely absent from deeper facies in Nevada. Orthids are, however, present in the extinction interval in the mixed carbonate–siliciclastic, storm-influenced succession of the Bison Creek Formation in Alberta (Westrop, Reference Westrop1989), where they occur in the Taenicephalus Zone (Fig 15.2). In similar facies of the Snowy Range Formation of Montana and Wyoming, brachiopods also appear in the Taenicephalus Zone, where species of Eoorthis, Billingsella, Huenella, and Otusia reach their epiboles and form shell beds in association with O. (Parabolinoides) (Grant, Reference Grant1965, p. 85, 90). Saltzman (Reference Saltzman1999, p. 929, 931, fig. 4B) also documented orthid shell beds (coquinite) from the same stratigraphic interval in Wyoming. In inner-shelf sandstone of the Lone Rock Formation, Eoorthis is abundant enough in the lower part of the Taenicephalus Zone to give its name to the Eoorthis Subzone (Berg, Reference Berg1953). On other Cambrian continents, the Orusia shell beds of Avalonia and Baltica appear at a similar stratigraphic level to the orthid beds of Laurentia (Landing and Westrop, Reference Landing and Westrop2015, p. 982).

Although absent from the extinction interval at the base of the Steptoean (Paibian) Stage, brachiopod shell beds are a feature of younger faunal changes in the Sunwaptan–Skullrockian stage boundary interval in the lower Survey Peak Formation of Alberta (Westrop, Reference Westrop1986a; Fig. 15.1) and the San Saba Member of the Wilberns Formation (Winston and Nicholls, Reference Winston and Nicholls1967, p. 92), in the Skullrockian–Stairsian boundary interval in shallow-water carbonates in the House Range of Utah (Adrain et al., Reference Adrain, Westrop, Karim and Landing2014, appendix 2, p. 206; Saltzman et al., Reference Saltzman, Edwards, Adrain and Westrop2015, fig. 3), and are bracketed by trilobite collections marking the pre- and post-extinction faunas at the base of the Tulean Stage in a subtidal, storm-influenced succession in the Garden City Formation of southeastern Idaho (Adrain et al., Reference Adrain, McAdams and Westrop2009, fig. 5).

Figure 15. Distribution of trilobite genera, trilobite shell beds and orthid brachiopod shell beds in extinction intervals at the bases of the Sunwaptan and Skullrockian stages at Wilcox Peak, Alberta (data from Westrop, Reference Westrop1984). (1) Lower Sunwaptan strata of the Bison Creek Formation. (2) Lower Skullrockian strata of the basal silty member, Survey Peak Formation.

These various brachiopod blooms were relatively short-lived episodes and at section RR, they correspond closely with turnover of trilobite genera and peak values of δ¹³C_carb (Fig. 16). Much like Irvingella, orthid brachiopods have some of the properties of opportunistic “disaster taxa” (Schubert and Bottjer, Reference Schubert and Bottjer1995; Petsios and Bottjer, Reference Petsios and Bottjer2016). This label does not singularly indicate causal factors for invasion of the shallow shelf by orthids. Using Lingula as an example, Petsios and Bottjer (Reference Petsios and Bottjer2016) suggested that disaster taxa invaded vacated ecospace briefly in the wake of extinction events, only to become competitively restricted in distribution as faunas recovered. However, it is far from obvious that this is the case for orthid brachiopods, which seem to be entering niche space that was essentially unoccupied prior to the extinction, and was not re-filled after the period of their high abundance and wide distribution. Orthids are generally rare between trilobite extinctions (e.g., see occurrence data in Freeman and Stitt, Reference Freeman and Stitt1996), but there is no other group filling the role of large, robustly calcified, sessile, epifaunal suspension feeder. Pelmatozoan debris is a major contributor to the grainstone and rudstone of the Honey Creek Formation (see also Stitt, Reference Stitt1971b; Donovan, Reference Donovan and Johnson2000) prior to, during, and after the extinction interval, indicating that invasion by orthids expanded the suspension-feeder niche, rather than simply entering vacated ecospace.

Figure 16. Data from the Royer Ranch section (RR), Murray County, Oklahoma (Stitt, Reference Stitt1971b, p. 64–66). Carbon isotope curve gives raw values without any form of moving-average smoothing. Black rectangles show occurrences and ranges of trilobite and agnostid genera from the collections listed by Stitt; white rectangles show occurrences and ranges of orthid brachiopod genera from Freeman and Stitt (Reference Freeman and Stitt1996). Stitt's (Reference Stitt1971b) measurements in feet were converted into meters. The gray-shaded band shows I. “major” Zone and Parabolinoides Subzone of the Taenicephalus Zone. The trajectory of the carbon isotope curve is more complex than at section CHC (Fig. 17), with a decline into the I. “major” Zone from a peak in the upper Elvinia Zone, followed by a steady rise of ~1‰ in the lower part of the Taenicephalus Zone. However, the upper Elvinia peak is pulled by a single extreme value, and if this outlier is ignored, the data arguably form a single rising trend from the pre-extinction to post-extinction intervals.

Diversity of trilobites does decline as brachiopods expand in shallow water, but trilobites remain abundant and occur as primary constituents of shell beds. At section KR1, trilobite- and brachiopod-dominated shell beds are interbedded (Fig. 11), which presumably reflects interplay of ecological processes, such as the vagaries of larval recruitment, and depositional processes. Trilobites are seemingly less abundant in brachiopod-rich beds (e.g., Fig. 10.6), but this is almost certainly a taphonomic signature, with thinner-shelled trilobites typically fragmented in accumulations of more-robust brachiopod valves.

It is difficult to imagine direct biological interactions between vagile and sessile benthos like trilobites and brachiopods that might explain the blooms. Rather, it seems more likely that they record the response to changes in the physical environment. Some modern species are tolerant of low-oxygen conditions and can reach high populations densities at depths of 50–700 meters in the fiords of coastal British Columbia (Tunnicliffe and Wilson, Reference Tunnicliffe and Wilson1988). Curry et al. (Reference Curry, Ansell, James and Peck1989) noted that metabolic rates of modern brachiopods are low, and reported that oxygen consumption is in the range of 10–50% of gastropods and bivalves housed under the same laboratory conditions. At first glance, these observations may seem to support hypotheses based on upwelling of dysoxic waters, but the environmental distribution of orthids in the extinction interval is problematic. They reach peak abundances in shallow water, with evidence in Oklahoma of tidal currents (e.g., Bucheit and Donovan, Reference Bucheit, Donovan and Johnson2000), suggesting well-mixed waters that were unlikely to be dysoxic. Moreover, in regions such as Nevada, where there is some sedimentary evidence to suggest a shift towards less-oxygenated environments, orthid brachiopods are absent.

In short, orthids are invasive taxa in the extinction interval, but the underlying causal factors are obscure. We also note that the diversification of rhynchonelliform brachiopods during the Ordovician Radiation records an expansion that, unlike the blooms of the late Cambrian and Early Ordovician, was sustained. An understanding of the ecological and evolutionary dynamics of Cambrian and Early Ordovician rhynchonelliforms may also throw light on the group's subsequent success.

Stratigraphic condensation influences the apparent tempo of faunal change.—Building on earlier work (e.g., Holland, Reference Holland2000), Holland and Patzkowsky (Reference Holland and Patzkowsky2015) explored the effect of stratigraphic biases on the record of mass extinctions. Their simulations showed that last occurrences of taxa will cluster at a variety of significant stratal surfaces, including flooding surfaces and forced regression surfaces, as well as various points of condensation within sequence architecture. Moreover, clusters of apparent extinction at a stratal surface may actually predate the time of peak extinction. Holland and Patzkowsky (Reference Holland and Patzkowsky2015) suggested that the record of Cambrian extinctions, which may be associated with flooding surfaces, is affected by such biases.

Faunal turnover within the condensed onlap shell bed at the top of the Barton Canyon Limestone Member (Fig. 9) indicates that our perception of the tempo of faunal change is indeed influenced by the sequence stratigraphic context. The data from Nevada and Oklahoma show that the onset of extinction and faunal change can be narrowed to stratal intervals measured in centimeters. Holland and Patzkowsky's (Reference Holland and Patzkowsky2015) simulations show that we cannot take this apparent abruptness of turnover at face value. At section RR (Fig. 16), the Irvingella “major” fauna was sampled at a single horizon, and bracketing by other collections shows that the zone can be no more than 0.75 m (2.5 feet) in thickness (Stitt, Reference Stitt1971b). As noted earlier, we have not been able to gain access to this or other classic localities in the Arbuckle Mountains region (e.g., Frederickson, Reference Frederickson1949; Stitt, Reference Stitt1971b) to make new observations, but condensation seems likely. Our thickest succession through the Irvingella “major” Zone at section KR1 expands to 2 m, and the base, which is not associated with a significant stratal surface, is localized to within 10 cm of a collection with a pre-extinction fauna (Fig. 11). At the same time, formation of shell beds implies some degree of sediment starvation (or an increase in shell production), so we cannot rule out some degree of condensation.

The evidence from Nevada and Oklahoma indicates that faunal replacement occurs across biostratigraphic boundaries that may correspond to significant stratal surfaces, but successive biofacies that occur between those boundaries (e.g., Westrop and Cuggy, Reference Westrop and Cuggy1999) are the products of the extinctions—assemblages that have been pruned by turnover but also augmented by immigration. Ironically, they are time-averaged records of stable paleocommunities (at the genus level, with replacement at the species level; Westrop, Reference Westrop1996) within a broader interval of taxonomic turnover. They represent the incumbents at different stages in what is a drawn-out process of faunal change. The actual extinctions will be difficult to resolve because they involved ecological processes that, as Schindel (Reference Schindel1980, table 2) recognized, occur too quickly to be identified in the sedimentary record, and that record may be further distorted by stratigraphic biases (Holland and Patzkowsky, Reference Holland and Patzkowsky2015).

Internally, the layers representing the Irvingella “major” Zone in the condensed shell bed at the top of the Barton Canyon Limestone Member are uniformly dominated by Irvingella, but there is variability in the associated taxa (Fig. 14). In Nevada, Anechocephalus and Elvinia, dominant genera in the pre-extinction fauna, are present at very low abundances in the basal layer (layer 3) of the I. major Zone at section CHC. At CHC, Comanchia is present in both sample intervals, and ranked second in abundance in the lower of the two, whereas O. (Parabolinoides) is absent. In contrast, O. (Parabolinoides) is present in both layers at STR, but Comanchia is absent. This could record ecological differences in coeval assemblages. Alternatively, in these time-rich beds, time is not sampled evenly between sites. Rather, individual layers at different sites provide snapshots of different time intervals. Pieced together, they may allow taxonomic turnover to be reconstructed in more detail.

A common pattern over much of Laurentian North America is the replacement of Irvingella-dominated biofacies by O. (Parabolinoides)-dominated biofacies (Westrop and Cuggy, Reference Westrop and Cuggy1999). Our data show that Irvingella and O. (Parabolinoides) coexisted for part of the Irvingella “major” Zone (see also Stitt, Reference Stitt1971b, p. 9, 57), with O. (Parabolinoides) becoming abundant and widespread following the disappearance of Irvingella. This replays the pattern noted above for Irvingella, which was rare in pre-extinction assemblages, but expanded shelf-wide after the first episode of turnover. This suggests that incumbency is important in maintaining biofacies composition during the interludes between episodes of turnover, perhaps mediated by processes such as competition, which is at least plausible for interactions between trilobite species.

Durations of the extinctions are difficult to estimate, but even the relatively expanded section through the I. “major” Zone at section KR1 suggests that they occurred over a geologically brief period. Minimum estimates from calibration of the SPICE δ¹³C excursion with dates from detrital zircons (Cothren et al., Reference Cothren, Farrell, Sundberg, Dehler and Schmitz2022) indicate that the Paibian Stage may have a duration of as little as 1.6 Myr (494.4–492.8 Ma). These dates bracket about 70 m of strata in northern Utah (Cothren et al., Reference Cothren, Farrell, Sundberg, Dehler and Schmitz2022). Leaping from dated tie points to estimates of durations of smaller packages of strata is difficult if not dangerous, but these new dates suggest that the extinctions occurred over a time scale that is likely measured in tens of thousands of years. Or, put another way, the period of relative stability recorded by the multiple shell beds of the Irvingella–Comanchia Biofacies at KR1 could easily be on par with the period of deglaciation and sea-level rise from the last glacial maximum to the present (e.g., Anderson and McBride, Reference Anderson and McBride1996).

Discussion.—During the extinctions at the base of the Sunwaptan Stage, the epicontinental seas of Laurentian North America were extraordinary from an ecological perspective. At various points in a succession of faunal replacements, a single trilobite genus is numerically dominant in all shelf environments for which samples are available (e.g., Westrop and Cuggy, Reference Westrop and Cuggy1999). Trilobite biofacies are a mixture of holdovers from underlying faunas and immigrants. Immigration and high abundance are also dominant themes in the appearance of orthid brachiopods in shallow-water sites over the shelf. In the late Cambrian and Early Ordovician, trilobite–orthid-dominated paleocommunities are known only in extinction intervals. The appearance of this paleocommunity type is one of the dominant biological signals of these events, and yet the causal factors behind the immigration and coordinated blooms of these taxa remain unknown.

CHC composite section

Measurements of δ¹³C came from samples collected from two sections measured through a 30-meter-thick interval comprising the Barton Canyon Limestone and lower Catlin members of the Windfall Formation on either side of Barton Canyon (Fig. 3; Table 2). These were compiled into a composite curve by aligning the sections using the top of the Barton Canyon Limestone Member as a tie point, which also allowed us to plot a composite range chart for trilobite genera (Fig. 17). Although we did not screen the samples petrographically for diagenetic alteration, a cross-plot of δ¹³C and δ¹⁸O (Fig. 18.1) did not reveal a significant correlation (Pearson's correlation coefficient [r]: 0.25, p = 0.23, r ² = 0.061), so there is no evidence for pervasive meteoric diagenetic alteration, albeit from just a single geochemical indicator.

Table 2. Stable isotope data from the composite section at Barton Canyon, Nevada (section CHC).

Figure 17. Composite section for CHC, combining data from sections CHC-1 and CHC-2 (Fig. 3), which were aligned using the top of the I. “major” Zone at the top of the Barton Canyon Limestone Member. Carbon isotope curve gives raw values without any form of moving-average smoothing. Black rectangles show occurrences and ranges of trilobite and agnostid genera. The gray-shaded band shows I. “major” Zone and Parabolinoides Subzone of the Taenicephalus Zone. In this interval, the carbon isotope curve shows a modest rise to a little more than 2‰.

Figure 18. Cross-plots of carbon and oxygen for sections CHC composite and RR. (1) CHC composite cross-plot showing no significant correlation between carbon isotope and oxygen isotope values; Pearson's correlation coefficient (r) = 0.25, p = 0.23, r ² = 0.061. (2) Section RR cross-plot, again showing no significant correlation between carbon isotope and oxygen isotope values; Pearson's correlation coefficient (r) = −0.17, p = 0.318; r ² = 0.03.

In the pre-extinction interval that comprises almost all of the Barton Canyon Limestone Member, there is no net trend in δ¹³C (Fig. 17), which oscillates around a mean value of 0.85‰ (median = 0.87‰) with a standard deviation of 0.32 (range = 0.12–1.25). Faunal replacement occurs across the lower and upper boundaries of the Irvingella “major” Zone in the condensed shell beds at the top of the Barton Canyon Limestone Member, with a sharp drop in genus diversity in collections immediately above this zone. Taking this succession of turnover and diversity decline as the primary extinction interval, the changes occur at the drowning unconformity in association with a trend towards higher values of δ¹³C that peak slightly above 2‰ (Fig. 17). Mean and median values in this interval are 1.44‰, with a standard deviation of 0.43 (range = 0.63–2.17). Mean values for pre-extinction and extinction intervals are significantly different (t-test; t = 3.579, p = 0.002). There are only three samples available above this interval, and δ¹³C drops, with a median value of 0.34‰.

The overall pattern is one of a modest positive excursion that begins near the onset of extinctions, continues beyond, and then peaks above the second phase of turnover (e.g., Westrop and Cuggy, Reference Westrop and Cuggy1999) at the top of the I. “major” Zone. This is congruent with the documented record of δ¹³C through this interval in Wyoming (Saltzman et al., Reference Saltzman, Davidson, Holden, Runnegar and Lohmann1995; Saltzman, Reference Saltzman1999) and at sites farther to the east of CHC in Nevada and Utah (Saltzman et al., Reference Saltzman, Runnegar and Lohmann1998).

Section RR

The Royer Ranch section, Murray County (Stitt, Reference Stitt1971b, p. 64) was selected for geochemical analysis because collections are closely spaced through the Steptoean–Sunwaptan boundary interval, and they yield abundant trilobites that are currently being revised by Westrop. The sample interval extends from the upper Elvinia Zone into the “Idahoia” lirae Zone. The samples are skeletal packstone, grainstone, and rudstone, and while micrite was sampled wherever possible, this could not be done with grainstone and rudstone lithologies. A cross-plot of δ¹³C and δ¹⁸O (Fig. 18.2) did not reveal a significant correlation (Pearson's correlation coefficient [r]: −0.17, p = 0.318; r ² = 0.03).

The δ¹³C curve (Fig. 16; Table 3) oscillates through much of the pre-extinction interval (Elvinia Zone), although there is a peak reaching 1.16‰ that is primarily the product of two data points located 1.2–0.91 m below the base of the Irvingella “major” Zone. The mean for the pre-extinction interval is 0.31‰ (median = 0.28‰), with a standard deviation of 0.37 and a range of −0.19 to 1.16‰. Values decline into the I. “major” Zone, reaching a minimum at the base of the Parabolinoides Subzone, and then rising steadily to a peak at 1.36‰ through 4.25 m of section into the Taenicephalus “shumardi” Subzone. This second rising trend has a mean value of 0.96‰ (median = 1.08‰), with a standard deviation of 0.35; the mean is significantly different from the pre-extinction values (t-test; t = 4.044; p = 0.0006). The second positive trend encompasses the diversity minimum of genera (Parabolinoides Subzone), and the “bloom” of orthid brachiopods.

Table 3. Stable isotope data from the section at Royer Ranch, Oklahoma (Stitt, Reference Stitt1971b; section RR).

The pattern of change in δ¹³C through the upper Elvinia Zone and into the Taenicephalus Zone could be partitioned into two trends of increasingly positive values separated by a reversal at the onset of extinctions. However, this is driven in part by the extreme value (1.16‰ at 42.37 m) in the upper Elvinia Zone. If this outlier is excluded, one could make a case for a single overall trend towards increasingly positive values, with a more “noisy” pattern in the Elvinia and I. “major” zones. The significant difference between mean values of the two trends is consistent with this alternative interpretation. A similar, essentially steady rise in δ¹³C from the Elvinia Zone into the Taenicephalus Zone was recorded by Saltzman (Reference Saltzman1999, fig. 8) in the Open Door Formation and correlative Snowy Range Formation in northwest Wyoming. Like section RR, the peak values of δ¹³C in Wyoming are lower (reaching 1.11‰ about 2 m above the base of the Taenicephalus Zone at Warm Springs Peak; Saltzman, Reference Saltzman1999, table 1) than maxima recorded in Nevada and Utah (Saltzman et al., Reference Saltzman, Davidson, Holden, Runnegar and Lohmann1995, Reference Saltzman, Runnegar and Lohmann1998).

Correlation and interpretation of the δ¹³C curves

Correlation of the curves for CHC, RR, and Saltzman et al.'s (Reference Saltzman, Runnegar and Lohmann1998) section at Little Horse Canyon, Orr Ridge, Utah, is shown in Figure 19. Given the nature of our data, we focus on general trends reproduced in all three curves rather than magnitudes of individual values. A shared feature of the curves is a trend towards increasingly positive δ¹³C values that begins in the vicinity the I. “major” Zone (either just below the base, as at Orr Ridge, or immediately above, as at RR). As noted above, a similar rise starting in the Elvinia Zone is apparent in the Steptoean–Sunwaptan boundary interval in Wyoming (Saltzman, Reference Saltzman1999). Below the Irvingella “major” Zone, there is no clear trend at CHC and Orr Ridge, whereas at RR, there is a “noisy” positive trend that could be interpreted as part of a longer-term rise into the Taenicephalus Zone.

Figure 19. Correlation of carbon isotope curves for sections RR (Oklahoma), CHC composite (Nevada), and Orr Ridge, Utah (plotted from data in Saltzman et al., Reference Saltzman, Runnegar and Lohmann1998, GSA Data Repository item 9804), using the base of the Irvingella “major” Zone as a datum; dashed lines indicate gaps in sampling. The overall pattern of a rising trajectory through the extinction interval can be identified at all three localities, and is consistent with the pattern reported from Nevada and Wyoming by Saltzman et al. (Reference Saltzman, Davidson, Holden, Runnegar and Lohmann1995, fig. 3; see also Saltzman, Reference Saltzman1999).

As discussed by Saltzman (Reference Saltzman1999), the trend towards increasingly more positive δ¹³C values across the top of the Steptoean reprises, at a shorter temporal scale, the early part of the SPICE excursion (Saltzman et al., Reference Saltzman, Runnegar and Lohmann1998) at the base of the stage. Saltzman et al. (Reference Saltzman, Davidson, Holden, Runnegar and Lohmann1995) interpreted the positive trend across the Steptoean–Sunwaptan boundary as the result of increased carbon burial due to upwelling of anoxic water, possibly in response to a sea-level rise. Gerhardt and Gill (Reference Gerhardt and Gill2016) offered a similar interpretation of the early part of the SPICE event near the base of the Steptoean. Landing (Reference Landing2012) provided stratigraphical and sedimentological evidence for an expanded oxygen minimum zone through much of the upper Cambrian.

Anoxia is supported by multiple lines of geochemical evidence. Gill et al. (Reference Gill, Lyons, Young, Kump, Knoll and Saltzman2011) documented a positive δ³⁴S_cas (carbonate-associated sulphate) excursion that paralleled the δ¹³C excursion of the SPICE in Laurentian North America and Australia, as well as a similar excursion in pyrite (δ³⁴S_pyr) in Sweden (Baltica). Gill et al. (Reference Gill, Lyons, Young, Kump, Knoll and Saltzman2011) interpreted these coupled excursions as the product of increased burial of organic carbon and pyritic sulfur under anoxic and euxinic conditions, and proposed that the extinction was the result of anoxic waters moving onto the shelf during a sea-level rise. Dahl et al. (Reference Dahl, Boyle, Canfield, Connelly, Gill, Lenton and Bizzarro2014) supported expansion of euxinic waters in the early part of the SPICE from a negative excursion of δ²³⁸U in the Mt. Whelan core, Australia, although these conditions were not sustained in the latter phases of the SPICE.

Changes in δ¹³C across the Steptoean–Sunwaptan boundary are consistent with upwelling of dysoxic waters, although the presence of trilobites in the Catlin Member indicates that oxygenation must have reached the threshold for respiration at least periodically. However, it remains to be seen exactly how anoxia may have caused these extinctions or, for that matter, the extinction at the onset of the SPICE at the base of the Steptoean Stage. There is independent evidence for sea-level rise and environmental change at section CHC, where the deep-water facies of the Catlin Member appear above the shallower carbonates of the Barton Canyon Limestone Member. However, anoxia is unlikely to have spread into shallow, wind- and tidally mixed waters of the shelf interior. Although extinction and immigration are associated with shifts in δ¹³C in the Honey Creek Formation, there is no independent sedimentary evidence for anoxia. This suggests a role for ecological and biogeographical effects (e.g., reduced population sizes and geographic ranges; interactions with invasive species) during the extinction in response to anoxic and dysoxic waters impinging on habitats in the outer part of the shelf (e.g., Westrop and Ludvigsen, Reference Westrop and Ludvigsen1987; Saltzman et al., Reference Saltzman, Edwards, Adrain and Westrop2015). Geographic restriction to isolated “relicts” and diachronous disappearances of taxa, which occurs at the end-Sunwaptan extinction, are consistent with ecological and biogeographic processes (Landing et al., Reference Landing, Westrop, Kröger and English2010, p. 541).