Apparent Competition, Lion Predation, and Managed Livestock Grazing: Can Conservation Value Be Enhanced?

Predator restorations often result in apparent competition, where co-occurring prey populations experience asymmetric predation pressure driven by predator preferences. In many rangeland ecosystems, livestock share the landscape with wildlife, including ungulates and the large carnivores that consume them. We examined whether apparent competition reorganized prey communities following restoration of lions (Panthera leo) to a savanna ecosystem, and whether and how livestock management could alter this indirect interaction between lions and their prey. Three lines of evidence supported the hypothesis that Jackson’s hartebeest (Alcelaphus bucelaphus lelwel; an ungulate of conservation concern) are suppressed via lion-mediated apparent competition. First, hartebeest exhibited an Allee effect where they were exposed to lions, but displayed negative density-dependent population growth where they were protected from lions. Second, spatial overlap between plains zebra (Equus burchelli; the primary prey of lions) and hartebeest further exacerbated lion predation on hartebeest. Finally, hartebeest were killed selectively by lions, whereas zebra were killed by lions in proportion to their abundance. We then tested whether glades (nutrient-rich hotspots created by abandoned cattle [Bos indicus] corrals) could be used to manipulate top-down control of hartebeest via their influence on the spatial distribution of zebra. Zebra aggregated at glades, and survival of hartebeest increased with increasing distance from glades, suggesting that corrals may be placed on the landscape away from hartebeest to create spatial refuges from lions. Our findings demonstrate how informed placement of livestock corrals can be used to manipulate the spatial distribution of primary prey (zebra), thereby reducing apparent competition suffered by hartebeest. Our work further provides an example of how integrating apparent competition theory with proactive livestock management can improve conservation efforts in multiple-use landscapes.

Predator restorations often result in apparent competition, where co-occurring prey populations experience asymmetric predation pressure driven by predator preferences. In many rangeland ecosystems, livestock share the landscape with wildlife, including ungulates and the large carnivores that consume them. We examined whether apparent competition reorganized prey communities following restoration of lions (Panthera leo) to a savanna ecosystem, and whether and how livestock management could alter this indirect interaction between lions and their prey. Three lines of evidence supported the hypothesis that Jackson's hartebeest (Alcelaphus bucelaphus lelwel; an ungulate of conservation concern) are suppressed via lion-mediated apparent competition. First, hartebeest exhibited an Allee effect where they were exposed to lions, but displayed negative density-dependent population growth where they were protected from lions. Second, spatial overlap between plains zebra (Equus burchelli; the primary prey of lions) and hartebeest further exacerbated lion predation on hartebeest. Finally, hartebeest were killed selectively by lions, whereas zebra were killed by lions in proportion to their abundance. We then tested whether glades [nutrient-rich hotspots created by abandoned cattle (Bos indicus) corrals] could be used to manipulate top-down control of hartebeest via their influence on the spatial distribution of zebra. Zebra aggregated at glades, and survival of hartebeest increased with increasing distance from glades, suggesting that corrals may be placed on the landscape away from hartebeest to create spatial refuges from lions. Our findings demonstrate how informed placement of livestock corrals can be used to manipulate the spatial distribution of primary prey (zebra), thereby reducing apparent competition suffered by hartebeest. Our work further provides an example of how integrating apparent competition theory with proactive livestock management can improve conservation efforts in multiple-use landscapes.

INTRODUCTION
The restoration of large carnivores can infuse optimism into conservation efforts that, historically, have been characterized by reactive, stopgap measures (Chapron et al., 2014;Wolf and Ripple, 2018). In the aftermath of such predator restoration, however, prey species often exhibit marked changes in numbers and behavior, such that composition and relative abundance of post-restoration prey communities may bear only slight resemblance to that of communities prior to predator extirpation (Lovari et al., 2009;DeCesare et al., 2010). Differences in carnivore-extirpated vs. carnivore-restored communities of prey often are attributed to apparent competition, in which preferences for prey by large carnivores lead to differences in the degree to which prey-some which are themselves conservation concerns-are suppressed via top-down control (Holt, 1977;Holt and Kotler, 1987;see DeCesare et al., 2010 for a review of cases in which apparent competition presents challenges for conservation; Holt and Bonsall, 2017).
Regardless of the mechanism by which they arise, prey preferences often result in secondary (less common) prey incurring greater risk of predation in proximity to primary prey. Secondary prey may decline while primary prey subsidize large carnivores, such that abundances of large carnivores and secondary prey are decoupled, driving the latter to rarity or even local extinction (e.g., Schmidt, 2004;Angulo et al., 2007;Hervieux et al., 2014;Serrouya et al., 2015;Chan et al., 2017). The likelihood that apparent competition generates these predatormediated Allee effects-the decrease in population growth at small population sizes, (Allee et al., 1949;Courchamp et al., 1999)-is greatest when large carnivores are restored after lengthy periods of extirpation, during which environmental conditions have changed, and prey abundance has increased (but see Berger et al., 2001;Ford and Goheen, 2015;Stier et al., 2016). Against the backdrop of shifted environmental conditions, effects of large carnivores can be stronger than expected, presenting a conservation challenge to species of secondary prey.
The factors that drive prey preferences-and thus predispose secondary prey to high rates of predation-include differential resilience of primary vs. secondary prey to predation, and differential space use between primary and secondary prey (DeCesare et al., 2010;Wittmer et al., 2013). This second feature of predator-prey dynamics can result in "gradients of consumption" (in which risk of predation varies monotonically and directionally; Orrock et al., 2008), creating refuges where encounters between large carnivores and their secondary prey are infrequent. In giving secondary prey a potential foothold for positive population growth, refuges are one of the few ways shown to negate predation-mediated Allee effects in nature, thereby providing a potential tool to ameliorate apparent competition (Sinclair et al., 1998).
African savannas hold promise and conservation importance for understanding how spatial refugia might be used to reduce apparent competition stemming from predator restoration. Almost invariably, communities of savanna ungulates are dominated by a single species that typically constitutes the primary prey for large carnivores, yet achieves sufficiently high abundances to escape top-down control. In contrast, populations of less abundant species (i.e., secondary prey) tend to be suppressed by large carnivores, thus creating strong potential for apparent competition (e.g., Sinclair, 1985;Harrington et al., 1999;Owen-Smith and Mason, 2005;Georgiadis et al., 2007a;Chirima et al., 2012). Ranching occurs alongside wildlife in many African savannas, and landscapes in these humanoccupied systems bear the imprint of livestock production in the form of glades: nutrient hotspots that attract wild ungulates and are derived from abandoned corrals or "bomas" (Augustine and McNaughton, 2006;Porensky and Veblen, 2015). Because livelihoods based purely on livestock production are becoming less reliably profitable, a changing mindset-to balance pastoralism with tourism, and potentially wildlife conservationis gaining traction in many areas (Prins et al., 2000;Odadi et al., 2011;Keesing et al., 2018). Viewing large predators consistently ranks as a top priority among tourists, potentially leading to financial benefits in ecotourism ventures (Lindsey et al., 2007). Restoring large predators along with diverse assemblages of wild ungulates may be key to economic viability in these regions (Cousins et al., 2008;Stein et al., 2010).
The goals of our study were 2-fold. First, we sought to test the hypothesis of apparent competition (Figure 1) in Laikipia County, Kenya. Here, multiple species of wild ungulatesmost notably Jackson's hartebeest (Alcelaphus buselaphus lelwel; hereafter simply "hartebeest")-have experienced recent declines following lion (Panthera leo) restoration in the late 1980s (Georgiadis et al., 2007b;Ng'weno et al., 2017). Restoration of lions has resulted from greater tolerance by ranch managers following decades of control via shooting and poisoning; however, declines in hartebeest populations are increasingly viewed as a particular conservation concern, leading some ranch managers to consider re-implementing lethal control of lions (Georgiadis et al., 2007b J. R. Goheen pers comm with Laikipia ranchers). Indeed, the impact of lion predation is sufficient to shift population growth of hartebeest from positive to negative (Ng'weno et al., 2017). In contrast to hartebeest and other declining species, populations of plains zebra (Equus quagga; hereafter simply "zebra")-the primary prey for lions in Laikipia-fluctuate in response to rainfall and density (Georgiadis et al., 2007a), but have not decreased with recovering lion numbers (Georgiadis, 2011;O'Brien et al., 2018). Consequently, we hypothesized that zebra populations suppressed those of hartebeest via apparent competition.
After assessing the hypothesis of apparent competition, we then tested whether livestock management could be used to manipulate the spatial distribution of predation risk to hartebeest, potentially enhancing coexistence of lions and their prey (and potentially removing a reason to lethally control lions). Under the hypothesis of apparent competition, we expected hartebeest in proximity to zebra to incur higher rates of predation. Because zebra (but not hartebeest) are attracted to glades (Veblen and Young, 2010;Augustine et al., 2011), we sought to quantify whether spatial separation between primary (zebra) and secondary (hartebeest) prey, driven by glade location, moderated apparent competition between zebra and hartebeest. Specifically, we predicted that (1) hartebeest occurring in areas FIGURE 1 | Hypothesized species interactions among lions, plains zebra, Jackson's hartebeest, and cattle in Laikipia, Kenya. (A) The hypothesis of apparent competition. Zebra populations are controlled bottom-up by rainfall, are the primary prey of lions, and support a growing population of lions since the cessation of lethal control in the late-1980s. Although they are less abundant than zebra, hartebeest are selectively killed by lions and are controlled top-down, such that hartebeest incur negative, indirect impacts from zebra via lions. (B) Cattle are corralled overnight in bomas where they deposit dung and urine for a period of ca. 1 month. (C) After abandonment, dung and urine break down over the course of 2-3 years, creating nutrient-rich forage in glades where zebra aggregate, but not hartebeest. This creates spatial separation between zebra and hartebeest, and potentially a refuge from predation if lions target hunting activity toward areas where zebra are abundant.
of high zebra density would incur greater risk of predation from lions; and (2) hartebeest survival would increase with increasing distance from glades. In the event that predation risk to hartebeest could be manipulated via spatial separation of hartebeest and zebra, this would provide a potential solution through which hartebeest numbers could be bolstered without lethal control of lions.

Estimation of Hartebeest Abundance and Population Growth
Between 2009 and 2015, we estimated hartebeest population growth and size within two zones under different management regimes: (1) a 294-km 2 area in which cattle production occurred alongside wildlife with the full complement of large carnivores (hereafter "control"); and (2) a 32-km 2 lion exclusion zone (hereafter "lion exclusion"), constructed with the intent of boosting numbers of declining ungulates, primarily hartebeest. The exclusion zone was demarcated from the adjacent control zone by a 2.5-m tall solar powered electrified (6,000-7,000volt) fence with nine strands spaced 0.2 m apart. The fence was fortified with chain-link 1.50 m above and 0.60 m beneath the ground, preventing lion incursion but remaining permeable to other large carnivores (Ng'weno et al., 2017). Stocking rates were maintained at equal densities of 20 cattle /km 2 in both zones. From 2009 to 2015, we conducted twice-monthly drive transects for hartebeest in both zones. During each sampling period, we systematically drove 17 (4-13 km) transects (12 in the control zone, 5 in the lion exclusion zone), based on a predetermined random starting point. We conducted surveys during 07:00-11:00 h, driving at a maximum speed of 10-15 km/h with two trained observers. At each hartebeest sighting, we recorded the group size and the distance and compass bearing to the group with binoculars and laser rangefinders following standard distance sampling methods (Buckland et al., 2015).
Hartebeest maintain small groups (Kingdon et al., 2013; x = 16 in this study), characterized by strong dominance hierarchies among females. Groups are relative sedentary and defend small (<5 km 2 ) territories against conspecifics; as a result, we were able to identify and track 179 unique individuals in 11 distinct groups over the latter 4 years of our study (2012-2015; see "Hartebeest survivorship and refugia analysis" below; Ng'weno et al., 2017).
Typically, apparent competition is characterized by one or more species of secondary prey that exhibit an Allee effect (Allee et al., 1949;Courchamp et al., 1999;Stephens et al., 1999). To test for an Allee effect, we calculated the instantaneous rates of increase (r) for each zone (control and lion exclusion) as r = where N is the estimate of population size from the ith survey at time t (Sinclair et al., 1998). For hartebeest populations in both zones, we related population size to r, and determined goodness-of-fit of curves from regression analysis (Zar, 2010).

Prey Selectivity of Lions
From 2012 to 2015, we used distance-sampling methods to conduct quarterly drive transects to estimate the abundance of 11 species of ungulates killed and consumed by lions (hereafter "lion prey"). We systematically drove the 12 transects in the control zone while two observers recorded species, group size, sighting distance, and bearing to all lion prey sighted along transects. To evaluate whether lions killed prey in proportion to their availability, we fit 5 female lions with Global Positioning System (GPS) collars (Vectronic Aerospace GmbH, Berlin, Germany) from March to May 2014. These lions belonged to 5 different prides, which collectively accounted for ca. 80% of lions on Ol Pejeta Ng'weno, 2017 in revision). Overlap between home ranges of prides was minimal (1-12%) throughout the course of our study (Ng'weno et al. in revision). All procedures were conducted with a veterinary team under the authority of the Kenya Wildlife Service.
From March 2014 to December 2015, we located lion kills from clustered locations of lions with GPS collars using an algorithm adapted from Knopff et al. (2009). Between August 2014 and December 2015, there were 246 instances in which prey carcasses were found to have been killed by lions Ng'weno, 2017 in revision). We then used Jacobs' index (Jacobs, 1974;Hayward and Kerley, 2005) to quantify seasonal selectivity for each of 11 species of lion prey, D = r−p 2r−2rp where r is the proportion of a given species among all kills, and p is its proportional abundance in the total prey population. Jacobs's index is bounded between −1 (highly avoided) and 1 (highly selected). Selectivity indices were calculated for each of the 11 species of lion prey, using carcasses from GPS clusters and estimates of prey abundance collected from August 2014 to December 2015.

Zebra Density and Risk of Predation to Hartebeest
To quantify the impact of primary prey (zebra) density on lion predation of hartebeest, we generated spatially-explicit density surfaces for zebra using resource selection functions (RSF) in a use-availability design (Manly et al., 2002). We constructed a minimum convex polygon (MCP) around all zebra herd locations from the quarterly surveys, and paired these with an equal number of random locations (n = 2,450) to achieve a 1:1 ratio of used to available locations. We then used a Rapideye satellite image (Digital Globe, Longmont, Co, USA) from May 2013 with 5-m spatial resolution to perform an unsupervised classification through isoclustering and maximum likelihood to group pixels with similar spectral reflectance into identify three habitat types: dense bushland (characterized by >50% cover of the tree Euclea divinorium), open bushland (characterized by 10-30% cover of the tree A. drepanolobium), and open grassland (characterized by <5% tree cover) (Birkett, 2002;Ng'weno et al., 2017). Following assignment of pixels, we ground-truthed our classification using 50-100 points in each habitat type. We performed all image processing using ERDAS Imagine, version 14 (Hexagonal Geospatial, Madison, Alabama) and ArcGIS version 10.3 (ESRI, Redlands, California).
Using the Euclidian distance and spatial join tools in ArcGIS 10.3, we extracted distances to the nearest water source and the nearest glade for each used and available location in each surveyspecific zebra MCP. We related locations of zebra herds to habitat type, distance to nearest water source, and distance to nearest glade in 30 × 30 m pixels. Collinearity between habitat variables (habitat type, distance to water, distance to glade) was minimal (r < 0.50; P > 0.20 for all possible pairwise combinations).
We used logistic regression to estimate RSF coefficients, with selection for or avoidance of a habitat variable indicated by coefficients >1.0 and < 1.0, respectively (Manly et al., 2002). We used the resultant coefficients to generate a zebra RSF, rescaled to create 16 continuous surfaces (one for each of the 16 quarterly surveys conducted during 2012-2015) with 30 × 30 m pixel values ranging between 0 (strongest avoidance) and 1 (strongest selection), divided into 5 bins of equal width following Morris et al. (2016). The first bin corresponded to the lowest probability of zebra selection (0.00-0.20) and the fifth bin corresponded to the highest probability of zebra selection (0.81-1.00). We combined survey-specific zebra densities with survey-specific RSFs to create spatially-explicit density estimates for zebra in each survey, which we validated independently using camera-trap grids (Appendix 1; Ng'weno et al.; Ng'weno, 2017 in revision). See Boyce et al. (2016) for a similar approach in using RSFs to estimate spatial variation in abundance.
Finally, we superimposed locations of hartebeest killed by lions (n = 27) obtained from GPS clusters to the spatially-explicit zebra density surface with which it was most closely associated in time. We used Chi-squared tests to assess differences in the proportion of hartebeest kills occurring within the highest (> 9.60 zebra/km 2 ) and lowest (<2.40 zebra/km 2 ) areas of estimated zebra density.

Glade Creation Experiment and Hartebeest Survival
For the period 2012 to 2015, we calculated survival rates of individual hartebeest, with biological years starting 1 October and ending 31 September of each subsequent year (corresponding to birth peaks). Survival of female adults, sub-adults, and calves collectively accounts for >70% of the variation in population growth of hartebeest at Ol Pejeta (Ng'weno et al., 2017). During our twice-monthly sight-resight drive surveys (see "Estimation of hartebeest abundance and population growth" above), we monitored survival of 179 adults, sub-adults, and calves within 11 groups in the control zone. Because they defend small (<5 km 2 ) territories, we were able to use unique marks (i.e., ear nicks, horn size, and shape, scars) to monitor survival through repeated surveys (Ng'weno et al., 2017). Because female hartebeest exhibit high site fidelity (Gosling, 1974), we equated disappearance with death. There were no instances in which an individual classified as "dead" during a particular survey was detected subsequently. Survival rates were calculated using the R package "survival" version 2.41-3 (Therneau, 2017).
Over the course of ∼3 years, abandoned boma sites transition into glades and become attractive to zebra because of high biomass of Cynodon and Pennisetum grasses (Veblen and Young, 2010;Porensky, 2011;Veblen, 2012). Glades vary in their sizes and shapes; we therefore selected and restricted our analyses to 37 glades derived from bomas that we established in 2009 and 2010 (Figure 4). Nineteen of these bomas were established along the edges of open plains as part of an experiment to understand how glades give rise to various edge effects (Porensky and Veblen, 2015;Porensky and Young, 2016). Subsequently, 18 additional bomas were established in the middle of open plains to further quantify resource selection of zebra and survival of hartebeest. All bomas were 17.2 ± 0.8 m SEM in diameter, and used by ∼200 cattle for 1 month.
To assess the influence of boma-derived glades on survival rates of hartebeest, we constructed 95% isopleths using fixed kernel density estimation to create utilization distributions (UD's), with least squares cross validation and a smoothing factor of 1,000 using R package "adehabitatHR" version 1.8.18 (Seaman and Powell, 1996;Powell, 2000;Calenge, 2006). We then overlaid glades onto hartebeest territories, and calculated the distance from the centroid of each sighting of each hartebeest group to the nearest glade.
To assess whether glades reduced hartebeest survival in accord with the hypothesis of apparent competition, we fit two Cox proportional hazards (PH) models (Cox, 1972) to a dataset on hartebeest mortality derived from two sources: the 27 kills that we detected from collared lion GPS clusters from May 2014 to December 2015, and 101 additional events in which individuals disappeared (and were equated with mortality events) from sight-resight surveys from January 2012 to May 2014. For every individual killed that was detected from GPS clusters, we noted absences through sight-resight surveys. No individuals disappeared in our sight-resight surveys between May 2014 and December 2015 that we did not independently discover through GPS kill-site clusters. In other words, 100% of the hartebeest in our study were killed by lions. Through the first PH model, we quantified how risk of mortality varied with the (categorical) presence of a glade within home ranges of hartebeest groups. Here, we assessed whether glade presence altered the hazard ratio (probability of death) of individual hartebeest over the course of 3 years. Because individuals belonged to one of 11 groups, and because groups were relatively sedentary over the course of our study, we incorporated frailties (random effects) associated with group identity (Fox et al., 2006; using R package "frailtypack" version 2.12-3 (Rondeau et al., 2012).
Finding random effects of group membership negligible (see "Results" below), we fit a second PH model in which distance to glade edge was used as a continuous covariate to explain hartebeest mortality (we did not attempt to combine frailties and continuous covariates within a single PH model, because both are challenging to integrate into PH models simultaneously; Lopez de Ullibarri et al., 2012). In calculating distances, we used the midpoint of the pair of sight-resight survey locations where an individual was last sighted with its group, and where it was first noted to have disappeared. All analyses were undertaken in R version 3.4.0 (R Development Core Team, 2017).

RESULTS
Exclusion of lions was sufficient to shift population growth of hartebeest from negative to positive, consistent with a lionmediated Allee effect (Figure 2). In the lion exclusion zone, population growth of hartebeest displayed a classic signal of negative density-dependence (Figure 2). In contrast, population growth of hartebeest increased with increasing population size in the control zone (Figure 2). Of the 11 species of ungulates killed by lions, zebra were consumed most frequently (40% of kills, n = 98) followed by warthog (15% of kills; n = 37), buffalo (14% of kills; n = 35), impala (14% of kills; n = 34), and hartebeest (13% of kills; n = 27). Lions consumed hartebeest and warthog disproportionately more than expected based on their abundance in both seasons, and zebra were consumed as frequently as expected based on their abundance in both seasons (Figure 3). The remaining 8 species of ungulates consumed by lions were avoided in at least one, and sometimes both, of the two seasons (Figure 3).
Zebra aggregated in and around glades (distance to glade: β = −2.30 ± 0.41 SEM, Z = −5.61, P < 0.0001), and the best supported RSF for zebra contained terms for distance to glade, distance to dense woodland, and distance to water (Supplemental Table 1). Hartebeest died more frequently than expected in areas of high zebra abundance (χ 2 4 = 42, P < 0.001; Figure 4). Consequently, and after nearly 3 years, survival among hartebeest lacking glades in their territories was over twice as high compared to those individuals whose territories contained glades (β = 2.35 ± 0.19 SEM, hazard ratio = 2.10, P < 0.001; Figure 5). Incorporating frailty terms to account for heterogeneity in survival among hartebeest groups did not significantly reduce PH model deviance (χ 2 10 = 11, P > 0.10). Using distance from the nearest glade's edge, we found evidence for spatial refugia around 600 m: with every 100 m closer to a glade, hartebeest incurred a roughly 15% greater chance of mortality (β = 0.14 ± 0.012 SEM, hazard ratio = 1.15, P < 0.001).

DISCUSSION
Multiple lines of evidence suggest that, over the course of 7 years, hartebeest were limited by apparent competition triggered by the restoration of lions on Ol Pejeta Conservancy. Lion exclusion erased Allee effects, and resulted in negative density dependence in hartebeest at populations below around 80 (Figure 2). Lions  consumed mostly zebra, and in proportion to their abundance, but selectively killed hartebeest (Figure 3). Zebra aggregated in and near glades, and mortality risk to hartebeest roughly doubled when glades occurred within hartebeest territories (Figure 5). Because ca. 80% of hartebeest mortality at Ol Pejeta is due to predation by lions (Ng'weno et al. in review), we attribute this change in patterns of mortality to a "gradient of consumption" (Orrock et al., 2008) that reduces survival of hartebeest in proximity to glades, and creates spatial refugia from lion predation beyond about 600 m from glades. Importantly, FIGURE 4 | Map with the locations of glades (≥3 years of age, n = 37) and hartebeest kill sites (n = 27) superimposed on a density surface for plains zebra. Estimates of zebra density (individuals/km 2 ) transitions from high (red) to low (green), and is driven largely by glade location.
FIGURE 5 | Survival rates and 90% confidence intervals from Cox proportional hazards models for hartebeest with territories containing vs. lacking glades. Survival for hartebeest with territories containing glades is roughly half that of those whose territories lack glades.
we did not study the mechanism underlying the Allee effect we observed, which reduces our ability to predict future hartebeest population trends at Ol Pejeta, or characterize the generality of our observations to other settings. Across vertebrates, Allee effects have been variously attributed to loss of genetic variability, reduced social facilitation, and difficulty finding mates. In our system, the effect appears linked to lion predation specifically: the higher preference of lions for hartebeest, or a reduction in shared vigilance in hartebeest herds. Lions may search for hartebeest groups, which tend to be in predictable places, and expend disproportionate effort to kill a hartebeest, once a group is detected.
Instances where apparent competition has been implicated in limiting ungulate populations are becoming numerous, including bighorn sheep (Ovis canadensis), mountain caribou (Rangifer tarandus), huemul (Hippocamelus bisulcus), and roan antelope (Hippotragus equinus; Harrington et al., 1999;DeCesare et al., 2010;McLellan et al., 2010;Wittmer et al., 2013;but see O'Brien et al., 2018 for an example contrary to the hypothesis of apparent competition involving Grevy's zebra [E. grevyi] in Laikipia). In our study system, spatial separation between zebra and hartebeest improved survival rates of hartebeest, probably by reducing encounters with lions hunting in areas with high zebra densities Ng'weno, 2017 in revision; see also Palmer et al., 2003;Forrester and Steele, 2004). Strategic placement of glades therefore offers a promising approach to creating refuges for hartebeest and perhaps other species of secondary prey. By siting bomas away from hartebeest, the establishment of glades was used to move zebras (and presumably lions) away from hartebeest territories.
The creation of spatial refugia and the reduction of primary prey alongside predators have previously been effective in releasing secondary prey from apparent competition in the past (Sinclair et al., 1998;Wittmer et al., 2013). Only the former is logistically feasible in our study system. Lions are Red-listed as Vulnerable by the IUCN and are at least as high a conservation priority as Jackson's hartebeest. At first glance, elimination of glades through reduction of cattle production might be considered another option for hartebeest conservation, but this also is impractical for at least two reasons. First, local ranchers are unlikely to reduce cattle numbers voluntarily. Second, reducing cattle would likely boost zebra numbers, and potentially lion numbers, increasing top-down control of hartebeest. This is because the diets of cattle and zebra overlap (Kartzinel et al., 2015) and both are likely limited by rainfall (as they are in the wider Laikipia region; Georgiadis et al., 2003Georgiadis et al., , 2007a and they probably compete in dry times and places (Odadi et al., 2011). Moreover, predation by lions on cattle is rare relative to predation on because zebra (Ogada et al., 2004;O'Brien et al., 2018). Alternative conservation interventions are required for the long-term persistence of lions and their prey not only on Ol Pejeta Conservancy, but more widely in Laikipia County and the whole of sub-Saharan Africa (Bauer et al., 2015).
The extent to which our results generalize and can be extended to other locales depends on several factors: whether the affinity of zebra for glades (a behavioral response) translates to a numerical response (population increase); the length of time over which glades are attractive to zebra, and how the creation and eventual disappearance of glades gives rise to a dynamic landscape of risk to hartebeest (Kohl et al., 2018). In turn, predation for hartebeest varies across the landscape in accord with the density of glades, and the degree to which the attraction of zebra to glades changes with glade density. If increasing density of glades increases zebra and thus lion abundance, predation pressure on hartebeest could also increase. Or, if glades only attract zebra but do not result in population increase, increasing glade density could further disperse zebras across landscapes, negating the ability of glades to aggregate zebra.
We showed that spatial and temporal dynamics of predators and prey on Ol Pejeta Conservancy are consistent with apparent competition theory, and suggested a practical application for how lions and their secondary prey can coexist. Using strategic livestock management to manipulate the spatial distribution of predation pressure, we have provided a possible way to alleviate top-down control on an ungulate of conservation concern (Jackson's hartebeest) without resorting to lethal control of their most important predator (lions). Such solutions for conservation of lions in non-protected areas are needed to complement traditional funding efforts for formally-protected areas (Lindsey et al., 2018).

AUTHOR CONTRIBUTIONS
CN, SB, NG, and JG conceptualized the research problem. CN and JG designed the study and performed the analyses. BG and AK implemented lion capture and assisted with analyses. LP implemented the boma experiment. DR collected and analyzed data on abundance and resource selection of plains zebra. CN, SB, and JG wrote the manuscript, with edits and input from all authors.  Fox. Cons. Biol. 21, 1082-1091. doi: 10.1111/j.1523-1739. 00721.x Augustine, D., and McNaughton, S. (2006. Interactive effects of ungulate herbivores, soil fertility, and variable rainfall on ecosystem processes in a semi-arid savanna. Ecosystems 9, 1242-1256. doi: 10.1007/s10021-005-0020-y