Recursive movement patterns have been used to detect behavioral structure within individual movement trajectories in the context of foraging ecology, home-ranging behavior, and predator avoidance. Some animals exhibit movement recursions to locations that are tied to reproductive functions, including nests and dens; while existing literature recognizes that, no method is currently available to explicitly target different types of revisited locations. Moreover, the temporal persistence of recursive movements to a breeding location can carry information regarding the fate of breeding attempts, but it has never been used as a metric to quantify recursive movement patterns. Here, we introduce a method to locate breeding attempts and estimate their fate from GPS-tracking data of central place foragers. We tested the performance of our method in three bird species differing in breeding ecology (wood stork ( Mycteria americana) , lesser kestrel ( Falco naumanni ), Mediterranean gull ( Ichthyaetus melanocephalus )) and implemented it in the R package ‘nestR’. We identified breeding sites based on the analysis of recursive movements within individual tracks. Using trajectories with known breeding attempts, we estimated a set of species-specific criteria for the identification of nest sites, which we further validated using non-reproductive individuals as controls. We then estimated individual nest survival as a binary measure of reproductive fate (success, corresponding to fledging of at least one chick, or failure) from nest-site revisitation histories during breeding attempts, using a Bayesian hierarchical modeling approach that accounted for temporally variable revisitation patterns, probability of visit detection, and missing data. Across the three species, positive predictive value of the nest-site detection algorithm varied between 87-100% and sensitivity between 88-92%, and we correctly estimated the fate of 86-100% breeding attempts. By providing a method to formally distinguish among revisited locations that serve different ecological functions and introducing a probabilistic framework to quantify temporal persistence of movement recursions, we demonstrated how the analysis of recursive movement patterns can be applied to estimate reproduction in central place foragers. Beyond avian species, the principles of our method can be applied to other central place foraging breeders such as denning mammals. Our method estimates a component of individual fitness from movement data and will help bridge the gap between movement behavior, environmental factors, and their fitness consequences.* GPS : Global Positioning System CART : Classification And Regression Trees MCMC : Markov Chain Monte Carlo JAGS : Just Another Gibbs Sampler GLM : Generalized Linear Model

A major aim of evolutionary biology is to understand why patterns of genomic variation vary among populations and species. Large-scale genomic studies of widespread species are useful for studying how the environment and demographic history shape patterns of genomic divergence, and with the continually decreasing cost of sequencing, such studies are now becoming feasible. Here, we carry out one of the most comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 28 populations, spanning almost the entire geographic range of the European great tit subspecies. We found that genome-wide variation was consistent with a recent colonisation across Europe from a single refugium in the Balkans and/or Turkey, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation' even among populations which are ostensibly panmictic. These islands of differentiation were consistently found in regions of low recombination, suggesting that background selection can rapidly promote population differentiation among even the most recently colonised populations. We also detected outlier regions that were unique to peripheral great tit populations, most likely as a result of recent directional selection at the range edges of this species. These 'unique' outlier regions contained candidate genes for morphology, stress response and colouration, supporting previous research in this species, and providing avenues for future investigation. Our study suggests that comprehensive screens of genomic variation in wild organisms can provide unique insights into evolution.

Abstract The effectiveness of agri‐environment schemes (AESs), the largest conservation‐related expenditure for farmland biodiversity conservation within the European Union, is often compromised by a limited spatial scale of implementation. We focused on multiannual forage crops, a surrogate habitat for grassland birds, to assess the scale‐dependent effects of mowing timing and frequency on the local population size of an iconic species, the skylark ( Alauda arvensis ). While there is much evidence for a negative impact of in‐field mowing activities on grassland birds, whether such effects occur also at broader spatial scales is largely unknown. We surveyed breeding skylarks in the Po Plain (northern Italy) to determine (1) the association between landscape composition/configuration and abundance and (2) how abundance is affected by forage crop mowing timing and frequency. We addressed both questions through scale optimisation, identifying the most influential spatial scales for each covariate. Forage crop mowing timing was assessed through a novel remote sensing algorithm based on high‐resolution Sentinel‐2 satellite images. We observed a strong scale dependence on the importance of different habitats in determining skylark abundance. Abundance increased with an increasing cover of forage crops locally (200 m) and of winter crops at a landscape scale (2600 m), suggesting that the species is favoured by heterogeneous agroecosystems. Locally (150–350 m), skylarks were more abundant when crops were aggregated, being negatively impacted by crop fragmentation caused by urbanization and by seminatural habitats. At the landscape scale (1150 m), the timing of mowing was consistent across years, with early‐mown areas supporting fewer skylarks. This is probably because, over longer temporal scales, early‐mown forage patches have limited or null productivity, eventually limiting local population size. Synthesis and applications . We provide a new perspective on the overarching influence of spatial scale in driving the abundance of a declining farmland bird species, supporting the urgency of designing landscape scale‐effective AESs. This should be framed within the new EU Common Agricultural Policy reform and operated by farmer collectives, whereby management interventions should be monitored by state‐of‐the‐art remote sensing techniques. These results suggest that implementing scale‐optimized AESs could be crucial for effective farmland biodiversity conservation.

Abstract Per- and polyfluoroalkyl substances (PFAS) are synthetic organofluorine compounds used in various products, which are highly durable in the environment and may pose risks to wildlife health. We investigated the blood cell concentrations of PFAS in breeding Scopoli's shearwaters (Calonectris diomedea) from three different colonies in the central and southern Mediterranean (Linosa, Malta, and La Maddalena). Shearwaters are flexible, high trophic level foragers, and foraging areas may differ according to sex and breeding stage. We examined inter- and intracolony differences in PFAS blood concentrations and compared them with exploited foraging areas and dietary tracers. Per- and polyfluoroalkyl substances were detected in all samples, with the major congeners detected in descending order being perfluoroctanesulfonic acid (PFOS), perfluoroundecanoic acid (PFuNA), perfluorododecanoic acid (PFDoDA), and perfluorotridecanoic acid (PFTriDA). The mean sum of PFAS during the chick-rearing phase was highest in the birds from Malta (145.1 ng/g dry wt, 95% confidence interval [CI] of the mean 106.8, 183.5) compared with Linosa (91.5 ng/g dry wt, 95% CI 72.9, 110.1) and La Maddalena (84.5 ng/g dry wt, 95% CI 61.7, 107.3), and the PFAS blood composition of shearwaters from La Maddalena and Malta differed. The PFAS concentrations in shearwaters from Linosa were higher during incubation than during chick-rearing, and males had higher PFAS concentrations than females during incubation. Some PFAS were associated with carbon and nitrogen stable isotope values. After baseline adjustment of stable isotope values, no differences were observed for adjusted δ15N and δ13C between the three colonies, suggesting that differences in PFAS levels attributed to diet were minor compared with regional differences. Our study highlights that shearwaters are useful biomonitors of PFAS exposure in remote marine areas.