Habitat selection by bluetongue lizards (Tiliqua, Scincidae) in tropical Australia: a study using GPS telemetry

Abstract

Background

Emerging global positioning system (GPS) technologies can clarify movement patterns of free-ranging animals in far more detail than has been possible with previous methods. We conducted long-term (mean, 65 days; maximum, 221 days) GPS radio-tracking of 41 northern bluetongue lizards (Tiliqua scincoides intermedia) and 8 centralian bluetongue lizards (T. multifasciata) at two study sites in northwestern Australia, close to the border between Western Australia and the Northern Territory.

Results

Individuals of both species spent long periods within small and distinctive habitat patches, interspersed with longer directional relocations from one patch to the next. Our sampling showed that these patches of core activity differed significantly from the surrounding landscape in several respects. The patches provided relatively shaded, cool, and damp conditions, with higher grass and more leaf-litter cover. The location of these patches in the landscape is probably determined by drainage patterns, soil moisture-holding ability, and stochastic recruitment of shade trees.

Conclusions

These scattered patches provide a critically important habitat for lizards (and probably, other taxa) within this hot dry landscape. Future conservation and management strategies need to prioritize the retention of such sites, at a spatial scale that allows animals to move between them.

Background

In all vertebrate species, organisms must move about for such functions as resource acquisition, mate location, predator avoidance, and thermoregulation [1]. Such movements are driven by complex processes that operate over a range of spatial and temporal scales [1–3]. Understanding the factors that influence movement patterns and habitat selection is central to managing species, particularly taxa that are affected by anthropogenic activities [4, 5].

Most landscapes are highly heterogeneous, and incorporate environmental variation in many abiotic and biotic conditions that can directly affect organismal fitness [6]. The tropical savanna of northern Australia represents one such region that has an extensive heterogeneous matrix imposed by both a highly seasonal climate [7] and anthropogenic influences, such as fire and grazing [8, 9]. These factors influence soil nutrients and water availability and the distribution of vegetation [10, 11]. In such a patchy habitat, the location that an individual inhabits can influence its physiology [12, 13], the availability of food resources [14, 15], its exposure to predation [16, 17], and its reproductive opportunities [18, 19].

For ectotherms such as reptiles, microhabitats often differ in the ease with which an individual can regulate its body temperature [20–24]. In temperate climates, the greatest challenge for an ectotherm may be to obtain body temperatures high enough to facilitate activity, prompting reptiles to shuttle between sheltered microhabitats and more exposed sun-basking sites to maintain suitable body temperatures [12, 20, 25, 26]. In the tropics, where ambient temperatures are high all year round [27], reptiles are able to maintain high, stable temperatures without the need for overt thermoregulation: that is, their selection of habitats is more likely to be driven by the need to avoid overheating than to find basking sites [28, 29]. Thus, tropical reptiles may select thermally distinctive microhabitats not to heat up, but to avoid overheating and water loss [30–32]. In the tropics, this ‘heat avoidance’ behavior probably predominates most of the year, whereas many temperate reptiles face similar thermal challenges for shorter periods (for example, during hot, dry summers [33]). The large temperate-zone scincid lizard Tiliqua rugosa predominately exhibits a ‘heat seeking’ thermoregulatory behavior, but switches to a ‘heat avoidance’ strategy to avoid dangerously high temperatures, and to conserve energy reserves and water during summer months [33].

To compare microhabitats from an organism’s perspective, we need to define those patches of microhabitat based on movement patterns of free-ranging organisms [34]. The movement ecology paradigm [3] provides objective techniques with which to differentiate between alternative modes of movement, potentially allowing us to compare areas that induce one movement mode (for example, long-distance, rapid dispersal) from areas that induce other movement paths (for example, sedentary, with frequent return to the original location). The logistical limitations of obtaining continuous, accurate locational data have largely precluded this approach with small vertebrates [35]. Recent developments in global positioning system (GPS) technology provide ways to overcome this obstacle [35–38].

The tropical Kimberley region of northern Australia supports a rich diversity of reptile taxa (39% of the Australian species total), including two species of bluetongue lizard (northern bluetongue Tiliqua scincoides intermedia and centralian bluetongue T. multifasciata [39]). Both are heavy-bodied, short-limbed omnivorous scincids [40, 41]. In this study, we use the movement ecology framework to examine mechanisms that influence the movement paths of free-ranging individuals of these two lizard species in tropical Australia, to clarify how these animals respond to habitat heterogeneity. In other analyses based on the same dataset, we have used a maximum likelihood approach to differentiate objectively between interpatch versus intrapatch movement phases by these lizards, and examined the influence of abiotic factors on dispersal rates [42]. Here, we extend this previous research to examine effects of biotic and abiotic characteristics on the lizards’ selection of core-area habitats, to better understand the factors that influence movement phases and habitat selection.

Results

Over two consecutive wet seasons (mid-November to mid-March), we radio-tracked a total of 49 bluetongue lizards (41 northern bluetongues, T. s. intermedia, and 8 centralian bluetongues, T. multifasciata) at two sites in the wet-dry tropics of northwestern Australia. Over the course of this study, we obtained >60,000 locations from the GPS system allowing us to distinguish clearly between core areas of high intensity use and areas that were used only infrequently by our radio-tracked lizards.

Point level

Overall, we conducted a total of 951 point-level habitat surveys (311 at the Northern Territory (NT) study site and 640 at the Western Australia (WA) study site). We found the lizards inactive (coiled within a refuge) at 180 (60%) of records at the NT study site, and 549 (88%) of records at the WA study site. Active lizards were generally found in relatively open habitats.

Lizards at the NT study site were predominately located (90% of locations) in scattered open forest (SF) with a native grass or spinifex understory (Figure 1a, Figure 2a). Lizards rarely occurred in scrubby areas dominated by Acacia leucophloea, or garden habitats (Figure 2a). The SF macrohabitat contained many flooded box or coolibah (Eucalyptus microtheca) with ghost gums (Corymbia bella), Cooktown ironwood (Erythrophleum chlorostachys) and Kimberley bauhinia (Lysiphyllum cunninghamii) also common. At the WA study site, lizards were predominately located in patches of riparian vegetation (R, 56%; Figure 1d), gardens (G, 14%), and built environments (BE, 12%; see Figure 2a). The dominant tree species in riparian areas was neem (Azadirachta indica, an exotic tree), with a midstory of Acacia species and an understory dominated by buffel grass (Cenchrus ciliaris, an introduced grass), exotic vine and mixed weeds.

Figure 1

Typical habitats predominantly utilized by bluetongue lizards at two sites within northern Australia. (a,b,c) Keep River National Park, NT and (d,e,f), Kununurra, WA. (a) Scattered open forests dominated by flooded box or coolibah, Eucalyptus microtheca. A typical core area is located at the base of the tree, typical non-preferred areas occur around the periphery. (d) A typical patch of riparian vegetation in the background with non-preferred open paddock habitat in the foreground. Smaller inserts display typical microhabitat characteristics that are associated with core areas: (b) native jasmine, (c) white currant bush, (e) buffel grass, (f) rough-leafed fig.

Figure 2

Percentage of locations in which a radio-tracked lizard was located: (a) by macrohabitat, (b) by microhabitat, (c) by ground-level structure. Descriptions of abbreviations are given in Table 4.

In relation to microhabitat, lizards at both study sites were predominately located within grass (G) and shrubs (S) (NT: 46% and 22%, respectively; WA: 53% and 13%, respectively; see Figures 1 and 2). At the NT study site, lizards (especially T. multifasciata) were also frequently located within spinifex (SX) (16%; see Figure 2b). Within the NT study site, shrubs typically occurred in discrete patches dominated by white currant bush (Fleuggea virosa melanthesoides), wild grape (Ampelocissus acetosa), native jasmine (Jasminum aemulum), stinking passionflower (Passiflora foetida), rough-leafed fig (Ficus opposita), and Acacia species. Often, these patches occurred around the base of a tree (Figure 1a). At the WA study site, shrubs primarily comprised Acacia sp. and exotic garden species. The dominant grass selected by lizards at the WA site was buffel grass (Figure 1e), a perennial grass up to 1 m high with a thick dome-shaped coverage (unlike grasses in the surrounding open areas, which were typically erect with a small basal area).

In relation to ground structure, at both study sites, lizards were found primarily within grass and leaf litter (LL) (NT: 35% and 27%, respectively; WA: 47% and 18%, respectively). At the WA study sites, the lizards often used burrows (B) (15%, probably excavated by varanid lizards), and artificial structures (13%, built structure (BS) and hard cover (HC) combined). At the NT study site, lizards sometimes used burrows (7%); spinifex (15%) was primarily utilized by T. multifasciata (Figure 2c).

Despite considerable differences in overall habitat characteristics of the two study sites, telemetered lizards chose broadly similar microhabitats (based on vegetation characteristics) at both sites (Table 1). Areas where lizards were located had a low proportion of bare ground (<8%), a high coverage of grass and leaf litter (50% and >34%, respectively) and usually one or two trees and three or four shrubs within 5 m of the site selected by the lizard (Table 1; Figure 1).

Table 1 Habitat variables at sites with located radio-tracked bluetongue lizards in tropical Australia

Effects of vegetation structure

At the NT study site, principal component analysis (PCA) based on the original eight vegetation structure variables identified four components that explained 83.27% of variation in habitat parameters (Table 2). The first and second axes from this PCA of vegetation structural data explained 42% and 18% of variation, respectively (eigenvalues 3.37 and 1.44). The first axis (PC1) was positively associated with number of trees and shrubs, and leaf-litter depth (Table 2). The second principal component (PC) axis (PC2) also described the number of trees (Table 2). The remaining PC axes described 23% of the variation, and again represented trees, shrubs and understory cover (PC3), and trees and shrubs (PC4) (Table 2).

Table 2 Eigenvectors from principal component analysis (PCA) of habitat variables

At the WA study site, the PCA also reduced the eight vegetation structural parameters to four PC scores, which together explained 77.37% of the variation in the dataset (Table 2). The primary axis (PC1) for the PCA explained 34.12% of the variation and was composed primarily of the number of trees and to a lesser extent the number of shrubs. PC2, which explained 15.74% of the variation, described the number of trees and shrubs, and amount of bare ground. PC3 (14.83% of the variation) also had a high loading on number of trees as well as leaf-litter depth. PC4 (12.67% of the variation) was associated with number of trees (Table 2).

We performed nominal logistic regressions to relate type of habitat (core vs. adjacent, less frequently used) to the four vegetation structure PC axes derived from our quadrat surveys. Habitat selection by both T. s. intermedia, and T. multifasciata at the NT study site was significantly positively related to PC1 and negatively related to PC2 (Table 3). Thus, lizards spent much of their time within core areas with dense vegetation, especially with more trees and shrubs and thus deeper leaf litter. The same pattern was seen for T. s. intermedia at the WA study site, with core areas containing dense vegetation (trees and shrubs) and minimal exposed ground (PC1 and PC2; see Tables 2 and 3).

Table 3 Nominal logistic regressions showing differences between ‘core areas’ and less frequently used areas

Biophysical factors

At the NT study site, all biophysical parameters that we measured were significantly linked to the PC axes derived from our data on vegetation characteristics. All PC components affected mean temperature (all P < 0.03); PC1 and PC2 affected soil-water-holding capacity and minimum temperature (both have P < 0.02), while PC4 also affected minimum temperature (P = 0.03). At the WA study site, PC1 influenced soil-water-holding capacity (P < 0.04), maximum temperature (P < 0.03), and minimum temperature (P < 0.02). Minimum and mean temperatures were also significantly influenced by PC3 (P = 0.002, 0.04 respectively).

At the NT study site, core areas differed significantly from adjacent, less frequently used parts of the lizards’ home ranges for all biophysical variables except minimum temperature (soil-water-holding capacity: F _1,192 = 30.79, P < 0.0001; mean temperature: F _1,114 = 107.29, P < 0.0001; maximum temperature: F _1,114 = 99.74, P < 0.0001; Figure 3a). There was no significant difference between minimum temperature at used versus unused microhabitats in the NT study site (F _1,114 = 0.72, P = 0.398). On average, minimum temperatures were 24°C in both core and adjacent patches. Mean and maximum temperatures were far lower in core patches than in nearby areas (31 ± 0.5°C vs. 34 ± 0.5°C and 43 ± 1.0°C vs. 57 ± 1.0°C, respectively). The water-holding capacity of the soil in core areas was greater than in nearby areas where lizards rarely spent time (60.91 ± 1.45% vs. 51 ± 1.45%; Figure 3a). These biophysical variables also differed between core areas and random sites at the WA study site (soil-water-holding capacity: F _1,233 = 12.91, P = 0.0004; mean temperature: F _1,108 = 27.80, P < 0.0001; maximum temperature: F _1,108 = 105.78, P < 0.0001; minimum temperature: F _1,108 = 12.82, P = 0.0005; Figure 3b). In core patches, mean and maximum temperatures were lower than at nearby random sites (23 ± 0.5°C vs. 25 ± 0.5°C and 30 ± 1.0°C vs. 41 ± 1.0°C, respectively), whereas minimum temperatures were marginally lower in random sites (17 ± 0.5°C vs. 19 ± 0.5°C; Figure 3b). Minor differences in thermal variables between sites, although statistically significant, are of dubious importance, given that the magnitude of difference is close to the accuracy of the measuring devices. The soil in core areas of lizard home ranges had a greater capacity to retain moisture than did that from adjacent areas that were rarely used by the lizards (66.66 ± 0.99% vs. 61 ± 1.14%; Figure 3b).

Figure 3

Comparisons between biophysical factors and main vegetation structures of core areas and non-preferred areas used by bluetongue lizards within tropical Australia. (a) NT study site, (b) WA study site.

Discussion

The northwestern tropics of Australia experience variable and challenging climatic conditions. For most of the year the area is hot (daytime maximum air temperatures often >40°C), with little or no rain during the extended (April to November) dry season. Any diurnally active ectotherm thus faces a high risk of lethal overheating if it remains in sun-exposed habitats; and may also be under hydric stress during prolonged dry periods [30, 43–45]. Presumably reflecting those challenges, a high proportion of large reptiles and amphibians living in this region are nocturnally rather than diurnally active [29, 46, 47]. Although bluetongue lizards are active by day [42, 48, 49], our radio-tracking studies have shown that the lizards escape lethally high temperatures in exposed sites by spending much of their time inactive within shaded retreat sites. Those sites are non-randomly distributed across the landscape, with the lizards spending on average 95% of their time in small core areas that constitute <40% of the home range [32].

Analyses in this paper show that those core areas possess a distinctive suite of attributes; basically, they are oases of cool moist sheltered conditions within the otherwise inhospitable landscape. At both study areas, most of those oases consisted of small concentrations of relatively dense vegetation; often, a few trees with a dense shrub or grass understory beneath them, creating a fairly thick ground cover of leaf litter in which the lizards could readily conceal themselves. At the highly disturbed WA site, some of the oases were created by artificial materials (concrete, wood, and so on) rather than by live trees. Despite this difference, both natural and artificially created core areas enabled the lizards to shelter in similar abiotic conditions (cool, moist, protected).

The causal bases for this heterogeneity in shelter opportunities are clear for the artificial shelters (which were constructed in places that people live or work), but obscure for the natural landscape. Small differences in elevation, and thus drainage patterns, as well as increased organic matter [50, 51] may have contributed to differences in soil type (and thus, moisture-holding ability), and these soil differences might have been exacerbated by the more profuse growth of vegetation in these small areas [11]. Stochastic events of tree germination and survival were probably also important; once a tree becomes established, its shade may allow the growth of numerous shrubs, and ultimately produce a small patch of moister, cooler conditions at ground level than can be found anywhere in the surrounding landscape.

It is likely that bluetongue lizards benefit from these conditions in a number of ways. Most obviously, they can avoid otherwise-lethal ground temperatures, and maintain water balance [31]. In addition, they are hidden from such predators as snakes, varanid lizards, raptors, and dingoes [52]. Although we did not quantify food availability, we observed more edible fruit, berries, invertebrates and small vertebrates (such as frogs) in the moist leaf litter than would be available elsewhere in the landscape. The differences in soil-water-holding capacity might also translate into a lower probability of inundation; after heavy rain, we often saw extensive pooling of water in areas between but not within core areas. Moister conditions of vegetation and soil within the patches may also provide refuges from fire, at least during low-intensity burns (for example, the early dry and wet seasons [53, 54]), particularly given the limited ability of bluetongues for rapid movement.

One of the most striking results of our study is the similarity in the attributes of core areas used by both species of bluetongue lizard, and by lizards at both sites. The two lizard species are broadly similar in overall size and morphology, but differ strongly in their geographic distributions and thus, in the range of abiotic conditions that they encounter over most of their extensive ranges. The northern bluetongue (T. s. intermedia) is restricted to the wet-dry tropics, primarily in open woodland habitat, whereas the centralian bluetongue (T. multifasciata) ranges widely across the arid interior, especially in spinifex grassland. The wet-dry tropics inhabited by northern bluetongues experience higher and more seasonally predictable rainfall than do the desert habitats of centralian bluetongues [55, 56]. Nonetheless, the broad abiotic challenges experienced by the two species are probably fairly similar: an open landscape that is lethally hot for long periods during daylight hours, with scattered patches of denser vegetation that provide shade, food, moisture, and protection from predators. At the site where the two species were sympatric (Keep River), the centralian bluetongues used spinifex clumps (rather than trees) for shelter more often than did the northern bluetongues (Figure 2b,c). This difference is consistent with the habitat types experienced by the two lizard species over the rest of their (allopatric) ranges.

The lizards’ frequent use of small isolated habitat patches, and the relative scarcity of such patches within the broader landscape, resulted in considerable overlap between individuals in patch use. For example, a single small (5 ha, or 5 × 10⁴ m²) patch of dense riparian habitat at the WA study site was used by at least 33 adult northern bluetongues over the course of our 9-month study at that site [32]. Such cohabitation, whether simultaneous or sequential, may substantially increase rates of parasite transmission among individuals [57–60], especially given that the cool, moist conditions in such sites would enhance survival of ectoparasites such as ticks, as well as the larvae of endoparasitic organisms [61]. These focal sites for reptile activity may also play a significant role in the social systems of the species involved, because many tiliquine species display complex sociality [62–69].

The lizards’ disproportionate use of small isolated patches of distinctive habitat also has strong implications for their conservation and management. Most obviously, these scattered patches of dense vegetation are critical for population persistence; destruction of such patches (for example, by feral pigs, or by cattle that use the shaded area as a diurnal refuge from sun exposure) would seriously reduce the carrying capacity of the broader landscape for these lizards. Future work could usefully examine issues of core-area connectedness in greater detail. Previous work has suggested that the visual perceptual range of a congeneric species (T. rugosa) is about 20 m, allowing lizards to execute direct movements towards refuge sites within this distance [70]. Thus, the perceptual range of a species is likely to play an important role in how lizards interact and move within the surrounding environment [71]. If habitat degradation reduces connectedness among core-area patches, such that many are isolated by distances too great for the lizards to detect from the nearest remaining patch, even otherwise-suitable remnant habitat patches may remain unused.

Another conservation implication of the lizards’ frequent usage of cool moist core areas is a negative one: the same patches will attract highly toxic cane toads (Rhinella marina), an invasive species currently spreading into this region [72–75]. Both bluetongues and toads are physiologically dependent on habitats that provide buffers against environmental extremes. Bluetongues are dependent on habitats that allow them to maintain the appropriate body temperature [23, 76, 77]. Similarly, cane toads have no morphological or physiological mechanisms to prevent evaporative water loss [78] and select shelter sites in relation to the protection they provide from desiccation [73–76, 79]. Toads at the invasion front are highly mobile, and spend daylight hours inactive in relatively cool moist areas [72]. Hence, toads congregate within refuge sites that are likely to be key foraging areas [74]. Bluetongues rely on relatively slow-moving prey, and are able to ingest very large items [80, 81]; thus, if encountered, a toad would be an easy target.

Conclusions

Cane toads pose a major threat to bluetongue populations, which have crashed in tropical Australia as the toad front has expanded westwards [82]. During this study, high bluetongue mortality was observed subsequent to the establishment of the first generation of these new invaders within both study sites (Price-Rees et al. unpublished data). Bluetongues from tropical Australia are highly sensitive to the toads’ distinctive toxins, and are likely to die if they attempt to consume even a single toad [82–84]. Thus, one unfortunate consequence of the lizards’ preference for moist cool diurnal retreat sites is an increased probability of encounter with a toxic toad. The kinds of habitat-use bias that we have documented in bluetongue lizards are probably widespread in tropical reptiles: in a landscape that is lethally hot for long periods, the optimal refuge habitat will be an area that is cool and moist. The scientific literature on reptile habitat use is dominated by studies of cool-climate species, where behavioral thermoregulation is facilitated by open, sun-exposed patches in a cooler moister habitat matrix [13]. However, most reptiles live in the tropics, where the challenge is typically to avoid overheating, rather than to attain high temperatures [29, 46, 85]. Thus, reliance on scattered patches of denser-than-usual vegetation that provide protection from predation, as well as from abiotic extremes, is probably a widespread feature of the ecology of squamate reptiles in many parts of the world.

Methods

Study site and species

Although most species from northwestern Australia are poorly known ecologically [86], this biota is under threat from several anthropogenic processes. Recently, the threats posed by overgrazing, agriculture, forestry, altered fire regimes [87, 88] have been exacerbated by the arrival of a highly toxic introduced species, the cane toad Rhinella marina [82]. Toad invasion has decimated bluetongue lizard populations in tropical areas of the NT, and may well have the same effect in WA [82, 83].

Our research was conducted at two sites within the Kimberley region of the wet-dry tropics. Ambient temperatures remain high all year round, but there is a pronounced seasonal variation in rainfall (90% of the 852 mm average annual rainfall falls during the four-month wet season, from mid-November to mid-March: Bureau of Meteorology 1968 to 2011 [89]). In the first year of our study we worked in Keep River National Park, in the far western NT near the border of WA (15°45′S, 129°6′E). Between November 2009 and March 2010, we tracked both Tiliqua multifasciata and T. scincoides intermedia. This area is largely unaffected by anthropogenic disturbance; the topography within this study area varied from sandstone plateaus and floodplains, in a matrix that is dominated by open eucalypt savanna woodland with an understory of native grasses [90]. The following year, we worked 6 km west of the town of Kununurra (15°46′S, 128°44′E), 50 km west of Keep River National Park. Between September 2010 and April 2011, we tracked the sole bluetongue species that occurs at this site (T. s. intermedia). The habitat was dominated by agricultural use, and included established gardens and built structures, roadside verges, and remnant riparian and swamp fragments.

Capture and GPS-tracking of animals

At the NT study site, we GPS-tracked nine T. s. intermedia (4 male, 5 female) and eight T. multifasciata (3 male, 5 female) for between two and 121 days each. At the WA study site, we GPS-tracked 32 T. s. intermedia (18 male, 14 female) for between three and 221 days each. Lizards were located by driving slowly along roads, and captured by hand. They were fitted with GPS transmitters using a vest-style backpack, and marked with a Trovan Unique™ ID-100B implantable transponder chip (Microchips Australia Pty Ltd., Keysborough, Victoria), and their sex determined by hemipene eversion. We also recorded morphology (mass, head width, snout to vent length, and length of hind leg, front leg, tail, and head). For details of capture methodology and transmitter attachment, see [32, 37].

Each lizard was fitted with a Sirtrack μGPS data logger with a custom-made attachment plate (mass, 24 g; only fitted to lizards >300 g; Sirtrack Ltd., Havelock North, New Zealand) and a single-stage VHF transmitter (battery life at 30 ppm, ~10 months; mass = 7 g) so that animals could be relocated and the GPS units retrieved if they became detached. Each GPS unit was programmed with four duty cycles that logged the location of individual lizards every 15 min during two periods of high activity (0600 to 1300 h; 1700 to 2300 h) and every hour during two periods of low activity (2300 to 0600 h; 1300 to 1700 h). Because of the high frequency of fixes (up to 63 per day), we re-caught lizards every 7 to 10 days and recharged the device (while it was still attached to the lizard) for up to two hours, after which lizards were released at their point of capture. At this time, GPS loggers were downloaded and the condition of the backpacks assessed; the backpack was replaced or repaired if necessary. Following release, lizards were visually relocated every 1 to 5 days using VHF telemetry. At each re-location, we recorded the lizard’s location, activity, date, time of day, and habitat.

Manufacturer specifications for Sirtrack GPS units claim an accuracy of <5 m 50% of the time, and <8 m 90% of the time. However, such factors as sedentary behavior, heavy canopy cover and complex terrain can decrease accuracy [91]. Thus, the horizontal dilution of precision (HDOP) recorded for every stored location fix was used to maintain the accuracy of the dataset. The HDOP is a standard measure of GPS accuracy, based on the number of satellites used for each fix and their relative position in the sky (Sirtrack GPS Receivers Manual, Sirtrack, Havelock North, NZ). Lower values indicate the likelihood of more accurate locations (HDOP values <2 = most accurate locations; HDOP values >10 = low accuracy). We discarded all locations with HDOP ≥5 from our analyses. In addition, regular visual relocations of lizards verified that locations where the lizards were recaptured were the same as those identified by the GPS units.

Surveys of habitat features in ‘core areas’

At the NT study site, we surveyed 106 ‘core areas’ based on the analyses of lizard movements (with one ‘random’ area surveyed for each core area, so N = 212 surveys). Of these sites, we also recorded temperatures at hourly intervals over a 24-h period in 61 core-area quadrats and 54 adjacent quadrats in areas less frequently used by lizards; soil samples were collected from all sites (106 core, 106 random). At the WA study site, 249 (141 core; 108 random) core-level habitat surveys were conducted. Temperatures were also recorded at hourly intervals over 24 h in 59 core-area quadrats and 50 random quadrats; overall, 236 (136 core; 100 random) soil samples were taken. Differing sample sizes reflect the impossibility of taking soil samples from some built environments, and malfunction or loss of iButtons. Overall, between 2,669 and 2,865 temperature measurements were recorded from core and random quadrats at each of the study sites.

Habitat features

We sampled attributes of the vegetation as well as biophysical variables (soil-water-holding capacity and air temperature). At both study sites, we collected two types of data on habitat parameters. First, we took habitat data whenever we located an animal (at the physical location of the lizard itself) and recorded whether or not the animal was active.

Our second type of analysis was based on a comparison between core areas used by lizards (as defined by GPS data on lizard movements) relative to the surrounding landscape. We did this by setting out a 1 m × l m quadrat within each identified core area, and a paired quadrat of the same size located 20 m from each core-area quadrat, in a random direction. To delineate core areas within each radio-tracked lizard’s home range, we used Ranges8 v2.7 [92, 93] to conduct a hierarchical incremental cluster analysis using nearest-neighbor distances. This analysis defines multinuclear core areas by generating a series of isopleths of diminishing percentage use (that is, 99%, 95%, 90%, 85%, …, 20%) and plotting polygons around the clusters (ICPs: incremental cluster polygons [92, 93]). This method links clusters of locations and excludes outlying locations by minimizing sums of nearest-neighbor distances as locations with longer distances are added, to define the core areas objectively [92, 93]. We used clusters that included 90% of locations to define the core-use areas (see Figure 4). For each lizard, up to ten of the highest density cores were selected for follow-up surveys; if a lizard had fewer than ten high-density patches, we surveyed all of the core areas that it had used. To locate survey sites, we used the easting and northing of a point within the center of the core cluster, and used the ‘go to’ function in a handheld GPS (model eTrex® H; Garmin, Olathe, Kansas, USA) to locate sites for our 1 m × l m quadrat habitat surveys.

Figure 4

Internal structure of the home range of an individual bluetongue lizard ( T. s. intermedia ) that was GPS-tracked for 64 days at the NT study site. Hierarchical incremental cluster analysis was used to identify core-use areas within the broader home range, with core areas defined by ICP 90% isopleths. Numbers show the intensely used core areas that were selected for our habitat surveys.

To determine random locations for paired quadrat surveys, we paced 20 m in a randomly allocated direction (obtained by numbering eight compass directions, and randomly drawing a number). If a random location coincided with a core area of the same or neighboring lizards, we discarded the site and randomly selected a different value for the direction. All sites for an individual lizard were surveyed on the same day, and all surveys were completed within a two-week period. At some locations within the WA study site, extensive overlap between conspecifics resulted in large blocks of frequently used areas [32]. Consequently, random sites could not be allocated without coinciding with other core areas. To overcome this problem, random sites were distributed outwards from the periphery of such dense-usage areas.

The spatial scale of this sampling (20 m from each lizard location to the paired ‘random’ site) was selected for three reasons. First, a site 20 m away is close enough for a lizard to detect using visual cues; the animals may be unable to detect sites further than this [70]. Second, most daily displacements by radio-tracked lizards were short (<20 m: see [42]), suggesting that animals are unlikely to move to an apparently suitable site if it is >20 m away from a lizard’s current location. Third, the broader landscape in which the study was conducted was homogeneous open savanna woodland, within which were located the small and distinctive patches of thicker forest in which lizards were usually found. The great spatial extent of the broader savanna habitat meant that samples at any spatial scale from 20 to at least 100 m would almost inevitably fall within that homogeneous matrix habitat – and thus, the choice of specific sampling distance (20 m) would have negligible influence on our results.

In each quadrat, we recorded the number of trees and shrubs within 5 m of the quadrat center, the percentage cover of canopy, understory, and leaf litter, and the percentage of bare ground and grass. We used categorical variables to describe broad macrohabitat types, understory microhabitat types, and the understory structure (see Table 4). Our classification systems at the two study sites differed slightly because of differences in habitat type (for example, the NT site was drier, and contained spinifex grassland not present at the WA site). At the center of each quadrat, we took a soil sample to determine soil-water-holding capacity (see below), and placed a Thermochron iButton® (model DS1921G-F5; Maxim Integrated Products, Sunnyvale, California, USA; accuracy: ± 1°C, temperature resolution: 0.5°C, operating temperature range: -40°C to +85°C) in every second quadrat (core and paired random site) to obtain 24-h temperature data. The iButtons were programmed to read temperature at 60-min intervals over the 24-hr period. If a lizard had five or fewer core areas, an iButton was placed in every quadrat. To quantify temperatures that lizards would experience, iButtons were placed in a location within each quadrat that best represented a potential refuge site. These placements were based on visual observation of the locations of telemetered lizards. The iButtons were placed on the ground within structures (for example, leaf litter, grass clumps). If ground-level structures were absent from the quadrat, iButtons were placed on the bare ground. For each lizard, we calculated mean, minimum, and maximum temperature for every quadrat sampled.

Table 4 Categories of broad macrohabitat types, understory microhabitat types, and ground-level structures recorded during point-level surveys

To assess the influence of soil characteristics, we quantified the water-holding capacity of soils from core areas and nearby ‘random’ areas using methodology modified from Warner et al. [19]. Both rainfall regime and soil types play a critical role in shaping tropical savanna ecosystems. The soil’s capacity to hold water can determine the composition of specific plant communities that can grow and persist throughout the dry season [7]. For example, during the wet season, soil is inundated with a surplus of water from monsoon rains but in the dry season porous soils with high clay content quickly lose this moisture, whereas loamier soils retain more water [7]. Thus, we measured the capacity of each soil sample to retain water, rather than the amount of moisture contained in a soil sample at the time of its collection. Soil samples from the center of each quadrat were placed in sealed plastic bags, and weighed prior to drying. To ensure that no moisture remained before testing, all samples were then dried, by placing the opened bags in a closed air-conditioned room maintained at 18°C for at least two weeks (until consecutive weighings revealed no further decrease in mass). Three subsamples of approximately 10 ml volume (uncompacted) were removed from each sample and placed into pieces of circular filter paper (90 mm diameter; Whatman International, Catalogue No. 1002090) stapled into a cone shape. The soil within the filter paper was weighed to the nearest 0.01 g, and then held upright in a rack of plastic vials. Once samples were arranged on the rack, 10 ml of tap water was added to the surface of each sample using a syringe. A pilot study determined that this volume exceeded the maximum capacity of any soil to retain (that is, no soil sample of 10 ml volume could retain 10 ml of water). Samples were reweighed after 45 min. The percentage water-holding capacity of the soil samples was calculated using the equation:

$$u=\frac{m_{\mathit{wet}}-m_{\mathit{dry}}}{m_{\mathit{dry}}}*100$$

Where m _wet is the saturated mass of the soil sample and m _dry is the initial dry mass of the soil sample before wetting.

Data analysis

Differences between habitat parameters as recorded at the point level (that is, individual locations of radio-tracked lizards) versus the core level (that is, habitat surveys) were evaluated with Pearson’s chi-square (χ²) test using JMP version 8.0 (SAS Institute, Cary, North Carolina, USA).

We conducted PCA on the habitat parameters (using JMP 8.0) to reduce the number of separate variables included in our analyses. Because the landscape types differed between the two study sites, we subjected vegetation structure variables from each site to separate PCAs. A nominal logistic regression analysis was used to identify effects of the most important PCs on lizard habitat selection (core areas vs. adjacent sites).

We used multivariate regression to assess the association between biophysical variables and the PCs.

To test whether biophysical variables differed between the core areas used by lizards versus less frequently used adjacent sites, we performed an analysis of variance with area type (core area vs. adjacent) as the factor and soil-water-holding capacity, and mean, minimum, and maximum temperature as dependent variables.