Malayan kraits (Bungarus candidus) show affinity to anthropogenic structures in a human dominated landscape

Animal movement can impact human–wildlife conflict by influencing encounter and detection rates. We assess the movement and space use of the highly venomous and medically important Malayan krait (Bungarus candidus) on a suburban university campus. We radio-tracked 14 kraits for an average of 114 days (min: 19, max: 218), during which we located individuals an average of 106 times (min: 21, max: 229) each. Most individuals displayed some level of attraction to buildings (n = 10) and natural areas (n = 12); we identified a similar unambiguous pattern of attraction to buildings and natural areas at the population level (of our sample). Snakes remained in shelter sites for long durations (max: 94 days) and revisited sites on average every 15.45 days. Over 50% of locations were within human settlements and 37.1% were associated with buildings. We found generally seasonal patterns of activity, with higher activity in wet seasons, and lower activity in the hot season. These results show frequent proximity between Malayan kraits and humans at the university; thereby, suggesting a near constant potential for human-wildlife conflict. Despite the fact that no snakebites from this species occurred at the university during our study period, substantial education and awareness training should be considered to ensure continued coexistence on campus.


Methods
Study site. The study area covers the campus of Suranaree University of Technology (SUT) and its surrounding landscape in Muang, Nakhon Ratchasima, Thailand (14.879° N, 102.018° E; Fig. 1). The university campus covers about 11.2 km 2 , and comprises a matrix of human modified lands interspersed with mixed deciduous forest fragments (at the onset of this study we identified there were 37 mixed deciduous forest fragments on campus, mean = 7.36 ± 1.48 ha, range = 0.45-45.6 ha [note, "±" is used for standard error throughout the text]). More than 15,000 students are enrolled at SUT, and there are numerous multi-story classrooms, laboratory and workshop buildings, residential housing, parking areas, eating and sports facilities, an elementary school, and a large hospital on the university campus. During the first term of the 2019 school year, 7622 students, as well as numerous SUT staff, lived in on-campus residential areas. The landscape surrounding the university is primarily dominated by agriculture, though there are also patches of less-disturbed areas as well as several densely populated villages and suburban housing divisions among the monoculture plots of upland crops (e.g., cassava, maize, and eucalyptus).
The study site is located within the Korat Plateau region with an altitude range of 205-285 m above sea level. Northeast Thailand has a tropical climate, and the average daily temperature from 1 January 2018 to 31 December 2020 in Muang Nakhon Ratchasima was 28.29 °C, with daily averages ranging from 19.3 to 34.1 °C 38 . The region receives an average annual rainfall ranging from 1270 to 2000 mm 39 . There are three distinct seasons in northeast Thailand: cold, wet, and hot, each are classified by annual changes in temperature and rainfall. Cold season is typically between mid-October and mid-February, hot season is generally from mid-February to May, while the highly unpredictable rainfall of the wet season is predominantly concentrated between the months May to October 39,40 .
Due to the representation of agriculture, semi-urban, and suburban areas with patches of more natural areas all within a relatively small area, we determined the university campus provided an ideal setting to examine how land-use features and human activity influence the movements of B. candidus. Additionally, past studies have indicated northeast Thailand hosts the most bites by B. candidus in Thailand 29,33 , making sites like ours ideal. www.nature.com/scientificreports/ Study animals. We opportunistically sampled Malayan kraits captured as a result of notifications from locals and ad-hoc encounters during transit due to low detectability in visual encounter surveys, in addition to those discovered through unstandardized visual encounter surveys. Upon capture, we collected morphometric data, including snout-vent length (SVL), tail length (TL), mass, and sex (Table 1, Supp. Table 1). We measured body lengths with a tape measure, measured body mass with a digital scale, and determined sex via cloacal probing, all while the snakes were anesthetized via inhaling vaporized isoflurane. We then housed individuals with an SVL > 645 mm and mass > 50 g in plastic boxes (with refugia and water) prior to surgical transmitter implantation by a veterinarian from the Nakhon Ratchasima Zoo. We attempted to minimize the time snakes were in captivity awaiting implantation; however, delays arose due to the veterinarian's availability, the snake being mid-ecdysis, or the snake having a bolus that needed to pass through the digestive tract before implantation (n = 21 implantations, mean = 5.02 ± 0.61 days, range = 0.60-13.02 days). The Nakhon Ratchasima Zoo veterinarian implanted radio transmitters (1.8 g BD-2 or 3.6 g SB-2 Holohil Inc, Carp, Canada) into the coelomic cavity using procedures described by Reinert and Cundall 41 , while the snake was anesthetized. We assigned each individual an ID according to sex and individual detection number (e.g., M02 = a male was the second B. candidus individual documented during the study). We released the implanted individuals as close as possible to their capture locations (mean = 65.31 m ± 13.7 m, range = 0-226.42 m), though on six occasions we moved individuals ≥ 100 m because the individual came from either residential areas or a busy road (all but one were moved < 155 m; see Supp. Table 2 for full details on captures, surgeries, and releases). Cameron Hodges released all snakes within 12 h post-surgery (after nightfall), though on one occasion retained two individuals (M32 and M33) for an additional night post-surgery to avoid heavy rainfall the night of the implantation surgery. We began collecting location data the day after their release. We included all tracking data in analyses, as the animals appeared to move and behave as usual immediately following their release. We radio-tracked 14 individuals (13 males, 1 female) between 8 June 2018 and 24 March 2020 within the SUT study site (Table 1) and classed individuals as adults if the SVL was > 800 mm; thus, nine of the males were adults and four were juveniles (though two of the males had an SVL > 720 mm, and therefore likely sub-adults). The single telemetered female was an adult.
Individual tracking durations varied (mean = 106.46 ± 15.36 days, range = 28.5-222.77 days; Supp. Fig. 1), as many individuals were lost due to unexpected premature transmitter failures (n = 5) or unsuccessful recapture efforts due to individuals sheltering under large buildings as the transmitter reached the end of its battery life (n = 4). We only recorded one confirmed mortality in the study, M01, who was killed by a motorized vehicle when crossing a road (n = 1). Another three individuals were lost due to unknown reasons, which may have been due to premature transmitter failure, mortality, or the animal moving beyond radio signal despite extensive search efforts. Thus, we only successfully recaptured and re-implanted five individuals (M01 once, M02 twice, M07 once, M27 once, and M33 twice). Transmitter batteries generally lasted approximately 90-110 days, so we aimed to replace transmitters after ≥ 90 days of use. At the end of the study, only one individual was successfully recaptured to remove the transmitter. Data collection. We used very high frequency radio-telemetry to locate each telemetered individual on average every 24.20 h (SE ± 0.41, 0.17-410.0 h; see Supp. Fig. 2 for distribution of tracking time lags). We aimed to locate each individual's shelter locations once each day during the daylight (06:00-18:00 h); however, we were occasionally (n = 34 days) unable to locate a snake for several consecutive days when we were unable to www.nature.com/scientificreports/ obtain radio signal due to an individual having moved far away or deep underneath a large structure. There were also a few occasions where we were unable to track snakes due to prolonged and heavy rainfall (n = 4 days), as the moisture damages equipment, or other reasons (n = 4 days). We additionally located snakes nocturnally (18:00-06:00 h) ad hoc and in an attempt to observe nocturnal behaviors and movement pathways when animals were active. We defined fixes as any time a telemetered individual was located, and relocations (i.e., moves) as the occasions where we located an individual > 5 m from its previous known location. Each day we manually honed in on signal via a radio receiver to locate individuals (as described by Amelon et al. 42 , and recorded locations in Universal Transverse Mercator (UTM; 47 N World Geodetic System 84) coordinate reference system with a handheld global positioning system (GPS) unit (Garmin 64S GPS, Garmin International, Inc., Olathe, Kansas) directly above the sheltered snake. We generally approached within one meter of sheltering snakes during daylight to precisely record shelter locations and identify shelter type. Since we could not visually confirm snake locations, we methodically eliminated all possible locations where the snake could possibly be while at close range with the minimum possible gain on the radio receiver.
Telemetered kraits tended to be inactive and sheltering underground during the daylight, thus we were confident that our diurnal location checks would not affect their movements. However, in some cases we resorted to determining an individual's location via triangulation, where multiple lines cast from different vantage points towards the snake intersect on the snake's location on the GPS, allowing us to determine the animal's coordinate location from approximately 10-30 m away. This helped ensure that we recorded locations with greater accuracy when snakes sheltered underneath large buildings, as it allowed us to move away from large structures that hindered the GPS accuracy. This technique was also implemented during some nocturnal location checks when a snake was believed to be active among dense vegetation, in an attempt to prevent disturbance of the animals' natural behavior. While we did hope to gain visual observations of active individuals during the night, we exercised more caution during nocturnal location checks, typically maintaining a minimum distance of approximately 5 m in attempt to lessen the chances of disturbing an active individual's behavior. If the animal was active we recorded the animal's observed behavioral state (i.e., moving, feeding, or foraging). When the radio signal was stable and the individual was not visible, we recorded the animal's behavior as "sheltering". We strived for an accuracy of < 5 m GPS accuracy when feasible during each location check.
For every location fix we recorded the time (dd/mm/yyyy hh:mm), location (Universal Transverse Mercator, UTM; World Geodetic System 84), relocation distance (straight-line distance between the last known location to the new location, relocation/move defined as > 5 m difference), and land-use type (e.g., mixed deciduous forest, human-settlement, semi-natural area, agriculture, plantation; see Supp. Figs. 3 and 4 for photos of land-use types), behavior (e.g., sheltering, moving, foraging, or feeding), and shelter type (e.g., anthropogenic, burrow, or unknown, note we also recorded if we suspected the shelter to be part of a termite tunnel complex due to a close proximity to a visible termite mound; Supp. Fig. 5).
During each location check we recorded the straight-line distance between the current and previous locations (distance moved/step length) with the GPS device. We then used step-lengths to summarize their movements by estimating the mean daily displacement (MDD; the total distance moved divided by the number of days the snake was located) and mean movement distance (MMD; the mean relocation distance, excludes distances ≤ 5). In order to limit biases due to some snakes being located multiple times within a given day/night, we limited our sample for estimating MMD and MDD to only include a single location per day. This was accomplished by manually removing "extra" nocturnal location checks that occurred within the same day, making sure to have all shelter relocations present within the dataset. When calculating MDD, we used the total number of daily location checks rather than the number of days between the individual's tracking start and stop date since there were some days where individuals were not tracked. We also used the same one location check per day dataset to calculate movement/relocation probabilities and to examine each individual's MMD, MDD, and relocation probability for the overall tracking duration as well as for each season.
When feasible, we positioned a Bushnell (Bushnell Corporation, Overland Park, Kansas) time lapse field camera (Trophy Cam HD Essential E3, Model:119837) with infrared night capability on a tripod spaced 2-5 m from occupied shelter sites. We positioned the cameras so that we may gather photos of the focal snake as it exited the shelter site and/or behaviors exhibited near the shelter. We programmed the cameras using a combined setting, including field scan, which continuously captured one photo every minute, along with a motion sensor setting, which took photos upon movement trigger outside of the regular 1-min intervals.
Space use and site fidelity. All analyses and most visualizations were done in R v.4.0.5 using RStudio v.1.4.1106 43,44 . We attempted to estimate home ranges for the telemetered B. candidus individuals using autocorrelated kernel density estimates (AKDEs) using R package ctmm v.0.6.0 45,46 in order to better understand the spatial requirements of B. candidus. However, examination of the variograms revealed that the majority of the variograms had not fully stabilized (i.e., limited evidence of range stability in our sample), and many individuals had extremely low effective sample sizes (21.82 ± 9.75, range = 1.49-135.75; Supp. Table 4). Therefore, we do not report home ranges in this text, as the AKDE estimates would violate the assumption of range residency and either underestimate or misrepresent B. candidus spatial requirements. We also examined the speed estimates resulting from fitted movement models. Resulting variograms and tentative home range estimates are included in a supplementary file for viewing only (Supp. Fig. 6, Supp. Table 4). The original code is from Montaño et al. 47 .
Since our data was not sufficient to estimate home range size for the telemetered B. candidus, we instead used Dynamic Brownian Bridge Movement Models (dBBMMs) with the R package move v.4.0.6 48 to estimate within study occurrence distributions. We caution readers that these are not home range estimates but instead modeling the potential movement pathways animals could have traversed 49  www.nature.com/scientificreports/ pathways and resource use, but it also allows us to examine movement patterns through dBBMM derived motion variance 50,51 . We selected a window size of 19 and margin size of 5, to catch short resting periods with the margin, while the window size of 19 is long enough to get a valid estimate of motion variance when the animals exhibit activity/movement. Contours however are somewhat arbitrary; therefore, we used three different contours levels (90%, 95%, 99%) to estimate dBBMM occurrence distributions (using R packages adehabitatHR v.0.4.19, and rgeos v.0.5.5), and show the sensitivity to contour choice 52,53 . All movement data, either including initial capture locations or beginning with the first location check ~ 24 h post release, was used for production of both the AKDEs and dBBMMs for each individual. We also estimated dBBMM occurrence distributions for each telemetered individual with the exception of M29, which only made three small moves within a burrow complex during the short time he was radio-tracked before transmitter failure.
We compared space use estimates to two previously published B. candidus tracking datasets 34,36 , and one unpublished dataset shared on the Zenodo data repository 54 36 . The previous studies on B. candidus only tracked the movements of a single individual each, had coarser tracking regimes, and used traditional-fundamentally flawed methods 55,56 -to estimate space use 34,36 . Therefore, we ran dBBMMs with these previous datasets using the same window (19) and margin size (5).
To quantify site reuse and time spent at sites (residency time) we used recursive analysis with the R package recurse v.1.1.2 57 . We defined each site as a circular area with a radius of 5 m around each unique location (matching the targeted GPS accuracy). Then we calculated each individual's overall number of relocations, each individual's total number of relocations to each site, and each individual's site revisit frequency and residency time at each unique site. Then we plotted revisited locations on a land-use map with space use estimates (95% and 99% dBBMM) in an attempt to help identify and highlight activity centers for telemetered individuals (see Supp. Figs. 7-13). All maps were created using Quantum Geographic Information System (QGIS v.3.8.2).

Habitat selection. We used Integrated
Step Selection Function models (ISSF) to examine the influence of land-use features on the movements of B. candidus at both the individual and population levels. We included movement data from all male individuals that used more than one habitat feature in our ISSF analysis. Therefore, we excluded F16 and M29 who both only used settlement habitat. Excluding M29 was justified by the individual having been tracked for the shortest duration (19 days) and had the fewest number of moves (n = 3), thus there were not enough relocations for ISSF models to work effectively. Using modified code from Smith et al. 51 that used ISSF with Burmese python radio-telemetry data, we used the package amt v.0.1.4 58 to run ISSF for each individual, with Euclidean distance to particular land-use features (natural areas, agriculture, settlement, buildings, and roads) to determine association or avoidance of features. Cameron Hodges created all land-use shape files in QGIS by digitizing features from satellite imagery and verified all questionable satellite land-use types via on-ground investigation.
The semi-natural areas, plantations, mixed deciduous forest and water bodies (such as irrigation canals and ponds which have densely vegetated edges) were all combined into a single layer of less-disturbed habitats which we refer to as "natural areas". All feature raster layers were then converted into layers with a gradient of continuous values of Euclidean distances to the land-use features, and were inverted in order to avoid zero-inflation of distance to feature values and to make the resulting model directional effects easier to more intuitive. We were able to generate 200 random steps per each observed step (following Smith et al. 51 ), due to the coarse temporal resolution of manually collected radio-telemetry data (i.e., we were not computational limited when deciding the number of random locations). Higher numbers of random steps are preferable as they can aid in detecting smaller effects and rarer landscape features 59 .
To investigate individual selection, we created nine different models testing for association to habitat features, with one being a null model which solely incorporated step-length and turning angle to predict movement 60 , five examining land-use features individually (agriculture, buildings, settlement, natural areas, roads), and the other three being multi-factor models. Each model considers distance to a land-use variable, step-length, and turnangle as an aspect of the model. After running each of the nine models for each individual, we then examined the AIC for each model, point estimates (with lower and upper confidence intervals), and p-values in order to identify the best models for each individual and determine the strongest relationships and trends among the samples. We considered models with ∆ AIC < 2 as top performing models.
We then investigated habitat selection at the population level, including all radio-telemetered except for M29 and F16. We used code from Smith et al. 51 , which was modified code originally from Muff et al. 61 , using a mixed conditional Poisson regression model with stratum specific effects. This model was essentially the equivalent to an ISSF at the population level. Both the step (strata), and the individual (individual ID) are modelling using Gaussian processes. As used for the individual level ISSF, we generated 200 random steps from each location, with a Gamma distribution for step length and Von Mises distribution for turn angle. We created five single factor models using the same land-use features used with the individual level ISSF (i.e., agriculture, buildings, natural areas, roads, and settlements [via the same inverted distance to land-use feature rasters]) with individual random intercepts and slopes. Following Muff et al. 61 , we set a fixed prior precision of 0.0001 for the stratum-specific random effect (i.e., step). We used a Penalized Complexity prior, PC (1, 0.05), for the other random slopes (i.e., individual), and uninformative normal priors, Normal (0, 103), for the fixed effects, as was done by Smith et al. 51 63,64 , we use these metrics as we attempt to maintain standard daily tracks throughout the study period. We calculated MMD and MDD for each individual within the defined seasons, using the methodology described for gaining overall MMD and MDD, and movement probabilities for each individual during each season (the total number of fixes divided by the number of relocations), in order to examine the raw data for possible movement trends.
We also examined dBBMM derived motion variance of all individuals to examine potential variations in temporal activity, and documented observations which may help determine when breeding takes place (e.g., conspecific interactions, presence of sperm plugs in males).
Approval for animal use.   75), but were incalculable for eight individuals (due to lack of model fit). Compared to MDD and MMD, the existing speed estimates (n = 5) appear very weakly connected to MMD (Supp. Fig. 14).
Mean motion variance for all telemetered B. candidus was low, at 1.70 ± 0.18 m (5.52 × 10 -5 -89.73 m). Mean motion variance was lowest for the single telemetered female (mean 0.16 ± 0.04 m), who remained within the same shelter complex (under the F1 building near the northern entrance) for the majority of her tracking duration, including 85 consecutive days (25 January-20 April 2019) spent there.
Occurrence distributions and site fidelity. Individual dBBMM occurrence distributions varied greatly, with the smallest (excluding M29) being the female (F16) with a 99% confidence area of 0.42 ha, and the greatest being the 99% confidence area for M32, at 119.55 ha (Fig. 2). We removed one male (M29) from this summary because the tracking duration was only 21 days with a single small relocation prior to returning back to the previous site. The telemetered male B. candidus (n = 12) had a mean 90% dBBMM confidence area of 6 Habitat use. Habitat use varied across individuals; however, the most frequently used habitat type overall was human settlement, with 51.2% of all fixes (Fig. 3A). Semi-natural areas were the second most commonly used habitat (25.2%), closely followed by mixed deciduous forests (22.8%). The least used land-use types were agriculture (0.5%) and plantation forests (0.3%). Of the points among human settlement habitat, 558 (72.47%) were associated with buildings and 99 (12.86%) were associated with concrete drainage ditches, sidewalks, or other concrete structures.
Shelter use. We determined individuals to be sheltering during 1,443 fixes. The most commonly used shelter type, with 582 (40.3%) fixes, were burrows, however, use of anthropogenic shelters was nearly equal, with 559 (38.7%) fixes (Fig. 3B, Supp. Fig. 5, Supp. Table 6). Anthropogenic shelters included fixes where snakes were directly underneath buildings (507), concrete drainage ditches (37), sidewalks (2), or other anthropogenic structures, while burrows included burrow systems and tunnels excavated by animals, such as rodents, but did not include shelters which appeared to be part of termite mound tunnel systems. Termite mounds/tunnels made up www.nature.com/scientificreports/ 217 (15%) of the shelter fixes. We were unable to identify shelter types during 85 (5.9%) of the fixes (though in as much as 55% of these "unknown" shelter fixes the snake was suspected to be sheltering among dense vegetation). One-hundred and fourteen of the shelters classed as non-"anthropogenic" were within a short distance (< 5 m) of concrete structures, such as drainage ditches or buildings, with 37 of the "burrow" and two "unknown" shelters within a single meter of a building, nine in or within a meter of a concrete drainage ditch, and 19 within a meter of a paved sidewalk.
Trends in foraging sites. Since we attempted to limit nocturnal location checks, it was uncommon for us to track an individual when it was active. Of the 38 occasions where we observed telemetered snakes moving, foraging, feeding, or otherwise active, nine occurred within a single meter of a paved sidewalk, 10 were associated with buildings (either inside or < 1 m to a building edge), and eight occurred in or within a single meter of a concrete drainage ditch (note these instances include initial capture locations). In total, 22 of 37 observations were associated (< 1 m) with concrete structures of some kind. We also gained five observations of B. candidus foraging or moving within agriculture (two within cassava fields, one within a fallow field, one among a fishery, and one on a road-side among a grass field). Several of the other observations occurred near the edge of a body of water, and another two telemetered individuals were observed moving within a meter of a chicken coop.
We recorded two individuals feeding on snakes during the study; both were within anthropogenic land use types (fully reported in 37,65 ). The first occurred in a concrete gutter within 2 m of student housing. The second event occurred within an open atrium within a faculty office building.
Habitat selection. All twelve of the B. candidus individuals included in the ISSF analysis exhibited positive association with natural habitats, and all but two of the individuals (M27 and M35) showed positive association with anthropogenic structures (Supp. Fig. 16). Five models best explained habitat selection across individuals (Supp. Table 7). Top models included model4 (buildings), model6 (natural), model7 (agriculture, natural, and buildings), model8 (roads, buildings, and natural), and model9 (roads, agriculture, and natural).
Credible intervals were quite broad and sometimes overlapped zero for several of the individual's models, thus limiting our ability to draw inferences. Much of the model uncertainty is likely resulting from the coarse resolution tracking data and the few and infrequent movements by the study animals. Interactions between steplength and distance to land-use features appeared nonexistent. However, the population level model revealed a potentially broader trend while absorbing some of the individual heterogeneity (via the random effect), with weak but positive associations with natural habitats (mean estimate = 0.0235, 95% CrI 0.0114-0.0416), buildings (mean estimate = 0.0094, 95% CrI 0.0038-0.0165), and settlements (mean estimate = 0.0059, 95% CrI 0.0001-0.0125). A possible weak negative association with agriculture was also present (mean estimate = − 0.0018, 95% CrI − 0.0058-0.0020; Fig. 4). Estimates for the influence of distance to land-use features on step lengths resulting Seasonality and temporal activity patterns. Males moved furthest distances (MMD) in the cold season, but longer distances per day (MDD) and more frequently in the wet season. In contrast males moved least (MMD and MDD) in the hot season (Fig. 5B, Supp. Table 8). Males were more likely to move between any given daily datapoint in the wet season than in other seasons. Mean motion variance was highest during the cold season (2.579 ± 0.332 m), where four larger adult males peaked in motion variance resulting from particularly large movements coinciding within the first two months of the cold season (Fig. 5A). There notably were several particularly high peaks in motion variance from a few different individuals (M02, M07, M12, and M32) both years. These four highest peaks all coincide within the first two months of the cold season, October and November. Similarly, twelve of fourteen particularly large movements (≥ 395 m, by 6 individuals) documented occurred within the late wet season (n = 4) or early cold season (n = 8), with the remaining two occasions occurring within the late cold season (January). In contrast, motion variance was lowest in the hot season (0.481 ± 0.042 m), and the wet season average was roughly halfway between the other seasons (1.174 ± 0.105 m).
In contrast to large moves, six of the male Malayan kraits showed prolonged stationary periods, remaining inactive within the same shelter for ≥ 20 consecutive days (mean = 35.05 ± 7.76 days, range = 20-94 days, n = 9). Shelter emergence times. We gathered a total of 1,160,970 time-lapse camera photos. Of these, we could only identify focal animals in 75 photos, from six different individuals, on 14 different occasions (i.e., independent nights), with a mean of 5.36 (range = 1-18) photos each occasion. During different occasions, individuals generally either peaked their heads out and slowly exited shelter sites (n = 7), or simply exited the shelter site and immediately moved away (n = 4), not to return again. Two individuals (M28 and M36) were photographed active during the night, presumably foraging near to the shelter site, before returning to the same shelter. On another occasion an individual (M12) spent a few minutes lying just outside the shelter before moving off. On camera, individuals tended to exit shelters and begin moving and/or foraging shortly after sundown, with all photographed individuals moving outside the shelter sites between 19:00 and 22:30 h, and with shelter site activity peaking at 19:30 h (Fig. 6).

Discussion
This study yields insight into a secretive snake species of medical importance-Bungarus candidus. We use a variety of approaches (i.e., dBBMMs, ISSFs, and AKDEs) to analyze low-resolution animal positional data. Despite identifying site fidelity, we were unable to estimate home ranges (via AKDE) of B. candidus due to a surprising lack of evidence of range stability, likely due to limited number of relocations, resolution, and tracking duration. However, dBBMM movement pathway estimates suggest an average area of 11.22 ± 4.37 ha (95% confidence area) potentially used during the study period per male individual. Bungarus candidus in our study were seasonal and frequently sheltered, foraged, and remained in human settlements, with high site fidelity and long-term use of shelters in areas constantly co-habited by humans. Despite B. candidus consistently occupying areas near humans (and potentially being most active during times of relatively high human traffic) there were no incidents of snakebite attributable to B. candidus in our site during the study period. Space use. Space use is frequently dependent on the environment 66 ; however, we witnessed very little difference between our kraits' confidence areas (n = 13, 95% confidence area mean = 11.22 ± 4.37 ha), and the previous studies in much less urban areas (n = 3, 95% contour estimate mean = 8.49 ± 2.4 ha). This similarity appears to contrast with findings from Tucker et al. 's 67 review, that showed mammals moving less in more anthropogenic environments. While the comparison is limited by the sample size from less-developed areas and differences in tracking duration, the initial lack of difference may indicate krait movement patterns are insensitive to human presence, or that changes are occurring on a scale undetectable by our dataset (e.g., movement changes occurring in < 24 h, changes in activity times). A lack of flexibility may lead to heightened vulnerability to human-snake conflict 68 , and also presents a foundation for exploring how such conflict may be predicted by human behaviors.
Site fidelity and habitat selection. All but one of the 14 radio-tracked individuals exhibited site fidelity. Snakes revisited both natural and developed areas repeatedly, specifically reusing termite mounds, tunnel systems, and crevices under anthropogenic structures (Fig. 3B, Supp. Fig. 5). www.nature.com/scientificreports/ We found most individuals to be associated with natural areas, commonly moving through and sheltering among available mixed deciduous forest and semi-natural areas. But most individuals also tended to use settlements more so than less-disturbed habitats, and were positively associated with anthropogenic structures, such as houses and university buildings. Snakes commonly used these areas as shelter sites (likely leading to the building association in the population model). Several studies document habitat use and movement with snakes in landscapes associated with humans 69,70 , though relatively few evaluated or even documented sheltering under concrete anthropogenic structures 71,72 . In cooler temperate areas, some snakes frequently used anthropogenic structures as refugia, as they provided more suitable thermal qualities than the available natural shelters 73,74 . More relevant to our study, Wolfe et al. 75 reported that several telemetered Pseudonaja affinis, a large diurnal elapid species from Australia capable of surviving in urban areas, occasionally sheltered underneath housing and paving stones. Our study demonstrates that even in a warm tropical climate, potentially dangerous snakes may be attracted to buildings and may use them as shelters even more frequently than available natural refuges.
Why are kraits using human settlements? These sites may modulate temperature by increasing ambient surrounding landscape temperature 76 . So perhaps Malayan kraits are using these insulated concrete structures to thermoregulate. While some suggest thermoregulation is less important to reptiles in tropical environments 77,78 , the concrete structures may aid animals in avoiding heat 79 . Findings from several studies suggest that habitat and shelter selection by some snakes is influenced by prey availability 80,81 , though in temperate areas it appears that thermoregulation may outweigh prey availability when selecting habitats 82,83 . Concrete structures likely house prey species such as other snakes, anurans, lizards, and mammals. Rats for example often use these structures because they fulfill three criteria: access to food via refuse, access to water via runoff, and access to shelter as a central foraging location 84 . Thorough investigations into these dynamics have yet to be held in the tropics, where thermoregulation, though still important 79,85 , likely has less of an impact on snake habitat use. Our study does provide evidence that B. candidus do forage in areas proximal to shelter sites, sometimes even returning to the same shelter site following active foraging activities (as observed through several camera trap images, also see Hodges et al. 37 . We found no clear evidence for road avoidance and limited evidence of attraction to roads in some individuals-but not the population model. This may help explain the numerous documented B. candidus road mortalities in our study site (C.W. Hodges, personal observation), though we suspect culverts were used for at least some road crossings, as appears to be the case for Ophiophagus hannah 86 . Furthermore, while our ISSF models did not reveal clear evidence of avoidance of agriculture, despite agriculture being the most widespread land use type in the area, few individuals used it for shelter. Agricultural areas could have occasionally been used as foraging sites during the night (such as observed in M27 and M33), but not for sheltering due to the lack of suitable shelters sites. More commonly, B. candidus sheltered in unmanicured field margins, similar to findings from Knierim and colleagues, which observed this in a single telemetered B. candidus 36 as well as several radio-tracked B. fasciatus 27 living among agriculture. They attributed the non-use of monoculture plots to insufficient shelter availability resulting from frequent and regular disturbance to top soil by human activities. It is also possible that B. candidus may have tended to avoid agricultural areas due to the increased risk of mortality 35,87 (C.W. Hodges, unpublished data) but we lack direct evidence to support this.
Seasonality. Tracked male B. candidus had more frequent but shorter relocations in the wet season, and moved less frequently, but covered greater distances, in the cold season, with the least movement in the hot season. The wet season may present abundant resources, influencing movement patterns 10,88 . Prey availability can influence snake movement and activity in the temperate 89 and the tropics 90,91 .
The fewest resources likely exist in the hot season, coinciding with the lowest B. candidus activity. Other reptiles in northeast Thailand, including Indotestudo elongata 92 , Ophiophagus hannah 93 , Trimeresurus macrops 94 , and Python bivittatus 51 responded similarly to dry seasons. Many reptiles reduce activity during particularly hot and dry periods to conserve energy and water [95][96][97][98] . Bungarus candidus are capable of prolonged inactivity; as we documented six adult males which remained within a single refugia for extended periods (up to 94 days, and this could have continued even longer, as this individual, M35, was recaptured during this period). Snakes may have occasionally exited shelters for nocturnal foraging and returned to the same shelter, but given camera traps failed to reveal movement from shelters, we consider it unlikely (though other exits could exist). Snakes sheltering long-term may have also foraged and fed fossorially/opportunistically from within the shelter complexes. The long-term dormancy might alternatively be extended mating behavior, though we did not detect conspecifics present during recaptures or shelter site camera trapping.
Peaks in cold season motion variance could be attributed to male B. candidus searching for mates over great distances-max of 1130 m (M02 expelled a sperm plug during processing on 17 September 2018, just prior to cold season). We suspect-but did not confirm-that the telemetered female nested beneath a building foundation in a refuge system where she remained for 85  www.nature.com/scientificreports/ Implications for human safety. There is extreme spatial and temporal overlap between humans and kraits in our study, similar to Glaudas 25 who documented overlap of Daboia russelii, another tropical Asian medically significant species. Our findings suggest an increased need for education and awareness among Thai communities, especially when paired with insights provided by Hodges et al. 37 , demonstrating that short-distance translocation for B. candidus is ineffective for preventing long-term conflicts. Using these data combined with occurrence data may aid in predicting where snakebites are most likely to occur and perhaps elucidate further preventative methods 100 . Such need for increased data is particularly apparent for other Bungarus spp. also known to occupy human-modified areas; for example, B. caeruleus that is responsible for thousands of deaths across South Asia annually 23,101,102 . Continual use of buildings for foraging and shelter highlights the need for access to antivenom and precautionary measures to prevent snakebite. In Nepal, Samuel et al. 103 , described how an education program -highlighting flashlight use, wearing boots when working in fields, eliminating openings in housing-led to reduced snakebites and mitigated the consequences. Another study from Nepal identified using a bed net whilst sleeping offered strong protection against snake bite 104 . As tracked B. candidus tended to avoid open land features like agriculture, land managers could potentially limit conflict by reducing vegetation cover near buildings-but this is yet untested, and cannot come at the expense of broader biodiversity benefits of university campuses 3 .

Limitations
Although the study sample of B. candidus (n = 14 [adult males = 9, juvenile males = 4, adult female = 1]) somewhat lower than the average sample for snake spatial studies over the last 20 years; our tracking intensity was comparable to the average 56 , and we were likely more consistent than average snake telemetry studies.
We attempted appropriate methods for comparing movements across seasons 105 , but the models failed as a result of our data structure forcing us to resort to weaker proxies (MMD and MDD). We suspect the short tracking duration (mean = 105.56 ± 15.16 days), coarse tracking resolution (mean = 24.03 ± 0.43 h), limited relocations (25.71 ± 4.09) and unequal sampling across seasons (308 fixes during the wet season/1505 fixes total) led to low effective sample sizes and thus model failure. Our inferences are further limited to shelter site selection as opposed to foraging and finer scale movements.
Prompted by the STRANGE framework 106 we consider the following as potentially impacting the generalizability of our findings. Social background or rearing history likely had no impact on the study with limited evidence of sociality in kraits, and all study animals were free-ranging with limited periods in captivity (mean = 5.02 ± 0.61 days). Trappability likely influenced the individuals we captured as we could not implement systematic trapping. For example, we more likely tracked individuals nearer humans, evidenced by 10/14 individuals originated from notifications from human-snake interactions. Our results cannot be generalized to female B. candidus because we tracked a single tracked female. We doubt that individuals acclimated or habituated to observer presence, as individuals were exposed to near continuous presence of human disturbance at the study site. In addition, our focus on day time tracks limited our presence during krait's activity periods as evidenced by our camera trap data. It is difficult to gauge the impacts of natural changes in responsiveness caused by unmeasured biological cycles (e.g., diel cycles, stress response attributed to mating season, prey populations). Genetic make-up may have influenced observed behaviors because of the small spatial scale, and high likelihood of closely related individuals. We suspect that experience (i.e., learning to avoid humans) played a minor role due to the short duration of the tracking, but radio-transmitter replantation for five individuals may have constituted traumatic events to prompt greater human avoidance; however, this was not apparent in individual step-selection models.

Conclusion
We identified that tracked krait space use during the study consistently overlapped with human structures and that snakes avoided using open landscape features. We also highlight heightened activity during the wet season and during the early hours of the night, suggesting two peaks in possible human encounters and resulting conflict. But this marked co-occurrence with humans in the absence of recorded bites suggests an alternative. Redoubling efforts to raise awareness of the habits and likely locations of these medically significant snakes is likely important to maintaining this harmonious balance of human-snake coexistence. Further, investigations into finer scale movements (e.g., building entry routes) could yield insights into further drivers enabling the coexistence we observed.