Conservation of skinks in New Zealand cities

Urbanisation dramatically alters ecosystems, disrupting key processes and threatening species persistence. Despite growing knowledge of the impacts of cities on wildlife, particularly for birds and mammals, there remains a paucity of research on how urbanisation impacts other taxa, including reptiles. We examined how urbanisation affects skink communities across multiple New Zealand cities by undertaking pitfall trapping surveys in a range of habitats. Of the 20 endemic species thought to have been present historically in the cities, only four were captured during surveys, Oligosoma aeneum in Hamilton; O. polychroma , O. aeneum , and O. ornatum in Wellington; O. polychroma in Nelson; and O. aff. polychroma Clade 5 in Dunedin. Our results reveal highly variable patterns of occupancy, abundance, and complex associations with key threats and habitat modifications. Sites in Nelson and Wellington were much more likely to record skinks than Hamilton and Dunedin, and modelling showed abundance of O. polychroma was positively related to rat tracking rates when grass cover was low but showed a negative relationship when grass cover was high. O. polychroma body condition was negatively associated with the proportion of urban land cover (built up areas or transport infrastructure). The low number of captures of several species is concerning and long-term monitoring is needed to assess the trajectories of these vulnerable populations. Effective conservation management of urban wildlife requires knowledge of where populations exist and how urban processes affect persistence. For understudied urban-dwelling taxa such as reptiles, the lack of this information may be putting populations at risk.


Introduction
The effects of human activities on ecosystems are widereaching, long-lasting and seldom more intense than in cities where the majority of the world's population now lives (McKinney 2002;United Nations 2018).Landscape-scale changes inherent to urbanisation dramatically decrease vegetated land cover, substantially reducing the availability of resources for wildlife (Valiela and Martinetto 2007;Grimm et al. 2008;Dupras and Alam 2015).Additionally, wildlife faces ongoing challenges of anthropogenic disturbance (e.g., traffic, development, etc.), and exotic competitors and predators (McKinney 2006).These urban processes tend to have strong negative effects on native biodiversity (McKinney 2006), and in the United States, urbanisation has been identified as one of the leading causes of species extinction (Czech et al. 2000).
As a result, over the last 20 years there have been numerous calls for conservation in cities (Miller and Hobbs 2002;Dunn et al. 2006), both to reduce negative impacts on the often significant biodiversity that can exist in cities (Ives et al. 2016;Soanes and Lentini 2019) and to ensure that humans continue to experience nature and interact with native species in their daily lives (Miller 2005;Dearborn and Kark 2009).However, conservation in cities requires knowledge of what species are present, how populations are faring and how the urban environment impacts different species (French et al. 2018).Research into urban wildlife has typically been biased towards mammals and birds, with reptiles, amphibians, fish, and arthropods being under-represented (Magle et al. 2012).This is concerning as responses to urbanisation differ considerably among taxa (McKinney 2008).
Globally, reptile populations are in decline and one contributing factor is land cover change resulting from urbanisation (Cox et al. 2022).The few studies examining reptiles in cities illustrate the highly variable nature of their responses to urbanisation -likely due to the highly variable nature of cities, with differing degrees of disturbance and high heterogeneity, particularly at moderate levels of urbanisation (French et al. 2018).In Tucson, Arizona, while reptile abundance and diversity generally decreased across a rural to urban gradient, it peaked at intermediate residential housing density (Germaine and Wakeling 2001).Similar patterns have been observed in other taxa, possibly resulting from the wide range of land uses that take place in the transition zone between a city's centre and outskirts (McKinney 2006(McKinney , 2008)).It is also inevitable, however, that some 'urban adaptor' species are better adapted than others to cope with the challenges of urban life (McKinney 2002;Kark et al. 2007).For example, the eastern blue-tongued skink (Tiliqua scincoides) in suburban Sydney, Australia, has likely been successful in these habitats because it exhibits strong site fidelity, readily makes use of artificial refuges, and has relatively r-selected reproductive traits (i.e., early maturation and large litters; Koenig et al. 2001).
In New Zealand, where ecosystems evolved in the absence of nearly all terrestrial mammals, threats to native reptiles in cities may differ from those in most other parts of the world.Like other modified landscapes across the country, New Zealand's cities retain only a fraction of the species diversity present in their regions before human colonisation (Woolley et al. 2019).Although the declines of now extirpated species were not documented, patterns of species loss in cities appear similar to those across much of the New Zealand mainland and likely resulted from similar threats: anthropogenic habitat loss and the introduction of mammalian predators (Tingley et al. 2013;Woolley et al. 2019).While the most vulnerable species are probably already lost from cities, of those that remain, seven of the nine skink species and five of the seven gecko species are classified as 'Declining', 'Nationally Vulnerable', or 'Relict' under the New Zealand Threat Classification System (Woolley et al. 2019;Hitchmough et al. 2021).
Despite the challenges of urban living, cities do offer conservation opportunities for lizards (Woolley et al. 2019).Green spaces such as parks and urban bush reserves are a mainstay for biodiversity, including lizard populations (Melzer and Bell 2014;Bell 2019).In addition to reserves, cities also offer valuable habitat in backyards and 'informal green spaces', such as the edges of transport infrastructure or urban parks.Although the value of these areas for biodiversity has not been well studied (Rupprecht et al. 2015), some of the only published research about lizards in New Zealand cities comes from these habitats (Barwick 1959;van Heezik and Ludwig 2012;Bell et al. 2018).In general, however, there has been little effort invested in the study of urban lizard populations (Woolley et al. 2019).Few biodiversity surveys for lizards have taken place in cities (but see Melzer and Bell 2014;Bell 2019), and much of the knowledge of species occurrence is anecdotal, often from salvage and biosecurity operations.As a result, there is little information about where lizard populations exist in cities, how populations are faring, and whether species might be vulnerable.
Here we illustrate the importance of this information for urban conservation management using New Zealand skinks as a case study.We conducted skink surveys in four New Zealand cities to assess the state of their populations in an urban context.Specifically, we aimed: (1) to compare the current species diversity in New Zealand cities with that thought to be present historically, (2) to compare skink abundance among cities and between urban and non-urban habitats, and (3) to explore relationships between population parameters and variation in the urban environment.

Study sites
Twelve study sites were selected in each of four New Zealand cities: Hamilton, Wellington, Nelson and Dunedin (n = 47 in total, as one Dunedin backyard was withdrawn from the study; Fig. 1).To ensure sites were representative of the diversity of habitats available in urban environments, in each city the sites were evenly divided among three strata: bush reserve, amenity area and backyard.Within strata, site characteristics differed widely, however, strata did have internal similarities and defining features.Bush reserves comprised areas of remnant or replanted native woody vegetation that were often actively managed and valued as 'green spaces' (Fig. 2a).Amenity areas were less densely vegetated and usually comprised more grassland Traps were 4 L plastic buckets dug into the ground so that their openings were flush with the substrate.A lid was pegged 10-20 mm above each trap to provide shade.Traps had 6-8 holes (approx.4 mm in diameter) drilled in the bottom to allow water drainage, and contained a layer of soil and foliage (approx.30 mm deep) to provide cover, a moistened sponge (40 mm × 60 mm × 10 mm) to prevent dehydration, and a square of wire mesh (100 mm 2 with ~ 10 mm mesh) curved over the floor of the trap to help prevent predation by rodents (Hare 2012).Traps were baited with canned pear (approx. 2 cm 3 ) and checked daily for 7-10 days.Two trapping sessions (early and late summer, i.e., November -December and February -March) took place in each of 2017-18 and 2018-19 (four sessions in total).Trapping was done in all four sessions in Wellington, but in only three sessions in Hamilton, Nelson and Dunedin (see S1 for details of trapping effort).
Captured animals were identified to species level, weighed and measured (snout-vent length, total length, length of tail regeneration).In Wellington and Nelson, skinks were made individually identifiable by clipping one toe on each foot (natural toe loss was incorporated into combinations).Skinks were not marked in Hamilton and Dunedin as prior knowledge suggested that the number of or scrubland, compared with taller vegetation in reserves.Amenity areas were often located at the margins of sports fields or transport infrastructure (Fig. 2b).Backyards consisted of privately-owned residential gardens, with participants recruited by door-knocking.Backyards were highly variable: lawns and cultivated exotic species were common components, though many had native plantings.The degree of management ranged from untidy (Fig. 2c) to highly maintained (Fig. 2d).

Skink surveys
Pitfall trapping was used to sample skinks at each site, as all skink species thought to be present in the cities are groundactive (rather than arboreal) and thus able to be captured in ground-based traps.In bush reserves and amenity areas, trapping grids consisted of 25 traps in a 5 × 5 arrangement with 2 m spacing between traps as recommended by other studies (Lettink et al. 2011;Wilson et al. 2017).In backyards, grids of either 6 or 10 traps were used, depending on space available.These grids were installed in 2 × 3 or 2 × 5 arrangements comparable with the other grids or, where space did not allow, traps were distributed around the garden at distances greater than 2 m.At each site, a line of ten stations of monitoring equipment was deployed at 50 m spacing.Each station consisted of a plastic footprint tracking tunnel (Black Trakka, NZ) with pre-inked card (Gotchatraps, NZ) and a plastic chew card pre-baited with peanut-butter-flavoured possum dough (traps.co.nz).Tracking tunnels were pegged into the ground and baited at both ends of the tunnel with rabbit paste (Erayz).Chew cards were folded and nailed to a tree or other suitable object approximately 300 mm from the ground and 2-5 m from the tunnel.Both tracking tunnels and chew cards were deployed for 6 nights.
We analysed data from spring (when surveys were conducted in all cities) using a subset of the dataset focussing on three key predators of lizards: rats (Rattus rattus, R. norvegicus), mice (Mus musculus), and European hedgehogs (Erinaceus europaeus).To make the mammal indices more relevant to the immediate area of the trapping grid, only the five closest stations were used in analyses.The maximum distance between any grid of pitfall traps and the furthest mammal station used in analyses was c. 200 m.
As tracking tunnel and chew card indices were correlated for the three mammal types (rats, r = 0.63, p < 0.001; mice, r = 0.35, p < 0.044; and hedgehogs, r = 0.37 p = 0.032), only the tracking tunnel indices were used in analyses as these indices were more sensitive to detections of the targeted species (especially hedgehogs; Balls 2019).

Statistical analyses
Site characteristics (Table 1) were compared using one-way anova, and where significant effects were found, pairwise comparisons were made using Tukey's HSD post-hoc tests.
For each site, an index of abundance for each skink species was made by calculating catch per unit effort (CPUE; number of skinks captured per 100 trap days) across all sessions.
At sites with grids of 25 traps where O. polychroma was recaptured, we estimated population density during each of three trapping sessions, using spatially explicit capturerecapture (SECR) models (Efford and Fewster 2013).Populations were assumed to be closed for the duration of each session.We fitted spatial detection models to the capture data, combining data from all sites and trapping sessions, by maximising full likelihood with the 'secr' package v. 3.2.0(Efford 2019) in R v. 3.5.1 (R Core Team 2018).SECR uses mark-recapture data to estimate density along with two spatial detection parameters: g 0 , the probability of capture at a trap located at the centre of an animal's home range, and σ, the width of the detection function, which was assumed to be half-normal.σ describes the spatial scale over which captures would be insufficient to perform mark-recapture analyses.Tail loss and natural toe loss were recorded, and pregnant females were identified visually by the shape of their abdomen.

Habitat-and landscape-scale covariates
Environmental covariates were collected to characterise the study sites, both at the scale of the local habitat and the wider landscape.A convex densiometer (Model A, Forest Suppliers) was used to measure summer canopy cover at each site by recording the proportion of mirrored squares in which the sky was obscured by foliage.Measurements were taken at heights of 1.2 m at the four cardinal compass points at each corner of the 25-trap grids or both ends of the smaller backyard grids.These measurements from each site were then averaged to calculate a mean canopy cover.Proportion of grass cover within pitfall grids was estimated visually on a 1-5 scale (1 = 0-19%, 2 = 20-39%, 3 = 40-59%, 4 = 60-79%, 5 = 80-100%).
Hourly temperature at the substrate was recorded at each site using iButton data loggers (DS1921G-F5#, Maxim Integrated) covered in 400 mm 2 squares of Parafilm (Parafilm M, Bemis Company Inc.) to increase weather resistance.Two data loggers were deployed at each site (n = 94) in November-December 2018 and were collected in January-February 2019, yielding a period of 34 days when all loggers were in the field at the same time (19/12/18-21/1/19).Data loggers were tethered with nylon monofilament to traps at opposite corners of the grid, positioned 200 mm from the trap, and oriented north.They were nestled into soil, vegetation or leaf litter to simulate realistic locations of skinks and to avoid direct sunlight.
Spatial analysis was used to calculate the proportion of urban land cover within a 500 m radius of each site.The Land Cover Database (LCDB v.4.1, https://lris.Scinfo.org.nz/layer/48423-lcdb-v41-land-cover-database-version-41mainland-new-zealand/, accessed: 12/3/19) was clipped within 500 m radii of each of the site locations using the buffer and clip tools in Arcmap v.10.5.1 (ESRI 2017).The resulting polygons were exported into their own layer and the dissolve tool used to calculate single values for the proportional area of each land cover type.Urban cover consisted of 'Built-up Area (settlement)', 'Urban Parkland/ Open Space' and 'Transport Infrastructure'.

Pest mammal surveys
Concurrent with lizard surveys, pest mammal surveys took place in each city.In Hamilton, Wellington, and Dunedin, mammal surveys were carried out in spring (late October -early December) 2017 and 2018, and autumn (late May for all non-pregnant individuals of a species (as used in other studies, e.g., Grayson et al. 2014;Batson et al. 2015).Linear mixed effects models were used to assess the effect of different environmental factors on abundance and body condition of the most abundant skink, O. polychroma, in Nelson and Wellington.Before running models, multicollinearity was assessed by testing for correlations among predictors using Spearman's rank correlation analyses.Where correlations were strong (|r s |>0.6;Hosmer and Lemeshow 2005), the predictor that correlated most strongly with other predictors was dropped from the model.CPUE and BCI were fitted as dependent variables in separate models.CPUE was log 10 -transformed to meet the assumptions of the model (normally distributed residuals, homogeneity of variance).Canopy cover, grass cover, mammal tracking indices (rat, hedgehog and mouse), substrate temperature and proportion of urban cover within 500 m of the sites were fixed effects, and site nested within city and session were random effects (Table 1).Interactions between grass cover and the three mammal tracking indices were tested but were removed from the model if non-significant.Model goodness-of-fit was determined by calculating marginal R 2 (variance explained by fixed effects only) and conditional R 2 (variance explained by both fixed and random effects; Nakagawa & Schielzeth 2012).All analyses were carried out using R version 3.5.1 (R Core Team 2018).

Site variation
The study sites varied substantially in their environmental characteristics within and among the cities and habitat types.The average substrate temperature during the 34-day summer period ranged from 24.6℃ at one backyard in Nelson to 12.9℃ in a bush reserve in Dunedin.Temperature varied significantly between the cities (F 3,43 = 43.35,p < 0.001); temperatures at sites in Dunedin were significantly lower than at all other cities (Tukey HSD, p < 0.001; Fig. 3a), and temperatures at sites in Nelson were warmer than those in Wellington (Tukey HSD, p < 0.001; Fig. 3a).Proportion of canopy cover also differed between cities (F 3,43 = 9.833, p < 0.001), with sites in Dunedin having more canopy cover than Nelson and Wellington (Tukey HSD, p < 0.001; Fig. 3b).Higher rat tracking indices were recorded at sites in Hamilton than all other cities (F 3,37 = 7.514; Tukey HSD, Dunedin and Nelson p < 0.01, Wellington p < 0.05; Fig. 3e).Hedgehog tracking rates also varied significantly between cities (F 3,37 = 11.92), with higher rates detected in Dunedin and Nelson compared with Hamilton and Wellington (Tukey HSD, p < 0.001 and p < 0.05 repectively; Fig. 3f).
an animal is detected and can be interpreted as implying a bivariate normal model for movement, where 2.45σ is the 95% home range radius (Wilson et al. 2007;Jennrich and Turner 1969).Skink home range centres were assumed to be distributed according to a Poisson point process, with the density parameter estimated from the capture data.Model selection was used to compare models with different covariates for the capture probability (g 0 ) term.The parameter g 0 was either constant (i.e., the null model) or with one of four behavioural responses: b, bk, B, Bk; where b is a permanent behavioural response to capture (either trap-happiness or trap-shyness), bk is a permanent behavioural response to capture in relation to a particular trap (i.e., trap-specific), B is a transient response of trap-happiness or -shyness only if it had been captured on the previous capture occasion, and Bk is a transient, trap-specific response (Efford 2019).For all models, density was modelled as a function of trapping session (a trapping period at a particular site) and σ was fitted as a constant.Models were compared using Akaike's Information Criterion adjusted for small sample sizes (AIC c ; Burnham and Anderson 2002).When more than one model was supported (ΔAIC c <8), we used modelaveraging (Burnham and Anderson 2002).
A body condition index (BCI) was calculated using the ratio of log 10 -transformed mass to log 10 -transformed SVL The city in which sampling took place (n = 2, random effect) Site The site at which sampling took place (n = 47, random effect) Session The session during which sampling took place (n = 4, random effect) of any species were for O. polychroma in Nelson and Wellington.In Nelson, mean CPUE across all sessions (± standard error of mean; SEM) was 6.47 ± 2.19 captures per 100 trap days, ranging from 0.17 to 28.00 (Fig. 4a).In Wellington, similarly high CPUEs were recorded at some sites (mean = 3.56 ± 1.2, range = 0.00-15.25).In Dunedin, CPUE of O. aff.polychroma Clade 5 at one coastal site (4.32 captures per 100 trap days) was comparable with those of O. polychroma in Wellington and Nelson (Fig. 4a).However, the species was only caught once at one other Dunedin site during the study (mean = 0.41 ± 0.3).In Wellington, mean CPUE of O. aeneum was 0.92 ± 0.74 captures per 100 trap days ranging from 0.00 to 9.03 and in Hamilton mean CPUE of O. aeneum was 0.12 ± 0.08 (range = 0.00-1.04)(Fig. 4b).
O. ornatum was captured only twice in the course of this study and only at one site, a backyard in Wellington (1.04 captures per 100 trap days).Mean CPUE for this species in Wellington was 0.09 ± 0.09.
No differences were detected between cities in terms of proportion of grass cover, proportion of urban land cover or mouse tracking rates (Fig. 3c, d and g).
The proportion of urban land cover did differ among the three habitat types (F 2,40 = 6.18, p < 0.01).Sites in residential backyards on average had 76% more urban cover within a 500 m radius compared with bush reserve sites.No other variables differed significantly among site types.

Skink species diversity and abundance in the four cities
Across the four cities, a total of five species of skinks were captured (four endemic and one introduced), though species diversity and site occupancy varied considerably among the cities with high spatial turnover (Table 2).
Catch per unit effort of the four endemic species varied considerably both by city and by site.The largest CPUEs We caught 148 individual O. polychroma (180 total captures) across all sites (all in Wellington and Nelson) in February 2018, 68 (80 total captures) in November 2018, and 190 (247 total captures) in February 2019.None of the four SECR models which modelled g 0 by the different behavioural responses were substantially better supported than each other, although all had more support than the null model (Table 3).Density estimates calculated by averaging these models ranged from 1.04 (95% CI 0.13-8.33) to 61.55 (33.29-113.82)skinks per 100 m 2 (Fig. 5).

Effect of environmental factors on Oligosoma polychroma
Only O. polychroma was captured in high enough numbers to allow the assessment of environmental factors on its populations.At sites where this species was captured, mouse tracking rate was negatively correlated with grass cover (Spearman's r s = -0.75,p < 0.001), and substrate temperature (r s = -0.47,p < 0.001) and positively corelated with rat tracking rate (r = 0.34, p = 0.004; S2).As these relationships were moderate to strong, mouse tracking rate was dropped from the models to prevent multicollinearity.All other correlations were either non-significant or had |r s | < 0.3 (S2).In linear mixed effects models, although rat tracking rate and grass cover were both positively related to O. polychroma CPUE, these main effects interacted; rat tracking rate was positively related to O. polychroma CPUE at low proportions of grass cover, and negatively at high (Table 4; Fig. 6a).Proportion of urban cover had a negative relationship with O. polychroma BCI (Table 5; Fig. 6b).No other factors were significant in either model.A summary of CPUE and BCI is available in S3.

Diversity and abundance of skinks in New Zealand cities
Of the 20 skink species thought to have been present in the regions of our study cities prior to human colonisation (Woolley et al. 2019), only four (20%) were detected in this study.Furthermore, in only one city (Wellington) was more than one skink species detected.These results are not altogether unexpected but are highly concerning and illustrate the dramatic reduction in biodiversity that has taken place in New Zealand's urban environments.
The four urban-dwelling skink species share characteristics that may make them more adept at surviving in the altered habitats that cities present.Compared with other New Zealand lizards, they have relatively high reproductive early 1970s in suburban Wellington found that O. aeneum populations can reach densities greater than 2000 skinks per hectare in suitable habitat (Bell et al. 2018).We were unable to estimate density for O. aeneum in this study, and capture rates were very low at most sites where they were captured (mean = 2.02 captures per 100 trap days, range = 0.15-9.03;Fig. 4).These results suggest that present urban populations may be smaller than they were 50 years ago.However, longterm monitoring is necessary to ascertain if populations are declining.Recent skink surveys for Wellington City Council recorded the species at fewer sites than expected across the city and described this low site occupancy as a cause for concern (Melzer and Bell 2014;Bell 2019).O. ornatum was captured only in a single backyard in Wellington and was not observed in Hamilton (despite its presence being recorded subsequently in a reserve in Hamilton in February 2020; pers.comm.M. Nelson-Tunley, Waikato Regional Council, 2020).This result too is a cause for concern.Other recent skink surveys of parks and reserves in Wellington have identified very few individuals of this species, and at few sites (Melzer and Bell 2014;Bell 2019).In contrast, the species is frequently captured in pitfall trapping surveys at the predator-exclusion-fenced Zealandia ecosanctuary (Nelson et al. 2016), suggesting that predation by introduced mammals may be an important factor driving these patterns.
outputs (c.five offspring per year) and fast maturation (less than 2 years to sexual maturity; Cree and Hare 2016; summarised in Table 6).They are also all habitat generalists and exist in a wide variety of environments (van Winkel et al. 2018).Behavioural and life history traits such as these have been identified in other urban-dwelling reptiles and it is possible that these traits allow these species to survive in urban environments where others cannot (Koenig et al. 2001).Notwithstanding these similarities, capture rates and site occupancy of the four species varied considerably and while some may be secure in cities, it is possible that others are vulnerable.
The highest site occupancy and capture rates were of Oligosoma polychroma in both Wellington and Nelson, where their estimated densities ranged widely from a few hundred per hectare (e.g., Reserve 2, Wellington) to several thousand per hectare (e.g., Reserve 2, Nelson; Table 7).In general, these densities were comparable with those of other studies of O. polychroma and the related O. aff.polychroma Clade 5 in a range of habitats with varying predator control regimes (Table 7).
Conversely, the low site occupancy and capture rates of Oligosoma aeneum and O. ornatum are cause for concern.O. aeneum was captured at only three sites in each of Hamilton and Wellington, and generally only in low numbers.Markrecapture studies of this species are rare, but one study in the Our surveys failed to detect other species thought to be present in cities (Table 2).In particular, O. zelandicum was not captured at any site in Wellington despite being regularly recorded in surveys in the fenced Zealandia ecosanctuary (Romijn 2013).This species has also not been captured in other recent surveys of Wellington reserves (Melzer and Bell 2014;Bell 2019).It has a Threat Status Classification of 'Declining' (Hitchmough et al. 2021).

Impact of urban environments on skink populations
We found no clear evidence for a negative impact of rats or hedgehogs on CPUE of O. polychroma.While many mammalian predators consume lizards (Newman 1994;Gillies and Clout 2003;Jones et al. 2013) and the impact of rodents on some lizard populations has been inferred from patterns of loss on New Zealand's mainland and population recovery on islands following pest mammal eradication (Towns 1991;Towns and Daugherty 1994;Norbury et al. 2014), evidence for the negative impacts of pest mammals on O. polychroma populations is mixed.As discussed in Wilson et al. (2017; and adapted in Table 7), there is no clear pattern indicating that populations of this species respond numerically to predator control.While locations with predator control have recorded some of the highest densities of O. aff.polychroma Clade 5 (e.g., Eglinton Valley; Table 7), some populations of O. polychroma achieved similarly high densities in the absence of predator control (e.g., Kelburn Cemetery and Pukerua Bay; Table 7).
The effect of rat tracking rate on CPUE was, however, dependent upon the proportion of grass cover, with a positive relationship when grass cover was low and a negative relationship when high (Table 4; Fig. 6a).This result suggests that grass cover may benefit populations of O. polychroma when rat abundance is low but not when it is high.It is possible that, in cities, grassland may act as an ecological trap for skinks, being preferentially selected as habitat but not providing adequate protection from rodent predation (Rotem et al. 2013;Battin 2004).Interactions between vegetation, herbivores and predators can result in complex and sometimes unexpected impacts for lizard populations (Norbury 2001;Hoare et al. 2007;Knox et al. 2012;Norbury et al. 2013).Densities of jewelled geckos (Naultinus gemmeus) in Coprosma spp.shrubland on the Otago Peninsula were four times higher in grazed sites than ungrazed, likely due to the higher activity of rodents at ungrazed sites (Knox et al. 2012).
Body condition of O. polychroma was negatively related to urban land cover (Table 5), in both Wellington and Nelson (Fig. 6b).Possibly, as the proportion of urban land cover

Table 3
Spatially explicit capture-recapture (SECR) models for O. polychroma during the three trapping sessions used for both Nelson and Wellington.Models were run with all sites (n = 13) and sessions (up to 3) combined (i.e., each session of trapping at a site was treated as a separate 'session' in the model.In alternative models, the parameter g   (Magle et al. 2012).Yet, in order to conserve these taxa in cities, ecological knowledge about their urban populations is imperative.Indeed, the present lack of research into how the urban environment affects these groups may be putting some species at risk.
Our study provides some of this essential information for New Zealand skinks and, by doing so, establishes a platform for their management in cities.The future of some lizard species in these environments is far from secure and it seems likely that population monitoring and conservation action will be required to ensure their longevity.
While the range of species captured in our sampling was generally consistent with existing knowledge of lizards in New Zealand cities, the failure to detect O. ornatum in Hamilton in a reserve where it has subsequently been found, increases, patch size of suitable habitat decreases making habitat and food resources more limited.Several studies document a decline in invertebrate abundance and diversity with increasing urbanisation (Niemelä et al. 2002;Ishitani et al. 2003;McKinney 2008).Alternatively, increased levels of disturbance may cause changes in foraging behaviours.Disturbances could be in the form of human activity (e.g., green space maintenance) or increased predator activity in more urban environments (e.g., domestic cats, Felis catus; Woolley and Hartley 2019).For example, several studies report differences in flight initiation distances and sprint speeds between lizards living in modified urban habitats and those in more natural environments (Prosser et al. 2006;Chejanovski et al. 2017).

Conservation of urban lizards -future directions
Urban wildlife research is especially challenging when target species have small home ranges or patchy distributions due to the difficulty of sampling across a highly heterogeneous environment.estimates are averaged from the four supported spatially explicit capture-recapture (SECR) models (ΔAICc < 3.2; Table 3).All models assumed that the parameter σ was constant across sites.Error bars are 95% confidence intervals surveyed in this study, using the full gamut of lizard sampling tools (Hare 2012a(Hare , 2012b(Hare , 2012c;;Lettink 2012).Due to the cryptic, secretive nature of many of New Zealand's lizards, it remains possible that species thought to have been lost from the mainland could yet be discovered (as was the case for Hoplodactylus duvaucelii at Maungatautari Mountain Sanctuary; Morgan-Richards et al. 2016 ).For this reason, managers should be open minded in their use of survey techniques.However, due to the heterogeneity of the urban landscape and costs associated with formal lizard surveys, sampling all potential habitats at the necessary intensity is likely not possible.Backyards, though highly variable in the habitat they offer, often constitute large proportions of urban green space in a city (e.g., 36% in Dunedin; Mathieu et al. 2007).Furthermore, the backyards surveyed in this study were among the most important sites for lizards in Wellington City, supporting high diversity (three species in one backyard), high abundance (highest CPUE for O. aeneum), and the only records of O. ornatum.While more formal surveys of backyards are recommended, these may be supplemented using citizen science.Records of lizard observations by members of the public could provide valuable information about species' distributions in backyards and other habitats that may allow the targeting of more intensive surveys.Such illustrates a key limitation of this study: that detection of some species, particularly those that are secretive, at low densities, or have patchy distributions, requires sustained effort and a range of detection tools.Further surveys in cities are needed to identify the locations of populations of skink species for which records are scarce (in particular, O. ornatum, O. aeneum and O. zelandicum).Additionally, future surveys in cities should target the full range of species that may be present, including geckos and arboreal skinks not  a citizen science project may also serve to increase public awareness for lizards and threats to them in cities (Ives et al. 2018).
Identifying the presence and location of lizard populations in cities, however, is only the first step towards their management and protection.It is possible that populations of some species are in unnoticed decline.To provide the information needed to identify population trends, as well as to demonstrate the efficacy of strategies proposed to benefit lizard populations (e.g., predator control, including cat management and habitat manipulation; Woolley et al. 2019), there is a need to establish long-term monitoring sites where surveys are repeated at regular intervals.Such long-term studies are rare in the urban context (Magle et al. 2012) but may be critical to ensure the persistence of some species in cities, particularly those that are cryptic and secretive.
Finally, it is possible that some of New Zealand's urbandwelling lizards may be well adapted to survive in urban environments, in particular O. polychroma.Internationally, other studies have identified apparently urban-adapted reptiles (Germaine and Wakeling 2001;Koenig et al. 2001) and in New Zealand, endemic birds vary in their resilience to introduced predators depending on their body size and degree of endemicity (Fea et al. 2020).Future research should identify what life history or behavioural traits facilitate the survival of these species in New Zealand cities. Understanding the mechanisms by which some species are able to survive in these highly modified environments will allow these species to be managed through further environmental change (Shine and Fitzgerald 1996;Koenig et al. 2001).

Fig. 2
Fig. 2 Examples of the three habitat strata from Wellington City.(a) a bush reserve site that receives community planting and weeding and is adjacent to a native forest; (b) an amenity area site of scrub/grassland

0
varied by one of four behavioural responses to capture: b = a permanent behavioural response, bk = a permanent trap-specific behavioural response, B = a transient behavioural response, and Bk = a transient trap-

Fig. 5
Fig. 5 Density estimates of O. polychroma at sites with 25-trap grids where more than one animal was captured in (a) Wellington and (b) Nelson.(c) shows estimates for Reserve 2 in Nelson with an extended y-axis to accommodate the high estimate in February 2019.Density

Fig. 6
Fig. 6 (a) Relationship between rat tracking index and catch per unit effort of O. polychroma at sites of high and low grass cover.Points are jittered horizontally for clarity.(b) Body condition index (BCI) of O.

Table 1
Summary of the variables used in linear mixed effects models

Table 2
Species detected in pitfall traps within the four cities in this study, along with other species recorded as present in the last 50 years and inferred to be present before human arrival (from Woolley et al. 2019 with species names and statuses updated using Hitchmough et al. 2021).All New Zealand skinks are in the genus Oligosoma except the introduced skink Lampropholis delicata.New Zealand threat classification status 2021 (NZTC) Hitchmough et al. 2021).* O. ornatum were recorded at a reserve in Hamilton in February 2020 (Pers.Comm.M. Tunley-Nelson).** O. kokowai were reintroduced to Zealandia Ecosanctuary, Wellington in 2016)

Table 4
Additionally, as cities have not Results of linear mixed effects model showing relationships between the model's fixed effects and catch per unit effort (CPUE) of O. polychroma.CPUE was log 10 transformed to meet the assumptions of normality.Effects with p < 0.05 are in bold.Marginal R 2 = 0.368 and conditional R 2 = 0.754

Table 5
Results of linear mixed effects model showing relationships between model's fixed effects and body condition index (BCI) of O. polychroma.Effects with p < 0.05 are in bold.Marginal R 2 = 0.231 and conditional R 2 = 0.546

Table 7
Summary of density estimates of O. polychroma and O. aff.polychroma Clade 5 from diverse studies throughout New Zealand, adapted from Wilson et al. (2017)