Together, our results from both experimental and observational systems showed that parasite density was greatly elevated at water sources, but that this varied by herbivore species. Notably, cattle and elephants drove variation in parasite eggs from dung, due to their high biomass, high infection intensity, and high degree of aggregation at water (Figure 2: A, B, C, SI Appendix Fig. S6). Importantly, drier contexts drove increased aggregations and parasite density near water.
Q1: To what extent do water sources concentrate hosts, their feces, and fecal-oral parasites across herbivore species?
Water removal from five ~25,000 liter pans resulted in significantly reduced total, grazing, and drinking activity (measured via camera traps) at experimental water sources relative to filled water sources for elephants and cattle. While drinking activity was significantly reduced for all animals together, the 46% drop (i.e. a two-fold reduction) in total activity at experimental pans was only marginally significant (t = -1.62, p = 0.1 for the drained x during interaction, SI Appendix, Tables S7, S8, Fig. S4). Interestingly, total, grazing, and drinking activity increased when water was replenished such that it even trended higher compared to pre-experimental levels (t = 2.27, p = 0.02 for the drained x post interaction, SI Appendix, Tables S7, S8).The interaction between experimental status and treatment was not an important parameter for models of buffalo, zebra, giraffe, or impala activity. However, after pans were refilled, zebra and buffalo activity was significantly higher at experimental pans relative to other phases of the experiment (p = 0.03, p = 0.01 for total activity, p = 0.09 and p < 0.001 for grazing, p = 0.01 and p < 0.001 for zebra and buffalo respectively; SI Appendix, Tables S7, S8).
Dung density at filled water sources relative to experimental water sources increased when water was drained for all animals together and herded cattle and elephants separately (Table 2). This effect was largest for elephants: when experimental pans were drained, dung density was over six times higher at filled pans (in the area closest to water), while it was no different pre-draining or post-refilling (t = 2.88, p = 0.004 for increased probability of zeros for the drained x during interaction, Table 2). We found a similar pattern for cattle, as dung aggregation at filled pans was over three times higher during the drained period (at the 0m mark), but not during the ‘pre’ phase (t = 3.09, p = 0.002 for the increased probability of zeros for the drained x during interaction, Table 2). While cattle dung density was significantly lower at drained pans relative to filled pans after refilling, this effect was substantially reduced. Since cattle and elephants accounted for the largest proportion of dung density (> 50%, Figure 1B), this drove a similar pattern for total dung. However, we detected no effect of water draining for zebra, impala, or giraffe considered separately, and buffalo dung density was slightly higher at experimental sites after refilling.
Total parasite density (order Strongylida) – estimated as the product of dung volume in the environment, dung physical density, and median fecal egg count for each species – was three times higher at filled pans compared to experimental pans during the experiment but was no different before or after (t = -2.93, p = 0.003 for the drained x during interaction, Table 2, Figure 3, see SI Appendix Fig. S7 for species-specific responses).
Finally, parasite egg counts in dry soil near water pans were relatively low but consistent across treatments throughout the experiment (X22 for status x treatment = 0.14, p = 0.93), suggesting that water removal did not substantially affect dry soil parasite density. However, after including eggs found in wet soils, density was 16 times higher at filled pans during the experiment, but not significantly different before or after (based on dry weight, X22 for status x treatment = 9.30, p = 0.01), since wet soils were not present at drained pans. Additionally, soil egg density in dry soils around both pan treatments averaged 4.5 times that at non-water sites (X2 for treatment = 5.58, p = 0.02).
Q2: To what extent are water sources parasite exposure hotspots across herbivores?
Combining dung parasite density results (Figure 2: A, B, C) with herbivore grazing behaviors showed that total parasite exposures for cattle and elephants were more than an order of magnitude higher near water if parasite mortality was equivalent at water and matrix sites (Figure 2: D, E, F). For the observational system, parasites of buffalos and elephants showed the strongest increase in total estimated exposures near water (85 and 39 times respectively, p = 0.09 and 0.002), with parasites of all other species except impala also trending higher. The strong results for buffalo were driven by the very low levels of buffalo observed grazing at matrix sites. With increasingly more conservative assumptions about parasite mortality differences across site types (assuming higher mortality of parasites near water), this effect diminished, such that at our most conservative assumptions (10x higher parasite mortality near water) only buffalo, elephants, and cattle had elevated estimated parasite exposures at water, although this increase was not statistically significant (Figure 2F).
For the experimental system comparisons between permanently filled water sources and matrix sites, parasite exposure was elevated near water for cattle, elephants, zebra, and impala (143, 67, 20, and 8 times higher respectively, p < 0.001 for all but impala for which p = 0.005). This remained significantly elevated (14 times and 7 times higher) for cattle and elephants, even in the scenario when parasite mortality between egg and infection was 10 times greater near water (p < 0.001, p = 0.002 respectively; Figure 2F).
Q3: How does herbivore activity, herbivore dung, and parasite density at watering holes vary across rainfall contexts?
Total herbivore grazing activity measured from camera traps was twice as high at water sources compared to matrix sites (t = 2.51, p = 0.01 for all species together), and it was significantly elevated for elephants, and marginally elevated for buffalo and zebra (Figure 2E). We found no significant interaction between mean annual precipitation (MAP) and animal grazing at water although both predictors were important separately (SI Appendix Table S9). This was likely due to the short deployment duration and low statistical power, as activity declined with increasing annual precipitation at both matrix sites and water sources, and it tended to be further elevated near water in drier areas (SI Appendix Table S9, Fig. S5). For all species except impala, grazing activity trended higher at water compared to matrix sites.
MAP and prior rainfall were important parameters in cattle, elephant, zebra, and total dung density models (Table 2, SI Appendix Fig. S8). In dry locations (~ 460 mm/year) close to water and following no rainfall, cattle dung density was three orders of magnitude higher at water relative to matrix sites, but this elevated density decreased as MAP and outward distance increased. This pattern was also strong for elephants: in dry areas following periods of no rainfall, elephant dung was approximately ten times higher close to water, but this effect weakened as MAP and outward distance increased (Table 2, SI Appendix Fig. S8). Zebra dung density was no different between water and matrix sites when there was little prior rainfall or low MAP, and we even observed potential aversion to water during the wettest periods in high MAP areas. Impala dung density was also slightly elevated near water in low-rainfall locations but depressed near water in wet conditions. We observed slightly higher dung density levels at watering holes relative to matrix sites for buffalo and giraffe in low-rainfall conditions, and there was a significant interaction between MAP and prior rainfall for giraffe (Table 2, SI Appendix Fig. S8).
Critically, outward distance from water, MAP, and prior rainfall all modulated parasite density at water sources compared to matrix sites. In areas were MAP was lowest (450 mm/year) and prior 30-day rainfall was 0mm, parasite egg density from dung was estimated to be more than 150 times higher than matrix sites in the closest area to water. This effect decreased sharply as MAP, prior rainfall and outward distance increased (MAP: t = 3.40, p = 0.001; prior rainfall: t = 5.06, p < 0.001; distance: t = 3.01, p = 0.003, Table 2; Figure 4).
In our negative binomial model of eggs found in soil, parasite densities differed significantly based on water proximity (X23 = 135.7, p < 0.001). Parasites (eggs per 20 grams dry soil) were more than two orders of magnitude higher in damp soil (based on dry weight, mean ± SE = 31.1 ± 8.0), and four times higher in dry soil (1.02 ± 0.35), near the water’s edge compared to locations 1 km from water (0.26 ± 0.13) (SI Appendix Fig. S9).