Mapping areas of spatial-temporal overlap from wildlife tracking data
© Long et al. 2015
Received: 18 February 2015
Accepted: 20 September 2015
Published: 1 November 2015
The study of inter-individual interactions (often termed spatial-temporal interactions, or dynamic interactions) from remote tracking data has focused primarily on identifying the presence of such interactions. New datasets and methods offer opportunity to answer more nuanced questions, such as where on the landscape interactions occur. In this paper, we provide a new approach for mapping areas of spatial-temporal overlap in wildlife from remote tracking data. The method, termed the joint potential path area (jPPA) builds from the time-geographic movement model, originally proposed for studying human movement patterns.
The jPPA approach can be used to delineate sub-areas of the home range where inter-individual interaction was possible. Maps of jPPA regions can be integrated with existing geographic data to explore landscape conditions and habitat associated with spatial temporal-interactions in wildlife. We apply the jPPA approach to simulated biased correlated random walks to demonstrate the method under known conditions. The jPPA method is then applied to three dyads, consisting of fine resolution (15 minute sampling interval) GPS tracking data of white-tailed deer (Odocoileus virginianus) collected in Oklahoma, USA. Our results demonstrate the ability of the jPPA to identify and map jPPA sub-areas of the home range. We show how jPPA maps can be used to identify habitat differences (using percent tree canopy cover as a habitat indicator) between areas of spatial-temporal overlap and the overall home range in each of the three deer dyads.
The value of the jPPA approach within current wildlife habitat analysis workflows is highlighted along with its simple and straightforward implementation and interpretation. Given the current emphasis on remote tracking in wildlife movement and habitat research, new approaches capable of leveraging both the spatial and temporal information content contained within these data are warranted. We make code (in the statistical software R) for implementing the jPPA approach openly available for other researchers.
Through movement ecology, wildlife researchers continue to build a more detailed understanding of processes that shape wildlife movement patterns. The study of wildlife movement has been enhanced by advances in remote tracking (e.g., GPS, VHF, Argos) that continue to improve data quality, inference, and cost-effectiveness [1, 2]. Remote tracking offers unique opportunities for studying wildlife movement over broad spatial and temporal extents and at increasingly fine resolutions; these features help to address when, where, how, and why animals move . While defining home ranges, the area used for normal wildlife activities , remains a common approach to exploring wildlife space-use patterns, increasingly high resolution tracking data makes analyzing more complex and detailed spatial-temporal patterns in wildlife movement and behaviour possible.
Inter- and intra-species interactions can play a key role in the movement patterns of many wildlife species. Behaviour arising from inter- and intra-species interactions can occur from a number of different wildlife movement processes, such as the development of social networks , mating , and territoriality . For example, female white-tailed deer exhibit social behaviour in which they form matrilineal groupings for most of the year, except during parturition when females become solitary in preparation of, and immediately following, birth . The aforementioned processes are examples of dependent behaviour in wildlife movement, where the movements of one individual influence another, or inter-individual interaction (commonly termed dynamic or spatial-temporal interaction [9, 10]). Observed patterns of social or interactive behaviour have important implications in the management of spatially explicit wildlife processes; one such process is disease spread, which is related to the spatial and temporal patterns of contacts between individual animals, either directly or indirectly [11, 12].
Exogenous factors (e.g., environment, landscape) are also known to shape the movement patterns observed in wildlife . For example, wildlife movement processes and interactions are related to changing climatic and environmental conditions (e.g., range shifts, ), and increasing anthropogenic change (e.g., natural resource extraction, ). When the landscape changes, movement is often impacted and wildlife may be forced to shift their range and interact with other individuals or new species due to limited habitats and resource availability (e.g., water, ). Other landscape features can physically shape wildlife interactions by creating movement barriers and corridors [16–18]. In northern Alberta, Canada, Latham et al.  found that ungulates were utilizing anthropogenic cut-lines as movement corridors, which were in turn being used by predators to track prey, changing interacting patterns. The interplay between landscapes and wildlife interactions influence health and survival, and due to the emphasis on wildlife management through landscapes, have substantial implications [20, 21].
Studies examining interactions using wildlife tracking data typically aim to quantify one of two conceptually related yet unique joint movement processes: 1) spatial-only interaction, and 2) inter-individual interaction [9, 10]. Methods for exploring spatial-only interaction (often termed static interaction, ) involve quantifying the joint space use between two individuals, often through calculation of the overlap zone (OZ) of two home ranges. Extending the home range OZ is the volume of intersection measure, which delineates the joint probability of occurrence of two animals whose movements are characterized by two different utilization distributions . However, neither the OZ nor the volume of intersection incorporates timing of joint space use, or likelihood of wildlife encounters. While spatial-only interaction measures highlight areas utilized by multiple wildlife ; high temporal resolution tracking data enables the possibility of quantifying the likelihood of wildlife encounters [24, 25].
Inter-individual interactions can be analyzed using a suite of indices that test for the presence of interactions in tracking data (reviewed by ). These methods define inter-individual interaction associated with contacts measured using spatial (d c ) and temporal (t c ) thresholds to identify when location fixes co-occur in space and time. Contact based measures of inter-individual interaction test the observed number of contacts against an expectation or null model in order to identify attraction (higher than expected contacts) or avoidance (lower than expected contacts) behaviour [26, 27]. Other methods exist to study altogether different aspects of inter-individual interaction patterns in wildlife tracking data such as coordinated movement [28, 29], group-dynamics , and flocking or herding [31, 32].
A major limitation of currently available indices of inter-individual interaction is that they do not facilitate a spatially explicit measure of where interactive behaviour occurs on the landscape. Wildlife are known to select areas within their home ranges unevenly  and we expect interactive behaviour to also exhibit different spatial patterns. In this paper, we propose a new method for mapping areas of spatial-temporal overlap from wildlife tracking data. Mapping areas of spatial-temporal overlap will provide new avenues for research aimed at studying the linkages between interactive behaviour and environmental factors. The new approach draws upon previous work using time geography to estimate wildlife home ranges , and extends methods used to study interactions in human movement studies  to the study of wildlife movement. First, we introduce the theory of time geography and its current application to wildlife studies and then describe the potential path area (PPA) approach to delineating the home range and how we extend this to compute a new measure of spatial-temporal overlap. Following this, we demonstrate the new method with simulated and empirical data. We finish with discussion of our findings and future opportunities for studying inter-individual interactions from wildlife tracking data.
Background – time geography
The space-time prisms of two individuals can be intersected in order compute the joint accessibility space – termed the social interaction space . The social interaction space represents a measure of the areas, in space and time, where direct inter-individual interaction is possible (i.e., contact between the two individuals are only possible within the spatial-temporal boundaries of the social interaction space). The social interaction space can be projected onto geographical space in order to create a map of spatial-temporal overlap – termed the joint potential path area (jPPA – Fig. 1b). Note that the jPPA is fundamentally different from a measure of home range, or spatial overlap, as it explicitly shows only those regions where two individuals have the potential to ‘meet’ in space and time.
Calculating the PPA
In application, calculation of the PPA is straightforward because the PPA is a perfect ellipse, the parameters of which can be calculated directly from the fixes and the mobility parameter v max (see ).
Calculating the jPPA
By definition, the jPPA will always be a sub-region of the intersection of the two PPA ellipses (see Fig. 1b) because the intersection of two ellipses generated by the PPA is the measure of the spatial-only overlap, while the jPPA represents spatial-temporal overlap. In order to compute the jPPA for a larger tracking dataset, we must recursively compute G τ A,B for all overlapping τ within the temporal period when animals A and B were simultaneously tracked.
A simulation study was used to examine the jPPA method and contrast it with the commonly employed spatial overlap measure in order to evaluate the ability of each for identifying and mapping interactive movement behaviour. Correlated random walks (CRW – [38, 39]) and biased correlated random walks (BCRW – [40, 41]) were used to generate simulated tracking data where inter-individual interactions would be expected and unexpected (null case). Four scenarios were implemented representing different types of inter-individual interactions commonly encountered in wildlife systems: i) no interaction, ii) grouping, iii) leading/following, and iv) joint resource use. In the first scenario, no interaction is simulated via two independent CRW, the second originating within the minimum convex polygon of the first. In the second scenario, grouping, a CRW is used to simulate the dynamic location of a group centroid, and the two individuals movements are biased towards this location . In the third scenario, leading, the movements of the second individual are biased towards the current position of the first . In the fourth scenario, joint resource use, the movements of the two individuals are biased towards a patch collocated between two home ranges . We allowed the simulated animals to switch between CRW and BCRW with some fixed probability at each step, to emulate a behavioural switching between interactive and non-interactive phases . Simulations where the spatial overlap in home ranges was < 10 % of the combined home range area were discarded corresponding to the idea that spatial overlap is a pre-cursor for spatial-temporal interaction . Similarly, in the scenarios 2, 3, and 4 we discarded all simulations where < 10 % of the time was spent in the interactive phase to ensure these scenarios demonstrated the expectation of interaction. We ran the simulation process for each scenario until 1000 simulations were achieved. More details on the simulation procedure, and accompanying R code, can be obtained from the Additional file 1.
For each simulated animal, we first computed the PPA estimate of home range . We then computed a measure of spatial-only interaction, the overlap zone defined as OZ = A∩B, where A and B are the PPA home ranges of the two simulated individuals. Then we computed the jPPA. To facilitate straightforward comparisons, the areas of the OZ and jPPA were normalized by the total joint home range area (A⋃B), such that each measure ranged from 0 to 1; 0 indicating no spatial overlap, and 1 complete spatial agreement. Evaluating the simulations was done by considering the 1000 runs from the first scenario (no interaction) as a test distribution representing where interaction is unexpected. Then we would expect that the values of the OZ and jPPA in the three scenarios where interaction was expected to lie in the outer tails of this distribution. In such a randomization test, the one-sided test statistic is computed as the probability p = (n e + 1)/n of getting a value equal to or more extreme in the null distribution in comparison to the observed value, where n e is the count of these extreme values. Using a critical value of α = 0.05, we examine the ability of the jPPA against the naive OZ statistic for identifying expected interaction behaviour. We contrast the performance jPPA against the OZ statistic in comparison to the OZ. Further, for the three simulated scenarios where BCRW were used to simulate known interactive behavior, we further contrast the jPPA with the proportion of time where the two animals were in the BCRW interactive phase (termed pInt). The value for pInt can be considered an indicator of the known or true level of interaction in BCRW simulations.
Empirical data: GPS tracking of white-tailed deer
We captured 38 white-tailed deer (Odocoileus virginianus) from a study area in south-central Oklahoma, USA from 1998 to 2004. The study site was 1,214 ha, and was surrounded by a 15-strand, high-tensile electric fence (2.5 m tall), thus partially restricting movement across property boundaries . Vegetation was consistent with that of the Cross Timbers and Prairies ecoregion . Deer were captured during the winter months using modified drop-net systems  and fitted with GPS collars (ATS G2000 remote-release collars; Advanced Telemetry Systems, Inc., Isanti, MN) programmed to collect fixes at a 15 min sampling interval. Data were successfully retrieved from 32 of 38 GPS collars. All capture, handling, and marking procedures were consistent with the guidelines of the American Society of Mammalogists  and were approved by permit from the Oklahoma Department of Wildlife Conservation.
Sex, age, and the tracking period of six white-tailed deer that were part of the empirical analysis of jPPA; the 6 deer represented 3 unique dyads
21-Nov - 08-Feb
21-Nov - 06-Dec
23-Nov - 17-Jan
14-Dec - 17-Feb
15-Feb - 28-Mar
15-Feb - 21-Mar
For each individual deer, we computed the PPA estimate of individual home range . From the PPA home range estimates, we computed the spatial overlap zone (OZ), similar to previous home range overlap analysis [23, 47] as a measure of spatial-only interaction. Finally, we computed the jPPA for each dyad. Calculating the PPA and jPPA home range areas requires estimating the v max parameter, which represents an upper bound on mobility. We expect deer to show higher levels of mobility at dawn and dusk than during the night and day; therefore, we estimated the v max parameter for the PPA and jPPA dynamically  for four time periods throughout the day: dawn (05:00 – 09:00), day (09:00 – 17:00), dusk (17:00 – 19:00), and night (19:00 – 05:00) following temporal intervals from . In all cases, we used the van der Watt  method (with k = 5) for estimating v max from the tracking data following [34, 48].
We investigated differences in vegetation within each of the individual home ranges, the OZ, and the jPPA for each of the three deer dyads. We chose a single metric – percent canopy cover – as a representative indicator of habitat. Vegetative cover is an important habitat component for white-tailed deer because it helps regulate the local thermal environment and provides concealment against predators [51, 52]. Percent canopy cover data were obtained from the US National Land Cover Database (NLCD, ). The NLCD percent canopy cover data are derived from Landsat satellite imagery and are represented at a spatial resolution of 30 m. For each of the three deer dyads, we calculated the mean and standard deviation of percent canopy cover (i.e., the mean of all pixels) for the pixels associated with the original telemetry points, the area within each individual’s PPA home range, the OZ, and the jPPA. Further, we examined the distribution of the percent canopy cover values associated with the individual points and within each of the home range and joint areas using overlaid density plots (with a bandwidth of 10) to explore variability in use of canopy cover across the range of potential values (0 – 100 %).
Performance of the OZ and jPPA method at correctly identifying the presence of interaction from the three simulation scenarios: group behavior, leading/following, and joint resource use (1000 simulations were run for each scenario)
Sig.a interaction identified by only one method
No sig.a interaction identified
Joint resource use
Empirical data: GPS tracking of white-tailed deer
The time geography movement model
The time geographic approach to animal space-use analysis is aided by its simple implementation and interpretation. Unlike current methods for studying inter-individual interaction, the jPPA method focuses on mapping inter-individual interactions, defined as the area of spatial-temporal overlap. Output polygons from the jPPA are easily interpreted and can be readily integrated into existing workflows common to many wildlife movement studies (e.g., home range and habitat analysis). The jPPA method does not facilitate a statistical test for the presence of significant inter-individual interaction, thus, the jPPA compliments the existing suite of methods for studying spatial and inter-individual interactions currently available to wildlife researchers.
Conceptually, the jPPA extends the time-geographic approach for quantifying space use in wildlife tracking data [34, 54], which explicitly considers the temporal ordering of fixes and a parameter of mobility in order to quantify what is termed the accessibility space of an individual. The jPPA is then defined as the joint accessibility space of two individuals, and has been successfully applied to numerous studies of human mobility and transportation patterns [35, 55]. The jPPA should therefore be interpreted as the spatial region delineating the areas that are jointly accessible, in space and time, by two animals. Herein, we introduced the jPPA with k = 2 individuals (termed dyads), but the same calculation can be extended to k > 2 individuals to delineate the joint accessibility space of a larger group of animals. For example, local resources may be simultaneously utilized between multiple species or wildlife assemblages. Where tracking data of multiple individuals are recorded simultaneously, the jPPA method can be used to locate these areas across the landscape.
Space-time prisms (which are the basis for the jPPA) are straightforward to construct from tracking data (e.g., fix locations and times), along with a parameter describing the mobility of the animal (termed v max ), which is interpreted as maximum travelling velocity. It is best to consider v max as a function of the sampling interval; for example, v max will change depending on whether tracking data were collected with a 10 min or 2 h sampling interval. The estimation of v max can be based on expert knowledge of animal biology or using statistical procedures based on the tracking data . When estimating v max from the tracking data it is useful to consider the distribution of the observed segment velocities (v i ) given by v i = d i / t i , where d i is the distance and t i the time between consecutive fixes. Based on the full distribution of the v i , statistical estimation procedures used for estimating the upper-bound of a distribution [50, 56] can be applied to estimate v max . Mobility levels change over space and time, thus the v max parameter can be modeled dynamically to better reflect changes in individual behaviour or travel mode .
The jPPA is most appropriately applied where tracking data are collected with a moderately high resolution and regular fix interval (i.e., such as the white-tailed deer data here using a 15 min interval). This is because the jPPA method explicitly considers the sequential ordering and the time duration between fixes in its derivation (i.e., equations (1), (2)). When the sampling interval is extremely high (e.g., some species are now being tracked at ≤ 1 min intervals, ) the potential for movement between known fixes is limited, and thus the PPA will be small. Similarly, when the sampling interval is extremely coarse, the potential for movement between known fixes is substantial, and thus the PPA will be large. In the first case, the jPPA is likely to underestimate potential interaction because of a mis-match between the data resolution and the functional scale of interaction (e.g., defined by sight, smell, sound etc.). In the second case, the jPPA will over-estimate interaction because of the immense possibility of potential movement when fixes are taken infrequently (e.g., daily or less). Much like the PPA estimate of home range [34, 48], jPPA should be used with caution with much coarser sampling resolutions because of the uncertainty associated with the delineation of the potential path area when the time duration between fixes is long. However, with high-resolution tracking data, the jPPA may be computed at multiple resolutions via re-sampling of the original tracking data. Plotting jPPA results against different tracking resolutions may provide insights into the functional scale of movement interaction, considering both the mobility of the animal and the potential spatial range associated with interactions. For example, questions may begin by discovering at what sampling interval the jPPA becomes significant, and then begin to consider what movement processes might correspond to this interval (e.g., fine vs. coarse scale movement).
Studying inter-individual interactions from wildlife tracking data
The jPPA provides an alternative approach to studying inter-individual interactions from wildlife tracking data, one that focuses specifically on spatially explicit mapping of inter-individual interaction. Unlike previously developed indices of inter-individual interaction, the jPPA does not facilitate a formal statistical test of the presence (or absence) of inter-individual interaction. Statistical tests for inter-individual interaction from wildlife tracking data can be problematic, owing to the issue of generating appropriate null distributions from which to test against [26, 27, 58]. Rather, the jPPA is able to detect, and more importantly map, infrequent and/or random inter-individual interaction areas (e.g., chance encounters) across the landscape. Such infrequent or random interactions typically go undetected when using formal statistical tests because they do not constitute a statistically significant interaction . Because of their importance in shaping biological processes, such as the spread of disease and predator–prey dynamics, methods capable of identifying and mapping random or unexpected encounters, such as jPPA, offer new potential for studying infrequent or random interactions by wildlife using tracking data.
The jPPA approach can be compared to existing measures of static interaction; defined as the spatial area used jointly by two (or more) animals [9, 10]. We compared the jPPA to the most common measure of static interaction, the overlap zone, defined as the spatial intersection of two home ranges. When the PPA home range estimate is employed, the jPPA will be a spatial sub-region of the OZ and represents those areas where animals could have potentially ‘meet’ in both space and time. Recent research has suggested that the utilization distribution, which represent space use as an uneven probability surface [59, 60], may be a more useful spatial measurement as animals typically use different areas of their home ranges unevenly in space and time . The volume of intersection [22, 33] then represents the analogous spatial interaction measure for utilization distributions; delineating a map showing the joint space-use probabilities. However, the volume of intersection of two (or more) utilization distributions does not explicitly consider time (i.e., simultaneous joint space use). Thus, the jPPA can be considered the spatial-temporal extension of the home range overlap zone and represents a discrete map (i.e., a polygon) of where inter-individual interaction was possible. Further extending the jPPA using some probabilistic models [62, 63] could facilitate probabilistic statements similar to those from the volume of intersection measure for studying inter-individual interaction probabilities.
New types of tracking sensors (i.e., proximity loggers) are capable of directly estimating inter-individual interactions (i.e., contacts) by sensing when, and for how long, two tagged animals are within a defined distance threshold . However, proximity loggers are not designed to record where contacts occur, and the actual distance between two animals is not known exactly, which has led to development of new devices that will combine location-aware technology (e.g., GPS) with proximity loggers . The temporal resolution of proximity logger data (often programmed to record continuously, ) will typically far exceed that of location data (e.g., from GPS), predominantly due to battery limitations. When proximity logger data are combined with GPS, proximity data can be incorporated into the calculation of the jPPA by computing the jPPA for only those times where the proximity logger identifies a contact to have occurred (i.e., τ in equation (6)). Such an extension will allow researchers to more precisely delineate regions where contacts occur when using the jPPA and incorporate spatially aware technology (e.g., GPS) to study the landscape context (habitat, topography, connectivity) associated with interactive behavior.
Insights from the white-tailed deer examples
One of the most important outcomes of the development of new algorithms for studying inter-individual interaction is the application of these methods to discover previously undocumented biological insights of animal movement behavior. The results from the white-tailed deer example depict only three selected cases, selectively chosen to demonstrate different scenarios of interaction, and should not be used to make more general inferences about deer movement ecology. However, we can make some interesting conclusions about the specific individuals in our study to provide context for future analysis. Specifically, dyad 1 consisted of two adult males that were monitored concurrently over an approximate 2-week period, which fell within the breeding season of deer in this study area . The results from our jPPA analysis for Dyad 1 suggest that any interaction may be related to brief and sporadic encounters. Biologically, this scenario is realistic considering that male deer do not socialize during the breeding season, but may combat with each other only temporarily, or may use similar areas when an estrous female is encountered. Dyad 2 consists of an adult male and adult female during the winter months suggesting that the interactive behavior observed in this dyad may be related to mating behavior (i.e., courting), which typically occurs on the Oklahoma, USA study area from 4 November to 24 December, with most breeding occurring from 18 November to 2 December . In this case, jPPA analysis revealed that the potential courting behavior occurred in a specific wooded region of the study area (i.e., localized spatial extent), containing higher levels of percent canopy cover in comparison with other areas of the PPA ranges and OZ. The pair of males in dyad 3 are young adults that tend to form bachelor groups during spring and summer . With dyad 3, the concurrent tracking period occurred from 15 February to 21 March, and the large jPPA area throughout their individual movement ranges may be an indication of the initiation of a bachelor pair (bachelor groups could be identified when analyzing k > 2 individuals). Habitat analysis revealed very little difference between the percent canopy cover within the PPA ranges, the OZ, and the jPPA, which would be expected considering the strong level of inter-individual interaction observed, and the general habitat needs of male deer following the breeding season, and during the antler shedding and growing cycle.
In general, it was also found that percent canopy cover estimates obtained using the telemetry points, differed from those from polygon-based home ranges. The general trend observed was that mean percent canopy cover was greater when the estimates were obtained from the telemetry points. One insight from this is that deer may be preferentially using forest edges , which when PPA home ranges are computed, the home range will include some open areas nearby as well. Then, it is interesting that the jPPA percent canopy cover values from dyad 2 were associated with higher percent canopy cover values than from the telemetry points, suggesting that this interaction phase occurred within a core forest area, as opposed to edge habitat. The problem of obtaining different results from habitat analysis when using telemetry points as the spatial unit vs. home range polygons as the spatial unit is well documented  and remains an ongoing challenge to the study of habitat composition from wildlife tracking data.
We hypothesize that using the jPPA to examine habitat use can reveal social behaviour and interactions of animals that may be driven by underlying landscape features (e.g., use of a shared or rare resource such as when animals visit watering holes, which could result in greater contact and interaction) or vice versa (use of landscape features is driven by social behaviour or interactions; e.g., during mating, animals may seek certain habitat types). Specifically, in dyad 3 from this study, we identified how the use of dense canopy cover within the jPPA would indicate generally important resources to male deer during winter (i.e., post-rut), which require greater energy and nutrition to recover from the stresses of the rut and for development of antlers in the spring [70, 71]. Knowing where and when habitat types and resources are used can facilitate management targeted at social groups (e.g., herds), rare resources where animals may be in contact (e.g., water sources), or during important life history phases (e.g., during mating).
Further applications of jPPA analysis
We envision that the jPPA approach will be attractive to many wildlife ecologists because comparing jPPA polygons with ancillary geographic datasets is straightforward within a geographic information system (GIS) and accommodates similar spatial analyses as commonly applied to wildlife home ranges. The jPPA offers a new approach aimed specifically at mapping areas across the landscape where inter-individual interaction occurs. Wildlife researchers can use the jPPA polygons to test spatial hypotheses related to how interactive movement behaviour relates to underlying environmental variables. For example, researchers can estimate habitat composition associated with interactive behaviour using jPPA polygons in a similar manner to how composition is analyzed within home ranges, core areas, and utilization distributions . We demonstrate such a process using a simple example where a single variable (percent canopy cover) is mapped across our study area. However, there is clear potential to use jPPA alongside multiple mapped covariates (e.g., landcover, topography, or weather data) to develop more sophisticated spatial analyses of resource selection to uncover greater biological insight into interactive behaviour. When performing habitat analysis from tracking data, it is important to consider how the organism perceives the environment, relating not only to the composition of habitat but also spatial pattern and scale [72–74]. As an alternative to pixel-based raster data, landscape features can be represented as distinct habitat ‘patches’, using a polygon data format, linear features, like roads, represented as lines, and other features, such as oil and gas well-sites, represented as points. In many species it will be interesting to examine how interactive behaviour is associated with these other types of features on the landscape (e.g., ). Output from jPPA analysis can be easily integrated with other geographic datasets within a GIS allowing more sophisticated spatial analysis of interactive movement behaviour.
The spatial patterns associated with the jPPA areas can also be used to understand biological processes associated with different interactive behaviours. For example, the configurational properties of jPPA regions may provide important information on the type of social interaction occurring. When jPPA regions are small and patchy relative to the home range and OZ it may be evidence of random encounters occurring across the landscape. Conversely, when the jPPA region is small but contiguous, the interaction observed may be an indication spatially localized resources being used simultaneously. When the jPPA covers a large proportion of the overlap zone, it is evidence of a sustained high-level of interactive behaviour easily corroborated by one of the many interaction statistics currently available [25, 26]; as evidenced by deer dyad 3 in our analysis. Further, more detailed analysis offers the potential to quantify which properties of the jPPA are associated with different social behaviour patterns. Polygon-shape indices, widely applied in the study of landscape patterns , offer potentially valuable metrics which could be included into the analysis of jPPA polygon regions in order to quantify, for example, compact or patchy shapes.
The study of inter-individual interactions in wildlife is important to many population-level processes and is thus of special interest to wildlife managers. Here we have demonstrated a new approach for mapping areas of inter-individual interaction from wildlife tracking data. The new approach extends a previously developed home range estimator  in-order to delineate areas of joint accessibility (termed the jPPA) between two (or more) simultaneously tracked animals. Maps of the jPPA provide researchers with a new spatial unit from which habitat analysis can be easily conducted and directly related to inter-individual interactions. We also demonstrate the application of the jPPA approach in habitat analysis by exploring the percent canopy cover with three deer dyads tracked using GPS tracking data. Spatial patterns of the jPPA, especially in relation to the OZ and home range areas, can be useful indication of the type of inter-individual interaction occurring in a dyad (e.g., random encounters vs. joint spatial use of local resources). Extending current indices that test for the presence or absence of interaction behaviour and new developments capable of mapping inter-individual interactions in space and time (such as the jPPA) are essential for studying the complex and infrequent social encounters of wildlife using remote tracking data. Finally, in order to assist other researchers wishing to utilize the jPPA in their own research, we make openly available code for computing the jPPA in the statistical software R  (for access to R code and tools please see: http://jedalong.github.io/wildlifeTG/).
The collection of the deer data used in this study was supported by funding from The Samuel Roberts Noble Foundation, Inc.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Cagnacci F, Boitani L, Powell RA, Boyce MS. Animal ecology meets GPS-based radiotelemetry: a perfect storm of opportunities and challenges. Philos Trans R Soc Lond B Biol Sci. 2010;365:2157–62.View ArticleGoogle Scholar
- Tomkiewicz SM, Fuller MR, Kie JG, Bates KK. Global positioning system and associated technologies in animal behaviour and ecological research. Philos Trans R Soc Lond B Biol Sci. 2010;365:2163–76.View ArticleGoogle Scholar
- Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D, et al. A movement ecology paradigm for unifying organismal movement research. Proc Natl Acad Sci. 2008;105:19052–9.View ArticleGoogle Scholar
- Burt WH. Territoriality and home range concepts as applied to mammals. J Mammal. 1943;24:346–52.View ArticleGoogle Scholar
- Prange S, Gehrt SD, Hauver S. Frequency and duration of contacts between free-ranging raccoons: uncovering a hidden social system. J Mammal. 2011;92:1331–42.View ArticleGoogle Scholar
- Stenhouse GB, Boulanger J, Lee J, Graham K, Duval J, Cranston J. Grizzly bear associations along the eastern slopes of Alberta. Ursus. 2005;16:31–40.View ArticleGoogle Scholar
- Giuggioli L, Kenkre VM. Consequences of animal interactions on their dynamics: emergence of home ranges and territoriality. Mov Ecol. 2014;2:20.View ArticleGoogle Scholar
- Bertrand MR, DeNicola AJ, Beissinger SR, Swihart RK. Effects of parturition on home ranges and social affiliations of female white-tailed deer. J Wildl Manage. 1996;60:899–909.View ArticleGoogle Scholar
- Macdonald DW, Ball FG, Hough NG. The evaluation of home range size and configuration using radio tracking data. In: Amlaner CJ, MacDonald DW, editors. A Handbook on Biotelemetry and Radio Tracking. Oxford: Pergamon Press; 1980. p. 405–24.View ArticleGoogle Scholar
- Kernohan BJ, Gitzen RA, Millspaugh JJ. Analysis of Animal Space use and Movements. In: Millspaugh JJ, Marzluff JM, editors. Radio Tracking and Animal Populations. New York: Academic; 2001. p. 125–66.View ArticleGoogle Scholar
- Ramsey D, Spencer N, Caley P. The effects of reducing population density on contact rates between brushtail possums: implications for transmission of bovine tuberculosis. J Appl Ecol. 2002;39:806–18.View ArticleGoogle Scholar
- Tosa MI, Schauber EM, Nielsen CK. Familiarity breeds contempt: combining proximity loggers and GPS reveals female white-tailed deer avoiding close contact with neighbors. J Wildl Dis. 2015;51.
- Parmesan C, Yohe G. A globally coherent fingerprint of climate change impacts across natural systems. Nature. 2003;421:37–42.View ArticleGoogle Scholar
- Sawyer H, Nielson R, Lindzey F, McDonald LL. Winter habitat selection of mule deer before and during development of a natural gas field. J Wiildlife Manag. 2006;70:396–403.View ArticleGoogle Scholar
- Rosenstock SS, Ballard WB, Devos JC. Viewpoint: benefits and impacts of wildlife water developments. J Range Manag. 1999;52:302–11.View ArticleGoogle Scholar
- Sawyer H, Kauffman M, Nielson R, Horne JS. Identifying and prioritizing ungulate migration routes for landscape-level conservation. Ecol Appl. 2009;19:2016–25.View ArticleGoogle Scholar
- Laberee K, Nelson TA, Stewart BP, McKay T, Stenhouse GB. Oil and gas infrastructure and the spatial pattern of grizzly bear habitat selection in Alberta, Canada. Can Geogr / Le Géographe Can. 2014;58:79–94.View ArticleGoogle Scholar
- Stewart BP, Nelson TA, Wulder MA, Nielsen SE, Stenhouse G, Nelson TA, et al. Impact of disturbance characteristics and age on grizzly bear habitat selection. Appl Geogr. 2012;34:614–25.View ArticleGoogle Scholar
- Latham A, Latham M, Boyce M, Boutin S. Movement responses by wolves to industrial linear features and their effect on woodland caribou in northeastern Alberta. Ecol Appl. 2011;21:2854–65.View ArticleGoogle Scholar
- Dzialak MR, Webb SL, Harju SM, Winstead JB, Wondzell JJ, Mudd JP, et al. The spatial pattern of demographic performance as a component of sustainable landscape management and planning. Landsc Ecol. 2011;26:775–90.View ArticleGoogle Scholar
- Dzialak MR, Harju SM, Osborn RG, Wondzell JJ, Hayden-Wing LD, Winstead JB, et al. Prioritizing conservation of ungulate calving resources in multiple-use landscapes. PLoS One. 2011;6, e14597.View ArticleGoogle Scholar
- Millspaugh JJ, Gitzen RA, Kernohan BJ, Larson M, Clay CL, Larson MA, et al. Comparability of three analytical techniques to assess joint space use. Wildl Soc Bull. 2004;32:148–57.View ArticleGoogle Scholar
- Robert K, Garant D, Pelletier F. Keep in touch: does spatial overlap correlate with contact rate frequency? J Wildl Manage. 2012;76:6.View ArticleGoogle Scholar
- Schauber E, Storm D, Nielson C. Effects of joint space use and group membership on contact rates among white tailed deer. J Wildl Manage. 2007;71:155–63.View ArticleGoogle Scholar
- Long JA, Nelson TA, Webb SL, Gee KL. A critical examination of indices of dynamic interaction for wildlife telemetry studies. J Anim Ecol. 2014;83:1216–33.View ArticleGoogle Scholar
- Miller JA. Using spatially explicit simulated data to analyze animal interactions: a case study with brown hyenas in northern Botswana. Trans GIS. 2012;16:271–91.View ArticleGoogle Scholar
- Benhamou S, Valeix M, Chamaillé-Jammes S, Macdonald DW, Loveridge AJ. Movement-based analysis of interactions in African lions. Anim Behav. 2014;90:171–80.View ArticleGoogle Scholar
- Shirabe T. Correlation Analysis of Discrete Motions. In: Raubal M, Miller HJ, Frank AU, Goodchild MF, editors. GIScience 2006, LNCS. 4197th ed. Berlin: Springer; 2006. p. 370–82.Google Scholar
- Long JA, Nelson TA. Measuring dynamic interaction in movement data. Trans GIS. 2013;17:62–77.View ArticleGoogle Scholar
- Langrock R, Hopcraft JGC, Blackwell PG, Goodall V, King R, Niu M, et al. Modelling group dynamic animal movement. Methods Ecol Evol. 2014;5:190–9.View ArticleGoogle Scholar
- Benkert M, Gudmundsson J, Hubner F, Wolle T. Reporting flock patterns. Comput Geom. 2008;41:111–25.View ArticleGoogle Scholar
- Perna A, Grégoire G, Mann RP. On the duality between interaction responses and mutual positions in flocking and schooling. Mov Ecol. 2014;2:22.View ArticleGoogle Scholar
- Fieberg J, Kochanny CO. Quantifying home-range overlap: the importance of the utilization distribution. J Wildl Manage. 2005;69:1346–59.View ArticleGoogle Scholar
- Long JA, Nelson TA. Time geography and wildlife home range delineation. J Wildl Manage. 2012;76:407–13.View ArticleGoogle Scholar
- Farber S, Neutens T, Miller HJ, Li X. The social interaction potential of metropolitan regions: a time-geographic measurement approach using joint accessibility. Ann Assoc Am Geogr. 2013;103:483–504.View ArticleGoogle Scholar
- Hägerstrand T. What about people in regional science? Pap Reg Sci Assoc. 1970;24:7–21.View ArticleGoogle Scholar
- Miller HJ. A measurement theory for time geography. Geogr Anal. 2005;37:17–45.View ArticleGoogle Scholar
- Kareiva PM, Shigesada N. Analyzing insect movement as a correlated random walk. Oecologia. 1983;56:234–8.View ArticleGoogle Scholar
- Codling E, Hill N, Pitchford J, Simpson S. Random walk models for the movement and recruitment of reef fish larvae. Mar Ecol Prog Ser. 2004;279:215–24.View ArticleGoogle Scholar
- Barton KA, Phillips BL, Morales JM, Travis JMJ. The evolution of an “intelligent” dispersal strategy: biased, correlated random walks in patchy landscapes. Oikos. 2009;118:309–19.View ArticleGoogle Scholar
- Fronhofer EA, Hovestadt T, Poethke H-J. From random walks to informed movement. Oikos. 2013;122:857–66.View ArticleGoogle Scholar
- Morales J, Haydon D, Frair J, Holsinger KE, Fryxell JM. Extracting more out of relocation data: building movement models as mixtures of random walks. Ecology. 2004;85:2436–45.View ArticleGoogle Scholar
- Webb SL, Gee KL, Demarais S, Strickland BK, DeYoung RW. Efficacy of a 15-strand high-tensile electric fence to control white-tailed deer movements. Wildl Biol Pract. 2009;5:45–57.View ArticleGoogle Scholar
- Gee KL, Porter MD, Demarais S, Bryant FC, Vreede GV. White-Tailed Deer: Their Foods and Management in the Cross Timbers. 2nd ed. Ardmore, Oklahoma: Samual Roberts Noble Foundation Publication; 1994.Google Scholar
- Gee KL, Holman JH, Demarais S. A man-Power Efficient Drop-net System for Capturing White-Tailed Deer. In: Abstracts of the 22nd Annual Southeast Deer Study Group Meeting. 1999. p. 31.Google Scholar
- Gannon WL, Sikes RS, The Animal Care and Use Committee of the American Society of Mammalogists. Guidelines of the american society of mammalogists for the use of wild mammals in research. J Mammal. 2007;88:809–23.View ArticleGoogle Scholar
- Webb SL, Hewitt DG, Hellickson MW. Scale of management for mature male white-tailed deer as influenced by home range and movements. J Wildl Manage. 2007;71:1507–12.View ArticleGoogle Scholar
- Long JA, Nelson TA. Home range and habitat analysis using dynamic time geography. J Wildl Manage. 2015;79:481–90.View ArticleGoogle Scholar
- Webb SL, Gee KL, Strickland BK, Demaris S, DeYoung RW, Demarais S. Measuring fine-scale white-tailed deer movements and environmental influences using GPS collars. Int J Ecol. 2010;2010(459610):12.Google Scholar
- Van der Watt P. A note on estimation bounds of random variables. Biometrika. 1980;97:712–4.View ArticleGoogle Scholar
- Pollock MT, Whittaker D, Demarais S, Zaiglin RE. Vegetation characteristics influencing site selection by male white-tailed deer in Texas. J Range Manag. 1994;47:235–9.View ArticleGoogle Scholar
- Wiemers DW, Fulbright TE, Wester DB, Ortega-S JA, Rasmussen GA, Hewitt DG, et al. Role of thermal environment in habitat selection by male white-tailed deer during summer in Texas, USA. Wildlife Biol. 2014;20:47–56.View ArticleGoogle Scholar
- Homer C, Dewitz J, Fry J, Coan M, Hossain N, Larson C, et al. Completion of the 2001 national land cover database for the counterminous United States. Photogramm Eng Remote Sens. 2007;73:337–41.Google Scholar
- Downs JA, Horner MW, Tucker AD. Time-geographic density estimation for home range analysis. Ann GIS. 2011;17:163–71.View ArticleGoogle Scholar
- Neutens T, Witlox F, de Weghe N, De Maeyer P, Van de Weghe N. Human interaction spaces under uncertainty. Transp Res Rec. 2007;2021:28–35.View ArticleGoogle Scholar
- Robson DS, Whitlock JH. Estimation of a truncation point. Biometrika. 1964;51:33–9.View ArticleGoogle Scholar
- Shamoun-Baranes J, van Loon EE, Purves RS, Speckmann B, Weiskopf D, Camphuysen CJ. Analysis and visualization of animal movement. Biol Lett. 2012;8:6–9.View ArticleGoogle Scholar
- Miller JA. Towards a Better Understanding of Dynamic Interaction Metrics for Wildlife: a Null Model Approach. Trans GIS 2015;19:342–361.
- Worton B. Kernel methods for estimating the utilization distribution in home-range studies. Ecology. 1989;70:164–8.View ArticleGoogle Scholar
- Van Winkle W. Comparison of several probabilistic home-range models. J Wildl Manage. 1975;39:118–23.View ArticleGoogle Scholar
- Keating KA, Cherry S. Modeling utilization distributions in space and time. Ecology. 2009;90:1971–80.View ArticleGoogle Scholar
- Buchin K, Sijben S, Willems EP, Arseneau TJM. Detecting Movement Patterns Using Brownian Bridges. In: ACM SIGSPATIAL. Redondo Beach, CA, USA: ACM Press; 2012. p. 119–28.View ArticleGoogle Scholar
- Downs JA, Lamb D, Hyzer G, Loraamm R, Smith ZJ, O’Neal BM. Quantifying spatio-temporal interactions of animals using probabilistic space–time prisms. Appl Geogr. 2014;55:1–8.View ArticleGoogle Scholar
- Prange S, Jordan T, Hunter C, Gehrt SD. New radiocollars for the detextion of proximity among individuals. Wildl Soc Bull. 2006;34:1333–44.View ArticleGoogle Scholar
- Davis MJ, Thokala S, Xing X, Hobbs NT, Miller MW, Han R, et al. Testing the functionality and contact error of a GPS-based wildlife tracking network. Wildl Soc Bull. 2013;7.
- Webb SL, Riffell SK, Gee KL, Demaris S, Demarais S. Using fractal analyses to characterize movement paths of white-tailed deer and response to spatial scale. J Mammal. 2009;90:1210–7.View ArticleGoogle Scholar
- Hirth DH. Social behavior of white-tailed deer in relation to habitat. Wildl Monogr. 1977;53:3–55.Google Scholar
- Alverson WS, Waller DM, Solheim SL. Forests Too deer: edge effects in northern Wisconsin. Conserv Biol. 1988;2:348–58.View ArticleGoogle Scholar
- Aebischer NJ, Robertson PA, Kenward RE. Compositional analysis of habitat use from animal radio-tracking data. Ecology. 1993;74:1313–25.View ArticleGoogle Scholar
- Hewitt DG. Nutrition. In: Hewitt DG, editor. Biology and Management of White-tailed Deer. Boca Raton, FL: CRC Press; 2011. p. 75–105.Google Scholar
- Demarais S, Strickland B. Antlers. In: Hewitt D, editor. Biology and Management of White-tailed Deer. Boca Raton, FL: CRC Press; 2011. p. 107–45.Google Scholar
- McIntyre NE, Wiens JA. How does habitat patch size affect animal movement? An experiment with darkling beetles. Ecology. 1999;80:2261–70.View ArticleGoogle Scholar
- Kotliar NB, Wiens JA. Multiple scales of patchiness and patch structure: a hierarchical framework for the study of heterogeneity. Oikos. 1990;59:253–60.View ArticleGoogle Scholar
- With KA. Using fractal analysis to assess how species perceive landscape structure. Landsc Ecol. 1994;9:25–36.View ArticleGoogle Scholar
- McGarigal K, Marks B. FRAGSTATS: Spatial Pattern Analysis Program for Quantifying Landscape Structure. Portland, OR; 1995 (August).
- R Development Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2015.Google Scholar