Bayesian joint models with INLA exploring marine mobile predator–prey and competitor species habitat overlap

Understanding spatial physical habitat selection driven by competition and/or predator–prey interactions of mobile marine species is a fundamental goal of spatial ecology. However, spatial counts or density data for highly mobile animals often (1) include excess zeros, (2) have spatial correlation, and (3) have highly nonlinear relationships with physical habitat variables, which results in the need for complex joint spatial models. In this paper, we test the use of Bayesian hierarchical hurdle and zero‐inflated joint models with integrated nested Laplace approximation (INLA), to fit complex joint models to spatial patterns of eight mobile marine species (grey seal, harbor seal, harbor porpoise, common guillemot, black‐legged kittiwake, northern gannet, herring, and sandeels). For each joint model, we specified nonlinear smoothed effect of physical habitat covariates and selected either competing species or predator–prey interactions. Out of a range of six ecologically important physical and biologic variables that are predicted to change with climate change and large‐scale energy extraction, we identified the most important habitat variables for each species and present the relationships between these bio/physical variables and species distributions. In particular, we found that net primary production played a significant role in determining habitat preferences of all the selected mobile marine species. We have shown that the INLA method is well‐suited for modeling spatially correlated data with excessive zeros and is an efficient approach to fit complex joint spatial models with nonlinear effects of covariates. Our approach has demonstrated its ability to define joint habitat selection for both competing and prey–predator species that can be relevant to numerous issues in the management and conservation of mobile marine species.     

Predicting the distributions of predator (snow leopard) and prey (blue sheep) under climate change in the Himalaya

Future climate change is likely to affect distributions of species, disrupt biotic interactions, and cause spatial incongruity of predator–prey habitats. Understanding the impacts of future climate change on species distribution will help in the formulation of conservation policies to reduce the risks of future biodiversity losses. Using a species distribution modeling approach by MaxEnt, we modeled current and future distributions of snow leopard (Panthera uncia) and its common prey, blue sheep (Pseudois nayaur), and observed the changes in niche overlap in the Nepal Himalaya. Annual mean temperature is the major climatic factor responsible for the snow leopard and blue sheep distributions in the energy‐deficient environments of high altitudes. Currently, about 15.32% and 15.93% area of the Nepal Himalaya are suitable for snow leopard and blue sheep habitats, respectively. The bioclimatic models show that the current suitable habitats of both snow leopard and blue sheep will be reduced under future climate change. The predicted suitable habitat of the snow leopard is decreased when blue sheep habitats is incorporated in the model. Our climate‐only model shows that only 11.64% (17,190 km²) area of Nepal is suitable for the snow leopard under current climate and the suitable habitat reduces to 5,435 km² (reduced by 24.02%) after incorporating the predicted distribution of blue sheep. The predicted distribution of snow leopard reduces by 14.57% in 2030 and by 21.57% in 2050 when the predicted distribution of blue sheep is included as compared to 1.98% reduction in 2030 and 3.80% reduction in 2050 based on the climate‐only model. It is predicted that future climate may alter the predator–prey spatial interaction inducing a lower degree of overlap and a higher degree of mismatch between snow leopard and blue sheep niches. This suggests increased energetic costs of finding preferred prey for snow leopards – a species already facing energetic constraints due to the limited dietary resources in its alpine habitat. Our findings provide valuable information for extension of protected areas in future.     

Diversity in thermal affinity among key piscivores buffers impacts of ocean warming on predator–prey interactions

Asymmetries in responses to climate change have the potential to alter important predator–prey interactions, in part by altering the location and size of spatial refugia for prey. We evaluated the effect of ocean warming on interactions between four important piscivores and four of their prey in the U.S. Northeast Shelf by examining species overlap under historical conditions (1968–2014) and with a doubling in CO₂. Because both predator and prey shift their distributions in response to changing ocean conditions, the net impact of warming or cooling on predator–prey interactions was not determined a priori from the range extent of either predator or prey alone. For Atlantic cod, an historically dominant piscivore in the region, we found that both historical and future warming led to a decline in the proportion of prey species’ range it occupied and caused a potential reduction in its ability to exert top‐down control on these prey. In contrast, the potential for overlap of spiny dogfish with prey species was enhanced by warming, expanding their importance as predators in this system. In sum, the decline in the ecological role for cod that began with overfishing in this ecosystem will likely be exacerbated by warming, but this loss may be counteracted by the rise in dominance of other piscivores with contrasting thermal preferences. Functional diversity in thermal affinity within the piscivore guild may therefore buffer against the impact of warming on marine ecosystems, suggesting a novel mechanism by which diversity confers resilience.     

A spatial theory for emergent multiple predator–prey interactions in food webs

Predator–prey interaction is inherently spatial because animals move through landscapes to search for and consume food resources and to avoid being consumed by other species. The spatial nature of species interactions necessitates integrating spatial processes into food web theory and evaluating how predators combine to impact their prey. Here, we present a spatial modeling approach that examines emergent multiple predator effects on prey within landscapes. The modeling is inspired by the habitat domain concept derived from empirical synthesis of spatial movement and interactions studies. Because these principles are motivated by synthesis of short‐term experiments, it remains uncertain whether spatial contingency principles hold in dynamical systems. We address this uncertainty by formulating dynamical systems models, guided by core habitat domain principles, to examine long‐term multiple predator–prey spatial dynamics. To describe habitat domains, we use classical niche concepts describing resource utilization distributions, and assume species interactions emerge from the degree of overlap between species. The analytical results generally align with those from empirical synthesis and present a theoretical framework capable of demonstrating multiple predator effects that does not depend on the small spatial or temporal scales typical of mesocosm experiments, and help bridge between empirical experiments and long‐term dynamics in natural systems.     

Climate and Demography Dictate the Strength of Predator-Prey Overlap in a Subarctic Marine Ecosystem

There is growing evidence that climate and anthropogenic influences on marine ecosystems are largely manifested by changes in species spatial dynamics. However, less is known about how shifts in species distributions might alter predator-prey overlap and the dynamics of prey populations. We developed a general approach to quantify species spatial overlap and identify the biotic and abiotic variables that dictate the strength of overlap. We used this method to test the hypothesis that population abundance and temperature have a synergistic effect on the spatial overlap of arrowtooth flounder (predator) and juvenile Alaska walleye pollock (prey, age-1) in the eastern Bering Sea. Our analyses indicate that (1) flounder abundance and temperature are key variables dictating the strength of flounder and pollock overlap, (2) changes in the magnitude of overlap may be largely driven by density-dependent habitat selection of flounder, and (3) species overlap is negatively correlated to juvenile pollock recruitment when flounder biomass is high. Overall, our findings suggest that continued increases in flounder abundance coupled with the predicted long-term warming of ocean temperatures could have important implications for the predator-prey dynamics of arrowtooth flounder and juvenile pollock. The approach used in this study is valuable for identifying potential consequences of climate variability and exploitation on species spatial dynamics and interactions in many marine ecosystems.     

Integrating disparate datasets to model the functional response of a marine predator: A case study of harbour porpoises in the southern North Sea

1. Quantifying consumption and prey choice for marine predator species is key to understanding their interaction with prey species, fisheries, and the ecosystem as a whole. However, parameterizing a functional response for large predators can be challenging because of the difficulty in obtaining the required data on predator diet and on the availability of multiple prey species.
2. This study modeled a multi‐species functional response (MSFR) to describe the relationship between consumption by harbour porpoises (Phocoena phocoena) and the availability of multiple prey species in the southern North Sea. Bayesian methodology was employed to estimate MSFR parameters and to incorporate uncertainties in diet and prey availability estimates. Prey consumption was estimated from stomach content data from stranded harbour porpoises. Prey availability to harbour porpoises was estimated based on the spatial overlap between prey distributions, estimated from fish survey data, and porpoise foraging range in the days prior to stranding predicted from telemetry data.
3. Results indicated a preference for sandeels in the study area. Prey switching behavior (change in preference dependent on prey abundance) was confirmed by the favored type III functional response model. Variation in the size of the foraging range (estimated area where harbour porpoises could have foraged prior to stranding) did not alter the overall pattern of the results or conclusions.
4. Integrating datasets on prey consumption from strandings, predator foraging distribution using telemetry, and prey availability from fish surveys into the modeling approach provides a methodological framework that may be appropriate for fitting MSFRs for other predators.

     

Improved estimates of digestion correction factors and passage rates for harbor seal (Phoca vitulina) prey in the northeast Atlantic

Diet composition in pinnipeds is widely estimated using hard prey remains recovered from feces. To estimate the size and number of prey represented in fecal samples accurately, digestion correction factors (DCFs) must be applied to measurements and counts of fish otoliths and cephalopod beaks. In this study, 101 whole prey feeding trials were conducted with six harbor seals (Phoca vitulina) and 18 prey species. We derived species‐ and grade‐specific estimates of digestion coefficients (DCs) and species‐specific recovery rates (RRs) to account for partial and complete digestion, respectively. Greater than 98% of otoliths were passed within three days of consumption. RRs were smallest for Atlantic salmon smolts (RR = 0.306, SE = 0.031) and increasingly larger for sandeels (RR = 0.494, SE = 0.017), flatfish (RR = 0.789, SE = 0.033), and large gadoids (RR = 0.944, SE = 0.034). Species‐specific otolith width DCs were smallest for Trisopterus species (DC = 1.14, SE = 0.015) and increasingly larger for flatfish (DC = 1.27, SE = 0.045), large gadoids (DC = 1.32, SE = 0.067) and sandeels (DC = 1.57, SE = 0.035). RRs were similar to those from gray seals (Halichoerus grypus), but harbor seal species‐ and grade‐specific DCs were generally smaller. Differences in partial and complete digestion rates among prey species and between seal species highlight the importance of applying DCFs when reconstructing diet.     

Geometric factors influencing the diet of vertebrate predators in marine and terrestrial environments

Predator–prey relationships are vital to ecosystem function and there is a need for greater predictive understanding of these interactions. We develop a geometric foraging model predicting minimum prey size scaling in marine and terrestrial vertebrate predators taking into account habitat dimensionality and biological traits. Our model predicts positive predator–prey size relationships on land but negative relationships in the sea. To test the model, we compiled data on diets of 794 predators (mammals, snakes, sharks and rays). Consistent with predictions, both terrestrial endotherm and ectotherm predators have significantly positive predator–prey size relationships. Marine predators, however, exhibit greater variation. Some of the largest predators specialise on small invertebrates while others are large vertebrate specialists. Prey–predator mass ratios were generally higher for ectothermic than endothermic predators, although dietary patterns were similar. Model‐based simulations of predator–prey relationships were consistent with observed relationships, suggesting that our approach provides insights into both trends and diversity in predator–prey interactions.     

Climate change effects on behavioral and physiological ecology of predator–prey interactions: Implications for conservation biological control

Habitat management under the auspices of conservation biological control is a widely used approach to foster conditions that ensure a diversity of predator species can persist spatially and temporally within agricultural landscapes in order to control their prey (pest) species. However, an emerging new factor, global climate change, has the potential to disrupt existing conservation biological control programs. Climate change may alter abiotic conditions such as temperature, precipitation, humidity and wind that in turn could alter the life-cycle timing of predator and prey species and the behavioral nature and strength of their interactions. Anticipating how climate change will affect predator and prey communities represents an important research challenge. We present a conceptual framework—the habitat domain concept—that is useful for understanding contingencies in the nature of predator diversity effects on prey based on predator and prey spatial movement in their habitat. We illustrate how this framework can be used to forecast whether biological control by predators will become more effective or become disrupted due to changing climate. We discuss how changes in predator–prey interactions are contingent on the tolerances of predators and prey species to changing abiotic conditions as determined by the degree of local adaptation and phenotypic plasticity exhibited by species populations. We conclude by discussing research approaches that are needed to help adjust conservation biological control management to deal with a climate future.     

How population loss through habitat boundaries determines the dynamics of a predator–prey system

The increased persistence of predator–prey systems when interactions are distributed through the space has been acknowledged by both empirical and theoretical studies. One salient feature of predator–prey interactions in heterogeneous space, for example, is the existence of cycles with reduced amplitude when compared with a homogeneous landscape. Although the role of spatial interactions in shaping the dynamics of predator–prey systems has been extensively studied, still very few works have focused on the effects of habitat loss and fragmentation on these systems. In this work, we study the population dynamics of a predator–prey system in a single finite habitat with flux at the boundaries. Species movement and growth are described through a reaction–diffusion model with Rosenzweig–MacArthur type local interactions. Conforming with the existing literature, we find that the reduction of habitat size, or increasing of species movement rates equivalently, has the potential to decrease the amplitude of oscillations and even bring the system to a steady coexistence equilibrium above a threshold. We observe, however, situations in which this trend is reversed. This occurs when species movement rates and response at patch boundaries interact to induce non-trivial patterns of species distributions. These distributions are characterized by anti-correlation between predator and prey, creating then spatial refugia for prey. Our results highlight the role of population loss through habitat boundaries in determining the dynamics of predator–prey interactions.     

Inside‐container effects drive mosquito community structure in Brazilian Atlantic forest

Metacommunity theory is a convenient framework in which to investigate how local communities linked by dispersal influence patterns of species distribution and abundance across large spatial scales. For organisms with complex life cycles, such as mosquitoes, different pressures are expected to act on communities due to behavioral and ecological partitioning of life stages. Adult females select habitats for oviposition, and resulting offspring are confined to that habitat until reaching adult stages capable of flight; outside‐container effects (OCE) (i.e., spatial factors) are thus expected to act more strongly on species distributions as a function of adult dispersal capability, which should be limited by geographic distances between sites. However, larval community dynamics within a habitat are influenced by inside‐container effects (ICE), mainly interactions with conspecifics and heterospecifics (e.g., through effects of competition and predation). We used a field experiment in a mainland‐island scenario to assess whether environmental, spatial, and temporal factors influence mosquito prey and predator distributions and abundances across spatial scales: within‐site, between‐site, and mainland‐island. We also evaluated whether predator abundances inside containers play a stronger role in shaping mosquito prey community structure than do OCE (e.g., spatial and environmental factors). Temporal influence was more important for predators than for prey mosquito community structure, and the changes in prey mosquito species composition over time appear to be driven by changes in predator abundances. There was a negligible effect of spatial and environmental factors on mosquito community structure, and temporal effects on mosquito abundances and distributions appear to be driven by changes in abundance of the dominant predator, perhaps because ICE are stronger than OCE due to larval habitat restriction, or because adult dispersal is not limited at the chosen spatial scales.     

Modeling plant species distributions under future climates: how fine scale do climate projections need to be?

Recent studies suggest that species distribution models (SDMs) based on fine‐scale climate data may provide markedly different estimates of climate‐change impacts than coarse‐scale models. However, these studies disagree in their conclusions of how scale influences projected species distributions. In rugged terrain, coarse‐scale climate grids may not capture topographically controlled climate variation at the scale that constitutes microhabitat or refugia for some species. Although finer scale data are therefore considered to better reflect climatic conditions experienced by species, there have been few formal analyses of how modeled distributions differ with scale. We modeled distributions for 52 plant species endemic to the California Floristic Province of different life forms and range sizes under recent and future climate across a 2000‐fold range of spatial scales (0.008–16 km²). We produced unique current and future climate datasets by separately downscaling 4 km climate models to three finer resolutions based on 800, 270, and 90 m digital elevation models and deriving bioclimatic predictors from them. As climate‐data resolution became coarser, SDMs predicted larger habitat area with diminishing spatial congruence between fine‐ and coarse‐scale predictions. These trends were most pronounced at the coarsest resolutions and depended on climate scenario and species' range size. On average, SDMs projected onto 4 km climate data predicted 42% more stable habitat (the amount of spatial overlap between predicted current and future climatically suitable habitat) compared with 800 m data. We found only modest agreement between areas predicted to be stable by 90 m models generalized to 4 km grids compared with areas classified as stable based on 4 km models, suggesting that some climate refugia captured at finer scales may be missed using coarser scale data. These differences in projected locations of habitat change may have more serious implications than net habitat area when predictive maps form the basis of conservation decision making.     

Marine reserves indirectly affect fine‐scale habitat associations,but not overall densities,of small benthic fishes

Many large, fishery‐targeted predatory species have attained very high relative densities as a direct result of protection by no‐take marine reserves. Indirect effects, via interactions with targeted species, may also occur for species that are not themselves targeted by fishing. In some temperate rocky reef ecosystems, indirect effects have caused profound changes in community structure, notably the restoration of predator–urchin–macroalgae trophic cascades. Yet, indirect effects on small benthic reef fishes remain poorly understood, perhaps because of behavioral associations with complex, refuge‐providing habitats. Few, if any, studies have evaluated any potential effects of marine reserves on habitat associations in small benthic fishes. We surveyed densities of small benthic fishes, including some endemic species of triplefin (Tripterygiidae), along with fine‐scale habitat features in kelp forests on rocky reefs in and around multiple marine reserves in northern New Zealand over 3 years. Bayesian generalized linear mixed models were used to evaluate evidence for (1) main effects of marine reserve protection, (2) associations with habitat gradients, including complexity, and (3) differences in habitat associations inside versus outside reserves. No evidence of overall main effects of marine reserves on species richness or densities of fishes was found. Both richness and densities showed strong associations with gradients in habitat features, particularly habitat complexity. In addition, some species exhibited reserve‐by‐habitat interactions, having different associations with habitat gradients inside versus outside marine reserves. Two species (Ruanoho whero and Forsterygion flavonigrum) showed stronger positive associations with habitat complexity inside reserves. These results are consistent with the presence of a behavioral risk effect, whereby prey fishes are more strongly attracted to habitats that provide refuge from predation in areas where predators are more abundant. This work highlights the importance of habitat structure and the potential for fishing to affect behavioral interactions and the interspecific dynamic attributes of community structure beyond simple predator–prey consumption and archetypal trophic cascades.     

A cascade of warming impacts brings bluefin tuna to Greenland waters

Rising ocean temperatures are causing marine fish species to shift spatial distributions and ranges, and are altering predator‐prey dynamics in food webs. Most documented cases of species shifts so far involve relatively small species at lower trophic levels, and consider individual species in ecological isolation from others. Here, we show that a large highly migratory top predator fish species has entered a high latitude subpolar area beyond its usual range. Bluefin tuna, Thunnus thynnus Linnaeus 1758, were captured in waters east of Greenland (65°N) in August 2012 during exploratory fishing for Atlantic mackerel, Scomber scombrus Linnaeus 1758. The bluefin tuna were captured in a single net‐haul in 9–11 °C water together with 6 tonnes of mackerel, which is a preferred prey species and itself a new immigrant to the area. Regional temperatures in August 2012 were historically high and contributed to a warming trend since 1985, when temperatures began to rise. The presence of bluefin tuna in this region is likely due to a combination of warm temperatures that are physiologically more tolerable and immigration of an important prey species to the region. We conclude that a cascade of climate change impacts is restructuring the food web in east Greenland waters.     

Increased autumn rainfall disrupts predator–prey interactions in fragmented boreal forests

There is a pressing need to understand how changing climate interacts with land‐use change to affect predator–prey interactions in fragmented landscapes. This is particularly true in boreal ecosystems facing fast climate change and intensification in forestry practices. Here, we investigated the relative influence of autumn climate and habitat quality on the food‐storing behaviour of a generalist predator, the pygmy owl, using a unique data set of 15 850 prey items recorded in western Finland over 12 years. Our results highlighted strong effects of autumn climate (number of days with rainfall and with temperature <0 °C) on food‐store composition. Increasing frequency of days with precipitation in autumn triggered a decrease in (i) total prey biomass stored, (ii) the number of bank voles (main prey) stored, and (iii) the scaled mass index of pygmy owls. Increasing proportions of old spruce forests strengthened the functional response of owls to variations in vole abundance and were more prone to switch from main prey to alternative prey (passerine birds) depending on local climate conditions. High‐quality habitat may allow pygmy owls to buffer negative effects of inclement weather and cyclic variation in vole abundance. Additionally, our results evidenced sex‐specific trends in body condition, as the scaled mass index of smaller males increased while the scaled mass index of larger females decreased over the study period, probably due to sex‐specific foraging strategies and energy requirements. Long‐term temporal stability in local vole abundance refutes the hypothesis of climate‐driven change in vole abundance and suggests that rainier autumns could reduce the vulnerability of small mammals to predation by pygmy owls. As small rodents are key prey species for many predators in northern ecosystems, our findings raise concern about the impact of global change on boreal food webs through changes in main prey vulnerability.     

Species interactions and climate change: How the disruption of species co‐occurrence will impact on an avian forest guild

Interspecific interactions are crucial in determining species occurrence and community assembly. Understanding these interactions is thus essential for correctly predicting species' responses to climate change. We focussed on an avian forest guild of four hole‐nesting species with differing sensitivities to climate that show a range of well‐understood reciprocal interactions, including facilitation, competition and predation. We modelled the potential distributions of black woodpecker and boreal, tawny and Ural owl, and tested whether the spatial patterns of the more widespread species (excluding Ural owl) were shaped by interspecific interactions. We then modelled the potential future distributions of all four species, evaluating how the predicted changes will alter the overlap between the species' ranges, and hence the spatial outcomes of interactions. Forest cover/type and climate were important determinants of habitat suitability for all species. Field data analysed with N‐mixture models revealed effects of interspecific interactions on current species abundance, especially in boreal owl (positive effects of black woodpecker, negative effects of tawny owl). Climate change will impact the assemblage both at species and guild levels, as the potential area of range overlap, relevant for species interactions, will change in both proportion and extent in the future. Boreal owl, the most climate‐sensitive species in the guild, will retreat, and the range overlap with its main predator, tawny owl, will increase in the remaining suitable area: climate change will thus impact on boreal owl both directly and indirectly. Climate change will cause the geographical alteration or disruption of species interaction networks, with different consequences for the species belonging to the guild and a likely spatial increase of competition and/or intraguild predation. Our work shows significant interactions and important potential changes in the overlap of areas suitable for the interacting species, which reinforce the importance of including relevant biotic interactions in predictive climate change models for increasing forecast accuracy.     

From projected species distribution to food‐web structure under climate change

Climate change is inducing deep modifications in species geographic ranges worldwide. However, the consequences of such changes on community structure are still poorly understood, particularly the impacts on food‐web properties. Here, we propose a new framework, coupling species distribution and trophic models, to predict climate change impacts on food‐web structure across the Mediterranean Sea. Sea surface temperature was used to determine the fish climate niches and their future distributions. Body size was used to infer trophic interactions between fish species. Our projections reveal that 54 fish species of 256 endemic and native species included in our analysis would disappear by 2080–2099 from the Mediterranean continental shelf. The number of feeding links between fish species would decrease on 73.4% of the continental shelf. However, the connectance of the overall fish web would increase on average, from 0.26 to 0.29, mainly due to a differential loss rate of feeding links and species richness. This result masks a systematic decrease in predator generality, estimated here as the number of prey species, from 30.0 to 25.4. Therefore, our study highlights large‐scale impacts of climate change on marine food‐web structure with potential deep consequences on ecosystem functioning. However, these impacts will likely be highly heterogeneous in space, challenging our current understanding of climate change impact on local marine ecosystems.     

Density‐dependent and ‐independent effects on the joint use of space by predators and prey in terrestrial arthropod food‐webs

The spatial distribution of predators and their prey is affected by their joint use of space. While the formation of such spatial patterns may be driven by density‐dependent and ‐independent factors our knowledge on the contribution of different land‐use activities on the formation of spatial patterns between predators and prey remains very limited. Agriculture is one of the most prevailing land‐use activities with strong effects on invertebrate densities and structural habitat conditions. Here, we used replicated conventionally and organically managed winter wheat fields to investigate the effects of agricultural land‐use on the spatial patterns of generalist predators and decomposer prey. We then identified the explanatory power of density‐dependent (prey and predator activity density) and density‐independent (vegetation structure) predictors for the observed spatial patterns. Generalist predators were regularly distributed only in conventionally managed fields and this pattern intensified with decreasing Collembola prey availability and increasing spider activity density. Segregation between carabid and spider predators was strongest in fields with lowest wheat plant height, suggesting more intense intraguild interactions in structurally less complex habitats. Collembola were aggregated independent of management and aggregation was strongest in fields with highest Collembola and carabid activity density. Spiders and Collembola prey were associated, but higher aphid densities under conventional management weakened or interrupted this spatial relationship. We conclude that active control of crop plant physiognomy by growth hormones and herbicides in conventionally managed fields promotes predator–predator segregation and that a high availability of aphid prey seems to decouple predator–Collembola prey associations. Our results emphasise the need for a more mechanistic understanding of the effects of land‐use on the formation of spatial patterns and species interactions, especially under scenarios of environmental change and an ongoing loss of biodiversity.     

Habitat overlap between predatory benthic fish and their invertebrate prey in streams: the relative influence of spatial and temporal factors on predation risk

1. The spatial heterogeneity of ecosystems as well as temporal activity patterns of organisms can have far‐reaching effects on predator–prey relationships. We hypothesised that spatiotemporal constraints in mesohabitat use by benthic fish predators would reduce habitat overlap with benthic invertebrates and lead to mesohabitat‐specific predation risks. 2. We analysed the spatiotemporal activity patterns of two small‐bodied benthivorous fishes, gudgeon (Gobio gobio) and stone loach (Barbatula barbatula), and of benthic invertebrates in a small temperate stream during three 24‐h field experiments. By applying a novel method of field video observation, we monitored the spatiotemporal foraging behaviour of the fish in their natural environment. A parallel analysis of invertebrate mesohabitat use by means of small area Hess sampling allowed a direct estimation of habitat overlap at a pool–riffle scale. 3. Gudgeon showed a dominant spatial activity pattern preferring pools at all times of day, whereas stone loach used both mesohabitats but with a distinct temporal (nocturnal) activity pattern. The patterns of residence were not identical with those of active foraging. Invertebrate community composition differed significantly between mesohabitats but not between times of day. More than half of the total dissimilarity between pools and riffles was accounted for by six invertebrate taxa. Five of these were subject to higher fish predation in pools than in riffles. The total prey consumption of the two fish species together in pools was about three times as high as in riffles. Trophic niche breadth of stone loach and thus its predation range was broader than that of gudgeon. 4. These results indicate that the potential predation risk for stream invertebrates depends on the combination of spatial and temporal patterns of both predator and prey. Given the distinct differences in predation risk found between pools and riffles, we conclude that spatial heterogeneity at the mesohabitat scale can influence mechanisms and consequences of selective predation. We also suggest that the analysis of spatiotemporal predator–prey relationships should not be based on the premise that the main residence habitat and active foraging habitat of a predator are identical.     

Surreptitious sympatry: Exploring the ecological and genetic separation of two sibling species

Climate change is having profound impacts on animal populations, and shifts in geographic range are predicted in response. Shifts that result in range overlap between previously allopatric congeneric species may have consequences for biodiversity through interspecific competition, hybridization, and genetic introgression. Harbor seals (Phoca vitulina) and spotted seals (Phoca largha) are parapatric sibling species and areas of co‐occurrence at the edges of their range, such as Bristol Bay, Alaska, offer a unique opportunity to explore ecological separation and discuss potential consequences of increased range overlap resulting from retreating sea ice. Using telemetry and genetic data from 14 harbor seals and six spotted seals, we explored the ecological and genetic separation of the two species by comparing their utilization distributions, distance from haul‐out, dive behavior (e.g., depth, duration, focus), and evidence of hybridization. Firstly, we show that harbor and spotted seals, which cannot be visually distinguished definitively in all cases, haul‐out together side by side in Bristol Bay from late summer to early winter. Secondly, we observed subtle rather than pronounced differences in ranging patterns and dive behavior during this period. Thirdly, most spotted seals in this study remained close to shore in contrast to what is known of the species in more northern areas, and lastly, we did not find any evidence of hybridization. The lack of distinct ecological separation in this area of sympatry suggests that interspecific competition could play an important role in the persistence of these species, particularly if range overlap will increase as a result of climate‐induced range shifts and loss of spotted seal pagophilic breeding habitat. Our results also highlight the added complexities in monitoring these species in areas of suspected overlap, as they cannot easily be distinguished without genetic analysis. Predicted climate‐induced environmental change will likely influence the spatial and temporal extent of overlap in these two sibling species. Ultimately, this may alter the balance between current isolating mechanisms with consequences for species integrity and fitness.