What predisposes two individuals to form and maintain a relationship is a fundamental question. Using facial recognition to see whether couples' faces change over time to become more and more similar, psychology researchers have concluded that couples tend to be formed from the start between people whose faces are more similar than average [1]. As the saying goes, birds of a feather flock together.

And what about in nature? Are these rules of assembly valid for communities of different species?

In his seminal contribution, Robert MacArthur (1984) wrote ‘To do science is to search for repeated patterns’ [2]. Identifying the mechanisms that govern the arrangement of life is a hot research topic in the field of ecology for decades, and an absolutely essential prerequisite to answer the outstanding question of what shape ecological patterns in multi-species communities such as species-area relationships, relative species abundances, or spatial and temporal turnover of community composition; amid others [3]. To explain ecological patterns in nature, some rely on the concept that every species - through evolutionary processes and the acquisition of a unique set of traits that allow a species to be adapted to its abiotic and biotic environment - occupies a unique niche: Species coexistence comes as the result of niche differentiation [4,5]. Such a view has been challenged by the recognition of the key role of neutral processes [6], however, in which demographic stochasticity contributes to shape multi-species communities and to explain why congener species coexist much more frequently than expected by chance [7,8]. While the niche-based and neutral theories appear seemingly opposed at first sight [9], the dichotomy may be more philosophical than empirical [4,5]. Many examples have come to support that both concepts are not incompatible as they together influence the structure, diversity and functioning of communities [10], and are simply extreme cases of a continuum [11]. From this perspective, extrinsic factors, i.e., environmental heterogeneity, may influence the location of a given community along the niche-neutrality continuum. 

The walk of species in nature is therefore neither random nor ecologically predestined. In microbial assemblages, the co-existence of these two antagonistic mechanisms has been shown both theoretically and empirically. It has been shown that a combination of stabilising (niche) and equalising (neutral) mechanisms was responsible for the existence of groups of coexistent species (clumps) in a phytoplankton rich community [12]. Analysing interannual changes (2003-2009) in the weekly abundance of diatoms and dinoflagellates located in a temperate coastal ecosystem of the Western English Channel, Mutshinda et al. [13] found a mixture of biomass dynamics consistent with the neutrality-niche continuum hypothesis. While niche processes explained the dynamic of phytoplankton functional groups (i.e., diatoms vs. dinoflagellates) in terms of biomass, neutral processes mainly dominated - 50 to 75% of the time - the dynamics at the species level within functional groups [13]. From one endpoint to another, defining the location of a community along the continuum is all matter of scale [4,11].

In their study, testing predictions made by an emergent neutrality model, Graco-Roza et al. [14] provide empirical evidence that neutral and niche processes joined together to shape and drive planktonic communities in a riverine ecosystem. Body size - the 'master trait' - is used here as a discriminant ecological dimension along the niche axis. From their analysis, they not only show that the specific abundance is organised in clumps and gaps along the niche axis, but also reveal that different clumps exist along the river course. They identify two main clumps in body size - with species belonging to three different morphologically-based functional groups - and characterise that among-species differences in biovolume are driven by functional redundancy at the clump level; species functional distinctiveness being related to the relative biovolume of species. By grouping their variables according to seasons (cold-dry vs. warm-wet) or river elevation profile (upper, medium and lower course), they hereby highlight how environmental heterogeneity contributes to shape species assemblages and their dynamics and conclude that emergent neutrality models are a powerful approach to explain species coexistence; and therefore ecological patterns.

References

[1] Tea-makorn PP, Kosinski M (2020) Spouses’ faces are similar but do not become more similar with time. Scientific Reports, 10, 17001. https://doi.org/10.1038/s41598-020-73971-8.

[2] MacArthur RH (1984) Geographical Ecology: Patterns in the Distribution of Species. Princeton University Press.

[3] Vellend M (2020) The Theory of Ecological Communities (MPB-57). Princeton University Press.

[4] Wennekes PL, Rosindell J, Etienne RS (2012) The Neutral—Niche Debate: A Philosophical Perspective. Acta Biotheoretica, 60, 257–271. https://doi.org/10.1007/s10441-012-9144-6.

[5] Gravel D, Guichard F, Hochberg ME (2011) Species coexistence in a variable world. Ecology Letters, 14, 828–839. https://doi.org/10.1111/j.1461-0248.2011.01643.x.

[6] Hubbell SP (2001) The Unified Neutral Theory of Biodiversity and Biogeography (MPB-32). Princeton University Press.

[7] Leibold MA, McPeek MA (2006) Coexistence of the Niche and Neutral Perspectives in Community Ecology. Ecology, 87, 1399–1410. https://doi.org/10.1890/0012-9658(2006)87[1399:COTNAN]2.0.CO;2.

[8] Pielou EC (1977) The Latitudinal Spans of Seaweed Species and Their Patterns of Overlap. Journal of Biogeography, 4, 299–311. https://doi.org/10.2307/3038189.

[9] Holt RD (2006) Emergent neutrality. Trends in Ecology & Evolution, 21, 531–533. https://doi.org/10.1016/j.tree.2006.08.003. 

[10] Scheffer M, Nes EH van (2006) Self-organized similarity, the evolutionary emergence of groups of similar species. Proceedings of the National Academy of Sciences, 103, 6230–6235. https://doi.org/10.1073/pnas.0508024103.

[11] Gravel D, Canham CD, Beaudet M, Messier C (2006) Reconciling niche and neutrality: the continuum hypothesis. Ecology Letters, 9, 399–409. https://doi.org/10.1111/j.1461-0248.2006.00884.x.

[12] Vergnon R, Dulvy NK, Freckleton RP (2009) Niches versus neutrality: uncovering the drivers of diversity in a species-rich community. Ecology Letters, 12, 1079–1090. https://doi.org/10.1111/j.1461-0248.2009.01364.x.

[13] Mutshinda CM, Finkel ZV, Widdicombe CE, Irwin AJ (2016) Ecological equivalence of species within phytoplankton functional groups. Functional Ecology, 30, 1714–1722. https://doi.org/10.1111/1365-2435.12641.

[14] Graco-Roza C, Segura AM, Kruk C, Domingos P, Soininen J, Marinho MM (2021) Clumpy coexistence in phytoplankton: The role of functional similarity in community assembly. bioRxiv, 869966, ver. 6 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/869966

In this article, Saade et al. (2021) investigate how the rate of local extinctions and their spatial distribution affect recolonization dynamics in metacommunities. They use an elegant combination of microcosm experiments with metacommunities of freshwater ciliates and mathematical modelling mirroring their experimental system. Their main findings are (i) that local patch extinctions increase both local (α-) and inter-patch (β-) diversity in a transient way during the recolonization process, (ii) that these effects depend more on the spatial distribution of extinctions (dispersed or clustered) than on their amount, and (iii) that they may spread regionally.
Microcosm experiments are already quite cool just by themselves and have contributed largely to conceptual advances in community ecology (see Fraser and Keddy 1997, or Jessup et al. 2004 for reviews on this topic), but they are here exploited to a whole further level by the fitting of a metapopulation dynamics model. The model allows both to identify the underlying mechanisms most likely to generate the patterns observed (here, competitive interactions) and to assess the robustness of these patterns when considering larger spatial or temporal scales. This release of experimental limitations allows here for the analysis of quantitative metrics of spatial structure, like the distance to the closest patch, which gives an interesting insight into the functional basis of the effect of the spatial distribution of extinctions.

A major strength of this study is that it highlights the importance of considering the spatial structure explicitly. Recent work on ecological networks has shown repeatedly that network structure affects the propagation of pathogens (Badham and Stocker 2010), invaders (Morel-Journel et al. 2019), or perturbation events (Gilarranz et al. 2017). Here, the spatial structure of the metacommunity is a regular grid of patches, but the distribution of extinction events may be either regularly dispersed (i.e., extinct patches are distributed evenly over the grid and are all surrounded by non-extinct patches only) or clustered (all extinct patches are neighbours). This has a direct effect on the neighbourhood of perturbed patches, and because perturbations have mostly local effects, their recovery dynamics are dominated by the composition of this immediate neighbourhood. In landscapes with dispersed extinctions, the neighbourhood of a perturbed patch is not affected by the amount of extinctions, and neither is its recovery time. In contrast, in landscapes with clustered extinctions, the amount of extinctions affects the depth of the perturbed area, which takes longer to recover when it is larger. Interestingly, the spatial distribution of extinctions here is functionally equivalent to differences in connectivity between perturbed and unperturbed patches, which results in contrasted “rescue recovery” and “mixing recovery” regimes as described by Zelnick et al. (2019).
 
Furthermore, this study focuses on local dynamics of competition and short-term, transient patterns that may have been overlooked by more classical, equilibrium-based approaches of dynamical systems of metacommunities. Indeed, in a metacommunity composed of several competitors, early theoretical work demonstrated that species coexistence is possible at the regional scale only, provided that spatial heterogeneity creates spatial variance in fitness or precludes the superior competitor from accessing certain habitat patches (Skellam 1951, Levins 1969). In the spatially homogeneous experimental system of Saade et al., one of the three ciliate species ends up dominating the community at equilibrium. However, following local, one-time extinction events, the community endures a recolonization process in which differences in dispersal may provide temporary spatial niches for inferior competitors. These transient patterns might prove essential to understand and anticipate the resilience of natural systems that are under increasing pressure, and enduring ever more frequent and intense perturbations (IPBES 2019). Spatial autocorrelation in extinction events was previously identified as a risk for stability and persistence of metacommunities (Ruokolainen 2013, Kahilainen et al. 2018). These new results show that autocorrelated perturbations also have longer-lasting effects, which is likely to increase their overall impact on metacommunity dynamics. As spatial and temporal autocorrelation of temperature and extreme climatic events are expected to increase (Di Cecco and Gouthier 2018), studies that investigate how metacommunities respond to the structure of the distribution of perturbations are more necessary than ever.
 
References

Badham J, Stocker R (2010) The impact of network clustering and assortativity on epidemic behaviour. Theoretical Population Biology, 77, 71–75. https://doi.org/10.1016/j.tpb.2009.11.003
 
Di Cecco GJ, Gouhier TC (2018) Increased spatial and temporal autocorrelation of temperature under climate change. Scientific Reports, 8, 14850. https://doi.org/10.1038/s41598-018-33217-0
 
Fraser LH, Keddy P (1997) The role of experimental microcosms in ecological research. Trends in Ecology & Evolution, 12, 478–481. https://doi.org/10.1016/S0169-5347(97)01220-2
 
Gilarranz LJ, Rayfield B, Liñán-Cembrano G, Bascompte J, Gonzalez A (2017) Effects of network modularity on the spread of perturbation impact in experimental metapopulations. Science, 357, 199–201. https://doi.org/10.1126/science.aal4122
 
IPBES (2019) Summary for policymakers of the global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services. S. Díaz, J. Settele, E. S. Brondízio E.S., H. T. Ngo, M. Guèze, J. Agard, A. Arneth, P. Balvanera, K. A. Brauman, S. H. M. Butchart, K. M. A. Chan, L. A. Garibaldi, K. Ichii, J. Liu, S. M. Subramanian, G. F. Midgley, P. Miloslavich, Z. Molnár, D. Obura, A. Pfaff, S. Polasky, A. Purvis, J. Razzaque, B. Reyers, R. Roy Chowdhury, Y. J. Shin, I. J. Visseren-Hamakers, K. J. Willis, and C. N. Zayas (eds.). IPBES secretariat, Bonn, Germany. 56 pages. https://doi.org/10.5281/zenodo.3553579 
 
Jessup CM, Kassen R, Forde SE, Kerr B, Buckling A, Rainey PB, Bohannan BJM (2004) Big questions, small worlds: microbial model systems in ecology. Trends in Ecology & Evolution, 19, 189–197. https://doi.org/10.1016/j.tree.2004.01.008
 
Kahilainen A, van Nouhuys S, Schulz T, Saastamoinen M (2018) Metapopulation dynamics in a changing climate: Increasing spatial synchrony in weather conditions drives metapopulation synchrony of a butterfly inhabiting a fragmented landscape. Global Change Biology, 24, 4316–4329. https://doi.org/10.1111/gcb.14280

Levins R (1969) Some Demographic and Genetic Consequences of Environmental Heterogeneity for Biological Control1. Bulletin of the Entomological Society of America, 15, 237–240. https://doi.org/10.1093/besa/15.3.237
 
Morel-Journel T, Assa CR, Mailleret L, Vercken E (2019) Its all about connections: hubs and invasion in habitat networks. Ecology Letters, 22, 313–321. https://doi.org/10.1111/ele.13192

Ruokolainen L (2013) Spatio-Temporal Environmental Correlation and Population Variability in Simple Metacommunities. PLOS ONE, 8, e72325. https://doi.org/10.1371/journal.pone.0072325

Saade C, Kefi S, Gougat-Barbera C, Rosenbaum B, Fronhofer EA (2021) Spatial distribution of local patch extinctions drives recovery dynamics in metacommunities. bioRxiv, 2020.12.03.409524, ver. 4 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2020.12.03.409524
 
Skellam JG (1951) Random Dispersal in Theoretical Populations. Biometrika, 38, 196–218. https://doi.org/10.2307/2332328
 
Zelnik YR, Arnoldi J-F, Loreau M (2019) The three regimes of spatial recovery. Ecology, 100, e02586. https://doi.org/10.1002/ecy.2586

One of the main issues of forest ecosystems is rising tree mortality as a result of extreme weather events (Franklin et al., 1987). Eventually, tree mortality reduces forest biomass (Allen et al., 2010), although its effect on forest ecosystem fluxes seems not lasting too long (Anderegg et al., 2016). This controversy about the negative consequences of tree mortality is joined to the debate about the drivers triggering and the mechanisms accelerating tree decline. For instance, there is still room for discussion about carbon starvation or hydraulic failure determining the decay processes (Sevanto et al., 2014) or about the importance of mortality sources (Reichstein et al., 2013). Therefore, understanding and predicting tree mortality has become one of the challenges for forest ecologists in the last decade, doubling the rate of articles published on the topic (*). Although predicting the responses of ecosystems to environmental change based on the traits of species may seem a simplistic conception of ecosystem functioning (Sutherland et al., 2013), identifying those traits that are involved in the proneness of a tree to die would help to predict how forests will respond to climate threatens.
Modelling tree mortality is complex, involving multiple factors acting simultaneously at different scales, from tree genetics to ecosystem dynamics and from microsite conditions to global climatic events. Therefore, taking into account different approaches to reduce uncertainty of the predictions is needed (Bugmann et al., 2019). Petit-Cailleux et al. (2020) uses statistical and process-based models to detect the main mortality drivers of a drought- and frost-prone beech population. Particularly, they assessed the intra-individual characteristics of the population, that may play a decisive role explaining the differences in tree vulnerability to extreme weather events. Comparing the results of both analytical approaches, they find out several key factors, such as defoliation, leaf phenology and tree size, that were consistent between them. Even more, the process-based model showed the physiological mechanisms that may explain the individual vulnerability, for instance higher loss of hydraulic conductance may increase the mortality risk of trees with early budburst phenology and large stem diameter. The authors also successfully model annual mortality rate with a linear relationship including only three parameters: loss of conductance, biomass of reserves and late frost days.
This valuable study is a good example of the complexity in understanding and predicting tree mortality. The authors carried out the ambitious commitment of studying the inter-annual variation in mortality with 14-year dataset. However, it might be not enough time to control for the dependence of temporal data to soundly model mortality rate. The authors also acknowledge that the use of two approaches increases the knowledge from different perspectives, but at the same time comparing their results is difficult because the parameters used are not identical. Particularly, process-based models tend to consider the same microclimatic conditions for every tree in the population, and may produce inconsistences with statistical models. Alternatively, individual-based modelling might overcome some of the incompatibilities between the approaches (Zhu et al., 2019).

(*) Number (and percentage) of articles found in Web of Sciences after searching (December the 10th, 2020) “tree mortality”: from 163 (0.006%) in 2010 to 412 (0.013%) in 2020.

References

Allen et al. (2010). A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Forest ecology and management, 259(4), 660-684. doi: https://doi.org/10.1016/j.foreco.2009.09.001
Anderegg et al. (2016). When a tree dies in the forest: scaling climate-driven tree mortality to ecosystem water and carbon fluxes. Ecosystems, 19(6), 1133-1147. doi: https://doi.org/10.1007/s10021-016-9982-1
Bugmann et al. (2019). Tree mortality submodels drive simulated long‐term forest dynamics: assessing 15 models from the stand to global scale. Ecosphere, 10(2), e02616. doi: https://doi.org/10.1002/ecs2.2616
Franklin, J. F., Shugart, H. H. and Harmon, M. E. (1987) Death as an ecological process: the causes, consequences, and variability of tree mortality. BioScience, 37, 550–556. doi: https://doi.org/10.2307/1310665
Petit-Cailleux, C., Davi, H., Lefèvre, F., Garrigue, J., Magdalou, J.-A., Hurson, C., Magnanou, E. and Oddou-Muratorio, S. (2020) Comparing statistical and mechanistic models to identify the drivers of mortality within a rear-edge beech population. bioRxiv, 645747, ver 7 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/645747
Reichstein et al. (2013). Climate extremes and the carbon cycle. Nature, 500(7462), 287-295. doi: https://doi.org/10.1038/nature12350
Sevanto, S., Mcdowell, N. G., Dickman, L. T., Pangle, R., and Pockman, W. T. (2014). How do trees die? A test of the hydraulic failure and carbon starvation hypotheses. Plant, cell & environment, 37(1), 153-161. doi: https://doi.org/10.1111/pce.12141
Sutherland et al. (2013). Identification of 100 fundamental ecological questions. Journal of ecology, 101(1), 58-67. doi: https://doi.org/10.1111/1365-2745.12025
Zhu, Y., Liu, Z., and Jin, G. (2019). Evaluating individual-based tree mortality modeling with temporal observation data collected from a large forest plot. Forest Ecology and Management, 450, 117496. doi: https://doi.org/10.1016/j.foreco.2019.117496

In this study, the authors report a new method for estimating the abundance of the Pyrenean brown bear population. Precisely, the methodology involved aims to apply Pollock's closed robust design (PCRD) capture-recapture models to estimate population abundance and trends over time. Overall, the results encourage the use of PCRD to study populations' demographic rates, while minimizing biases due to inter-individual heterogeneity in detection probabilities.

Estimating the size and trends of animal population over time is essential for informing conservation status and management decision-making (Nichols & Williams 2006). This is particularly the case when the population is small, geographically scattered, and threatened. Although several methods can be used to estimate population abundance, they may be difficult to implement when individuals are rare, elusive, solitary, largely nocturnal, highly mobile, and/or occupy large home ranges in remote and/or rugged habitats. Moreover, in such standard methods,

• the population is assumed to be closed both geographically (no immigration nor emigration) and demographically (no births nor deaths) and
• all individuals are assumed to have identical detection probabilities regardless of their individual attributes (e.g., age, body mass, social status) and habitat features (home-range location and composition) (Otis et al. 1978).

However, these conditions are rarely met in real populations, such as wild mammals (e.g., Bellemain et al. 2005; Solbert et al. 2006), and therefore the risk of underestimating population size can rapidly increase because the assumption of perfect detection of all individuals in the population is violated.

Focusing on the critically endangered Pyrenean brown bear that was close to extinction in the mid-1990s, the study by Vanpe et al. (2022), uses protocol and opportunistic data to describe a statistical modeling exercise to construct mark-recapture histories from 2008 to 2020. Among the data, the authors collected non-invasive samples such as a mixture of hair and scat samples used for genetic identification, as well as photographic trap data of recognized individuals. These data are then analyzed in RMark to provide detection and survival estimates. The final model (i.e. PCRD capture-recapture) is then used to provide Bayesian population estimates. Results show a five-fold increase in population size between 2008 and 2020, from 13 to 66 individuals. Thus, this study represents the first published annual abundance and temporal trend estimates of the Pyrenean brown bear population since 2008.

Then, although the results emphasize that the PCRD estimates were broadly close to the MRS counts and had reasonably narrow associated 95% Credibility Intervals, they also highlight that the sampling effort is different according to individuals. Indeed, as expected, the detection of an individual depends on

• the intraspecific home range size variation that results in individuals that move the most being most likely to be detected and
• the mortality rate which is higher on cubs than on adults and subadults (due to infanticide by males, predation, death of the mother, or abandonment).

Overall, the PCRD capture-recapture modelling approach, involved in this study, provides robust estimates of abundance and demographic rates of the Pyrenean brown bear population (with associated uncertainty) while minimizing and considering bias due to inter-individual heterogeneity in detection probabilities.

The authors conclude that mark-recapture provides useful population estimates and urge wildlife ecologists and managers to use robust approaches, such as the RDPC capture-recapture model, when studying large mammal populations. This information is essential to inform management decisions and assess the conservation status of populations.

References

Bellemain, E.V.A., Swenson, J.E., Tallmon, D., Brunberg, S. and Taberlet, P. (2005). Estimating population size of elusive animals with DNA from hunter-collected feces: four methods for brown bears. Cons. Biol. 19(1), 150-161. https://doi.org/10.1111/j.1523-1739.2005.00549.x

Nichols, J.D. and Williams, B.K. (2006). Monitoring for conservation. Trends Ecol. Evol. 21(12), 668-673. https://doi.org/10.1016/j.tree.2006.08.007

Otis, D.L., Burnham, K.P., White, G.C. and Anderson, D.R. (1978). Statistical inference from capture data on closed animal populations. Wildlife Monographs (62), 3-135.

Solberg, K.H., Bellemain, E., Drageset, O.M., Taberlet, P. and Swenson, J.E. (2006). An evaluation of field and non-invasive genetic methods to estimate brown bear (Ursus arctos) population size. Biol. Conserv. 128(2), 158-168. https://doi.org/10.1016/j.biocon.2005.09.025

Vanpé C, Piédallu B, Quenette P-Y, Sentilles J, Queney G, Palazón S, Jordana IA, Jato R, Elósegui Irurtia MM, de la Torre JS, and Gimenez O (2022) Estimating abundance of a recovering transboundary brown bear population with capture-recapture models. bioRxiv, 2021.12.08.471719, ver. 4 recommended and peer-reviewed by PCI Ecology. https://doi.org/10.1101/2021.12.08.471719

Modelling the temporal dynamics of trophic networks has been a key challenge for community ecologists for decades, especially when anthropogenic and natural forces drive changes in species composition, abundance, and interactions over time. So far, most modelling methods fail to incorporate the inherent complexity of such systems, and its variability, to adequately describe and predict temporal changes in the topology of trophic networks. 

Taking benefit from theoretical computer science advances, Gaucherel and colleagues (2024) propose a new methodological framework to tackle this challenge based on discrete-event Petri net methodology. To introduce the concept to naïve readers the authors provide a useful example using a simplistic predator-prey model.

The core biological system of the article is a freshwater trophic network of western France in the Charente-Maritime marshes of the French Atlantic coast. A directed graph describing this system was constructed to incorporate different functional groups (phytoplankton, zooplankton, resources, microbes, and abiotic components of the environment) and their interactions. Rules and constraints were then defined to, respectively, represent physiochemical, biological, or ecological processes linking network components, and prevent the model from simulating unrealistic trajectories. Then the full range of possible trajectories of this mechanistic and qualitative model was computed.

The model performed well enough to successfully predict a theoretical trajectory plus two trajectories of the trophic network observed in the field at two different stations, therefore validating the new methodology introduced here. The authors conclude their paper by presenting the power and drawbacks of such a new approach to qualitatively model trophic networks dynamics.

Reference

Cedric Gaucherel, Stolian Fayolle, Raphael Savelli, Olivier Philippine, Franck Pommereau, Christine Dupuy (2024) Diagnosis of planktonic trophic network dynamics with sharp qualitative changes. bioRxiv 2023.06.29.547055, ver. 2 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2023.06.29.547055