How Individual Variation Shapes Ecological Niches in Two <i>Pipistrellus</i> Bat Species

Niche partitioning is a key mechanism explaining species coexistence and biodiversity. While different axes of niche partitioning can account for species coexistence, variation between and within individuals can significantly influence this process. This empirical study investigates niche partitioning in a comprehensive framework using two co-occurring bat species as model organisms: Nathusius’ pipistrelle <i>Pipistrellus nathusii</i> and Kuhl’s pipistrelle <i>P. kuhlii.</i> <i>Pipistrellus kuhlii</i> is a recent invader of the study area, rapidly expanding its geographic range. As human-induced alterations in land cover and climate are likely to accelerate this range expansion, understanding the mechanisms enabling the co-occurrence of these species is crucial. We examined niche partitioning across habitat, time, and resource use, considering variation between species, between-individual, and within-individual levels. Specifically, we analysed the movement patterns of 58 individuals of both species and employed metabarcoding on guano samples to assess resource use. Our results demonstrate that individual variation exceeded species-level differences and varied between species.<i>Pipistrellus nathusii</i> exhibited greater between-individual variation, while the expanding <i>P. kuhlii</i> displayed stronger within-individual variation, which may contribute to its range expansion. Our study suggests an important role of within-individual variation in the ongoing expansion of an invading bat species, reshaping animal communities in the context of global change.<br>Data and scripts of “The ecological niche and individual variation of a rapidly expanding bat”<br><b>a_TriangulationEstimation</b><br>This folder contains the R-scripts required to triangulate the pre-filtered data. First step is to pre-filter data based on the interval between consecutive signals. Then, the bearings are calculated using a cosine approach. Finally, the bearings are triangulated using Azimuthal Telemetry Models (ATM).The fits of the ATM have a size of 56.42 GB. Thus, they are not provided on Dryad, but can be shared on request.<br>Input: Raw dataoutput: Triangulated location estimates including inaccuracy estimates<br><br><b>b_filterRawPoints</b><br>This folder contains the scripts to filter the triangulated locations based on activity (moving vs roosting) and speed.<br>Input: Triangulated location estimates including inaccuracy estimatesOutput: Filtered location estimates<br><br><b>c_MovementModels</b><br>This folder contains the scripts to transform the location data to dynamic Brownian Bridge movement models and burst the data to 30 min, 10 min, and 1 hour intervals.<br>Input: Filtered location estimatesOutput: Movement models, burst to 30 min, 10 min, and 1 hour intervals<br><br><b>d_landscapePreparation</b><br>This folder contains the scripts to crop the land cover data and calculate selection ratios of the respective land cover types.<br>Input: land cover data (Corine, Land Burgenland), movement modelsOutput: Temporal explicit selection ratios<br><b>e_ nicheModels</b><br>This folder contains the scripts to build niche models based on the temporal explicit selection ratios<br>Input: Temporal explicit selection ratiosOutput: Individual niche models<br><b>f_metabarcoding</b><br>This folder contains the analysis of the metabarcoding.<br>Input: Output of metabarcoding analysisOutput: PERMANOVA of metabarcoding results, NMDS and Pianka’s niche overlap<br>