Elephant population growth in Kruger National Park , South Africa , under a landscape management approach

African savanna elephants (Loxodonta africana) significantly influence several conservation management objectives that authorities seek to achieve in protected areas. Managers of protected areas typically seek to conserve biodiversity – elephants act as ecological engineers that may alter habitats and various aspects of biological diversity (Valeix et al. 2011). Protected areas also serve as focal places where elephant presence forms a key component of tourist attractions and experiences (Kerley, Geach & Vial 2003). For people living next to parks, elephants leaving protected areas threaten livelihoods and lives (Osborn & Parker 2003). South African National Parks (SANParks) seek to manage these effects that elephants have on the ability of managers to achieve conservation management objectives.


Introduction
African savanna elephants (Loxodonta africana) significantly influence several conservation management objectives that authorities seek to achieve in protected areas. Managers of protected areas typically seek to conserve biodiversity -elephants act as ecological engineers that may alter habitats and various aspects of biological diversity (Valeix et al. 2011). Protected areas also serve as focal places where elephant presence forms a key component of tourist attractions and experiences (Kerley, Geach & Vial 2003). For people living next to parks, elephants leaving protected areas threaten livelihoods and lives (Osborn & Parker 2003). South African National Parks (SANParks) seek to manage these effects that elephants have on the ability of managers to achieve conservation management objectives.
The effects that elephants have on biodiversity, tourist experiences and rural people's livelihoods may associate primarily with their location rather than abundance or density (Van Aarde, Jackson & Ferreira 2006). In small parks, however, density of elephants may aggravate the effects they have on biodiversity. The way elephants use landscapes depends on the distribution of water (Chamaillé-Jammes, Valeix & Fritz 2007;Loarie, Van Aarde & Pimm 2009), food (Codron et al. 2006) and shade (Kinahan, Pimm & Van Aarde 2007) -resources that elephants select at different intensities during different seasons (Young, Ferreira & Van Aarde 2009a). In addition, elephants avoid people by choice (Douglas-Hamilton, Krink & Vollrath 2005;Graham et al. 2009). Within this context SANParks embraced a landscape approach to elephant management and seek to restore the landscape variability and limitations on elephant spatial and temporal dynamics (Harris et al. 2008). Within Kruger National Park (KNP), this translated to stopping the control of elephant numbers through culling (Van Aarde, Whyte & Pimm 1999), the closure of boreholes and destruction of several earthen dams (Gaylard, Owen-Smith & Redfern 2003), the removal of fences between Kruger National Park and adjacent conservation-friendly land-uses (Venter, Naiman, Biggs & Pienaar 2008), and the use of a fire management policy focusing on patch burns (Van Wilgen et al. 2004).
The creation of spatial and temporal heterogeneity in elephant resources predicts responses in how elephants use landscapes. Elephants use landscapes more variably (Smit & Ferreira 2010) South African National Parks (SANParks) manage landscapes rather than numbers of elephants (Loxodonta africana) to mitigate the effects that elephants may have on biodiversity, tourism and stakeholder conservation values associated with protected areas. This management philosophy imposes spatial variability of critical resources on elephants. Restoration of such ecological processes through less intensive management predicts a reduction in population growth rates from the eras of intensive management. We collated aerial survey data since 1995 and conducted an aerial total count using a helicopter observation platform during 2015. A minimum of 17 086 elephants were resident in the Kruger National Park (KNP) in 2015, growing at 4.2% per annum over the last generation of elephants (i.e. 12 years), compared to 6.5% annual population growth noted during the intensive management era ending in 1994. This may come from responses of elephants to density and environmental factors manifested through reduced birth rates and increased mortality rates. Authorities should continue to evaluate the demographic responses of elephants to landscape scale interventions directed at restoring the limitation of spatial variance in resource distribution on elephant spatiotemporal dynamics and the consequences that may have for other conservation values.
Conservation implications: Conservation managers should continue with surveying elephants in a way that allows the extraction of key variables. Such variables should focus on measures that reflect on how theory predicts elephants should respond to management interventions.

Read online:
Scan this QR code with your smart phone or mobile device to read online. http://www.koedoe.co.za Open Access following the restoration of various ecological processes, where previous management philosophies and approaches fostered homogeneity (Venter et al. 2008). Variability in landscape use also predicts population numerical responses -when daily distances travelled by cows increase, weaned calves have higher mortality rates (Young & Van Aarde 2010). The closure of boreholes increased distances and configuration (Gaylard 2015) required for travel between water and good quality food. In addition, when SANParks stopped culling, elephant densities increased (Van Aarde, Whyte & Pimm 1999).
Higher densities associated with a less clumped distribution of elephants across the Kruger landscape (Young, Ferreira & Van Aarde 2009b). These mechanisms predict decreased population growth rates with increasing density and decreasing water provisioning.
SANParks seek to evaluate the consequences of reduced management intensity directed at mitigating the effects of elephants on various conservation management objectives. Monitoring focuses on spatial and temporal responses of elephants. As part of this monitoring, SANParks survey elephants to derive estimates of population variables. Here we report on the KNP elephant survey of 2015, trends extracted from recent surveys (SANParks, Unpublished data 1 ) and check population growth rates post 1994 against those at the end of the intensive elephant management era in the mid-1990s (Whyte 2001).

Methods
The KNP, 19 485 km² in size, spans the low-lying savannas of the eastern parts of the Limpopo and Mpumalanga Provinces, South Africa. Mozambique abuts the Park in the east, while Zimbabwe forms the boundary to the north. The Park has a mean annual rainfall that varies from 750 mm in the south to 450 mm in the north, with ± 80% of rain falling during summer (Gertenbach 1980). Granite and gneiss soils dominate KNP's western half and nutrient-rich basalt soils dominate the eastern half, with a narrow band of Karoo sediment occurring at the junction of the granite and basalt (Schutte 1986). The vegetation on the southern basalts is largely wooded savanna, with Sclerocarya caffra and Vachellia nigrescens dominating the tree canopy. Mixed Combretum spp. and Vachellia spp. dominate the southern granites. In the north Colophospermum mopane dominates both the basalt and granite substrates (Gertenbach 1983). This produces a diverse range of landscapes resulting in spatial and temporal heterogeneity of resources that elephants use.
Water is an essential resource provided through several perennial (i.e. Crocodile, Sabie, Olifants and Letaba) and large seasonal rivers (i.e. Biyamiti, Nwaswitshaka, Nwatswitsontso, Timbavati, Tsendze, Shingwedzi, Levhuvhu and Limpopo). The provision of additional water, through drilling of boreholes and construction of several dams that started in 1933, disrupted spatial and temporal heterogeneity From these surveys, we summed all observations to find the total number of elephants observed. We collated the total number of elephants observed in surveys since 1995. These values are estimates of minimum numbers of elephants known to be alive and do not account for availability, detection probability, observer bias or statistical sampling error (Caughley 1974).
We fitted population models that included exponential: and heterogeneity: models where N t is minimum number of elephants alive at time t, r is exponential growth rate, r m is maximum growth rate, K is an estimate of equilibrium population size, R is mean annual rainfall and α is a constant. We estimated r, r m , K and α using maximum likelihood approaches (Edwards 1972;Microsoft Excel macro provided by Hood 2005). We used Akaike's Information Criteria (AIC c ) to select the most appropriate model (Johnson & Omland 2004). When differences between AIC c values were less than 2, we used model averaging (Claeskens & Hjort 2008) http://www.koedoe.co.za Open Access to obtain a population growth estimate and associated confidence intervals.
The recording of calves allowed us to derive a recruitment rate as the fraction of elephants less than 1 year old. Such a recruitment rate is the consequence of calves born and their subsequent survival during the first year. To estimate overall elephant survival (s) we calculated an age distribution parameter a, the rate at which age-specific frequencies decay with age in stable age distributions from: where n x is the frequency of individual x years of age and w is elephant longevity (Ferreira & Van Aarde 2008). We set w at 55 years (Lee et al. 2016). We estimated a using maximum likelihood approaches (Edwards 1972; Microsoft Excel macro provided by Hood 2005). This allowed us to calculate survival from a = λs (Eberhardt 1988) where λ is finite population growth (λ = 1+r). To minimise the consequences of variation in age distributions, we recalculated population growth using only estimates spanning an elephant generation defined as the age at which member of a given cohort is expected to reproduce (Charlesworth 1994), that is, 12 years (Robson 2015) from 2004 until 2015 and used this as input to calculate survival rate. We could then estimate birth rate by dividing recruitment rate with the average estimated survival rate across all age classes.
SANParks' intensive management philosophy ceased in 1994 that is, culling stopped, water provisioning decreased, some fences were removed and patch burns were implemented since then. We thus extracted population growth rates as well as estimates of survival and recruitment rates prior to 1994 (Whyte 2001) to check against recent estimates of these variables. Note that mean annual rainfall prior to 1994 was 510.1 mm (SE = 38.2, n = 18) similar to that after culling stopped (565.1, SE = 44.5, n = 20) and from 2004 to 2015 (529.4, SE = 34.0, n = 11). We thus did not consider rainfall driven influences in our comparisons.

Results and discussion
We encountered elephants throughout KNP (Figure 1). Two of the large perennial rivers dissecting the Park -the Olifants and Sabie Rivers served as boundaries for regional estimates. We provide estimates for areas north of the Olifants River, between the Olifants and Sabie Rivers and south of the Sabie River separately. Observers encountered 7107 and 1163 elephants in mixed (n = 527) and bull herds (n = 546), respectively, north of the Olifants River at a density of 0.84 n.km -2 . Between the Sabie and Olifants Rivers observers noted 3526 and 462 elephants living in mixed (n = 289) and bull herds (n = 229), respectively, at a density of 0.81 n.km -2 . Highest density of 1.15 n.km -2 was south of the Sabie River and came from observers recording 4232 and 596 elephants in mixed (n = 306) and bull herds (n = 302), respectively. These regional counts resulted in a minimum of 17 086 elephants living in KNP during 2015.
Models describing the trends in minimum elephants known to be alive after culling stopped were equivalent (r 2 values ranged from 0.96 to 0.98) ( Table 1). Model averaging (Figure 2) thus recorded annual exponential growth rate since 1995 at 4.9% (2.1% -7.6%) per annum, similar to that reported previously (Young et al. 2009b). The point estimates of 6.5% reported prior to 1995 (Whyte 2001), the period of intensive management, is located at the 86.7th percentile of growth since 1995. Population growth since 2003, an elephant generation, was 4.2% (1.1% -7.3%) with 6.5% located at the 92.7th percentile of growth since 2004. These results may reflect a likely reduction in elephant population growth rate over time and as densities increased supporting the prediction that population growth rates should be declining. Note that poaching had little effect on the elephants living in KNP with between 0 and 22 elephants poached during any year (SANParks, unpublished data 3 ). It is uncertain, however, how poaching of elephants could change in future and what the impact would be on the population in KNP.
We acknowledge that movements of elephants to areas abutting KNP can also account for reduction in population growth rates that we noted. For this reason, we also estimated vital rates given that previous studies implied variance in birth rates (Robson 2015) and declines in survival rates as mechanisms associated with increase landscape limitations on elephant resources (Young et al. 2009b;Young & Van Aarde 2010). The estimated population growth since 2003 and an estimated age distribution parameter over that same period (a = 0.94 ± 0.02, n = 10) allowed us to estimate survival during the first year at 0.982 (0.941-0.999), while other age classes had a survival rate of 0.983 (0.954-0.999). This estimate ignores agespecific survival rates for elephants older than 1 year. The overall recruitment rate (5.67%, 5.23% -6.02%) predicted an annual birth rate of 5.78% (5.33% -6.41%). Note that during the period of intensive management estimates of annual birth rates (expressed as potential rates of increasei.e. rate of increase in the absence of mortality) ranged from 7.8% to 16.8% (Whyte 2001). Demographic responses of elephants to restoration of ecological processes may thus provide some potential mechanisms explaining how elephant growth rates declined since 1995. Recent analyses, however, highlight that variability reflecting on ecological productivity during the conception and gestation phases of elephant reproduction drives variability in elephant birth rates (Robson 2015).
The estimated survival rate translates to an overall mortality rate of 1.8% (1.72% -1.88%) higher than the point estimated of 1.5% reported for the period of intensive  management (Whyte 2001). We are mindful of confounding influences of different methods, but our results could support the predictions derived from mechanism of how changes in resource distribution and density influence survival schedules of elephants (Young & Van Aarde 2010

Conclusion and recommendations
The temporal responses of elephants following the implementation of a philosophy of managing landscapes rather than managing elephants directly (   Minimum number of elephants FIGURE 2: Minimum number of elephants living in the Kruger National Park since 1995 after intensive management of elephants stopped in 1994. We also illustrate the average model (line) and its 95% confidence intervals (thin lines) since then.