We sought to test the hypothesis that shrubs facilitated lizards by providing thermoregulatory opportunity. We predicted that lizards would associate with shrubs for a meaningful proportion of their daily activity cycle; that shrub association would increase in the afternoon when daytime temperatures peak (Filazzola, Sotomayor, Sotomayor, & Lortie, 2017 ); and that lizard association with shrubs would be correlated with thermoregulatory behaviors. The results of our study confirm the application of radiotelemetry to ecological facilitation studies and the application of such studies to the description of beneficial interactions between shrubs and vertebrate ectotherms.

Ecological facilitation theory provides a roadmap for describing and predicting the beneficial interactions of shrubs with other organisms within their communities (Bruno, Stachowicz, & Bertness, 2003 ; Bulleri, Bruno, Silliman, & Stachowicz, 2016 ; Filazzola & Lortie, 2014 ; Filazzola, Westphal, et al., 2017 ; Mcintire & Fajardo, 2014 ). Using facilitation theory, Filazzola, Westphal, et al. ( 2017 ) extended the exploration of the beneficial interactions between desert shrubs and vertebrates and found that one species of shrub provided facilitative benefits to a target species of lizard. We sought to confirm and add depth to their findings using radiotelemetry tracking of the same target species (Figure 1 ). Radio telemetry is a well‐tested and powerful tool that allows the longitudinal tracking of individual animals throughout their daily behavioral cycles (McGowan et al., 2017 ) and enables the direct observation of habitat interactions and behaviors. We used radiotelemetry to test and refine our understanding of the beneficial interaction of shrubs with lizards. To our knowledge, incorporating radiotelemetry into a facilitation study is a novel use of the method.

Deserts are highly distinct ecosystems that contribute significantly to global biodiversity and global ecosystem function. The conversion and loss of desert habitat is therefore a global biodiversity crisis requiring immediate intervention, including conservation of remaining undisturbed habitat and restoration of degraded desert (Bachelet, Ferschweiler, Sheehan, & Strittholt, 2016 ; Cook, 2004 ; Hannah, Carr, & Lankerani, 1995 ; Hoekstra, Boucher, Ricketts, & Roberts, 2005 ; Kéfi et al., 2007 ; Mouat & Lancaster, 2008 ; Westphal, Stewart, Tennant, Butterfield, & Sinervo, 2016 ). Identifying the drivers of ecological health in desert communities will be a crucial component of such interventions. Shrubs can maintain the diversity of desert plant communities (Flores & Jurado, 2003 ) and are predicted to play significant roles in the thermal ecology of desert ectotherms (Basson, Levy, Angilletta, & Clusella‐Trullas, 2017 ; Sears et al., 2016 ). Shrubs can also facilitate ectotherm populations in the face of climate change (Adolph, 1990 ; Kearney, Shine, & Porter, 2009 ; Sears & Angilletta, 2015 ; Sears et al., 2016 ; Sinervo et al., 2010 ).

Shrub density was calculated by visually counting individual shrubs within each lizard's MCP using aerial photographs (Google Earth, image taken December 20, 2016, accessed November 2017) and dividing that number by the area in square meters of the MCP. We calculated a standardized measure of shrub association zone area using on‐the‐ground measurements of a randomly chosen sample of shrubs in the study area ( n = 61), from which we calculated an average radius for each shrub following the method of Filazzola, Westphal, et al. ( 2017 ) and to which we added the 0.5 m association criterion described above. We calculated the area of each shrub association zone using the formula πr 2 and then took the average across the sample. We multiplied this standardized shrub association zone area by the number of shrubs counted in each MCP to obtain an estimate of the percent area of an MCP subsumed by shrub association zones.

Analyses were conducted in R (version 3.3.2). Habitat association was analyzed using a generalized linear model (Bolker et al., 2009 ) with the multcomp package (Hothorn, Bretz, & Westfall, 2008 ). Behavioral data were analyzed with a multinomial logistic regression using the nnet package that accounts for the multiple levels of nominal outcomes of the observations (Venables & Ripley, 2002 ). Home range size was calculated using a 95% minimum convex polygon (MCP) estimation (Mohr, 1947 ) using the adehabitatHR package (Calenge, 2006 ). MCPs were visualized in two dimensions in R.

Gambelia sila were relocated using a 3‐element Yagi antenna and Model R‐100 telemetry receiver (Communications Specialists, Inc., Orange, CA, USA). Once found, a location was taken for each lizard using a Garmin 64st GPS unit (Garmin Ltd., Olathe, KS, USA) and a laser rangefinder (Bushnell Outdoor Products, Overland Park, KS, USA). Habitat was categorized as whether a lizard was within 0.5 m of a shrub (shrub) or not (open) (henceforth, the “shrub association zone”), and behavior was scored from a suite of predetermined behavioral syndromes (Supporting Information Table S1 ). Disturbance from the observer to the lizard was kept to a minimum for each observation to avoid influencing behavior and habitat selection. At the completion of the study, all collars were removed from the lizards.

In the first 2 days following release, all captured G. sila individuals were relocated (i.e., repeatedly sighted using radio telemetry) several times to ensure that the lizards were successfully adjusting to the collars and that impacts to their behavior and survival were minimal. We looked for any negative effects the collar had on the lizards, such as impacts on movement or any other deviation from normal lizard behaviors. Gambelia sila were then formally surveyed for 24 consecutive days. Surveys were conducted on each lizard 3 times a day. Two of these daily surveys were conducted during daylight hours, when lizards were typically active above ground. One survey was conducted before noon, and one was conducted after noon. The third survey was conducted during the night when lizards are inactive below ground. The “night survey” was conducted before 7:30 a.m. or after 7:30 p.m. on each day.

Gambelia sila individuals were located during foot and vehicle surveys and captured using a pole and noose. Individuals were collared following the method of Germano and Rathbun, ( 2016 ). VHF radio transmitters (Holohil model BD‐2, frequency 151–152 MHz, battery life 8–16 weeks, Holohil Systems Ltd., Carp, ON, Canada) were attached to a small beaded chain collar using jewelry wire and epoxy, and the collars were then fastened around the lizard's neck. Gambelia sila were kept overnight to ensure the collar was fitted correctly and did not irritate or harm the animal and were then released at their capture site. Collars weighed 1.6–2.2 g (depending on the size of chain needed for the lizard's neck), and we ensured that the weight of the collar did not exceed between 5% and 10% of the body mass of the individual.

Ephedra californica , a basal gymnosperm in the Gnetophyta division, is a large, slow‐growing woody shrub restricted to arid environments in western North America (Sawyer, Keeler‐Wolf, & Evens, 2009 ). Although the genus has a worldwide distribution and is represented by over a dozen species in the desert southwest of North America, E. californica is the only species that occurs in the San Joaquin Valley, where it is locally considered rare and sensitive (Sawyer et al., 2009 ) and has been documented to be a foundation species in the San Joaquin Desert community (Hawbecker, 1951 ; Lortie, Filazzola, & Westphal, 2017 ; Lortie, Liczner, et al., 2017 ). Ephedra californica is the only large shrub represented at our study site. Gambelia sila is a state and federally listed endangered species endemic to the San Joaquin Valley and restricted to San Joaquin Desert habitat (Germano & Rathbun, 2016 ; Germano & Williams, 1992 ; Germano et al., 2011 ; U.S. Fish & Wildlife Service, 1998 ; Warrick, Kato, & Rose, 1998 ). Gambelia sila are diurnal and mainly insectivorous though they may eat smaller lizard species on occasion (Germano et al., 2007 ; Warrick et al., 1998 ). Though G. sila can bury themselves and will occasionally dig primitive burrows, they mostly utilize abandoned burrows of other animals such as D. ingens (Fields, Coffin, & Gosz, 1999 ; Prugh & Brashares, 2012 ). Adult G. sila are inactive in burrows for much of the year, emerging only from late March or April through July (Germano & Rathbun, 2016 ; U.S. Fish & Wildlife Service, 1998 ; Warrick et al., 1998 ). During the active season, G. sila will also spend the night underground in burrows and may return to a burrow during the day if the temperature becomes too hot or cold (Germano & Rathbun, 2016 ; Warrick et al., 1998 ).

The study was conducted on the Elkhorn Plain within Carrizo Plain National Monument (San Luis Obispo County, California, USA, 35.1914°N, 119.7929°W; Figure 2 ) within the San Joaquin Desert ecosystem (Germano et al., 2011 ). Average annual precipitation within the Monument ranges from 15 cm in the southeast to 25 cm in the northwest (Hijmans, Cameron, Parra, Jones, & Jarvis, 2005 ). The Elkhorn Plain is located within the Monument on an elevated plain separated from the main valley floor of the Carrizo Plain by the San Andreas Fault (Germano & Williams, 2005 ). The area has been heavily invaded by non‐native annual grasses including Bromus madritensis, Erodium cicutarium , and Hordeum murinum (Gurney, Prugh, & Brashares, 2015 ; Schiffman, 1994 ; Stout, Buck‐Diaz, Taylor, & Evens, 2014 ) but still provides habitat for endemic keystone species such as the giant kangaroo rat Dipodomys ingens (Bean et al., 2014 ). California jointfir, Ephedra californica was the dominant shrub at our study site. A much smaller woody perennial, Gutierrezia californica , can be found in some portions of the site at low frequency. The blunt‐nosed leopard lizard, Gambelia sila , was well documented on the study site (Germano, Smith, & Tabor, 2007 ).

Shrub use by individual lizards did not vary significantly as a function of shrub density within that lizard's home range (Figure 4 ). Percent of MCP areas subsumed by shrub association zones ranged from 1% to 15% with an average of 5% of total surface area, and frequency of shrub use by lizards was significantly higher than predicted by the percent of MCP area subsumed by shrubs ( Z = −4.714 from a Wilcoxon Signed ranks test, p < 0.001).

Behavior differed significantly between habitat types (Figure 3 , Table 2 , p < 0.0001). Lizards were observed cooling under shrubs significantly more than in the open (Figure 3 , Table 2 , p < 0.0001). Because simple presence under shrubs may not necessarily imply cooling, we used cues from individual lizard postures when scoring their behavior as “cooling” (Supporting Information Table S1 ). Lizards were also observed avoiding predators under shrubs more frequently than at other microhabitat types (Table 2 , p < 0.0001). The predators that lizards were observed avoiding in this study were all aerial predators (either ravens or raptor species). Predator interactions were all indirect and based on the observer's intuition; therefore, this result should be regarded as preliminary data. Burrowing and interacting occurred significantly less often under shrubs ( p < 0.0001). Other types of behavior such as sunning, hunting, or active observation did not differ significantly between habitat types. Observed behavior also differed significantly between different times of day, for example, lizards were more frequently observed sunning in the morning in both habitat types compared to the afternoon and more often burrowing and avoiding predators in the afternoon (Figure 2 , Table 2 , p < 0.001).

The frequency of lizard observations differed significantly between shrub and open (Figure 3 , Table 1 , p < 0.01). Shrub association frequencies of individual lizards ranged from 0 to 0.63 with a mean of 0.23 ± 0.035 SE (Supporting Information Table S2 ). Observations of lizards within open habitat did not differ between different times of day, but observations of lizards associating with shrubs differed significantly between morning and afternoon with lizards being found more frequently at shrubs in the afternoon (Table 1 , p = 0.0252).

On average the home ranges of the lizards overlapped with only two other individuals within a population throughout the entire sampling period (mode = 2 overlapping mcp polygons, one‐sample t test for µ = 2, t = −0.22535, df = 26, p ‐value = 0.8235), and there were no significant differences between the two genders in the extent of overlapping number of home ranges (GLM, family = poisson with total area per individual as covariate, χ 2 = 42.416, p = 0.39806). Our results were thus consistent with Tollestrup, (1983), Warrick et al., ( 1998 ) and Germano and Rathbun, ( 2016 ).

4 DISCUSSION

Shrubs are foundation species in many ecosystems due to the facilitative benefits that they provide to both plant and animal species (Filazzola & Lortie, 2014; Lortie, Filazzola, & Sotomayor, 2016). We hypothesized that E. californica facilitates G. sila by providing thermoregulatory benefits. Our finding that individual lizards were associated with shrubs more than predicted by shrub area within their home ranges as well as the observed significant shift toward shrub association in the afternoon during peak daytime temperatures supports our hypothesis. Our behavioral observations provided a suite of potential behaviors that lizards were likely to perform during the observation period; thus, our finding that cooling predominated under shrubs compared to other behaviors further supports our hypothesis. Our observation that shrub use was not correlated with shrub density suggests that lizards are actively choosing shrubs over open habitats rather than as a consequence of shrubs being more densely distributed in their home ranges. The observed association of G. sila with shrubs is consistent with results of studies of thermoregulatory behavior of lizards (Sears et al., 2016, Vickers, Manicom, & Schwarzkopf, 2011, Vickers & Schwarzkopf 2016, Basson et al., 2017, Grimm‐Seyfarth, Mihoub, & Henl 2017) and suggests that shrubs may facilitate G. sila by providing shade. Although facilitation as an ecological process does not necessarily include lifetime fitness as a component (Stachowicz, 2001, Bruno et al., 2003, Michalet et al., 2011, Michalet & Pugnaire, 2016), our inference that E. californica facilitates G. sila would take on additional relevance, particularly to the potential for community structuring and the promotion of resilience in lizard populations, if the effects of E. californica facilitation on G. sila individual fitness were quantified. In the case of diurnal lizards, the link between thermal habitat quality and individual fitness has been firmly established by both theory and empirical testing (Kirchhof et al., 2018, Ortega, Mencía, & Pérez‐Mellado, 2016, Vickers et al., 2016, Pontes‐da‐Silva et al., 2018, Sinervo et al., 2018, Camacho et al., 2018), but we believe further study on the E. californica—G. sila relationship with respect to individual fitness is warranted.

Shrubs can buffer the extremes of multiple environmental conditions such as temperature, wind, and solar radiation, creating a moderate microclimate under their canopy (Kerr & Bull, 2004; Pugnaire, 2010). At the landscape scale, the presence of shrubs and their pattern of distribution (i.e., clumped vs. dispersed) will affect lizard thermoregulatory behavior and can be crucial to an ectotherm's thermoregulatory efficiency (Basson et al., 2017; Sears et al., 2016). Sources of shade are particularly important for ectotherms, which must maintain body temperature through behavior (Díaz & Cabezas‐Díaz, 2004; Huey, 1974; Huey & Slatkin, 1976; Kerr & Bull, 2004). Visual concealment from predators and physical protection is also important (Fields et al., 1999, Filazzola, Westphal, et al., 2017). Shrubs may therefore provide important mechanisms of facilitation for G. sila. Our results suggest an important mechanism (shrub restoration) for the management of desert ectotherms such as G. sila and provide support for radiotelemetry as a viable method for studying ecological facilitation.

Shrub use by G. sila was addressed in one previous paper that also used radiotelemetry. Germano and Rathbun, (2016) employed post hoc tests to answer the question of whether shrubs are important components of G. sila habitat. One test depended on an assumption based on Schoepf, Schmol, Keonig, Pillay, and Schradin (2015) that home ranges were resource‐based (i.e., shrub‐limited) and would thus be smaller in the presence of high‐quality habitat (=shrubs), while another test sought to bound the amount of shrub habitat present in lizard home ranges away from a null expectation. The authors found no effect of shrubs on home range size but did find more shrubs present within lizard home ranges than predicted. Our a priori approach (i.e., taking direct observations of association with shrubs) provided evidence that lizards actively seek out shrubs rather than randomly encountering them during their daily activity. The lack of a correlation between individual shrub use and shrub density suggests that a threshold presence of shrubs may be sufficient to provide thermoregulatory opportunity, therefore a strong correlation between absolute number of shrubs within a home range, and home range size would not be predicted. This conclusion is further supported by the Germano and Rathbun's (2016) result that home ranges tended to include more shrub habitat than predicted by the study‐site wide prediction, that is, it is likely beneficial that some shrubs be available within the daily activity theater of individual lizards. Our results therefore confirm and are consistent with the results from Germano and Rathbun, (2016).

Germano and Rathbun, (2016) also provide a caveat against overestimating the importance of shrubs to G. sila by noting that G. sila occurs in places that lack shrubs. Given the variation that we observed in lizard shrub association within one population is not surprising that entire populations can persist in relatively shrubless areas. It should be pointed out that the possibility that lizards in shrubless areas may be using alternative strategies to effectively thermoregulate (Germano & Rathbun, 2016 suggest rodent burrows may substitute for shrubs) the fact they may do so does not negate our findings that shrubs provide thermoregulatory benefits to lizards.

Although heritability of thermoregulatory response in lizards has been found to be low in species where it has been estimated (Logan et al., 2018; Paranjpe, Bastiaans, Patten, Cooper, & Sinervo, 2013), heritable variation in propensity to use shrubs could allow a population to adapt to the loss of shrubs at the landscape scale (presuming that shrubs were primordially present; Logan, Cox, & Calsbeek, 2014). However, where population‐scale variation exists in the predisposition to use shrubs, such as we found in this study, it would be reasonable to propose that shrubs be made available to those lizards that are predisposed to associate with shrubs. The net effect would be to optimize the habitat available for that population. Such optimization may be crucial to impart population resilience to climate change (Sears et al., 2016; Sinervo et al., 2010). Additionally, structured and/or heterogeneous habitats are becoming increasingly recognized as important to achieve individual‐scale thermoregulatory optimization for lizards (Basson et al., 2017; Clusella‐Trullas & Chown, 2014; Goller, Goller, & French, 2014; Sears et al., 2016).