Vegetation cover and configuration drive reptile species distributions in a fragmented landscape

Keywords: fragmentation, habitat loss, landscape structure, MaxEnt, reptiles, species distribution modelling, squamate.

Habitat loss and fragmentation occur globally due to processes such as the expansion of agriculture and plantation forestry, and as a result, 79% of the world’s forested area is currently considered fragmented (Potapov et al. 2017). Landscape structure, defined as the composition and configuration of landscape elements (Cushman et al. 2008), influences species distributions, with many species persisting entirely in structurally modified landscapes (Bennett et al. 2006; Fletcher et al. 2018). A major challenge for ecologists and land managers is understanding (1) which components of landscape structure drive species’ responses, and (2) how landscape structure influences species distributions relative to other known drivers such as climate (Cushman et al. 2008).

Landscape modification and fragmentation can influence species in various ways, including increased risk of extinction (Crooks et al. 2017). Increases in the amount of habitat edge have been associated with increased rates of predation (Weldon and Haddad 2005; Hansen et al. 2019), and changes in land cover may influence road kills (Rincón-Aranguri et al. 2019). Further, changes in the spatial configuration and distance between habitat patches can alter connectivity and affect the capacity of species to move through the landscape and access key resources (Driscoll et al. 2013; Nimmo et al. 2019). In modified landscapes, species may alter their movement patterns to avoid unfavourable habitat (Hansen et al. 2020). This may in turn influence migration patterns and rates of recolonisation and dispersal, leading to reduced gene flow and consequently reduced genetic diversity (Neuwald and Templeton 2013; Herrmann et al. 2017). Genetic diversity is necessary for population persistence because it improves fitness and helps animals persist in the face of both stochastic events and shifts in environmental conditions (Hughes et al. 2008; Banks et al. 2013).

Although reptiles are one of the most diverse vertebrate groups, far less is known about their responses to landscape structure compared with birds and mammals (McGarigal and Cushman 2002; Böhm et al. 2013). This is concerning, because reptiles worldwide are considered particularly vulnerable to habitat loss and landscape fragmentation (Rodrigues 2005; Araújo et al. 2006; Cabrera-Guzmán and Reynoso 2012; Munguia-Vega et al. 2013). Indeed, a global meta-analysis of the impact of habitat modification on squamate reptiles found that modified habitat predominantly had a negative impact on species abundance (Doherty et al. 2020). Reptiles are known to respond differently to climatic variables compared with other vertebrate taxa, in part due to their thermoregulatory requirements (Powney et al. 2010; Buckley et al. 2012), and therefore may also respond differently to landscape structure. For example, in an urban environment jacky dragons (Amphibolurus muricatus) used open areas of lawn for basking (a thermoregulatory behaviour), but preferred to forage in nearby native vegetation that afforded them more cover (Burgin et al. 2011). Reptiles are primarily found in native vegetation and have a limited capacity for dispersal; therefore, loss of native vegetation cover is likely to influence their occurrence (Dubey and Shine 2010; Munguia-Vega et al. 2013).

Australia has the highest levels of reptile species diversity and endemism in the world, with around 1000 described species in 18 families (Webb et al. 2014). Since European colonisation around half of Australia’s forests have been significantly modified or cleared (Bradshaw 2012). Victoria is Australia’s most cleared state with around 34% of native vegetation remaining intact, primarily due to the growth of the wheat, sheep and cattle industries in the 1890s (Bradshaw 2012). The state is currently experiencing unprecedented population growth and is forecast to grow by at least 1.5% per year, reaching 8 million people by 2031 (Commissioner for Environmental Sustainability Victoria 2018a). This has led to record levels of housing development and urban land use, particularly at the fringes of the capital city of Melbourne, putting increased pressure on biodiversity (Commissioner for Environmental Sustainability Victoria 2018b). Fragmentation has already been implicated in the decline of Australian reptile species (Driscoll 2004; Brown et al. 2008), including the threatened striped legless lizard (Delma impar) (Scroggie et al. 2019) and the pygmy blue-tongue lizard (Tiliqua adelaidensis), which was considered extinct until it was rediscovered in 1992 (Armstrong et al. 1993).

Species distribution models are important tools for understanding relationships between environmental variables and species occurrence, and have been used widely in conservation planning and management (Araújo and Peterson 2012; Guisan et al. 2013; Guillera-Arroita et al. 2015). However, although fragmentation is commonly implicated in reptile species declines, most species distribution models are developed using only climate variables. In this study, our aim was to model the distributions of reptile species in Victoria using a combination of biophysical (i.e. climate, soils) and landscape structure variables. We quantified five landscape structure metrics, two representing habitat cover and three representing configuration. We modelled each species using the same set of predictor variables to allow us to compare the relative contribution of structural and biophysical predictors, and to identify the extent to which landscape structure influences reptile species distributions across structurally modified landscapes.

The study was conducted in Victoria, the southernmost state of mainland Australia (Fig. 1). We chose Victoria as our study system because in addition to a diverse reptile community it has a substantial gradient in both biophysical states and landscape structure; high levels of land clearing and urbanisation since European colonisation have resulted in marked landscape change. Victoria has a total land area of 227 496 km² and is Australia’s most densely populated state, with a total population of 6.3 million people (2017 figure) (Australian Bureau of Statistics 2018). Just over three-quarters (77%) of the population are concentrated in the state’s capital city of Melbourne (Australian Bureau of Statistics 2018). The total extent of native vegetation is estimated to have decreased by about half since European colonisation (i.e. pre 1750) (Eigenraam et al. 2013). Nevertheless, Victoria has a rich diversity of ecosystems including grasslands, wet and dry forest, coastal habitats, woodlands, and alpine and mallee communities (Cheal 2010). In most areas of the state the climate is temperate, with mild/warm-hot, dry summers (December to February) and cold, wet winters (Bureau of Meteorology 2016). The average minimum temperature for Victoria is 4.5°C in winter, with an average maximum of 27.4°C in summer. The annual average rainfall is 637 mm.

Fig. 1.  Map showing the distribution of the five land cover types (cropping, grazing pasture, native vegetation, plantation forest and other) across Victoria, Australia.

We obtained reptile species occurrence records from the Victorian Biodiversity Atlas, a database that collates data of species observations across Victoria (State Government of Victoria 2018). We included records between 2005 and 2017 to increase the congruence between the sample data and the date of spatial layers (i.e. 2016/2017) used to represent landscape structure. Turtles and sea snakes were excluded from our dataset, as were any records that did not identify individuals by both the genus and species names. Records were also removed if they included descriptions suggesting that an individual may have been sighted outside its range (for example, records that were described as ‘beach washed’ or ‘found in predator scat’). To improve independence and spatial accuracy, data were thinned to one occurrence point per 100-m cell, and records with an accuracy of >100 m were excluded. Finally, we excluded occurrence points within 5 km of the state boundary where landscape structure could not be evaluated accurately (see below), resulting in a dataset of 5749 species occurrence records. Records with additional naming information were renamed to only include the genus and species (for example, ‘Liopholis whitii PLAIN BACK MORPH’ was renamed ‘Liopholis whitii’, and ‘Eulamprus tympanum marnieae’ and ‘Eulamprus tympanum tympanum’ were renamed ‘Eulamprus tympanum’).

To account for spatial sampling bias, we generated background locations using the target group sampling method. This involved using the species occurrence locations described above as background data, operating under the assumption that the same sampling biases will apply to taxonomically related species sampled using similar methods as those used to observe the focal species (Merow et al. 2013). Because the 5749 occurrence records that we extracted were fewer than the recommended number of background points for MaxEnt (Elith et al. 2011), an additional 5000 points were generated within the core distribution of the occurrence records. We achieved this by using the 5749 data points to calculate a 50% kernel density polygon and then selecting the additional points at random from within this space. This resulted in a background sample of approximately 10 000 locations that was used as background in all models. Kernel density estimation was conducted using the adehabitatHR package (Calenge 2006) in the R statistical environment (R Core Team 2020).

Predictor variables were selected within two broad categories: landscape structure and biophysical variables.

A map of land cover types for Victoria was obtained from the DataVic website (State of Victoria (Agriculture Victoria) 2018). The original layer included 15 land cover types, which we reclassified into five classes: (1) cropping; (2) grazing pasture; (3) native vegetation; (4) plantation forests; and (5) other (e.g. horticulture and urban areas) (Fig. 1). The layer was clipped to the study area boundary and converted to a raster layer with a pixel size of 250 × 250 m using ArcGIS Desktop 10.5. A moving window analysis of the land cover layer was performed using FRAGSTATS (ver. 4.2; McGarigal et al. 2012) to calculate five landscape structure metrics: (1) Contagion Index; (2) Euclidean nearest neighbour (mean); (3) Native Vegetation Cover; (4) Shape Index (arithmetic mean); and (5) Shannon’s Diversity Index (Table 1) (McGarigal 2015). We used a square window with a side length of 5 km, and accepted samples with a maximum 10% of border/no data for all metrics. This analysis calculates a specified metric within a 5-km window for each pixel defined by the base land cover layer, resulting in a spatial layer of the specified metric with the same pixel size (250 × 250 m). In preliminary analysis we calculated a range of metrics, but chose the five listed above because they were effective at characterising different aspects of landscape structure and were minimally correlated (Pearson’s correlation: −0.2 < r < 0.6).

Table 1.  Description of landscape structure variables used in reptile species distribution modelling (Source: McGarigal 2015).

We selected climate variables known to influence reptile species distributions: annual mean temperature; annual precipitation; maximum temperature of the warmest month; and minimum temperature of coldest month. These variables have been successfully used to model distributions of Australian elapid snakes (Cabrelli et al. 2014) and skinks (Cabrelli and Hughes 2015). Bioclimatic data at a spatial resolution of ~1 km² was obtained from the WorldClim ver. 2 database (Fick and Hijmans 2017).

A spatial layer of soil surface texture for Victoria was obtained from Agriculture Victoria (Department of Jobs, Precincts and Resources). Soil surface texture is a classification system that estimates the relative amounts of sand, silt and clay particles, and influences some physical properties of soil (National Committee on Soil and Terrain 2009). The 19 categories in the original soil layer were reclassified into six broader soil categories (sands, sandy loam, loams, clay loams, light clays, medium clay), plus riverine (riverine areas where surface texture is not categorised) using information from Agriculture Victoria (State of Victoria (Agriculture Victoria) 2020).

We used MaxEnt to develop species distribution models in the R statistical environment (R Core Team 2020) using the dismo package (Hijmans et al. 2020). MaxEnt uses presence–background data to model species distributions, and performs well compared to other algorithms (Elith et al. 2006) and on datasets with a range of sample sizes (Wisz et al. 2008). We used MaxEnt to develop species distribution models for 40 reptile species with ≥30 occurrence records. We chose a minimum of 30 records as the cut off for analysis, because MaxEnt performs well for samples of this size or greater (Wisz et al. 2008).

Our aim was to compare the relative influence of different predictor variables on species distributions (rather than to identify the parsimonious model), so distribution models were generated for all 40 species using the same set of predictors, features and regularisation parameters. We used hinge features only to simplify the model while enabling the representation of non-linear responses, and chose the default regularisation value as a practical option for modelling multiple species (Merow et al. 2013). We selected the logistic output to allow for comparison of species responses to the different predictor variables (Merow et al. 2013). Area under the curve (AUC) values were used to assess the predictive ability of the models. AUC measures the area under the Receiver Operating Characteristic curve, and in the context of presence–background data represent the ability of a model to discriminate between species presences and background sample points. We followed the thresholds set out in Thuiller et al. (2003), where AUC values greater than 0.9 represent excellent predictive ability, 0.8–0.9 represent good predictive ability, 0.7–0.8 represent fair predictive ability, 0.6–0.7 represent poor predictive ability, and less than 0.6 represents models that perform little better than random. We used five-fold cross-validation to estimate the uncertainty associated with predictive ability, presenting the mean, minimum and maximum AUC values from the five runs. Maps of the relative likelihood of occurrence predicted by MaxEnt were produced for all 40 species in R using the raster (Hijmans 2020), tidyverse (Wickham et al. 2019) and cowplot (Wilke 2020) packages. Graphs representing species responses to individual predictors were produced in R using tidyverse. Relative likelihood of occurrence was represented on scale from 0 to 1 to allow for comparison between predictor variables and species.

Permutation importance values were used to assess the relative influence of the environmental variables to the species distribution models. Permutation importance is calculated by randomly permuting the values of each variable between the presence and background points, and measuring the resulting decrease in AUC. The bigger the decrease, the greater importance of that variable to the model. Permutation importance is considered a strong measure of variable influence because it does not depend on the pathway used to obtain the final model (Phillips 2017). We considered predictors influential if permutation importance was ≥5%.

Models were run using an average of 97 occurrence records per species (range 30–357). Models for 37 of the 40 species had an AUC score above 0.8, and 17 models had an AUC score above 0.9, indicating good to excellent predictive capacity (Table 2).

Table 2.  Results from modelling the distributions of reptile species across Victoria, Australia, showing the number of records (n) used in each model, model performance measured as area under the curve (AUC) and permutation importance of landscape structure and biophysical predictor variables.

The permutation importance of landscape structure variables was ≥5% for 80% (32 of 40) of species (Table 2, Fig. 2). Extent of native vegetation cover was the most influential of the five structure variables and resulted in the highest importance values for two skink species, Lerista bougainvillii and Tiliqua rugosa (29.4 and 22.8 respectively; Table 2). We selected 12 species with notable responses to landscape structure and display their relative likelihood of occurrence maps in Fig. 3. Despite the strong influence of landscape structure on the distributions of some species, climate variables were the strongest predictors of species distributions overall (Table 2, Fig. 2).

Fig. 2.  Box plot showing the range of permutation importance values for variables used in MaxEnt models to predict relative likelihood of occurrence of 40 reptiles across Victoria, Australia. Midlines of the boxes represent medians, and bottom and top of the boxes represent the 25th and 75th percentiles, with whiskers extending to 1.5 times the interquartile range. Blue circles represent raw data for biophysical variables and purple circles represent raw data for landscape variables. P, annual precipitation; TM warm, maximum temperature of warmest month; T, annual mean temperature; VEG, native vegetation cover; TM cold, minimum temperature of coldest month; SOIL, soil texture; SHA, Shape Index; SHDI, Shannon’s Diversity Index; CON, Contagion Index; ENN, Euclidean nearest neighbour.

Fig. 3.  MaxEnt predictions of relative likelihood of occurrence across Victoria, Australia, for 12 of the 30 reptile species with notable responses to landscape structure variables. Purple represents areas of high relative suitability; blue represents areas of low relative suitability.

Shape index was the most influential measure of landscape configuration, and influenced the distributions of seven species (permutation importance ≥5%). The permutation importance was higher for Pseudonaja textilis, Lampropholis guichenoti and Lerista punctatovittata (range 14.07–19.49) than the other species (range 5.03–5.57) (Table 2) (Fig. 3). The direction of responses to shape index varied among species; three species responded negatively and four species responded positively (Fig. 4). Contagion index and Euclidean nearest neighbour were the least influential configuration metrics (Fig. 2). Contagion index was an important variable for L. bougainvillii and Notechis scutatus (permutation importance 11.22 and 11.72 respectively) (Table 2) (Fig. 3). Both species responded negatively to contagion index, that is, relative likelihood of occurrence decreased with increasing aggregation of the landscape (Fig. 4). Only one species, D. impar, responded to Euclidean nearest neighbour; however, the permutation importance of this variable was low (5.22) (Table 2) (Fig. 3), with relative likelihood of occurrence increasing with distance between landscapes of the same type (Fig. 4).

Fig. 4.  Selected reptile species responses to landscape structure variables (native vegetation cover, Shape, Shannon’s Diversity Index, Contagion Index, and Euclidean nearest neighbour) predicted by MaxEnt models. Relative likelihood of occurrence values are presented on a scale from 0 to 1 to allow for comparison among species and predictor variables. Black lines represent the mean response predicted by the MaxEnt model; blue bands represent minimum and maximum values predicted from five-fold cross validation.

Shannon’s Diversity Index was an important predictor for three skinks, T. rugosa, Tiliqua scincoides and Eulamprus tympanum (permutation important 11.45, 10.80 and 6.37 respectively) (Table 2) (Fig. 3). The responses among species varied; T. rugosa responded positively to Shannon’s Diversity Index, T. scincoides had a hump-shaped response, and E. tympanum had a negative response (Fig. 4).

Climate variables were key predictors in all species distribution models. Annual precipitation was important for 37 species, and was the variable with the highest permutation importance for 17 species (ranging from 24.7 to 84.5) (Table 2). Maximum temperature of the warmest month was influential in 34 species models, and was the variable with the highest permutation importance for 15 species models (ranging between 81.55 and 24.57) (Table 2). Annual mean temperature and minimum temperature of the coldest month were both influential in 27 species’ models (Table 2). Soil texture was an important predictor variable for only eight species (Table 2); the soil texture permutation importance was highest for Liopholis whitii (11.54), which was more likely to occur on clay soils (i.e. Clay loams, Light clay and Medium clay).

Understanding the responses of animals to landscape structure is an important part of managing species in modified landscapes, but few studies have used landscape structure variables to model species distributions (but see Rodrigues and Lima-Ribeiro 2018; Cable et al. 2021; Senior et al. 2021). We modelled the distributions of 40 reptile species across a fragmented landscape in Victoria, Australia to compare the relative influence of landscape structure, climate and soil. We found that although climate was generally the strongest driver of species distributions, most species also respond to landscape structure. For L. bougainvillii and T. rugosa, the extent of native vegetation had the greatest influence on their modelled distribution.

Our results indicate that landscape structure plays an important role in driving the distributions of some reptile species in Victoria. The extent of native vegetation was the most influential landscape structure variable, and was important for species within all seven reptile families included in the study. Extent of native vegetation cover was an important driver of the distributions of three-quarters of the species in this study. Most of these species had a positive response to native vegetation cover, such as Diplodactylus vittatus, which appears to occur more commonly in regions with warm, dry climates and high vegetation cover. Likewise, Ctenophorus fordi was found to have a strong positive relationship with extent of native vegetation, which may be due to their association with hummock grasses (Triodia sp.) (Verdon et al. 2020). In semiarid parts of Australia, many reptiles, including C. fordi, are associated with the extent cover of hummock grasses, which are a dominant ground cover and an important resource for thermoregulation and protection from predators (Verdon et al. 2020). Habitat availability is critical to the persistence of most species, and habitat loss and modification are a clear driver of reptile declines in eastern Australia (Driscoll 2004; Brown et al. 2008) and worldwide (Doherty et al. 2020). Most reptiles, and particularly skinks (which represent most species in this study), are ground dwelling and therefore sensitive to loss and disturbance of ground and understorey layers including leaf litter (Brown 2001).

Interestingly, several species had negative responses to the extent of vegetation cover. For example, D. impar, a small, threatened legless lizard (~300 mm long), preferred areas with less native vegetation cover, potentially reflecting this species’ reliance on threatened native grassland habitat, rather than a preference for low levels of vegetation cover per se (Kutt 1993). Since European colonisation, grasslands have suffered extensive loss and degradation, and now remain primarily in isolated remnants (Scroggie et al. 2019). Disturbances such as fire and grazing may promote the persistence of D. impar when applied in isolation by promoting structural and floristic diversity of grasslands and increasing connectivity (Scroggie et al. 2019). However, when applied simultaneously, fire and grazing were associated with a high risk of local extinction, and any further loss of native vegetation cover is likely to exacerbate the decline of D. impar (Scroggie et al. 2019).

Shape index was the most influential configuration metric considered in this study, and reflects the perimeter-to-area ratio and therefore the amount of edge environment. Four species responded positively to shape index, including P. textilis which was more likely to occur in irregularly shaped habitat patches and where there was an intermediate level of native vegetation cover. More irregular, elongated patches have a higher perimeter-to-area and are therefore more exposed to surrounding land cover types, which may advantage species that engage in landscape complementation and require resources in multiple land cover types (Nimmo et al. 2019). Pseudonaja textilis are a widely distributed species known to exploit agricultural areas where prey species (specifically Mus musculus) are abundant (Whitaker and Shine 2003). Similarly, for two species (N. scutatus and L. bougainvillii), relative likelihood of occurrence was higher when the contagion index was low (representing low clumping and high intermixing of patch types), yet had hump-shaped responses to vegetation cover extent, such that occurrence was more likely when vegetation cover was intermediate. This suggests that although vegetation cover is important, these species may be able to take advantage of resources associated with increased intermixing of different land use types.

In contrast, several species, including L. punctatovittata, were more likely to occur in more regular-shaped habitat. Irregular and elongated patches may also have undesirable implications for species that inhabit them, including increased exposure to introduced plants and animals (Brown et al. 2006) and disease (Jellinek et al. 2014; Uriostegui-Velarde et al. 2018). Negative relationships have been found between exotic plant cover and both the abundance and richness of lizard species (Jellinek et al. 2004). Because weed cover is expected to be higher in irregular shaped habitat, such areas may represent less suitable habitat for reptiles. Additionally, more elongated patches have a higher edge-to-area ratio; edges are subject to more to extreme climatic conditions and more extreme variations in temperature, which may impact the ability of reptiles to thermoregulate (Jellinek et al. 2004). Shape index has also been found to negatively influence plant species richness, possibly due to the greater availability of resources in more regular-shaped patches (Li et al. 2020), which may in turn afford more resources for animals.

Shannon’s Diversity Index was important to the distributions of three reptile species, but the directions of relationships were variable. Eulamprus tympanum responded negatively to increasing diversity, and positively to native vegetation cover. This is probably at least partly due to this species’ preference for shelter sites with logs (Langkilde et al. 2003), which are likely to be more common in native vegetation than other land cover types. In general, the presence of other land cover types is thought to negatively impact reptile species persistence, and a global meta-analysis found most types of habitat modification (including mining, grazing, agriculture, plantations and patch size reduction) have an overall negative impact on reptile species abundance (Doherty et al. 2020). In south-eastern Australia, agriculture is thought to negatively affect reptiles (Driscoll 2004), and pine plantations have been shown to act as a barrier to connectivity in some species and increase connectivity for others (Mortelliti et al. 2015). In our study, T. rugosa, a habitat generalist, was more likely to occur in areas with moderate levels of diversity and native vegetation cover. This species is known to favour open habitat and commonly found on roadsides, and may utilise both natural and artificial shelter (Norval and Gardner 2020), characteristics that may provide it with the flexibility to exploit multiple land cover types.

Although our study shows that native vegetation cover and configuration influences the distributions of some species at a statewide scale, it is important to note other reptiles may respond to aspects of landscape structure at finer spatial scales. Habitat cover and configuration, as well as other aspects of landscape structure such as the edge-to-area ratio, influence microclimate and the availability of shade, which is required for thermoregulation. South-eastern Australia is also subject to low-frequency, high-intensity fires (Murphy et al. 2013) that create mosaics of vegetation age classes in the landscape. Where species select different age classes, fire may represent an important additional source of habitat fragmentation, which is likely to influence distributions and threaten species persistence (Sitters and Di Stefano 2020). Furthermore, it is anticipated following future changes to the climate, many species worldwide will alter their present distributions (Bezeng et al. 2018). However, such range expansions will not be possible if the landscape structure does not facilitate this movement. Incorporating landscape structure, as well as dynamic predictor variables such as fire (Senior et al. 2021; Swan et al. 2021), into models may improve our ability to predict future species distributions and enable more informed management strategies.

We found climate variables, especially annual precipitation and the maximum temperature of the warmest month, were the strongest drivers of species distributions for almost all reptile species. Multiple studies have shown climate is a key factor determining both the current and future distributions of reptiles (Ihlow et al. 2012; Cabrelli et al. 2014; Hosseinzadeh et al. 2020) and their nesting habitat (Pike 2013). Environmental temperature plays a particularly important role in regulating the metabolic rates and activity patterns of ectotherms (Buckley et al. 2012), and regulates the reproductive success and sex determination in many species (Witt et al. 2010; Ihlow et al. 2012; Castelli et al. 2021).

Models of both reptile and amphibian species distributions have predicted that, on average, species are likely to experience increases in their climatic range due to climate change (Araújo et al. 2006; Hosseinzadeh et al. 2020). However, both groups are generally considered to be poor dispersers, and this factor, combined with extensive fragmentation, means that we should not assume that species have the ability and pathways to disperse and colonise new habitat (Araújo et al. 2006). Indeed, modelling of European amphibians and reptiles has shown that when models assume no dispersal, virtually all species are projected to experience range contractions (Araújo et al. 2006). Therefore, it is important that other non-climatic factors like dispersal (Sahlean et al. 2014; Rodrigues and Lima-Ribeiro 2018) and changes to landscape structure (Uriostegui-Velarde et al. 2018) are incorporated into models if we are to properly understand the factors driving current and future species distributions.

Soil texture was an important driver of the distributions of a fifth of the species in this study, including common species such as T. rugosa and L. guichenoti, which both showed an aversion to clay loams. Clay soils experience seasonal changes, becoming waterlogged during the wet season and cracked during the dry season (Woinarski et al. 1999). Although clay soils may provide an important temporary water source for frogs (Ferreira et al. 2018), this seasonality may be unsuitable for many reptiles (Woinarski et al. 1999). Notable differences have also been reported among reptile communities on sandy and clay soils that support different plant communities in arid and semiarid parts of Australia (Driscoll et al. 2012). Amphisbaenians, which likely share traits with other burrowing reptiles with absent or reduced limbs, prefer sandy, loose soils and avoid compact, heavy soils – presumably due to the reduced energetic costs of moving underground in sandy soils (Martín et al. 2013). Interestingly, however, we found no relationships between sandy textured soils and either Lerista species, which are semifossorial and known for their ability to ‘swim’ through sand (Morinaga and Bergmann 2020).

Scale is an important factor in ecological studies because species respond to different drivers at different spatial scales (Cushman and McGarigal 2004). In the context of presence–background modelling, the area within which background locations are sampled (landscape extent) is an important consideration because it can influence model outputs (Elith et al. 2011; Merow et al. 2013). We chose the state of Victoria in south-eastern Australia as the landscape extent for all modelled species because it has a substantial gradient in both biophysical characteristics and landscape structure – the two groups of predictor variables of primary interest. Although using a common landscape extent to model the distributions of multiple species is a common approach (Guerin and Lowe 2013; Liu et al. 2013; Hageer et al. 2017; Swan et al. 2021), the importance of particular predictor variables is likely to be influenced by the range size of each species. For example, in our study, the relative influence of climate compared with landscape structure may be reduced for species whose range is larger than the landscape extent. Results should be interpreted with this in mind because the geographic range size of the species we modelled is variable.

Although native vegetation is generally considered important for reptiles, most studies that model species distributions at large spatial scales focus on species responses to climate variables. Nevertheless, landscape structure is expected to be important to reptiles because the composition and configuration of habitat influence access to resources including shelter, food and mates (Fischer and Lindenmayer 2007). Our findings show that the extent of native vegetation had an important influence on the distribution of many terrestrial reptile species in our data set, including endangered species such as D. impar, E. tympanum and V. varius. Therefore, management that focuses on minimising the removal of native vegetation and increasing revegetation is fundamental to maintaining species persistence. Our study also showed that the configuration of native vegetation was important for some reptiles, including common species such as L. guichenoti and P. textilis. This suggests there is potential to use landscape structure metrics such as configuration to help identify patches of native vegetation where habitat restoration and revegetation are likely to have the greatest conservation outcomes (Jellinek 2017). Strategic revegetation may be especially pertinent in regions such as Victoria, where increases in human population and urban development over the next three decades (Commissioner for Environmental Sustainability Victoria 2018a) are likely to have further negative consequences for biodiversity.

The data that support this study are available in FigShare at https://doi.org/10.26188/61dbdefb9427d.

The authors declare no conflicts of interest.

This research was supported by the Holsworth Wildlife Research Endowment and the Parks Victoria Applied Park Management Research Scholarship.