Hostname: page-component-788cddb947-nxk7g Total loading time: 0 Render date: 2024-10-09T17:30:33.295Z Has data issue: false hasContentIssue false

Coupling of geographic range and provincialism in Cambrian marine invertebrates

Published online by Cambridge University Press:  09 December 2022

Lin Na
Affiliation:
State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology, Chinese Academy of Sciences, Beijing East Road 39, 210008 Nanjing, China. E-mail: linna@nigpas.ac.cn
Ádám T. Kocsis
Affiliation:
GeoZentrum Nordbayern, Department of Geography and Geosciences, University of Erlangen-Nuremberg, Loewenichstraße 28, 91054 Erlangen, Germany. E-mail: adam.kocsis@fau.de, wolfgang.kiessling@fau.de
Qijian Li*
Affiliation:
State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology, Chinese Academy of Sciences, Beijing East Road 39, 210008 Nanjing, China; and GeoZentrum Nordbayern, Department of Geography and Geosciences, University of Erlangen-Nuremberg, Loewenichstraße 28, 91054 Erlangen, Germany. E-mail: qjli@nigpas.ac.cn
Wolfgang Kiessling
Affiliation:
GeoZentrum Nordbayern, Department of Geography and Geosciences, University of Erlangen-Nuremberg, Loewenichstraße 28, 91054 Erlangen, Germany. E-mail: adam.kocsis@fau.de, wolfgang.kiessling@fau.de
*
*Corresponding author.

Abstract

The Cambrian saw a dramatic increase in metazoan diversity and abundance. Between-assemblage diversity (beta diversity) soared in the first three Cambrian stages, suggesting a rapid increase in the geodisparity of marine animals during the Cambrian radiation. However, it remains unclear how these changes scale up to first-order biogeographic patterns. Here we outline time-traceable provinces for marine invertebrates across the Cambrian period using a compositional network based on species-level fossil occurrence data. Results confirm an increase in regional differences of faunal composition and a decrease in by-species geographic distribution during the first three stages. We also show that general biogeography tends to be reshaped after global extinction pulses. We suggest that the abrupt biogeographic differentiation during the Cambrian radiation was controlled by a combination of tectonics, paleoclimate, and dispersal capacity changes.

Type
Articles
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), 2022. Published by Cambridge University Press on behalf of The Paleontological Society

Introduction

The Cambrian witnessed a noteworthy radiation in biodiversity (Sepkoski Reference Sepkoski1997; Maloof et al. Reference Maloof, Porter, Moore, Dudas, Bowring, Higgins, Fike and Eddy2010; Na and Kiessling Reference Na and Kiessling2015) and a remarkable divergence in Bauplans (Erwin et al. Reference Erwin, Laflamme, Tweedt, Sperling, Pisani and Peterson2011) of marine metazoans. Although the temporal and spatial patterns of Cambrian biodiversity have been widely documented (e.g., Zhuravlev and Riding Reference Zhuravlev and Riding2001; Jago et al. Reference Jago, Zang, Sun, Brock, Paterson and Skovsted2006; Li et al. Reference Li, Steiner, Zhu, Yang, Wang and Erdtmann2007; Peng et al. Reference Peng, Babcock, Cooper, Ogg, Ogg and Gradstein2012; Na and Kiessling Reference Na and Kiessling2015; Kröger et al. Reference Kröger, Franeck and Rasmussen2019; Rasmussen et al. Reference Rasmussen, Kröger, Nielsen and Colmenar2019; Fan et al. Reference Fan, Shen, Erwin, Sadler, MacLeod, Cheng, Hou, Yang, Wang, Wang, Zhang, Chen, Li, Zhang, Shi, Yuan, Chen, Zhang, Li and Zhao2020), less attention has been paid to the fundamental biogeographic structure and its changes through time. Biogeographic patterns may provide insights into evolutionary processes during large-scale diversifications at different temporal scales (Valentine et al. Reference Valentine, Foin and Peart1978; Fortey et al. Reference Fortey, Briggs and Wills1996; Crame and Owen Reference Crame and Owen2002; Huang et al. Reference Huang, Jin and Rong2018).

Cambrian provinciality has been recorded among various fossil groups at different taxonomic levels (Theokritoff Reference Theokritoff1979; Lu Reference Lu1981; Zhuravlev Reference Zhuravlev1986; Steiner et al. Reference Steiner, Li, Qian, Zhu and Erdtmann2007; Yang et al. Reference Yang, Steiner and Keupp2015). A two-part or three-part provincial scheme has been suggested to represent early Cambrian trilobite biogeography (Lu et al. Reference Lu, Zhu, Qian, Lin, Zhou and Yuan1974; Palmer and Repina Reference Palmer and Repina1993; Fortey et al. Reference Fortey, Briggs and Wills1996): the Redlichiid Province in Gondwana; the Olenellid Province in Laurentia, Baltica, and Siberia; and sometimes an intermediate province representing overlap domains (Lu Reference Lu1981; Pillola Reference Pillola1989). The identification of trilobite biogeography relied largely on distinctive endemic trilobites among major continents and has proved unreliable in light of phylogenetic analyses (Lieberman Reference Lieberman1998; Meert and Lieberman Reference Meert and Lieberman2004). Taking phylogenetic history of trilobites into account, more complex biogeographic schemes among trilobites were advocated for the early and the late Cambrian (Meert and Lieberman Reference Meert and Lieberman2004; Álvaro et al. Reference Álvaro, Ahlberg, Babcock, Bordonaro, Choi, Cooper, Ergaliev, Gapp, Pour and Hughes2013), and many attempts have been made to compare phylogenetic units among tectonic blocks to reconstruct plate tectonic configurations (Jell Reference Jell1974; Lieberman Reference Lieberman1997, Reference Lieberman2005; Álvaro et al. Reference Álvaro, Ahlberg, Babcock, Bordonaro, Choi, Cooper, Ergaliev, Gapp, Pour and Hughes2013; Steiner et al. Reference Steiner, Hohl, Yang, Huang and Li2021). However, traditional approaches based on taxon- and time slice–specific analysis (e.g., Theokritoff Reference Theokritoff1979; Yang et al. Reference Yang, Steiner and Keupp2015) can hardly provide a panoramic perspective on secular trends in biogeography and their link to diversification patterns (Kocsis et al. Reference Kocsis, Reddin and Kiessling2018a,b). Moreover, the heterogenous fossil coverage distorts our understanding of deep-time biogeographic patterns (Alroy et al. Reference Alroy, Aberhan, Bottjer, Foote, Fürsich, Harries, Hendy, Holland, Ivany, Kiessling, Kosnik, Marshall, McGowan, Miller, Olszewski, Patzkowsky, Peters, Villier, Wagner, Bonuso, Borkow, Brenneis, Clapham, Fall, Ferguson, Hanson, Krug, Layou, Leckey, Nürnberg, Powers, Sessa, Simpson, Tomašových and Visaggi2008; Close et al. Reference Close, Benson, Saupe, Clapham and Butler2020), which has necessitated multi-taxon approaches and a more reproducible quantitative framework in delineating biogeographic units (Kocsis et al. Reference Kocsis, Reddin, Scotese, Valdes and Kiessling2021).

Here, we derive the first-order biogeographic patterns of macroinvertebrates across the entire Cambrian period. We apply a network analysis of species occurrence data while accounting for uneven sampling. The network-based approach has been widely applied to fossils to reconstruct biogeographic patterns in deep time (Sidor et al. Reference Sidor, Vilhena, Angielczyk, Huttenlocker, Nesbitt, Peecook, Steyer, Smith and Tsuji2013; Vilhena et al. Reference Vilhena, Harris, Bergstrom, Maliska, Ward, Sidor, Strömberg and Wilson2013; Rojas et al. Reference Rojas, Patarroyo, Mao, Bengtson and Kowalewski2017; Dunne et al. Reference Dunne, Close, Button, Brocklehurst, Cashmore, Lloyd and Butler2018; Huang et al. Reference Huang, Jin and Rong2018; Rong et al. Reference Rong, Harper, Huang, Li, Zhang and Chen2020), but no study has focused on the Cambrian period. Our methodology traces community structures and outlines time-traceable biogeographic units (bioregions; Kiel Reference Kiel2017; Kocsis et al. Reference Kocsis, Reddin and Kiessling2018a,Reference Kocsis, Reddin and Kiesslingb). We employ a multi-fauna approach, but for comparison with earlier works, we also apply our method to only trilobite occurrences. By assessing the relationship between biogeographic changes and species’ geographic ranges, we further explore the timing of provincial divergence/convergence and derive the potential drivers of biogeographic changes during the Cambrian radiation.

Material and Methods

Data

We downloaded occurrence data of Cambrian marine invertebrate species from the Paleobiology Database (https://paleobiodb.org) on September 3, 2021. The raw data comprised 9214 collections and 41,912 fossil occurrences. For temporal binning, we updated the framework for Cambrian biostratigraphic correlation following Na and Kiessling (Reference Na and Kiessling2015) (Supplementary File F1), with time binning created with the divDyn package (Kocsis et al. Reference Kocsis, Reddin, Alroy and Kiessling2019) in R (R Core Team 2021). Ages refer to the Geologic Time Scale 2020 (Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2020). Only collections that could be assigned to a single stage were retained. We also removed all occurrences that could not be assigned to a species with confidence as well as form species and ichnospecies. Species occurrences with the qualifiers “cf.”, “?”, “aff.”, “n. sp.”, “n. gen.”, or “indet.” were disregarded. The final dataset comprised 4397 collections and 15,622 species occurrences of 2644 species (Supplementary Table S1; also see Supplementary Material available through Zenodo). Nearly half of the original collections were dropped due to the uncertainty in stratigraphic binning.

Paleo-coordinates were obtained by rotating the present-day coordinates of fossil collections with the GPlates v. 2.2.0 software (Müller et al. Reference Müller, Cannon, Qin, Watson, Gurnis, Williams, Pfaffelmoser, Seton, Russell and Zahirovic2018). We used the rotation file supplied by C. R. Scotese in the PALEOMAP PaleoAtlas for GPlates package (Scotese Reference Scotese2016). Accurate Cambrian reconstructions of the configuration of Earth's continents come with substantial spatial and temporal uncertainties (Landing et al. Reference Landing, Westrop and Bowring2013). With a different paleogeographic model that reconstructs Cambrian geography differently, our results might be somewhat different. We rotated all collections to their positions at 510 Ma (the middle of the Cambrian) to (1) make the variation among cells comparable through time without largely influenced by tectonic variations, (2) better determine how provincial assignments evolved in a particular region, and (3) better relate the source data to our paleobiogeographic reconstructions.

We aggregated the collections based on their paleo-coordinates to equal-area hexagonal grid cells using a tessellated icosahedral geographic grid from the icosa package in R (Kocsis Reference Kocsis2017). To decrease the mismatch between the areas covered by a geographic cell in two subsequent time slices due to plate tectonic movements while retaining within-continent variation as much as we could, we used a rather coarse spatial resolution that gives roughly 252 cells per time interval with one grid cell being roughly 2 million km2 in area (edge length = 8°). In this way, 2520 spatiotemporal cells (10 geological stages × 252 geographic cells) were outlined and used in the network analysis to construct a time-uniform biogeographic partitioning scheme. Following a previous analysis (Kocsis et al. Reference Kocsis, Reddin and Kiessling2018b), we employed a minimum quota of at least 10 occurrences in each spatiotemporal unit to be included in the bipartite network.

Network Analysis

The collections in a time interval were assigned to equal-area geographic grid cells, and species were tabulated in every spatiotemporal cell. From that matrix, a bipartite occurrence network (Vilhena and Antonelli Reference Vilhena and Antonelli2015) was constructed based on shared species between the 2520 spatiotemporal units, wherein every spatiotemporal cell was represented with a vertex, and units were connected if the spatiotemporal cells shared species. The strength of connections was determined by the number of shared species and sampling. We performed a weighted projection from the bipartite network onto the node subset (Rojas et al. Reference Rojas, Patarroyo, Mao, Bengtson and Kowalewski2017) and applied the Infomap community detection algorithm (Rosvall and Bergstrom Reference Rosvall and Bergstrom2008) to partition the spatiotemporal grid cells within the projected network (Edler et al. Reference Edler, Guedes, Zizka, Rosvall and Antonelli2016; Rojas et al. Reference Rojas, Patarroyo, Mao, Bengtson and Kowalewski2017). The defined modules (i.e., bioregions) were then projected into geographic space.

To investigate the biogeographic dynamics, we used two indices (Kocsis et al. Reference Kocsis, Reddin and Kiessling2018a): (1) by-cell turnover and (2) by-bioregion turnover. By-cell turnover represents total biogeographic turnover and is described as the changes of bioregion assignments of each grid cell between consecutive time bins. This metric ranges from 0 (no change) to 1 (complete turnover) and is insensitive to changes in sampling intensity (Kocsis et al. Reference Kocsis, Reddin and Kiessling2018a). By-bioregion turnover estimates the proportions of emergence and disintegration of bioregions to the total number of outlined bioregions in the time interval. In other words, this index traces the origination and extinction of bioregions through time. By-bioregion turnover series are based on subsampling (classical rarefaction) to 10 geographic cells in every age with 100 iterations.

Although the dataset includes a large number of Cambrian data (more than 40,000 fossil occurrences), it remains incomplete. When spatial sampling is sparse, as in our case, the biogeographic membership of a cell can be sensitive to the exact location of a cell boundary when there is considerable turnover over distances that one cell covers. To overcome the biasing effect of this phenomenon, we report time series based on the biogeographic partitioning as averages of 400 replicates with the grid rotated randomly in each iteration. To further determine whether biogeographic patterns change with spatial resolution, the analyses were repeated with different geographic resolutions (edge lengths of 10°, 6.7°, 5.7°, 5°, and 4.4°, respectively).

Provinciality and Geographic Ranges

Following Kocsis et al. (Reference Kocsis, Reddin, Scotese, Valdes and Kiessling2021), we employed Hurlbert's probability of interspecific encounter (PIE) to estimate provinciality. PIE measures the probability that two randomly sampled grid cells belong to different provinces. This provinciality metric ranges from values of 0 (perfectly uniform) to 1 (perfectly uneven). Although direct calculations of species’ ranges (e.g., Stigall and Lieberman Reference Stigall and Lieberman2006; Hendricks et al. Reference Hendricks, Lieberman and Stigall2008; Na and Kiessling Reference Na and Kiessling2015) can provide information about species’ maximum possible geographic distribution, they can be biased. To circumvent the problem of heterogenous spatial sampling, we employed proportional occupancy as a proxy of species’ geographic ranges. Geographic occupancy is here defined as the proportion of all sampled grid cells in which a particular species occurs (Kiessling and Kocsis Reference Kiessling and Kocsis2016). We also calculated species’ stratigraphic durations by the number of stratigraphic intervals during which a species existed. All calculations were implemented in the R programming language and environment (R Core Team 2021).

Results

The first three Cambrian stages, that is, the main phase of diversification during the Cambrian radiation (Na and Kiessling Reference Na and Kiessling2015), saw a prominent increase in the number of bioregions (Fig. 1). The number of bioregions peaked in Stage 4. The total biogeographic turnover increased substantially from Stage 2 until the Wuliuan (Fig. 2). Biogeographic turnover was highest in the Wuliuan and the Paibian stages. Patterns are robust to changes of spatial resolution and minimum occurrence quota (Supplementary Figs. S1, S2). The rate of emergence of bioregions decreased through the first four stages and fluctuated at a lower level afterward, while the rate of disintegration of bioregions increased slightly from Stage 3 to the Guzhangian (Fig. 2).

Figure 1. Paleogeographic positions of sampled bioregions for the 10 Cambrian stages. Numbers and colors indicate time-traceable bioregions. N/A entries denote cells with inadequate sampling (fewer than 10 occurrences). Maps are the PALEOMAP raster series of Scotese (Reference Scotese2016).

Figure 2. By-cell biogeographic turnover and by-bioregion turnover (emergence and disintegration) throughout the Cambrian. The estimates are the averages of 400 random grid rotations. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

The extent of biogeographic changes varied geographically (Fig. 1). The increase in bioregions during the Cambrian radiation was prominent in tropical areas of Laurentia, Siberia, and Gondwana. Tropical areas were also dominated by single-interval species (singletons; Table 1). The distribution of species’ stratigraphic durations was right-skewed in both tropical and nontropical areas, but the mean species’ stratigraphic longevity was slightly shorter in tropical areas than in nontropical areas (Table 1). Similar patterns are also evident in trilobite biogeography (Table 1, Fig. 3). Overall bioregions spanned across continents throughout the Cambrian, but most trilobite bioregions were constrained to a single continent from Stage 3 to the Wuliuan (Fig. 3).

Figure 3. Trilobite biogeographic distribution from Stage 3 to Stage 10. Numbers and colors indicate time-traceable bioregions. N/A entries denote cells with inadequate sampling (fewer than 10 occurrences). Maps are the PALEOMAP raster series of Scotese (Reference Scotese2016).

Table 1. Comparisons of the number in singleton and species median/mean stratigraphic duration between tropical and nontropical areas. Median stratigraphic stages/duration value is calculated by median absolute deviation, whereas mean stratigraphic stages/duration value is standard error of the mean. Tropical and nontropical areas are separated at 30° rotated paleolatitude. Durations are expressed as the number of stages.

The disintegration and subsequent reorganization of bioregions since Stage 3 showed no distinct pattern, but most bioregions seemed to follow expansion–contraction trajectories with respect to their geographic distribution. For example, the fauna of bioregion 1 (red grids in Fig. 1) that persisted from the Fortunian to Stage 4 achieved worldwide distribution (from the tropics to the southern polar region) during Stages 2 and 3 until its disintegration during Stage 4. Only the Wuliuan, the Guzhangian, and the Jiangshanian featured some short-lived (single-interval) bioregions (Fig. 1).

The PIE provinciality metric indicates a profound increase in provinciality from Stage 2 to Stage 3 (Fig. 4). Although the increase of provinciality at the beginning of the Cambrian is concordant with the main diversification phase (Fig. 5), we found no correlation between changes of sampling-standardized species diversity and provinciality (rho = 0.18, p = 0.64), suggesting that the changes in provinciality were not related to changes in species diversity. On the contrary, we found a strong correlation between the provinciality and the number of singletons (rho = 0.78, p = 0.01) over the whole period studied, suggesting that the singletons exerted a potential role in driving Cambrian provincialism.

Figure 4. Trajectory of provinciality from the Fortunian to Stage 10 of the Cambrian, based on Hurlbert's probability of interspecific encounter (PIE). Note the high values during the period from Stage 3 to the Wuliuan. Shading indicates the quantiles of distribution with the 400 random grid rotations. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Figure 5. Global species diversity of marine invertebrates through the Cambrian based on shareholder quorum subsampling with 60% frequency coverage per stage. Error bar is standard deviation of 100 trials. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Species’ geographic occupancy decreased visibly from the Fortunian to Stage 3 (Fig. 6), despite the profound increase in the diversity of mobile species during the same period (Fig. 7). This contraction of geographic occupancy during the critical diversification phase is also evident in subsets of mobile/non-mobile and singleton/non-singleton species (Supplementary Fig. S3).

Figure 6. Box plots of geographic occupancy of species for 10 Cambrian stages. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Figure 7. Trajectory of sample-standardized diversity of actively mobile occurrences in the Cambrian based on shareholder quorum subsampling with 60% frequency coverage per stage. Error bar is standard deviation of 100 trials. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Discussion

Our results confirm a profound increase in provinciality during the Cambrian radiation. This increase of provinciality is paralleled by a decrease of geographic ranges, which is consistent with earlier findings (e.g., Hendricks et al. Reference Hendricks, Lieberman and Stigall2008). A rise of global diversity can be fueled by an increase in provincialism (Valentine et al. Reference Valentine, Foin and Peart1978; Lieberman Reference Lieberman2003; Stigall et al. Reference Stigall, Bauer, Lam and Wright2017). However, we found no correlation between provinciality and diversity over the entire Cambrian.

Our trilobite biogeographic partitioning is largely concordant with traditional trilobite demarcations between Laurentia and Gondwana, but failed to delineate subrealms among tectonic blocks and microcontinents, such as Australia, Avalonia, India, Kazakhstan, and Mongolia (Fig. 3). Moreover, a biogeographic boundary between northern Laurentia and southern Laurentia, as suggested by earlier studies (Meert and Lieberman Reference Meert and Lieberman2004; Álvaro et al. Reference Álvaro, Ahlberg, Babcock, Bordonaro, Choi, Cooper, Ergaliev, Gapp, Pour and Hughes2013), is not discernible during Series 2 and the Miaolingian in our trilobite partitioning outcome, although it is evident in our overall partitioning (Figs. 1, 3). This suggests that the dissimilarity in multi-taxonomic composition within Laurentia was not strongly associated with faunal changeover in local olenelloid trilobites (Pates et al. Reference Pates, Daley, Edgecombe, Cong and Lieberman2019). The differences between the traditional and network-based biogeographic models may reflect limitations of our approach, which is ignorant of phylogenetic histories.

Abrupt biogeographic differentiation during the Cambrian radiation has been suggested by vicariance associated with the breakup of supercontinents (Lieberman Reference Lieberman2005; Meert and Lieberman Reference Meert and Lieberman2008; Na and Kiessling Reference Na and Kiessling2015) such as Pannotia (Scotese Reference Scotese2009). The within-continent variation in bioregions was prevalent since Cambrian Stage 3 (Fig. 1), suggesting that the vicariance process may have occurred not only at continental scales but also at regional and local scales. Within-continent bioregionalization is also evident with different spatial resolution of grid cells (Supplementary Fig. S2) and supported by previous analyses (Lieberman Reference Lieberman1998; Brock et al. Reference Brock, Engelbretsen, Jago, Kruse, Laurie, Shergold, Shi and Sorauf2000; Meert and Lieberman Reference Meert and Lieberman2004; Skovsted Reference Skovsted2004; Paterson et al. Reference Paterson, García-Bellido, Jago, Gehling, Lee and Edgecombe2016).

Across the entire Phanerozoic, both climate and continental configuration have been important drivers of marine provinciality. There is evidence that climate might have been a stronger driver (Kocsis et al. Reference Kocsis, Reddin, Scotese, Valdes and Kiessling2021), with cooling leading to increased provincialism. Determining the relationship between climate and Cambrian provinciality is compromised by limited climate proxy data in the Cambrian. An early Cambrian greenhouse and subsequent cooling is suggested by a few mid-Cambrian geochemical data (Hearing et al. Reference Hearing, Harvey, Williams, Leng, Lamb, Wilby, Gabbott, Pohl and Donnadieu2018; Wotte et al. Reference Wotte, Skovsted, Whitehouse and Kouchinsky2019) and modeling (Hearing et al. Reference Hearing, Pohl, Williams, Donnadieu, Harvey, Scotese, Sepulchre, Franc and Vandenbroucke2021). Goldberg et al. (Reference Goldberg, Present, Finnegan and Bergmann2021) demonstrated a late Cambrian warming, but there were no data points below the Miaolingian. It is thus unclear whether the increasing provinciality across the Cambrian was linked to a long-term cooling trend. The reduced temperature gradients in a warmer world, however, might in part explain the slight decrease in provincialism after the Wuliuan (Fig. 4).

The increase in the number of bioregions during the Cambrian radiation may have been facilitated by morphological and ecological innovations such as dispersal and predation, but underlying mechanisms may vary with time (Nürnberg and Aberhan Reference Nürnberg and Aberhan2013). The global distribution of a single bioregion (Fig. 1) and maximum geographic ranges of species in the Terreneuvian (Fig. 6) could result from a low-competition system (Na and Kiessling Reference Na and Kiessling2015), but could also be due to the lack of diagnostic characters in the earliest Cambrian faunas (Steiner et al. Reference Steiner, Li, Qian, Zhu and Erdtmann2007). The recent discovery of a Cambrian bryozoan in traditional small shelly beds (Ernst and Wilson Reference Ernst and Wilson2021; Zhang et al. Reference Zhang, Zhang, Ma, Taylor, Strotz, Jacquet, Skovsted, Chen, Han and Brock2021) proves a hidden diversity in the earliest Cambrian, which could also account for higher geographic occupancy and lower provinciality in the first two Cambrian stages than afterward. The development of novelties in trophic networks later in Stage 3, such as pervasive predation, would have caused niche overlap and contraction (Na and Kiessling Reference Na and Kiessling2015). This process would ultimately have decreased geographic ranges of species and increased biogeographic differentiation.

The metazoan larvae reported from abundant phosphatized microfossils in strata at the basal Cambrian (e.g., Chen et al. Reference Chen, Oliveri, Li, Zhou, Gao, Hagadorn, Peterson and Davidson2000; Steiner et al. Reference Steiner, Zhu, Li, Qian and Erdtmann2004) indicate a larval revolution (Raff Reference Raff, Maximilian and Littlewood2009). Cosmopolitanism in the Fortunian and Stage 2 could have been achieved by larval dispersal through passive transport (De Bie et al. Reference De Bie, De Meester, Brendonck, Martens, Goddeeris, Ercken, Hampel, Denys, Vanhecke and Van Der Gucht2012). Short larval durations under warm temperature (Álvarez-Noriega et al. Reference Álvarez-Noriega, Burgess, Byers, Pringle, Wares and Marshall2020) might also be associated with the pattern that biogeographic boundaries occurred principally in low latitudes during Stages 3 and 4 (Figs. 1, 3, Supplementary Fig. S2).

Despite an increase in average mobility and perhaps in dispersibility of marine invertebrate species from the Fortunian to Stage 3 (Fig. 7), biogeographic patterns became more pronounced. This may suggest that geographic distribution of species is not dependent on their dispersal ability but on taxonomic duration during which they could persist (Kiessling and Aberhan Reference Kiessling and Aberhan2007). Finding a direct link between dispersibility and geographic distribution could be compromised by our long-term timescale, which is in millions of years. But the biogeographic structure could be partly regulated by some widespread species. For example, the uniform biogeographic structures that were featured since the Wuliuan were likely related to the expanded distribution of the cosmopolitan agnostoid trilobites (Babcock et al. Reference Babcock, Peng and Ahlberg2017), such as Ptychagnostus atavus and Glyptagnostus reticulatus (Zhu et al. Reference Zhu, Babcock and Peng2006).

Another possible explanation for the increase in provinciality might be the large number of species with short stratigraphic duration (e.g., singletons) in the Cambrian (Supplementary Fig. S4), as species might not have attained their potential geographic range simply because they had short life spans (Kiessling and Aberhan Reference Kiessling and Aberhan2007; Foote et al. Reference Foote, Crampton, Beu and Cooper2008). The high proportion of singletons in low latitudes and shorter mean durations (Table 1) suggest a more rapid evolutionary turnover of tropical species, which would reduce geographic ranges and accelerate biogeographic differentiation.

Extinction events strongly affect biogeographic patterns (Kocsis et al. Reference Kocsis, Reddin and Kiessling2018a), and our results suggest that the Cambrian is no exception. Several biological crises are evident in the Cambrian. The best known is the extinction of archaeocyath sponges in Stage 4 (Sinsk event), the end of the Marjumiid Biomere at the Guzhangian/Paibian boundary, and the end of the Pterocephaiid Biomere at the Paibian/Jiangshanian boundary (e.g., Babcock et al. Reference Babcock, Peng, Brett, Zhu, Ahlberg, Bevis and Robison2015). In the aftermath of the extinction of archaeocyath sponges (Zhuravlev and Wood Reference Zhuravlev and Wood1996), a major biogeographic turnover occurred in the Wuliuan (Fig. 2), suggesting that the collapse of reefs exerted a strong influence on biogeographic structures. The other two events paved the way for a more biologically uniform world in the subsequent stages (Paibian, Jiangshanian, and Stage 10), when expansions of oxygen minimum zones, with their corresponding habitat disruption, were probably recurrent across the shallow shelf (Kröger et al. Reference Kröger, Franeck and Rasmussen2019).

Although our network-based approaches have removed some biases arising from an incomplete fossil record and uneven spatial coverage, a series of problems still remain, such as taxonomic issues and stratigraphic resolution. For example, the long stage duration (e.g., Fortunian) might lead to an underestimation of bioregions and an overestimation of the frequency of singletons compared with short stage duration (e.g., Paibian). If evolutionary rates were constant through time, differences in stage duration would severely distort biogeographic structures. However, in our case, given that the Terreneuvian is characterized by low taxonomic turnover rates and high genus survivorships (Kröger et al. Reference Kröger, Franeck and Rasmussen2019), more time slices cut into the Terreneuvian are unlikely to change our biogeographic patterns. On the other hand, Cambrian fossil groups may have been taxonomically oversplit, such as in small shelly faunas (Steiner et al. Reference Steiner, Li, Qian, Zhu and Erdtmann2007), trilobites (Geyer Reference Geyer2019; Zhao et al. Reference Zhao, Liu and Bicknell2020), archaeocyath sponges. (Zhuravlev and Wood Reference Zhuravlev and Wood2020). To what extent this taxonomic issue would impact our biogeographic patterns through time will rely on the variation in proportion of the oversplit taxa among stages, which is worth exploring in the future.

In summary, an increase in regional differences of faunal composition occurred during the Cambrian radiation and was linked with a decrease in by-species geographic distribution. We show that the various ways in which faunal composition was partitioned among regions may have been determined by geography, biotic interaction, climatic changes, and tectonics, providing further evidence for the evolution of biogeographic patterns during and after the Cambrian radiation.

Acknowledgments

We thank M. Hopkins, B. S. Lieberman, and C. M. Ø. Rasmussen for their insightful and helpful reviews, which greatly improved this article. L.N. and Q.L. acknowledge funding from the Youth Innovation Promotion Association of the Chinese Academy of Sciences (CAS; 2019310), and the State Key Laboratory of Paleobiology and Stratigraphy (Nanjing Institute of Geology and Palaeontology, CAS) (20202103), and the CAS (XDB26000000). A.T.K. acknowledges funding from the Deutsche Forschungsgemeinschaft (Ko 5382/1-2 and Ko 5382/2-1). Q.L. is grateful for the invitation to be a guest researcher at GeoZentrum Nordbayern. This is a contribution to IGCP Project 668: “The Stratigraphic and Magmatic History of Early Palaeozoic Equatorial Gondwana and Its Associated Evolutionary Dynamics and IGCP 735 “Rocks and the Rise of Ordovician Life: Filling Knowledge Gaps in the Early Palaeozoic Biodiversification.” This is Paleobiology Database publication no. 436.

Declaration of Competing Interests

The authors declare no competing interests.

Data Availability Statement

Raw data and scripts required to replicate the analysis are accessible through Zenodo: https://doi.org/10.5281/zenodo.6898357. Supplementary File F1: Cambrian biostratigraphic correlation; Supplementary File F2: Supplementary Figs. S1–S4 and Supplementary Table S1.

References

Literature Cited

Alroy, J., Aberhan, M., Bottjer, D. J., Foote, M., Fürsich, F. T., Harries, P. J., Hendy, A. J. W., Holland, S. M., Ivany, L. C., Kiessling, W., Kosnik, M. A., Marshall, C. R., McGowan, A. J., Miller, A. I., Olszewski, T. D., Patzkowsky, M. E., Peters, S. E., Villier, L., Wagner, P. J., Bonuso, N., Borkow, P. S., Brenneis, B., Clapham, M. E., Fall, L. M., Ferguson, C. A., Hanson, V. L., Krug, A. Z., Layou, K. M., Leckey, E. H., Nürnberg, S., Powers, C. M., Sessa, J. A., Simpson, C., Tomašových, A., and Visaggi, C. C.. 2008. Phanerozoic trends in the global diversity of marine invertebrates. Science 321:97100.CrossRefGoogle ScholarPubMed
Álvarez-Noriega, M., Burgess, S. C., Byers, J. E., Pringle, J. M., Wares, J. P., and Marshall, D. J.. 2020. Global biogeography of marine dispersal potential. Nature Ecology and Evolution 4:11961203.CrossRefGoogle ScholarPubMed
Álvaro, J. J., Ahlberg, P., Babcock, L. E., Bordonaro, O. L., Choi, D. K., Cooper, R. A., Ergaliev, G. K., Gapp, I. W., Pour, M. G., and Hughes, N. C.. 2013. Global Cambrian trilobite palaeobiogeography assessed using parsimony analysis of endemicity. Geological Society of London Memoir 38:273296.CrossRefGoogle Scholar
Babcock, L. E., Peng, S.-C., Brett, C. E., Zhu, M.-Y., Ahlberg, P., Bevis, M., and Robison, R. A.. 2015. Global climate, sea level cycles, and biotic events in the Cambrian Period. Palaeoworld 24(1–2):515.CrossRefGoogle Scholar
Babcock, L. E., Peng, S.-C., and Ahlberg, P.. 2017. Cambrian trilobite biostratigraphy and its role in developing an integrated history of the Earth system. Lethaia 50:381399.CrossRefGoogle Scholar
Brock, G., Engelbretsen, M., Jago, J., Kruse, P., Laurie, J., Shergold, J., Shi, G., and Sorauf, J.. 2000. Palaeobiogeographic affinities of Australian Cambrian faunas. Memoirs of the Association of Australasian Palaeontologists 23:161.Google Scholar
Chen, J.-Y., Oliveri, P., Li, C.-W., Zhou, G.-Q., Gao, F., Hagadorn, J. W., Peterson, K. J., and Davidson, E. H.. 2000. Precambrian animal diversity: putative phosphatized embryos from the Doushantuo Formation of China. Proceedings of the National Academy of Sciences USA 97:44574462.CrossRefGoogle ScholarPubMed
Close, R., Benson, R., Saupe, E., Clapham, M., and Butler, R.. 2020. The spatial structure of Phanerozoic marine animal diversity. Science 368:420424.CrossRefGoogle ScholarPubMed
Crame, J. A., and Owen, A. W.. 2002. Palaeobiogeography and biodiversity change: the Ordovician and Mesozoic–Cenozoic radiations. Geological Society of London Special Publication 194.Google Scholar
De Bie, T., De Meester, L., Brendonck, L., Martens, K., Goddeeris, B., Ercken, D., Hampel, H., Denys, L., Vanhecke, L., and Van Der Gucht, K.. 2012. Body size and dispersal mode as key traits determining metacommunity structure of aquatic organisms. Ecology Letters 15:740747.CrossRefGoogle Scholar
Dunne, E. M., Close, R. A., Button, D. J., Brocklehurst, N., Cashmore, D. D., Lloyd, G. T., and Butler, R. J.. 2018. Diversity change during the rise of tetrapods and the impact of the “Carboniferous rainforest collapse.” Proceedings of the Royal Society of London B 285:20172730.Google ScholarPubMed
Edler, D., Guedes, T., Zizka, A., Rosvall, M., and Antonelli, A.. 2016. Infomap bioregions: interactive mapping of biogeographical regions from species distributions. Systematic Biology 66:197204.Google Scholar
Ernst, A., and Wilson, M. A.. 2021. Bryozoan fossils found at last in deposits from the Cambrian period. Nature 599:203204.CrossRefGoogle ScholarPubMed
Erwin, D. H., Laflamme, M., Tweedt, S. M., Sperling, E. A., Pisani, D., and Peterson, K. J.. 2011. The Cambrian conundrum: early divergence and later ecological success in the early history of animals. Science 334:10911097.CrossRefGoogle ScholarPubMed
Fan, J.-X., Shen, S.-Z., Erwin, D. H., Sadler, P. M., MacLeod, N., Cheng, Q.-M., Hou, X.-D., Yang, J., Wang, X.-D., Wang, Y., Zhang, H., Chen, X., Li, G.-X., Zhang, Y.-C., Shi, Y.-K., Yuan, D.-X., Chen, Q., Zhang, L.-N., Li, C., and Zhao, Y.-Y.. 2020. A high-resolution summary of Cambrian to Early Triassic marine invertebrate biodiversity. Science 367:272277.CrossRefGoogle ScholarPubMed
Foote, M., Crampton, J. S., Beu, A. G., and Cooper, R. A.. 2008. On the bidirectional relationship between geographic range and taxonomic duration. Paleobiology 34:421433.CrossRefGoogle Scholar
Fortey, R. A., Briggs, D. E. G., and Wills, M. A.. 1996. The Cambrian evolutionary “explosion”: decoupling cladogenesis from morphological disparity. Biological Journal of the Linnean Society 57:1333.Google Scholar
Geyer, G. 2019. A comprehensive Cambrian correlation chart. Episodes: Journal of International Geoscience 42:321332.Google Scholar
Goldberg, S. L., Present, T. M., Finnegan, S., and Bergmann, K. D.. 2021. A high-resolution record of early Paleozoic climate. Proceedings of the National Academy of Sciences USA 118:e2013083118.CrossRefGoogle ScholarPubMed
Gradstein, F. M., Ogg, J. G., Schmitz, M. D., and Ogg, G. M.. 2020. Geologic time scale 2020. Elsevier, Amsterdam.Google Scholar
Hearing, T. W., Harvey, T. H. P., Williams, M., Leng, M. J., Lamb, A. L., Wilby, P. R., Gabbott, S. E., Pohl, A., and Donnadieu, Y.. 2018. An early Cambrian greenhouse climate. Science Advances 4:eaar5690.CrossRefGoogle ScholarPubMed
Hearing, T. W. W., Pohl, A., Williams, M., Donnadieu, Y., Harvey, T. H., Scotese, C. R., Sepulchre, P., Franc, A., and Vandenbroucke, T. R.. 2021. Quantitative comparison of geological data and model simulations constrains early Cambrian geography and climate. Nature Communications 12:111.Google Scholar
Hendricks, J. R., Lieberman, B. S., and Stigall, A. L.. 2008. Using GIS to study palaeobiogeographic and macroevolutionary patterns in soft-bodied Cambrian arthropods. Palaeogeography, Palaeoclimatology, Palaeoecology 264(1–2):163175.CrossRefGoogle Scholar
Huang, B., Jin, J., and Rong, J.-Y.. 2018. Post-extinction diversification patterns of brachiopods in the early–middle Llandovery, Silurian. Palaeogeography, Palaeoclimatology, Palaeoecology 493:1119.CrossRefGoogle Scholar
Jago, J. B., Zang, W.-L., Sun, X., Brock, G. A., Paterson, J. R., and Skovsted, C. B.. 2006. A review of the Cambrian biostratigraphy of South Australia. Palaeoworld 15(3–4):406423.CrossRefGoogle Scholar
Jell, P. A. 1974. Faunal provinces and possible planetary reconstruction of the Middle Cambrian. Journal of Geology 82:319350.CrossRefGoogle Scholar
Kiel, S. 2017. Using network analysis to trace the evolution of biogeography through geologic time: a case study. Geology 45:711714.Google Scholar
Kiessling, W., and Aberhan, M.. 2007. Geographical distribution and extinction risk: lessons from Triassic–Jurassic marine benthic organisms. Journal of Biogeography 34:14731489.CrossRefGoogle Scholar
Kiessling, W., and Kocsis, Á. T.. 2016. Adding fossil occupancy trajectories to the assessment of modern extinction risk. Biology Letters 12:20150813.CrossRefGoogle Scholar
Kocsis, Á. T. 2017. The R package icosa: coarse resolution global triangular and penta-hexagonal grids based on tessellated icosahedra, R package version 0.9.81. https://CRAN.R-project.org/package=icosa, accessed 13 December 2019.Google Scholar
Kocsis, Á. T., Reddin, C. J., and Kiessling, W.. 2018a. The biogeographical imprint of mass extinctions. Proceedings of the Royal Society of London B 285:20180232.Google ScholarPubMed
Kocsis, Á. T., Reddin, C. J., and Kiessling, W.. 2018b. The stability of coastal benthic biogeography over the last 10 million years. Global Ecology and Biogeography 27:11061120.CrossRefGoogle Scholar
Kocsis, Á. T., Reddin, C. J., Alroy, J., and Kiessling, W.. 2019. The R package divDyn for quantifying diversity dynamics using fossil sampling data. Methods in Ecology and Evolution 10:735743.CrossRefGoogle Scholar
Kocsis, Á. T., Reddin, C. J., Scotese, C. R., Valdes, P. J., and Kiessling, W.. 2021. Increase in marine provinciality over the last 250 million years governed more by climate change than plate tectonics. Proceedings of the Royal Society of London B 288:20211342.Google ScholarPubMed
Kröger, B., Franeck, F., and Rasmussen, C. M. Ø.. 2019. The evolutionary dynamics of the early Palaeozoic marine biodiversity accumulation. Proceedings of the Royal Society of London B 286:20191634.Google ScholarPubMed
Landing, E., Westrop, S. R., and Bowring, S. A.. 2013. Reconstructing the Avalonia palaeocontinent in the Cambrian: a 519 Ma caliche in South Wales and transcontinental middle Terreneuvian sandstones. Geological Magazine 150:10221046.CrossRefGoogle Scholar
Li, G.-X., Steiner, M., Zhu, X.-J., Yang, A.-H., Wang, H.-F., and Erdtmann, B.. 2007. Early Cambrian metazoan fossil record of South China: generic diversity and radiation patterns. Palaeogeography, Palaeoclimatology, Palaeoecology 254:229249.CrossRefGoogle Scholar
Lieberman, B. S. 1997. Early Cambrian paleogeography and tectonic history: a biogeographic approach. Geology 25:10391042.2.3.CO;2>CrossRefGoogle Scholar
Lieberman, B. S. 1998. Cladistic analysis of the Early Cambrian olenelloid trilobites. Journal of Paleontology 72:5978.CrossRefGoogle Scholar
Lieberman, B. S. 2003. Taking the pulse of the Cambrian radiation. Integrative and Comparative Biology 43:229–37.CrossRefGoogle ScholarPubMed
Lieberman, B. S. 2005. Geobiology and paleobiogeography: tracking the coevolution of the Earth and its biota. Palaeogeography, Palaeoclimatology, Palaeoecology 219(1–2):2333.CrossRefGoogle Scholar
Lu, Y.-H. 1981. Provincialism, dispersal, development, and phylogeny of trilobites. Geological Society of America Special Paper 187:143152.Google Scholar
Lu, Y.-H., Zhu, Z.-L., Qian, Y.-Y., Lin, H.-L., Zhou, Z.-Y., and Yuan, K.-X.. 1974. Bio-environmental control hypothesis and its application to the Cambrian biostratigraphy and palaeozoogeography. Memoir of the Nanjing Institute of Geology and Palaeontology, Academia Sinica 5:27110.Google Scholar
Maloof, A. C., Porter, S. M., Moore, J. L., Dudas, F. O., Bowring, S. A., Higgins, J. A., Fike, D. A., and Eddy, M. P.. 2010. The earliest Cambrian record of animals and ocean geochemical change. Geological Society of America Bulletin 122(11–12):17311774.CrossRefGoogle Scholar
Meert, J. G., and Lieberman, B. S.. 2004. A palaeomagnetic and palaeobiogeographical perspective on latest Neoproterozoic and early Cambrian tectonic events. Journal of the Geological Society 161:477487.CrossRefGoogle Scholar
Meert, J. G., and Lieberman, B. S.. 2008. The Neoproterozoic assembly of Gondwana and its relationship to the Ediacaran–Cambrian radiation. Gondwana Research 14:521.CrossRefGoogle Scholar
Müller, R. D., Cannon, J., Qin, X., Watson, R. J., Gurnis, M., Williams, S., Pfaffelmoser, T., Seton, M., Russell, S. H., and Zahirovic, S.. 2018. GPlates: building a virtual Earth through deep time. Geochemistry, Geophysics, Geosystems 19:22432261.CrossRefGoogle Scholar
Na, L., and Kiessling, W.. 2015. Diversity partitioning during the Cambrian radiation. Proceedings of the National Academy of Sciences USA 112:47024706.CrossRefGoogle ScholarPubMed
Nürnberg, S., and Aberhan, M.. 2013. Habitat breadth and geographic range predict diversity dynamics in marine Mesozoic bivalves. Paleobiology 39:360372.CrossRefGoogle Scholar
Palmer, A. R., and Repina, L. N.. 1993. Through a glass darkly: taxonomy, phylogeny, and biostratigraphy of the Olenellina. University of Kansas, Lawrence.Google Scholar
Paterson, J. R., García-Bellido, D. C., Jago, J. B., Gehling, J. G., Lee, M. S., and Edgecombe, G. D.. 2016. The Emu Bay Shale Konservat-Lagerstätte: a view of Cambrian life from East Gondwana. Journal of the Geological Society 173:111.CrossRefGoogle Scholar
Pates, S., Daley, A. C., Edgecombe, G. D., Cong, P., and Lieberman, B. S.. 2019. Systematics, preservation and biogeography of radiodonts from the southern Great Basin, USA, during the upper Dyeran (Cambrian Series 2, Stage 4). Papers in Palaeontology 7:235262.CrossRefGoogle Scholar
Peng, S.-C., Babcock, L., and Cooper, R.. 2012. The Cambrian Period. Pp. 3746 in Ogg, J. G., Ogg, G. M., and Gradstein, F. M., eds. The concise geological timescale. Cambridge University Press, Cambridge.Google Scholar
Pillola, G. L. 1989. Trilobites du Cambrien inférieur du SW de la Sardaigne, Italie. Palaeontographia Italica 78:1173.Google Scholar
Raff, R. A. 2009. Origins of metazoan body plans: the larval revolution. Pp. 4351 in Maximilian, J. T. and Littlewood, D. T. J., eds. Animal evolution. Oxford University Press, New York.CrossRefGoogle Scholar
Rasmussen, C. M. Ø., Kröger, B., Nielsen, M. L., and Colmenar, J.. 2019. Cascading trend of Early Paleozoic marine radiations paused by Late Ordovician extinctions. Proceedings of the National Academy of Sciences USA 116:72077213.CrossRefGoogle ScholarPubMed
R Core Team. 2021. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.Google Scholar
Rojas, A., Patarroyo, P., Mao, L., Bengtson, P., and Kowalewski, M.. 2017. Global biogeography of Albian ammonoids: a network-based approach. Geology 45:659662.CrossRefGoogle Scholar
Rong, J.-Y., Harper, D., Huang, B., Li, R.-Y., Zhang, X.-L., and Chen, D.. 2020. The latest Ordovician Hirnantian brachiopod faunas: new global insights. Earth-Science Reviews 208:103280.CrossRefGoogle Scholar
Rosvall, M., and Bergstrom, C. T.. 2008. Maps of random walks on complex networks reveal community structure. Proceedings of the National Academy of Sciences USA 105:11181123.CrossRefGoogle ScholarPubMed
Scotese, C. R. 2009. Late Proterozoic plate tectonics and palaeogeography: a tale of two supercontinents, Rodinia and Pannotia. Geological Society of London Special Publication 326:6783.CrossRefGoogle Scholar
Scotese, C. R. 2016. PALEOMAP PaleoAtlas for GPlates and the PaleoData Plotter Program. https://www.earthbyte.org/paleomap-paleoatlas-for-gplates, accessed 20 July 2021.Google Scholar
Sepkoski, J. J. Jr. 1997. Biodiversity: past, present, and future. Journal of Paleontology 71:533539.CrossRefGoogle ScholarPubMed
Sidor, C. A., Vilhena, D. A., Angielczyk, K. D., Huttenlocker, A. K., Nesbitt, S. J., Peecook, B. R., Steyer, J. S., Smith, R. M., and Tsuji, L. A.. 2013. Provincialization of terrestrial faunas following the end-Permian mass extinction. Proceedings of the National Academy of Sciences USA 110:81298133.CrossRefGoogle ScholarPubMed
Skovsted, C. B. 2004. Mollusc fauna of the Early Cambrian Bastion Formation of north-east Greenland. Bulletin of the Geological Society of Denmark 51:1137.CrossRefGoogle Scholar
Steiner, M., Zhu, M.-Y., Li, G.-X., Qian, Y., and Erdtmann, B. D.. 2004. New early Cambrian bilaterian embryos and larvae from China. Geology 32:833836.CrossRefGoogle Scholar
Steiner, M., Li, G.-X., Qian, Y., Zhu, M.-Y., and Erdtmann, B. D.. 2007. Neoproterozoic to early Cambrian small shelly fossil assemblages and a revised biostratigraphic correlation of the Yangtze Platform (China). Palaeogeography Palaeoclimatology Palaeoecology 254(1–2):6799.CrossRefGoogle Scholar
Steiner, M., Hohl, S. V., Yang, B., Huang, X.-T., and Li, D.. 2021. Rewriting the Cambrian biogeography of the Central Asian Orogenic Belt using combined faunal cluster, zircon age and C isotope analysis. Geophysical Research Letters 48:e2021GL093133.CrossRefGoogle Scholar
Stigall, A. L., and Lieberman, B. S.. 2006. Quantitative palaeobiogeography: GIS, phylogenetic biogeographical analysis, and conservation insights. Journal of Biogeography 33:20512060.CrossRefGoogle Scholar
Stigall, A. L., Bauer, J. E., Lam, A. R., and Wright, D. F.. 2017. Biotic immigration events, speciation, and the accumulation of biodiversity in the fossil record. Global and Planetary Change 148:242257.CrossRefGoogle Scholar
Theokritoff, G. 1979. Early Cambrian provincialism and biogeographic boundaries in the North Atlantic region. Lethaia 12:281295.CrossRefGoogle Scholar
Valentine, J. W., Foin, T. C., and Peart, D.. 1978. A provincial model of Phanerozoic marine diversity. Paleobiology 4:5566.CrossRefGoogle Scholar
Vilhena, D. A., and Antonelli, A.. 2015. A network approach for identifying and delimiting biogeographical regions. Nature Communications 6:6848.CrossRefGoogle ScholarPubMed
Vilhena, D. A., Harris, E. B., Bergstrom, C. T., Maliska, M. E., Ward, P. D., Sidor, C. A., Strömberg, C. A., and Wilson, G. P.. 2013. Bivalve network reveals latitudinal selectivity gradient at the end-Cretaceous mass extinction. Scientific Reports 3:1790.CrossRefGoogle Scholar
Wotte, T., Skovsted, C. B., Whitehouse, M. J., and Kouchinsky, A.. 2019. Isotopic evidence for temperate oceans during the Cambrian Explosion. Scientific Reports 9:19.CrossRefGoogle ScholarPubMed
Yang, B., Steiner, M., and Keupp, H.. 2015. Early Cambrian palaeobiogeography of the Zhenba–Fangxian Block (South China): independent terrane or part of the Yangtze Platform? Gondwana Research 28:15431565.CrossRefGoogle Scholar
Zhang, Z.-L., Zhang, Z.-F., Ma, J.-Y., Taylor, P. D., Strotz, L. C., Jacquet, S. M., Skovsted, C. B., Chen, F.-Y., Han, J., and Brock, G. A.. 2021. Fossil evidence unveils an early Cambrian origin for Bryozoa. Nature 599:251255.CrossRefGoogle ScholarPubMed
Zhao, W.-Y., Liu, J.-N, and Bicknell, R. D.. 2020. Geometric morphometric assessment of Guanshan trilobites (Yunnan Province, China) reveals a limited diversity of palaeolenid taxa. Palaeontologia Electronica 23:115.Google Scholar
Zhu, M.-Y., Babcock, L. E., and Peng, S.-C.. 2006. Advances in Cambrian stratigraphy and paleontology: integrating correlation techniques, paleobiology, taphonomy and paleoenvironmental reconstruction. Palaeoworld 15(3–4):217222.CrossRefGoogle Scholar
Zhuravlev, A. Y. 1986. Evolution of archaeocyaths and palaeobiogeography of the Early Cambrian. Geological Magazine 123:377385.CrossRefGoogle Scholar
Zhuravlev, A. Y., and Riding, R.. 2001. The ecology of the Cambrian radiation. Columbia University Press, New York.Google Scholar
Zhuravlev, A. Y., and Wood, R. A.. 1996. Anoxia as the cause of the mid-early Cambrian (Botomian) extinction event. Geology 24:311314.2.3.CO;2>CrossRefGoogle Scholar
Zhuravlev, A. Y., and Wood, R.. 2020. Dynamic and synchronous changes in metazoan body size during the Cambrian Explosion. Scientific Reports 10:18.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Paleogeographic positions of sampled bioregions for the 10 Cambrian stages. Numbers and colors indicate time-traceable bioregions. N/A entries denote cells with inadequate sampling (fewer than 10 occurrences). Maps are the PALEOMAP raster series of Scotese (2016).

Figure 1

Figure 2. By-cell biogeographic turnover and by-bioregion turnover (emergence and disintegration) throughout the Cambrian. The estimates are the averages of 400 random grid rotations. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Figure 2

Figure 3. Trilobite biogeographic distribution from Stage 3 to Stage 10. Numbers and colors indicate time-traceable bioregions. N/A entries denote cells with inadequate sampling (fewer than 10 occurrences). Maps are the PALEOMAP raster series of Scotese (2016).

Figure 3

Table 1. Comparisons of the number in singleton and species median/mean stratigraphic duration between tropical and nontropical areas. Median stratigraphic stages/duration value is calculated by median absolute deviation, whereas mean stratigraphic stages/duration value is standard error of the mean. Tropical and nontropical areas are separated at 30° rotated paleolatitude. Durations are expressed as the number of stages.

Figure 4

Figure 4. Trajectory of provinciality from the Fortunian to Stage 10 of the Cambrian, based on Hurlbert's probability of interspecific encounter (PIE). Note the high values during the period from Stage 3 to the Wuliuan. Shading indicates the quantiles of distribution with the 400 random grid rotations. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Figure 5

Figure 5. Global species diversity of marine invertebrates through the Cambrian based on shareholder quorum subsampling with 60% frequency coverage per stage. Error bar is standard deviation of 100 trials. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Figure 6

Figure 6. Box plots of geographic occupancy of species for 10 Cambrian stages. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.

Figure 7

Figure 7. Trajectory of sample-standardized diversity of actively mobile occurrences in the Cambrian based on shareholder quorum subsampling with 60% frequency coverage per stage. Error bar is standard deviation of 100 trials. Dru, Drumian; For, Fortunian; Guz, Guzhangian; Jia, Jiangshanian; Pai, Paibian; Wul, Wuliuan; St2, Stage 2; St3, Stage 3; St4, Stage 4; St10, Stage 10.