Spatial structure of woody cover affects habitat use patterns of ocelots in Texas

About 80% of the known breeding population of ocelots (Leopardus pardalis) in the USA occurs exclusively on private ranches in northern Willacy and Kenedy counties in South Texas. These private ranches support several large contiguous undisturbed patches of thornscrub, which is preferred by ocelots. Past studies have indicated ocelots in South Texas select for woody patches that contain extremely dense thornscrub (i.e., 95% canopy cover and 85% vertical cover) and require large patches of woody cover to survive. Landscape metrics have been used to explain ocelot habitat use in fragmented areas, but their application in less-fragmented rangelands is lacking. From 2011 to 2018, we used camera traps on the East Foundation’s El Sauz Ranch to assess seasonal habitat use of ocelots relative to landscape structure, configuration, and complexity and other site-level factors in South Texas. Seasonal habitat use and detection were positively influenced by larger mean patch area and lower landscape shape index values. We also observed ocelots were less likely to be detected during periods of drought and exhibited a seasonal trend in detection. Ocelots used woody patches that were larger and more regularly shaped, indicating a preference for areas with a lower degree of fragmentation across the study area. As patches become larger, they will coalesce over time and form larger woody aggregates, which will promote ocelot habitat use. Brush management needs to be strategic as patch area and shape index are a limiting factor to promote ocelot habitat use on working rangelands in South Texas. These results demonstrate the ability to use landscape metrics to discern the effects of spatial structure of vegetation communities relative to ocelot occupancy parameters.

Assessing landscape-level changes in landscape structure and determining the effects of historic or current fragmentation can aid in defining conservation strategies needed for recovering carnivore populations (Jackson et al. 2005;Zemanova et al. 2017). Throughout much of the world, human disturbance such as oil and gas exploration, deforestation, large-and small-scale agriculture and livestock production, and urbanization has had extensive effects on mammalian carnivores and other species (Sanderson et al. 2002;McKinney 2008). Unlike protected tracts of land such as national parks and wildlife refuges, private lands (i.e., rangelands, farms, forest plantations) may not be immune to changes in landscape structure, composition, or complexity and could become an impediment to wildlife use over time (Lombardi et al. 2018). Ninety percent of threatened and endangered species exist on private lands in the USA (Scott et al. 2001).
Among mammals, medium-large carnivores face the highest potential extinction risks due to higher trophic levels, lower population densities, and large home ranges (Cardillo et al. 2004). Within the family Felidae, small-medium-sized felids (e.g., ocelot (Leopardus pardalis), margay (Leopardus weidii), and marbled cat (Pardofelis marmorata)) are generally Communicated by: Krzysztof Schmid Electronic supplementary material The online version of this article (https://doi.org/10.1007/s13364-020-00501-2) contains supplementary material, which is available to authorized users. understudied compared to the big cats (Grassman et al. 2005;Nagy-Reis et al. 2017;Anile et al. 2019). Through the use of camera traps, researchers have now been able to peel back the curtain and gain a stronger understanding of how these lesser felids respond to natural and anthropogenic factors that may impact patterns of habitat use, population density, and occupancy (Jackson et al. 2005;Singh et al. 2014;Gómez-Ramírez et al. 2017;Anile et al. 2019). Statistical approaches like occupancy modeling have become a useful technique for discerning species occurrence or habitat use at different temporal and spatial scales (Wang et al. 2019). This statistical approach combined with robust long-term camera trap monitoring datasets, landscape metrics, and other site-level factors can aid in elucidating potential thresholds of habitat use and help guide recovery strategies, land stewardship, and management decisionmaking (Crooks 2002;Zemanova et al. 2017).
Ocelots are listed on Appendix I on the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) due to habitat loss and fragmentation, vehicle collisions, loss of genetic diversity, killing due to livestock depredation, and heavy hunting in the early-mid 1900s (Haines et al. 2005a;Janečka et al. 2011Janečka et al. , 2014Janečka et al. , 2016IUCN 2020). Ocelots are listed as Endangered in the USA and Vulnerable in Brazil, Colombia, and Argentina (Hunter 2019;IUCN 2020).
Ocelots inhabit a variety of vegetation communities including but not limited to Pine (Pinus spp.)-oak (Quercus spp.) and oak woodlands, semi-arid thornscrub, subtropical and tropical broadleaf forests, and gallery forests (Horne et al. 2009;Gomez-Ramirez et al. 2017;Paolino et al. 2018;Satter et al. 2019;Wang et al. 2019). Ocelots have been documented in open grasslands and flood plains, but studies have linked highest indices of ocelot occurrence to areas of dense canopy cover (Emmons 1988;Shindle and Tewes 1998;Harveson et al. 2004;Horne et al. 2009;Wang et al. 2019). Ocelots are solitary, semi-arboreal, and nocturnal predators that prey on small mammals, birds, and reptiles across their range (Booth-Binczik et al. 2013;Hunter 2019). Ocelot ecology has increasingly become more understood over the last 35 years with the advent of camera trapping, especially in Central and South America (Massara et al. 2015(Massara et al. , 2018Paolino et al. 2018;Satter et al. 2019;Wang et al. 2019); however, there is a paucity in similar data in the USA (Haines et al. 2006;Lombardi et al. 2020a).
In the USA, < 80 ocelots exist in two isolated breeding populations in Kenedy, Willacy, and Cameron counties in southern Texas (Haines et al. 2006;Janečka et al. 2014Janečka et al. , 2016Leonard et al. 2020). The larger ranch population exists on private working lands in Kenedy and northern Willacy counties and is home to about 80% of ocelots in Texas. These working lands also have the largest patches of thornscrub and live oak (Quercus virginiana) forest remaining in the region (Navarro-López 1985;Tewes 2017;Leonard et al. 2020;Lombardi et al. 2020b). The smaller Refuge population occurs in and around the Laguna Atascosa National Wildlife Refuge in eastern Cameron County (Haines et al. 2005a, b;Tewes 2017).
The majority of ocelot habitat studies in Texas to date have occurred on Laguna Atascosa National Wildlife Refuge, an area surrounded by a matrix of agriculture, highspeed roadways, and encroaching urban areas (Harveson et al. 2004;Haines et al. 2005a, b;Laack et al. 2005;Jackson et al. 2005;Horne et al. 2009). Jackson et al. (2005) reported that ocelots in these areas strongly select for areas that have greater shape index values, contain more edge, and are more isolated, indicating a greater degree of fragmentation in these areas. Unfortunately, however, this higher degree of fragmentation has contributed to severe genetic erosion over the last 60 years (Janečka et al. 2011(Janečka et al. , 2014(Janečka et al. , 2016. Though studies have been conducted examining ocelot habitat use patterns in fragmented landscapes, habitat use patterns in the larger working rangelands where there are likely larger patches and a greater degree of spatial heterogeneity in the structure, configuration, and complexity of woody cover are poorly known (Lombardi et al. 2020b).
We conducted this present study to provide valuable information regarding ocelot habitat relationships (e.g., effects of composition, configuration, and landscape complexity) within woody communities on working rangelands in South Texas. The results will help formulate recovery and conservation strategies for ocelots and guide habitat management plans. For this study, we monitored camera traps over 7 years to examine seasonal patterns of ocelot habitat use relative to landscape structure, configuration, and complexity and other site-level factors. Due to the absence of urban development in our study site, we expected ocelots to use large forested habitats on working rangelands differently than compared to other fragmented areas of South Texas. We propose the following hypotheses: (1) site-level metrics (e.g., increased in canopy cover, horizontal cover) will positively influence habitat use patterns of ocelots, (2) spatial structure of woody cover will positively influence what areas ocelots are likely to use, and (3) ocelots are less likely to be detected during periods of drought than during wet periods.

Camera surveys
From 1 May 2011 to 31 March 2018, we conducted camera surveys on the El Sauz Ranch as a part of a long-term ocelot population monitoring project. We used a systematic, gridbased sampling method with a randomized point within each grid cell (Meek et al. 2014;Lombardi et al. 2020a). We maintained a minimum of 1 km spacing between all sampling locations (i.e., camera stations), as required for ocelot camera surveys by the USA Fish and Wildlife Service (Permit Number permit TE822908-0). This spacing was based on historic ocelot telemetry data collected in the early 2000s (Lombardi et al. 2020a). We focused camera grids in areas of live oak and thornscrub forest in the northwestern (15 camera stations; 15 km 2 ) and southwestern (13 camera stations; 13 km 2 ) areas of El Sauz. These areas were selected due to ocelot occurrence being strongly tied to these vegetation communities (Navarro-López 1985; Shindle and Tewes 1998;Horne et al. 2009). At each randomized location, we established a camera station with two cameras in areas within or immediately adjacent to thornscrub or live oak patches. To avoid biasing our detection rates, we did not use bait or lures at the camera station to ensure ocelot behavior was not affected (Lombardi et al. 2020a). At each camera station, there were two Cuddeback® Expert Scouting Cameras and Cuddeback® X-Change Color cameras (Non-Typical Inc., Green Bay, Wisconsin, USA) with a 30-s delay between photographic events. Cameras were attached to trees or wooden stakes 0.5 m above ground and offset 1-2 m to enable individual identification of individuals for population monitoring (Satter et al. 2019).

Landscape structure and site-level metrics
Landscape and site-level variables were quantified to examine their influence on seasonal habitat use (Table 1). Site-level biotic variables examined included percent canopy cover, percent horizontal cover, canopy height (m), and shrub height (m). We marked the center of each camera station and estimated percent canopy cover using a Geographic Resource Solutions® (Geographic Resource Solutions, Arcata, CA) convex densitometer at 5 m in four cardinal directions and the center; the five values were then averaged for each station. We defined a categorical canopy cover variable based on open (< 25%), mixed (25-75%), and closed canopies (> 75%) at each camera station (Lombardi et al. 2020a). We used the point-intercept method to estimate mean horizontal cover at three different height profiles (0.0-0.5, 0.5-1.0, and 1.0-2.0 m) using a polyvinyl-chloride cover pole (2 × 3 m) marked at 0.5 m heights. We measured the height of each the lower shrub canopy and upper overhead canopy at each camera station at the center of each camera station (Horne et al. 2009).
To assess the influence of landscape structure, configuration, and complexity on ocelot habitat use, we conducted an unsupervised land cover classification of the El Sauz Ranch using 2014 1-m imagery from the National Agriculture Imagery Program in ERDAS IMAGINE (Hexagon Geospatial, Norcross, GA, USA. Following Lombardi et al. (2020a), land cover classes were based on four broad categories: woody cover (i.e., mesquite, huisache, thornscrub, and live oak forest, and woody areas within palustrine emergent wetlands), bare ground (i.e., inland dunes, barren ground, levees, paved and caliche roads), herbaceous (i.e., grasslands, Texas Gulf Coast Prairie, and herbaceous cover within palustrine emergent wetlands and along caliche roads), and water (i.e., lagunas, canals, and anthropogenic ponds). To assess our classification, we collected 629 ground-truth points in June and September 2016 using a Trimble® Geo 7 Series Handheld Computer with 1 m precision and Trimble Nomad® 1050 Series Handheld Computer with GBSS 1 m precision (Trimble Navigation, Ltd., Sunnyvale, CA, USA). We then used a confusion matrix to compare visually observed field data and the classified imagery (Congalton 1991) achieving an overall accuracy of 85.1%. Using the classified imagery, we placed a 500-m buffer around each camera site (sampling unit) and we examined six class-level metrics to quantify the spatial structure of woody and herbaceous cover types using

Occupancy modeling
We examined the influence of landscape structure, configuration, and complexity on detectability and multi-season habitat use patterns of ocelots using a simple multi-season occupancy model in program PRESENCE (v 12.1, MacKenzie et al. 2018). Because ocelot's home ranges in this region exceed that of the sampling unit (Leonard 2016), we will interpret true occupancy in this context as habitat use (MacKenzie et al. 2005(MacKenzie et al. , 2018. Before initial model construction, we used Pearson's correlation analysis to examine correlation(s) among site-level and landscape-level covariates. Non-correlated variables (|r| < 0.70) were grouped and combined based on 12 a priori hypotheses, which were based on the ecology of ocelots in Texas (Supplementary material Table 1). Continuous variables were normalized within the program PRESENCE (Murphy et al. 2017;Lombardi et al. 2018).
To build our habitat use models, we defined a 7-day (weekly) occasion, where the first ocelot detection within a given week was marked as detected that week. We defined seasons based on ocelot biology and climate criterion. Unlike bobcats (Lynx rufus), ocelots in Texas are a-seasonal breeders, having been documented giving birth to young kittens all year and young are not weaned until 18 months ). There are also no studies to support seasonal trends in mortality rates in this region. South Texas is also marked by a semiarid subtropical climate with summer temperatures often > 30°C and cooler periods of < 20°C. Due to the biology of ocelots (Hunter 2019;IUCN 2020) and the difference in seasonal temperatures, we partitioned seasons based on the monthly average temperature for the region (National Oceanic Atmospheric Administration: National Climate Database, https://www.ncdc.noaa.gov/cdo-web/datasets, accessed 17 March 2017) over the sampling period (hotseason [29 o C]: 1 May to 30 September; cool-season [18.4 o C]: 1 November to 31 March). We excluded the transitional months of April and October because regional average temperatures vary between hot and cool. We generated a binary capture history containing 14 seasons with 22 weekly occasions each season. The binary capture history was thoroughly inspected and was not found to violate the assumptions of closed populations or independence between occasions and between seasons (D. MacKenzie, Proteus, pers. com).
We used the fourth parameterization in the simple multiseason model framework using the seasonal occurrence and detection in PRESENCE. Habitat use (ψ) in this context was defined as the probability that ocelots will use a given site and movement in and out of that site is random (MacKenzie et al. 2018). Detection (ρ) is defined as the probability of detecting an ocelot at a site, at a given occasion or season (MacKenzie et al. 2018). For model selection, we used a two-step process following methodologies described by Noigere et al. (2013), in which we first modeled detection and then used the topperforming detection covariates to derive habitat use (ψ) models. Detection was modeled as a function of 15 biologically relevant models that included a drought variable (i.e., Palmer Drought Severity Index [drought: < − 1.0; wet > 1.0]) and two temporal variables (between-and withinseason variable), for a total of 30 total candidate detection models (Supplementary material Table 2). We used the standard temporal models available in program PRESENCE to examine whether detection changed between survey periods (between-season) or if detection of ocelots changed from occasion to occasion during each season (within-season model). Using the top-performing detection model, we modeled habitat use patterns (ψ) as a function of 11 biologically relevant candidate models based on our a priori hypotheses (Supplementary material Table 1). We assessed model fit by running 1000 bootstraps and did not consider models that failed to converge, as this may be due to overdispersion (b c > 1.0) (Murphy et al. 2017). Models with ΔAIC < 2.00 were considered to have substantial support (Burnham and Anderson 2002). Following the recommendation by Arnold (2010), we estimated 85% confidence intervals for coefficients in each top model.

Discussion
Landscape-level metrics are critical to the understanding of habitat use patterns for ocelots. We found that the mean patch area and landscape shape indices of woody patches were stronger predictors of ocelot habitat use and detection than canopy and horizontal cover on large working rangelands. The use of landscape metrics in occupancy modeling enables researchers to discern the effect of the spatial structure of vegetation communities and land cover and reflect levels of fragmentation on the landscape, which is important for examining ocelot population size and distribution. The research also highlights the importance of factoring in regional climatic conditions such as drought in occupancy models for ocelots and other mammals that occur in semi-arid areas. To date, only six studies have used occupancy modeling to examine various aspects of ocelot ecology (Massara et al. 2015(Massara et al. , 2018Paolino et al. 2018;Wang et al. 2019;Lombardi et al. 2020a). However, only this study and Wang et al. (2019) examine ocelot habitat use over 7 years. This study contributes new information on ocelot habitat relationships in the extreme northern edge of their geographic range and on working rangelands. From 2011 to 2018, habitat use of ocelots on El Sauz increased relative to increasing mean patch area and decreasing landscape shape index. Mean patch area is one of the strongest predictors of species occurrence or susceptibility for extirpation (Harrison 1991;Hanski et al. 1995;Crooks 2002;Michalski and Peres 2005). Bender et al. (1998) reported that forest interior species tend to decrease in population with smaller patch sizes due to habitat fragmentation. Janečka et al. (2011) reported that ocelots in South Texas do not respond well to fragmentation and often are unable to disperse between populations. Isolation of patches and subsequently reduced dispersal between the remaining populations has led to lower genetic diversity within the two ocelot populations in Texas. Habitat fragmentation can affect forests by not only decreasing size and area but also the shape and amount of edge available (Williams and Pearson 1997). As landscape shape index increases, patches become more irregularly shaped, and perimeter to area ratio increases, leading to a change in vegetation composition and structure (Arroyo-Rodríguez et al. 2008).
Our study is the opposite of Jackson et al. (2005), who found ocelots in the more isolated refuge population (2 9 km south of our study area) strongly selected areas with larger shape index values and inherently more edge, indicating a greater degree of fragmentation. As compared to our study area which contains the largest undisturbed patches of woody cover occurring in coastal South Texas, 95% of the native woody cover surrounding the Laguna Atascosa NWR has been lost since the 1930s and remaining unprotected patches are continued to be fragmented by rapid urbanization and rowcrop agriculture (Tremblay et al. 2005;Leslie 2016;Tewes 2017;Lombardi et al. 2020b).
We found that ocelots on El Sauz are more likely to use large areas of woody cover, which corroborates Leonard (2016) who observed ocelots strongly selecting larger patches of woody cover in proportion with open areas in South Texas. Ocelots have also been linked to increased forested cover in Peru (Emmons 1988), Argentina (Di Bitetti et al. 2006), Brazil (Wang et al. 2019), and Texas (Shindle and Tewes 1998;Harveson et al. 2004;Haines et al. 2006;Horne et al. 2009). Like these studies, larger patches of woody cover may be linked to higher availability of preferred prey (Booth-Binczik et al. 2013;Wang et al. 2019) or greater access to potential mates. As woody patches increase in size over time, adjacent patches will merge and coalesce over time and form larger patches, which may promote ocelot habitat use and potentially assist in dispersal and movement to unperturbed habitat patches.
We were unable to find strong evidence that canopy cover or horizontal cover affected ocelot habitat use. Ocelots have been strongly linked to dense vegetation communities (> 75% canopy cover) possibly for dense screening cover, ample prey, and secure den sites these vegetation communities may provide (Shindle and Tewes 1998;Harveson et al. 2004;Laack et al. 2005;Horne et al. 2009;Booth-Binczik et al. 2013). Heterogeneity in canopy cover (x = 56.7, 2.5-97%) and related foliage across our camera grids and areas used by ocelots may have limited our ability to discern these effects on ocelot occurrence during each survey season. Consequently, future research is needed to understand the effects these site-level metrics have on ocelots occurring in predominantly mixed canopy forests.
Drought may have negatively affected our ability to detect ocelots during this study. Although it is not clear why detection was higher during a severe drought (PDSI < − 4.0) than during moderate-low levels of drought (− 4.0 to − 1.0). Episodic drought is a common occurrence in the semi-arid climate of South Texas (Norwine and Kuruvilla 2007). This study began during one of the most severe droughts in Texas history and ended during another period of drought. Periods of drought also likely coincide with reductions in rodent populations, which ocelots and bobcats (Lynx rufus) in South Texas rely on (Tewes and Hornocker 2007). Potential reductions in rodent populations have been linked to bobcat home ranges in the region expanding 100%, which likely may also occur with ocelots (Blankenship 2000). Ocelots may be less detectable during periods of episodic drought due to increased daily movements or expansion of their home ranges to compensate for the reduction in rodent abundance (Tewes and Hornocker 2007). Future research may be warranted to understand the ecological effects of episodic drought on ocelots in South Texas.

Management implications
The contrasting results of this study compared to a previous study detailing ocelot habitat use in more fragmented areas Fig. 4 Influence of mean patch area (ha) and landscape shape index of woody cover on ocelot (Leopardus pardalis) detection the East Foundation's El Sauz Ranch in Willacy and Kenedy counties in South Texas underscores the importance of large private working rangelands for ocelot recovery and conservation in South Texas. Historic brush clearing 50 years ago on the southwestern edge of our study area and nearby working rangelands followed by subsequent woody cover encroachment has resulted in smaller patches of woody cover in some areas, which may not provide enough structural cover to support continued use by ocelots. However, due to the land stewardship ethic of these landowners, these patches will likely increase in size, thus increasing ocelot habitat use. Private working rangelands in South Texas will need to maintain large patches of woody cover to help facilitate ocelot habitat use and recovery in Texas.
Funding information The East Foundation, BF Vaughn III, Tim and Karen Hixon Foundation, Brown Foundation, Travis and Bettina Mathis, Anova LNG and their parent company Exelon Generation, and Wild Cat Conservation, Inc. provided financial support for this research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.