Introduction
Human-induced climate change is currently viewed as one of the principal threats to biodiversity conservation (Stott et al. Reference Stott, Tett, Jones, Allen, Mitchell and Jenkins2000, Thomas et al. Reference Thomas, Cameron, Green, Bakkenes, Beaumont, Collingham, Erasmus, de Siqeira, Grainger, Hannah, Hughes, Huntley, van Jaarsveld, Midgley, Miles, Ortega-Huerta, Peterson, Phillips and Williams2004). Mean global temperature has increased by 0.6°C since 1950 and is forecast to rise by around 3°C at the end of the current century (IPCC 2001). This warming at the global scale is predicted to result in significant changes (both increase and decrease) in precipitation and temperature at the local scale affecting habitats (Hannah et al. Reference Hannah, Midgley and Millar2002), with impacts on the distribution and abundance of species (Erasmus et al. Reference Erasmus, Van Jaarsveld, Chown, Kshatriya and Wessels2002, Keith et al. Reference Keith, Akcakaya, Thuiller, Midgley, Pearson, Phillips, Regan, Araujo and Rebelo2008).
Ecological responses of organisms to climate change have now been documented throughout the world (see Parmesan Reference Parmesan2006 and references therein), with at least one species having become extinct as a direct result of climate change (Pounds et al. Reference Pounds, Fogden and Campbell1999). Impacts on birds have also been predicted and documented. For example, mean range shifts of more than 500 km to the north-east have been predicted for European nesting birds, with concomitant decrease in distributional range (Huntley et al. Reference Huntley, Collingham, Green, Hilton, Rahbek and Willis2006, Reference Huntley, Green, Collingham and Willis2007, Reference Huntley, Collingham, Willis and Green2008). However, only a few studies have investigated the effect of climate change on the conservation status of birds in Africa (Erasmus et al. Reference Erasmus, Van Jaarsveld, Chown, Kshatriya and Wessels2002, Simmons et al. Reference Simmons, Barnard, Dean, Midgley, Thuiller and Hughes2004, Coetzee et al. Reference Coetzee, Robertson, Erasmus, van Rensburg and Thuiller2009, Holt et al. Reference Holt2009). These studies typically predict a decrease in distributional range (Thomas et al. Reference Thomas, Cameron, Green, Bakkenes, Beaumont, Collingham, Erasmus, de Siqeira, Grainger, Hannah, Hughes, Huntley, van Jaarsveld, Midgley, Miles, Ortega-Huerta, Peterson, Phillips and Williams2004) and/or a decline in population size (Wichmann et al. Reference Wichmann, Jeltsch, Dean, Moloney and Wissel2003) over the next century, with some possible extinctions (Erasmus et al. Reference Erasmus, Van Jaarsveld, Chown, Kshatriya and Wessels2002). The severity of these impacts has been considered to be related to ecological and life-history traits of the species; long-lived, territorial species with low reproductive rates, restricted global distributions and reliance on vulnerable habitats are especially at risk (Simmons et al. Reference Simmons, Barnard, Dean, Midgley, Thuiller and Hughes2004). Birds of prey characteristically exhibit such features and may be prime candidates for studies on climate change (Wichmann et al. Reference Wichmann, Jeltsch, Dean, Moloney and Wissel2003).
The Sokoke Scops Owl Otus ireneae is a small, range-restricted species endemic to coastal Kenya and north-eastern Tanzania (Evans Reference Evans1997, Virani et al. Reference Virani, Njoroge and Gordon2010) that is currently listed as globally “Endangered” (IUCN 2012). The listing was justified by the very small, severely fragmented range with declining habitat quality and the removal of nesting trees (IUCN 2012). Its entire global distribution consists of three small, isolated populations: 1) the adjoining Manga and Kwamgumi forest reserves covering 97 km2 in the lowlands of the East Usambara mountains, Tanzania (Evans et al. Reference Evans, Watson, Hipkiss, Kiure, Timmins and Perkin1994, Virani Reference Virani1995, Evans Reference Evans1997); 2) the Arabuko-Sokoke Forest Reserve (ASF) containing 220 km2 of suitable habitat, coastal Kenya (Ripley Reference Ripley1966, Virani et al. Reference Virani, Njoroge and Gordon2010); and 3) the Dakatcha woodlands 30 km north of ASF covering an unknown area, Kenya (Virani et al. Reference Virani, Njoroge and Gordon2010). In total, therefore, its global range is unlikely to exceed 500 km2, and may well be less than 400 km2.
Two additional factors exacerbate the situation further; the species’s restricted habitat requirements and the loss and fragmentation of this habitat. Despite the availability of three different forest habitats at ASF, the Sokoke Scops Owl occurs only in Brachylaena/Cynometra-dominated forest (Virani et al. Reference Virani, Njoroge and Gordon2010), avoiding Brachystegia and mixed forests. Forest cover outside ASF has disappeared over the past half-century in this region (Kelsey and Langdon Reference Kelsey and Langdon1984), and illegal harvesting of timber is affecting the forest quality within ASF (Thompson et al. Reference Thompson, Eshiamwata, Githiru, Matiku and Ayiemba2007). The habitat of the Sokoke Scops Owl has not been quantified in the East Usambaras or at Dakatcha, but in both areas this species inhabits some form of lowland forest (Evans Reference Evans1997). Furthermore, in logged parts of the Usambaras, Sokoke Scops Owl densities were half that of unlogged forests, emphasising the dependence of this species on undisturbed lowland forest. Not surprisingly, this species has been flagged as being at high risk of extinction in the medium term (Collar et al. Reference Collar, Crosby and Stattersfield1994). Disturbingly, the population of this species has been documented to be in decline at ASF from an estimated 1,025 pairs in 1993 down to 800 pairs in 2008 (Virani et al. Reference Virani, Njoroge and Gordon2010).
In this paper, we investigate the possible effects of climate change on the distribution of the Sokoke Scops Owl. Our primary aim is to assess the likelihood of long-term persistence of this species. Our approach is to use niche modelling to describe the current distribution of the species, and then project its distribution into the future. Our specific objectives are: 1) to model the current distribution of the Sokoke Scops Owl using Maxent (Phillips et al. Reference Phillips, Anderson and Chapire2006), and 2) to model its distribution in 2080 under the IPCC Special Report on Emission Scenarios (SRES) scenario A2. The rationale for our approach is that modelling will allow us to predict the future distribution of this endangered species and therefore develop a conservation management plan for this species in light of the ongoing change in global climate.
Methods
Study sites
The known global distribution of the Sokoke Scops Owl ranges patchily from approximately 3° to 5° S and from 38° to 40° E, making it a narrow-ranged East African endemic. Hence, we designated our study area from 1.5° to 7.5° S and from 38° to 42° E (Figure 1). Much of the natural habitat of the Sokoke Scops Owl (i.e. coastal lowland forest) within this region has been transformed, degraded or disturbed (Thompson et al. Reference Thompson, Eshiamwata, Githiru, Matiku and Ayiemba2007), the existing tracts being fragmented and isolated (Figure 1).
Data collection
The distribution of the Sokoke Scops Owl was based on published (Evans Reference Evans1997, Virani et al. Reference Virani, Njoroge and Gordon2010) and unpublished (C. Jackson unpubl. data) sources, totalling 42 point locations (Figure 1). Additionally, surveys were conducted in apparently suitable habitat in southern coastal Kenya (between Mombasa and the Tanzanian border) where this species has not been recorded before, but where suitable climatic conditions are predicted to occur (see Results below). A total of nine Kaya forests (comprising pristine coastal forest habitat) were surveyed for the presence of Sokoke Scops Owl. These nine forests, in order of increasing size, are: Kaya Waa (20 ha), Diani (80 ha), Gandini (150 ha), Mtswakara (247 ha), Mrima Hill (250 ha), Dzombo Hill (295 ha), Mwache (345 ha), Marenje Hill (1,480 ha), and Shimba Hills (9,500 ha) lying 30–300 m above sea level.
Each of the Kaya forests was visited once in either August-September 2007 or August 2008. Each forest was surveyed over two nights, with an early evening 19h00–22h00 and an early morning session 03h00-06h00. Survey methods followed Virani et al. (Reference Virani, Njoroge and Gordon2010) and involved 10-min owl counts at points separated by 200 m. At the beginning of each point, owl calls were elicited for 30 seconds using pre-recorded Sokoke Scops Owl calls recorded from ASF.
Niche modelling
Niche modelling involves predicting the geographical distribution of a species based on environmental conditions at locations where the species is known to occur (presence only data), and has become a popular and effective tool in conservation biology, ecology, evolution and phylogeography (Townsend and Robins Reference Townsend and Robins2003, Graham et al. Reference Graham, Ron, Santos, Schneider and Moritz2004, Martinez-Meyer et al. Reference Martinez-Meyer, Peterson and Hargrove2004, Ficetola et al. Reference Ficetola, Thuiller and Miaud2007, Raxworthy et al. Reference Raxworthy, Ingram, Rabibisoa and Pearson2007, Lamb et al. Reference Lamb, Ralph, Goodman, Bogdanowicz, Fahr, Gajewska, Bates, Eger, Benda and Taylor2008). A number of different modelling techniques have been developed for such “presence only” datasets, of which Maxent (maximum entropy; Phillips et al. Reference Phillips, Anderson and Chapire2006) has consistently performed well (Elith et al. Reference Elith, Graham and Anderson2006, Phillips and Dudik Reference Phillips and Dudik2008, Wisz et al. Reference Wisz, Hijmans, Peterson, Graham and Guisan2008). Furthermore, it has been suggested that Maxent performs well even when dealing with a small number of occurrence records (Hernandez et al. Reference Hernandez, Graham, Master and Albert2006, Pearson et al. Reference Pearson, Raxworthy, Nakamura and Peterson2007, Wisz et al. Reference Wisz, Hijmans, Peterson, Graham and Guisan2008).
One categorical and eight continuous environmental variables were used as predictors in the Maxent models; these variables are listed in Table 1. The continuous variables included a digital elevation model (SRTM30) or were extracted from WorldClim (Hijmans et al. Reference Hijmans, Cameron, Parra, Jones and Jarvis2005; www.worldclim.org). The only categorical variable included was a soil map extracted from the harmonised world soil database (FAO/IIASA/ISRIC/ISSCAS/JRC 2009). All the environmental variables were resampled to a grid resolution of 30 arc sec (roughly 1 x 1 km) and clipped to eastern Kenya and north-eastern Tanzania using ArcView 3.2 and the extension Garp Datasets. All models were run with the seven BioClim variables, the elevation model and the soil map. The 42 presence records correspond to 35 pixels (1 x 1 km).
The Maxent model was run with 75% training and 25% random test data, and the regularisation multiplier was set to 1 (Lamb et al. Reference Lamb, Ralph, Goodman, Bogdanowicz, Fahr, Gajewska, Bates, Eger, Benda and Taylor2008). All other Maxent settings were left on default. Model performance was evaluated by examining the area under the receiver operating characteristic curve (AUC), and by a jack-knife test which examines the importance of each environmental variable, first by removing one variable at a time and then each variable in isolation (Phillips et al. Reference Phillips, Anderson and Chapire2006).
Maxent provides an output of the relative suitability for a species as a continuous variable. We converted this into a binary variable by applying a threshold, produced by Maxent. Hence, values below the threshold indicate areas that were unsuitable for the species, whereas values above the threshold were those that were suitable. The logistic threshold “Equal Training Sensitivity and Specificity” gave the most realistic current distribution prediction for the Sokoke Scops Owl and we therefore used it for all the models.
We projected the future distribution of the Sokoke Scops Owl by using the HadCM3 general circulation model (Gordon et al. Reference Gordon, Cooper, Senior, Banks, Gregory, Johns, Mitchell and Wood2000). We extracted the same bioclimatic data as was used above to model the predicted current distribution from the WorldClim site (www.worldclim.org) for 2080 under the “A2” scenario. This scenario depicts a heterogeneous future world based on a high human population growth rate coupled with increased energy expenditure and large land-use changes (i.e. higher CO2 emissions). We used the “clamp” feature in Maxent to this future prediction which restricts the effect of projecting onto environmental features not present in the training dataset (Pyron et al. Reference Pyron, Burbrink and Guiher2008). The protected areas were extracted from the 2012 World Database on Protected Areas (WDPA), which is a joint project of UNEP and IUCN (www.wdpa.org).
Results
The potential current distribution of the Sokoke Scops Owl is shown in Figure 2a. This model has a good fit (training AUC = 0.996, test AUC = 0.990), but produces a potential distribution that is only partially realised. Some areas (e.g. to the south of ASF) showing as highly suitable have in fact been transformed and the native forest has been replaced with an agricultural landscape unsuitable for Sokoke Scops Owl. The jack-knife test for the model indicated that soil type and precipitation of the coldest quarter (Bio 19), both with regularised training gains of > 1.0, were the most important variables; they contributed 57.3% and 30.3% of the explanatory power of the model, respectively. The response curves showing the way the seven bioclimatic variables affect the predicted distributions are presented in Figure 3. The variables isothermality (the mean diurnal range divided by the temperature annual range; Bio 03) and precipitation seasonality (Bio 15) are inversely related while mean temperature of the driest month (Bio 09) and annual rainfall (Bio 12) are directly related. In contrast, the other variables tended to increase to a peak and then decline again. The bioclimatic suitability model for the Sokoke Scops Owl, therefore, emphasises high, year-round rainfall with low seasonality.
The surveys of Sokoke Scops Owl in the Kaya forests of southern Kenya failed to record this species in any of the nine forests that were surveyed.
The predicted distribution of the Sokoke Scops Owl in 2080 is shown in Figure 2b. The areas with the highest suitability are shifted southwards and eastwards, away from the current core distribution of the species in Kenya. In the case of Tanzania, the future predicted distribution does not include any suitable climatic conditions for the Sokoke Scops Owl. Concomitant with a dramatic shift, the actual area of suitable climatic conditions is also predicted to shrink by 64% in the case of the pixels with highly suitable conditions (Table 2). A more relevant measure to the future conservation of the Sokoke Scops Owl is the predicted number of pixels with highly suitable conditions that fall within protected areas. The current model indicates 840 pixels (c.840 km2) falling within protected areas, whereas by 2080 this declines to 298 pixels (298 km2), a decline of 65%. It is important to note that the Kwamgumi forest reserve area (and all surrounding suitable areas in the East Usambaras) is lost under the future model.
Discussion
Our model predicts that the distribution of Sokoke Scops Owl will shift southward and eastward by 2080. In addition to the shift in the predicted distribution, there are also changes in the amount and location of the suitable environmental space within the distribution. Areas predicted to have high probability-of-presence for the species, or high environmental suitability, decrease by 2080. Given the already restricted range of the Sokoke Scops Owl, decreases in the amount of area with high environmental suitability is likely to severely increase the extinction risk of the species. Southward and eastward shifts of the whole distribution, or areas of high environmental suitability, may have dire consequences for its long-term persistence. To the south and east of the ASF population of Sokoke Scops Owl are vast areas of highly degraded or entirely transformed landscapes without any suitable habitat for the owls (Thompson et al. Reference Thompson, Eshiamwata, Githiru, Matiku and Ayiemba2007). Even more dramatic is the predicted total loss of suitable conditions in the East Usambara population. How this predicted change in future climate will affect the forests currently persisting at these two sites is unknown. However, due to the extremely restricted distribution of the Sokoke Scops Owl, any reduction in suitable habitat is of high conservation concern. Species with high habitat specificity and narrow ranges have been highlighted as the most vulnerable to climate change (Manne and Pimm Reference Manne and Pimm2001, Thomas et al. Reference Thomas, Cameron, Green, Bakkenes, Beaumont, Collingham, Erasmus, de Siqeira, Grainger, Hannah, Hughes, Huntley, van Jaarsveld, Midgley, Miles, Ortega-Huerta, Peterson, Phillips and Williams2004, Williams et al. Reference Williams, Bolitho and Fox2003). Restricted-range species are more vulnerable to extinction overall, due to stochastic events (McKinney Reference McKinney1997), as recognised by the IUCN Red List criteria (IUCN 2001).
Small distributions are often linked to lower abundance and higher habitat specificity (Brown Reference Brown1984). While we have not directly measured abundance in this study, some predicted patterns may be inferred. The environmental suitability predicted by Maxent has been shown to correlate with species’ abundance: where environmental suitability is low, species abundance is generally also low; where environmental suitability is high, species’ abundance can range from high to low (VanDerWal et al. Reference VanDerWal, Shoo, Johnson and Williams2009). The decrease in suitability within a species’s distribution may indicate that the species is likely to decrease in abundance by 2080.
Many studies have shown that the synergy between climate change and other forms of habitat degradation is driving species towards extinction faster than either factor operating alone (Millenium Ecosystem Assessment 2005). This has been shown for birds in particular at a global scale (Jetz et al. Reference Jetz, Wilcove and Dobson2007). For the Sokoke Scops Owl, forest destruction in combination with shifting climate space may lead to highly limited suitable habitat. In addition, it is shown that current conservation efforts will decrease in efficacy with climate change. Current protected areas may decrease in species richness, as species distributions change with changing climate (Coetzee et al. Reference Coetzee, Robertson, Erasmus, van Rensburg and Thuiller2009, Hole et al. Reference Hole, Willis, Pain, Fishpool, Butchart, Collingham, Rahbek and Huntley2009). Very little of the future distribution of the Sokoke Scops Owl will overlap with current protected areas. This is the most likely future situation, as current global greenhouse gas emissions are already exceeding the most severe of the IPCC future scenarios (Raupach et al. Reference Raupach, Marland, Ciais, Le Quéré, Canadell, Klepper and Field2007), making the A2 scenario used in this study very likely.
The future of the Sokoke Scops Owl looks bleak. With predicted decrease in Sokoke Scops Owl distribution, potential decrease in abundance (corresponding to the predicted decrease in suitability), and the lack of suitable protected areas, this owl is facing a very real possibility of extinction within the next 100 years. In addition, the 65% decline in suitable habitat falling within protected areas and ongoing habitat degradation further emphasises the potential for significant population declines in this species. Given this scenario, and the fact that a 25% decline in population over the past two decades has already been documented (Virani et al. Reference Virani, Njoroge and Gordon2010), we strongly recommend that owl surveys to estimate population size continue at least every two years following the protocol outlined by Virani et al. (Reference Virani, Njoroge and Gordon2010).
Some conservation interventions may be possible for the Sokoke Scops Owl. Our study has identified the potential suitable climatic space available to this species in the future. It is theoretically possible that reforestation or habitat restoration could create suitable areas for its future distribution (Shoo et al. Reference Shoo, Storlie, VanDerWal, Little and Williams2011). Re-evaluation of protected areas, and ideally creating new protected areas, is likely to be crucial to the persistence of many species. In addition to the Sokoke Scops Owl, other threatened or localised species that are likely to benefit from habitat restoration and the proclamation of future protected areas in this region include the following birds: Sokoke Pipit Anthus sokokensis, Clarke’s Weaver Ploceus golandi, and Amani Sunbird Hedydipna pallidigaster; and mammals: Sokoke bushy-tailed mongoose Bdeogale omnivora, Aders Duiker Cephalophus adersi and the easternmost race of the African golden cat Felis aurata (Virani Reference Virani1993).
We suggest that exploratory surveys north of the ASF and the Dakatcha Woodlands, coupled with habitat assessments to identify suitable habitat for the Sokoke Scops Owl should be a priority goal. In addition, a species action plan for the Sokoke Scops Owl urgently needs to be drafted with the possibility of considering translocation of the species as a way of expanding its range (Hoegh-Guldberg et al. Reference Hoegh-Guldberg, Hughes, McIntyre, Lindenmayer, Parmesan, Possingham and Thomas2008). The Kaya forests are the obvious candidate for any such translocation since they fall within the predicted future climatic conditions suitable for this species. The absence of the Sokoke Scops Owl from these Kaya forests is surprising since they form a natural link between the forests in the East Usambaras and the ASF. The possibility of this species having recently gone extinct in these Kaya forests needs to be investigated before any translocations are conducted.
Although the use of niche modelling (in particular Maxent) in conservation has increased steadily in the past few years (Ficetola et al. Reference Ficetola, Thuiller and Miaud2007, Recio and Virgos Reference Recio and Virgos2010, Jimenez-Valverde et al. Reference Jimenez-Valverde, Lira-Noriega, Peterson and Soberon2010), its effectiveness as a tool for predicting future distributions and status of species has been challenged (Jimenez-Valverde et al. Reference Jimenez-Valverde, Lobo and Hortal2008). There are a number of reasons for this including: possible lack of relationship between probability of distribution and abundance (Jimenez-Valverde et al. Reference Jimenez-Valverde, Diniz, de Azevedo and Borges2009); confusion of potential and realised distributions (Jimenez et al. Reference Jimenez-Valverde, Lobo and Hortal2008); inappropriate use of AUC (Lobo et al. Reference Lobo, Jimenez-Valverde and Real2007); biases due to small sample size (Wisz et al. Reference Wisz, Hijmans, Peterson, Graham and Guisan2008); lack of distinction between absence and pseudo-absence data; and obviously the reliance on future climate models that cannot be tested in the present. Hence, the interpretation of our data needs to be carefully considered in light of these observations.
Conclusion
The combination of the current restricted range of the Sokoke Scops Owl, and the future projections of range shifts and loss into the future, suggest that this species has a high risk of extinction. Niche modelling has provided a useful mechanism in which to predict these likely changes in distribution. We strongly recommend close monitoring of this species and suggest the investigation of potential intervention in the form of translocations to other suitable habitat patches currently not known to support populations of the Sokoke Scops Owl.
Acknowledgements
Financial support for this study was provided by The Peregrine Fund through its Africa programme. We thank the Kenya Wildlife Service and the National Museums of Kenya for their support. We are indebted to David Ngala, Gabriel Katana, Simon Musila, Bernard Cheruiyot Soi and Shiv Kapila for assistance in the field.