Management factors affecting densities of common grassland birds of high elevation grasslands of eastern South Africa: Ingula as a case study

Globally, the grassland biome and associated species are threatened and require management interventions to prevent further habitat loss and species extinction (With et al. 2008; Ceballos et al. 2010; Richardson et al. 2014). In southern Africa, grasslands are also one of the most threatened biomes, with 23% under cultivation, 60% irreversibly transformed, 2% formally protected and most of the remaining natural areas used as rangeland for livestock (Fairbanks et al. 2000; Reyers et al. 2001; O’Connor and Kuyler 2009). An assessment of conservation priorities in the Grassland Biome in southern Africa identified 36.7% of the biome as being important for biodiversity conservation (Egoh et al. 2011). In particular, South Africa’s moist grasslands harbor globally significant biodiversity, supply essential ecosystem services, support crop and livestock, forestry and settlement, yet are poorly conserved (O’Connor and Kuyler 2009). This area also coincides with high avifaunal species richness, diversity and endemism (Maphisa et al. 2016). As a result, a number of Important Bird and Biodiversity Areas (IBAs) have been proclaimed in the area (e.g. Barnes 1998). Despite this, only a few studies (e.g. Mentis and Bigalke 1981; Jansen et al. 1999) have examined abundance or densities of bird species in relation to habitat features that can be managed. Regionally, fire is an important conservation management tool used both on private and state lands to make habitats suitable for birds (Parr and Chown 2003). In addition to heavy livestock grazing and other agricultural activities, this area is a target for energy development projects to meet growing human demand for electricity and water for consumption (Maphisa et al. 2016).

Basing conservation decisions on reliable density estimates is critical to halt current species habitat loss and species extinction (Wintle et al. 2005; Post van der Burg et al. 2011). However, estimating population densities of animals is notoriously difficult because some animals remain undetected during surveys (Diefenbach and Brauning 2003). Many grassland bird species look alike and are therefore difficult to identify in the field. Furthermore, the detection probability varies across methods and different observers, changes with habitats and weather conditions (Petitot et al. 2014). Hierarchical distance-sampling methods account for the detection process and are therefore commonly used to study animal density (Marques et al. 2007; Oedekoven et al. 2013). Improvements in distance sampling allow for the incorporation of environmental and habitat variables to fully explore variation in density of species across sites (Royle et al. 2004; Thomas et al. 2010).

Hierarchical distance sampling models consist of two components (Fiske and Chandler 2011; Sillett et al. 2012): the first component models the observation process (detection probability as a function of distance from the transect line) and the second models the biological process (density of the focal bird species). The observation model involves the choice of a distance function whose parameters will be estimated from the data. This approach has recently been used to relate density of birds to habitat variables within Fynbos Vegetation in South Africa (Lee et al. 2015) but has not been used within high-elevation grasslands of southeastern Africa. We present density estimates based on hierarchical distance sampling models for common bird species in relation to vegetation attributes at Ingula, a high-elevation grassland in eastern South Africa which will soon be declared a nature reserve (Maphisa 2012).

Few studies examine how habitat and management actions affect bird densities at high-elevation grasslands of eastern South Africa. The few studies that exist, examine the effect of burning and grazing on a single species (Parr and Chown 2003) and do not account for the observation process (e.g. Jansen et al. 1999; Little et al. 2013). Burning is the main management tool in agricultural lands and conservation areas in this region (Parr and Chown 2003). Some studies indicate that inappropriate use of fire and grazing leads to degradation of grasslands in this region (Maphisa et al. 2009, 2016; Little et al. 2015a).

We relate grass height and cover in summer to eight small grassland birds that we found to be the most common during preliminary data analysis. Summers coincide with breeding for most birds and is the period with the highest species richness in this area (Maphisa et al. 2016). We used the following eight species: Cape Longclaw (Macronyx capensis), Wing-snapping Cisticola (Cisticola ayresii), African Pipit (Anthus cinnamomeus), African Quail Finch (Ortygospiza atricollis), Yellow-breasted Pipit (Anthus chloris), Red-capped Lark (Calandrella cinerea), Zitting Cisticola (Cisticola juncidis) and Common Quail (Coturnix coturnix). Of these, only the Yellow-breasted Pipit is threatened nationally (Taylor et al. 2015). Conservation managers often rely on surrogates to estimate biodiversity to manage habitat for other species whose density cannot be estimated by distance sampling because they are rare (Nally and Fleishman 2004; Mattsson and Marshall 2009; Lewandowski 2010). Because the suite of species we chose to study represents a range of different habitat requirements in terms of their responses to grass height and cover, we suggest that they can be used as indicators to monitor habitat suitability for rare threatened species which are of conservation concern.

The number of detections (shown in brackets) per species differed amongst species. We found the following species: Cape Longclaw (184), African Pipit (109), Wing-snapping Cisticola (106), Yellow-breasted Pipit (70), Zitting Cisticola (65), Common Quail (60), African Quailfinch (36) and Red-capped Lark (33). For most species, the best detection function depended on habitat variables, implying that grass height and cover affected our ability to detect individual birds (Table 1). In particular, the model that combined grass height plus cover best explained the detection of four out of eight species (Table 1). Correcting for these effects, we then examined how grazing, burning, grass height and grass cover affected the density of the eight species differently (Table 2; Figs. 1, 2, 3, 4, 5, 6).

	Detection models	CLC	WSC	AP	AQF	YBP	RCL	ZC	CQ
1.	dens(.)p(.)	4.08	7.35	6.00	14.41	9.71	13.36	10.26	1.49
2.	dens(.)p(Burning)	0.00	5.47	7.84	16.37	6.50	2.71	6.68	0.97
3.	dens(.)p(Grazing)	6.41	9.13	9.24	3.81	10.22	0.00	4.80	4.26
4.	dens(.)p(cover)	4.15	8.63	7.58	9.94	10.00	7.34	12.26	2.46
5.	dens(.)p(avh)	6.08	0.00	1.70	16.19	5.13	13.45	4.24	0.00
6.	dens(.)p(avh + cover)	5.37	1.58	0.00	0.00	0.00	6.08	0.00	2.64

The table shows delta AIC values, i.e., the difference in AIC between a particular model and the best. The fitted models are hierarchical distance sampling models and full model selection details are given in Additional file 1: Appendix S1. The best detection model for each species based on lowest delta AIC is highlighted in italics. Density dens(.) was held constant while considering the effects of transect habitat on the detection function p(.). The habitat variables were grass cover (‘cover’), grass height (‘avh’) and whether the area had been grazed (‘Grazing’) or burned (‘Burning’) prior to the survey. Species are sorted from the most common (left) to least common (right); Cape Longclaw (CLC), Wing-snapping Cisticola (WSC), African Pipit (AP), African Quailfinch (AQF), Yellow-breasted Pipit (YBP), Red-capped Lark (RCL), Zitting Cisticola (ZC) and Common Quail (CQ)

	Density models	CLC	WSC	AP	AQF	YBP	RCL	ZC	CQ
1.	dens(.)	2.72	18.36	26.61	15.93	24.24	8.81	22.78	8.03
2.	dens(Burning)	2.20	0.00	28.40	11.87	0.00	5.86	24.08	7.12
3.	dens(Grazing)	0.34	19.34	30.46	13.56	24.26	6.38	18.26	8.29
4.	dens(cover)	0.11	20.13	24.86	6.03	26.23	8.25	0.00	8.09
5.	dens(avh)	0.00	20.27	11.93	17.83	16.84	5.58	43.39	7.05
6.	dens(avh + cover)	0.55	22.09	0.00	0.00	16.55	0.00	12.60	8.29
7.	dens(s(cover))	1.71	8.24	30.53	4.82	23.17	21.58	24.94	0.00
8.	dens(s(avh))	0.73	3.59	21.01	35.61	22.31	20.08	3.93	6.95

For each species the detection function was modeled as in the best model in Table 1. The relationship between density and grass height and cover was modeled as linear (on the log scale) relationships (avh and cover) and using regression splines with two knots [s(avh) and s(cover)]. See notes to Table 1 for more detail

For Zitting Cisticola, Cape Longclaw, Red-capped Lark and African Quailfinch, the model allowing for different densities between transects that differed in how they were grazed (Model 3, Table 2) was slightly better supported by the data than the constant model (Model 1, Table 2), even though it was not the best supported model for any species. Zitting Cisticolas tended to occur at higher densities on heavily grazed transects than on ungrazed ones (Fig. 1), while the Cape Longclaw occurred at its lowest density on heavily grazed transects (Fig. 1). Red-capped larks occurred in their highest densities on lightly grazed transects where African Quailfinches had the lowest densities (Fig. 1).

For the Wing-snapping Cisticola and Yellow-Breasted Pipit, the model allowing densities to vary between burned and unburned transects (Model 2, Table 2) was the best model, while for Cape Longclaw, African Quailfinch, Red-capped Lark and Common Quail, this model was better than the constant model (Model 1, Table 2). These species were more abundant on burned transects, except for Red-capped Larks, which were more abundant on unburned transects (Fig. 2).

For all species, one of the models involving grass height or cover (Models 4–8, Table 2) was better supported than the constant model (Model 1, Table 2). Model 6, allowing density to vary with both grass height and cover was clearly the best model for the African Pipit, African Quailfinch and Red-capped Lark. For Cape Longclaws Model 5, including only grass height, was the best even though several of the other models were close competitors for this species. For the Zitting Cisticola, the best model was Model 4 assuming a linear effect of grass cover on density. The Common Quail was the only species for which density was best described by a model that included a nonlinear relationship, Model 8, with a regression spline for cover. Overall, all eight species responded differently to the effect of increasing grass height and grass cover (Figs. 3, 4, 5, 6). Of these species, Cape Longclaws apparently benefitted from little grazing with their density increasing with increases in grass height and grass cover (Fig. 3). The Wing-snapping Cisticola, African Quailfinch, Yellow-breasted Pipit, Red-capped Lark and Common Quail densities decreased with increasing grass height and grass cover (Figs. 3, 4, 5, 6). African Pipit densities increased with increasing grass cover (Fig. 4). On the one hand the densities of Zitting Cisticolas increased with an increase in grass cover and Common Quail densities clearly peaked at intermediate levels of grass cover (Fig. 6).

Our results support previous findings that sometimes grassland birds co-occur with differing habitat requirements (Maphisa et al. 2009, 2016). Fire and grazing affect aboveground biomass differently (Richardson et al. 2014). Depending on the intensity of fire and grazing, heterogeneity of micro niches is created across the landscape that affects bird species densities differently. The impact of fire and grazing on landscape and habitat attributes has been well documented (Jacobs et al. 2012). Our study is the first in the region to examine how vegetation variables affect the density of birds in a hierarchical way that incorporates detection. Grazing and burning affect species differently depending on their functional groups or the preferences of individual species (Richardson et al. 2014; Mero et al. 2015; Maphisa et al. 2016). Under managed burning and grazing (e.g. Klug et al. 2010; Richardson et al. 2014; Mero et al. 2015) we expected species that prefer short grass to be more abundant on burned transects (e.g. Apfelbaum and Haney 1981) and thereafter with little grazing as it was at Ingula the density was expected to decrease with an increase in grass height. Most of the species we studied exhibited this with only two exceptions: both Cape Longclaws and Zitting Cisticolas increased in density with an increase in grass height (Figs. 3, 6). Because burning and grazing were not managed at our study site, we therefore expected grass height and cover to be better predictors of species density because of the ambiguity of visual classification of grazing in the field. For example, higher densities of Red-capped Larks and Common Quails (Fig. 2) on burned and unburned transects contradict our knowledge of habitat preferences for these two species. The same is true for the Common Quail which occurs in similar densities across the three categories of grazing (Fig. 1).

For three of the eight species we examined, we found that a combination of grass height and cover is more important than just grass height alone or grass cover alone. This agrees with previous findings that fire and burning complement each other to create habitats for birds (e.g. Fuhlendorf et al. 2009). Our study demonstrates that grass height and cover or each one alone, depending on the habitat requirements of individual species, determines habitat use by birds (Winter et al. 2005; Devereux et al. 2006; Donald et al. 2010). Regional studies that link species abundance to habitat do not take into account detection (e.g. Jansen et al. 1999; Little et al. 2013, 2015b). The same is true about the studies from north eastern African high-elevation grasslands (Mamo et al. 2016).

Grasslands habitat in our study area is mostly in its natural state other than for degradation from annual burning and heavy grazing. A similar study from northern Uruguay found that vegetation structure has variable effects on the density of grassland birds too (Azpiroz and Blake 2016). However, this study did not account for variability in detection because of insufficient sample size. To deal with problems of insufficient data some studies estimate relative densities (Post van der Burg et al. 2011). The study by Azpiroz and Blake (2016) mostly resembles our study because it includes natural grasslands, managed with cattle, game or sheep. Overall, Azpiroz and Blake (2016) found that bird numbers were more correlated with grass cover than with grass height in contrast to observations of Henslow’s Sparrow (Ammodramus henslowii) at modified Missouri grasslands (e.g. Jacobs et al. 2012). But this is because the comparison of Azpiroz and Blake (2016) was made with croplands where vegetation height tends to be higher but sparse. Also the definition of grass cover for these two studies includes forbs or shrubs. In our study we excluded forbs in our analysis because of sparse forb data.

As expected in our study area, the African Quailfinch indicated steeper density declines with increasing grass cover than with grass height. This species is associated more with croplands than with natural grasslands (Kosicki et al. 2014; Maphisa 2015). Cape Longclaws, the most common species in a group, benefitted from little grazing in our study area with their density rising with increasing grass height and grass cover. Our results contradict Mamo et al. (2014) where a similar but endemic Abyssinian Longclaw (Macronyx flavicollis) occurred with higher densities in more open habitats.

Our results for the African Pipit (Fig. 4), a grassland generalist, agree with our knowledge of how this species reacts to changes in grass height and cover. The rising density of Zitting Cisticolas with increasing grass height agrees with the observation that this species requires high grass to breed (Maphisa 2015). The lowest densities of the threatened Yellow-breasted Pipit on heavily grazed transects (Fig. 1) may explain why this species is rare, given the widespread heavy grazing in the region (Maphisa et al. 2016). On the one hand the same species occurs in much higher densities on transects that were burned (Fig. 2) indicating a need for visual inspection of habitat (e.g. Maphisa et al. 2016). This makes sense given that in the case of our study area there was little grazing. Thick tall grass is worse than heavily grazed grass because of the possible effects of visual obstruction (e.g. Devereux et al. 2006). This is important for ground nesting birds that breed and feed on the ground.

The Common Quail is a regional local migrant with birds easily identified by call and sight early in the breeding season, when birds establish their territories. Birds are easy to see in the field at this time given widespread burning. But it will be harder to see or hear later in the season even when they are present because of increase in grass height and grass cover (Fig. 2) and that bird calls are less frequent. Alternatively, our study site could be a stopover with birds moving to areas that are intensively ploughed. Agricultural fields are a more suitable habitat for this species (Kosicki et al. 2014).

For most species in our study, grazing had no strong effect on density. However, this could be because there was little grazing on site to compare amongst categories of grazing. In contrast, a study with replicated control plots from mountainous grasslands in Norway (Loe et al. 2007) managed with sheep, found that the density of birds was much higher on farms with high sheep densities. But high bird densities on heavily grazed grasslands do not necessarily indicate habitat suitability because this might simply indicate food availability (Loe et al. 2007) but may not support breeding.

Land set aside for conservation in our region is increasing in the form of private reserves and state owned lands. On these lands, we recommend that researchers implement the recommendations by Parr and Chown (2003) regarding experimental research on the effect of fire intensity, timing and variability across the landscape on fauna. A further consideration that needs to be taken into account is the number of replicates compared between sites that are managed differently (e.g. Parr and Chown 2003).

Most species of conservation concern are found on privately owned farms where application of fire and grazing is mostly variable (Howell et al. 2009; Hanley et al. 2012; Jacobs et al. 2012; Ciuzio et al. 2013). In these differently managed lands, we suggest that measurements of vegetation variables be limited to grass height and cover because these are the key habitat variables that are mostly associated with habitat use by birds (Fisher and Davis 2010; Hovick et al. 2014) and doing so also allows survey efforts to cover a large area in a short space of time. Finally we suggest that investigation of factors that affect birds density or suitability be based on surveys from inside the farms, rather than counting along roads (e.g. Wellicome et al. 2014). Habitat on road verges often does not resemble habitat within farms and birds nesting along road verges may be more exposed to predation (Wellicome et al. 2014) and may suffer more from disturbances.

Our study supports previous suggestions that worldwide fire and grazing complement each other to make habitats suitable for birds (Fuhlendorf et al. 2009). However, the most important proximate habitat variables in grasslands are grass height and cover (Fisher and Davis 2010). These two variables can be managed with fire and grazing to improve habitat suitability for grassland bird communities (Fuhlendorf et al. 2009; Maphisa et al. 2009).

Grass height and cover were important predictors of density for all eight common grassland species we examined. However, the species reacted to variation in grass height and cover in different ways. Our results therefore suggest that management should aim to provide a mosaic of patches with different grass height and cover to promote habitat suitability for grassland birds in general in the grasslands of eastern South Africa. Combined use of fire and grazing could be used to create such a mosaic, which is likely to benefit rare species, in addition to the common ones we studied.