The geography of metapopulation synchrony in dendritic river networks

Heterogeneous dispersal networks to improve biodiversity science

Dispersal has a key role in shaping spatial patterns of biodiversity, yet its spatial heterogeneity is often overlooked in biodiversity analyses and management strategies. Properly parameterised heterogeneous dispersal networks capture the complex interplay between landscape structure and species-specific dispersal capacities. However, this heterogeneity is recurrently neglected when studying the processes underlying biodiversity variation. To address this gap, we introduce a conceptual framework detailing the fundamental processes driving dispersal heterogeneity and its effects on biodiversity dynamics. We propose methods to parameterise heterogeneous dispersal networks, facilitating their integration into commonly used quantitative frameworks for biodiversity analyses. By considering the architecture of heterogeneous dispersal networks, we demonstrate their critical role in guiding biodiversity management strategies.

Climatic effects on the synchrony and stability of temperate headwater invertebrates over four decades

Important clues about the ecological effects of climate change can arise from understanding the influence of other Earth-system processes on ecosystem dynamics but few studies span the inter-decadal timescales required. We, therefore, examined how variation in annual weather patterns associated with the North Atlantic Oscillation (NAO) over four decades was linked to synchrony and stability in a metacommunity of stream invertebrates across multiple, contrasting headwaters in central Wales (UK). Prolonged warmer and wetter conditions during positive NAO winters appeared to synchronize variations in population and community composition among and within streams thereby reducing stability across levels of ecological organization. This climatically mediated synchronization occurred in all streams irrespective of acid-base status and land use, but was weaker where invertebrate communities were more functionally diverse. Wavelet linear models indicated that variation in the NAO explained up to 50% of overall synchrony in species abundances at a timescale of 4-6 years. The NAO appeared to affect ecological dynamics through local variations in temperature, precipitation and discharge, but increasing hydrochemical variability within sites during wetter winters might have contributed. Our findings illustrate how large-scale climatic fluctuations generated over the North Atlantic can affect population persistence and dynamics in inland freshwater ecosystems in ways that transcend local catchment character. Protecting and restoring functional diversity in stream communities might increase their stability against warmer, wetter conditions that are analogues of ongoing climate change. Catchment management could also dampen impacts and provide options for climate change adaptation.

The effects of longitudinal fragmentation on riverine beta diversity are modulated by fragmentation intensity

The loss of longitudinal connectivity affects river systems globally, being one of the leading causes of the freshwater biodiversity crisis. Barriers alter the dispersal of aquatic organisms and limit the exchange of species between local communities, disrupting metacommunity dynamics. However, the interplay between connectivity losses due to dams and other drivers of metacommunity structure, such as the configuration of the river network, needs to be explored. In this paper, we analyzed the response of fish communities to the network position and the fragmentation induced by dams while controlling for human pressures and environmental gradients. We studied three large European catchments covering a fragmentation gradient: Upper Danube (Austrian section), Ebro (Spain), and Odra/Oder (Poland). We quantified fragmentation through reach-scaled connectivity indices that account for the position of barriers along the dendritic network and the dispersal capacity of the organisms. We used generalized linear models to explain species richness and Local Contributions to Beta Diversity (LCBD) and multilinear regressions on the distance matrix to describe Beta Diversity and its Replacement and Richness Difference components. Results show that species richness was not affected by fragmentation. Network centrality metrics were relevant drivers of beta diversity for catchments with lower fragmentation (Ebro, Odra), and fragmentation indices were strong beta diversity predictors for the catchment with higher fragmentation (Danube). We conclude that in highly fragmented catchments, the effects of network centrality/isolation on biodiversity could be masked by the effects of dam fragmentation. In such catchments, metapopulation and metacommunity dynamics can be strongly altered by barriers, and the restoration of longitudinal connectivity (i.e. the natural centrality/isolation gradient) is urgent to prevent local extinctions.

Non-perennial segments in river networks

Non-perennial river segments - those that recurrently cease to flow or frequently dry - occur in all river networks and are globally more abundant than perennial (always flowing) segments. However, research and management have historically focused on perennial river segments. In this Review, we outline how non-perennial segments are integral parts of river networks. Repeated cycles of flowing, non-flowing and dry phases in non-perennial segments influence biodiversity and ecosystem dynamics at different spatial scales, from individual segments to entire river networks. Varying configurations of perennial and non-perennial segments govern physical, chemical and ecological responses to changes in the flow regimes of each river network, especially in response to human activities. The extent of non-perennial segments in river networks has increased owing to warming, changing hydrological patterns and human activities, and this increase is predicted to continue. Moreover, the dry phases of flow regimes are expected to be longer, drier and more frequent, albeit with high regional variability. These changes will likely impact biodiversity, potentially tipping some ecosystems to compromised stable states. Effective river-network management must recognize ecosystem services (such as flood risk management and groundwater recharge) provided by non-perennial segments and ensure their legislative and regulatory protection, which is often lacking.

Modelling environmental DNA transport in rivers reveals highly resolved spatio-temporal biodiversity patterns

The ever-increasing threats to riverine ecosystems call for novel approaches for highly resolved biodiversity assessments across taxonomic groups and spatio-temporal scales. Recent advances in the joint use of environmental DNA (eDNA) data and eDNA transport models in rivers (e.g., eDITH) allow uncovering the full structure of riverine biodiversity, hence elucidating ecosystem processes and supporting conservation measures. We applied eDITH to a metabarcoding dataset covering three taxonomic groups (fish, invertebrates, bacteria) and three seasons for a catchment sampled for eDNA at 73 sites. We upscaled eDNA-based biodiversity predictions to approximately 1900 reaches, and assessed α- and β-diversity patterns across seasons and taxonomic groups over the whole network. Genus richness predicted by eDITH was generally higher than values from direct eDNA analysis. Both predicted α- and β-diversity varied depending on season and taxonomic group. Predicted fish α-diversity increased downstream in all seasons, while invertebrate and bacteria α-diversity either decreased downstream or were unrelated to network position. Spatial β-diversity mostly decreased downstream, especially for bacteria. The eDITH model yielded a more refined assessment of freshwater biodiversity as compared to raw eDNA data, both in terms of spatial coverage, diversity patterns and effect of covariates, thus providing a more complete picture of freshwater biodiversity.

Destabilizing Effects of Environmental Stressors on Aquatic Communities and Interaction Networks across a Major River Basin

Human-driven environmental stressors are increasingly threatening species survival and diversity of river systems worldwide. However, it remains unclear how the stressors affect the stability changes across aquatic multiple communities. Here, we used environmental DNA (eDNA) data sets from a human-dominated river in China over 3 years and analyzed the stability changes in multiple communities under persistent anthropogenic stressors, including land use and pollutants. First, we found that persistent stressors significantly reduced multifaceted species diversity (e.g., species richness, Shannon's diversity, and Simpson's diversity) and species stability but increased species synchrony across multiple communities. Second, the structures of interaction networks inferred from an empirical meta-food web were significantly changed under persistent stressors, for example, resulting in decreased network modularity and negative/positive cohesion. Third, piecewise structural equation modeling proved that the persistent stress-induced decline in the stability of multiple communities mainly depended upon diversity-mediated pathways rather than the direct effects of stress per se; specifically, the increase of species synchrony and the decline of interaction network modularity were the main biotic drivers of stability variation. Overall, our study highlights the destabilizing effects of persistent stressors on multiple communities as well as the mechanistic dependencies, mainly through reducing species diversity, increasing species synchrony, and changing interaction networks.

The color of environmental noise in river networks

Despite its far-reaching implications for conservation and natural resource management, little is known about the color of environmental noise, or the structure of temporal autocorrelation in random environmental variation, in streams and rivers. Here, we analyze the geography, drivers, and timescale-dependence of noise color in streamflow across the U.S. hydrography, using streamflow time series from 7504 gages. We find that daily and annual flows are dominated by red and white spectra respectively, and spatial variation in noise color is explained by a combination of geographic, hydroclimatic, and anthropogenic variables. Noise color at the daily scale is influenced by stream network position, and land use and water management explain around one third of the spatial variation in noise color irrespective of the timescale considered. Our results highlight the peculiarities of environmental variation regimes in riverine systems, and reveal a strong human fingerprint on the stochastic patterns of streamflow variation in river networks.

Geographic patterns of taxonomic and phylogenetic β-diversity of aquatic angiosperms in China

China covers a vast territory harbouring a large number of aquatic plants. Although there are many studies on the β-diversity of total, herbaceous or woody plants in China and elsewhere, few studies have focused on aquatic plants. Here, we analyse a comprehensive data set of 889 aquatic angiosperm species in China, and explore the geographic patterns and climatic correlates of total taxonomic and phylogenetic β-diversity as well as their turnover and nestedness components. Our results show that geographic patterns of taxonomic and phylogenetic β-diversity are highly congruent for aquatic angiosperms, and taxonomic β-diversity is consistently higher than phylogenetic β-diversity. The ratio between the nestedness component and total β-diversity is high in northwestern China and low in southeastern China. The geographic patterns of taxonomic and phylogenetic β-diversity of aquatic angiosperms in China are obviously affected by geographic and climatic distances, respectively. In conclusion, the geographic patterns of taxonomic and phylogenetic β-diversity of aquatic angiosperms are consistent across China. Climatic and geographic distances jointly affect the geographic patterns of β-diversity of aquatic angiosperms. Overall, our work provides insight into understanding the large-scale patterns of aquatic angiosperm β-diversity, and is a critical addition to previous studies on the macroecological patterns of terrestrial organisms.

Scale of population synchrony confirms macroecological estimates of minimum viable range size

Global ecosystems are facing a deepening biodiversity crisis, necessitating robust approaches to quantifying species extinction risk. The lower limit of the macroecological relationship between species range and body size has long been hypothesized as an estimate of the relationship between the minimum viable range size (MVRS) needed for species persistence and the organismal traits that affect space and resource requirements. Here, we perform the first explicit test of this assumption by confronting the MVRS predicted by the range-body size relationship with an independent estimate based on the scale of synchrony in abundance among spatially separated populations of riverine fish. We provide clear evidence of a positive relationship between the scale of synchrony and species body size, and strong support for the MVRS set by the lower limit of the range-body size macroecological relationship. This MVRS may help prioritize first evaluations for unassessed or data-deficient taxa in global conservation assessments.

The effect of population model choice and network topologies on extinction time patterns in dendritic ecological networks

Models of populations in habitat networks are vital for understanding and linking processes and patterns across individuals, environments, ecological interactions, and population structures. River ecosystem models combine the physical structure of the networks with the biological processes of the organisms using structural and functional models, respectively. Previous studies on dendritic river networks have employed different functional (population) models and either directly claimed or implied that the results illustrate general properties of actual river systems. However, these studies have used different approaches and assumptions when modeling population characteristics and behavior, and it is possible that inferences regarding a system may vary based on the combination of functional model and the spatial structure of a network. This study aims to understand if different functional models in river systems produce substantially different model results and, therefore, whether conclusions are model-dependent. We compare variation in extinction time and occupancy proportion of river networks with linear, trellis, dendritic and ring-lattice topologies, using three population models (uniform, age-class and individual based) and one metapopulation-based (patch-occupancy) model. Dendritic, linear, and trellis structures did not show notable differences among extinction times for any of the four models. The difference between topologies was higher for the patch-occupancy model compared to the three population models. There were significant differences in the variations of patch-occupancy between the metapopulation and the population models, but the three population models of differing complexity produced broadly similar results. Therefore, if the occupancy data is obtained based on local subpopulations, spatial arrangement and connectivity does not appear to be the sole predictor of single-species metapopulation responses. We conclude that the outputs from functional models are robust to assumptions and varying levels of detail as long as they contain at least some detail at the level of individuals within habitat nodes. Also, if we are modeling network-scale populations, models that include at least some detailed information on individuals are a far better choice than considering populations implicitly.

Tail-dependent spatial synchrony arises from nonlinear driver-response relationships

Spatial synchrony may be tail-dependent, that is, stronger when populations are abundant than scarce, or vice-versa. Here, 'tail-dependent' follows from distributions having a lower tail consisting of relatively low values and an upper tail of relatively high values. We present a general theory of how the distribution and correlation structure of an environmental driver translates into tail-dependent spatial synchrony through a non-linear response, and examine empirical evidence for theoretical predictions in giant kelp along the California coastline. In sheltered areas, kelp declines synchronously (lower-tail dependence) when waves are relatively intense, because waves below a certain height do little damage to kelp. Conversely, in exposed areas, kelp is synchronised primarily by periods of calmness that cause shared recovery (upper-tail dependence). We find evidence for geographies of tail dependence in synchrony, which helps structure regional population resilience: areas where population declines are asynchronous may be more resilient to disturbance because remnant populations facilitate reestablishment.

Asynchrony, density dependence, and persistence in an amphibian

Understanding drivers of metapopulation dynamics remains a critical challenge for ecology and conservation. In particular, the degree of synchrony in metapopulation dynamics determines how resilient a metapopulation is to a widespread disturbance. In this study, we used 21 years of egg mass count data across 64 nonpermanent freshwater ponds in Connecticut, USA to evaluate patterns of abundance and growth and to assess regional as well as local factors in shaping the population dynamics of wood frogs (Rana sylvatica = Lithobates sylvaticus). In particular, we asked whether a species known to undergo metapopulation dynamics exhibited spatial synchrony in abundances. With the exception of a single year when breeding took place during severe drought conditions, our analyses revealed no evidence of synchrony despite close proximity (mean minimum distance <300 m) of breeding ponds across the 3213 ha study area. Instead, local, pond-scale conditions best predicted patterns of abundance and population growth rate. We found negative density dependence on population growth rate within ponds as well as evidence that larger neighboring pond populations had a negative effect on focal ponds. Beyond density, pond depth was a critical predictor; deeper ponds supported larger populations. Drought conditions and warm winters negatively affected populations. Overall, breeding ponds vary in critical ways that either support larger, more persistent populations or smaller populations that are not represented by breeding pairs in some years. The infrequency of spatial synchrony in this system is surprising and suggests greater resilience to stressors than would have been expected if dynamics were strongly synchronized. More generally, understanding the characteristics of systems that determine synchronous population dynamics will be critical to predicting which species are more or less resilient to widespread disturbances like land conversion or climate change.

Spatial patterns in the contribution of biotic and abiotic factors to the population dynamics of three freshwater fish species

BACKGROUND

Population dynamics are driven by a number of biotic (e.g., density-dependence) and abiotic (e.g., climate) factors whose contribution can greatly vary across study systems (i.e., populations). Yet, the extent to which the contribution of these factors varies across populations and between species and whether spatial patterns can be identified has received little attention.

METHODS

Here, we used a long-term (1982-2011), broad scale (182 sites distributed across metropolitan France) dataset to study spatial patterns in the population's dynamics of three freshwater fish species presenting contrasted life-histories and patterns of elevation range shifts in recent decades. We used a hierarchical Bayesian approach together with an elasticity analysis to estimate the relative contribution of a set of biotic (e.g., strength of density dependence, recruitment rate) and abiotic (mean and variability of water temperature) factors affecting the site-specific dynamic of two different size classes (0⁺ and >0⁺ individuals) for the three species. We then tested whether the local contribution of each factor presented evidence for biogeographical patterns by confronting two non-mutually exclusive hypotheses: the "range-shift" hypothesis that predicts a gradient along elevation or latitude and the "abundant-center" hypothesis that predicts a gradient from the center to the edge of the species' distributional range.

RESULTS

Despite contrasted life-histories, the three species displayed similar large-scale patterns in population dynamics with a much stronger contribution of biotic factors over abiotic ones. Yet, the contribution of the different factors strongly varied within distributional ranges and followed distinct spatial patterns. Indeed, while abiotic factors mostly varied along elevation, biotic factors-which disproportionately contributed to population dynamics-varied along both elevation and latitude.

CONCLUSIONS

Overall while our results provide stronger support for the range-shift hypothesis, they also highlight the dual effect of distinct factors on spatial patterns in population dynamics and can explain the overall difficulty to find general evidence for geographic gradients in natural populations. We propose that considering the separate contribution of the factors affecting population dynamics could help better understand the drivers of abundance-distribution patterns.

Biological traits can mediate species-specific, quasi-extinction risks of macroinvertebrates in streams experiencing frequent extreme floods

Most research on the ecological responses to extreme floods examines impacts at short time scales, whereas long-term datasets combining hydrological and biological information remain rare. Using such data, we applied time-series analysis to investigate simultaneous effects of a biotic factor (density dependence), an abiotic factor (extreme floods), and spatial synchrony in the population dynamics of three riverine insects. Spatial synchronization of population dynamics by extreme floods varied among species. These different responses to extreme floods could be explained by species-specific biological traits. Moreover, density dependence influenced the population dynamics under the context of extreme floods. Accordingly, quasi-extinction risks were highest for species that were simultaneously influenced by biotic and abiotic factors. An understanding of ecological responses to increasing hydrological extremes may be enhanced by recognizing long-term, climatic non-stationarity.

Size, connectivity and edge effects of stream habitats explain spatio-temporal variation in brown trout (Salmo trutta) density

Ecological theory postulates that the size and isolation of habitat patches impact the colonization/extinction dynamics that determine community species richness and population persistence. Given the key role of lotic habitats for life-history completion in rheophilic fish, evaluating how the distribution of swift-flowing habitats affects the abundance and dynamics of subpopulations is essential. Using extensive electrofishing data, we show that merging island biogeography with meta-population theory, where lotic habitats are considered as islands in a lentic matrix, can explain spatio-temporal variation in occurrence and density of brown trout (Salmo trutta). Subpopulations in larger and less isolated lotic habitat patches had higher average densities and smaller between-year density fluctuations. Larger lotic habitat patches also had a lower predicted risk of excessive zero-catches, indicative of lower extinction risk. Trout density further increased with distance from the edge of adjacent lentic habitats with predator (Esox lucius) presence, suggesting that edge- and matrix-related mortality contributes to the observed patterns. These results can inform the prioritization of sites for habitat restoration, dam removal and reintroduction by highlighting the role of suitable habitat size and connectivity in population abundance and stability for riverine fish populations.

Two classes of functional connectivity in dynamical processes in networks

The relationship between network structure and dynamics is one of the most extensively investigated problems in the theory of complex systems of recent years. Understanding this relationship is of relevance to a range of disciplines-from neuroscience to geomorphology. A major strategy of investigating this relationship is the quantitative comparison of a representation of network architecture (structural connectivity, SC) with a (network) representation of the dynamics (functional connectivity, FC). Here, we show that one can distinguish two classes of functional connectivity-one based on simultaneous activity (co-activity) of nodes, the other based on sequential activity of nodes. We delineate these two classes in different categories of dynamical processes-excitations, regular and chaotic oscillators-and provide examples for SC/FC correlations of both classes in each of these models. We expand the theoretical view of the SC/FC relationships, with conceptual instances of the SC and the two classes of FC for various application scenarios in geomorphology, ecology, systems biology, neuroscience and socio-ecological systems. Seeing the organisation of dynamical processes in a network either as governed by co-activity or by sequential activity allows us to bring some order in the myriad of observations relating structure and function of complex networks.

Drought effects on invertebrate metapopulation dynamics and quasi-extinction risk in an intermittent river network

Ecological communities can remain stable in the face of disturbance if their constituent species have different resistance and resilience strategies. In turn, local stability scales up regionally if heterogeneous landscapes maintain spatial asynchrony across discrete populations-but not if large-scale stressors synchronize environmental conditions and biological responses. Here, we hypothesized that droughts could drastically decrease the stability of invertebrate metapopulations both by filtering out poorly adapted species locally, and by synchronizing their dynamics across a river network. We tested this hypothesis via multivariate autoregressive state-space (MARSS) models on spatially replicated, long-term data describing aquatic invertebrate communities and hydrological conditions in a set of temperate, lowland streams subject to seasonal and supraseasonal drying events. This quantitative approach allowed us to assess the influence of local (flow magnitude) and network-scale (hydrological connectivity) drivers on invertebrate long-term trajectories, and to simulate near-future responses to a range of drought scenarios. We found that fluctuations in species abundances were heterogeneous across communities and driven by a combination of hydrological and stochastic drivers. Among metapopulations, increasing extent of dry reaches reduced the abundance of functional groups with low resistance or resilience capacities (i.e. low ability to persist in situ or recolonize from elsewhere, respectively). Our simulations revealed that metapopulation quasi-extinction risk for taxa vulnerable to drought increased exponentially as flowing habitats contracted within the river network, whereas the risk for taxa with resistance and resilience traits remained stable. Our results suggest that drought can be a synchronizing agent in riverscapes, potentially leading to regional quasi-extinction of species with lower resistance and resilience abilities. Better recognition of drought-driven synchronization may increase realism in species extinction forecasts as hydroclimatic extremes continue to intensify worldwide.

Network ecology in dynamic landscapes

Network ecology is an emerging field that allows researchers to conceptualize and analyse ecological networks and their dynamics. Here, we focus on the dynamics of ecological networks in response to environmental changes. Specifically, we formalize how network topologies constrain the dynamics of ecological systems into a unifying framework in network ecology that we refer to as the 'ecological network dynamics framework'. This framework stresses that the interplay between species interaction networks and the spatial layout of habitat patches is key to identifying which network properties (number and weights of nodes and links) and trade-offs among them are needed to maintain species interactions in dynamic landscapes. We conclude that to be functional, ecological networks should be scaled according to species dispersal abilities in response to landscape heterogeneity. Determining how such effective ecological networks change through space and time can help reveal their complex dynamics in a changing world.