Oikos

Agroecosystems are dynamic ecosystems, constituted by patches of vegetation and agricultural use, where biodiversity is shaped by spatial and temporal variability. While most studies have focused on spatial composition and configuration, the role of temporal variability remains poorly understood. Yet, temporal dynamics can strongly modify species composition, abundance, and persistence in ecological communities. Temporal variability is particularly relevant in agroecosystems with rainfed agriculture where environmental conditions shift dramatically between rainy and dry seasons. In this paper, we assess the role of temporal variability on biodiversity maintenance in an agricultural matrix using a metacommunity model that simulates an agricultural landscape under rainfed conditions, that is, with abrupt seasonal changes in the agricultural patches. This model couples a local community network dynamic with a migration dynamic and is based on empirically documented features of rainfed agricultural matrices. Our results show that temporal variability provides new opportunities for species to recover from low densities. However, the effect of temporal variability is not straightforward. It depends on the initial and final conditions, the migration and mortality rates and the intensity of temporal variability. Overall, our findings highlight the need to further investigate temporal variability to better understand its role in shaping biodiversity in agricultural landscapes.

One of the most intuitive ideas in ecology is that diversity at lower trophic levels in food webs provides niches to support diversity at higher trophic levels. This accumulation of diversity can be limited by survival of species in the landscape, but revealing these limits has been challenging. We analyze spatial turnover in a hyperdiverse parasitoid-caterpillar-plant food web across 75,000 km2 of continuous lowland rainforest in Papua New Guinea. Species turnover across sites is higher in parasitoids than in their caterpillar hosts. Furthermore, turnover of interactions is also higher in parasitoid-caterpillar than caterpillar-plant networks. Spatial turnover thus amplifies upwards across trophic levels, forcing parasitoids to live closer to spatial persistence limits. Consequently, progressing rainforest fragmentation can especially endanger parasitoids.

The diversity of species responses to environmental change (response diversity) is a key mechanism of ecological stability. However, anticipating where strong or weak stabilizing responses emerge is challenging because species responses can depend on the local community and the specific stability metric. Whether species traits can consistently inform on how species respond to disturbances, enabling less context-dependent predictions, remains an open question. To address this gap, we use microcosm experiments on marine phytoplankton to test how response diversity supports multiple aspects of community stability under pulse temperature changes, testing both an increase (heatwave) and a decrease in temperature (coldspell). We then map species traits to their responses in a community to identify which traits modulate and predict species' sensitivities. Fundamental response diversity, based on the diversity of species responses to temperature measured in isolation, was a weak predictor of community stability, and relationships differed between disturbances (i.e., heatwave and coldspell). Instead, species traits were consistent predictors of species responses in communities. Small, fast-growing species were more tolerant and benefited from the disturbance, while large, slow-growing species were less tolerant and decreased in proportion - these patterns were consistent across disturbances and community compositions. These results suggest that strong trait-performance relationships might reduce the importance of response diversity for stability. But these findings also show that general species traits, such as size and growth rate, can predict which species, and how, contribute to community responses, providing an empirical basis to relate species traits to stability outcomes under climate change.

Under complex perturbation regimes, biodiversity dynamics show temporal variability in species and community abundance around long-term population trends. Many species indeed show long-term declines while other species increase, putting natural communities far from stationary regimes, while variability is often studied near equilibrium. We contribute to bridging this gap by investigating population and community variability during long-term trends caused by press perturbations in stochastic models of population dynamics. By estimating the deterministic changes in mean and variance during the transient regime, we show that population variability deviates from stationary expectations. Moreover, the deviation strongly depends on the sign of the population trends: increases generate excesses of variability while declines generate deficits. Scaling up to community variability, we propose a decomposition of community variability deviation, allowing to highlight that community variability in the transient regime depends on how the press perturbation is distributed within species relative abundances and growth rates. These results challenge the equilibrium assumption and open new perspectives for the study of the variability of ecological systems under multiple perturbation types.

In order to feed a growing global population without silencing nature, conceiving agricultural management strategies reconciling yield and conservation goals is key. Using numerical simulations of a metacommunity model, we explore the possibilities for compromise offered by spatial management strategies of farmed areas. Each strategy is characterized by its farming intensity, the proportion of farmed lands and their spatial aggregation. We show that achieving equitable yield-biodiversity compromise is difficult. While conciliatory strategies offering top yield and biodiversity are typically not possible, accepting slightly lower yields (ie, "Pretty Good Yield strategies") allows to recover substantial biodiversity. Such reconciliation possibilities are limited for species with small dispersal. Yield increases mainly through farmland expansion, whereas farming intensity strongly influences biodiversity, increasing it at low intensity before decreasing with further intensification. Finally, we demonstrate that reconciliation is easier if agricultural production relies on biodiversity through ecosystem services.

Understanding how habitat connectivity shapes biodiversity remains a major ecological challenge. In particular, the roles of connectivity and ecological heterogeneity on co-variation in plant species diversity and intraspecific genetic diversity is not understood. We combined species distribution modelling, resistance-to-movement mapping, landscape connectivity analysis and population genomics to investigate diversity patterns in three wet meadow herbs, Scorzonera humilis, Oenanthe peucedanifolia and Lychnis flos-cuculi, and their surrounding plant communities. Genetic diversity patterns differed strongly among co-occurring species. Connectivity metrics explained genetic diversity only in O. peucedanifolia, and environmental drivers of genetic diversity were highly species specific. Importantly, genetic diversity changed with the presence of some species in the community, but it was consistently unrelated to indicators of local plant community diversity. Overall, the processes shaping within-species biodiversity may differ fundamentally from those structuring habitat connectivity and plant species communities, with important implications for conservation.

The impacts of a changing abiotic environment on fitness and performance arise not only from low tolerance to new environmental conditions, but also from changes in the abundance and interaction intensity with other species. The strength of the interaction may itself depend on how well each species performs across environments, but there is a dearth of studies investigating how intrinsic fitness and interaction intensity covary across an abiotic environmental gradient. We addressed this question in a hypersaline consumer-resource system: the microalga Dunaliella spp. grazed by the brine shrimp Artemia franciscana. We exposed four Dunaliella strains to a range of salinities above seawater, with or without brine shrimps, and tracked their population sizes over time and the survival of their predators, to estimate basic parameters of a Lotka-Volterra model. We found that the intrinsic growth rate of algae, the survival rate of predators, and the per-capita predation rate, all varied with salinity and algal strain. Significant interactions between strain and salinity further revealed that these ecological responses to salinity are evolvable. Together with correlations between demographic parameters across salinity, this suggests that predation may influence the evolution of salinity tolerance curves, blurring the line between the fundamental and realized niches.

Protected areas cannot halt biodiversity loss in isolation; integrating them with surrounding human-dominated landscapes is critical. However, this integration is challenged by substantial landscape heterogeneity at their borders, hindering our understanding of cross-border changes in ecosystem service provision. We introduce a novel framework for characterizing these dynamics by analyzing ecosystem service gradients along protected area borders. For 16 protected areas in the French Alps, we assessed 12 ecosystem services using a mix of established biophysical models and novel connectivity-based models for mobile species. These were aggregated into three stakeholder-driven domains reflecting respectively rural, cultural, and urban management priorities. Automated polynomial regression analysis classified borders into five gradient types. The most common were 'Decreasing Gradients', representing a decline in ecosystem services outside the protected area, and 'Increasing Gradients', with the opposite pattern. Our analysis reveals these patterns are driven by specific landscape configurations, uncovering frequent trade-offs between the three management priorities, where, for instance, landscapes supporting rural priorities often degrade cultural and urban ones. We also identify key opportunities for synergies, by identifying areas where ecosystem services for all three priority domains increase simultaneously outside the protected area. This spatially explicit typology provides a powerful diagnostic tool for designing targeted interventions, such as prioritizing habitat restoration where ecosystem services decline or managing agricultural landscapes to mitigate conflicts across management priorities, supporting a more effective integration of protected areas into the wider landscape.

Climate change is expected to alter plant populations not only through direct environmental shifts but also via changes in biotic interactions, such as with herbivores and pathogens. As plant species are also expected to differ in their responses to both climate and antagonists, plant responses to both factors are expected to be variable and species-specific. To assess whether interactive effects of climate and antagonists on plant population dynamics are common and whether the strength and direction of plant responses vary across species, we conducted a multi-year field experiment that manipulated realistic climate change and experimentally reduced insect herbivores and fungal pathogens. We measured responses of plant vital rates, such as survivorship, growth, and reproduction across six grassland species. Using Integral Projection Models (IPMs) and Life Table Response Experiments (LTREs), we quantified changes in population growth rate across experimental treatments and the contribution of each vital rate to that observed change. Two of the study species declined so drastically over the course of the experiment that demographic quantification of population growth rates was not possible. From the remaining species, Bromus erectus and Plantago lanceolata show significant interactive responses of climate and antagonist reduction on population growth rates. In contrast, Dianthus carthusianorum and Tragopogon orientalis showed limited responses to experimental treatments. Notably, our results indicate that in some species biotic interactions may amplify the effects of climate change: the presence of plant antagonists exacerbates the negative effects of the future climate treatment on plant population dynamics. Our findings highlight the complexity in predicting plant population responses to climate change and provide insights for grassland management under future environmental conditions.

Despite the ecological importance of ungulate migrations, we lack a complete understanding of why some ungulates migrate and others do not. Though progress has been made towards understanding differences across species and between populations, migratory behavior varies even within populations: in many populations, some individuals remain behind as residents (partial migration). Theoretical population-level work has suggested that these different migratory tactics can coexist, but such approaches stop short of providing insights into how individuals make the decision to stay or go each year. Using long-term data from three ungulate populations, we find that individuals probabilities of migrating are highly variable across years, which points to a non-trivial context-dependent decision-making process, whose underlying mechanisms must be probed via individual-level modeling. Drawing on existing knowledge, we propose a decision-making model of ungulate migration onset wherein individuals probabilistically decide to start migrating based on the local intensity of environmental and/or social cues. Residents arise as a robust collective organization phenomenon in our model. At sufficiently large population sizes, the number of residents is invariant with total population size, consistent with empirical patterns. Instead, resident numbers are influenced by the severity of the bad season, by relevant character differences among individuals, and by how individuals contribute and respond to environmental and/or social cues; for instance, when social cues contribute to decision-making in addition to environmental ones, fewer residents result, and migration is more likely to be complete. Overall, our model provides a potential mechanistic explanation for how residents might emerge within migratory ungulate populations.

1. Ruby-throated hummingbirds (Archilochus colubris) exhibit flower preferences and readily visit human subsidies like feeders. However, foraging behavior and plant-animal interactions may vary across spatial, temporal, and landscape gradients. Despite the popularity and ubiquity of Ruby-throated Hummingbirds in the eastern U.S., there has never been a quantitative assessment of their flower preferences or feeder use across broad scales. 2. We investigated how feeder visitation, flower visitation, and flower trait preferences vary by latitude, season, and land use in the northeastern United States using over 6 million occurrence records of Ruby-throated Hummingbirds and flowering plants, >2,100 annotated hummingbird-flower interactions, and >2,700 feeder visit occurrences. 3. We found that hummingbird feeder use declined over the year, increased with latitude, and was higher in developed landscapes. Flower visitation increased over the year across all latitudes, with higher visitation in developed landscapes. Finally, we found that native plant use diverged between landscapes, such that the probability of visiting a native flower increased over time in non-developed land uses but declined over time in developed ones, demonstrating that hummingbirds track the advancement of native floral phenology and use non-native, cultivated flowers as a human subsidy due to either availability or preferences. 4. Our preference and network models revealed that while hummingbird-plant network structure was similar across landscapes, the composition of important taxa shifted from native, wild species like Monarda and Impatiens to non-native, cultivated species like Salvia. 5. Using trait-based models of flower visitation, we found that hummingbirds preferred native, tubular, and red/orange flowers fitting the hummingbird pollination syndrome despite visiting >260 different plant species. Red and orange flowers were preferred across all seasons, suggesting color may be a reliable signal of nectar availability across species and contexts. Native and tubular flowers were strongly preferred during the breeding season; however, preferences relaxed during spring and fall migration. 6. These findings reveal the consistent preferences of Ruby-throated Hummingbirds for native, tubular, and red/orange flowers, and underscore how spatial and temporal factors reshape foraging behavior and trait preferences. Our results also highlight the value of community-collected data in characterizing plant-pollinator interactions across broad spatial and temporal scales.

Environmental variability can strongly alter coexistence among competing species and their extinction risk, particularly when population dynamics are shaped by behavioral interactions, such as fear. In this work, we develop a novel stochastic differential equation competition model that incorporates both non-consumptive fear effects and environmental variability to investigate how behavioral interactions influence species coexistence under random fluctuations. Our result reveals that environmental stochasticity can drive species to extinction even when the corresponding deterministic system admits coexistence. In particular, under an explicit stability condition on the fear and competition parameters and sufficiently strong averaged noise intensities, we prove that both competing species become extinct exponentially almost surely. Conversely, we derive a stochastic persistence criterion in terms of fear, competition, and noise-induced suppression parameters for the fearful species. We further demonstrate that environmental noise may reverse classical competition-exclusion outcomes, leading to qualitatively different long-term dynamics from those predicted deterministically. These results provide rigorous thresholds separating stochastic extinction from persistence and highlight the critical role of environmental variability in fear-mediated competitive ecosystems. From an applied perspective, these results provide insight into how behavioral interactions and environmental variability influence species survival, with potential applications in ecological management and conservation.

AO_SCPLOWBSTRACTC_SCPLOWO_LICommunity assembly graphs (CAGs) summarize which species combinations can coexist and how single-species invasions drive transitions between them, encoding the pathways, alternative endpoints, and cycles that make up a communitys assembly history. Constructing CAGs from dynamical models requires methods that are both computationally tractable and faithful to the underlying ecological dynamics. However, existing methods rely on restrictive assumptions, such as global stability, that exclude alternative stable states and non-equilibrium dynamics known to occur in empirical systems. C_LIO_LIWe develop a computational pipeline that constructs CAGs from any generalized Lotka-Volterra model. Building on the invasion graph framework and its connection to permanence, the pipeline verifies that community dynamics are bounded, identifies which subsets of species coexist in the sense of permanence, determines which single-species invasions are dynamically realized, and assigns each community a topographic height equal to the length of the longest assembly path leading to it. We also provide a numerical algorithm to simulate the dynamics of community assembly. C_LIO_LIWe prove several general properties of the resulting graphs, including that a successful invader is never subsequently excluded and that, in the absence of assembly cycles, permanent communities can be reassembled by introducing their species one at a time in the right order. We prove that the CAG faithfully reproduces the compositional shifts seen in the numerically simulated dynamics of assembly. Applying the pipeline to three empirically based models (a New Zealand grassland, a European pasture, and a Puerto Rican ant community), we show how competition strength and mutualistic feedbacks reshape the assembly landscape and how intransitive competition generates assembly cycles. C_LIO_LIOur approach accommodates alternative stable states and non-equilibrium dynamics without requiring global stability, and it turns the long-standing landscape metaphor into a quantitative, mechanistically grounded object by resolving what "height" means. More broadly, it makes the topography of the assembly pathways measurable, providing a way to compare the historical contingency and predictability of the assembly in ecological systems. C_LI

This study investigates changes in habitat connectivity and meta-ecosystem resilience between 1817 and 2022 along a 150 km section of a large European river (Danube) and its adjacent floodplains. The analysis was based on a time series of functional habitat networks (graphs) constructed from historical records and remote sensing data, integrating habitat suitability and dispersal modes of functional organism groups. The results indicate that overall habitat availability declined by about 50% since 1817, leading to the near-complete loss of functional connectivity among dynamic habitats by 1910. Less dynamic habitats persisted or expanded but became functionally less connected. Connectivity-based habitat classification further revealed four distinct functional connectivity clusters and the near-complete loss of an originally dominant, dynamic, and highly connected habitat type. Furthermore, the results indicated a fundamental loss of meta-ecosystem resilience. This development was reflected by increased spatial modularity among functionally similar habitats and by reduced layer dissimilarity and structural robustness in multilayer networks representing the spatial habitat structure and functional habitat connectivity across different functional organism groups.

The exchange of energy and organisms across habitat boundaries links aquatic and terrestrial ecosystems and sustains ecosystem functioning. Although disturbance may disrupt these linkages, the mechanisms at play remain poorly understood. Here, we investigated the extent to which flow intermittency may disrupt riparian-aquatic ecosystem linkages by altering consumer communities in the recipient ecosystem or by altering resource quality in the donor ecosystem. We ran an experiment in an intermittent river network in California, focusing on a critical forest-to-river subsidy (organic matter in the form of leaf litter), its transformation, and its reciprocal benefit (aquatic insect production). Using three riparian species (willow, cottonwood, and oak) at sites spanning a gradient of flow permanence, we quantified intraspecific plasticity in leaf traits (specific leaf area, nitrogen and phosphorus concentrations, and {delta}13C), measured decomposition rates, and estimated the secondary production of aquatic shredders (Plecoptera). Across all leaf species, decomposition rates were 16-36% lower at intermittent than perennial sites, an effect largely driven by intraspecific leaf trait plasticity rather than changes in consumer abundance. At high flow intermittency, willow experienced water stress (enriched {delta}13C) and reduced specific leaf area, while cottonwood showed primarily stoichiometric responses (reduced leaf nitrogen and phosphorus). Despite these divergent strategies, all species produced lower-quality litter at intermittent sites. Variance partitioning confirmed that initial litter quality uniquely explained 51.5% of variation in decomposition rates, more than double the contribution of invertebrate community metrics; and structural equation modeling revealed that both leaf traits and stonefly (Plecoptera) secondary production significantly predicted decomposition rates, with leaf traits exerting the stronger effect. Notably, stonefly secondary production was 37-98% lower at intermittent sites across leaf species. Because these insects later emerge as terrestrial adults, they provide a significant energy flux to riparian predators, and, thus, impoverished litter quality suppresses the reciprocal transfer of energy back to terrestrial food webs. As drought intensifies globally, the decoupling of terrestrial-aquatic linkages may begin in the riparian canopy.

As the worlds oceans change in response to climate change, phytoplankton communities will adapt to warmer, more stratified surface waters via plasticity, evolution, and range shifts. Current global ocean models assume that size structured phytoplankton communities have fixed trait relationships, and as a result generally predict that smaller size classes will become more dominant globally. However, this general expectation fails to consider how intra-species trait tradeoffs may operate orthogonally from large-scale inter-species tradeoffs--allowing for alternative evolutionary pathways given the limits and/or possibilities available to ancestral populations. To identify evolutionary pathways phytoplankton populations might take, we develop a novel modeling framework that combines a trait-based phytoplankton quota model with stochastic evolution (ecoTRACE). EcoTRACE explicitly decouples key phytoplankton traits from interspecific allometric relationships, allowing for novel phenotypes to emerge. We validated ecoTRACE against a long-term artificial size selection experiment on Dunaliella tertiolecta. We show that ecoTRACE captures multi-dimensional evolved phenotypes that quota models based on interspecific relationships fail to reproduce. Under fluctuating multi-stressor growth, model populations evolve phenotypic plasticity that deviates from predicted interspecific allometric relationships. EcoTRACE provides a framework for generating hypotheses as to the evolutionary trajectories that phytoplankton will experience in a warmer, more variable ocean.

Arctic and Boreal raptor communities will continue to be affected by borealization and other climate change related processes, providing a challenge for ecologists predicting future sates. However, by using community assembly theory and species traits, future communities may be predictable. In this study, we analyzed variation in reproduction traits as a consequence of diet specialization for 29 raptors, 2 skuas and 3 corvids. We assessed and implemented foraging traits for specialists and generalists into predator-prey models from which successful invasion conditions were derived. Specialist raptors produced larger clutch sizes, had a higher proportion of fledged per clutch and also expressed more variation compared to generalist raptors. These results suggest a relationship between diet specialization and reproductive traits which was also observed within phylogenetic orders. Specialist owls (Strigiformes) produced higher clutch sizes with a larger clutch range compared to generalist owls. The same pattern was observed for falcons (Falconiformes). No clear difference in reproduction was observed for specialist and generalist hawks, kites and eagles (Accipitriformes). Corvids expressed clutch sizes similar to that of specialist raptors while having the lowest proportion of fledged per clutch. Differences in foraging traits between specialists and generalists could be distinguished using functional response curves. A predator-prey model parameterized with foraging trait data showed that a generalist can coexist with a resident specialist if it has access to prey unavailable to the resident specialist. Otherwise, the native specialist outcompetes the invading generalist due to foraging efficiency. The combined empirical and theoretical findings in this study show how diet specialization affects both reproduction and the potential invasion success of raptors.

Understanding the relationships between soil macrofauna and decomposition is crucial for predicting how land-use change impacts ecosystem function in fragmented systems. This is because soil macrofauna affect decomposition and also respond to the changes in abiotic conditions across habitat gradients. This study investigates edge effects on the macrofauna contributions to decomposition across forest-field ecotones. We used bait lamina assay to quantify aboveground and belowground feeding activity of soil macrofauna in Autumn 2025 in three deciduous forest-old field ecotones and one coniferous forest-old field ecotone, in Southwestern Ontario, Canada. Vegetation diversity and composition, LAI and soil characteristics (i.e., soil organic matter, pH, temperature and moisture) were measured at each plot along the ecotone. Pitfall trap data collected in Summer 2025 at the same sites were used to characterize macrofauna communities. We used generalized linear mixed effects models to estimate the effect of distance to edge, site, and depth into the soil on bait lamina consumption and soil macrofauna, with transect nested within site as random effects. Consumption activity increased with distance into the forest from the field, with the edge representing an intermediate; and, decreased with increasing depth into the soil. In contrast, soil macrofauna abundance, especially isopods, decrease with distance into the forest from the field. These trends varied significantly across sites, so that consumption activity and abundance sometimes remained constant across the ecotone (i.e., site x distance interaction). The results demonstrate that macrofaunal contributions to bait consumption varied along the ecotone, shaped by interacting environmental gradients and shifts in community composition unique to each site.

Urban green spaces are increasingly recognised as important refuges for biodiversity, yet their ecological value depends strongly on design and management. Here, we investigate how fine-scale structural and microhabitat components shape urban ant assemblages, using ants as indicators of broader arthropod responses to urbanisation. Ant communities were sampled in twelve urban green spaces in Cordoba (southern Spain) over a ten-year period (2004 to 2013) using pitfall traps, alongside detailed characterisation of vegetation structure and ground-layer microhabitats. In total, 38 species and 25,578 individuals were recorded. Microhabitat variables explained 58% of the variation in species occurrence. Community differences among microhabitats were driven primarily by nestedness, with dense herbaceous cover acting as a core habitat and edge-related components contributing disproportionately to beta diversity. Tree abundance showed a unimodal relationship with species richness, with maximum diversity at intermediate densities, while shrub and lawn cover had weak or inconsistent effects. Fine-scale elements such as leaf litter, stones, woody debris, and small bare-ground patches strongly influenced species occurrence by providing thermal refugia, nesting substrates, and foraging opportunities. The invasive Argentine ant (Linepithema humile) exhibited strong but spatially restricted dominance and species-specific negative effects on native ants, emphasising the role of habitat context in mediating invasion impacts. Our results demonstrate that urban biodiversity is maximised by enhancing fine-scale habitat heterogeneity rather than increasing green cover alone. We highlight practical design principles for urban green infrastructure that prioritise structural diversity and ground-layer complexity to support resilient arthropod communities.

Artificial light at night (ALAN) can alter the movement of nocturnal pollinators such as moths, with potential consequences for pollination services. However, most studies have focused on the presence or properties of lights, while the role of spatial lighting configuration remains poorly understood. We conducted a small-scale field experiment to test how different lighting configurations affect pollination success in moth-pollinated plants, using low-intensity LED garden lights. Potted phytometer plants of three moth-pollinated species were exposed to one of three treatments: no-light control, isolated LED point lights, or 25 m linear arrays of multiple LED lights. We quantified pollination success as both the probability of seed set and reproductive output through seed and capsule production. Experimental lighting affected plant reproductive success, with light arrays increasing seed-set probability, seed and seed capsule number relative to unlit controls, whereas point lights showed no clear effect. These results provide preliminary evidence that the spatial configuration of artificial lights may influence nocturnal pollination outcomes. Future work combining phytometer assays with direct tracking of moth movement is needed to assess whether light arrays facilitate or redirect the dispersal of nocturnal pollinators.