A new general index of home range overlap and segregation: the Lesser Kestrel in Southern Italy as a case study

There is increasing interest in evaluating home-range overlap (or, otherwise, segregation) between bird species, and between or within bird populations, to inform spatial planning. So far, studies of home-range overlap typically make use of comparisons between pairs of individuals, populations or species, and return a matrix of pairwise overlaps (e.g., percent overlaps). However, when the number of individuals, populations or species to be compared is elevated, an overlarge overlap matrix is difficult to interpret from an ecological viewpoint. We propose here a new, conceptually simple and computationally efficient index (general overlap index; GOI) for the ready computation within GIS of home range overlap of an arbitrarily large number (i.e., n ≥ 2) of individuals, populations or species. Whatever the number of home ranges to be compared, GOI always returns a single score between 0 and 100. As a case study, we applied our index to 24,074 GPS points of 10 Lesser Kestrels (Falco naumanni) in order to estimate within-colony and between-colony overlaps in two neighboring colonies in Southern Italy. Within-colony overlap was elevated for both colonies (96.41% at Cassano delle Murge, n = 5 individuals; 81.38% at Santeramo in Colle, n = 5 individuals), while between-colony overlap was low (19.12%; n = 2 colonies) and, after a randomization procedure, more spatially-segregated than expected by chance. Modern biotelemetry offers huge amounts of data describing the space use of animal species. The use of intuitive and straightforward indices, like GOI, can be useful to promptly extract ecological information from such an amount of data (e.g. detecting change in space use over successive years, evaluating the reliability of various home-range estimators).


Background
Intraspecific and interspecific interactions impact the extent and spatial patterning of animal home ranges (Adams 2001). Animals compete for resources (e.g. food, shelter and mates) and it has long been known (e.g., Brown and Orians 1970;Davies 1978;Newton 1998) that one way to compete is to exclude potential competitors from the area containing resources. Overlap, or otherwise segregation, between bird species (Warning and Benedict 2015;Zhao et al. 2015), and between or within bird populations (Yang et al. 2011;Clay et al. 2016) has been progressively recorded in recent years due to the increased availability of telemetry data (Wang et al. 2010). The increasing popularity of multi-species studies in the context of spatial management (e.g. Lascelles et al. 2016) has enhanced the need to calculate home range overlap for a large number of individuals, populations or species.

Open Access
Avian Research *Correspondence: sgtpm@libero.it Conservation Department, LIPU -BirdLife Italy (Lega Italiana Protezione Uccelli), Via Udine 3/a, 43122 Parma, Italy Home-range overlap indices have several important applications to wildlife research and management. Overlap indices can be useful for assessing the degree of interaction among individuals as well as site fidelity for a particular individual. In addition, overlap measures may be used to measure the reliability of various home-range estimators (Fieberg and Kochanny 2005). There are several overlap indices available in the literature (Kernohan et al. 2001). The available approaches make use of comparisons between pairs of individuals, populations or species, and return a matrix of pairwise overlaps. The most common and intuitive approach is percent overlap, i.e. the proportion of animal i's home range that is overlapped by animal j's home range (Kernohan et al. 2001). However, when the number of individuals, populations or species to be compared is elevated, an overlarge overlap matrix is difficult to interpret from an ecological viewpoint. For example, with only 10 individuals (populations, or species), a 10 × 10 pairwise overlap matrix is produced, whose ecological interpretation could be not that simple. Researchers have interpreted multiple pairwise comparisons using mean overlaps (e.g., Macias-Duarte and Panjabi 2013), however mean values can be scarcely representative of the pairwise overlap matrix if the dispersion of overlap values around the mean is elevated.
Accordingly, we propose here a new, conceptually simple and computationally efficient index (general overlap index; GOI hereafter) for the ready computation within GIS of home range overlap of an arbitrarily large number (i.e., n ≥ 2) of individuals, populations or species. As a case study, we applied our index to ten Lesser Kestrels (Falco naumanni) in order to estimate within-colony and between-colony overlap/segregation in two neighboring urban colonies (Cassano delle Murge and Santeramo in Colle; Apulia region) in Southern Italy. From an ecological point of view, it was a good study system to investigate competition during the breeding season as the two colonies represent the most elevated density of Lesser Kestrels in urban areas worldwide (Gustin et al. 2013).
The Lesser Kestrel is a colonial, small falcon breeding in steppe-like grasslands and cultivated landscapes with short vegetation and extensive crops (BirdLife International 2017). It is present among Annex I species of EU Wild Birds Directive (2009/147/EEC), and its important breeding habitats have been designated as Special Protection Areas of the Natura 2000 Network. In Southern Italy, the Lesser Kestrel has been widely studied in the nearby colonies of Gravina in Puglia and Altamura (Gustin et al. 2014a(Gustin et al. , b, 2017aGustin et al. 2014bGustin et al. , 2017bGustin et al. , 2018Ferrarini et al. 2018a, b;Ferrarini et al. 2018bFerrarini et al. , 2019. The study area ( Fig. 1) is an agricultural landscape located within the SPA (Special Protection Area) "Murgia Alta" IT9120007, and also included within the IBA (Important Bird Area) "Murge".

Methods
We monitored ten birds (5 at Cassano delle Murge and 5 at Santeramo in Colle; Table 1) during the nestling period (Additional file 1: Fig. S1) in the two urban colonies of Cassano delle Murge (between June 22th and July 6th 2017; 11,993 GPS points) and Santeramo in Colle (between June 13th and June 29th 2017; 12,081 GPS points) where Lesser Kestrels have their artificial nests (Additional file 1: Fig. S2). We tracked birds using TechnoSmart GiPSy-4 and GiPSy-5 data loggers (23 × 15 × 6 mm, 5 g weight; Additional file 1: Fig. S3), which provided information about date, time, latitude, longitude, altitude (m a.s.l.) and instantaneous speed (m/s). GPS sampling frequency was one fix every three minutes. We fitted birds with data loggers at their nest boxes when they were delivering food to nestlings. All devices were tied dorsally to the base of two central tail feathers (Additional file 1: Fig. S4). The weight of the devices in relation to that of the birds was less than 4% for all individuals. The attachment of transmitters did not take more than 15 min, and had no visible deleterious effects on the birds. To download the data from the data-loggers, we recaptured birds at their nest boxes.
We imported GPS data into GIS and estimated the individual home ranges using the fixed-mean minimum convex polygons (Kenward 1987) which calculates the arithmetic mean of all X (longitude) and Y (latitude) coordinates, then selects the requested percentage of points closest to that arithmetic mean point. We also estimated colony-specific home ranges (i.e. home ranges calculated after pooling the locations of all individuals of each colony). We chose the 95% isopleth to represent home range as this value is widely used in the literature (White and Garrott 1990).
In order to quantify home range overlaps, we first used the most common method (percent overlap; j is the proportion of home-range i that is overlapped by home-range j, A i is the area of home-range i, and A i,j is the area of overlap between the two homeranges. As HR i,j ≠ HR j,i (i.e., directional indices), we quantified the degree of overlap using both HR i,j and HR j,i . In addition, we also employed our general overlap index (GOI). In the case of perfectly disjoint (i.e. non-overlapping) home ranges (Fig. 2a), the total area (A T ) covered by the home ranges is simply the sum of their extents (i.e. A i ). In the case of perfectly nested (i.e. overlapping) home ranges (Fig. 2b), A T is simply the extent of the largest home range (i.e. max(A i )). In the intermediate case (i.e., partially overlapping home ranges; Fig. 2c), A T corresponds to the union of the home range polygons (i.e. A i ). Unioning a set of (partially) overlapping polygons is a standard GIS procedure with the effect of merging their areas (Fig. 2d). Therefore, the difference between A i and max(A i ) represents the maximum distance possible (Dist MAX ) from a perfectly  non-overlapping situation. The difference between A i and A i is the observed distance (Dist OBS ) from the perfectly disjoint situation. GOI was calculated as (Eq. 1): where n is the number of home ranges under study. GOI thus measures the distance of the observed overlaps from two extremes (perfect overlap and perfect non-overlap). If Dist OBS = 0, then GOI = 0 (perfect non-overlap); if Dist OBS = Dist MAX , then GOI = 100 (perfect overlap). If home ranges partially overlap, then 0 < GOI < 100. The pseudocode of the algorithm used to calculate GOI is described in the Additional file 2: Text S1.
Our overlap index corresponds, in essence, to the linear equation ) where a is the extent of the largest home range polygon, b is the sum of the home (1) range extents and X is the extent of the union of home range polygons, which varies depending upon the degree of overlap. Since b ≥ X and b > a, then the denominator is always positive while the numerator can be positive or null. In addition, since X ≥ a, then b -X is always less than, or equal to, ba. Thus GOI is constrained in the interval [0, 100], independently of the number of home ranges under study. In the case of perfectly non-overlapping home ranges, X = b then GOI = 0. In the case of perfectly overlapping home ranges, X = a then GOI = 100. Finally, a general segregation index (GSI) was computed as the complement to 100 of GOI (Eq. 2): The first derivatives of GOI and GSI (Eqs. 3, 4) show their rate of change with respect to A i : . 2 a Perfectly disjoint (i.e., non-overlapping) home ranges, b perfectly nested (i.e., overlapping) home ranges, c partially-overlapping home ranges, d union (on the right) of partially-overlapping home range polygons (on the left) Ferrarini et al. Avian Res (2021) 12:4 and Therefore, every unitary increase/decrease (e.g. 1 ha if home ranges are expressed in hectares, 1 km 2 if they are expressed in km 2 ) of A i determines a decrease/ increase in GOI, and a correspondent increase/decrease of GSI, equal to We applied GOI to the individual and the colonyspecific home ranges, in order to estimate withincolony (GOI W ) and between-colony (GOI B ) overlaps respectively. As suggested by the "diplomacy" hypothesis (Grémillet et al. 2004), spatial segregation among nearby colonies may mitigate intraspecific competition for resources. In order to test this hypothesis, we used a randomization procedure to determine if GOI B was greater than expected by chance. Under the null hypothesis of no spatial segregation between the two colonies, GOI B should not be significantly different from the size of the overlap if the GPS points of each colony were randomly and independently assigned. As Lesser Kestrels are central place foragers, distance is highly relevant and thus we could not assume they were free to visit all locations within the study area. Thus, we generated our null expectation by using a rotation with a random angle of the observed GPS points (by anchoring points to the coordinates of the correspondent urban colony), therefore randomly positioning the GPS points of each colony while keeping distances from the colony equal (Ferrarini et al. 2018a, b). In order to apply a rotational resampling of the colony data, all the data from each individual were randomly rotated around its nest location, independently from the other individuals. Mathematically, we used the standard algorithm (Eq. 5 and Eq. 6) for rotating points around a centre of rotation: where (x 0 , y 0 ) is the point to be rotated, (x c , y c ) are the coordinates of the nest location, θ is the angle of rotation (positive counterclockwise), (x 1 , y 1 ) are the coordinates of point after rotation. ( We then computed the randomly created home ranges (HR Rand ) for each colony, and overlaps (GOI Rand hereafter) between the two colonies. We repeated our randomizations 9999 times. The P-value for each colony was determined by the proportion of randomly created overlaps GOI Rand that were smaller than the observed overlap GOI B .
The five Lesser Kestrels of Santeramo in Colle scored an average home range size equal to 6132.58 ha (± 4026.29 std. dev.). The smallest and largest home ranges were  (Table 1). Pairwise percent overlaps (Table 3) ranged from 20.42% to 100%, with an average value equal to 58.09% (± 27.09 std. dev.). Individual home ranges were mostly nested within the home range of the second largest home range (individual M24 S), except for individual M4 S (Fig. 3). GOI W was equal to 81.38% (i.e., A i = 30,662.86 ha; max(A i ) = 11,147.12 ha; A i = 14,779.66 ha), thus GSI W was equal to = 18.62% (Fig. 3).
Colony-specific home ranges were 17,652.74 ha at Cassano and 13,228.31 ha at Santeramo, respectively. Between-colony overlap was 2529.46 ha. GOI B was equal to 19.12% (i.e., HR i = 30,881.05 ha; max(A i ) = 17,652.74 ha; A i = 28,351.58 ha; Fig. 4), thus GSI B was 80.88%. The proportion of randomly created overlaps GOI Rand that were smaller than the observed  Table 1 for the GPS ID of the individuals. The red squares represent the towns of Cassano and Santeramo where the Lesser Kestrels have their nests. GOI and GSI stand for general overlap index and general segregation index, respectively overlap GOI B was 4.36% (436 randomizations out of 9999; P = 0.0436), therefore the null hypothesis of no spatial segregation was rejected (P < 0.05).

Discussion
In this study, we have proposed a non-pairwise metric of home range overlap/segregation, and have applied it to two neighboring Lesser Kestrel's colonies in Southern Italy.
Our overlap index follows a simple idea: given n home ranges, it is always possible to calculate the extent of two spatial configurations, perfect segregation and perfect overlap. In the former case (Fig. 2a), the extent covered by the home ranges is simply the sum of their areas, and in the latter case (Fig. 2b) it is equal to the area of the largest home range. Our index simply measures the distance of the observed overlaps from these two extremes. In doing so, our overlap index does not require calculating pairwise overlaps among individual home ranges. The uniquely non-pairwise nature of the metric leads to two interesting properties: first, it is computationally fast as it just requires the union of home range polygons (Fig. 2d) to be calculated within GIS; second, the overlap score provided by GOI is semantically different from the overlap scores provided by pairwise overlap indices. In fact, GOI provides an estimate of how nested different home ranges areas are. This explains why in both colonies GOI did not duplicate the information provided by pairwise overlap measures (Tables 2, 3), and not even some statistical properties (e.g. mean or median) of such pairwise measures. In fact, mean pairwise overlap was 65.51% ± 27.64 (mean ± std. dev.) at Cassano, and it was 58.09% ± 27.09 at Santeramo. Therefore, GOI was outside the mean ± std. dev. interval at Cassano, and almost outside the right tail of the same interval at Santeramo. In addition, GOI was much easier to interpret in comparison to the 5 × 5 pairwise overlap matrices (Tables 2, 3).
Our overlap index also has several other desirable properties: (1) GOI can be applied to an arbitrarily large number of home ranges (i.e., n ≥ 2) belonging to individuals, populations (colonies) or species; (2) whatever the number of home ranges under study, GOI returns a single overlap measure; (3) in the case of perfectly disjoint home ranges, GOI is equal to 0; (4) in the case of perfectly nested (overlapping) home ranges, GOI is equal to 100; (5) in any other case, GOI returns a value between 0 and 100; (6) GOI varies linearly between 0 and 100, independently of the number of home ranges under study. In fact, Eqs. 1-4 ensure that GOI and GSI and their rates of change are independent of (a) the number of observations and (b) the initial value assumed by A i , but only depend on the geometric and positional properties of the home ranges. The linear nature of these metrics also ensures that small/big changes to the home range overlaps proportionally determine small/big changes to GOI (Additional file 2: Text S2).
In this study, we have applied GOI to 2D home ranges, however our overlap index can be readily applied to 3D home ranges as well (Tracey et al. 2014;Ferrarini et al. 2018b). In the case of volumetric home ranges, the 2D home range size should be simply replaced by 3D estimation, but GOI (and also GSI) would maintain the same properties described above. We estimated home ranges through the minimum convex polygons algorithm, however the application of GOI (and also GSI) is successive, and thus independent, of the type of algorithm (e.g. low convex hull; Getz et al. 2007) employed to assess birds' home ranges. Thus, both GOI and GSI can be applied to home range polygons derived from any type of home range estimator (Signer et al. 2015). We have applied GOI to a central-place forager, as this type of species presents elevated within-colony overlap thus making the use of an overlap index very appropriate. In the case of bird species with low overlaps, the alternative GSI index could be more suitable to readily assess the degree of home range segregation. The populations studied here showed elevated intracolony overlap and between-colony segregation. These results are in agreement with findings from the nearby colonies of Gravina in Puglia and Altamura (Ferrarini et al. 2018a, b), although in that case segregation was computed using a standard pairwise overlap index. During the chick rearing interval the demand for food is the highest, thus this might affect spatial segregation between the colonies. By foraging in spatially segregated areas, individuals from different colonies may avoid interference competition for food (Grémillet et al. 2004). It is therefore plausible that spatial segregation is relaxed in other periods when food demand is lower.

Conclusions
Our overlap index addresses the question of generalizing pairwise measures of home range overlap to a single measure of overlap within or across populations or species. It is not intended to replace the commonly used Fig. 4 Between-colony overlap. We first estimated colony-specific home ranges (i.e. home ranges calculated after pooling the locations of all individuals of each colony), then we calculated the between-colony overlap. The red squares represent the towns of Cassano and Santeramo where the Lesser Kestrels have their nests. GOI and GSI stand for general overlap index and general segregation index, respectively