Abstract Introduction Study area Methods Results Discussion Conclusions and management implications Acknowledgements References Appendix 1 About the Author(s) Helga van der Merwe South African Environmental Observation Network, Kimberley, South Africa Plant Conservation Unit, University of Cape Town, Cape Town, South Africa Noel van Rooyen South African Environmental Observation Network, Kimberley, South Africa Hugo Bezuidenhout South African National Parks, Kimberley, South Africa Applied Behavioural Ecology and Ecosystem Research Unit, University of South Africa, Florida, South Africa Jacobus du P. Bothma Centre for Wildlife Management, University of Pretoria, Pretoria, South Africa Margaretha W. van Rooyen South African Environmental Observation Network, Kimberley, South Africa Department of Plant and Soil Science, University of Pretoria, Pretoria, South Africa Citation Van der Merwe, H., Van Rooyen, N., Bezuidenhout, H., Bothma, J. du P. & Van Rooyen, M.W., 2019, ‘Vachellia erioloba dynamics over 38 years in the Kalahari Gemsbok National Park, South Africa’, Koedoe 61(1), a1534. https://doi.org/10.4102/koedoe.v61i1.1534 Original Research Vachellia erioloba dynamics over 38 years in the Kalahari Gemsbok National Park, South Africa Helga van der Merwe, Noel van Rooyen, Hugo Bezuidenhout, Jacobus du P. Bothma, Margaretha W. van Rooyen Received: 16 Mar. 2018; Accepted: 21 Nov. 2018; Published: 29 Apr. 2019 Copyright: © 2019. The Author(s). Licensee: AOSIS. This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Vachellia erioloba is a keystone tree species in the southern Kalahari. This long-term study over nearly four decades tracks two populations in different landscapes (the interior sandy duneveld versus the clayey Nossob riverbed) of a large conservation area and offers valuable data on this species under natural soil moisture conditions and with limited anthropogenic influences. In 1978, 18 trees were permanently marked in a 1 ha plot in the interior duneveld of the Kalahari Gemsbok National Park (Dankbaar site). In the Nossob riverbed all trees in a 1 ha plot were surveyed in 1979 (Grootkolk site). At both sites, tree height and stem circumference were subsequently measured at irregular intervals until 2016 in order to investigate growth rates and population structure. Of the 18 marked trees at Dankbaar, six died and three showed coppice regrowth following substantial dieback after a fire. A mean height increase of 60 mm/year was recorded and the mean height of the remaining uncoppiced trees was 6.8 m in 2016. Stem diameter growth rate per year varied widely between trees and between years with a mean value of 2.5 mm/year over the 38-year period. Growth rate calculated for three 10-year intervals varied. Using the mean growth rate derived in the current study and stem size of the dead trees, the mean age of the trees when they died was estimated. At the Grootkolk site, the position of the centroid in relation to the midpoint of the diameter class range suggests that this population is gradually becoming a mature to old population with limited recruitment. This was supported by the size class distribution curves. However, no differences between slopes or intercepts of the stem diameter size class distributions were found. Conservation implications: This study was conducted in a large conservation area, that is, a natural ecosystem excluding most of the anthropogenic threats that are present outside of the park. The study illustrated that in the duneveld the population studied was self-sustaining, with recruitment occurring and large individuals presumably dying of old age. Although fire caused a few individuals to coppice, no fire-related deaths were reported. In the Nossob riverbed, surveys started in a stand of predominantly young trees and the size class distribution at that stage already showed a lack of recruitment. This stand is ageing and will likely disappear at this site; however, new young stands are appearing at other sites in the Nossob riverbed. Under the current conditions with negligible anthropogenic influences, it therefore appears that some V. erioloba populations in the park are increasing in size while others are decreasing, but that overall the species will persist. The impact of global climate change on this species is, however, unknown. Keywords: conservation area; growth rate; Kgalagadi Transfrontier Park; long-term study; size class distribution; tree age. Introduction Ecological patterns generated from a single snapshot in time should be used carefully when trying to identify the underlying processes driving the system (Jeltsch, Moloney & Milton 1999), and long-term data can be very useful to assist with the interpretation of ecosystem patterns and processes. However, appropriate long-term data are difficult to find in the case of long-lived organisms such as trees, whose lifetimes are usually considerably longer than those of researchers (Menges 2000). In arid systems, this is further exacerbated by low rates of biotic change (Cody 2000; Lawley et al. 2013; Van Rooyen et al. 2015). In such systems, to accommodate climatic fluctuations, the monitoring of woody species should take place over a timespan in excess of 20 years (O’Connor 1985). In arid savannas, the tree component may seem to be constituted by even-aged cohorts (Midgley & Bond 2001; Wiegand, Ward & Saltz 2005); however, what appear to be cohorts among adult trees may not be even-aged individuals (Seymour 2008). Even-aged stands are interpreted as evidence of the importance of abiotic factors, particularly rainfall, to produce pulsed recruitment. Apparent even-aged stands signify the release from fire or herbivory whereby different-aged saplings are allowed to recruit simultaneously into larger height classes, thus appearing of equal age (Midgley & Bond 2001). In the southern Kalahari, Vachellia erioloba is regarded as a keystone species because it provides food, shelter, perches, nests and roost sites to many types of animals and harbours a distinct assemblage of plant species below its canopy (Barnes 2001; Leistner 1967; Milton & Dean 1995; Seymour 2008; Steenkamp et al. 2008). It is listed as a protected tree species in South Africa (Government Notice no. 1602, 2016). There are, however, growing concerns over increased mortality in mature V. erioloba in the southwestern Kalahari (Schachtschneider & February 2013). The reasons proposed for the increase in mortality include an increase in fire frequency (Seymour & Huyser 2008), increasing groundwater abstraction (Powell 2005; Shadwell & February 2017), harvesting (Seymour 2008), diseases (Slippers et al. 2014), arboricides (Powell 2005) and competition for resources from alien invasive species such as Prosopis (Robertson & Woodborne 2002; Schachtschneider & February 2013; Zachariades, Hoffman & Roberts 2011). These threats are especially relevant outside of conservation areas, emphasising the importance of protected environments in which V. erioloba individuals, exposed only to natural threats such as fire, can be studied. Vachellia erioloba is a long-lived tree with a lifespan in excess of 240 years (Steenkamp et al. 2008). Mature individuals of this species are reportedly not easily affected by short-term climatic variations (Barnes 2001), possibly because individuals have some of the deepest roots of any known species (the maximum recorded is 68 m), allowing them access to deep groundwater sources (Jennings 1974 in Moustakas et al. 2006). Quantitative information on growth rates and lifespan of the species are, however, still inadequate (Steenkamp et al. 2008). From 1978, surveys were initiated in the Kalahari Gemsbok National Park (KGNP) to gather information on the growth rate, lifespan and population dynamics of V. erioloba in the interior duneveld (Dankbaar) and in the Nossob riverbed habitats (Grootkolk). The objectives of this paper were (1) to report on the mean annual increase in stem diameter and tree height of the 18 marked V. erioloba trees at the Dankbaar site; and (2) to document the mean age of those trees that died in the interior duneveld of the park; and (3) at the Grootkolk site to describe changes in the size class distribution of a V. erioloba population in the Nossob riverbed. The publication of these findings contributes to our understanding of V. erioloba population dynamics within large conservation areas without many of the anthropogenic factors proposed to increase mortality in this species. Study area The KGNP is situated in the southern Kalahari (24°15’S–26°30’S and 20°00’E–20°45’E) and covers approximately 9600 km² (Figure 1). It forms part of the Kgalagadi Transfrontier Park, together with the Gemsbok National Park in Botswana. FIGURE 1: Location of the two Vachellia erioloba study sites, Dankbaar situated in the interior duneveld and Grootkolk located in the Nossob riverbed, within the Kalahari Gemsbok National Park. The mean annual rainfall recorded from 1976 to 2015 at Nossob Rest Camp in the park is 194 mm (Figure 2), with most rain falling between January and April. The annual rainfall has a high coefficient of variation and the rain often falls as short-duration, high-intensity thunderstorms. Temperatures show a large range, with winter lows reaching -10.3 °C and summer highs reaching 45.4 °C (Van Rooyen & Van Rooyen 1998). FIGURE 2: Total annual rainfall (January to December) at Nossob Rest Camp from 1976 until 2016. The orange line represents the mean annual rainfall of 194 mm. Sampling was conducted at two sites (Figure 1). The first site is situated 2.7 km east of Dankbaar (1040 m above sea level) in the interior duneveld in the Acacia (Vachellia) erioloba – Schmidtia kalahariensis low duneveld landscape (Van Rooyen et al. 2008). This open tree savanna, dominated by V. erioloba (tree) and S. kalahariensis (grass), occurs in the interior duneveld in the northern part of the park. The aeolian sands of the Low Duneveld landscape are piled into a gently undulating landscape. Beneath the sand lies a vast sheet of calcareous or silicified sand or sandstone that contains grits and minor conglomerates (Malherbe 1984). Van Rooyen (1984) describes the duneveld soils as deep, red, eutrophic, fine and medium sandy soils of the Hutton soil form, or yellow-brown, eutrophic and calcareous, fine to medium sandy soils of the Clovelly soil form (Soil Classification Working Group 1991). The second site, at Grootkolk (1041 m above sea level), is situated in the Nossob riverbed in the Acacia (Vachellia) erioloba – Acacia (Vachellia) haematoxylon riverbeds landscape (Van Rooyen et al. 2008). This landscape is characterised by large scattered V. erioloba trees with a dense grass stratum dominated by Panicum coloratum. Van Rooyen (1984) described the Nossob riverbed’s soil as deep, calcareous, reddish-brown and grey-brown loam with clay soils of the Oakleaf, Dundee and Valsrivier soil forms. Upstream from the Grootkolk site, shallow, brown and yellow-brown calcareous sands and loams of the Mispah soil form occur (Soil Classification Working Group 1991). Fine-grained soils in arid areas show the so-called inverse texture effect (Noy-Meir 1973), with the result that there is less available soil water and consequently lower production on fine-textured soils than on coarse-textured soils (Van Rooyen et al. 1984). Methods Field surveys Two study sites (Dankbaar site and Grootkolk site) were surveyed on numerous occasions; however, survey methods used at these two sites differed. In 1978, all V. erioloba trees were marked in a 1 ha plot at the Dankbaar site. Tree height and stem circumference were recorded for each of the 18 permanently marked trees at irregular intervals from 1978 to 2016 (survey years 1978, 1980, 1983, 1985, 1988, 1989, 1991, 1994, 1998, 1999, 2008 and 2016). In 2016, for the first time since the original survey in 1978, all additional trees within the plot were measured and permanently marked. At the 1 hectare (ha) Grootkolk site, height and stem circumference of all V. erioloba individuals were measured at irregular intervals from 1979 to 2016 (survey years 1979, 1984, 1986, 1989, 1994, 1996, 1999, 2008 and 2016). The stem circumferences of all dead trees were also recorded at Grootkolk in the year of the survey when the tree was found to have died. Analysis Because of the different survey methods used at the two sites; the data collected could not be analysed in the same manner for both the sites. For example, only the 18 marked trees at Dankbaar were surveyed throughout the 38-year period; thus an analysis of the population structure is inappropriate as the entire population in the 1 ha plot was not surveyed. An estimate of the growth rate per year over the entire study period (38 years) was calculated for each marked tree at the Dankbaar site by subtracting the first measured stem diameter (1978 survey) from the last live measured stem diameter and dividing by the number of years between these two measurements. Similarly, this method was used to calculate changes in height of a tree. A linear regression using GraphPad 7.02 (GraphPad Software, www.graphpad.com) of stem diameter versus annual growth rate was used to investigate the relationship between these two variables. Additionally, stem diameter growth rates over the first 10-year period (1978–1988), second 10-year period (1988–1998) and third 10-year period (1988–2008) were calculated for each marked individual. This was calculated in the same way as for the entire period, but using the first and last measured diameters for the appropriate periods and calculating a mean across all the individuals. Dead and coppiced individuals were excluded from the calculations from the survey year in which the death or coppicing was recorded for the first time. The mean stem diameter in 1978 and 2016 of all marked trees was used to estimate the average age using the growth rate determined in the current study. The average age of all the marked trees in 1978 and in 2016 were calculated by dividing the mean stem diameter for the particular year by the growth rate determined in the current study. For comparative purposes the mean age of the trees was also calculated by the same method using the mean diameter growth rate for only the interior duneveld (3.2 mm/year) and for the interior duneveld and riverbed combined (3.6 mm/year) reported by Steenkamp et al. (2008). Age was also determined using the growth rate for the interior duneveld and riverbed combined (5.1 mm/year) based on the equation derived for the relationship between radiocarbon age and stem circumference (Steenkamp et al. 2008). These calculations were made using the mean stem diameter of the live trees in 2016 as well as the last ‘alive’ measured stem diameter for the dead trees. These calculations were also performed for the Grootkolk population using the mean stem diameter of the stem circumference classes. The stem circumferences for trees at the Grootkolk 1 ha plot were classified into 12 equal stem circumference classes (0–20, > 20–40, > 40–60, > 60–80, > 80–100, > 100–120, > 120–140, > 140–160, > 160–180, > 180–200, > 200–220, > 220–240 starting with class 1 and ending with class 12) for each of the survey years. The size class distribution was examined visually to classify the population into one of four population structure types (Gaugris & Van Rooyen 2007; Figure 3). In this model, Type I size class distribution represents a growing population and follows an inverse J-shaped curve, which represents the ideal curve shape for a growing population in a natural environment. There is a constant decrease in the number of individuals as the size classes increase, and the curve spans the whole range of size classes described for the species. Populations with a Type II curve have a near ideal population structure that differs from a Type I only by having a lower abundance of smaller size classes. Populations with Type III size class distributions are typically associated with disturbance. The Type IIIa curve broadly resembles a Type I size class distribution but applies to populations where a section of the curve spanning one or several size classes is missing or where abundance is lower than expected. The typical large pool of small individuals is generally missing in Type IIIb and the size class distribution shows a succession of minor peaks, even with some gaps where one or several size classes may be missing altogether. In the Type IV curve, one or two peaks are observed over a restricted range of size classes, and this curve thus represents an abnormal size class distribution. The curve is either bell-shaped with a peak more or less in the middle of the size class range, or it shows a high frequency of large size classes, with little or no representation of smaller classes. FIGURE 3: Idealised shapes of the four different population curve types. For the 1 ha Grootkolk site the number of individuals recorded per stem circumference size class was calculated in order to obtain the density of trees per size class. Regressions were performed on the density (Di) of individuals against size class (Mi) (Gaugris & Van Rooyen 2007). The diameter class midpoint (Mi) was calculated as the mean of the upper and lower limit of each size class (Condit et al. 1998). Logarithmic transformations, ln (Di + 1) and ln (Mi), were used to standardise the data (Lykke 1998; Niklas, Midgley & Rand 2003) before performing least square linear regressions, that is, ln (Di + 1) for the y-axis and ln (Mi) for the x-axis. The value of +1 was added because some size class bins were empty (Lykke 1998). The linear regressions were evaluated in terms of their slopes and y-axis intercepts (Gaugris, Vasicek & Van Rooyen 2012; Van der Merwe & Geldenhuys 2017). Strongly negative slopes typically represent regenerating populations (Condit et al. 1998), and a high y-axis intercept indicates many small individuals in the population. A y-axis intercept near to zero indicates that few small individuals are present (Everard, Midgley & Van Wyk 1994). Both slopes and y-axis intercepts of regressions were compared statistically across years by means of an analysis of covariance using GraphPad 7.02 (GraphPad Software, www.graphpad.com. Two additional mean values were calculated for the stem diameter distribution for each survey year: (1) the ‘centroid’ was calculated as the arithmetic mean of all stem diameters recorded; and (2) the ‘midpoint of stem diameter range’ was calculated as the mean of the upper and lower limits of the diameter range included in the size class distribution. Niklas et al. (2003) found a relationship between the position of the centroid and the diameter class range midpoint for a species. If a centroid value is smaller than the midpoint value, the population is young and growing. When the centroid value is larger than the midpoint value, the population is mature to old. In general, as the mean stem diameter increases in a population, self-thinning occurs and the density decreases, in such a way that an increasing part of the biomass is found in fewer individuals (Niklas et al. 2003; Silvertown & Charlesworth 2001). As self-thinning occurs the centroid of a population shifts towards the right within the size class distribution. Results Growth rate Over the period of 38 years, 6 of the 18 originally marked trees at the Dankbaar site died (33%) and three showed coppice regrowth following substantial dieback because of a fire (Table 1-A1). In 1978 the mean height of the 18 trees was 5.1 m and mean stem diameter 235.85 mm. In 2016 the mean height of the remaining nine originally marked trees that had not experienced dieback was 6.8 m and mean stem diameter 282.7 mm. The mean height for the six dead individuals at the onset of the study was 6.1 m and at the time of last measurement before death 6.8 m. A mean height increase of 60 mm/year was calculated. The 2016 survey measured the remaining 12 originally marked trees as well as 12 additional trees now present; thus, after 38 years, the population increased by six trees in the 1 ha plot. TABLE 1: Age calculations for Vachellia erioloba trees in the interior duneveld at the Dankbaar site using mean stem diameter of the 18 originally marked trees in 1978, the nine remaining live trees recorded in 2016 and the six dead trees last measured stem diameter when last surveyed as alive, derived by different methods. A linear regression of growth rate against stem diameter (Figure 4) indicated that growth rate decreases significantly as stem diameter increases, that is, as the individuals age (p = 0.022). The mean growth rate for the first 10-year period (1978–1988) was 2.99 mm/year, the second 10-year period (1988–1998) 1.01 mm/year and the third 10-year period (1998–2008) 2.84 mm/year, indicating a variable growth rate over time. Stem diameter growth rate per year varied widely between the nine remaining trees (0.91 mm – 4.86 mm/year), with a mean value of 2.5 mm/year. Depending on the choice of mean annual growth rate or application of the regression between stem circumference and age, the mean age of the six dead trees could be estimated at between 71 and 142 years (Table 1). FIGURE 4: Linear regression of stem diameter versus annual stem diameter growth rate for the 18 marked Vachellia erioloba trees at the Dankbaar site in the Kalahari Gemsbok National Park. Population structure At the Grootkolk site, 61 live and 14 dead individuals were recorded in the initial survey in 1979. The number of live individuals at Grootkolk in the Nossob riverbed gradually decreased over the survey period (Figure 5). Initially the size class distribution of the trees from 1979 to 1999 could visually be classified as Type IIIb curves; however, there was a gradual transition to a Type IV curve from 2008 to 2016 (Figures 3 and 5). In spite of the decreasing tree numbers and the change of curve shape, linear regressions on the density of individuals against the diameter class midpoint (Figure 1-A1) revealed no significant differences between the slopes (Table 2-A1) or intercepts (Table 3-A1) of the fitted lines. FIGURE 5: Size class distribution curves of Vachellia erioloba individuals at the Grootkolk site in the various survey years. (a) 1986; (b) 1984; (c) 1986; (d) 1989; (e) 1994; (f) 1996; (g) 1999; (g) 2008; (i) 2016. Using calculations from Steenkamp et al. (2008) for the interior duneveld and riverbed, the largest live tree at Grootkolk was aged between 195 and 212 years, while the relationship between stem circumference and radiocarbon age for the interior duneveld and riverbed combined aged the tree between 136 and 148 years (Table 2). When using the same two methods of calculation, most of the dead trees fell in the age categories 25–53, 37–71 and 49–89 years. TABLE 2: Age calculations, derived by different methods, for Vachellia erioloba trees at the Grootkolk site in the Nossob riverbed using mean stem diameter of the circumference class for the largest live tree; and for circumference size classes that recorded the majority of the dead trees. In all survey years the position of the centroid (arithmetic mean of all stem diameters recorded per site) was to the left of the midpoint of the size class distribution. However, the difference between these two values decreased over time and in 2016 the two values were nearly the same (Table 4-A1), indicating that the population was changing from a young and growing population to a mature to old population. Discussion Growth rate Marked V. erioloba tree stem diameters at Dankbaar in the interior duneveld increased at a mean rate of 2.50 mm/year (range 0.91 mm – 4.86 mm/year) over the monitored period. This mean annual diameter increase is less than the 3.6 mm/year established by Steenkamp et al. (2008) and substantially lower than the 12 mm – 14 mm/year previously reported for Acacia species growing in arid and semi-arid environments (Gourlay 1992). The Steenkamp et al. (2008) study used annual growth ring and radiocarbon ageing to determine growth rates. A mean growth rate was derived from individuals in the interior duneveld as well as in the Nossob riverbed. If only the individuals in the interior duneveld in the Steenkamp et al. (2008) study are considered, then the mean annual diameter increment is 3.2 mm/year. Although Steenkamp et al. (2008) concluded that the growth rate did not differ significantly between trees growing in the dry Nossob riverbed and those growing in the interior duneveld, growth rate varied substantially between individual trees. Variability in the growth rates of other tree species has been reported and the suggested underlying causes for this phenomenon include microsite conditions and/or genetic potential of individual trees (Stahle et al. 1996). Depending on the equation used to determine the age of the dead trees there could be a twofold or greater difference in calculated age at both the Dankbaar (Table 1) (71 years as opposed to 142 years) and Grootkolk (Table 2) (25 years as opposed to 127 years) sites. The high variability in estimated age at which 6 of the 18 marked trees died (48, 103, 163, 178, 181 and 191 years, using the mean annual diameter increase of 2.5 mm/year calculated for this current study) was unexpected. Moustakas et al. (2008), using a combination of aerial photographs and satellite images covering a period of 61 years to provide spatial data on mortality, found that V. erioloba mortality risk was greatest for trees of intermediate size. This does not seem to be the case for the Dankbaar population, where mortality was found to be predominantly among the large trees. At Grootkolk, the data seem to indicate that mortality was highest in smaller size classes and gradually declined for larger size classes. The density of trees at Grootkolk was low, and consequently self-thinning (Silvertown & Charlesworth 2001) is unlikely to have occurred. Population structure On the basis of the position of the centroid in relation to the midpoint of the stem diameter class range, the Grootkolk population could be considered as a population gradually changing from young and growing to mature to old. This is supported by a visual assessment of the size class distribution, which showed a gradual transition from a Type IIIb to a Type IV curve shape over the 38 years. The large peak of young individuals in size class 4 that was evident at the first survey gradually became smaller and shifted to size class 5, and a prominent second peak emerged at size class 8 at the last survey. At Grootkolk, few seedlings were recorded and few individuals appeared to be recruiting into the population. The low survival of size class 2, 3 and 4 individuals over the monitored period could possibly be ascribed to individuals progressing to larger size classes over the monitored period without being replenished. This is supported by the fact that size classes 5 and 8 had more individuals at the end of the monitoring period than when monitoring commenced. Overall, the survival of individuals in the larger size classes (e.g. size classes 6, 7 and 9) exceeded that in the smaller size classes (e.g. size classes 1, 2 and 3). In spite of the fact that the Grootkolk population was almost halved over the 38-year period, the slopes and intercepts of the size class distribution regressions showed no differences between years (Figure 1-A1). In their study, Steenkamp et al. (2008) reported an ageing size class distribution with poor recruitment at some sites in the Nossob riverbed. At other sites, many young individuals were found with only a few mature or dying trees in these populations. The Grootkolk site situated in the same riverbed illustrated an ageing population with few young individuals. It has been suggested that flooding of the Nossob River and/or consecutive years of above-average rainfall were necessary for regeneration of V. erioloba. Steenkamp et al. (2008) were, however, unable to find support to substantiate the need for flooding. Seed germination of V. erioloba is best after good rain events, and seedling survival depends on the availability of moisture during the next few seasons. Mortality rates during seedling establishment are generally higher than during the adult stages (Silvertown & Charlesworth 2001). Seedlings and saplings suffer high mortality in dry years (Barnes 2001; Theron, Van Rooyen & Van Rooyen 1985; Van Rooyen & Van Rooyen 1998; Van Rooyen et al. 1984). The high density of wildlife in the Nossob riverbed could also have contributed to the ageing of the Grootkolk population, with losses in young individuals because of trampling and browsing by herbivores and occasional toppling and breaking of larger trees by wind. Likewise, insects and rodents that impact on the seed pool may be more numerous in the riverbed habitat (Steenkamp et al. 2008). Furthermore, competition from a well-established grass layer may inhibit tree seedling establishment and growth in the riverbed (Steenkamp et al. 2008); however, this was not evident at the Grootkolk site. Fire also significantly impacts V. erioloba populations. After good rainfall, sufficient fuel accumulates and enables fire to spread through the landscape. In general, up to one-third of individuals in an area that experiences a fire are killed by the fire, with most of the individuals subsequently resprouting from the base (Van der Walt & Le Riche 1984). Large old trees, especially those with hollow stems, are most sensitive to fire (Milton & Dean 1995). During the study period, fire only occurred at the Dankbaar study site (large fires in 1995 and 2012). No tree mortality resulting from the fire was evident at the site; however, three of the individuals displayed coppiced growth subsequent to fire. Vachellia erioloba population dynamics differed at the two study sites over the last 38 years. Considering all tree mortalities and recruitments, the Dankbaar population can be considered a self-sustaining population with an increase of 12 new individuals (net increase of six trees) at the site over the 38-year period. The density (number of individuals per the 1 ha plot) of the Grootkolk population is declining and the size classes of the remaining individuals increasing as the population ages. Conclusions and management implications There is growing concern that mortality rates in V. erioloba, a keystone species of the Kalahari, are increasing as a result of threats such as increasing fire frequency, groundwater abstraction, competition from alien invasive species, diseases, harvesting and the effect of arboricides. This study in the KGNP provided an ideal opportunity to investigate V. erioloba dynamics in a natural environment, under different soil types and consequently water availability, where threats were negligible or absent. The different survey methods used at the two sites and alternative data analysis techniques have provided findings on the growth of individual trees as well as on population dynamics. Additionally, the study provided evidence that the population studied in the interior duneveld was self-sustaining and although the population in the Nossob riverbed was declining, new young populations were being established elsewhere in the riverbed (Steenkamp et al. 2008). The continued monitoring of the population at the two sites, where baseline data are already available, is highly recommended. However, monitoring sites should also be established at additional sites outside of park borders to incorporate sites in which various threats to this species can be investigated. It would also be imperative to determine growth rates in the riverbeds. Rainfall stations should be erected at all sites to enable the investigation of the influence of significant rainfall events and drought on seed germination, seedling survival and adult mortality. Acknowledgements The authors gratefully acknowledge the financial assistance of the University of Pretoria, Centre for Wildlife Management and the National Research Foundation (GUN 2053522, grant numbers 61277, 81089 and 91434). The numerous honours students in Wildlife Management and Botany are thanked for their assistance in the field. South African National Parks (SANParks) are thanked for the opportunity to conduct the research in the Kalahari Gemsbok National Park. The field surveys of 2016 were successfully conducted with the assistance of SANParks and South African Environmental Observation Network employees. The anonymous reviewers and editor are also thanked for their valuable contributions to this manuscript. Competing interests The authors declare that they have no financial or personal relationships that may have inappropriately influenced them in writing this article. Authors’ contributions N.v.R. donated the data set to South African Environmental Observation Network. The authors provided financial and/or institutional support and conducted numerous field surveys over the years. H.v.d.M. and M.W.v.R. compiled the manuscript with valuable input from the remaining co-authors. Funding information Funding for this study was provided by the University of Pretoria, Centre for Wildlife Management and the National Research Foundation. References Barnes, M.E., 2001, ‘Effects of large herbivores and fire on the regeneration of Acacia erioloba in Chobe National Park, Botswana’, African Journal of Ecology 39, 340–350. https://doi.org/10.1046/j.1365-2028.2001.00325.x Cody, M.L., 2000, ‘Slow-motion population dynamics in Mojave Desert perennial plants’, Journal of Vegetation Science 11, 351–358. https://doi.org/10.2307/3236627 Condit, R., Sukumar, R., Hubbel, S. & Foster, R., 1998, ‘Predicting population trends from size distributions: A direct test in a tropical tree community’, American Naturalist 152, 495–509. https://doi.org/10.1086/286186 Everard, D.A., Midgley, J.J. & Van Wyk, G.F., 1994, ‘Disturbance and the diversity of forests in Natal, South Africa: Lessons for their utilization’, Strelitzia 1, 275–285. Gaugris, J.Y. & Van Rooyen, M.W., 2007, ‘The structure and harvesting potential of the sand forest in Tshanini Game Reserve, South Africa’, South African Journal of Botany 73, 611–622. https://doi.org/10.1016/j.sajb.2007.06.004 Gaugris, J.Y., Vasicek, C.A. & Van Rooyen, M.W., 2012, ‘Herbivore and humans impact on woody species dynamics in Maputaland, South Africa’, Forestry 4, 497–512. https://doi.org/10.1093/forestry/cps046 Gourlay, I.D., 1992, ‘Age and growth rate determination of some African Acacia species’, MSc dissertation, University of Oxford, London. Government Notice no. 1602, 2016, ‘Notice of the list of protected tree species under the National Forests Act (1998) (Act no. 84 of 1998)’, Government Printer, Republic of South Africa. Jeltsch, F., Moloney, K.A. & Milton, S.J., 1999, ‘Detecting processes from snap-shot pattern – Lessons from tree spacing in the southern Kalahari’, Oikos 85, 451–466. https://doi.org/10.2307/3546695 Lawley, V., Parrot, L., Lewis, M., Sinclair, R. & Ostendorf, B., 2013, ‘Self-organization and complex dynamics of regenerating vegetation in an arid ecosystem: 82 years of recovery after grazing’, Journal of Arid Environments 88, 156–164. https://doi.org/10.1016/j.jaridenv.2012.08.014 Leistner, O.A., 1967, ‘The plant ecology of the southern Kalahari’, Memoirs of the Botanical Survey of South Africa 38, 1–172. Lykke, A.M., 1998, ‘Assessment of species composition change in savanna vegetation by means of woody plants’ size class distributions and local information’, Biodiversity and Conservation 7, 1261–1275. https://doi.org/10.1023/A:1008877819286 Malherbe, S.J., 1984, ‘The geology of the Kalahari Gemsbok National Park’, Koedoe 33–44. Menges, E.S., 2000, ‘Population viability analyses in plants: Challenges and opportunities’, Trends in Ecology & Evolution 15, 51–56. https://doi.org/10.1016/S0169-5347(99)01763-2 Midgley, J.J. & Bond, W.J., 2001, ‘A synthesis of the demography of African Acacias’, Journal of Tropical Ecology 17, 871–886. https://doi.org/10.1017/S026646740100164X Milton, S.J. & Dean, W.R.J., 1995, ‘How useful is the keystone species concept, and can it be applied to Acacia erioloba in the Kalahari Desert?’, Zeitschrift für Ökologie und Naturschutz 4, 147–156. Moustakas, A., Guenther, M., Wiegand, K., Mueller, K-H., Ward, D., Meyer, K.M. et al., 2006, ‘Long-term mortality patterns of the deep-rooted Acacia erioloba: The middle class shall die!’, Journal of Vegetation Science 17, 473–480. https://doi.org/10.1658/1100-9233(2006)17[473:LMPOTD]2.0.CO;2 Moustakas, A., Wiegand, K., Getzina, S., Ward, D., Meyer, K.M., Guenther, M. et al., 2008, ‘Spacing patterns of an Acacia tree in the Kalahari over a 61-year period: How clumped becomes regular and vice versa’, Acta oecologica 33, 355–364. https://doi.org/10.1016/j.actao.2008.01.008 Niklas, K.J., Midgley, J.J. & Rand, R.H., 2003, ‘Tree size frequency distributions, plant diversity, age and community disturbance’, Ecological Letters 6, 405–411. https://doi.org/10.1046/j.1461-0248.2003.00440.x Noy-Meir, I., 1973, ‘Desert ecosystems: Environment and producers’, Annual Review of Ecology and Systematics 4, 25–51. https://doi.org/10.1146/annurev.es.04.110173.000325 O’Connor, T.G., 1985, A synthesis of field experiments concerning the grass layer in the savanna regions of southern Africa, South African National Programmes Report No. 114, F.R.D., CSIR, Pretoria. Powell, E., 2005, What is happening to the camel thorns of Kathu?, Veld and Flora, Cape Town, pp. 34–35. Robertson, I. & Woodborne, S., 2002, Carbon isotopes confirm the competitive advantages of Prosopis over Acacia erioloba. Study of environmental change using isotope techniques, International Atomic Energy Agency, Vienna. IAEA-CSP-13/P. Schachtschneider, K. & February, E.C., 2013, ‘Impact of Prosopis invasion on a keystone tree species in the Kalahari Desert’, Plant Ecology 4, 597–605. https://doi.org/10.1007/s11258-013-0192-z Seymour, C., 2008, ‘Grass, rainfall and herbivores as determinants of Acacia erioloba (Meyer) recruitment in an African Savanna’, Plant Ecology 197, 131–138. https://doi.org/10.1007/s11258-007-9366-x Seymour, C.L. & Huyser, O., 2008, ‘Fire and the demography of camelthorn (Acacia erioloba Meyer) in the southern Kalahari—evidence for a bonfire effect?’, African Journal Ecology 46, 594–601. https://doi.org/10.1111/j.1365-2028.2007.00909.x Shadwell, E. & February, E., 2017, ‘Effects of groundwater abstraction on two keystone tree species in an arid savanna national park’, Peer J 5, e2923. https://doi.org/10.7717/peerj.2923 Silvertown, J. & Charlesworth, D., 2001, Introduction to plant population biology, 4th edn., Blackwell Science, Oxford, London. Slippers, B., Roux, J., Wingfield, M.J., Van der Walt, F.J.J., Jami, F., Mehl, J.W.M. et al., 2014, ‘Confronting the constraints of morphological taxonomy in the Botryosphaeriales’, Persoonia 33, 155–168. https://doi.org/10.3767/003158514X684780 Soil Classification Working Group, 1991, ‘Soil classification: A taxonomic system for South Africa’, Memoirs on the Agricultural Natural Resources of South Africa 15, 1–262. Stahle, D.W., Cleaveland, M.K., Haynes, G.A., Klimowicz, J., Ngwenya, P. & Mushove, P., 1996, ‘Preliminary tree-ring investigations of indigenous trees in Zimbabwe’, in P.T. Mushove, E.M. Shumba & F. Matose (eds.), Sustainable management of indigenous forests in the dry tropics. Proceedings of an International Conference in Kadoma, Zimbabwe, May 28–June 01, 1996, Zimbabwe Forestry Commission, Harare, pp. 57–61. Steenkamp, C.J., Vogel, J.C., Fuls, A., Van Rooyen, N. & Van Rooyen, M.W., 2008, ‘Age determination of Acacia erioloba trees in the Kalahari’, Journal of Arid Environments 72, 302–313. https://doi.org/10.1016/j.jaridenv.2007.07.015 Theron, G.K., Van Rooyen, N. & Van Rooyen, M.W., 1985, ‘Vegetation structure and vitality in the lower Kuiseb’, in B.J. Huntley (ed.), The Kuiseb environment: The development of a monitoring baseline, South African Scientific Programmes Report 106, pp. 81–91, Council for Scientific and Industrial Research, Pretoria. Van der Merwe, H. & Geldenhuys, C., 2017, ‘Proposed long-term monitoring protocol and applications for Aloidendron dichotomum populations’, South African Journal of Botany 109, 253–262. https://doi.org/10.1016/j.sajb.2017.01.008 Van der Walt, P.T &, Le Riche, E.A.N., 1984, ‘The influence of veld fire on an Acacia erioloba community in the Kalahari Gemsbok National Park’, Koedoe 103–106. Van Rooyen, T.H., 1984, ‘The soils of the Kalahari Gemsbok National Park’, Koedoe 27(suppl.) 45–61. https://doi.org/10.4102/koedoe.v27i2.568 Van Rooyen, M.W., Le Roux, A., Geldenhuys, C., Van Rooyen, N., Broodryk, N.L. & Van der Merwe, H., 2015, ‘Long-term vegetation dynamics (40yr) in the Succulent Karoo, South Africa: Effects of rainfall and grazing’, Applied Vegetation Science 18, 311–322. https://doi.org/10.1111/avsc.12150 Van Rooyen, N., Van Rensburg, D.J., Theron, G.K. & Bothma, J.d.P., 1984, ‘A preliminary report on the dynamics of the vegetation of the Kalahari Gemsbok National Park’, Koedoe 27(suppl.) 83–102. https://doi.org/10.4102/koedoe.v27i2.570 Van Rooyen, N. & Van Rooyen, M.W., 1998, ‘Vegetation of the South-western arid Kalahari: An overview’, Transactions of the Royal Society of South Africa 53, 113–140. https://doi.org/10.1080/00359199809520381 Van Rooyen, M.W., Van Rooyen, N., Bothma, J.d.P. & Van den Berg, H., 2008, ‘Landscapes in the Kalahari Gemsbok National Park, South Africa’, Koedoe 50, 99–112. https://doi.org/10.4102/koedoe.v50i1.154 Wiegand, K., Ward, D. & Saltz, D., 2005, ‘Multi-scale patterns and bush encroachment in an arid savanna with a shallow soil layer’, Journal of Vegetation Science 16, 311–320. https://doi.org/10.1111/j.1654-1103.2005.tb02369.x Zachariades, C., Hoffman, J.H. & Roberts, A.P., 2011, ‘Biological control of Mesquite (Prosopis species) (Fabaceae) in South Africa’, African Entomology 19, 402–415. https://doi.org/10.4001/003.019.0230 Appendix 1 TABLE 1-A1: Changes in stem diameter (mm) from 1978 to 2016 of 18 marked Vachellia erioloba trees at the Dankbaar site in the interior duneveld of the Kalahari Gemsbok National Park. TABLE 2-A1: Statistical comparison of slope values for linear regressions fitted on the density of individuals against the diameter class midpoint of live trees at Grootkolk in the Kalahari Gemsbok National Park using an analysis of covariance. TABLE 3-A1: Statistical comparison of intercept values for linear regressions fitted on the density of individuals against the diameter class midpoint of live trees at Grootkolk in the Kalahari Gemsbok National Park using an analysis of covariance. TABLE 4-A1: Calculated centroid (arithmetic mean of all stem diameters recorded per site in mm) and the midpoint (mm) of the size class distribution for each year surveyed at the Grootkolk site in the Nossob riverbed. FIGURE 1-A1: Linear regressions fitted on the ln (Mi) against ln (Di + 1) of live trees at Grootkolk in the Kalahari Gemsbok National Park in (a) 1979, (b) 1984, (c) 1986, (d) 1989, (e) 1994, (f) 1996, (g) 1999, (h) 2008 and (i) 2016, where Mi, midpoint of the size class and Di, density of individuals.