Empirical evidence for source-sink populations: a review on occurrence, assessments and implications

A new threshold selection method for species distribution models with presence-only data: Extracting the mutation point of the P/E curve by threshold regression

Selecting thresholds to convert continuous predictions of species distribution models proves critical for many real-world applications and model assessments. Prevalent threshold selection methods for presence-only data require unproven pseudo-absence data or subjective researchers' decisions. This study proposes a new method, Boyce-Threshold Quantile Regression (BTQR), to determine thresholds objectively without pseudo-absence data. We summarize that the mutation point is a typical shape feature of the predicted-to-expected (P/E) curve after reviewing relevant articles. Analysis based on source-sink theory suggests that this mutation point may represent a transition in habitat types and serve as an appropriate threshold. Threshold regression is introduced to accurately locate the mutation point. To validate the effectiveness of BTQR, we used four virtual species of varying prevalence and a real species with reliable distribution data. Six different species distribution models were employed to generate continuous suitability predictions. BTQR and nine other traditional methods transformed these continuous outputs into binary results. Comparative experiments show that BTQR has advantages in terms of accuracy, applicability, and consistency over the existing methods.

Higher Mortality Is Not a Universal Cost of Dispersal: A Case Study in African Wild Dogs

AbstractMortality is considered one of the main costs of dispersal. A reliable evaluation of mortality, however, is often hindered by a lack of information about the fate of individuals that disappear under unexplained circumstances (i.e., missing individuals). Here, we addressed this uncertainty by applying a Bayesian mortality analysis that inferred the fate of missing individuals according to information from individuals with known fate. Specifically, we tested the hypothesis that mortality during dispersal is higher than mortality among nondispersers using 32 years of mark-resighting data from a free-ranging population of the endangered African wild dog (Lycaon pictus) in northern Botswana. Contrary to expectations, we found that mortality during dispersal was lower than mortality among nondispersers, indicating that higher mortality is not a universal cost of dispersal. Our findings suggest that group living can incur costs for certain age classes, such as limited access to resources as group density increases, that exceed the mortality costs associated with dispersal. By challenging the accepted expectation of higher mortality during dispersal, we urge for further investigations of this key life history trait and propose a robust statistical approach to reduce bias in mortality estimates.

Chemical Mixtures and Multiple Stressors: Same but Different?

Ecosystems are strongly influenced by multiple anthropogenic stressors, including a wide range of chemicals and their mixtures. Studies on the effects of multiple stressors have largely focussed on nonchemical stressors, whereas studies on chemical mixtures have largely ignored other stressors. However, both research areas face similar challenges and require similar tools and methods to predict the joint effects of chemicals or nonchemical stressors, and frameworks to integrate multiple chemical and nonchemical stressors are missing. We provide an overview of the research paradigms, tools, and methods commonly used in multiple stressor and chemical mixture research and discuss potential domains of cross-fertilization and joint challenges. First, we compare the general paradigms of ecotoxicology and (applied) ecology to explain the historical divide. Subsequently, we compare methods and approaches for the identification of interactions, stressor characterization, and designing experiments. We suggest that both multiple stressor and chemical mixture research are too focused on interactions and would benefit from integration regarding null model selection. Stressor characterization is typically more costly for chemical mixtures. While for chemical mixtures comprehensive classification systems at suborganismal level have been developed, recent classification systems for multiple stressors account for environmental context. Both research areas suffer from rather simplified experimental designs that focus on only a limited number of stressors, chemicals, and treatments. We discuss concepts that can guide more realistic designs capturing spatiotemporal stressor dynamics. We suggest that process-based and data-driven models are particularly promising to tackle the challenge of prediction of effects of chemical mixtures and nonchemical stressors on (meta-)communities and (meta-)food webs. We propose a framework to integrate the assessment of effects for multiple stressors and chemical mixtures. Environ Toxicol Chem 2023;42:1915-1936. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Coevolution of Patch Selection in Stochastic Environments

AbstractSpecies interact in landscapes where environmental conditions vary in time and space. This variability impacts how species select habitat patches. Under equilibrium conditions, evolution of this patch selection can result in ideal free distributions where per capita growth rates are zero in occupied patches and negative in unoccupied patches. These ideal free distributions, however, do not explain why species occupy sink patches, why competitors have overlapping spatial ranges, or why predators avoid highly productive patches. To understand these patterns, we solve for coevolutionarily stable strategies (coESSs) of patch selection for multispecies stochastic Lotka-Volterra models accounting for spatial and temporal heterogeneity. In occupied patches at the coESS, we show that the differences between the local contributions to the mean and the variance of the long-term population growth rate are equalized. Applying this characterization to models of antagonistic interactions reveals that environmental stochasticity can partially exorcize the ghost of competition past, select for new forms of enemy-free and victimless space, and generate hydra effects over evolutionary timescales. Viewing our results through the economic lens of modern portfolio theory highlights why the coESS for patch selection is often a bet-hedging strategy coupling stochastic sink populations. Our results highlight how environmental stochasticity can reverse or amplify evolutionary outcomes as a result of species interactions or spatial heterogeneity.

Catch Effectiveness Revealed by Site-Related Differences in Capture-Mark-Recapture Methods: A Butterfly Metapopulation Study

Understanding metapopulation structures is very important in the context of ecological studies and conservation. Crucial in this respect are the abundances of both the whole metapopulation and its constituent subpopulations. In recent decades, capture-mark-recapture studies have been considered the most reliable means of calculating such abundances. In butterfly studies, individual insects are usually caught with an entomological net. But the effectiveness of this method can vary for a number of reasons: differences between fieldworkers, in time, between sites etc. This article analyses catch effectiveness data with respect to two subpopulations of the Apollo butterfly (Parnassius apollo) metapopulation in the Pieniny National Park (Polish Carpathians). The results show that this parameter varied significantly between sites, probably because of differences in microrelief and plant cover. In addition, a method is proposed that will include information on catch effectiveness for estimating the sizes of particular subpopulations and will help to elucidate the structure of the entire metapopulation.

Population dynamics of western gorillas at Mbeli Bai

Long-term studies of population dynamics can provide insights into life history theory, population ecology, socioecology, conservation biology and wildlife management. Here we examine 25 years of population dynamics of western gorillas at Mbeli Bai, a swampy forest clearing in Nouabalé-Ndoki National Park, the Republic of Congo. The Mbeli population more than doubled from 101 to 226 gorillas during the study. After adjusting for a net influx of gorillas into the study population, the increase represents an inherent growth rate of 0.7% per year, with 95% confidence limits between -0.7% and 2.6%. The influx of gorillas mainly involved immigration of individuals into existing study groups (social dispersal), but it also included the appearance of a few previously unknown groups (locational dispersal). The average group size did not change significantly during the study, which is consistent with the possibility that western gorillas face socioecological constraints on group size, even when the population is increasing. We found no significant evidence of density dependence on female reproductive success or male mating competition. The distribution of gorillas among age/sex categories also did not change significantly, which suggests that the population had a stable age structure. Our results provide evidence of population stability or growth for some western gorillas (albeit within a small area). The results highlight the value of law enforcement, long-term monitoring, and protected areas; but they do not diminish the importance of improving conservation for this critically endangered species.

Habitat amount or landscape configuration: Emerging HotSpot analysis reveals the importance of habitat amount for a grassland bird in South Dakota

Habitat loss and fragmentation are two important drivers of biodiversity decline. Understanding how species respond to landscape composition and configuration in dynamic landscapes is of great importance for informing the conservation and management of grassland species. With limited conservation resources, prescribed management targeted at the appropriate landscape process is necessary for the effective management of species. We used pheasants (Phasianus colchicus) across South Dakota, USA as a model species to identify environmental factors driving spatiotemporal variation in population productivity. Using an emerging Hotspot analysis, we analyzed annual count data from 105 fixed pheasant brood routes over a 24-year period to identify high (HotSpot) and low (ColdSpot) pheasant population productivity areas. We then applied classification and regression tree modeling to evaluate landscape attributes associated with pheasant productivity among spatial scales (500 m and 1000 m). We found that the amount of grassland at a local spatial scale was the primary factor influencing an area being a HotSpot. Our results also demonstrated non-significant or weak effects of fragmentation per se on pheasant populations. These findings are in accordance with the habitat amount hypothesis highlighting the importance of habitat amount in the landscape for maintaining and increasing the pheasant population. We, therefore, recommend that managers should focus on increasing the total habitat area in the landscape and restoring degraded habitats. Our method of identifying areas of high productivity across the landscape can be applied to other species with count data.

The Role of Source-Sink Dynamics in the Assessment of Risk to Nontarget Arthropods from the Use of Plant Protection Products

The concept of source-sink dynamics as a potentially important component of metapopulation dynamics was introduced in the 1980s. The objective of the present review was to review the considerable body of work that has been developed, to consider its theoretical implications as well as to understand how source-sink dynamics may manifest under field conditions in the specific case of nontarget arthropods in the agricultural environment. Our review concludes that metapopulation dynamics based on field observations are often far more complex than existing theoretical source-sink models would indicate, because they are dependent on numerous population processes and influencing factors. The difficulty in identifying and measuring these factors likely explains why empirical studies assessing source-sink dynamics are scarce. Furthermore, we highlight the importance of considering the spatial and temporal heterogeneity of agricultural landscapes when assessing the population dynamics of nontarget arthropods in the context of the risk from the use of plant protection products. A need is identified to further develop and thoroughly validate predictive population models, which can incorporate all factors relevant to a specific system. Once reliable predictive models for a number of representative nontarget arthropod species are available, they could provide a meaningful tool for refined risk evaluations (higher tier level risk assessment), addressing specific concerns identified at the initial evaluation stages (lower tier level risk assessment). Environ Toxicol Chem 2021;40:2667-2679. © 2021 ERM, FMC, Syngenta, Bayer AG, BASF SE, Corteva agriscience. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

High turn-over rates at the upper range limit and elevational source-sink dynamics in a widespread songbird

The formation of an upper distributional range limit for species breeding along mountain slopes is often based on environmental gradients resulting in changing demographic rates towards high elevations. However, we still lack an empirical understanding of how the interplay of demographic parameters forms the upper range limit in highly mobile species. Here, we study apparent survival and within-study area dispersal over a 700 m elevational gradient in barn swallows (Hirundo rustica) by using 15 years of capture-mark-recapture data. Annual apparent survival of adult breeding birds decreased while breeding dispersal probability of adult females, but not males increased towards the upper range limit. Individuals at high elevations dispersed to farms situated at elevations lower than would be expected by random dispersal. These results suggest higher turn-over rates of breeding individuals at high elevations, an elevational increase in immigration and thus, within-population source-sink dynamics between low and high elevations. The formation of the upper range limit therefore is based on preference for low-elevation breeding sites and immigration to high elevations. Thus, shifts of the upper range limit are not only affected by changes in the quality of high-elevation habitats but also by factors affecting the number of immigrants produced at low elevations.

A sink host allows a specialist herbivore to persist in a seasonal source

In seasonal environments, sinks that are more persistent than sources may serve as temporal stepping stones for specialists. However, this possibility has to our knowledge, not been demonstrated to date, as such environments are thought to select for generalists, and the role of sinks, both in the field and in the laboratory, is difficult to document. Here, we used laboratory experiments to show that herbivorous arthropods associated with seasonally absent main (source) habitats can endure on a suboptimal (sink) host for several generations, albeit with a negative growth rate. Additionally, they dispersed towards this host less often than towards the main host and accepted it less often than the main host. Finally, repeated experimental evolution attempts revealed no adaptation to the suboptimal host. Nevertheless, field observations showed that arthropods are found in suboptimal habitats when the main habitat is unavailable. Together, these results show that evolutionary rescue in the suboptimal habitat is not possible. Instead, the sink habitat functions as a temporal stepping stone, allowing for the persistence of a specialist when the source habitat is gone.

Source-sink dynamics assists the maintenance of a pollinating wasp

Dispersal that unites spatially subdivided populations into a metapopulation with source-sink dynamics is crucial for species persistence in fragmented landscapes. Understanding such dynamics for pollinators is particularly urgent owing to the ongoing global pollination crisis. Here, we investigated the population structure and source-sink dynamics of a pollinating wasp (Wiebesia sp. 3) of Ficus pumila in the Zhoushan Archipelago of China. We found significant asymmetry in the pairwise migrant numbers for 22 of 28 cases on the historical timescale, but only two on the contemporary timescale. Despite a small population size, the sole island not colonized by a superior competitor wasp (Wiebesia sp. 1) consistently behaved as a net exporter of migrants, supplying large sinks. Comparable levels of genetic diversity, with few private alleles and low genetic differentiation (total F_st : 0.03; pairwise F_st : 0.0005-0.0791), were revealed among all the islands. There was a significant isolation-by-distance pattern caused mainly by migration between the competition-free island and other islands, otherwise the pattern was negligible. The clustering analysis failed to detect multiple gene pools for the whole region. Thus, the sinks were most probably organized into a patchy population. Moreover, the estimates of effective population sizes were comparable between the two timescales. Thus the source-sink dynamics embedded within a well-connected population network may allow Wiebesia sp. 3 to persist at a competitive disadvantage. This study provides evidence that metapopulations in the real world may be complicated and changeable over time, highlighting the necessity to study such metapopulations in detail.

Dispersal in a house sparrow metapopulation: An integrative case study of genetic assignment calibrated with ecological data and pedigree information

Dispersal has a crucial role determining ecoevolutionary dynamics through both gene flow and population size regulation. However, to study dispersal and its consequences, one must distinguish immigrants from residents. Dispersers can be identified using telemetry, capture-mark-recapture (CMR) methods, or genetic assignment methods. All of these methods have disadvantages, such as high costs and substantial field efforts needed for telemetry and CMR surveys, and adequate genetic distance required in genetic assignment. In this study, we used genome-wide 200K Single Nucleotide Polymorphism data and two different genetic assignment approaches (GSI_SIM, Bayesian framework; BONE, network-based estimation) to identify the dispersers in a house sparrow (Passer domesticus) metapopulation sampled over 16 years. Our results showed higher assignment accuracy with BONE. Hence, we proceeded to diagnose potential sources of errors in the assignment results from the BONE method due to variation in levels of interpopulation genetic differentiation, intrapopulation genetic variation and sample size. We show that assignment accuracy is high even at low levels of genetic differentiation and that it increases with the proportion of a population that has been sampled. Finally, we highlight that dispersal studies integrating both ecological and genetic data provide robust assessments of the dispersal patterns in natural populations.

Empirical tests of habitat selection theory reveal that conspecific density and patch quality, but not habitat amount, drive long-distance immigration in a wild bird

Individuals that disperse long distances from their natal site must select breeding patches with no prior knowledge of patch suitability. Despite decades of theoretical studies examining which cues dispersing individuals should use to select breeding patches, few empirical studies have tested the predictions of these theories at spatial scales relevant to long‐distance dispersal in wild animal populations. Here, we use a novel assignment model based on multiple intrinsic markers to quantify natal dispersal distances of Wood Thrush (Hylocichla mustelina) breeding in forest fragments. We show that long‐distance natal dispersal in this species is more frequent than commonly assumed for songbirds and that habitat selection by these individuals is driven by density‐dependence and patch quality but not the amount of habitat surrounding breeding patches. These results represent an important contribution to understanding habitat selection by dispersing individuals, especially with regards to long‐distance dispersal.

Why we should care about movements: Using spatially explicit integrated population models to assess habitat source-sink dynamics

Assessing the source–sink status of populations and habitats is of major importance for understanding population dynamics and for the management of natural populations. Sources produce a net surplus of individuals (per capita contribution to the metapopulation > 1) and will be the main contributors for self‐sustaining populations, whereas sinks produce a deficit (contribution < 1). However, making these types of assessments is generally hindered by the problem of separating mortality from permanent emigration, especially when survival probabilities as well as moved distances are habitat‐specific.

To address this long‐standing issue, we propose a spatial multi‐event integrated population model (IPM) that incorporates habitat‐specific dispersal distances of individuals. Using information about local movements, this IPM adjusts survival estimates for emigration outside the study area.

Analysing 24 years of data on a farmland passerine (the northern wheatear Oenanthe oenanthe), we assessed habitat‐specific contributions, and hence the source–sink status and temporal variation of two key breeding habitats, while accounting for habitat‐ and sex‐specific local dispersal distances of juveniles and adults. We then examined the sensitivity of the source–sink analysis by comparing results with and without accounting for these local movements.

Estimates of first‐year survival, and consequently habitat‐specific contributions, were higher when local movement data were included. The consequences from including movement data were sex specific, with contribution shifting from sink to likely source in one habitat for males, and previously noted habitat differences for females disappearing.

Assessing the source–sink status of habitats is extremely challenging. We show that our spatial IPM accounting for local movements can reduce biases in estimates of the contribution by different habitats, and thus reduce the overestimation of the occurrence of sink habitats. This approach allows combining all available data on demographic rates and movements, which will allow better assessment of source–sink dynamics and better informed conservation interventions.

When to stay and when to leave? Proximate causes of dispersal in an endangered social carnivore

Reliable estimates of birth, death, emigration and immigration rates are fundamental to understanding and predicting the dynamics of wild populations and, consequently, inform appropriate management actions. However, when individuals disappear from a focal population, inference on their fate is often challenging. Here we used 30 years of individual-based mark-recapture data from a population of free-ranging African wild dogs Lycaon pictus in Botswana and a suite of individual, social and environmental predictors to investigate factors affecting the decision to emigrate from the pack. We subsequently used this information to assign an emigration probability to those individuals that were no longer sighted within their pack (i.e. missing individuals). Natal dispersal (i.e. emigration from the natal pack) showed seasonal patterns with female dispersal peaking prior to the mating season and male dispersal peaking at the beginning of the wet season. For both sexes, natal dispersal rate increased in the absence of unrelated individuals of the opposite sex in the pack. Male natal dispersal decreased with increasing number of pups in the pack and increased in larger packs. Female natal dispersal decreased with increasing number of pups in larger packs, but increased with increasing number of pups in smaller packs. Individuals of both sexes were less likely to exhibit secondary dispersal (i.e. emigration from a pack other than the natal pack) if they were dominant and if many pups were present in the pack. Our models predicted that 18% and 25% of missing females and males, respectively, had likely dispersed from the natal pack, rather than having died. A misclassification of this order of magnitude between dispersal and mortality can have far-reaching consequences in the evaluation and prediction of population dynamics and persistence, and potentially mislead conservation actions. Our study showed that the decision to disperse is context-dependent and that the effect of individual, social and environmental predictors differs between males and females and between natal and secondary dispersal related to different direct and indirect fitness consequences. Furthermore, we demonstrated how a thorough understanding of the proximate causes of dispersal can be used to assign a dispersal probability to missing individuals. Knowledge of causes of dispersal can then be used within an integrated framework to more reliably estimate mortality rates.

Low wolverine (Gulo gulo) density in a national park complex of the Canadian Rocky Mountains

Large carnivores are sensitive to human-caused extirpation due to large home ranges, low population densities, and low reproductive rates. Protected areas help maintain populations by acting as sources, but human-caused mortality, habitat displacement, and edge effects occurring at protected area boundaries may reduce that function. The national parks Banff, Yoho, and Kootenay in the Canadian Rocky Mountains are refugia for large carnivores, including wolverines (Gulo gulo (Linnaeus, 1758)). Despite growing conservation concern, empirical baseline population data for wolverines remain scarce throughout their range, including most of Canada. We hypothesized (i) that in these national parks, wolverine density matched values expected for high-quality habitat, and (ii) that edge effects decreased density towards park boundaries. We conducted systematic non-invasive genetic sampling surveys covering >7000 km2 (2011 and 2013). Using spatial capture–recapture models, we estimated mean (±SE) female (1.5 ± 0.3 and 1.4 ± 0.3 wolverine/1000 km2), male (1.8 ± 0.4 and 1.5 ± 0.3 wolverine/1000 km2), and combined (3.3 ± 0.5 and 3.0 ± 0.4 wolverine/1000 km2) densities for 2011 and 2013, respectively. These estimates were lower than predictions based on density extrapolation from nearby high-quality habitat, and density decreased towards park boundaries. To benefit the population, we recommend creating buffer zones around parks that protect female habitat and prohibit harvest.

Contrasting demographic trends and asymmetric migration rates in a spatially structured amphibian population

Natural populations often persist at the landscape scale as metapopulations, with breeding units (subpopulations) experiencing temporal extinction and recolonization events. Important parameters to forecast population viability in these systems include the ratio of the effective number of breeders (N_b ) to the total number of adults (N_a ) and migration rates among subpopulations. Here, we present the results of a 10-year integrative monitoring program of a metapopulation of the Iberian green frog (Pelophylax perezi) in central Spain. We characterized population dynamics at two main breeding ponds (Gravera and Laguna) using capture-mark-recapture data to estimate N_a in each breeding season, and multilocus genotypes to estimate the effective population size (N_e ), N_b , individual breeding success, and migration rates. Both ponds experienced population decline after a dry season, with Gravera subsequently recovering and Laguna suffering a bottleneck associated with genetic impoverishment. In this subpopulation, average allelic richness and private alleles decreased from 2010 (10.87 and 1.67, respectively) to 2018 (8.0 and 0.20). The N_b /N_a ratio in Laguna in 2018 was twice as high (0.95) than in Gravera (0.41) or in pre-bottleneck Laguna (0.50), suggesting plasticity or genetic compensation through increased individual breeding success. Migration rates were asymmetric between ponds, with a stronger contribution from Gravera to Laguna (29.9% vs. 16.2% in the opposite direction) that may result in a rescue effect. This study emphasizes the importance of integrative demographic approaches for the monitoring of natural populations based on a better understanding of their spatio-temporal dynamics, which provides valuable information for conservation actions.

The spatial scale of dispersal revealed by admixture tracts

Evaluating species dispersal across the landscape is essential to design appropriate management and conservation actions. However, technical difficulties often preclude direct measures of individual movement, while indirect genetic approaches rely on assumptions that sometimes limit their application. Here, we show that the temporal decay of admixture tracts lengths can be used to assess genetic connectivity within a population introgressed by foreign haplotypes. We present a proof-of-concept approach based on local ancestry inference in a high gene flow marine fish species, the European sea bass (Dicentrarchus labrax). Genetic admixture in the contact zone between Atlantic and Mediterranean sea bass lineages allows the introgression of Atlantic haplotype tracts within the Mediterranean Sea. Once introgressed, blocks of foreign ancestry are progressively eroded by recombination as they diffuse from the western to the eastern Mediterranean basin, providing a means to estimate dispersal. By comparing the length distributions of Atlantic tracts between two Mediterranean populations located at different distances from the contact zone, we estimated the average per-generation dispersal distance within the Mediterranean lineage to less than 50 km. Using simulations, we showed that this approach is robust to a range of demographic histories and sample sizes. Our results thus support that the length of admixture tracts can be used together with a recombination clock to estimate genetic connectivity in species for which the neutral migration-drift balance is not informative or simply does not exist.

Causes of maladaptation

Evolutionary biologists tend to approach the study of the natural world within a framework of adaptation, inspired perhaps by the power of natural selection to produce fitness advantages that drive population persistence and biological diversity. In contrast, evolution has rarely been studied through the lens of adaptation's complement, maladaptation. This contrast is surprising because maladaptation is a prevalent feature of evolution: population trait values are rarely distributed optimally; local populations often have lower fitness than imported ones; populations decline; and local and global extinctions are common. Yet we lack a general framework for understanding maladaptation; for instance in terms of distribution, severity, and dynamics. Similar uncertainties apply to the causes of maladaptation. We suggest that incorporating maladaptation-based perspectives into evolutionary biology would facilitate better understanding of the natural world. Approaches within a maladaptation framework might be especially profitable in applied evolution contexts - where reductions in fitness are common. Toward advancing a more balanced study of evolution, here we present a conceptual framework describing causes of maladaptation. As the introductory article for a Special Feature on maladaptation, we also summarize the studies in this Issue, highlighting the causes of maladaptation in each study. We hope that our framework and the papers in this Special Issue will help catalyze the study of maladaptation in applied evolution, supporting greater understanding of evolutionary dynamics in our rapidly changing world.

Transient local adaptation and source-sink dynamics in experimental populations experiencing spatially heterogeneous environments

Local adaptation is determined by the strength of selection and the level of gene flow within heterogeneous landscapes. The presence of benign habitat can act as an evolutionary stepping stone for local adaptation to challenging environments by providing the necessary genetic variation. At the same time, migration load from benign habitats will hinder adaptation. In a community context, interspecific competition is expected to select against maladapted migrants, hence reducing migration load and facilitating adaptation. As the interplay between competition and spatial heterogeneity on the joint ecological and evolutionary dynamics of populations is poorly understood, we performed an evolutionary experiment using the herbivore spider mite Tetranychus urticae as a model. We studied the species's demography and local adaptation in a challenging environment that consisted of an initial sink (pepper plants) and/or a more benign environment (cucumber plants). Half of the experimental populations were exposed to a competitor, the congeneric T. ludeni. We show that while spider mites only adapted to the challenging pepper environment when it was spatially interspersed with benign cucumber habitat, this adaptation was only temporary and disappeared when the populations in the benign cucumber environment were expanding and spilling-over to the challenging pepper environment. Although the focal species outcompeted the competitor after about two months, a negative effect of competition on the focal species's performance persisted in the benign environment. Adaptation to challenging habitat in heterogeneous landscapes thus highly depends on demography and source-sink dynamics, but also on competitive interactions with other species, even if they are only present for a short time span.

Advances in population ecology and species interactions in mammals

The study of mammals has promoted the development and testing of many ideas in contemporary ecology. Here we address recent developments in foraging and habitat selection, source–sink dynamics, competition (both within and between species), population cycles, predation (including apparent competition), mutualism, and biological invasions. Because mammals are appealing to the public, ecological insight gleaned from the study of mammals has disproportionate potential in educating the public about ecological principles and their application to wise management. Mammals have been central to many computational and statistical developments in recent years, including refinements to traditional approaches and metrics (e.g., capture-recapture) as well as advancements of novel and developing fields (e.g., spatial capture-recapture, occupancy modeling, integrated population models). The study of mammals also poses challenges in terms of fully characterizing dynamics in natural conditions. Ongoing climate change threatens to affect global ecosystems, and mammals provide visible and charismatic subjects for research on local and regional effects of such change as well as predictive modeling of the long-term effects on ecosystem function and stability. Although much remains to be done, the population ecology of mammals continues to be a vibrant and rapidly developing field. We anticipate that the next quarter century will prove as exciting and productive for the study of mammals as has the recent one.

Ants Response to Human-Induced Disturbance in a Rain Tropical Forest

A high rate of human-induced disturbance of tropical ecosystems results in enormous loss of biodiversity due to local extinctions. Yet, mechanisms at the population level that lead to the extinction are still poorly understood. Here we tested the hypothesis that human-induced disturbance results in smaller amount of nesting sites for wood-dwelling arthropods that leads to smaller population size and diminished reproduction, and therefore, may promote local extinctions. We completed censuses in less-disturbed and human-disturbed secondary rain forest plots in Puerto Rico. We measured population size and brood production in wood-nesting ants and examined whether these parameters differ between less-disturbed and more-disturbed habitats. In addition, we measured volume of wood parts of all inhabited and potential nesting sites to assess nest site availability. We found that more human-disturbed forests furnish smaller nest sites, resulting in diminished population size and lowered brood production. Our study shows that human-induced disturbance decreases volume of available nesting sites that leads to decreased population size and lowered reproduction. Thus, in addition to the well-documented loss of species richness in human-disturbed tropical habitats, we demonstrated the direct effect of the disturbance that may promote vulnerability of local populations.

A multispecies test of source-sink indicators to prioritize habitat for declining populations

For species at risk of decline or extinction in source-sink systems, sources are an obvious target for habitat protection actions. However, the way in which source habitats are identified and prioritized can reduce the effectiveness of conservation actions. Although sources and sinks are conceptually defined using both demographic and movement criteria, simplifications are often required in systems with limited data. To assess the conservation outcomes of alternative source metrics and resulting prioritizations, we simulated population dynamics and extinction risk for 3 endangered species. Using empirically based habitat population models, we linked habitat maps with measured site- or habitat-specific demographic conditions, movement abilities, and behaviors. We calculated source-sink metrics over a range of periods of data collection and prioritized consistently high-output sources for conservation. We then tested whether prioritized patches identified the habitats that most affected persistence by removing them and measuring the population response. Conservation decisions based on different source-sink metrics and durations of data collection affected species persistence. Shorter time series obscured the ability of metrics to identify influential habitats, particularly in temporally variable and slowly declining populations. Data-rich source-sink metrics that included both demography and movement information did not always identify the habitats with the greatest influence on extinction risk. In some declining populations, patch abundance better predicted influential habitats for short-term regional persistence. Because source-sink metrics (i.e., births minus deaths; births and immigrations minus deaths and emigration) describe net population conditions and cancel out gross population counts, they may not adequately identify influential habitats in declining populations. For many nonequilibrium populations, new metrics that maintain the counts of individual births, deaths, and movement may provide additional insight into habitats that most influence persistence.

Constraints Imposed by a Natural Landscape Override Offspring Fitness Effects to Shape Oviposition Decisions in Wild Forked Fungus Beetles

Oviposition site decisions often maximize offspring fitness, but costs constraining choice can cause females to oviposit in poor developmental environments. It is unclear whether these constraints cumulatively outweigh offspring fitness to determine oviposition decisions in wild populations. Understanding how constraints shape oviposition in natural landscapes is a critical step toward revealing how maternal behavior influences fundamental phenomena like the evolution of specialization and the use of sink environments. Here, we used a genetic capture-recapture technique to reconstruct the oviposition decisions of individual females in a natural metapopulation of a beetle (Bolitotherus cornutus) that oviposits on three fungus species. We measured larval fitness-related traits (mass, development time, survival) on each fungus and compared the oviposition preferences of females in laboratory versus field tests. Larval fitness differed substantially among fungi, and females preferred a high-quality (high larval fitness) fungus in laboratory trials. However, females frequently laid eggs on the lowest-quality fungus in the wild. They preferred high-quality fungi when moving between oviposition sites, but this preference disappeared as the distance between sites increased and was inconsistent between study plots. Our results suggest that constraints on oviposition preferences in natural landscapes are sufficiently large to drive oviposition in poor developmental environments even when offspring fitness consequences are severe.

Spatially varying density dependence drives a shifting mosaic of survival in a recovering apex predator (Canis lupus)

Understanding landscape patterns in mortality risk is crucial for promoting recovery of threatened and endangered species. Humans affect mortality risk in large carnivores such as wolves (Canis lupus), but spatiotemporally varying density dependence can significantly influence the landscape of survival. This potentially occurs when density varies spatially and risk is unevenly distributed. We quantified spatiotemporal sources of variation in survival rates of gray wolves (C. lupus) during a 21‐year period of population recovery in the Upper Peninsula of Michigan, USA. We focused on mapping risk across time using Cox Proportional Hazards (CPH) models with time‐dependent covariates, thus exploring a shifting mosaic of survival. Extended CPH models and time‐dependent covariates revealed influences of seasonality, density dependence and experience, as well as individual‐level factors and landscape predictors of risk. We used results to predict the shifting landscape of risk at the beginning, middle, and end of the wolf recovery time series. Survival rates varied spatially and declined over time. Long‐term change was density‐dependent, with landscape predictors such as agricultural land cover and edge densities contributing negatively to survival. Survival also varied seasonally and depended on individual experience, sex, and resident versus transient status. The shifting landscape of survival suggested that increasing density contributed to greater potential for human conflict and wolf mortality risk. Long‐term spatial variation in key population vital rates is largely unquantified in many threatened, endangered, and recovering species. Variation in risk may indicate potential for source‐sink population dynamics, especially where individuals preemptively occupy suitable territories, which forces new individuals into riskier habitat types as density increases. We encourage managers to explore relationships between adult survival and localized changes in population density. Density‐dependent risk maps can identify increasing conflict areas or potential habitat sinks which may persist due to high recruitment in adjacent habitats.

More than 75 percent decline over 27 years in total flying insect biomass in protected areas

Global declines in insects have sparked wide interest among scientists, politicians, and the general public. Loss of insect diversity and abundance is expected to provoke cascading effects on food webs and to jeopardize ecosystem services. Our understanding of the extent and underlying causes of this decline is based on the abundance of single species or taxonomic groups only, rather than changes in insect biomass which is more relevant for ecological functioning. Here, we used a standardized protocol to measure total insect biomass using Malaise traps, deployed over 27 years in 63 nature protection areas in Germany (96 unique location-year combinations) to infer on the status and trend of local entomofauna. Our analysis estimates a seasonal decline of 76%, and mid-summer decline of 82% in flying insect biomass over the 27 years of study. We show that this decline is apparent regardless of habitat type, while changes in weather, land use, and habitat characteristics cannot explain this overall decline. This yet unrecognized loss of insect biomass must be taken into account in evaluating declines in abundance of species depending on insects as a food source, and ecosystem functioning in the European landscape.