Search

Artificial light at night (ALAN) is increasing in extent and intensity across the globe. It has been shown to interfere with animal sensory systems, orientation and distribution, with the potential to cause significant ecological impacts. We analysed the locations of 102 mountain lions (Puma concolor) in a light-polluted region in California. We modelled their distribution relative to environmental and human-disturbance variables, including upward radiance (nearby lights), zenith brightness (sky glow) and natural illumination from moonlight. We found that mountain lion probability of presence was highly related to upward radiance, that is, related to lights within approximately 500 m. Despite a general pattern of avoidance of locations with high upward radiance, there were large differences in degree of avoidance among individuals. The amount of light from artificial sky glow was not influential when included together with upward radiance in the models, and illumination from moonlight was not influential at all. Our results suggest that changes in visibility associated with lunar cycles and sky glow are less important for mountain lions in their selection of light landscapes than avoiding potential interactions with humans represented by the presence of nearby lights on the ground. This article is part of the theme issue Light pollution in complex ecological systems.

Dispersal is a key life-history trait for most species and is essential to ensure connectivity and gene flow between populations and facilitate population viability in variable environments. Despite the increasing importance of range shifts due to global change, dispersal has proved difficult to quantify, limiting empirical understanding of this phenotypic trait and wider synthesis. Here, we introduce a statistical framework to estimate standardised dispersal kernels from biased data. Based on this, we compare empirical dispersal kernels for European breeding birds considering age (average dispersal; natal, before first breeding; and breeding dispersal, between subsequent breeding attempts) and sex (females and males) and test whether different dispersal properties are phylogenetically conserved. We standardised and analysed data from an extensive volunteer-based bird ring-recoveries database in Europe (EURING) by accounting for biases related to different censoring thresholds in reporting between countries and to migratory movements. Then, we fitted four widely used probability density functions in a Bayesian framework to compare and provide the best statistical descriptions of the different age and sex-specific dispersal kernels for each bird species. The dispersal movements of the 234 European bird species analysed were statistically best explained by heavy-tailed kernels, meaning that while most individuals disperse over short distances, long-distance dispersal is a prevalent phenomenon in almost all bird species. The phylogenetic signal in both median and long dispersal distances estimated from the best-fitted kernel was low (Pagel's lambda < 0.25), while it reached high values (Pagel's lambda >0.7) when comparing dispersal distance estimates for fat-tailed dispersal kernels. As expected in birds, natal dispersal was on average 5 km greater than breeding dispersal, but sex-biased dispersal was not detected. Our robust analytical framework allows sound

Dispersal is a key life-history trait for most species and is essential to ensure connectivity and gene flow between populations and facilitate population viability in variable environments. Despite the increasing importance of range shifts due to global change, dispersal has proved difficult to quantify, limiting empirical understanding of this phenotypic trait and wider synthesis. Here, we introduce a statistical framework to estimate standardised dispersal kernels from biased data. Based on this, we compare empirical dispersal kernels for European breeding birds considering age (average dispersal; natal, before first breeding; and breeding dispersal, between subsequent breeding attempts) and sex (females and males) and test whether different dispersal properties are phylogenetically conserved. We standardised and analysed data from an extensive volunteer-based bird ring-recoveries database in Europe (EURING) by accounting for biases related to different censoring thresholds in reporting between countries and to migratory movements. Then, we fitted four widely used probability density functions in a Bayesian framework to compare and provide the best statistical descriptions of the different age and sex-specific dispersal kernels for each bird species. The dispersal movements of the 234 European bird species analysed were statistically best explained by heavy-tailed kernels, meaning that while most individuals disperse over short distances, long-distance dispersal is a prevalent phenomenon in almost all bird species. The phylogenetic signal in both median and long dispersal distances estimated from the best-fitted kernel was low (Pagel's lambda < 0.25), while it reached high values (Pagel's lambda >0.7) when comparing dispersal distance estimates for fat-tailed dispersal kernels. As expected in birds, natal dispersal was on average 5 km greater than breeding dispersal, but sex-biased dispersal was not detected. Our robust analytical framework allows sound

Poaching is driving many species toward extinction, and as a result, lowering poaching pressure is a conservation priority. This requires understanding where poaching pressure is high and which factors determine these spatial patterns. However, the cryptic and illegal nature of poaching makes this difficult. Ranger patrol data, typically recorded in protected area logbooks, contain information on patrolling efforts and poaching detection and should thus provide opportunities for a better understanding of poaching pressure. However, these data are seldom analyzed and rarely used to inform adaptive management strategies. We developed a novel approach to making use of analog logbook records to map poaching pressure and to test environmental criminology and predator–prey relationship hypotheses explaining poaching patterns. We showcase this approach for Golestan National Park in Iran, where poaching has substantially depleted ungulate populations. We digitized data from >4800 ranger patrols from 2014 to 2016 and used an occupancy modeling framework to relate poaching to (1) accessibility, (2) law enforcement, and (3) prey availability factors. Based on predicted poaching pressure and patrolling intensity, we provide suggestions for future patrol allocation strategies. Our results revealed a low probability (12%) of poacher detection during patrols. Poaching distribution was best explained by prey availability, indicating that poachers target areas with high concentrations of ungulates. Poaching pressure was estimated to be high (>0.49) in 39% of our study area. To alleviate poaching pressure, we recommend ramping up patrolling intensity in 12% of the national park, which could be achievable by reducing excess patrols in about 20% of the park. However, our results suggest that for 27% of the park, it is necessary to improve patrolling quality to increase detection probability of poaching, for example, by closing temporal patrolling gaps or expanding informant networks. O, Peer Reviewed

Effective conservation of animal populations depends on the availability of reliable data derived from rigorous monitoring protocols, which allows us to assess trends and understand the processes they are governed by. Nevertheless, population monitoring schemes are hampered by multiple sources of errors resulting from specific logistical and survey constraints. Two common complications are the non-visitation of some sites in certain years and preferential sampling (PS), that is, the tendency to survey “better” sites disproportionately more often. Both factors can lead to serious biases unless accommodated into models. We used 22 yr of nest-monitoring data to develop a dynamic multistate occupancy model, including a PS component to investigate occupancy and reproduction dynamics in a peripheral Black Stork (Ciconia nigra) population in Spain. We analyzed the effects of climate and nesting substrate (tree vs. cliff) on population dynamics and accounted for PS and non-visitation biases using a model that distinguished three territorial states: unoccupied, occupied without, or occupied with successful reproduction. We found strong evidence for positive PS, and when accounting for this bias, lower population size estimates were generated. Black stork nests had a high probability of remaining in the same state from one year to the next, with successful nests more likely to be occupied again and to be successful the following year than occupied but unsuccessful or unoccupied nests. Nesting substrate and spring precipitation did not influence state transition probabilities or the probability of reproductive success; nevertheless, cliff nest occupancy was overall higher than tree nest occupancy. Our results highlight the importance of correcting for non-visitation and PS in habitat occupancy models. If these potential biasing effects are not accounted for, inferences of population size may be overestimate. Multistate occupancy models with correction for PS offer a powerful a, Peer Reviewed

This article was supported by the German Research Foundation (DFG) and the Open Access Publication Fund of Humboldt-Universität zu Berlin., Effective conservation of animal populations depends on the availability of reliable data derived from rigorous monitoring protocols, which allows us to assess trends and understand the processes they are governed by. Nevertheless, population monitoring schemes are hampered by multiple sources of errors resulting from specific logistical and survey constraints. Two common complications are the non-visitation of some sites in certain years and preferential sampling (PS), that is, the tendency to survey “better” sites disproportionately more often. Both factors can lead to serious biases unless accommodated into models. We used 22 yr of nest-monitoring data to develop a dynamic multistate occupancy model, including a PS component to investigate occupancy and reproduction dynamics in a peripheral Black Stork (Ciconia nigra) population in Spain. We analyzed the effects of climate and nesting substrate (tree vs. cliff) on population dynamics and accounted for PS and non-visitation biases using a model that distinguished three territorial states: unoccupied, occupied without, or occupied with successful reproduction. We found strong evidence for positive PS, and when accounting for this bias, lower population size estimates were generated. Black stork nests had a high probability of remaining in the same state from one year to the next, with successful nests more likely to be occupied again and to be successful the following year than occupied but unsuccessful or unoccupied nests. Nesting substrate and spring precipitation did not influence state transition probabilities or the probability of reproductive success; nevertheless, cliff nest occupancy was overall higher than tree nest occupancy. Our results highlight the importance of correcting for non-visitation and PS in habitat occupancy models. If these potential biasing effects are not accounted for, inferences of population size may be overestimate. Multistate occupancy models with correction for PS offer a powerful a, Peer Reviewed

[Background and Aims]: It is widely accepted that changes in the environment affect mean trait expression, but little is known regarding how the environment shapes intra-individual and intra-population variance. Theory suggests that intra-individual variance might be plastic and under natural selection, rather than reflecting developmental noise, but evidence for this hypothesis is scarce. Here, we experimentally tested whether differences in intrinsic environmental predictability affect intra-individual and intra-population variability of different reproductive traits, and whether intra-individual variability is under selection., [Methods]: Under field-conditions, we subjected Onobrychis viciifolia to more and less predictable precipitation over four generations and four years. We analysed effects on the coefficient of intra-individual (CVi-i) and the coefficient of intra-population variation (CVi-p), assessed whether the coefficients of intra-individual variation (CsVi-i) are under natural selection and tested for transgenerational responses (ancestor environmental effects on offspring)., [Key Results]: Less predictable precipitation let to higher CsVi-i and CsVi-p, consistent with plastic responses. The CsVi-i of all studied traits were under consistent stabilizing selection, and precipitation predictability affected the strength of selection and the location of the optimal CVi-i of a single trait. All CsVi-i differed from the optimal CVi-i and the maternal and offspring CVi-i were positively correlated, showing that there was scope for change. Nevertheless, no consistent transgenerational effects were found in any of the three descendant generations, which contrasts with recent studies that detected rapid transgenerational responses in trait means of different plant species. This suggests that changes in intra-individual variability take longer to evolve than changes in trait means, which may explain why high intra-individual variability is maintained, despite the stabilizing selection., [Conclusions]: The results indicate that plastic changes of intra-individual variability are an important determinant of whether plants will be able to cope with changes in environmental predictability induced by the currently observed climatic change.

Global biodiversity is under high and rising anthropogenic pressure. Yet, how the taxonomic, phylogenetic, and functional facets of biodiversity are affected by different threats over time is unclear. This is particularly true for the two main drivers of the current biodiversity crisis: habitat destruction and overexploitation. We provide the first long-term assessment of multifaceted biodiversity changes caused by these threats for any tropical region. Focussing on larger mammals in South America's 1.1 million km Gran Chaco region, we assessed changes in multiple biodiversity facets between 1985 and 2015, determined which threats drive those changes, and identified remaining key areas for all biodiversity facets. Using habitat and threat maps, we found, first, that between 1985 and 2015 taxonomic (TD), phylogenetic (PD) and functional (FD) diversity all declined drastically across over half of the area assessed. FD declined about 50% faster than TD and PD, and these declines were mainly driven by species loss, rather than species turnover. Second, habitat destruction, hunting, and both threats together contributed ~57%, ~37%, and ~6% to overall facet declines, respectively. However, hunting pressure increased where TD and PD declined most strongly, whereas habitat destruction disproportionally contributed to FD declines. Third, just 23% of the Chaco would have to be protected to safeguard the top 17% of all three facets. Our findings uncover a widespread impoverishment of mammal species richness, evolutionary history, and ecological functions across broad areas of the Chaco due to increasing habitat destruction and hunting. Moreover, our results pinpoint key areas that should be preserved and managed to maintain all facets of mammalian diversity across the Chaco. More generally, our work highlights how long-term changes in biodiversity facets can be assessed and attributed to specific threats, to better understand human impacts on biodiversity and to guide conserva

This paper explores how environmental traits shape the winter distribution of passerine (O. Passeriformes) richness in the south-western Palearctic, a major wintering ground for the European avifauna. We apply three complementary approaches. First, we assess the way climate, landscape and habitat affect species richness by means of field counts across the study area. Second, we model the spatial distribution of six common passerines using ring recoveries as presence data. Finally, where the resulting models predict the actual distribution of birds, we employ the models to forecast the future distribution of richness according to predictions of climate change, i.e. of increasing temperature and decreasing precipitation. The results support an effect of landscape, habitat structure, temperature and precipitation on bird richness. Distribution models fit the actual distribution of bird richness and their predictions suggest that winter species richness will decrease in lowlands and increase in highlands. This pattern could be related to direct effects of temperature on thermoregulatory costs and indirect effects on winter primary productivity affecting food resources. These results also indicate that, in a context of climate warming, species-rich sectors will shift from lowlands to highlands, suggesting that uplands will have a regulatory role in the future winter distribution of birds. - Tellería, J.L., Fandos, G. & Fernández-López, J. (2021). Winter bird richness distribution in the south-western Palearctic: current patterns and potential changes. Ardeola, 68: 17-32., Este trabajo explora el efecto de los factores ambientales que determinan la distribución invernal de la riqueza de especies de paseriformes (O. Passeriformes) en el Paleártico suroccidental, un importante área de invernada para la avifauna europea. Para ello, aplicamos tres aproximaciones complementarias. Primero, evaluamos el efecto del clima, el hábitat y la estructura del paisaje sobre la riqueza de especies evaluada sobre el terreno en el área de estudio. Segundo, modelamos la distribución de seis especies utilizando las recuperaciones de aves anilladas. Finalmente, si los modelos predicen su distribución real, serán utilizados para predecir los cambios en la distribución de la riqueza de especies según las predicciones del cambio climático (aumento de la temperatura y descenso de la precipitación). Los modelos generados por las recuperaciones de aves anilladas se ajustan bien a la distribución real de la riqueza de especies y sus predicciones sugieren que este parámetro disminuirá en las zonas bajas y aumentará en las más elevadas. Este patrón puede relacionarse con el efecto de la temperatura sobre el coste de termorregulación y la productividad primaria invernal que, a su vez, afecta a la disponibilidad de alimento para las aves. Los resultados sugieren que, en un contexto de cambio climático, los sectores más ricos en especies se desplazarán de las zonas más bajas hacia las altas lo que apunta un papel regulador de las grandes mesetas del Paleártico suroccidental sobre los cambios en la distribución invernal de las aves.

Gaps in the forest canopy appear to play a major role in local species richness as they increase habitat heterogeneity and border effects. In the case of bats, these small clearings seem to play an important role as foraging and commuting sites. However, further research is required to set the actual role of forest gaps on bat biology in order to consider them as conservation targets in forest management. In this study, we test whether bat species richness and activity are higher in gaps within a large Scot pine (Pinus sylvestris) forest than in the adjacent tree covered control sites. We also test these differences in the forest specialists since this group is particularly sensitive to forest management. To investigate this, we used bat detectors in 9 gaps and their adjacent control points during ten nights in July - August 2016 and 2017. The trends resulting from the analysis of 228,108 bat calls support that species richness and activity of bat assemblages were higher in gaps than in the adjacent control points. A similar pattern was detected in the most sensitive tree dwelling bat species. These results stress the importance of gaps as useable areas for the forest bat assemblages and suggest the interest of managing gaps for bat conservation within large areas of continuous tree cover., Comunidad de Madrid/Fondo Social Europeo, Depto. de Biodiversidad, Ecología y Evolución, Fac. de Ciencias Biológicas, TRUE, pub

Global biodiversity is under high and rising anthropogenic pressure. Yet, how the taxonomic, phylogenetic, and functional facets of biodiversity are affected by different threats over time is unclear. This is particularly true for the two main drivers of the current biodiversity crisis: habitat destruction and overexploitation. We provide the first long-term assessment of multifaceted biodiversity changes caused by these threats for any tropical region. Focussing on larger mammals in South America's 1.1 million km2 Gran Chaco region, we assessed changes in multiple biodiversity facets between 1985 and 2015, determined which threats drive those changes, and identified remaining key areas for all biodiversity facets. Using habitat and threat maps, we found, first, that between 1985 and 2015 taxonomic (TD), phylogenetic (PD) and functional (FD) diversity all declined drastically across over half of the area assessed. FD declined about 50% faster than TD and PD, and these declines were mainly driven by species loss, rather than species turnover. Second, habitat destruction, hunting, and both threats together contributed ~57%, ~37%, and ~6% to overall facet declines, respectively. However, hunting pressure increased where TD and PD declined most strongly, whereas habitat destruction disproportionally contributed to FD declines. Third, just 23% of the Chaco would have to be protected to safeguard the top 17% of all three facets. Our findings uncover a widespread impoverishment of mammal species richness, evolutionary history, and ecological functions across broad areas of the Chaco due to increasing habitat destruction and hunting. Moreover, our results pinpoint key areas that should be preserved and managed to maintain all facets of mammalian diversity across the Chaco. More generally, our work highlights how long-term changes in biodiversity facets can be assessed and attributed to specific threats, to better understand human impacts on biodiversity and to guide conserv, Peer Reviewed

Species distribution models (SDMs) constitute the most common class of models across ecology, evolution and conservation. The advent of ready‐to‐use software packages and increasing availability of digital geoinformation have considerably assisted the application of SDMs in the past decade, greatly enabling their broader use for informing conservation and management, and for quantifying impacts from global change. However, models must be fit for purpose, with all important aspects of their development and applications properly considered. Despite the widespread use of SDMs, standardisation and documentation of modelling protocols remain limited, which makes it hard to assess whether development steps are appropriate for end use. To address these issues, we propose a standard protocol for reporting SDMs, with an emphasis on describing how a study's objective is achieved through a series of modeling decisions. We call this the ODMAP (Overview, Data, Model, Assessment and Prediction) protocol, as its components reflect the main steps involved in building SDMs and other empirically‐based biodiversity models. The ODMAP protocol serves two main purposes. First, it provides a checklist for authors, detailing key steps for model building and analyses, and thus represents a quick guide and generic workflow for modern SDMs. Second, it introduces a structured format for documenting and communicating the models, ensuring transparency and reproducibility, facilitating peer review and expert evaluation of model quality, as well as meta‐analyses. We detail all elements of ODMAP, and explain how it can be used for different model objectives and applications, and how it complements efforts to store associated metadata and define modelling standards. We illustrate its utility by revisiting nine previously published case studies, and provide an interactive web‐based application to facilitate its use. We plan to advance ODMAP by encouraging its further refinement and adoption by the sc, Peer Reviewed

This article was supported by the German Research Foundation (DFG) and the Open Access Publication Fund of Humboldt-Universität zu Berlin., Species distribution models (SDMs) constitute the most common class of models across ecology, evolution and conservation. The advent of ready‐to‐use software packages and increasing availability of digital geoinformation have considerably assisted the application of SDMs in the past decade, greatly enabling their broader use for informing conservation and management, and for quantifying impacts from global change. However, models must be fit for purpose, with all important aspects of their development and applications properly considered. Despite the widespread use of SDMs, standardisation and documentation of modelling protocols remain limited, which makes it hard to assess whether development steps are appropriate for end use. To address these issues, we propose a standard protocol for reporting SDMs, with an emphasis on describing how a study's objective is achieved through a series of modeling decisions. We call this the ODMAP (Overview, Data, Model, Assessment and Prediction) protocol, as its components reflect the main steps involved in building SDMs and other empirically‐based biodiversity models. The ODMAP protocol serves two main purposes. First, it provides a checklist for authors, detailing key steps for model building and analyses, and thus represents a quick guide and generic workflow for modern SDMs. Second, it introduces a structured format for documenting and communicating the models, ensuring transparency and reproducibility, facilitating peer review and expert evaluation of model quality, as well as meta‐analyses. We detail all elements of ODMAP, and explain how it can be used for different model objectives and applications, and how it complements efforts to store associated metadata and define modelling standards. We illustrate its utility by revisiting nine previously published case studies, and provide an interactive web‐based application to facilitate its use. We plan to advance ODMAP by encouraging its further refinement and adoption by the sc, Peer Reviewed

Global biodiversity is under high and rising anthropogenic pressure. Yet, how the taxonomic, phylogenetic, and functional facets of biodiversity are affected by different threats over time is unclear. This is particularly true for the two main drivers of the current biodiversity crisis: habitat destruction and overexploitation. We provide the first long-term assessment of multifaceted biodiversity changes caused by these threats for any tropical region. Focussing on larger mammals in South America's 1.1 million km2 Gran Chaco region, we assessed changes in multiple biodiversity facets between 1985 and 2015, determined which threats drive those changes, and identified remaining key areas for all biodiversity facets. Using habitat and threat maps, we found, first, that between 1985 and 2015 taxonomic (TD), phylogenetic (PD) and functional (FD) diversity all declined drastically across over half of the area assessed. FD declined about 50% faster than TD and PD, and these declines were mainly driven by species loss, rather than species turnover. Second, habitat destruction, hunting, and both threats together contributed ~57%, ~37%, and ~6% to overall facet declines, respectively. However, hunting pressure increased where TD and PD declined most strongly, whereas habitat destruction disproportionally contributed to FD declines. Third, just 23% of the Chaco would have to be protected to safeguard the top 17% of all three facets. Our findings uncover a widespread impoverishment of mammal species richness, evolutionary history, and ecological functions across broad areas of the Chaco due to increasing habitat destruction and hunting. Moreover, our results pinpoint key areas that should be preserved and managed to maintain all facets of mammalian diversity across the Chaco. More generally, our work highlights how long-term changes in biodiversity facets can be assessed and attributed to specific threats, to better understand human impacts on biodiversity and to guide conserv, Peer Reviewed

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (i) drift in geographic isolation, (ii) differentiation along continuous selective gradients, and (iii) isolation by adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

[Background] Space use by animals is determined by the interplay between movement and the environment, and is thus mediated by habitat selection, biotic interactions and intrinsic factors of moving individuals. These processes ultimately determine home range size, but their relative contributions and dynamic nature remain less explored. We investigated the role of habitat selection, movement unrelated to habitat selection and intrinsic factors related to sex in driving space use and home range size in Iberian ibex, Capra pyrenaica. We used GPS collars to track ibex across the year in two different geographical areas of Sierra Nevada, Spain, and measured habitat variables related to forage and roost availability., [Results] By using integrated step selection analysis (iSSA), we show that habitat selection was important to explain space use by ibex. As a consequence, movement was constrained by habitat selection, as observed displacement rate was shorter than expected under null selection. Selection-independent movement, selection strength and resource availability were important drivers of seasonal home range size. Both displacement rate and directional persistence had a positive relationship with home range size while accounting for habitat selection, suggesting that individual characteristics and state may also affect home range size. Ibex living at higher altitudes, where resource availability shows stronger altitudinal gradients across the year, had larger home ranges. Home range size was larger in spring and autumn, when ibex ascend and descend back, and smaller in summer and winter, when resources are more stable. Therefore, home range size decreased with resource availability. Finally, males had larger home ranges than females, which might be explained by differences in body size and reproductive behaviour., [Conclusions] Movement, selection strength, resource availability and intrinsic factors related to sex determined home range size of Iberian ibex. Our results highlight the need to integrate and account for process dependencies, here the interdependence of movement and habitat selection, to understand how animals use space. This study contributes to understand how movement links environmental and geographical space use and determines home range behaviour in large herbivores.