Connections between Community, Ecosystem, Linnaean, “Geological Time Scale,” and Bretskyan Hierarchies

The Linnaean or genealogical hierarchy has received considerable interest in debates on the levels of selection, and evolution more generally (Vrba and Gould Reference Vrba and Gould1986; Lloyd Reference Lloyd1994; Gould and Lloyd Reference Gould and Lloyd1999; Gould Reference Gould2002; Okasha Reference Okasha2006, Reference Okasha2012; Simpson Reference Simpson2011). One of the major discussion points is where we should find the focus of selection processes, and which entities are selected and which are sorted or receiving side effects of selection events (Eldredge Reference Eldredge, Jablonski, Erwin and Lipps1996; Jablonski Reference Jablonski2008b). In any case, the Linnaean hierarchy, even though reticulated by myriad hybridization and lateral gene transfer events (Arnold Reference Arnold2006), shows a distinct treelike structure and represents a hierarchy of informational inheritance. The economic hierarchy is more complicated, because (1) it can be represented as a hierarchy of interacting populations and communities (Eldredge Reference Eldredge2008), or (2) alternatively, it can be seen as a hierarchy of ecosystems (O'Neill et al. Reference O'Neill, DeAngelis, Waide and Allen1986). The ecosystem view is even more complex, because it also encompasses abiotic energy and matter transfer, which are tied to geological processes that span many time scales and interact not only horizontally (in geographic space of longitude [λ] and latitude [φ] but also in the third dimension [z] across all planetary spheres—crust, mantle, and core). There is no doubt that, at the long time scales of thousands to billions of years, life is an integral part of the Earth's dynamics, with feedbacks working in both directions and most clearly apparent on time scales longer than millions of years (Lovejoy Reference Lovejoy2015; Wignall Reference Wignall2015; Knoll and Nowak Reference Knoll and Nowak2017; Spiridonov et al. Reference Spiridonov, Venckutė-Aleksienė and Radzevičius2017c, Reference Spiridonov, Balakauskas and Lovejoy2022; Eichenseer et al. Reference Eichenseer, Balthasar, Smart, Stander, Haaga and Kiessling2019; Anderson Reference Anderson, Foulger, Hamilton, Jurdy, Stein, Howard and Stein2022; Sánchez-Baracaldo et al. Reference Sánchez-Baracaldo, Bianchini, Wilson and Knoll2022; Spencer et al. Reference Spencer, Davies, Gernon, Wang, McMahon, Morrell, Hincks, Pufahl, Brasier and Seraine2022; Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022, Reference Spiridonov and Lovejoy2023). Therefore, the hierarchy of ecosystems, and especially the dynamics of their spatial and compositional structure on long time scales, should be understood as a part of the geological/intra-planetary hierarchy of processes and structures. The Bretskyan hierarchy traces fates and spatial configurations of component populations of interacting functional and genealogical collectives. It combines approaches of all major biological hierarchies—community hierarchy of biotic functional interactions, ecosystem hierarchy of Earth–life interactions, and genealogical hierarchy of patterns of information transfer.

To distinguish the functional hierarchy of biocoenoses (Tëmkin Reference Tëmkin, Brooks, DiFrisco and Wimsatt2021), which we here call the Vernadskyan hierarchy, from the historical hierarchy of eco-genealogical groupings, we present here a reassessed Bretskyan hierarchy (Table 1). The Bretskyan units are formed by proximately interacting biota co-occurring in contingently determined time and space. The fact that most Bretskyan units are polyphyletic stems from Herbert Spencer's principle of the “instability of the homogeneous”—the basis of the so-called biology's first law (McShea and Brandon Reference McShea and Brandon2010; Brandon and McShea Reference Brandon and McShea2020), and also the second law of thermodynamics as applied to collectives—the random movement of individuals in space from different groups results in progressive mixing (increase in disorder of configurations). The unmixing requires a directional effort (e.g., spatial sorting of some sort [Shine et al. Reference Shine, Brown and Phillips2011]). It should be noted that Bretskyan units can be composed of monophyletic units, but this requires very specific mechanisms of aggregation, such as formation of kin groups such as human families (but see Rosenberg [2022] for a different view). The Bretskyan units at smaller spatial scales can be highly integrated and thus described as holobionts or even holobiont-derived organisms in a strict sense (such as eukaryotes). The Bretskyan units, which form at larger scales and are composed of thousands to millions of species, are separated from other such units by externally forced isolation by means of geographic barriers. The large-scale Bretskyan units can be subunits of still larger units up to the ultimate scale of the largest geobiomic unit that also has a phylogenetic significance—Gaia herself. To reiterate, here the term “Gaia” is used as the given name for the largest Bretskyan unit on our planet, which is characterized by both ecosystem and genealogical properties.

Geobiomes have one significant difference from holobionts, in that their biotic parts are usually sparsely dispersed in the abiotic environment and a significant number of interactions at all scales is dominated by interactions with that physical environment. Therefore, the physical environment is at the same time the embedding matrix for a geobiome and a causal constituent part of it. Geobiomes are principal geobiological individuals. At larger scales, they have a potential to change physical conditions to a significant degree and thus are self-modifying vectors of their own evolution (Kiessling Reference Kiessling2009; Knoll and Nowak Reference Knoll and Nowak2017)—they, to a significant degree, are causes and effects of their own change.

The central tenet of the Bretskyan hierarchy is the integration of polyphyletic individuals, by means of isolation, which can be sustained by means of biotic interactions or by means of externally enforced formidable physical barriers. The large geological structures of greater size or areal extent, belonging to the geological time scale hierarchy, are longer lasting than smaller geological structures—oceans last longer than seas, which last longer than lakes; and continents last longer than large islands, and large islands last longer than small islands due to works of erosion. While cores of continents last eons—billions of years (Nance Reference Nance2022), oceanic islands usually last millions to tens of millions of years (Clague and Dalrymple Reference Clague and Dalrymple1987). Similarly, ocean basins, which are larger than epicontinental seas by orders of magnitude, last tens of millions to hundreds of millions of years (Björck Reference Björck1995; Woodcock Reference Woodcock2004). This is a hierarchy of sizes and, therefore, significance of effects. Also, smaller structures are parts of larger structures, and therefore represent a case of compositional scalar hierarchy (smaller things inside larger things). For example, a pond could be a part of a terrestrial valley, which is a part of a plateau. A pond exists for a shorter period of time than a valley, which exists for a shorter period of time than a plateau. The geographic area and the duration of Bretskyan units are therefore congruently determined by scaling of sizes and durations of geological/geomorphic structures. The phylogenetic congruence of biota inside Bretskyan units can be measured by proportion (and ranks) of endemic taxa. At the scale of the planet, the entire clade “Earth's life” becomes spatially located, and therefore is equal to the largest Bretskyan individual (Fig. 1). At the singularity of the largest scale, contiguous assemblage of all taxa becomes completely monophyletic.

Figure 1. The congruence of the Bretskyan hierarchy of geobiomes with the Linnaean hierarchy of taxa and clades. At the smallest scales, random mixing dominates and the significance and duration of physical barriers of dispersal of organisms are low. The probability of encountering barriers and the duration of these barriers increases as a function of scale—larger physical structures have a higher chance of persistence. Isolated islands, lakes, epeiric seas, or whole continents will impose multimillion year congruence in phylogenies of clades that occupy these areas. The larger the Bretskyan units are, the longer their probable persistence time and the higher their phylogenetic significance. At the largest scale of the planet, the Bretskyan and Linnaean hierarchies converge into singularity of the largest clade and the largest geobiome—Gaia. Because taxa at high taxonomic levels are spatially well mixed on our planet at smaller scales, when we approach the size of the whole planet, the congruence of Bretskyan units and Linnaean taxa should rise sharply, thus forming a concave upward relationship.

Geological time scales are constructed by recognizing the hierarchy of significance of events subdividing Earth's history into time intervals typified by unique features. These geochronological subdivisions have their stratigraphic equivalents represented as sets of all geological bodies formed at given time intervals, from the largest rank to the smallest: eonothems, erathems, systems, series, and stages (Ogg et al. Reference Ogg, Ogg and Gradstein2008). Time and structure are intrinsically linked in geology and reflect the hierarchical nature of planetary change. No scientific discipline better reflects this than sequence stratigraphy, which ties patterns of structure to concrete physical events and uses them to recognize hierarchically organized patterns, repeated at many scales, for geological time subdivision and correlation of isolated strata (Catuneanu Reference Catuneanu2006; Miall Reference Miall2010). The progress of geological mapping has enabled the restoration of paleogeography (paleotectonic configurations) in 4D, thus recognizing the hierarchical composition, assembly, and disassembly of continents, oceans, and smaller terrains (Torsvik and Cocks Reference Torsvik and Cocks2016).

The topographic barriers caused by compositional differences and physical forces affecting the solid Earth create not only barriers for the dispersal of biota but also compartmentalize geological and other geophysical processes in space as well as time, thus creating a hierarchy of sedimentary basins and regional to local climates. Therefore these ≈2D features continue for significant periods of time, thus forming branching and merging 3D anastomosing (reticulated treelike) networks of ever-changing geological structures and processes. The geobiological Bretskyan hierarchy essentially inherits this structure. Currently, the closest approach that tackles such structures is so-called macrostratigraphy, which quantifies gap-bound or hiatus-bound sedimentary “packages”—archives of relatively homogenous paleogeographies at sub-continental scales (Meyers and Peters Reference Meyers and Peters2011). This information can be used in studying spatiotemporal changes in Earth systems (Peters et al. Reference Peters, Quinn, Husson and Gaines2022). Lithosomes—geological formations, or even higher-order lithologically defined bodies such as groups or supergroups (Murphy and Salvador Reference Murphy and Salvador1999), which represent time intervals of tens to hundreds of millions of years at still larger spatial scales—are the archives of compartments of Bretskyan units that preserve information on physical conditions and sizes of geobiomes.

Longer-lasting geobiomes experience more rounds of microevolution and macroevolution. The exact shape of the spatial scaling relation between the geographic sizes/areas of Bretskyan units and the coherence between Bretskyan units and Linnaean taxonomy will depend on the geophysical state of dynamics of the planet. For example, if the plate tectonics is intermittent with the “single lid” tectonics regime (as appears to be the case [Stern Reference Stern2023], then we should expect extreme differences in levels of geobiomic individuality and consequently rates and patterns of coevolution between the Mesoproterozoic and Neoproterozoic–Phanerozoic eons.

This view of connection between spatial and temporal scales resonates well with the “sloshing bucket” model of evolution (Eldredge Reference Eldredge1999, Reference Eldredge, Crutchfield and Schuster2003), wherein large-scale changes in taxonomic composition and ecological functioning need proportionally larger external perturbations (Fig. 2). Large-scale perturbations are either rare and their recurrence time is long (Raup and Sepkoski Reference Raup and Sepkoski1982, Reference Raup and Sepkoski1984; Lieberman and Melott Reference Lieberman and Melott2007) or of long duration (Lovejoy Reference Lovejoy2015; Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022).

Figure 2. Conceptual diagram that frames the relations between the economic hierarchy of functional communities, the Linnaean hierarchy of clades, and the Bretskyan hierarchy of holobionts and geobiomes. Community units have high chances for evolutionary integration into discrete individual-like modules at the very smallest and at the largest scales. At the smallest scales, the unity and integration is achieved by means of cooperation or other forms of symbiosis between limited numbers of individuals at that level. The dotted line shows the zone where classical “major transitions in individuality” are probable: (1) egalitarian in the Bretskyan hierarchy and (2) fraternal in the Linnaean hierarchy. As systems become larger, the number of components (interacting species) grows, as does the variance in their traits, which makes integration of polyphyletic and polymorphic individuals harder and harder. Then, at still larger spatial and temporal scales, another and unrelated factor starts to play the prominent role in integration of polyphyletic ecological communities into evolutionary units—the force external to the biota itself—spatial isolation, that is, geology in a broad sense. Because the size of spatial barriers is directly related to their durations, this means larger and more diverse ecological entities are more likely to be isolated for longer periods of time from any such similarly sized (hierarchically comparable) entities. At larger sizes of geobiomes, effective integrating coevolution will be more and more visible if we compare biogeographic units of increasing size. The changes in such units can be apparent, due to their size inertia, either through a long time or by means of rare but exceedingly powerful events—the “sloshing bucket” mechanism of larger entities needing larger/longer perturbation for initiating evolutionary change (Eldredge Reference Eldredge1999, Reference Eldredge, Crutchfield and Schuster2003). As in the case of water level in the metaphorical “bucket,” we need sufficiently large in magnitude perturbation in order to permanently change the volume of the liquid or, in the case of geobiomic evolution, the composition of biota in a province.

The Bretskyan hierarchy of holobionts and geobiomes connects in the conceptual understanding of the major transitions in individuality (Buss Reference Buss1987; Maynard Smith and Szathmáry Reference Maynard Smith and Szathmáry1998) with the existence of biological hierarchies (Eldredge Reference Eldredge, Jablonski, Erwin and Lipps1996; Tëmkin and Eldredge Reference Tëmkin, Eldredge, Serrelli and Gontier2015). In the presented view, borrowing David Queller's terminology (Queller Reference Queller2000), the “egalitarian” major transitions (integration of polyphyletic communities into higher-level individuals) as opposed to the “fraternal transitions” of genetically related individuals in forming higher-level integrated evolutionary entities can be understood as self-isolation and self-integration processes working in the lower tiers of the Bretskyan hierarchy. Because the most species-rich clade of life, the eukaryotes, are a result of one such historical individuation event, we can posit that Bretskyan and Linnaean kinds of hierarchies can be nontrivially connected and supply individuals in both directions in a rank-free way (Bretskyan units ⇄ Linnaean units). When we are describing Bretskyan and Linnaean kinds of hierarchies, we mean that the same logic of subdivision into ancestry-based treelike patterns (Linnaean kind) and spatially and temporally isolated and coevolution and/or physical isolation–bound patterns (Bretskyan kind) can also be applied to other evolving systems: for example, cornets could be seen (Tëmkin and Eldredge Reference Tëmkin and Eldredge2007) as units of cultural evolution forming phylogenies, but they also can be seen as being parts of semi-isolated temporal and spatial evolving cultural units having both economic and genealogical components (Bretskyan-like hierarchy or “maker hierarchy” in the cultural realm [Eldredge Reference Eldredge, Prentiss, Kuijt and Chatters2009]). If viruses are a grade of organization—self-replicating coalitions of sequences that appeared many times by means of individuation from more complex cellular ancestors or their plasmids (as was inferred by [Kazlauskas et al. Reference Kazlauskas, Varsani, Koonin and Krupovic2019])—each separate independently originating monophyletic branch of “virus-level” evolutionary units could be understood as a separate hierarchy of the Linnaean kind.

The dual hierarchy view (first introduced by Eldredge and Salthe [Reference Eldredge, Salthe, Dawkins and Ridley1984]) accentuates the distinction of germ and soma (the Weismann barrier) or “functional” and “reproductive” in biological hierarchies. There is basic ontological truth to this generalization. On the other hand, the specialization of cells in organisms to function as reproducers of organisms or as servers of the organism's current function (with similar specializations at lower levels of cells and higher levels of superorganisms or eusocial colonies) is an evolved feature of subdivision of labor (Simpson Reference Simpson2011). Here, using the population genetics definition of fitness as a relative contribution of a given type to the next generation (or “probability of reproductive success”), we can think of somatic cells (non-reproducing individuals at the cellular level) as cells having zero long-term fitness (w ≈ 0), and germ cells as having larger than zero fitness (w > 0); therefore, we can think about them as heterogenous dimorphic populations with respect to fitness. The same principle applies to animal superorganisms, such as eusocial hymenopterans. On the other hand, somatic cells, as well as workers in hymenopteran colonies (Bourke Reference Bourke2011), have a potential to “darwinize back” (using the terminology of Godfrey-Smith [Reference Godfrey-Smith2009])—in the first case, cancers emerge, in the second case, disintegration of the colony through parasitism (which is dynamically equivalent to cancer) ensues, workers in a sense become parasites of their own colony. Darwinization is understood here (following Godfrey-Smith [Reference Godfrey-Smith2009]) as reacquisition of autonomy in reproduction of individuals at a given level, and the transition from epigenetic reproduction control (Simpson Reference Simpson2011), which is dependent on organismal (as in the case of cells) or colony-wise (as in the case of superorganisms) context, to genetically controlled reproduction; in other words, there is a reacquisition of capacity to participate in dynamics driven by mutation, selection, and statistical drift. We should emphasize that the same conceptualization can be repeated at many levels of economic, genealogical, and Bretskyan hierarchies: eukaryotic cells can be conceptualized as economic powerhouses, as coevolving heterogenous/polyphyletic collectives, or as single replicating units, as can organisms (unicellular and multicellular), colonies, and so on. The fluidity in exact definition and fuzziness of boundaries of the paradigmatic individuals—the organisms—which is very apparent in plants (White Reference White1979) or fungi (Lakovic and Rillig Reference Lakovic and Rillig2022), strengthens the case that in principle, we can treat most spatially and temporally contiguous (if not all?) biological entities as Bretskyan units of differing inclusivity levels. Consequently, we can conceptualize cases of emergence of cancers and their competitive dynamics with native cells inside organisms as also being of the Bretskyan unit type. There have been numerous attempts at understanding the sub-microevolutionary dynamics of cancers in ecological paradigms (Walther et al. Reference Walther, Hiley, Shibata, Swanton, Turner and Maley2015; DeGregori and Eldredge Reference DeGregori and Eldredge2020)—coevolution of spontaneously evolved (darwinized) unicellular parasites usually results in the death of an organismic Bretskyan unit. But sometimes a given cancer becomes successful and evolves into a transmissible form (Murchison et al. Reference Murchison, Tovar, Hsu, Bender, Kheradpour, Rebbeck, Obendorf, Conlan, Bahlo and Blizzard2010) or acquires a completely new niche (Van Valen and Maiorana Reference Van Valen and Maiorana1991), while in other cases, neoplasm cells become able to “steal” organelles from their hosts (Rebbeck et al. Reference Rebbeck, Leroi and Burt2011)—the inside-community large-scale gene transfer resulting in long-term survival of a parasite by means of escaping the boundaries of its ancestral multicellular organism.

A major difference between the Linnaean and the Vernadskyan hierarchies is the problem of formulating heritability in ecosystems. For evolution to happen, the progeny should resemble the parents (Lewontin Reference Lewontin1970), or more generally there should be significant temporal correlation between the states. Ecosystems are constituted of two components—the living populations, which are characterized by genetic heritability, and the nonliving one, which is harder to conceptualize with respect to heritability. Here we argue that both of these components have heritability, and the heritability of a nonliving system is well recognized in the geological time scale hierarchy—geologists routinely trace fission and the coalescence of continents and terrains at all time scales (Table 1). What enables the determination of this geodynamic phylogeny of geological bodies? The essence of heritability is the inertia or tendency for persistence—it can be achieved by a number of ways. In biology, the key role is played by replication and the adaptive error correction. In geology, persistence emerges in a number of forms, for example, (1) mechanical/thermal/chemical durability, such as in the case of zircons, which can persist as components of rock for eons; (2) dynamic stability enabled by buoyancy, such as in the case of low subduction potential of continental crust, which also has potential for indefinite persistence; (3) and the size itself—due to locality of interactions, larger structures are harder to modify (e.g., by erosion, which works on surfaces). All these factors determine time scales of persistence and therefore heritability of geological features. Persistence-based natural selection is apparently a universal feature across all kinds of systems (Lieberman and Melott Reference Lieberman and Melott2013). The persistence of geological structures causes the persistence of both living and nonliving components of ecosystems while explaining the nonliving matter part of the ecosystem heritability equation. Therefore Bretskyan units are characterized by two equally important kinds of heritability: (1) replication of genealogical units of all ranks and (2) geologically caused persistence of the physical space and material properties of matter they occupy and are immersed in.

We should mention that in the original definition of the dual biological hierarchy theory (Eldredge and Salthe Reference Eldredge, Salthe, Dawkins and Ridley1984; Eldredge Reference Eldredge1985, Reference Eldredge1986), organisms were considered to belong to both genealogical and economic hierarchies, thus being the crucial connection between functional and genetic multilevel processes. In this sense, by recognizing the special status of the organism level, the dual hierarchy approach already suggested the need for the intermediary hierarchy of eco-genealogical entities. Our present work reveals that the ecosystemic interactions and genealogical significance to varying degrees can be found in many spatially and temporally bounded biological structures. The Bretskyan hierarchy simply explicitly acknowledges all these properties and frames their relative significance in biological ontology as a function of size and duration of the studied entities located in real space, time, and material matrix (Table 1).

Hierarchy of Physical Barriers, Forced Individuality of Large-Scale Bretskyan Units, the Concept of Multiscale Allopatry, and the Kinds of Bretskyan Units

Isolating mechanisms are enablers of the emergence of individuality of higher-level evolutionary units: speciation is a paradigmatic case in which populations split due to internal incompatibilities or due to external barriers that prevent gene flow, which can result in quick extirpation of incipient species (Allmon et al. Reference Allmon, Morris, McKinney, McKinney and Drake1998; Stanley Reference Stanley2008), or alternatively in often punctuated change and subsequent stasis of a lineage (Stanley Reference Stanley1979; Futuyma Reference Futuyma1987; Eldredge et al. Reference Eldredge, Thompson, Brakefield, Gavrilets, Jablonski, Jackson, Lenski, Lieberman, McPeek and Miller2005; Hunt Reference Hunt2007). Another important set of aspects includes spatial structures of populations and locality of their interactions (Dieckmann et al. Reference Dieckmann, Law and Metz2000) and the multiscale heterogeneity of environments (Jablonski Reference Jablonski2008a). For example, even in the earliest stages of evolution before the emergence of cells, the locality of self-replicating chemical reactions on 2D surfaces should have significantly contributed to the emergence of organism-like reproducers (Szathmáry and Smith Reference Szathmáry and Smith1997; Czárán and Szathmáry Reference Czárán, Szathmáry, Dieckmann, Law and Metz2000; Szathmáry Reference Szathmáry2015). The spatial separation itself can be understood as a form of quantitative/graded isolation—the interactions are exponentially more plausible in proximity than with distance (Edelstein-Keshet Reference Edelstein-Keshet2005). Spatial separation individuates the constituent populations of a species, while forming clinal variation (Gould and Johnston Reference Gould and Johnston1972; Yablokov Reference Yablokov1986) and complex spatiotemporal patterns of microevolution (Lieberman et al. Reference Lieberman, Brett and Eldredge1995; Brombacher et al. Reference Brombacher, Wilson, Bailey and Ezard2023).

Even though the raw separation in space is an important factor in isolation, life-forms evolved and were subjected to a range of more effective isolating mechanisms that played a role in generating and sustaining evolutionary individuality. The acquisition of cellular structures of autonomous life-forms—the lipid membrane, separating the internal molecular network from most of the external genetic material and metabolically active molecules, was a key step in organizing replicating molecules into cooperative evolutionary units (Maynard Smith and Szathmáry Reference Maynard Smith and Szathmáry1998). The recognition of members of your own colony is an essential integrating mechanism in maintaining superorganismal individuality (Wilson and Hölldobler Reference Wilson and Hölldobler2005; Nowak et al. Reference Nowak, Tarnita and Wilson2010). The strongest cases of in-group and out-group recognition can be found in the most integrated individuals—for example, the CRISPR-Cas immune system in archaea and bacteria at the single-cell level, defending the cell from foreign genetic parasites/influences (Koonin and Makarova Reference Koonin and Makarova2019), and the more familiar but far more complex immunity and cell recognition mechanisms inside multicellular organisms (Gilbert and Tauber Reference Gilbert and Tauber2016). Also in the paradigmatic cases of individuality at the species level, active reproductive isolating mechanisms are crucial in maintaining individuality of species by means of Paterson's “species mate recognition system” (SMRS; Vrba Reference Vrba1980; Vrba and Eldredge Reference Vrba and Eldredge1984; Paterson Reference Paterson and Vrba1985).

Allopatric speciation and its para- and peripatric variations (Gavrilets Reference Gavrilets2004)—a kind of speciation due to geographic separation of ranges—is thought to be the major factor responsible for the individuation of populations and, consequently, generation of new species (Dobzhansky Reference Dobzhansky1937, Reference Dobzhansky1940; Mayr Reference Mayr1940; Eldredge Reference Eldredge1971; Eldredge and Gould Reference Eldredge, Gould and Schopf1972; Allmon and Sampson Reference Allmon, Sampson, Allmon and Yacobucci2016; Hernández-Hernández et al. Reference Hernández-Hernández, Miller, Román-Palacios and Wiens2021; Anderson and Weir Reference Anderson and Weir2022). There is a growing amount of evidence that the Earth's geomorphic structures play the major role in generating biodiversity (Rahbek et al. Reference Rahbek, Borregaard, Antonelli, Colwell, Holt, Nogues-Bravo, Rasmussen, Richardson, Rosing and Whittaker2019), following the idea that the complexity of the planet is the main modulator of diversity and stability of habitats (Archibald et al. Reference Archibald, Bossert, Greenwood and Farrell2010) and isolation opportunities (Cracraft Reference Cracraft1982). Therefore, speciation and extinction processes are intrinsically linked to the structure and dynamics of the Earth and its climate (Stanley Reference Stanley1990; Vrba Reference Vrba1993). The same recurring principle of environmental heterogeneity in maintaining and generating diversity apparently scales down to even short-term ecological processes, which helps to explain the so-called Hutchinson's paradox of the plankton (Hutchinson Reference Hutchinson1961; Descamps-Julien and Gonzalez Reference Descamps-Julien and Gonzalez2005). If the system is scaling and multifractal, such as topography or the geology of the planet (Gagnon et al. Reference Gagnon, Lovejoy and Schertzer2006), then we should expect non-uniformity (much more extreme than Gaussian) at all spatial and temporal scales—for example, some regions can be topographically very rough, while others have very low relative variability in any of the measured parameters, therefore creating very quickly changing biomes, as in mountain ranges, or the opposite, the very monotonic features such as the Great Eurasian Steppe or the Amazon rain forest. And even here heterogeneity plays a role, such that the former biome is desolate and species-poor, while the latter is species-rich. The complexity in macroevolution is usually meant as a single-level quantity—a scalar—although it is (in the case of scaling systems such as Earth) a multiscale property, which can be represented by a spectrum or structure function of the scale. Here we use the definition in a latter meaning—the Earth is heterogenous and diverse at multiple resolutions.

Speciation and extinction dynamics exhibit significant coherence between clades—X-genealogical turnovers—at both the regional (Vrba Reference Vrba1985, Reference Vrba1992) and global scales (Sepkoski Reference Sepkoski1981, Reference Sepkoski1984; Bambach Reference Bambach2006; Alroy Reference Alroy2010; Cuthill et al. Reference Cuthill, Guttenberg and Budd2020). This pattern signals the importance of common external drivers of biodiversity dynamics usually tied to changes in climate states, sea level, extraterrestrial causes, and tectonics (Valentine and Moores Reference Valentine and Moores1970; Alvarez et al. Reference Alvarez, Alvarez, Asaro and Michel1980; Finnegan et al. Reference Finnegan, Heim, Peters and Fischer2012; Stigall et al. Reference Stigall, Edwards, Freeman and Rasmussen2019; Mathes et al. Reference Mathes, van Dijk, Kiessling and Steinbauer2021; Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022). Some of the best examples relevant for the establishment of the picture of hierarchical emergence and persistence of geobiomes come from studies of Paleozoic marine faunas or of so-called coordinated stasis of multiple lineages and community patterns (Brett and Baird Reference Brett, Baird, Erwin and Anstey1995; Brett et al. Reference Brett, Ivany and Schopf1996).The majority of lineages formed recurring communities, and the taxonomic and ecological turnover was concentrated within relatively short periods of time (Ivany et al. Reference Ivany, Brett, Wall, Wall and Handley2009). The application of dynamical systems tools to Silurian Period faunas revealed that the patterns of species dominance were modulated by oceanic circulation and sea-level states and sharp transitions between them (Spiridonov et al. Reference Spiridonov, Brazauskas and Radzevičius2015, Reference Spiridonov, Brazauskas and Radzevičius2016, Reference Spiridonov, Kaminskas, Brazauskas and Radzevičius2017a,Reference Spiridonov, Stankevič, Gečas, Šilinskas, Brazauskas, Meidla, Ainsaar, Musteikis and Radzevičiusb, Reference Spiridonov, Stankevič, Gečas, Brazauskas, Kaminskas, Musteikis and Kaveckas2020b; Venckutė-Aleksienė et al. Reference Venckutė-Aleksienė, Radzevičius and Spiridonov2016; Spiridonov Reference Spiridonov2017; Crampton et al. Reference Crampton, Meyers, Cooper, Sadler, Foote and Harte2018; Whittingham et al. Reference Whittingham, Radzevičius and Spiridonov2020), with a pattern showing coherence in space (Spiridonov et al. Reference Spiridonov, Samsonė, Brazauskas, Stankevič, Meidla, Ainsaar and Radzevičius2020a) and repeating itself at several time scales (Rinkevičiūtė et al. Reference Rinkevičiūtė, Stankevič, Radzevičius, Meidla, Garbaras and Spiridonov2022). The spatial evolutionary processes are modulated by sea-level change, climate, and tectonic perturbations—these factors play a crucial role in modulating diversity at large temporal and spatial scales (Lieberman and Eldredge Reference Lieberman and Eldredge1996; Lieberman Reference Lieberman2003; Radzevičius et al. Reference Radzevičius, Spiridonov, Brazauskas, Dankina, Rimkus, Bičkauskas, Kaminskas, Meidla and Ainsaar2016; Stigall et al. Reference Stigall, Bauer, Lam and Wright2017).

From such paleobiogeographic and macroevolutionary examples discussed earlier, we implicitly can see the hints of a general pattern suggesting what we call hierarchical or multilevel allopatry, which should be directly related to the nature of geophysical processes and the hierarchical structure of the planet. The classical examples of coordinated stasis and species factories show the importance of sea-level and climate perturbations, which generate repeated episodes of isolation (= individuation of species and whole local geobiomes) and dispersal of newly generated taxa (Stigall et al. Reference Stigall, Bauer, Lam and Wright2017). Epeiric seas are shallow and thus susceptible to relatively small-magnitude sea-level/climate changes. Such changes could be generated by relatively fast “grand” Milankovitch cycles with period lengths of 1.3–2.6 Myr (Van Dam et al. Reference Van Dam, Aziz, Sierra, Hilgen, van den Hoek Ostende, Lourens, Mein, Van Der Meulen and Pelaez-Campomanes2006; Crampton et al. Reference Crampton, Meyers, Cooper, Sadler, Foote and Harte2018), and possiby by even faster precession cycles at scales ≤ 25 kyr (Crampton et al. Reference Crampton, Cooper, Foote and Sadler2020). These time scales of millions of years are commensurable with durations of separate species and do not allow a deep repeated multigenerational (measuring in species generations) coevolution. Therefore, the scale of species-level multigenerational coevolution should correspond to the continental or oceanic spatial scales and the time scales of tens of millions of years. The concept of a continent and its distinction from an island is an arbitrary one. The same goes for oceans and seas. Tectonic plates exhibit power law scaling or a “scale-free” pattern, which means plates do not have characteristic sizes (Sornette and Pisarenko Reference Sornette and Pisarenko2003). The same applies to the fragmentation index of continents—meaning that we should expect fragmentation to happen on all time scales (Spiridonov et al. Reference Spiridonov, Balakauskas and Lovejoy2022). Due to positive scaling of fragmentation indices with increasing time scales, we find the largest differences between continental configurations at the longest time scales. Therefore, larger spatial structures on average have more formidable barriers, which also in turn last proportionately longer. The connection of large geographic sizes and long durations, which are far longer than typical species durations (typically of millions of years [Raup Reference Raup1991; Cooper and Sadler Reference Cooper and Sadler2010; Žliobaitė and Fortelius Reference Žliobaitė and Fortelius2022]), implies the significance of these planetary structures in keeping forced individuality of geobiomes (Fig. 3).

Figure 3. Allopatric space and the hierarchies of allopatries. Species-level allopatry requires relatively short time scales—higher end estimated from the duration of species, which is on the order of several millions of years, although it can potentially occur in a large space, depending on the dispersal abilities of taxa. Here the species allopatric space is a subset of Bretskyan units’ allopatric space—where whole biotas (“provinces”) diverge and experience coherent coevolution and integration in externally isolated or semi-isolated ecosystem units. The size of a province need not have a particular scale and can be a fragment of a biota of any size larger than a certain threshold determined by the characteristics of the planet's geodynamics, which determines the intensity, duration, and other physical parameters crucial for biotic isolation and persistence of the isolation for a sufficient amount of time.

Bretskyan units at larger scales—although more fuzzy than, for example, species, which also have a property of fuzzy/non–sharply separated sets (Van Valen Reference Van Valen1978)—are true ontological individuals, because they are spatiotemporally bounded, with origins, historical durations, and terminations. The apparent fuzziness of geobiomes stems in part from our scale-boundedness as observers. There is a growing appreciation of such geobiomic ontology (Ung and Buttiegieg Reference Ung and Buttiegieg2023). Geobiomic units constitute a rank-free scalar hierarchy in space and a network-like structure in time. As seen in time, the Bretskyan hierarchy resembles a braided river with anastomosing units emerging, merging together, separating, or changing hierarchical significance in space depending on the degree of isolation as barrier intensity and spatial extent fluctuates (Fig. 4). This flow of geobiomic individuals at long time scales achieves global mixing, when all previously separated units (through contacts with neighbors) exchange/percolate evolved biotas. It was estimated that these mixing or synchronizing processes led by plate tectonics, such as “Noah's Ark”/“docking events” of continents or islands (McKenna Reference McKenna1972), start to dominate diverging local and regional tendencies at time scales longer than 40 Myr (Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022, Reference Spiridonov and Lovejoy2023). Such mixing events, which work in synchronizing and homogenizing biota at long time scales, are now well documented. One example is the connection event where trilobites of genera Kootenia and Bailiella from the previously distinct western Newfoundland “Pacific Province” and eastern Newfoundland “Atlantic Province” were found together showing the merger event and the origin of a new larger fused Bretskyan paleobiogeographic unit during the middle Cambrian (Kay and Eldredge Reference Kay and Eldredge1968). The great American biotic interchange (GABI) is another great example (Marshall et al. Reference Marshall, Webb, Sepkoski and Raup1982) that will be discussed later in the paper.

Figure 4. Conceptual figure of spatiotemporal dynamics of large Bretskyan units (geobiomes) at two levels. On the left side of the figure at Time 1, we start with two large-scale geobiomes (bioprovinces) that are also composed of many smaller persisting, splitting, fusing, and disappearing geobiomes. At Time 2, we have a merger event of two provinces that become one. At Time 3, we have the fission of province C into three new provinces D, E, and F. During this time period, province E goes extinct (e.g., during disappearance of an isolated oceanic basin during collision of continents). On the right-hand side of the figure, we can see the strict (although fuzzy to different levels) scalar hierarchical structure of bioprovinces at any given time. Nested patterns show biotic similarities between bioprovinces.

Larger and more formidable geological structures serve as barriers for dispersal of biota for longer periods of time. Thus, we should expect a hierarchy of coevolutionary processes in space and time that integrates geobiomes into hierarchical symbiotic collectives—larger geobiomes with more formidable dispersal barriers will experience proportionately more integration with more constituent taxa (Fig. 5). It shows allopatry as a scale-free concept—the same repeating principles of isolation and contingent divergence are applicable not only at the species level but also for much larger and far more complex entities. The dynamics of geobiomes, though, have a significant difference with allopatric speciation, because Bretskyan unit isolation is inevitably followed by coalescence (yet hybridization between closely related sister taxa is not at all uncommon). In this way, the dynamics of Bretskyan units or geobiomes topologically can be described as representing a process akin to hierarchical “sexual reproduction.” This concept perhaps sounds rather exotic, but it can be described in a clear way—in a multilevel process of fragmentation and merging of geobiomes, smaller geobiomes diverge and merge with other smaller geobiomes, thus sharing evolved features and coevolved relations with each other at short time scales; while at the same time these collectives of isolating and coalescing geobiomes constitute still larger geobiomes that repeat the process at still longer time scales and larger space scales, and so to infinitum or in the real world to the limit of the size of a planet with active geodynamics. Here individuality of geobiomes emerges and disappears in a scaling manner: there is a continuum and proportionality of their sizes and durations. The structures on our planet are characterized by multifractal variability (Spiridonov et al. Reference Spiridonov, Balakauskas and Lovejoy2022), that is, they are scaling, self-similar, and show strong variability at allscales.Footnote *

Figure 5. Spatiotemporal evolution of isolating barriers, for the ease of understanding presented in 1D space (x-axis) and 1D time (y-axis), as a function of dispersal abilities of clades (z-axis). Weak barriers have shorter duration and spatial extent (due to positive scaling of durations of geological bodies/barriers in relation to their size). Smaller barriers are sufficiently long-lasting to work as species-generating barriers (species-level allopatric barriers), and larger and longer-lasting barriers work as barriers for large geobiomes (larger-scale Bretskyan units). These very significant barriers (e.g., oceans for terrestrial non-volant species) isolate in tandem many different clades into integrative units. Note that many generations of lower-level barriers (= many generations of species) are present; this makes micro- and macroevolutionary coadaptation all the more likely for the whole biota.

A similar observation was made by Charles Darwin almost 200 years ago, when he, during the voyage on the Beagle, compared floras of two sides of the Andean mountains (Darwin Reference Darwin1839: p. 399): “I was very much struck with the marked difference between the vegetation of these eastern valleys and that of the opposite side: yet the climate, as well as the kind of soil, is nearly identical, and the difference of longitude very trifling.” Later he continues, implicitly linking the contrasts in floras and faunas to the great size, height, and consequently long duration of the Andes, which he understood as an epochal barrier: “This fact is in perfect accordance with the geological history of the Andes; for these mountains have existed as a great barrier, since a period so remote that whole races of animals must subsequently have perished from the face of the earth.” Although the concept of scaling did not exist when Darwin wrote these words, the implicit link between the dimensions of geological structures and their geological longevity, and as a consequence their significance in producing biological diversity, was hinted at in these earliest observations of modern science.

Barriers to the spread of holobionts exist not only in geography, but also in time—sharp environmental changes of short duration and large magnitude cause the demise and consequent turnover of the global biota (Raup and Sepkoski Reference Raup and Sepkoski1982; Bambach Reference Bambach2006). Similarly, as in the case of topography, mass extinctions, and presumably their causes, are the extreme right tails of continuous distributions (Marshall Reference Marshall2023). The majority of extinctions and originations at the global scale are probably caused by many such impulse-like perturbations of varying sizes (Foote Reference Foote2005). And indeed the statistics of global originations and extinctions at scales of millions of years and also populational processes at shorter time scales of tens of thousands of years show the signature of multiplicative multifractal variability (Plotnick and Sepkoski Reference Plotnick and Sepkoski2001; Yacobucci Reference Yacobucci2005; Spiridonov et al. Reference Spiridonov, Brazauskas and Radzevičius2016; Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022), which implies sharp, intermittent (typified by extremes), and hierarchically modulated changes in time that, as seen from the Bretskyan units perspective, create sharp temporal boundaries in the composition of the biota (up to time scales of tens of thousands of years [Crampton et al. Reference Crampton, Cooper, Foote and Sadler2020]). These findings support the idea that the global biota—Gaia—experiences and transitions through sharp separation events from its previous compositional states. The marine animal “communities” originally studied and distinguished at time scales of geological periods (tens of millions of years) by Bretsky (Reference Bretsky1968) were essentially representations of Gaia. The geological time scale, which is itself hierarchical, reflects the levels of severity of the perturbations, or the degree of transformativeness more generally, in the evolution of Gaia.

It was suggested that, although of a very different nature, mass extinctions represent a kind of “major transition in evolution” (Currie Reference Currie2019) similar in importance to traditional “major transitions” instances of purely creative eukaryogenesis or the origins of the human language (Maynard Smith and Szathmáry Reference Maynard Smith and Szathmáry1998; Bickerton Reference Bickerton2009; Calcott and Sterelny Reference Calcott and Sterelny2011). In the light of the presented theory, mass extinctions indeed could quite literally be seen as the major transitions in evolutionary individuality between compositionally and functionally differing states of Gaia. These temporal barriers are, of course, different from spatial barriers due to the presence of “time's arrow” or the irreversibility of the time dimension, although similar in other aspects. As biotas can be starkly distinct on continents isolated by strong environmental barriers in space, similarly biotas can be starkly distinct by means of their sharp environmental isolation in time—they differ in ecological functioning (Wagner et al. Reference Wagner, Kosnik and Lidgard2006; Aberhan and Kiessling Reference Aberhan and Kiessling2015) and the states of macroevolutionary dynamics (Miller and Foote Reference Miller and Foote2003; Jablonski Reference Jablonski2005; Brayard et al. Reference Brayard, Escarguel, Bucher, Monnet, Brühwiler, Goudemand, Galfetti and Guex2009; Alroy Reference Alroy2010), including in geographic space itself (Krug et al. Reference Krug, Jablonski and Valentine2009). Therefore, Gaias (i.e., differing instantiations of Gaia in different times) as the largest units in the Bretskyan hierarchy, separated by mass extinctions, show significant structural and functional differences and thus can be seen as sequential individuals with ancestor–descendant relations occupying approximately the same space. The same arguments work for mass originations and mixtures of high rate of extinction and origination events (Marshall Reference Marshall2006; Cuthill et al. Reference Cuthill, Guttenberg and Budd2020).

In a similar vein, the International Chronostratigraphic Chart represents a hierarchy of changes in the composition and functioning of global Bretskyan units; regional and local stratigraphic scales, especially when based on biostratigraphy, should to a significant extent represent transitions in states of smaller in areal extent, regional and subregional Bretskyan units. Therefore, geological time scales reflect ontologies of change in individuality of geobiomes. One such example is the NALMA (North American land mammal ages) time scale subdividing Late Cretaceous to Cenozoic terrestrial sequences of North America based on the compositional changes in mammalian assemblages at the continental scale (Flynn et al. Reference Flynn, MacFadden and McKenna1984). The apparent diachrony of

first and last appearance events in different regions of North America (Alroy Reference Alroy1998b) can be understood as a reflection of the presence of finer-scale geobiomic substructures in the North American geobiome.

We can conceptualize some important evolutionary events as representing additional major transitions of individuality—the major transitions in kinds of sub-Gaian Bretskyan units—examples of which are terrestrialization and the expansion of life into the open ocean or Earth's deep interior (some of the “megatrajectories” of Knoll and Bambach [Reference Knoll and Bambach2000]). For example, terrestrialization at the beginning of the Phanerozoic led to formation of a “parallel biogeography” weakly interacting with its marine geobiomic ancestor (interactions inside the realm are much more important than those outside it). The new characteristics of previously unoccupied environments impose new sets of physical and geochemical constraints and feedbacks on these constraints (Kenrick et al. Reference Kenrick, Wellman, Schneider and Edgecombe2012; Lenton et al. Reference Lenton, Dahl, Daines, Mills, Ozaki, Saltzman and Porada2016). Such transitions result in qualitatively different drivers, rates, and scales for the coevolution inside Bretskyan units of new type due to radically different ranges of physical parameters of occupied environment. The largest difference between the marine biosphere and the terrestrial one is the topology—the marine biosphere at first approximation is always contiguous. It topologically resembles a plane with holes. On the other hand, continents are essentially islands isolated by hardly penetrable marine barriers. In the future, percolation network approaches (e.g., Palamara et al. Reference Palamara, Rozenfeld, de Santana, Klecka, Riera, Eguíluz and Melián2023) will give us a comprehensive understanding of the nature of differences in temporal topological connectivity of oceanic and terrestrial geobiomes at geological time scales. The example of emergence of qualitatively different geobiomes by means of break-through (escape) into qualitatively new environments of a number of clades, and their later ecological integration, shows that there is a qualitative difference in modes of major transitions at different levels of the Bretskyan hierarchy. Traditional (small inclusivity scale) major transitions (e.g., eukaryogenesis) are dominated by the pure fusion of lower-level polyphyletic units, while the largest scale transition and the origins of the quasi-autonomous terrestrial biosphere is a case of fission of an already formed (and eons old) marine sub-Gaian unit, resembling a higher-order “cladogenetic” event. At the same time, it is also an instance of the growth or expansion (Van Valen Reference Van Valen1989) and differentiation of Gaia.

Time Scales of Coevolution inside Bretskyan Units, Fitness Landscapes of Geobiomes, and Spatiotemporality of Contingency, Adaptation, and Progress in Macroevolution

Here we present evidence that the structure of a biota is not an epiphenomenon, as it is sometimes presented (Hoffman Reference Hoffman1979; Bambach and Bennington Reference Bambach, Bennington, Jablonski, Erwin and Lipps1996), but an essential ontological feature of biology, necessary for the causal understanding of evolution. The major obstacle in understanding evolution of communities stems from previous lack of a coherent conceptual framework that could explain the mechanisms of individuation of communities at the whole range of time and space scales. The spatially and temporally defined hierarchy of Bretskyan units is the closest representation of the living world we know of. All interactions (metabolic and genealogical) happen in ecosystems, which are separated to varying degrees from other such systems. The mechanisms and the typical magnitudes of separation and thus individuation vary as a function of scale. At larger scale, the barriers individuating biota have a geological (in a broad sense of the word) origin, and their space–time characteristics are determined by geodynamics. Therefore, the knowledge of the statistical time and space scaling characteristics of geological phenomena provides a unified theory on the patterns of individuation, fusion, and scales of coevolution inside and among geobiomes.

The evolution of holobionts at lower spatial and temporal levels of the Bretskyan hierarchy is a well-developed field (e.g., Gilbert and Tauber Reference Gilbert and Tauber2016). Here we concentrate on the much less understood large spatial and long temporal scales—regions, provinces, and Gaia herself—which respectively can last from thousands to up to billions of years. The main concept here is coevolution and, more precisely, the geobiomic depth of coevolution, which is directly linked to the physical characteristics of Bretskyan units. The next concept of importance is so-called fitness landscapes of geobiomes—they are analogues of classical fitness or adaptive landscapes (Gavrilets Reference Gavrilets2004; McGhee Reference McGhee2006). Characters of geobiomes are statistical organismal-level, populational-level, species-level, or functional characteristics of constituent taxa. Fitness is here understood as the probability of survival/proliferation of evolved characteristics at a given geobiomic level.

Coevolution implies the reciprocal/mutual evolutionary causation of change in two or more groups of individuals (Levins and Lewontin Reference Levins and Lewontin1987; Thompson Reference Thompson2013). The basic prerequisite for the emergence of coevolution is proximity or overlap of ranges of evolutionary individuals in time and space. Speciation, which is the result of splitting and isolation of one lineage into at least two, entails the separation of microevolutionary dynamics inside these lineages. After the split, lineages follow contingent evolutionary paths while having little genetic exchange of newly acquired adaptations and other traits. Isolated lineages start contingently adapting to the surrounding biotic and abiotic environments, including experiencing coevolution inside the lineage itself (Prum Reference Prum2017). We can add that even simple models of evolution, with contingent separation and isolation events, result in wide exploration of the morphospace (Niklas Reference Niklas2004; Solé Reference Solé2022). This is an example of the evolutionary “butterfly effect,” a metaphor based on the short story “A Sound of Thunder” written by Ray Bradbury in 1952. This metaphor says that even small differences in initial conditions of evolving lineages can subsequently cause significant divergence in traits or sequences of selective events (Niklas Reference Niklas2004).

The critical role in defining, and making species coherent entities is played by the SMRS (Paterson Reference Paterson and Vrba1985); species individuality is maintained by the “glue” of coevolution, which includes gene exchange and a mutual selection and other “economic” effects of organisms in a Red Queen evolutionary mode (Van Valen Reference Van Valen1973). Therefore, the individuality of species mechanistically is maintained by more intense interaction between organisms of the same group than with other such groups. These interactions are only possible if the lower-level individuals in question occupy the same geographic space, and if they reside in the same time interval. The same arguments work for more inclusive spatiotemporal eco-genealogical individuals that constitute the Bretskyan hierarchy.

The genome of a population subject to natural selection maximizes so-called Fisher information (about the environment), here understood as an information measure that reflects an unknown (environmental) probability distribution (Frank Reference Frank2009)—for example, the presence of thick fur reflects cold climate, fat stores reflect unevenness of food supply, and so on. It was suggested that evolution represents a multilevel learning process (Vanchurin et al. Reference Vanchurin, Wolf, Katsnelson and Koonin2022); here we show how this learning is organized in the real-world hierarchies. It should be noted that the information maximized within a population of evolutionary individuals in their genomes is only about “effective” factors that contribute to the survival and expansion of lineages. To understand the importance of maximization of Fisher information in evolutionary dynamics of geobiomes, we should remember that the effective environment (space in which the holobiont evolves) is characterized by size and duration. Therefore we should expect scaling in acquired information about the Earth system by any holobiont or geobiome. The distribution, size, and duration of barriers determine the geodiversity. Geobiomes of larger sizes and durations, and occupying more diverse geological settings, are expected to be more “learned.” Because only part of the environment is a physical one, and the other is the biological one (Levins and Lewontin Reference Levins and Lewontin1987), we should also expect acquisition of more Fisher information about biota in larger areas. This follows from species–area curves (Rosenzweig Reference Rosenzweig1995) and the law of mass action (the number of possible interactions is proportional to diversity). More species means more ecological interactions from which all species could mutually “learn.” In a larger geobiome, there will be greater coevolutionary depth than in smaller, less diverse, and shorter-lived geobiomes. The rates of biological progress—competitiveness level in as wide a range of biotic and abiotic environments as possible—are directly related to hierarchical dynamics of planetary geology in a broad sense.

Physical isolation works as an integrative external force for isolated monophyletic clades, because it causally divorces some members of the clade from others (Fig. 6). The coevolution of species of a single clade inside a large-scale geobiome (e.g., bioprovince) should structurally resemble the “game against a version of yourself.” Species belonging to the same family or genus show higher similarity in all their characteristics between each other than to other less related species due to phylogenetic constraints (Raup and Gould Reference Raup and Gould1974; Gould Reference Gould2002). Therefore, species are more likely to be exploiting a similar resource base and acting as (from the ecosystem point of view) approximate functional equivalents of each other. This high initial similarity of competing interactors results in stiff competition and escalation that can either result in extinction or evolution of specialized morphologies due to so-called character displacement (Pfennig and Pfennig Reference Pfennig and Pfennig2012). The described dynamics are not unlike those in artificial intelligence games, which self-master by playing against copies of themselves without any external supervision, as was implemented in the case of AlphaGo Zero, which achieved superhuman mastery of the ancient game of Go (Silver et al. Reference Silver, Schrittwieser, Simonyan, Antonoglou, Huang, Guez and Hubert2017). Again, the greatest depth of coevolution and the greatest progress are expected in the largest geobiomes/Bretskyan units.

Figure 6. Clades separated for long periods of times will experience multiple rounds of competitive and other biotic interaction–mediated evolution inside the bioregions. This will increase the coadaptability of species in comparison to randomly drawn members of a clade from other regions.

Integration by means of multiscale coevolution inside geobiomes is also an expected pattern for any set of taxa occupying the same bioregion (Fig. 7). Because selection and cross-optimization between different functional units should be more powerful in larger, longer-lasting, and more diverse bioregions, the geobiomes of these regions should be of higher resilience to external biotic and abiotic perturbations. Because it is expected that selection on ecosystems results in increased stability of ecosystems at all scales (Borrelli et al. Reference Borrelli, Allesina, Amarasekare, Arditi, Chase, Damuth and Holt2015), we should expect a greater level of stability to perturbations to evolve in larger and longer-lasting Bretskyan units. The largest—the Gaian-level Bretskyan unit—is potentially eternal (up to the astrophysical limits); therefore, in the long run, it should achieve the maximal possible stability level. Empirical studies apparently confirm this pattern (Van Valen Reference Van Valen1984; Foote Reference Foote2000; Miller and Foote Reference Miller and Foote2003; Markov Reference Markov2009; Lieberman and Melott Reference Lieberman and Melott2013).

Figure 7. What is essentially the same process as in Fig. 6 works for mixtures of clades. Different regions separated from one another can receive invasions of different unrelated clades that will interact with the descendants of a given clade. This will create contingently different coevolutionary complexes (“splendid isolations” of George Gaylord Simpson's “South Americas” [Reference Simpson1983]). Here we have strong compartments that serve as robust individuating boundaries, not unlike reproductive isolation in the case of sexually reproducing species. Species evolution (selection and sorting) inside Bretskyan units should work analogous to natural selection inside species, but here we will have directional changes in the distribution of features of whole geobiomes—analogous to changes in average phenotypes of individual organisms inside their respective species, as happens at shorter time scales.

Many aspects of evolution inside Bretskyan units can be conceptualized and understood using a generalization of fitness landscapes for geobiomes. In the case of geobiomes, instead of characters specific for Linnaean units such as organisms or species, we use characters that can describe geobiomes and holobionts in a meaningful way. These geobiome-scale phenotype variables can be, for example, average herbivore or average carnivore sizes in large vertebrate guilds (Fig. 8), variables describing ecological communities. These are essentially axes of ecospace (Novack-Gottshall Reference Novack-Gottshall2016a,Reference Novack-Gottshallb) with quantitative gradations and fitness values assigned to them. The measurement of fitness in geobiomic fitness landscapes is a scaling concept, as are geobiomes themselves: different levels of inclusivity of traits are needed as we compare geobiomes of different scales. The many measurements of such combinations over long time periods should produce an estimate of geobiomic fitness landscape and reveal the directionality in evolution. A fine example of such an approach is the study that determined that Cenozoic North American mammals (or the mammalian component of the North American Cenozoic geobiome) have two body mass attractors (Alroy Reference Alroy1998a).

Figure 8. The concept of the fitness landscapes of large-scale Bretskyan units—geobiomes. In the case of grand spatial scales and grand geographic barriers that last tens to hundreds of millions of years, we are entering into the arena of grand evolutionary contingency. Sets of distinct taxa that populate bioregions after major extinction or origination events or as a result of prolonged isolation create very distinct geobiome-level fitness landscapes, attracting evolutionary change to highly distinct sets of parameter values. Shown here are two hypothetical examples: Bretskyan unit 1 could be imagined as representing the case of non-avian dinosaurs. Their anatomic characteristics, such as bone pneumaticity and egg laying, enabled an r strategy of population dynamics, which increases survivability of large animals in uncertain environments, so non-avian dinosaurs evolutionarily achieved much larger sizes than any terrestrial mammals (Sander and Clauss Reference Sander and Clauss2008; Sander et al. Reference Sander, Christian, Clauss, Fechner, Gee, Griebeler, Gunga, Hummel, Mallison and Perry2011; Botha et al. Reference Botha, Choiniere and Benson2022). All these factors create boundary conditions for the emergence of radically different ecosystem fitness landscapes than in mammal-dominated ecosystems (represented by Bretskyan unit 2). Dinosaurs would evolve toward much larger predators and prey, and mammals would have optima at smaller body sizes. This difference in the real world happened in time and not in space, but this does not change the main conclusion. If non-avian dinosaurs, by any chance, had survived on some highly isolated continent, we currently would have had two very different geobiome fitness landscapes—one for mammals and the other for non-avian dinosaurs. The closest we currently have in the real world is the case of Australia vs. the rest of the inhabited continents. Here placental and marsupial mammals are highly distinct due to the presence of different developmental constraints (and their nature and importance for evolution can also change with time [Salazar-Ciudad et al. Reference Salazar-Ciudad, Jernvall and Newman2003])—marsupials probably will never develop truly marine species or flying or hoofed species, because their embryos need to have functional grabbing arms in order to get to the mother's pouch (Sears Reference Sears2004); moreover, the developmental constrains in marsupials have much wider effects (Fabre et al. Reference Fabre, Dowling, Portela Miguez, Fernandez, Noirault and Goswami2021). The noted spatially contingent features constrain the ranges and directions of diversification inside the Bretskyan units and thus define their geobiome-level fitness landscapes.

One of the major themes of evolution is the role of chance and contingency in the history of life (Gould Reference Gould, Briggs and Crowther2001; Erwin Reference Erwin, Ramsey and Pence2016; Ramsey and Pence Reference Ramsey and Pence2016). In his book Wonderful Life, Stephen Jay Gould argued that the extinction of body plans that were represented in the Cambrian by low numbers of species should have had a major consequence for all subsequent evolution, as it closed some major windows of evolutionary radiations and enabled some others represented by surviving phyla (Gould Reference Gould1990). Whether or not the early extinctions in animal evolution had such a profound effect on macroevolution and ecology in the Phanerozoic, the temporal events arguably are of great importance in forming the composition and patterns of evolutionary dynamics (Jablonski Reference Jablonski1986; Raup Reference Raup1992, Reference Raup1994; Brayard et al. Reference Brayard, Escarguel, Bucher, Monnet, Brühwiler, Goudemand, Galfetti and Guex2009). The Bretskyan hierarchy worldview presents a generalization of evolutionary contingency in space–time.

Different states of a given biota can be separated by physical barriers of catastrophic conditions (mass extinctions events) in time, so too can parts of a biota be separated by physical barriers in space. The composition of clades before and after an extinction event can radically differ; also the composition of clades can differ in space. Depending on the composition of clades that are parts of geobiomes, any given geobiome could differ significantly in the shape of its fitness landscape, and therefore its vectors of evolutionary dynamics. For example, the fitness landscape of average masses of apex terrestrial vertebrate predators in the South American Miocene should have been very different from the fitness landscape of the North American geobiome. In the former case, South America was inhabited by members of four clades that filled that niche—metatherian sparassodont mammals, terrestrial crocodiles, phorusrhacids or “terror birds,” and the giant snakes Madtsoiidae (Prevosti et al. Reference Prevosti, Forasiepi and Zimicz2013); while in the latter case, all top predators were placental mammals. Many other examples of the role of the space–time contingency in evolution could be given (Fig. 8).

The spatial structure of biodiversity implies the existence of scaling of contingency of evolution in space as well as in time as well as its limits. Endemicity and provincialism are always present, but at long time scales, species from different provinces invade most environments suitable for their existence. Here tectonics and also climate fluctuations help the dispersal of taxa with different environmental tolerances and set limits on how different local and regional biotas can be and how long they can be effectively isolated (Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022). The same kind of processes of merging of lower-level Bretskyan units after some time of isolation and forced individuality should be expected to proceed at all spatial and temporal scales—smaller spatial scales imply shorter time scales of merging and therefore shorter time scales of contingent evolution.

The last point of importance is evolutionary interaction and the question of mechanisms of selection of Bretskyan units. In the case of small-scale Bretskyan units (endosymbiotic cells, lichen, etc.), holobionts could be targets of Type II multilevel selection (MLS2; Damuth and Heisler Reference Damuth and Heisler1988) when there are differences in reproductive capacities of holobionts themselves (“groups” produce more “groups”). In the case of large-scale Bretskyan units, which should have significant influence on the fitness of taxa residing and being part of their causal network, we should expect Type I multilevel selection (MLS1) when the differences in fitness reside in components of the groups (Damuth and Heisler Reference Damuth and Heisler1988; Okasha Reference Okasha2006; Lean and Jones Reference Lean and Jones2023). The special case here is the global holobiont-like unit Gaia (zu Castell et al. Reference zu Castell, Lüttge, Matyssek, Wegner and Lüttge2019), which is singular, the largest and indivisible in the long-term Bretskyan individual. The Gaian unit is modulating fitness of all lower-level units (Bourrat Reference Bourrat2023) and also possesses its own fitness, defined here as resilience to its ultimate extinction—therefore it combines features of both Type I and Type II multilevel fitness. Contrary to popular conception and models that Gaia increases its chances of survival by means of regulation of environment (Lenton Reference Lenton1998; Arthur and Nicholson Reference Arthur and Nicholson2022), current evidence shows that such regulation is absent (Lovejoy Reference Lovejoy2015; Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022), and life resists this unchecked change at the longest time scales (Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022, Reference Spiridonov and Lovejoy2023).

The effects of belonging to a given Bretskyan unit can be measured by differential survival or proliferation of taxa that originated in a geobiome. In this way, geobiomes do not produce more of the same geobiomes, as is the case for organism-level selection, but “successful” geobiomes make future geobiomes resemble them compositionally (Fig. 9). The merging of Bretskyan units during collision results in multiple invasions (Stigall Reference Stigall2019) that are natural experiments that can be used in testing of hierarchical scaling laws of spatiotemporal evolution of geobiomes at large scales. One of the finest and best-researched examples of such mergers of large Bretskyan units is GABI, which happened in the late Cenozoic between North America and South America (Marshall et al. Reference Marshall, Webb, Sepkoski and Raup1982). Apparently South American mammalian faunas experienced higher extinction rates in comparison to North American mammalian faunas (Carrillo et al. Reference Carrillo, Faurby, Silvestro, Zizka, Jaramillo, Bacon and Antonelli2020), which is in line with expectations of larger and more environmentally diverse geobiomes having higher competitiveness.

Figure 9. Interactions between higher-level geobiomes, here representing hypothetical bioprovinces experiencing sequential fusion events; these could be large islands or terrains isolated by a sea that merge sequentially during eustatic regression. Compositions of bioprovinces are represented by colors. Type I multilevel selection (MLS1) Bretskyan fitness of provinces B and C is equal (competitive ability of component taxa is the same); therefore, the newly merged unit BC has exactly intermediate composition. Bretskyan unit A has much higher fitness than unit BC; therefore its taxa dominate unit ABC, which is formed by the final merger.

Apparently, as was shown by the case studies of modern South and North American freshwater fishes, the diversity history and the distribution of species richness is significantly affected by dynamics of sub-continental structures, even at scales of millions to tens of millions of years, thus showing the great influence of Bretskyan units in molding some of the most recent evolutionary radiations (Cassemiro et al. Reference Cassemiro, Albert, Antonelli, Menegotto, Wüest, Cerezer, Coelho, Reis, Tan and Tagliacollo2023; Stokes et al. Reference Stokes, Kim, Gallen, Benavides, Keck, Wood, Goldberg, Larsen, Mollish and Simmons2023). Therefore, we can see that fragmentation and dynamics of landscapes at multiple spatiotemporal scales are crucial in modulating small-scale taxonomic changes, as well as changes in whole faunas—the mosaic of evolution and its coordination is sustained and reshaped by branching (multiscale allopatry) and merging (or the Geo-Red Queen [Spiridonov and Lovejoy Reference Spiridonov and Lovejoy2022]) of landscapes by means of tectonic geodynamics and climate change.

Many more such studies across many clades, areas, times, and time scales (Patzkowsky Reference Patzkowsky2017; Stigall et al. Reference Stigall, Bauer, Lam and Wright2017) are needed, including cases of geographic differences in recoveries after mass extinctions and other significant events (Jablonski Reference Jablonski2005), in order to generalize the patterns of evolution and to derive scaling laws applicable to different levels of the Bretskyan hierarchy.

The major purpose of this paper is the explication of the basic biological ontology across all scales. Nonetheless, the described properties of the Bretskyan units and the processes that determine and form the Bretskyan hierarchy at large scales of geobiomes present a range of empirical predictions that, as in the case of the “first law of biology” (McShea and Brandon Reference McShea and Brandon2010), should be understood as idealized patterns, which can be best detected at ceteris paribus conditions. Predictions listed here are related to the greater coevolutionary depth of taxa and greater match to a wider range of environments with which taxa interacted in space and time while being part of larger, older, and more structured geobiomes (by possessing more of the Fisherian information and consequently being more “learned” and more integrated with each other) than in smaller, younger, and less complex geobiomes. These are, upon fusion of geobiomes: (1) taxa from larger/older geobiomes will on average occupy a wider range of environments; (2) taxa from larger/older geobiomes will have higher chances of invasion and establishment; (3) taxa from larger/older geobiomes will be found in ecological associations more frequently and for longer geological time; (4) taxa from geodynamically more active/complex regions (more spatiotemporally structured geobiomes) should exhibit higher fitness given the same area, climate, and the age of a geobiome; (5) taxa from larger/older geobiomes should have longer survivorship half-lives; and (6) average duration of a taxon should increase as time passes in a given geobiome, consequently if sampled at arbitrary time (e.g., recent), the average taxon ages should be longer in an older/larger geobiome (e.g., marine vs. terrestrial geobiomes).

Regarding prediction 2: The modern ecosystems that are highly influenced by anthropogenic factors are affected by massive influx of invasive species (Stigall Reference Stigall2019), the most disruptive of which is arguably humans. The Bretskyan hierarchy approach gives a predictive framework on how this increased dispersal (or fusion of multiscale geobiomes) event can proceed. Predictions span the whole range of time scales, from immediate establishment, to the environmental breadth of dispersal, to ultimately the differential evolutionary success of the invaders. The same principles could be cross-tested in deep-time settings, thus providing deeper insights into strength of geobiomic factors in ecology and evolution.