Habitat selection by free-roaming domestic dogs in rabies endemic countries in rural and urban settings

Study sites and study design

The study was performed in the frame of a dog ecology research project, with details on the study locations published elsewhere^15,42,43. For the current study, five study sites located in Indonesia and Guatemala were included. Site selection was carried out by each country’s research team, taking into consideration rural and urban settings, as well as differing expected number of dogs present at each location. The Indonesian study sites were semi-urban Habi and rural Pogon, in the Sikka regency, at the eastern area of Flores Island (Supplementary Fig. 6). In Guatemala, the study sites were Poptún (urban setting), Sabaneta and La Romana (both rural settings), located in the Guatemalan department of Péten, in the northern part of the country (Supplementary Fig. 7). Data were collected during May to June 2018 in Guatemala and from July to September 2018 in Indonesia.

In each location, a 1 km² area was predefined using Google Earth within which the study took place. The 1 km² area was chosen because of the research goals of another part of the project, investigating the contact network of the dogs¹⁵. Within these areas, the teams visited all dog-owning households. In each household, the study was presented to an adult of the family, who was then asked if they owned a dog and if they were willing to participate in the study. After the dog owner’s oral or written consent was granted, a questionnaire was answered, and the dogs collared. The handling of the dogs was performed by a trained veterinarian or a trained veterinary paramedic of the team.

The questionnaire data was collected through interviews with the dog owners. Multiple dogs per household could be included as multiple entries in the questionnaire. The detailed questionnaire contains information on the household location, dog demographics (age, sex, reproductive status) and management (dog’s purpose, origin, confinement, vaccination status, feeding and human-mediated transportation within and outside the pre-determined area).

All dogs of a household fulfilling the inclusion criteria were equipped with a geo-referenced contact sensor (GCS) developed by Bonsai Systems (https://www.bonsai-systems.com), containing a GPS module and an Ultra-High-Frequency (UHF) sensor for contact data recording^43,44. GCS devices report a 5-m maximum accuracy, a run-time of up to 10 years, can store up to 4 million data points and carry a lithium-polymer-battery (LiPo). For this study, only GPS data were analysed. The GCS were set to record each dog’s geographical position at one-minute intervals. Dogs remained collared for 3 to 5 days with the duration of the data collection being limited by the device’s battery capacity, as batteries were not re-charged or changed during the study. Throughout the time of recording, date, hour, GPS coordinates and signal quality (HDOP) raw data were collected by the GPS module and amassed into the workable databases.

Exclusion criteria were dogs of less than four months of age (since they were not big enough to carry a collar), sick dogs and pregnant bitches (to avoid any risk of stress-induced miscarriages). Reasons for non-participation of eligible dogs included dog owner’s absence, dog’s absence, inability to catch the dog, and refusal of participation by the dog owner. In addition, dogs foreseen for slaughtering within the following four days were excluded in Indonesia to ensure data collection for at least four to five days. All dogs included in this study were constantly free roaming or at least part-time (day only, night only and for some hours a day).

Human and/or animal ethical approval were obtained depending on the country-specific regulations. All the procedures were carried out in accordance with relevant guidelines. Ethical clearance was granted in Guatemala by the UVG’s International Animal Care and Use Committee [Protocol No. I-2018(3)] and the Community Development Councils of the two rural sites, which included Maya Q’eqchi’ communities⁴⁵. In Indonesia, the study was approved by the Animal Ethics Commission of the Faculty of Veterinary Medicine, Nusa Cendana University (Protocol KEH/FKH/NPEH/2019/009). In addition, dogs that participated in the study were vaccinated against rabies and/or dewormed to acknowledge the owners for their participation in the study.

Data cleaning

Data were stored in an application developed by Bonsai Systems compatible with Apple operating system (iOS iPhone Operating Systems), downloaded as individual csv file for each unit, and further analysed in R (version 3.6.1)⁴⁶.

The GPS data were cleaned based on three automatised criteria. First, the speed was calculated between any two consecutive GPS fixes, and fixes with speed of > 20 km/h were excluded, given the implausibility of a dog running at such speed over a one-minute timespan⁴⁷. It is noteworthy that car travel causes speeds over 20 km/h. However, as we were interested in analysing the dog’s behaviour outside of car transports, removing these fixes was in line with our objectives. Second, the Horizontal Dilution of Precision (HDOP), which is a measure of accuracy⁴⁸ and automatically recorded by the devices for each GPS fix, was used to exclude fixes with low precision. According to Lewis et al.⁴⁹, GPS fixes with HDOP higher than five were excluded, which deleted 1.3% of data in Habi, 2.2% in Pogon, 3.3% in Poptún, 1.8% in La Romana and 2.1% in Sabaneta. Third, the angles built by three consecutive fixes were calculated for each dog. When studying animals’ trajectories as their measure of movement, acute inner angles are often connected to error GPS fixes⁵⁰. The fixes having the 2.5% smallest angles were excluded, to target those fixes with highest risks of being errors, while balancing against the loss of GPS fixes due to the cleaning process. With the exclusion of the smallest angles, 2.6% of data were deleted in Habi, 3% in Pogon, 2.9% in Poptún, 2.6% in La Romana and 2.7% in Sabaneta. After the automatised cleaning was concluded, 18 obvious error GPS fixes (unachievable or inexplicable locations by dogs) still prevailed in the Habi dataset and were manually removed.

Habitat resource identification and calculation of terrain slope

To analyse habitat selection of the collared FRDD, resources were delimited by a 100% Minimum Convex Polygon (MCP) including all cleaned GPS fixes per study site, using QGIS⁵¹ (Fig. 1).

Resources were defined by taking into consideration the following criteria: resources are (i) likely to impact upon movement patterns of dogs, (ii) identifiable by landscape satellite topography, and (iii) chosen considering information on relevant gathering places for FRDD observed by the field teams. Three resources were disclosed in all study sites: buildings, roads and vegetation coverage. All habitat relevant resources were manually identified within the available area (MCP) in QGIS using satellite imagery. All building-like structures were identified using vector polygons and summed under the layer “buildings”. Roads were identified and manually traced using vector lines in all sites, except in Poptún where the roads were automatically traced using an OpenStreetMap road layer of the area (https://www.openstreetmap.org/export). A buffer vector polygon was generated to encompass the full potential width of the roads, with a 5 m width in Habi and Poptún (semi-urban and urban site) and a 2 m width in Pogon, La Romana and Sabaneta (rural sites). In Habi, a “beach” layer was defined by generating a five-meter buffer from the shoreline in both directions using a vector polygon. The layer “sea” was defined as the vector polygon resulting from the difference between the MCP sea outer limit and the beach buffer polygon. Vegetation coverage was distinct between study sites with sparse vegetation and bushes present in all sites except Pogon, and dense forest-like vegetation present in La Romana and Pogon. These two types of vegetation were defined as “low” and “high vegetation”, respectively. In Habi and La Romana, “low” and “high vegetation”, respectively, were manually identified using vector polygons and summarised under the respective layers. Finally, open field in Habi, high vegetation in Pogon and low vegetation in Poptún, La Romana and Sabaneta were the last vector layers to be established since they represented the difference between all other polygon vector layers and the MCP total area. After all resource vector polygons had been created, an encompassing vector layer was generated by merging all resource polygon vectors for final resource classification (Fig. 2). As part of the resource classification in Habi, the airport terminal and runaway as well as waterways enclosed in the MCP area were identified but excluded from the analysis.

After the construction of the habitat resource layers, all GPS fixes were assigned to the respective resource they were located, using the QGIS join attributes by location algorithm. Fixes located exactly on the MCP border in Indonesia were not classified automatically and had to be manually classified to the respective resource.

In non-flat topographies (all locations expect Habi) we tested the hypothesis of whether the steepness would influence the dogs’ movement patterns. The degrees of slope were calculated using a 30-m raster-cell resolution (STRM 1-Arc Second Global, downloaded from the United States Geological Survey (USGS) Earth Explorer, https://earthexplorer.usgs.gov/). The slope was assigned by the QGIS join attributes by location algorithm to each GPS fix.

Statistical analysis

To quantify habitat selection in each study site, we compared resources used by the dogs with the resources available, according to Freitas et al.⁵². Adapting the methodology applied by O’Neill et al.¹⁸, the observed number of GPS fixes for each dog was used to generate an equivalent number of locations that were randomly distributed within the MCP area using the Random points in layer bound vector tool from QGIS. For example, if dog “D300” had 100 recorded GPS fixes, 100 random points were generated within the MCP of the respective study site and assigned to “D300”. Random points were then assigned to the respective resources and slope of that location, as previously done with the observed GPS fixes. Using this approach, the habitat resources used by each dog could be compared to the available resources in the respective study site, using a regression model.

Observation independence is a fundamental presupposition of any regression model. However, the spatial nature of the point-referenced data permits perception of spatial dependence. In our dataset, spatial autocorrelation was proven for all study sites using the Moran’s I test. Therefore, we applied a spatial regression model, which takes into consideration spatial autocorrelation while exploring the effects of the study variables. A mixed effects logistic regression model accounting for spatial autocorrelation was created to quantify the effect of variables on used (i.e. observed GPS fix) versus available (i.e. randomly generated GPS fixes) resources, using the fitme function in the spaMM package in R^53,54. The model’s binary outcome variable was defined as either observed (1) or random (0) GPS fix, i.e. the dog being present or absent from a position. The explanatory variable was the resource classification with “buildings”, “roads”, “low vegetation”, “beach”, “sea” and “open field” as levels in Habi; “buildings”, “roads” and “high vegetation” in Pogon; “buildings”, “roads”, “low vegetation” in Poptún and Sabaneta; and “buildings”, “roads”, and “high” and “low vegetation” in La Romana. Different habitat resources were used interchangeably as reference level. In all study sites except Habi, the slope was included as an additional explanatory variable. As observations were not evenly distributed in time, with less observations recorded towards the end of the study, a variable ”hour” was added as an additional continuous fixed effect.Each observed GPS fix was assigned to the hour of its record, with the earliest timestamp registered in each study site being assigned the hour zero. The randomly generated points were randomly assigned to an hour within the determined time continuum of the observed GPS fixes. As our focus was investigating habitat selection at a population-level, we assumed there was no within-dog autocorrelation (space/time) and each dog was independent and exhibited no group behaviour³⁸. Still, to partially account for spatial autocorrelation of each dog’s household, the random effects included in models were defined as each dog’s household geographical location recorded during fieldwork by a GPS device. The restricted maximum likelihood (REML) through Laplace approximations, which can be applied to models with non-Gaussian random effects⁵⁵, and the Matérn correlation function were used to fit the spatial models with the Matérn family dispersion parameter ν, indicator of strength of decay in the spatial effect, was set at 0.5⁵⁴.

Source: Ecology - nature.com