Landscape requirements of a primate population in a human-dominated environment
© Hoffman and O'Riain; licensee BioMed Central Ltd. 2012
Received: 1 November 2011
Accepted: 23 January 2012
Published: 23 January 2012
As urban and rural land development become widespread features of the global landscape so an understanding of the landscape requirements of displaced and isolated wildlife species becomes increasingly important for conservation planning. In the Cape Peninsula, South Africa, rapid human population growth, and the associated urban and rural land transformation, threatens the sustainability of the local chacma baboon population. Here we analyse spatial data collected from nine of the 12 extant troops to determine their population-level landscape requirements. We use hurdle models to ascertain the key landscape features influencing baboon occurrence and abundance patterns on two hierarchical spatial scales.
Both spatial scales produced similar results that were ecologically reliable and interpretable. The models indicated that baboons were more likely to occur, and be more abundant, at low altitudes, on steep slopes and in human-modified habitats. The combination of these landscape variables provides baboons with access to the best quality natural and anthropogenic food sources in close proximity to one another and suitable sleeping sites. Surface water did not emerge as an influential landscape feature presumably as the area is not water stressed.
The model results indicate that land development in the Cape Peninsula has pushed baboons into increasingly marginal natural habitat while simultaneously providing them with predictable and easily accessible food sources in human-modified habitats. The resultant spatial competition between humans and baboons explains the high levels of human-baboon conflict and further erosion of the remaining land fragments is predicted to exacerbate competition. This study demonstrates how the quantification of animal landscape requirements can provide a mechanism for identifying priority conservation areas at the human-wildlife interface.
The primary goal of ecologists is to understand the ecological factors that determine species distribution and abundance patterns . Furthermore, as urban expansion and rural land development become more widespread on the global landscape , so the understanding of the spatial requirements of species becomes increasingly important for conservation planning and management .
Habitat domination by humans , and the concomitant compression, fragmentation and conversion of primate habitats , are the driving forces behind human-primate conflict and one of the greatest threats to primate survival . The use of space has thus become a central theme in primate studies , with conservationists relying on patterns of habitat use and minimum resource requirements for the effective conservation and management of various primate populations . This is particularly true for those inhabiting small, isolated and fragmented habitats .
Within primates, baboons (genus Papio) are among the genera exhibiting the greatest degree of spatial overlap with humans . This success is attributed to their agility, dexterity, high levels of sociality and co-operation, combined with dietary and behavioural flexibility . Like the Macaca and Cercopithecus genera [9, 11], baboons can survive, and even appear to thrive, in human-modified habitats . However, as human populations expand and more land is developed, so the benefits afforded to baboons by habitat alteration are likely to be exceeded by the deleterious consequences of competition for space .
The chacma baboon population in the Cape Peninsula, South Africa provides one of the best examples of primate commensalism with humans. The earliest records of baboon and human co-existence in the Cape Peninsula date back to the 15th century with the arrival of Dutch settlers in South Africa . Over the last century the human population has grown rapidly , and all Cape Peninsula troops now have contact with humans, albeit to different degrees, in both residential and tourist-frequented areas. Furthermore, nearly half of the Cape Peninsula landscape has been transformed by a combination of urbanisation, farming and invasions by self-sown alien vegetation. Despite reductions in available land, geographical isolation caused by urban sprawl, historical extirpation of troops  and current high levels of human-induced injury and mortality , the local baboon population has persisted and has shown a steady increase in size over the last decade from 365 [1998; ] to 475 [2011; Beamish E.K., University of Cape Town, unpublished information].
Of concern to the sustainability of the Cape Peninsula baboon population are the continued expansion of the human population and the increasing spatial extent of the city of Cape Town, both of which have doubled over the last 30 years . In addition to exacerbating already high levels of human-baboon conflict, further encroachment of humans into natural areas of the Cape Peninsula may compromise baboon conservation in several ways. First, the spatial concentration of the baboon population may make it susceptible to intra-specific infectious diseases . Second, increasing overlap of baboons and humans could heighten the probability of bidirectional interspecies disease transmission [19, 20]. Third, the survivability and ecological role of future generations of baboons could be compromised if young baboons that grow up in troops heavily reliant on human food sources do not learn the necessary skills for finding and processing indigenous food . Finally, as contact between humans and baboons increases so baboons may become increasingly aggressive towards humans, as observed for vervet monkeys [Cercopithecus aethiops: ] and macaques [Macaca thibetana: ], and in extreme cases may be euthanized to protect public health and safety.
Scale 1 model results
Occurrence model coefficients
Abundance model coefficients
Scale 2 model results
Occurrence model coefficients
Abundance model coefficients
p (> |z|)
The Scale 1 and Scale 2 models found the same relationships between baboon abundance and the predictor variables (Tables 1 and 2). Habitat had the greatest influence on baboon abundance, followed by slope and then altitude. Among the topographic variables, the predicted abundance of baboons increased significantly with increasing slope and decreasing altitude. Relative to habitat, and compared to natural habitat, the predicted baboon abundance increased significantly in invasive alien vegetation and urban habitat. Distance to water and agricultural habitat had no significant influence on baboon abundance.
Post-hoc investigation results
Altitudinal vegetation characteristics
Shrubs 1-2 m
Shrubs 1-2 m
Shrubs 1-2 m
Shrubs 1-2 m
Shrubs > 2 m
Shrubs > 2 m
Shrubs > 2 m
Shrubs > 2 m
Low trees < 10 m
Shrubs > 2 m
Low trees < 10 m
Low trees < 10 m
Shrubs > 2 m
Low trees < 10 m
Low trees < 10 m
Shrubs > 2 m
Low trees < 10 m
Shrubs 1-2 m
Cumulative and total area of remaining natural habitat
Undeveloped area (km2)
Cumulative area (km2)
Undeveloped area (km2)
Cumulative area (km2)
Ecological reliability and interpretability of models
While the complexity of biological systems inhibits the ability of ecological models to reflect all reality, a model that suitably approximates the information contained in empirical data allows interesting inferences about ecology to be made . On account of their generalist nature , baboons are likely to have occurred throughout the Cape Peninsula prior to urbanisation, but with abundance being higher in more favoured habitat. The models reflect this pattern with the evaluation results and output maps showing that, despite being poorly correlated, the models provided qualitative predictions of baboon ranging patterns that can be considered accurate for four reasons. First, both the Scale 1 and Scale 2 models accurately predicted the current distribution of all studied troops. Second, both models predicted a higher abundance of baboons on land that is currently being used by troops that were not included in this study (n = 3; Figure 5) as well as historically recorded locations of extirpated troops (n = 3). Third, in only small, non-contiguous patches did the models predict baboon occurrence or abundance patterns that are not supported by either historic records or current baboon distribution patterns. Finally, the model predictions deviated most notably for the two troops most actively herded within their home ranges by baboon monitors . Here, both the Scale 1 and Scale 2 models predicted a lower probability of occurrence and a lower predicted abundance than we would expect given the troops ranging patterns [SK and DG troops: ]. This would suggest that the herding of baboons by monitors has impacted on their habitat use.
The low levels of positive spatial auto-correlation in the model residuals resulted from underlying landscape patterns, and may explain the average calibration of the models. Spatial correlation is almost always present in grid-datasets [31, 32] but as it does not bias regression coefficients  it does not affect the ecological interpretability of the models.
Key landscape features
It is possible, based on the inherent properties of the landscape, that the mountainous spine that runs the length of the Cape Peninsula has never provided sufficient food resources to support a large, spatially continuous baboon population. Baboons may have always been reliant on access to low land to obtain sufficient food, with high altitude areas acting as a demographic sink for the expanding population [e.g., ]. Thus, the rapid growth of the human population and extensive urbanisation over the last two centuries  has not only isolated the local baboons from all other populations but has also annexed most of the low lying and more productive foraging areas. That baboons have continued to persist in this environment despite these challenges is testament to their ability to modify their foraging behaviour and coexist with humans. This coexistence, however, has come at a severe cost with whole troops having been extirpated  and frequent cases of human-induced injury and mortality .
Both models indicated that baboons are more likely to occur, and be more abundant, at low altitudes, on steep slopes and in some human-modified habitats. These patterns are congruent with predictions based on baboon ecology (see below) as the combination of these variables provides baboons with access to food and sleeping sites - two resources critical to their survival.
Optimal foraging strategies for primates simultaneously maximise nutrient gain  and the efficient use of available time . Accordingly, patterns of primate distribution and abundance across the landscape can be explained primarily by the distribution of the most lucrative foraging sites . Baboon occurrence and abundance in the Cape Peninsula converged with the areas of the landscape that have the most profitable food sources, namely lower altitudes (see Table 3) and human-modified habitats.
Altitude and food
In the Cape Peninsula the benefits of foraging at low altitudes are threefold for baboons. Firstly, not only do the lower altitudes appear to contain larger plants than the rocky mountaintops - as revealed by landscape surveys - but the productivity of vegetation also decreases as altitude increases [A.R. Rebelo, South African National Biodiversity Institute, unpublished information]. Secondly, baboons gain access to high protein food resources along the coastline by consuming a variety of marine intertidal organisms. Not all local troops have access to the latter food source but troops that spend the majority of their time in the protein-poor indigenous vegetation  routinely include marine-food sources in their diet [Lewis M.C., University of Cape Town, unpublished information]. Thirdly, most human-modified habitat occurs at lower altitudes and offers highly concentrated and predictable food resources.
Human-modified habitats and food
Anthropogenic habitat alteration can dramatically affect the quality, availability and distribution of food resources and the addition of anthropogenic food sources into primate diets can have a positive effect on both their abundance and fecundity . In addition to the food sources available in urban (e.g., fruit trees in gardens, garbage in refuse bins, food items in houses) and agricultural habitats (e.g., grapes in vineyards, pine nuts in Pinus plantations, ostrich feed in livestock farms), humans in the Cape Peninsula have introduced many species of invasive alien plants [e.g., Pinus, Acacia and Eucalyptus spp., ] which have both higher seed production and standing biomass than indigenous vegetation . Human-modified habitats thus offer abundant, accessible and calorie-rich food sources that baboons favour over the low quality forage of local indigenous vegetation [29, 42]. That the models detected, on average, that baboons preferred human-modified habitats to natural habitat is thus unsurprising. There were, however, some interesting exceptions to this trend.
Baboons were less likely to occur in urban than in natural habitat despite the abundance of high quality food sources available in both houses and gardens. This is almost certainly a consequence of conflict with humans, with baboons suffering from harassment, injury and mortality when foraging in urban habitat . However, in time baboons are able to adapt to profitable foraging conditions by improving their raiding success while simultaneously minimising the costs associated with foraging in high risk habitats . Indeed baboons in the Cape Peninsula manage to mitigate against human threats while maximizing nutrient gain by spending minimal time raiding in urban habitat, acquiring human food quickly and returning thereafter to the relative safety of other habitats . This raiding strategy would explain the low probability of occurrence predicted for urban habitat.
However, an additional explanation for these results may lie in the statistical procedure of the modelling process. In both the Scale 1 and Scale 2 models the proportion of urban habitat used was much lower than what was available. When the amount of available urban habitat was decreased relative to the habitat use values (Scale 2 model) the negative effect of urban habitat on baboon occurrence was reduced. This suggests a sensitivity of the occurrence models to large discrepancies between use and availability values. With these discrepancies controlled, as they were in the abundance models, urban habitat was found to be favourable to natural habitat at both spatial scales.
A similar pattern was found for agricultural habitat where the occurrence models determined that baboons preferred this habitat to natural habitat, but the abundance models found the preferences for the habitats to be similar. There are two possible - perhaps interacting - explanations for this pattern. First, some crops (e.g., vineyards) have a distinctly seasonal growth cycle and are consequently not used consistently by baboons on an annual basis . Second, agricultural habitats are the source of income for farmers who, across sub-Saharan Africa consider baboons to be pests, capable of more crop damage than any other primates [43, 44] or indeed any other wildlife species . Of the types of agricultural habitat in the Cape Peninsula, baboons are tolerated in pine and eucalyptus plantations, but vineyard owners actively chase baboons when the vines are in fruit , and livestock farmers routinely chase baboons off their property throughout the year . Seasonal differences in crop use and farm policing would not affect the presence-absence analyses of the occurrence models, but would reduce the significance of the overall patterns of baboon abundance in agricultural habitat.
Invasive alien vegetation
Both the occurrence and abundance models revealed that baboons found invasive alien vegetation to be significantly preferential to natural habitat. This result is unsurprising given the 3-10 fold increase in above-ground biomass associated with invasive alien vegetation  and the resultant preference for foraging in this habitat . Furthermore, unlike urban habitat, baboons are able to exploit invasive alien vegetation without the cost of human harassment.
While food resources play a crucial role in determining primate spatial distributions, their availability to the animals is constrained by their proximity to other critical resources. In affording baboons safety from predators  and providing them with suitable vantage points for area surveillance , sleeping-sites fundamentally affect baboon ranging patterns and dictate the intensity to which they use the landscape . Baboons use a variety of sleeping-sites including cliffs, trees, and caves . In the Cape Peninsula baboons sleep primarily in trees and on cliffs [[29, 30]].
The baboons' use of trees as sleeping-sites provides an additional explanation for the preference shown for invasive alien vegetation and agricultural habitat over natural habitat. Natural habitat in the Cape Peninsula is characterised by shrublands, grasslands and low trees . By contrast both self-sown invasive alien vegetation and cultivated plantation trees are suitably sized for baboon sleeping-site requirements [29, 42]. Furthermore, plantation trees have also been cultivated alongside vineyards and urban areas and, because of the high levels of disturbance at the urban/natural habitat interface and the significant source of alien propagules presented by suburban gardens , self-sown alien plants tend to invade and establish in close proximity to urban habitat. Thus the spatial distribution of tall trees in the Cape Peninsula provides baboons with suitable sleeping-sites in close proximity to favoured foraging areas.
The importance of cliff sleeping-sites was also detected by the models. Both models indicated that baboons are more likely to occur, and be more abundant, on steep slopes - the inaccessibility of which provides them with a safe night-time refuge. The coincident preference for steep slopes and low altitudes once again represents the importance of proximity of sleeping-sites to favoured foraging areas. In the Cape Peninsula, slope steadily increases from the 100 m contour line (Figure 4). Because urban development is constrained by the exposure and inaccessibility of high altitudes and steep slopes  the spatial extent of urban habitat is restricted to the flat land below the 80 m contour line. Consequently steep cliffs, suitable as baboon sleeping-sites, occur directly above favourable urban habitat in many areas.
Given that water is a critical resource for baboons it is surprising that there was a significantly higher probability of baboons occurring far from permanent surface water sources than near to them. An explanation for this pattern emerges when results are considered in context with the hydrological attributes of the Cape Peninsula. With permanently flowing surface waters, the presence of freshwater wetlands and vegetation prone to seasonal water-logging , the Cape Peninsula is not a water-stressed environment. These factors explained the relative lack of importance of water in determining the ranging patterns of one of the local troops [TK troop: ]. For this troop the high water content of vegetation explained the lack of 'drinking sessions' [sensu] where many baboons converge at a waterhole and drink simultaneously. These explanations seem equally applicable to all troops in the Cape Peninsula, especially given the results of the abundance models that indicated no significant relationship between baboon land use and surface water. Thus, rather than revealing an interesting ecological phenomenon where animals avoid water, the model results should rather be interpreted as an indication that permanent surface water need not be considered as a key landscape feature for baboons in the Cape Peninsula during years of good rainfall. However, permanent water sources may indeed become a good predictor of baboon occurrence during years of drought or below par rainfall, particularly during the dry summer season.
Currently, the most widely used method of baboon management is the employment of baboon monitors who attempt to reduce levels of human-baboon conflict by preventing troops from crossing the urban edge . Because of the spatial attributes of urbanised areas, monitors must typically herd baboons away from favourable low lying land and up the mountain, into increasingly marginal habitat (Figure 5). That monitored troops still suffer extraordinary high levels of human-induced injury  is not attributable to atypical or errant baboon behaviour, but rather indicative of the intensity of the competition between baboons and humans for low lying land and the associated high quality natural and anthropogenic food resources. The inevitable persistence of this competition through time, will forever compromise the ability to effectively manage baboons at the interface of natural and urban habitats.
The implications of extensive land development for the future of baboon management and conservation are equally concerning. Currently, minimal amounts of the land considered to be ecologically suitable for baboons (probability of occurrence > 0.5) remain undeveloped. This is most pronounced in the northern half of the Cape Peninsula (Figure 5) where suitable natural habitat comprises non-contiguous islands made inaccessible to baboons by the sea of urbanisation surrounding them. The landscape of the southern half of the Cape Peninsula, where three-quarters of the current baboon population range, holds more promise with continuous stretches of ecologically suitable natural habitat still available. Of the natural habitat available in the southern Cape Peninsula, 87% is conserved as part of the Table Mountain National Park, with 13% potentially subject to urban or agricultural transformation. Transformation of this land will force baboons farther into the increasingly marginalised habitat of higher altitudes and will ultimately exacerbate levels of human-baboon conflict, with inevitable deleterious consequences for baboons. If attempts to conserve the Cape Peninsula baboon population are to be effective, then mitigating against the development of the remaining natural fragments of the landscape must be addressed urgently.
Conserving this undeveloped land will require the involvement of various landowners , including both private owners (62% of land) and national and provincial government (38%). The local government departments responsible for approving developments have already accepted the GPS data from this research and included then within new zoning schemes. While this will not automatically prevent further erosion of natural habitats, it will trigger baboon-specific Environmental Management Plans and initiate assessments of the possible impacts of such development on the affected troop(s), with input from interested and affected parties. However, because most undeveloped land is owned by private parties, it is vital that the responsibility of landscape conservation is not left entirely to environmental authorities. Rather, a great deal of effort should be placed on encouraging and or providing incentives for private landowners to participate in the conservation process. This can be achieved using a variety of methods including one or a mix of voluntary, property-based, price-based or regulatory mechanisms [51, 52]. Ultimately, however, the effectiveness and sustainability of any landscape conservation strategy relies on the a combination of 'top-down rigour' and 'bottom-up participation' , where partnerships between management authorities, non-governmental organisations and private parties will ensure that land development plans are not only ecologically sustainable, but socially sustainable [54–56]. In line with this, a baboon liaison group was established in the Cape Peninsula in 2010 to serve as an intermediary between baboon management authorities, researchers and local communities. In providing a platform for communication, the establishment of this liaison group will improve the chances of the local baboon population being sustainable on the long-term.
This study highlights the complexities of wildlife management and conservation at the interface of natural and human-modified habitats. This is particularly true for wildlife whose landscape requirements are concurrent with those of humans, and whose ecological flexibility allows them to thrive in human-modified habitats - qualities which may increase their probability of experiencing conflict with humans. However, by enhancing our understanding of the fundamental drivers of human-wildlife conflict, the quantification of animal landscape requirements can support and validate wildlife conservation efforts. Furthermore, an understanding of animal spatial ecology can assist in discriminating between the relative importance of landscape quality and landscape quantity, and can provide a mechanism for identifying priority conservation areas at the human/wildlife interface.
Located at the south-western most point of the African continent (latitude 33°55'-34°21' S; longitude 18°25'-18°28 E), the Cape Peninsula (Figure 1) comprises a combination of natural and human-modified habitats bounded by the Atlantic Ocean. The topography is characterised by the Peninsula Mountain chain and has an altitudinal range of 0-1100 m. Lower elevations are predominantly urbanised, mid-elevations are used for agriculture (vineyards, plantations and livestock farming) and higher elevations are almost exclusively indigenous fynbos vegetation. Fynbos is a species-rich but nutrient poor, sclerophyllous shrubland that is a key component of the Cape Floristic Region . More than half of the Cape Peninsula remains undeveloped and is conserved within the Table Mountain National Park. Chacma baboons are the only non-human primates present, and no natural baboon predators remain.
Study population and data collection
At the time of data collection the Cape Peninsula baboon population comprised 12 troops ranging in size from 16-115 baboons . Urban habitat and neighbouring troops served as the only barriers to troop movement and troops were able to range freely in approximately 250 km2 of natural habitat. All nine study troops were accustomed to the presence of people and were habituated to close (≤ 10 m from baboons) observation at the start of the study. Four troops were managed by baboon monitors - people employed by the local management authorities to minimise human-baboon conflict by preventing troops from entering urban areas . Troop home ranges show low levels of overlap (mean ± 1SE: 7.3 ± 4.9%)  and consequently troops seldom interact with each other. However, sub-adult and adult male baboons remain able to disperse to neighbouring troops in accordance with typical baboon behaviour.
Between 2006 and 2009 we collected one full year of spatial data for each of the study troops. Logistical and financial constraints prevented us from sampling all troops simultaneously, but there was no interannual variation in mean rainfall (Kruskal-Wallis test: H(5, n = 2191) = 6.6317, p = 0.250), mean maximum (ANOVA: F5,2185 = 0.8219, df = 2185, p = 0.0534) or mean minimum temperatures (ANOVA: F5,2185 = 1.5676, df = 2185, p = 0.166) across the study years. We recorded Global Positioning System (GPS) data points for each troop using (a) handheld devices (Garmin eTrex) operated by field researchers (n = 5 troops), (b) tracking collars (n = 3 troops) and (c) a combination of both methods (n = 1 troop). Field researchers recorded the GPS location of the centre point of the troop (visually estimated) at 20-minute intervals between sunrise and sunset for a mean of 109 days (± 28 days 1SE, range: 71-170 days, n = 6 troops) per troop. Tracking collars recorded the GPS point of a single troop member at 3-hourly intervals between sunrise and sunset for a mean of 302 days (± 54 days 1SE, range: 247-334 days, n = 3 troops) per troop. The terrain within the study area was easily traversable on foot and visibility of baboons within all habitat types was excellent. Only GPS data points that had an estimated level of accuracy of ≤ 10 m were included in our analyses. A total of 24,618 GPS data points was recorded for the population, with a mean of 2735 ± 768 GPS data points 1SE (range: 1668 - 5018, n = 9 troops) recorded per troop.
Model study areas and datasets
To produce tabular datasets for the models we assigned a matrix of grid cells to the Scale 1 and Scale 2 study areas. Each grid cell was 0.023 km2 in area (150 m × 150 m) set to reduce the likelihood of any troop being spread through more than one cell simultaneously . Consequently, the total area of each cell was sufficient to encompass the mean troop spread of the largest troop (mean ± 1SE: 0.021 ± 0.011 km2, n = 5 spreads). We merged GPS data from the baboon population to the Scale 1 and Scale 2 grids. The GPS data comprised 1000 GPS data points from each of the nine study troops, randomly selected to control for intertroop differences in sample sizes and sampling regimes, and the effects of group size and season on baboon ranging patterns. We pooled the GPS data together to generate a population-level dataset (Figure 1). For both study areas we determined a use value for each entered grid cell by counting the number of GPS data points within it. We assigned use-values of zero to non-entered cells.
Cell count details and predictor variable attributes of the model study areas
Scale 1 study area
Scale 2 study area
Cell count details
Mean ± SEM
Mean ± SEM
155.1 ± 2.4 m
0 - 1069.3 m
154.5 ± 2.5 m
0 - 911.81 m
0 - 61.8°
0 - 57.4°
0.8 ± 0.02 km
0 - 10.2 km
1.1 ± 0.03 km
0 - 10.2 km
We analysed cell use as a function of the following predictor variables for both the Scale 1 and Scale 2 models: altitude (m), slope (°), the distance to permanent surface water sources (km) and habitat. We assigned explanatory variables (altitude, slope, distance to permanent surface water sources and habitat) to each cell in the Cape Peninsula grid matrix using the following methods. We calculated Altitude using a 30 m digital elevation model [DEM; Environmental Systems Research Institute. 1998. CSDGM FGDC Metadata DTD 3.0.0 19981217]. We used the same DEM to determine the slope of each cell, and calculated distances to permanent surface water sources using a shapefile [glcrveg; South African National Parks, unpublished information] that details drainage systems in the region. We used GIS maps containing landscape information specific to the Cape Peninsula (details below) to categorise habitat for each grid cell. In addition, we used observer records of location-specific habitats, and information gleaned from digitisation of the Cape Peninsula using Google Earth imagery to confirm and identify habitat in areas where the extent and or the detail of the GIS layers was insufficient or inaccurate in its descriptions. The percentage cover of each habitat within every cell was calculated using the Intersect Function of the Geoprocessing Wizard in ArcView 3.3. We converted the percentage values to categorical variables, and assigned habitat categories based on the dominant habitat (> 50% of cover) within each cell. We categorised habitat as natural, urban, agricultural and invasive alien vegetation. Natural habitat, categorised using a shapefile from , included indigenous vegetation, rocky shores and beaches. Urban habitat consisted of all urban areas delineated in the City of Cape Town's Generalised Zoning shapefile, cells with human-made structures such as buildings, gardens and grass patches adjacent to buildings, and cells dominated in cover by roads and sports fields. Agricultural habitat included plantations - delineated using a shapefile [lease_2006_06_28; South African National Parks, unpublished information] - and vineyards and an ostrich farm that were mapped digitally using Google Earth. We also used Google Earth to digitally map the range and extent of invasive alien vegetation (Pinus, Acacia and Eucalyptus spp.).
Multicollinearity of predictor variables
The study areas differed in their overall composition of habitats but were similar in their topographic profiles (Table 5). We used both broad-scale and fine-scale habitat variables. The broad-scale variable (Broad Habitat) categorised habitat as being natural or human-modified. The fine-scale variable (Fine Habitat) included the broad-scale natural category and the human-modified sub-categories, namely urban habitat, agricultural habitat and invasive alien vegetation. To determine which of these habitat variables were most suitable for inclusion in the final models for each dataset, we evaluated each habitat variable in turn at the model building stage using the methods described under 'Model selection and evaluation' below. Once the final models were selected we evaluated each predictor in terms of its overall contribution to each respective model . We used the habitat variable of 'natural' as the intercept category for both models.
Modelling algorithms and model fitting
For all datasets we used hurdle models  to analyse cell use as a function of the predictor variables. Potts and Elith  found that relative to four other regression models (Poisson, negative binomial, quasi-Poisson and the zero-inflated Poisson) the hurdle model had the greatest predictive performance when assessing the relationship between the abundance of an organism and its environment. Hurdle models are modified count models that separate data into two parts: one containing zero values and one containing positive counts . As such, hurdles account for two ecological processes: the first is the process that causes an animal to be present at a site (occurrence) , and the second is the process that influences the numbers of animals found at a site, given that they occur there (abundance). Hurdle models model occurrence using binary (presence/absence) models with a binomial probability and model the abundance (positive counts) using zero-truncated count models . Quasi-likelihood removes the effect of zero-inflation in the binary models and the effect of over-dispersion in the zero-truncated models .
Count data is typically modelled using the Poisson distribution, which assumes equality between means and variances . However, actual count data are often over-dispersed (mean < variance) relative to the Poisson distribution [e.g., [67, 61]]. In this instance it becomes more appropriate to use the negative binomial distribution, which allows for a quadratic relationship between the mean and the variance [e.g., [66, 67]]. On account of the over-dispersion in both the Scale 1 and Scale 2 datasets (Table 5) we fitted the positive count models with negative binomial distributions.
We ran occurrence and abundance models for the Scale 1 and Scale 2 study areas. For the occurrence models we used a binomial distribution with logit link, and for the abundance models we used a negative binomial distribution with log link to ensure that the predicted values were always positive . We conducted all statistical analyses using the R language and environment, an integrated software suite for statistical computing . We fitted all models using the R package 'pscl' .
Model selection and evaluation
AIC values for candidate models
Scale 1 candidate models
Scale 2 candidate models
FH and ALT and SL and WAT*
FH and ALT and SL and WAT
BH and ALT and SL and WAT
BH and ALT and SL and WAT
ALT and SL and WAT
ALT and SL and WAT
FH and ALT and SL
FH and ALT and SL
FH and ALT
FH and ALT
FH and SL
BH and ALT and SL
FH and SL
BH and ALT and SL
ALT and SL
BH and ALT
ALT and SL
BH and ALT
BH and SL
BH and SL
ALT:SL (interaction term)
ALT:SL (interaction term)
Spatial autocorrelation deals with a lack of independence of data points and measures the degree to which a variable is correlated to itself in space . This phenomenon is pervasive in ecological datasets  and can stem from movement patterns of the study subject or underlying patterns of the landscape. Spatial autocorrelation can be problematic in analyses as it can lead to Type 1 statistical errors (false positives) and can result in inflated probabilities for predictor variables .
The random selection of GPS data points for inclusion in the model datasets accounted for any spatial autocorrelation attributable to animal movement patterns. During the modelling process, once the models are fitted to the data, predictor variables should account for any autocorrelation caused by landscape patterns. If this is not the case, then spatial autocorrelation should be evident in the model residuals . We used GeoDa 0.9.5-i [GeoDa Center for Geospatial Analysis and Computation, Arizona] to test for spatial autocorrelation in the Scale 1 and Scale 2 model residuals using Monte Carlo simulation (999 permutations) of Moran's I. Moran's I ranges from -1.0 - 1.0, with non-zeroes indicating that the abundance values produced for spatially connected grid cells are either more similar (positive autocorrelation) or more different (negative autocorrelation) than would be expected given a random association among the cells . We calculated Moran's I using a weight matrix defined by k-nearest neighbours, with the value for k defined by the number of cells within a 1 km radius from each cell (k = 224).
Model selection and model checking
To enhance our ability to interpret the model results in the context of the Cape Peninsula landscape, we conducted two post-hoc investigations to establish the spatial relationships between altitude and slope, and altitude and vegetation biomass. First, we used ArcView 3.3 to delineate the Cape Peninsula landscape into 100 m altitudinal belts increasing from sea level to 1100 m (regional maximum), and to calculate the mean (± 1SE) slope of each belt. Second, to analyse the relationship between altitude and vegetation biomass, we surveyed a section of the Cape Peninsula landscape. We assessed changes in biomass along three altitudinal transects running from sea level to 600 m above sea level. At 100 m intervals along each transect we visually determined the growth form and canopy cover of the dominant stratum within 10 × 5 m quadrates [sensu]. These particular transects were selected because they were located in the only area of the Cape Peninsula that contained all of the following three attributes: (1) they included an extensive and traversable altitudinal range relative to other regions of the Cape Peninsula, (2) they covered an area of land stretching from sea level to mountain top that was undeveloped and dominated by natural habitat, and (3) they fell within a large enough area to allow three replicate surveys to be conducted while controlling for geology, hydrology and invasive alien vegetation.
Finally, using the Scale 1 model results, we calculated the absolute and cumulative areas of natural habitat remaining in the Cape Peninsula for each level of probability and each category of predicted abundance.
This research was funded by the Table Mountain Fund (WWF-SA), Table Mountain National Parks, the National Research Foundation and Back to Africa. The research complied with protocols approved by the ethics committees of the University of Cape Town, South African National Parks and the Society for the Prevention of Cruelty to Animals, and adhered to South African legal requirements. We thank South African National Parks for the permission to conduct fieldwork in the Table Mountain National Park and for their assistance with data collection and analysis. We would like to thank Esme Beamish, Alta de Vos, Bentley Kaplan and Matthew Lewis of the Baboon Research Unit at the University of Cape Town for their field assistance and insightful discussion, Angela van Doorn for her contributions to the spatial dataset, and Andrew King for his comments on earlier versions of this manuscript. We would like to thank Anthony Rebelo for providing botanical data. We are grateful to Anneli Hardy, Birgit Erni and Christien Thiart from the Department of Statistics at the University of Cape Town for their extensive help with the statistical analyses, and to Hamish Currie, David Zimmerman, Mark Dittberner and Justin Benade for their veterinary assistance. Finally, without the volunteer efforts of our long-term field assistants (Caroline Booth, Nelly Boyer, Shamini Bundell, Shahrina Chowdhury, Caroline Feinberg, Cody Freas, Alistair Fyfe, Mandy Gibson, Amy Harrington, Tracy Herigstad, Jolle Jolles, Dana Kirschmann, Susan Lin, Cody Massing, Liliana Mendel, Ania Paluch, Lindsay Patterson, Darby Proctor, Rachel Schill, Reinhardt Scholtz, Nitin Sekar, Leslie Steele, Kathryn Tarr, Nadia Timol, Denise Thompson, Simon Wood, Nicolina van der Merwe) and short-term field assistants (Jonathan Aronson, Brionie Benchley, Chloe Botha, Charine Collins, Stacey Jordaan, Abby Joustra, Anne Ketley, Megan Laird, Christopher Mills, Glenn Moncrieff Jacqui Stephenson, Alicia Thomas, Liezel Tolmay) this study would not have been possible.
- McLoughlin PD, Morris DW, Fortin D, Vander Wal E, Contasti AL: Considering ecological dynamics in resource selection functions. J Anim Ecol. 2010, 79: 4-12. 10.1111/j.1365-2656.2009.01613.x.View ArticlePubMedGoogle Scholar
- Western D: Human-modified ecosystems and future evolution. P Natl Acad Sci USA. 2001, 98: 5458-5465. 10.1073/pnas.101093598.View ArticleGoogle Scholar
- Beckmann JP, Berger J: Rapid ecological and behavioural changes in carnivores: the responses of black bears (Ursus americanus) to altered food. J Zool, London. 2003, 261: 207-212. 10.1017/S0952836903004126.View ArticleGoogle Scholar
- Paterson JD, Wallis J: Special Topics in Primatology. Volume 4: Commensalism and Conflict: The Human-Primate Interface. 2005, Norman, Oklahoma: American Society of PrimatologistsGoogle Scholar
- Strum SC: The development of primate raiding: implications for management and conservation. Int J Primatol. 2010, 31: 133-156. 10.1007/s10764-009-9387-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Laurance WF, Lovejoy TE, Vasconcelos HL, Bruna EM, Didham RK, Stouffer PC, Gascon C, Bierregaard RO, Laurance SG, Sampaio E: Ecosystem decay of Amazonian forest fragments: A 22-year investigation. Conserv Biol. 2002, 16: 605-618. 10.1046/j.1523-1739.2002.01025.x.View ArticleGoogle Scholar
- Arrowood BHC, Trevest A, Mathews NE: Determinants of day-range length in the black howler monkey at Lamanai, Belize. J Trop Ecol. 2003, 19: 591-594. 10.1017/S026646740300364X.View ArticleGoogle Scholar
- Robbins RM, McNeilage A: Home range and frugivory patterns of mountain gorillas in Bwindi Impenetrable National Park, Uganda. Int J Primatol. 2003, 24: 467-491. 10.1023/A:1023741328338.View ArticleGoogle Scholar
- Hill CM: People, crops, and primates: A conflict of interests. Commensalism and Conflict: The Human-Primate Interface. Edited by: Paterson JD, Wallis J. 2005, Norman, Oklahoma: The American Society of Primatologists, 40-59.Google Scholar
- Swedell L: African Papionins: diversity of social organization and ecological flexibility. Primates in Perspective. Edited by: Campbell CJ, Fuentes A, MacKinnon KC, Panger M, Bearder Slangkop. 2011, New York: Oxford University Press, 241-277. 2Google Scholar
- Lee PC, Brennan EJ, Else JG, Altmann J: Ecology and behaviour of vervet monkeys in a tourist lodge habitat. Primate Ecology and Conservation. Edited by: Else JG, Lee PC. 1986, Cambridge: Cambridge University Press, 229-236.Google Scholar
- Strier KB: Primate Behavioral Ecology. 2007, Boston: Pearson Education Inc, 3Google Scholar
- Goodwin AJH: Jan van Riebeeck and the Hottentots 1652-1662. S Afr Archaeol Bull. 1952, 7: 2-6. 10.2307/3887530.View ArticleGoogle Scholar
- Sinclair-Smith K: The expansion of the urban Cape Town. Metropolitan Spatial Planning Branch. 2009, City of Cape Town: Spatial Planning and Urban Design DepartmentGoogle Scholar
- Skead CJ: Historical Mammal Incidence in the Cape Province. 1980, Cape Town: Department of Nature and Environmental Conservation, 1:Google Scholar
- Beamish EK: Causes and consequences of mortality and mutilation in the Cape Peninsula baboon population, South Africa. MSc thesis. 2010, University of Cape Town, Zoology DepartmentGoogle Scholar
- Kansky R, Gaynor D: Baboon Management Strategy for the Cape Peninsula. Final Report. 2000, Table Mountain Fund Project number ZA 568, Cape Town, South AfricaGoogle Scholar
- Wrangham RW: Artificial feeding of chimpanzees and baboons in their natural habitat. Anim Behav. 1974, 22: 83-94. 10.1016/S0003-3472(74)80056-4.View ArticleGoogle Scholar
- Ravasi DFC: Gastrointestinal parasite infections in chacma baboons (Papio h. ursinus) of the Cape Peninsula, South Africa: the influence of individual, troop and anthropogenic factors. PhD thesis. 2009, University of Cape Town, Zoology DepartmentGoogle Scholar
- Drewe J, O'Riain MJ, Beamish EK, Currie H, Parsons S: A survey of infections transmissible between baboons and humans in Cape Town, South Africa. Emerg Infect Dis.Google Scholar
- Asquith PJ: Provisioning and the study of free-ranging primates: History, effect, and prospects. Yearb Phys Anthropol. 1989, 32: 129-158. 10.1002/ajpa.1330320507.View ArticleGoogle Scholar
- Zhao Q-K: Birth timing shift with altitude and its implications for Macaca thibetana at Mt.Emei, China. Oecolog Montana. 1994, 3: 24-26.Google Scholar
- Johnson DH: The comparison of usage and availability measurements for evaluating resource preference. Ecology. 1980, 61: 65-71. 10.2307/1937156.View ArticleGoogle Scholar
- Barton RA, Whiten A, Strum SC, Byrne RW, Simpson AJ: Habitat use and resource availability in baboons. Anim Behav. 1992, 43: 831-844. 10.1016/S0003-3472(05)80206-4.View ArticleGoogle Scholar
- Crook JH, Aldrich-Blake P: Ecological and behavioural contrasts between sympatric ground-dwelling primates in Ethiopia. Folia primatol. 1968, 8: 192-227. 10.1159/000155154.View ArticlePubMedGoogle Scholar
- Whiten A, Byrne RW, Henzi SP: The behavioral ecology of mountain baboons. Int J Primatol. 1987, 8: 367-389. 10.1007/BF02737389.View ArticleGoogle Scholar
- Altmann SA, Altmann J: Baboon Ecology: African Field Research. 1970, Chicago and London: University of Chicago PressGoogle Scholar
- Burnham KP, Anderson DR: Model selection and multimodel inference: A practical information-theoretic approach. 2nd edition. 2002, New York: Springer-VerlagGoogle Scholar
- van Doorn A: The interface between socioecology and management of chacma baboons (Papio ursinus) in the Cape Peninsula, South Africa. PhD thesis. 2009, University of Cape Town, Zoology DepartmentGoogle Scholar
- Hoffman TS: The spatial ecology of chacma baboons (Papio ursinus) in the Cape Peninsula, South Africa: Towards improved management and conservation strategies. PhD thesis. 2011, University of Cape Town, Zoology DepartmentGoogle Scholar
- Rahbek C, Graves GR: Multiscale assessment of patterns of avian species richness. P Natl Acad Sci USA. 2001, 98: 4534-4539. 10.1073/pnas.071034898.View ArticleGoogle Scholar
- Diniz-Filho JAF, Bini LM, Hawkins BA: Spatial autocorrelation and red herrings in geographical ecology. Global Ecol Biogeogr. 2003, 12: 53-64. 10.1046/j.1466-822X.2003.00322.x.View ArticleGoogle Scholar
- Hawkins BA, Diniz-Filho JAF, Bini LM, De Marco P, Blackburn TM: Red herrings revisited: spatial autocorrelation and parameter estimation in geographical ecology. Ecography. 2007, 30: 375-384.View ArticleGoogle Scholar
- Henzi SP, Dyson ML, Deenik A: The relationship between altitude and group size in mountain baboons (Papio cynocephalus ursinus). Int J Primatol. 1990, 11: 319-325/. 10.1007/BF02193003.View ArticleGoogle Scholar
- Oates JF: Food distribution and foraging behavior. Primate Societies. Edited by: Smuts BB, Cheney DL, Seyfarth RM, Wrangham RW, Struhsaker TT. 1987, Chicago: University of Chicago Press, 197-209.Google Scholar
- Dunbar RIM, Korstjens AH, Lehmann J: Time as an ecological constraint. Biol Rev. 2009, 84: 413-429.View ArticlePubMedGoogle Scholar
- Clutton-Brock TT, Harvey PH: Primate ecology and social organization. J Zool, London. 1977, 183: 1-39.View ArticleGoogle Scholar
- Cowling RM, MacDonald IAW, Simmons MT: The Cape Peninsula, South Africa: physiographical, biological and historical background to an extraordinary hot-spot of biodiversity. Biodiver Conserv. 1996, 5: 527-550. 10.1007/BF00137608.View ArticleGoogle Scholar
- Higham JP, Warren Y, Adanu J, Umaru BN, MacLarnon AM, Sommer V, Ross C: Living on the edge: life-history of olive baboons at Gashaka-Gumti National Park, Nigeria. Am J Primatol. 2009, 71: 293-304. 10.1002/ajp.20651.View ArticlePubMedGoogle Scholar
- Alston KP, Richardson DM: The roles of habitat features, disturbance, and distance from putative source populations in structuring alien plant invasions at the urban/wildland interface on the Cape Peninsula, South Africa. Biol Conserv. 2006, 132: 183-198. 10.1016/j.biocon.2006.03.023.View ArticleGoogle Scholar
- Honig MA, Cowling RM, Richardson DM: The invasive potential of Australian Banksias in South African fynbos: A comparison of the reproductive potential of Banksia ericifolia and Leucadendron laureolum. Aust J Ecol. 1992, 17: 305-14. 10.1111/j.1442-9993.1992.tb00812.x.View ArticleGoogle Scholar
- Hoffman TS, O'Riain MJ: The spatial ecology of chacma baboons (Papio ursinus) in a human-modified environment. Int J Primatol. 2010, 32: 308-328.View ArticleGoogle Scholar
- Hill CM: A conflict of interest between people and baboons: crop raiding in Uganda. Int J Primatol. 2000, 21: 299-315. 10.1023/A:1005481605637.View ArticleGoogle Scholar
- Tweheyoa M, Hill CM, Obua J: Patterns of crop raiding by primates around the Budongo Forest Reserve, Uganda. Wildlife Biol. 2005, 11: 237-247. 10.2981/0909-6396(2005)11[237:POCRBP]2.0.CO;2.View ArticleGoogle Scholar
- Naughton-Treves L: Predicting Patterns of Crop Damage by Wildlife around Kibale National Park, Uganda. Conserv Biol. 1998, 12: 156-168. 10.1046/j.1523-1739.1998.96346.x.View ArticleGoogle Scholar
- Versfeld DB, van Wilgen BW: Impacts of woody aliens on ecosystem properties. The Ecology and Control of Biological Invasions in South Africa. Edited by: Mac-donald IAW, Kruger FJ, Ferrar AA. 1986, Cape Town: Oxford University Press, 239-246.Google Scholar
- DeVore I, Hall KRL: Baboon Ecology. Primate Behavior: Field Studies of Monkeys and Apes. Edited by: DeVore I. 1965, New York: Holt, Rinehart and Winston, 20-52.Google Scholar
- Anderson JR: Ethology and ecology of sleep in monkeys and apes. Adv Stud Behav. 1984, 14: 165-229.View ArticleGoogle Scholar
- Richardson DM, van Wilgen BW, Higgins SI, Trinder-Smith TH, Cowling RM, McKell DH: Current and future threats to plant biodiversity on the Cape Peninsula, South Africa. Biodivers Conserv. 1996, 5: 607-647. 10.1007/BF00137612.View ArticleGoogle Scholar
- Norton DA: Editorial: Conservation Biology and Private Land: Shifting the Focus. Conserv Biol. 2000, 14: 1221-1223. 10.1046/j.1523-1739.2000.01451.x.View ArticleGoogle Scholar
- Gunningham N, Young MD: Toward optimal environmental policy: The case of biodiversity conservation. Ecol Law Quart. 1997, 24: 243-298.Google Scholar
- Elmendorf CS: Ideas, incentives, gifts, and governance: Toward conservation stewardship of private land, in cultural and psychological perspective. U Illinois Law Rev. 2003, 423-505.Google Scholar
- Lochner P, Weaver A, Gelderbloma C, Peart R, Sandwith T, Fowkes S: Aligning the diverse: the development of a biodiversity conservation strategy for the Cape Floristic Region. Biol Conserv. 2003, 112: 29-43. 10.1016/S0006-3207(02)00394-4.View ArticleGoogle Scholar
- Hannah L, RakotosaMimanana B, Ganzhorn J, Mittermeier RA, Olivier S, Iyer L, Rajaobelina S, Hough J, Andriamialisoa F, Bowles I, Tilkin G: Participatory planning, scientific priorities, and landscape conservation in Madagascar. Environ Conserv. 1998, 25: 30-36. 10.1017/S0376892998000071.View ArticleGoogle Scholar
- Younge A, Fowkes S: The Cape Action Plan for the Environment: overview of an ecoregional planning process. Biol Conserv. 112: 15-28.Google Scholar
- Lombard AT, Cowling RM, Vlok JHJ, Fabricius C: Designing conservation corridors in production landscapes: assessment methods, implementation issues, and lessons learned. Ecology and Society. 2010, 15: 7-Google Scholar
- Manly BFJ, McDonald LL, Thomas DL: Resource selection by animals. 1993, New York: Chapman and HallView ArticleGoogle Scholar
- Dickson BG, Beier P: Home range and habitat selection by adult Cougars in southern California. J Wildlife Manage. 2002, 66: 1235-1245. 10.2307/3802956.View ArticleGoogle Scholar
- Henzi SP, Byrne RW, Whiten A: Patterns of movement by baboons in the Drakensberg Mountains: Primary responses to the environment. Int J Primatol. 1992, 13: 601-629. 10.1007/BF02551256.View ArticleGoogle Scholar
- Potts JM, Elith J: Comparing species abundance models. Ecol Model. 2006, 199: 153-163. 10.1016/j.ecolmodel.2006.05.025.View ArticleGoogle Scholar
- Zuur AF, Ieno EN, Walker NJ, Savaliev AA, Smith GM: Mixed Effects Models and Extensions in Ecology with R. 2009, New York: Springer-VerlagView ArticleGoogle Scholar
- Graham MH: Confronting multicollinearity in ecological multiple regression. Ecology. 2003, 84: 2809-2815. 10.1890/02-3114.View ArticleGoogle Scholar
- Graham MH: Factors determining the upper limit of giant kelp Macrocystis pyrifera Agardh along the Monterey Peninsula, central California, USA. J Exp Biol. 1997, 218: 127-149. 10.1016/S0022-0981(97)00072-5.View ArticleGoogle Scholar
- Hilbe JM: Negative Binomial Regression. 2008, Cambridge: Cambridge University PressGoogle Scholar
- Cragg JG: Some statistical models for limited dependent variables with application to the demand for durable goods. Econometrica. 1971, 39: 829-844. 10.2307/1909582.View ArticleGoogle Scholar
- Maunder MN, Punt AE: Standardizing catch and effort data: a review of recent approaches. Fish Res. 2004, 70: 141-159. 10.1016/j.fishres.2004.08.002.View ArticleGoogle Scholar
- Punt AE, Walker TI, Taylor BL, Pribac F: Standardization of catch and effort data in a spatially-structured shark fishery. Fish Res. 2000, 45: 129-145. 10.1016/S0165-7836(99)00106-X.View ArticleGoogle Scholar
- Venables WN, Smith DN: An introduction to R. Notes on R: a programming environment for data analysis and graphics version 2.11.1. 2010, [http://www.r-project.org/]Google Scholar
- Jackman S, Tahk A, Zeileis A, Maimone C, Fearon J: pscl: Political Science Computational Laboratory, Stanford University. R package version 1.03.5. 2010, [http://pscl.stanford.edu/]Google Scholar
- Heinänen S, Rönkä M, von Numers M: Modelling the occurrence and abundance of a colonial species, the arctic tern Sterna paradisaea in the archipelago of SW Finland. Ecography. 2008, 31: 601-611. 10.1111/j.0906-7590.2008.05410.x.View ArticleGoogle Scholar
- Pearce J, Ferrier S: Evaluating the predictive performance of habitat models developed using logistic regression. Ecol Model. 2000, 133: 225-245. 10.1016/S0304-3800(00)00322-7.View ArticleGoogle Scholar
- Cliff AD, Ord JK: Spatial processes: Models and applications. 1981, London: PionGoogle Scholar
- Legendre P: Spatial autocorrelation: Trouble or new paradigm?. Ecology. 1993, 74: 1659-1673. 10.2307/1939924.View ArticleGoogle Scholar
- Boyce MS, Vernier PR, Nielsen SE, Schmiegelow FKA: Evaluating resource selection functions. Ecol Model. 2002, 157: 281-300. 10.1016/S0304-3800(02)00200-4.View ArticleGoogle Scholar
- Tyre AJ, Possingham HP, Lindenmayer DB: Matching observed pattern with ecological process: can territory occupancy provide information about life history parameters?. Ecol Appl. 2001, 11: 1722-1738. 10.1890/1051-0761(2001)011[1722:IPFPCT]2.0.CO;2.View ArticleGoogle Scholar
- Campbell BM: A classification of the mountain vegetation of the fynbos biome. Mem Bot Surv S Afr. 1985, 50: 7-11.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.