Ecosphere

Ecologists remotely track movement steps of animals (e.g., via global positioning systems) and use step selection functions to study the effect of environmental factors upon their movement decisions. Constructing such functions requires pairing each observed step with some number of unobserved but feasible comparison steps. Larger numbers of comparison steps generally yield better estimates but also incur potentially challenging computational demands. Thus, it is important to determine an appropriate number of comparison steps. No established guidance exists for this decision. Here, we use simulated tracks to assess how many comparison steps are needed, fitting each set of steps to a conditional logistic regression model. We monitor errors in estimated effects for several classes of tracks, identifying the number of comparison steps for which mean relative absolute error in estimated effects is consistently low. By this criterion, 32 comparison steps per observed step are needed for our primary class of simulated tracks. Tracks in more homogeneous landscapes, tracks with shorter mean step lengths, or shorter tracks generally require more comparison steps (ranging from 64 to 128 per observed step) to achieve the same level of accuracy. Longer tracks generally require fewer comparison steps (16 per observed step). These results clearly demonstrate that the number of comparison steps influences how well step selection functions estimate covariate effects and provides initial direction in a research area that currently lacks quantitative guidance. Movement ecologists should take care when selecting the number of comparison steps paired with each observed step because those decisions matter.

Many Sockeye salmon (Oncorhynchus nerka) populations have declined over recent decades, and climate change is likely to exacerbate these declines through direct and indirect ecological effects. The response to the associated environmental changes is likely to vary among life stages, populations, and regions. Quantitative estimates of climate change driven impacts that account for this variability could fill a critical gap and provide forward-looking insights into how sockeye are expected to respond to future climate-driven change across their lifecycle. To address this need we developed a hierarchical population dynamics model parameterized with juvenile, adult return and spawner abundance data from 13 sockeye salmon populations from Washington State to northern British Columbia. We used a formal causal inference framework that paired salmon abundance data with a suite of environmental covariates hypothesized to represent ecological conditions across the lifecycle. We used the model to estimate population-specific responses to each environmental driver, then combined parameter estimates with projections from down-scaled climate change models to estimate productivity responses to anticipated environmental change. We found that historical sockeye productivity was strongly associated with environmental covariates, which explained more interannual variability in return abundance than spawner abundance in most populations. However, the life stages and specific environmental covariates with the largest impacts differed among populations and regions, often displaying a latitudinal gradient. Increases in coastal ocean temperatures and mixed layer depth generally had negative effects though they varied among regions. Increased freshwater summer rearing and return migration temperatures had weaker but consistently negative effects. Under future climate conditions, projected changes in these environmental covariates are expected to result in substantial declines in productivity across most populations. Sockeye salmon display varying degrees of sensitivity to climate change across life stages, populations, and regions. Effective future management will require explicitly accounting for these life stage and population-specific responses.

The reintroduction of extinct or endangered species to restore ecosystem function is an essential aspect of rewilding. The Wilder Blean Project at West Blean and Thornden Woods in Canterbury, UK, is committed to rewilding natural processes and enhancing biodiversity in one of England's oldest and largest areas of ancient woodland. The introduction of European bison (Bison bonasus) is an important part of the project. However, how the reintroduction of large herbivores influences local biodiversity and ecosystem functions during the early stages of rewilding remains poorly understood. Soil samples were collected from the same sampling sites before and two years after bison were reintroduced and profiled by metagenomic sequencing using Oxford Nanopore Technologies sequencing platforms. The results showed that the alpha diversity of soil organisms did not change significantly before and after the introduction of European bison, while beta diversity showed modest shifts in community composition. The relative abundance of some nitrogen-fixing and photosynthetic microbial genera showed declines in the 2024 Bison Area, while the mycorrhizal fungus genus Rhizophagus was significantly less abundant than in the 2024 Control Area. Despite relatively stable taxonomic diversity, functional composition differed significantly between the 2022 and 2024 Bison areas and among the 2024 rewilding treatments, revealing a decoupling between taxonomic diversity and functional composition. Amino acid synthesis pathways and carbon metabolism pathways were significantly enriched. These findings highlight the potential of long-read Oxford Nanopore metagenomics to reveal functional shifts that may not be apparent from taxonomic diversity alone. Although these early-stage responses cannot yet predict long-term rewilding trajectories, continued longitudinal monitoring integrating microbial, soil physicochemical, and ecosystem-level measurements will be essential to determine the persistence and ecological significance of these functional shifts.

Replication is central to most experimental and sampling designs, increasing inferential power and capturing fine-scale data heterogeneity. However, its importance remains poorly evaluated in some ecological and evolutionary settings. This is the case of metabarcoding studies using DNA recovered from sedimentary archives, in which biological signals may integrate ecological information through depositional and burial processes, and are often inferred from a single sediment core per site. Here, we evaluated the effect of different types of replication using sedimentary DNA (sedaDNA) metabarcoding data from two genetic markers (mitochondrial COI and nuclear 18S), under a nested sampling design. The design included three intertidal sites, three spatially separated sediment cores per site (biological replicates), two sediment depth horizons per core, and eight PCR (technical) replicates per sediment sample. Variance partitioning showed that site identity and sediment age group together explained >70% of the variation in beta diversity, indicating that among-site spatial variation and stratigraphic variation were the dominant drivers of community composition. In contrast, variation among different cores within sites was small and non-significant (<5%). Among PCR replicates from the same sediment sample, richness varied substantially, whereas Shannon diversity was more consistent. Despite this variability, differences in community composition among technical replicates remained smaller than among biological replicates and site identity, indicating limited influence on broader ecological patterns. Community composition was highly similar among replicate cores within sites, consistent with stratigraphic coherence. These results indicate limited within-site heterogeneity and suggest that, under stratigraphically coherent conditions, increasing biological replication may yield limited additional information, whereas enhancing technical replication and stratigraphic resolution can improve ecological inference from sedaDNA metabarcoding datasets.

Root-associated fungal (RAF) communities can influence tree nutrient acquisition and performance, yet their structuring factors and relationship with tree growth remain poorly understood, particularly near environmental treelines. We characterized root-associated fungal (RAF) communities on fine roots of white spruce (Picea glauca [Moench] Voss) in paired forest and treeline plots across two elevational and one moisture-limited treeline ecotones in Alaska. Using ITS2 DNA metabarcoding of fine root tips, we characterized individual RAF community alpha and beta diversity and tree growth based on basal area increment (BAI) over 5-30 years. As a result, sampling site was the strongest predictor of RAF composition, explaining 19.6% of variation, while soil pH explained an additional 11.7%. Treeline effects on RAF communities were weaker and context-dependent, with treeline trees showing distinct community composition in the Alaska Range, lower alpha diversity across all Hill numbers in Interior Alaska, and higher ECM relative abundance in the Brooks Range. RAF composition did not differ between fast- and slow-growing trees within sites. In contrast, alpha diversity was negatively associated with tree growth: OTU richness, Shannon and Inverted Simpson indices significantly predicted BAI over the previous 5, 10 and 15 years, with fast-growing trees supporting less diverse RAF and ECM communities. This relationship was strongest for recent growth and weakened with longer BAI averaging periods. Our results suggest that high growth of mature P. glauca is not necessarily associated with greater RAF or ECM diversity or specific taxa but may instead rely on fewer dominant or functionally effective fungal partners.

Understanding the spatial distribution of species and their patterns of endemism is necessary for establishing effective conservation priorities. Despite the vital pollination services bees provide, Californias native bee distribution patterns remain largely unexplored relative to plants and butterflies. We analyzed bee species richness and endemism across California and their concordance with plant distributions. Richness was high across areas of the California Floristic Province, including the Sierra Nevada, San Francisco Bay Area and Central Coast, South Coast Ranges, and the Transverse and Peninsular ranges. Bee endemism was more localized, concentrated in the San Joaquin Valley, eastern Sierra Nevada and adjacent Great Basin, Sierra Nevada foothills, and California deserts. Because richness and endemism appear to operate at different spatial scales and likely respond to different environmental drivers, effective conservation strategies must address both. Additionally, conservation plans that incorporate both plant and bee diversity are needed to achieve more comprehensive biodiversity protection.

Managers must make decisions in the face of uncertainty, especially when available resources are limiting. Identifying thresholds when certain conditions are met or exceeded enable the potential to mitigate risks. In 2005, sustainable harvest levels of willow ptarmigan were identified to avoid harvest efforts exceeding three hunter days km2. Here we evaluate these recommendations by analyzing line transect counts and harvest data from six areas forming three open/closed pairs in a region of state managed willow ptarmigan harvest. We developed three sets of Bayesian hierarchical models, one static distance model, and two dynamics models. One mechanistic hazard model and a Gompertz phenomenological model. Adult and juvenile density showed pronounced year-to-year variation that was largely synchronous across all six sites regardless of hunting status. The harvest effort parameter shows a striking difference between the two models. In the Hazard model, is positive, and excludes zero with near certainty, but in the Gompertz model, the parameter is highly uncertain. However, the two models do not contradict each other but answer complementary questions with different sensitivity to the harvest signal, harvest mortality is additive at the individual level, but this additive mortality is masked at the level of population abundance. The demographic cost of harvest is therefore real and quantifiable through the survival chain, but bounded in the long run by the stabilizing dynamics. A fixed limit anchored to monitored effort and bag is not a crude substitute for adaptive management but the appropriate design under the information commonly at hand. It will be a precautionary instrument grounded in the one relationship this study establishes firmly, the translation of hunter effort into harvest mortality.

Gray fox (Urocyon cinereoargenteus) populations in the Midwestern USA have suffered precipitous declines in recent decades, yet they are relatively understudied. However, understanding survival and cause-specific mortality is vital for declining populations and the limited existing survival studies have been performed outside of the Midwest. We equipped 13 gray foxes in southern Illinois with GPS radio collars to investigate their survival and cause-specific mortality. We calculated the Kaplan-Meier 6- and 12-month survival rates to be 0.79 (95% CI: 0.57-1.0) and 0.53 (95% CI: 0.27-1.0), respectively. We recorded 4 mortalities: 1 disease, 1 gunshot, and 2 unknown causes. While our study has a small sample size, it contributes key information on a data-deficient mesocarnivore suffering from a population decline driven by undefined causes. We recommend further research into the survival and mortality of this elusive mesocarnivore.

Wildflower plantings are an important conservation measure for supporting wild bee diversity. They aim to enhance floral resources for nutrition, but do not explicitly consider that bees also require suitable nesting and overwintering resources. Wildflower plantings may provide nesting habitat for ground-nesting bees, but the effects of soil management, e.g. ploughing, on ground-nesting bees are hardly known. To study how ploughing and age of wildflower plantings affect ground-nesting bees, we sampled bees on ploughed and unploughed wildflower plantings aged 0-4 years. We used emergence traps to sample bees directly after emergence from the ground, allowing inference on nest numbers. We found that ploughing and wildflower planting age negatively affected overwintering bee nest numbers. Nesting peaked in the year of wildflower planting establishment and declined thereafter, indicating lower habitat suitability at later successional stages. Annual ploughing caused the greatest reduction in nest numbers (-72 %), providing evidence for an ecological trap.

Flow regulation for hydropower production affects stream ecosystems through decreased connectivity and changed timing and magnitude of flows. Hydropeaking, a form of regulation in which infrequent large peak flows are replaced with frequent small peaks, increases riparian erosion and causes water and drought stress for riparian vegetation. Hydropeaking is likely to affect soil seed bank (SSB) formation and composition, while SSBs are important sources of self-restoration should a systems flow regulation be relaxed. We tested how hydropeaking intensity affects the size and composition of SSBs, including the functionally important group of large graminoids, and by calculating Ellenberg values for Moisture, Light and Soil disturbance. SSB samples were taken at 15 riparian zones across central and northern Sweden. Each site was regulated, but sites differed in their hydropeaking intensities. SSBs were subjected to a seedling emergence experiment, from which over 700 seedlings from 53 taxa emerged. We found that hydropeaking intensity affects the composition of soil seed banks on multiple levels. Seedling density was negatively correlated with hydropeaking intensity at the sites where samples were taken. SSB richness varied (two to eighteen species per site) and was not affected by hydropeaking intensity. The proportion of large graminoids in the seed bank showed a near-significant decrease with increasing hydropeaking intensity, and community-weighted means for Moisture, Light and Soil disturbance increased (non-significantly) with increasing intensity. Our results suggest that riparian SSBs, should flow regulation be relaxed or ceased, are not sufficient for self-restoration of functional riparian vegetation. Seeds of large graminoids and species that are tolerant to drought in the germination stage are less present in riparian SSBs of heavily-regulated streams. Supply of seeds of these groups, or even planting, may need to be considered when changes in flow management are implemented. HighlightsO_LIHydropeaking negatively affects riparian soil seed banks (SSBs) in Sweden C_LIO_LISSB size slightly decreases with hydropeaking intensity, but richness does not change C_LIO_LIThe proportion of large graminoid seeds in SSBs decreases with hydropeaking intensity C_LIO_LIRiparian SSBs from less-impacted sites have most potential for self-restoration C_LI

1. Natal dispersal is a key process in population ecology because it links local demographic processes to broader-scale population dynamics by redistributing individuals. When using capture-recapture data, multistate capture-recapture models using discrete spatial units as states are the gold standard for estimating natal dispersal among spatial units while accounting for spatial variation in survival, recruitment and imperfect detection. However, because their computational cost increases rapidly with the number of spatial units, applications have been limited to a small number of units. Therefore, in practice, these models cannot provide spatially detailed inference on natal dispersal across large landscapes. 2. We develop a computationally efficient Bayesian capture-recapture model, called the efficient natal dispersal (END) model, to estimate natal dispersal among discrete spatial units jointly with spatial variation in demographic parameters and detection probabilities. The END model relies on two key structural features: juveniles and breeders are separated into two arrays, and resightings outside the natal spatial unit are aggregated over time for individuals released as juveniles. 3. Using simulations, we show that the END model is considerably (up to 30 times) more computationally efficient than a conventional multistate model, while maintaining comparable parameter accuracy. We then apply the END model to white stork (Ciconia ciconia) capture-recapture data from Germany across 101 hexagonal spatial units, a spatial resolution at which a conventional multistate model is computationally infeasible. We estimate natal dispersal among units jointly with spatial variation in survival and recruitment. This allows us to identify areas of lower or higher survival, earlier or delayed recruitment, and dispersal probabilities among all units. By combining estimated dispersal probabilities with existing data on the number of juveniles born in each spatial unit, we estimate natal dispersal in terms of numbers of individuals and identify units with positive or negative net migration, sources and sinks. 4. Overall, our approach moves capture-recapture analyses from estimating natal dispersal among a few spatial units to inferring dispersal networks and assessing their demographic consequences across large domains. Our approach is applicable to many spatially structured capture-recapture datasets, opening new opportunities for studying spatial population dynamics.

Freshwater bioassessment relies on assessing aquatic assemblages to infer ecological conditions, yet conventional surveys require extensive field sampling, specimen processing, and specialized taxonomic expertise. Existing environmental DNA (eDNA) methods have not yet provided a practical alternative to conventional macroinvertebrate assays in part because current approaches cannot feasibly recover broad taxonomic diversity at sufficient taxonomic resolution. Here, we evaluated targeted hybridization capture of mitochondrial cytochrome oxidase I (COI) target sequences as a unified molecular approach for cross-phylum freshwater bioassessment. Environmental DNA was collected at 18 sites along 63 km of Boulder Creek spanning nearly 1,500 m of elevation from forested headwaters to agricultural plains. COI targets were enriched using custom RNA bait panels designed to target regional freshwater arthropods, annelids, and molluscs. Hybridization capture increased recovery of COI sequences [~]1,760-fold relative to unenriched shotgun libraries, generating Folmer-region COI contigs that averaged [~]400 bp. Across the watershed, we recovered sequences for approximately 450 macroinvertebrate genera across 8 phyla. Detected macroinvertebrate richness averaged 56 genera per site and increased down Boulder Canyon before declining downstream of the city. Macroinvertebrate assemblage composition from hybridization capture paralleled patterns observed with past conventional bioassessment. These results demonstrate that targeted hybridization capture enables robust, cross-phylum detection of species used for freshwater bioassessment from environmental DNA.

Understanding how closely related species coexist in highly seasonal and unpredictable environments is central to studies of ecological differentiation and niche partitioning. We investigated the feeding ecology of sympatric populations of Microcebus murinus and M. griseorufus within a contact zone in Andohahela National Park, southeastern Madagascar, across dry and wet seasons. Using a combination of direct behavioral observations (1,611 feeding records), fecal sample analyses (n = 56), and vegetation phenology surveys, we quantified dietary composition, seasonal shifts in resource use, and habitat-related variation. Seasonal changes in diet were pronounced, with dry-season feeding dominated by exudates and wet-season diets incorporating greater proportions of fruit and flowers, closely tracking phenological patterns at both sites. Diets of both species were dominated by plant resources, but consistent interspecific differences in dietary strategy were evident. Although both species consumed comparable proportions of insect prey, M. murinus showed pronounced wet-season increases in the use of high-sugar, carbohydrate-rich floral resources (15.9%) and hemipteran-associated honeydew (30.6%). In contrast, M. griseorufus relied more consistently on predictable exudates throughout the year. Fecal analyses supported observational data but revealed differences in the detectability of dietary components, with increased representation of invertebrates in the wet season and seeds in the dry season. These results indicate substantial dietary overlap but consistent differences in resource use, suggesting that coexistence is facilitated by fine-scale trophic differentiation within a broadly shared omnivorous niche. Such subtle but persistent differences in feeding strategy likely reduce competitive overlap and enable continued sympatry in a climatically variable and resource-limited system.

Investigating the factors that stabilize biological communities is a central topic in ecology. Response diversity, defined as variation in species responses to environmental change, has been proposed as a key mechanism underlying the biodiversity-ecosystem functional stability (BEFS) relationship, whereby greater species diversity enhances ecological stability. Previous studies have shown that response diversity promotes ecological stability by generating asynchronous population fluctuations and the resulting compensatory dynamics. Although several metrics have been proposed to quantify response diversity, they do not explicitly consider the presence of insensitive species whose performance is unaffected by current environmental conditions. To examine how insensitive species influence response diversity, species persistence, and ecological stability, we conducted numerical simulations of a generalized Lotka-Volterra model under environmental forcing. We first confirmed that increasing variation among sensitive species increased the response diversity index and stabilized community dynamics. We then examined a scenario in which response diversity depended solely on the proportion of sensitive and insensitive species, assuming that all sensitive species responded identically to environmental change. Under this assumption, the response diversity index was maximized when sensitive and insensitive species occurred in equal proportions, whereas increasing the number of sensitive species monotonically destabilized community dynamics. Consequently, the relationship between response diversity and community stability depended on how response diversity was generated, such that higher response diversity could even be associated with lower community stability. These findings demonstrate that overlooking environmentally insensitive species can obscure the mechanisms linking response diversity and ecological stability. More broadly, our results reveal that response diversity comprises at least two distinct biological components--species sensitivity and response variation among sensitive species--that can have contrasting consequences for community stability. We therefore highlight the need to quantify sensitive species empirically and to develop response diversity metrics that distinguish these components. Author SummaryUnderstanding why some communities remain stable despite environmental change is a longstanding goal in ecology. Response diversity, which refers to differences in how species respond to environmental change, has been proposed as a key mechanism explaining why greater biodiversity (species richness) can promote ecological stability. Because species respond differently to changing environments, declines in some species can be compensated by increases in others, helping to stabilize community dynamics. However, previous studies have rarely considered species that are insensitive to current environmental changes. Using a mathematical model, we show that response diversity can arise from two distinct biological components--the number of sensitive species and variation in their responses--and that these components can have contrasting effects on ecological stability. When response diversity reflects variation among sensitive species, greater response diversity stabilizes community dynamics, as expected. In contrast, when response diversity changes only because of the proportions of sensitive and insensitive species, higher response diversity can be associated with lower community stability. Our findings highlight the importance of quantifying the number of sensitive species and developing response diversity metrics that distinguish species sensitivity from variation in responses among sensitive species.

After introduction, many non-native marine species are dispersed planktonically. Secondary spread within the non-native range has been shown to prevent the establishment of the introduced species if the advection of larvae prevents sufficient return of larvae to maintain the population in the face of competition with native species. However, those studies have largely neglected the effects of spatial variation in alongshore larval transport. We examine the introduction of a novel species with planktonic dispersal into a more realistic coastal environment which includes spatial variation in larval transport estimated from the Mercator Ocean 1/12th degree global circulation model. The introduction may either be from a distant habitat, or through range expansion. We find that there are locations in the global coastal ocean where introduced species are more likely to persist because of spatial variation of coastal currents. These include regions where alongshore larval transport diverges, such as estuaries. The location where a non-native species is introduced may not be where it flourishes - it cannot be assumed that the region where invading species are first noticed to be abundant is the region where it was introduced. We extend closed-population theory to open coastal systems to estimate persistence as a function of local circulation, habitat extent, and the competitive advantage of the introduced species. Software is provided which allows the estimations of regions where introduced species are more likely to persist and flourish as a function of larval depth behavior, planktonic duration and release timing.

O_LISeasonal migration in birds involves a substantial spatial redistribution of avian biodiversity each year and drives seasonal changes in community composition. Migrants experience different combinations of species interactions over space and time, generating regular disassembly and reassembly of bird communities throughout their annual cycles. However, the effects of seasonal migration on phylogenetic community structure remain poorly understood. C_LIO_LIWe assess spatiotemporal variation in phylogenetic community structure of North American passerines to test how seasonal migration restructures the evolutionary relatedness and dominant assembly mechanisms in bird communities throughout the annual cycle. Using distributional projections, we calculated metrics describing the phylogenetic dispersion of passerine communities each week of the year. We then tested the relationship between seasonal turnover in community phylogenetic dispersion and seasonal variation in species richness and proportion of migratory species. C_LIO_LISeasonal migration, by changing spatial patterns of avian diversity, simultaneously drives a complex continental redistribution of phylogenetic community structure. We find evidence of taxonomic scale dependency to our results, wherein throughout North America, the seasonal influx of migrant passerines yields communities that are overall more phylogenetically clustered, yet also exhibit greater phylogenetic overdispersion at smaller taxonomic scales. C_LIO_LISeasonal shifts in phylogenetic dispersion, though complex, track changes in diversity, manifesting as fluctuations in phylogenetic dispersion between northern and southern regions as seasonal migrants move between these regions. Our findings reveal a dynamic continental landscape of phylogenetic community structure directed by the movements of seasonal migrants. C_LI

Wildlife diseases cause well-documented and often dramatic reductions in host survival. However, the impact of infectious diseases on host reproduction remains understudied, especially with respect to effects of prior and/or current pathogen exposure on reproductive development. Here we experimentally tested how prior and/or current infection with a common bacterial pathogen, Mycoplasma gallisepticum ( MG), alters reproductive development for female versus male house finches (Haemorhous mexicanus). Finches were inoculated with either MG or sterile media while in wintering condition and subsequently received one of these treatments while in breeding condition. In females, MG exposure had both immediate and carry-over effects on reproduction: controls had higher odds of laying eggs compared to females inoculated with MG in spring only, higher odds than females inoculated in both winter and spring, and higher odds than females given MG in the winter only. Conversely, breeding-condition males inoculated with MG in spring had higher testosterone levels than males receiving only control inoculations, and there were no carryover effects of winter MG inoculation or inoculations during both seasons on testosterone. Sex bias in the reproductive impacts of infectious diseases may have important knock-on effects on the epidemiology and population-regulating effects of pathogens, thereby warranting further study.

The tidal wetlands of Louisiana comprise about 25% of those found throughout the conterminous United States yet estimates of wetland loss rates in the region between 1932 and 2016 have exceeded 60 km2 yr-1. To mitigate further degradation and wetland loss in the region, a globally unprecedented $50B, 50-year plan for coastal Louisiana is driving restoration efforts, and demand exists from multiple stakeholders for regularly updated, regional-scale, accurate land cover information. We used machine learning (random forests; RF) and cloud computing to develop a new Landsat-based, marsh vegetation community geospatial dataset. The dataset depicts wetland vegetation community types defined in a previous study at annual (1985-2025) time steps at 30-m resolution. An RF algorithm was used to integrate training samples with feature variables derived from Landsat imagery, and the resulting geospatial data product achieved an overall correct classification rate of 78%. The approach for development of the land cover dataset presented here has potential for application in other coastal wetland habitats throughout the world.

The red fox, Vulpes vulpes, is abundant in England and can exert limiting effects on their avian and mammalian prey. Large-scale gamebird releases in the UK may be sustaining high predator numbers, leading to greater predation on other prey species, especially once gamebird stocks are depleted. However, little is known how gamebird release affects the broader diet of the red fox. Here we used DNA metabarcoding to assess the diet of foxes from 18 agricultural estates in Southern England, 10 of which released large numbers of gamebirds (red-legged partridges - Alectoris rufus, and common pheasants - Phasianus colchicus) and 8 which did not. Scats were collected over one year, allowing for seasonal investigation of the foxes' diet. We investigated the vertebrate species consumed and compared the non-gamebird dietary diversity and composition between release and non-release estates and across seasons. The field vole (Microtus agrestis) was found to be the most frequently predated species overall. Brown hares and field voles were detected more on release sites, while bank voles and dog faeces were detected more on non-release sites. We found little evidence that foxes predate ground nesting birds or other species of concern. The dietary diversity was significantly lower on estates that released gamebirds, and this difference was most notable during the post-shoot, spring months (February to April). On both estate types, diversity was highest in the summer months. The altered predatory behaviour due to gamebird release is likely to affect the populations of the non-game prey of the red fox and this should be considered when designing policies regarding gamebird management and biodiversity conservation. Alternative predation combined with predator control may be reducing predation pressure on non-game prey where gamebirds are released, but high fox density elsewhere likely results in higher predation pressure on a wide range of species.

As ecosystems face a growing number of threats, coextinctions (resultant extinctions following a primary extinction) are expected to proliferate. However, less is known about the conditions under which coextinctions could outpace primary extinctions. Because coextinctions often occur through lost species interactions, we posit that aspects of food web structure and complexity can help predict differences in vulnerability to coextinction across ecosystems. To test this, we leverage Bayesian network models to assess the extent to which variation in ecosystem vulnerability to coextinction varies with food web structure. We find that food webs with high maximum trophic level are most vulnerable to coextinction, and that maximum trophic level is a better predictor than other aspects of food web structure, such as species richness or trophic connectance. Extending this approach, we also find that maximum trophic level uncovers the relative vulnerability of ecosystem services to species coextinction across 12 empirical food webs.