Next Article in Journal
Genetic and Haplotype Diversity of Manila Clam Ruditapes philippinarum in Different Regions of China Based on Three Molecular Markers
Next Article in Special Issue
Evaluating Different Factors That Affect the Nesting Patterns of European and Algerian Hedgehogs in Urban and Suburban Environments
Previous Article in Journal
Effects of Salinity Stress on Histological Changes, Glucose Metabolism Index and Transcriptomic Profile in Freshwater Shrimp, Macrobrachium nipponense
Previous Article in Special Issue
Are Tourists Facilitators of the Movement of Free-Ranging Dogs?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Feral Donkey Distribution and Ecological Impacts in a Hyper-Arid Region

1
The Royal Commission for AlUla, Riyadh 12512, Saudi Arabia
2
Plant Production Department, College of Food & Agriculture Sciences, King Saud University, Riyadh 11451, Saudi Arabia
3
Department of Botany, Faculty of Science, Port Said University, Port Said 42511, Egypt
*
Author to whom correspondence should be addressed.
Animals 2023, 13(18), 2885; https://doi.org/10.3390/ani13182885
Submission received: 6 August 2023 / Revised: 4 September 2023 / Accepted: 7 September 2023 / Published: 11 September 2023
(This article belongs to the Special Issue Ecological Impact of Feral Animals)

Abstract

:

Simple Summary

We estimated the population density and the population size of feral donkeys in northwest Saudi Arabia. The estimated overall population density was 1.03 (0.19 SE) donkey/km2 and an abundance at ~1135 individuals. The negative impact of feral donkeys on natural resources included overgrazing, habitat fragmentation, and competition for resources with native species. Our study adds evidence to the detrimental impacts of feral donkeys and calls for urgent actions to control the presence of feral donkeys in the region. We recommend humane eradication for controlling the feral donkeys as it would be the most efficient and applicable action that can significantly abate the feral donkeys’ negative impacts.

Abstract

The feral donkey (Equus asinus L.) is an invasive species in Saudi Arabia and can cause severe damage to natural and cultural heritage. Over the last 30 years, feral donkeys have become a serious problem, as their abundance and geographic distribution has increased drastically. The impacts of feral donkeys are not well documented, and information about their abundance and distribution is lacking, certainly in Saudi Arabia, which hampers the implementation of effective management plans. Accordingly, we used the minimum population number approach (MPN) to determine the number of feral donkeys in this part of northwest Saudi Arabia. A total of 1135 feral donkeys were encountered in the region. The area around Khaybar harbors ~25% (n = 338) of the feral donkey population, whereas Tayma and AlGhrameel nature reserves were the least-inhabited sites (almost absent). The average population density of feral donkeys was estimated as 1.03 (0.19 SE) donkey/km2. We documented the negative ecological impact of feral donkeys on natural resources, which constituted overgrazing that resulted in habitat fragmentation and competition for resources with native species. We propose urgent actions to control the presence of feral donkeys in the region and suggest humane eradication as the most efficient and applicable to significantly reduce the negative impacts of feral donkeys.

1. Introduction

For more than 4000 years, donkeys have been domesticated as beasts of burden for packing, transporting goods, and tillage [1]. Donkeys (Equus asinus L.), also called “feral burros” [2], are believed to be evolved from the African wild ass (Equss africanus) and lived and worked with humans for 1000s of years in Arabia [3]. They are an equine species with distinct large heads and long ears [4]. The long ears allow donkeys to hear distant calls of other donkeys and are used for cooling in hot weather [5]. They vary in color from white to black, while dun grey is the most common color pattern [4]. Most often, feral donkeys have dorsal and shoulder strips that are primitive markings [6]. Donkeys have a specialized digestive system that allows the consumption of large amounts of low-quality forage and the subsequent extraction of nutrients [2]. Donkeys’ weight and height range between 110 and 142 kg and between 102 and 142 cm, respectively [6]. These variations are a result of the quantity and quality of the available food [7,8]. Donkeys were used for their strength, agility, and adeptness in navigating terrain that was unsuitable for both horses and camels [3]. Since the 1970s, the establishment of an extensive road network has reduced the need for donkeys, and they were set free into the wild, to become feral and now breed and graze in unmanaged herds [9].
The IUCN Invasive Species Specialist Group (ISSG) classified feral donkeys in Saudi Arabia as an invasive species in 2011 [10]. The FAO Statistics Division estimated a population of 55,000 donkeys in 1961 and 100,000 by 2020 in Saudi Arabia [11]. The population size of feral donkeys presents significant ecological problems with consequential socioeconomic implications. For instance, previous studies documented the negative impact of donkey grazing behavior on the availability of forage for other herbivores [12,13,14]. Further, the presence of feral donkeys has various effects on native wildlife, including competition for food and shelter, habitat destruction, and disease transmission [15]. In arid regions, feral equines can spend up to two-thirds of their day at water sources, using aggressive behavior to prevent native wildlife from using these resources [16,17].
Despite the negative ecological impacts of feral donkeys on ecosystem processes and functions, their associated impacts have not been well documented in the Arabian Desert and particularly in Saudi Arabia. Moreover, the numbers of feral donkeys and their distribution are still unknown and no measures have yet been taken to reduce their impact.
The Royal Commission for AlUla (RCU) has committed to restoring and conserving the natural heritage of the Arabian Desert, specifically in the AlUla region. To this end, RCU has established five nature reserves in the region and implemented several restoration and reintroduction projects within these reserves. However, there was a concern that the presence of feral animals, particularly donkeys, could challenge the success of these projects. Accordingly, it is important for the RCU to determine the current spatial distribution and the relative abundance of the feral donkeys, and to evaluate their associated impacts on the natural resources. Such information can be used to inform conservation managers to implement an effective control plan for feral donkeys. Therefore, this study was aimed at assessing the current distribution of feral donkeys and estimating their population size in the region. To this end, we developed a stratified sampling design to estimate the population size of feral donkeys and model the association between the estimated number and the environmental variables.

2. Materials and Methods

2.1. Study Area

The study area is located in the northwest of Saudi Arabia (Figure 1) and under the authority of the Royal Commission for AlUla (RCU). The study area covers a total area of ~25,000 km2 that enclosed AlUla County and the surrounding nature reserves, AlGhrameel Nature Reserve (GNR), Harrat Uwayrid Biosphere Reserve (HUBR), Sharaan Nature Reserve (SNR), Harrat AlZabin Nature Reserve (HZNR), and Wadi Nakhlah Nature Reserve (WNNR). Three other areas were also covered during the survey, Jabal ElWard (120 km west of AlUla County), Tayma (70 km northeast of AlUla County), and Khaybar (200 km southeast of AlUla County) (Figure 1). All reserves are unfenced, with the exception of SNR, which was not included in the survey. The study area is an inland desert habitat that varies in altitude (200–2000 m.a.s.l.). It is characterized by heterogeneous habitats of lava fields, plateaus, sandstone outcrops, sand dunes, and several permanent springs that support wetland and oasis habitats. The RCU’s nature reserves were designed to represent key desert habitats and to facilitate structural and functional connectivity with other reserves. For instance, King Salman bin Abdulaziz Royal Natural Reserve is connected with RCU’s nature reserves through GNR and Tayma, while Prince Mohammed bin Salman Royal Reserve is connected with RCU’s nature reserves through HUBR.

2.2. Data Collection and Analysis

Initially, we developed a stratified sampling design to estimate the population size of feral donkeys within the study area, based on a distance sampling approach [18]. However, applying this approach was very challenging, due to the donkey’s behavior and uneven use of the study area. Feral donkeys are highly sociable and form large groups of around 15 individuals [19]. When humans are present, they emit loud braying sounds as a warning to others, prompting them to flee and disappear into the surroundings [20]. This behavior poses a challenge to distance sampling assumptions, which rely on the detection of objects at their original location and precise measurements [18]. As a result, estimating the detection probability of this species becomes difficult. Additionally, since this study is a baseline to inform the feral animal control plan, rather than a population estimate for long-term monitoring, we decided to use an alternative simpler approach. Hence, we used the minimum population number approach (MPN) to determine the number of feral donkeys in the study area [14,21]. A recent study showed that MPN can calculate the population number with high accuracy and precision and also facilitates temporal comparison [21]. We divided the study area including Jabal ElWard into 900 grids, 5 km × 5 km each. A total of 109 grids were selected to represent the study area. Since the distribution of feral donkeys was uneven across the landscape, the selection of these grids was based on the result of the interview with the local community and shepherds. Additionally, we considered the distance between the selected grids to minimize spatial autocorrelation and to ensure the independence of the observations. For Khaybar and Tayma, due to their small spatial scales (~50 and 20 km2, respectively), each was sampled by a single 5 km × 5 km grid. We used a daytime survey to collect field data from the grids and to calculate the MPN in the study area. The survey was carried out between February and May 2022. This involved both driving transect (~25 km/h) along dirt roads and walking transect through inaccessible areas (e.g., narrow valleys, outcropped habitats, wetlands, and abandoned farms). For the whole study area, excluding GNR and Tayma, transect design (length and the direction) was unsystematic and determined by the topographic feature of the grids, whereas in GNR and Tayma, with a dominant open habitat such as gravel plain and sand sheets, the grids were surveyed using two parallel transects (2 km apart) with a southeast–northwest orientation. The survey was carried out with a team of four observers. For each sighting of a feral donkey, the location (using Garmin GPSMAP 64s, Taiwan) and the number of animals were recorded. The MPN was calculated as the total number of all recorded animals across the study area, while feral donkey population density was calculated as the number of recorded donkeys/km2.
The landscape covariates associated with each grid were also recorded, which include human presence, vegetation cover, presence of water source (yes/no), and topographic roughness. Human presence was calculated as a perpendicular distance to the nearest settlement or farm. The vegetation cover was estimated visually along the sampled transects on a continuous scale between “0” (no vegetation) and “100” (complete cover). The topographic roughness was used here to represent the surface evenness and was calculated from the elevation data using the “terrain” function implemented in “raster” R package [22].
We used the Generalized Linear Mixed Models (GLMMs) to test whether the MPN of feral donkeys was driven by the associated landscape covariates [23]. GLMMs allow for fitting ecological data that are not normally distributed while accounting for random effects by including the sites as a random factor [24,25]. We fitted GLMMs with all possible combinations including interaction terms using the Poisson distribution using the package glmmTMB [26]. The best-fitted model was selected using Akaike’s Information Criterion (AIC) [27].
The spatial distribution of the feral donkeys was mapped using Inverse Distance Weighted interpolation (IDW), a deterministic spatial interpolation approach [28]. The fundamental principle behind IDW is that the distance between a known point and an unknown point influences the unknown property being estimated (i.e., the impact decreases as the distance increases) [28].

3. Results

A total survey effort of 2225 km covered the 109 grids across the study area. The calculated MPN derived was ~1135 feral donkeys (Table 1). The Khaybar area was the most inhabited site with 338 individuals, whereas Tayma and GNR were the least inhabited sites (almost absent). The average population density was estimated as 1.03 (0.19 SE) feral donkeys/km2. This density was estimated over the study area including Jabal ElWard, Tayma, and Khaybar.
The model that included both the linear and quadratic terms of the presence of human and vegetation cover covariates was the best supported (lowest AIC = 207.1) (Table 2). GLMM revealed an association between the MPN of feral donkeys and both the presence of humans (z = −3.690, p = 0.0001; Figure 2(a1)) and the extent of vegetation cover (z = −4.359, p = 0.0001; Figure 2(a2)). Although the best-fitted model showed no statistically significant impact of topographic roughness, feral donkeys were observed as tending to occupy sites that are dominated by rocky outcrop habitats and avoid open plains. The spatial distribution of the calculated MPN showed a high probability of occurrence for feral donkeys on the western side of the study area compared to the eastern side (Figure 2b).

4. Discussion

This study assesses the distribution and population size of the feral donkey in the AlUla region and shows that at a minimum there are approximately 1135 feral donkeys occurring in the study area. The average population density was estimated as 1.03 (0.19 SE) donkeys/km2, with a group size between 5 and 28 individuals. We acknowledge that our estimate was conservative and certainly underestimated the number of feral donkeys in the study area. Nevertheless, this underestimated figure strongly indicates that feral donkeys are overabundant and present a management issue.
Khaybar area harbors the highest number, ~338, of feral donkeys. This could be explained by the presence of ~350 natural water points and several abandoned farms, which allow survival in the extreme desert conditions. Conversely, the near absence of feral donkeys at GNR and Tayma could be attributed to water scarcity and the degradation of the vegetation cover (pers. obs.). Additionally, both GNR and Tayma are open habitats, mostly sand and gravel plains with no outcrops, which makes donkey survival in these areas difficult due to the direct exposure to harsher conditions. Despite these conditions, grazing is still possible in the area where the shepherds provide food (e.g., Alfalfa) and water for their flocks. Such practices might increase the possibility of feral donkeys occupying both GNR and Tayma in the very near future.
Our model shows a positive association between feral donkey population size and human activity. This finding is consistent with a previous study that found a high population density of donkeys in proximity to human presence in South Africa [14]. This is expected as donkeys are taking advantage of the resources provided by humans, such as fodder for their livestock. This association between feral donkeys and humans could increase their negative impacts, such as road traffic accidents and disease transmission, due to the continuous interaction with livestock. The best-fitted model showed no statistical significance for topography. Although this finding is consistent with a previous study [20], we observed a tendency for donkeys to dwell on site with relatively high topographic variability.
During our fieldwork, we recorded a multitude of adverse impacts for feral donkeys, including ecological and socio-economic impacts. Unlike native herbivore species such as Arabian oryx (Oryx leucoryx), Sand gazelle (Gazella marica), Arabian gazelle (Gazella arabica), and Nubian ibex (Capra nubiana), feral donkeys were observed to uproot native plant species during their grazing. Further, they also negatively affect the native plants through overgrazing, selective bark stripping of woody species, and trampling due to their tendency to form herds and their high population size. These detrimental effects of feral donkeys on floristic composition and diversity have been reported by other studies [9,29,30], and they can result in deterioration of ecosystem processes and functions of arid regions [31,32] and contribute to soil erosion [33].
In the deserts of North America, bark damage induced by feral donkeys resulted in significant deterioration in the recruitment and the reproduction of native trees [34]. In our study, feral donkeys were also observed to damage Acacia trees (Vachellia spp.; keystone species in the study area) by girdling the bark of the trees. This behavior could be a result of drought, lacking some minerals or salts in donkeys’ diets, and they might find moisture and minerals in Acacia trees’ bark. However, girdling was very common across the study area and would increase the possibility of Acacia dieback [34].
Living in large groups allows feral donkey to dominate available resources, and their aggressive behavior can drive native wildlife from an area [35]. For instance, feral donkeys were reported to pollute and foul water resources and prevent other native species from drinking [36]. A recent study showed that the presence of feral donkeys in the study area resulted in a ~50% decline in the number of the endemic Arabian Partridges Alectoris melanocephala [37]. In addition, feral donkeys can also facilitate disease transmission between wildlife and domestic animal [15,38].
Feral donkeys also have a socio-economic impact, as they do become associated with human settlements and are reported as destroying local people’s assets and properties. It is a very common behavior now for donkeys to approach human settlements and overturn garbage bins looking for food, which, in turn, increases the dispersal of litter and plastic. Crossing highways and roads was also documented in this study, which increases the likelihood of vehicle–animal collisions in the region. Although there is no authentic documentation of human–donkey accidents in the Arabian peninsula and particularly in Saudi Arabia, it is documented in other countries that feral donkeys are responsible for almost 300 human deaths/year [39].

5. Conclusions

Our study adds evidence to the detrimental impacts of feral donkeys and calls for urgent actions to control the presence of feral donkeys in the region. These actions include castration, adoption/removal, and eradication [39]. However, castration and removal might not be efficient in controlling the donkeys because castrated individuals can still roam large area and continue to destroy habitats for at least 10 years before they die. Meanwhile, the removal of animals would expand the spatial scale of the donkey, without solving the issue. Humane eradication would be the most efficient and applicable action that can significantly abate the feral donkey’s negative impacts. Humane lethal population control, eradication, is a very common practice and is applied in several parts of the world (Australia, Africa, and America) [40,41], and it was prioritized to control invasive species without conflicting animal welfare standards [42]. Applying a non-poisonous eradication program can have a positive secondary effect on other native species, where it increases food availability for many species, particularly regionally endangered species, such as the Striped hyena (Hyaena hyaena) and the Griffon vulture (Gyps fulvus) [43,44].

Author Contributions

Conceptualization, A.S., S.B., A.A. and M.E.-B.; methodology, A.S., M.E.-B., M.D., A.A., N.A.-J. and Y.A.; software, A.S.; validation, A.S., A.M.A., M.E.-B. and S.B.; formal analysis, A.S.; investigation, A.S., M.E.-B., M.D., A.A., N.A.-J. and Y.A.; resources, A.S. and S.B.; data curation, A.S.; writing—original draft preparation, A.S.; writing—review and editing, A.S., S.B., A.M.A., M.D. and M.E.-B.; visualization, A.S.; supervision, A.S., S.B., A.M.A. and M.E.-B.; project administration, A.S. and S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Primary data not presented in this study are available on request from the corresponding author.

Acknowledgments

The authors wish to acknowledge the RCU ranger team and Mohammed Sulayem for facilitating the fieldwork and providing us with a detailed description of the site’s accessibility. Special thanks for Lourens Du Preez–Van Essen for his comments and suggestions on the revised MS.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Meutchieye, F.; Kwalar, N.N.; Nyock, R.A.F. Donkey Husbandry and Production Systems. In Current Donkey Production and Functionality: Relationships with Humans. Book 1; González, F.J.N., Delgado, J.V., Burgos, J.C.V., Eds.; UcoPress: Córdoba, Spain, 2016. [Google Scholar]
  2. Schoenecker, K.A.; King, S.R.B.; Nordquist, M.K.; Nandintsetseg, D.; Cao, Q. Habitat and Diet of Equids. In Wild Equids: Ecology, Management, and Conservation; Ransom, J., Kaczensky, P., Eds.; Johns Hopkins University Press: Baltimore, MD, USA, 2016; pp. 41–57. [Google Scholar]
  3. Turnbull, A.; Wernery, U.; Wernery, R.; Anandh, J.P.; Kinne, J. Survey of Six Infectious Diseases of Feral Donkeys in the United Arab Emirates. Equine Vet. Educ. 2010, 14, 33–38. [Google Scholar] [CrossRef]
  4. Global Invasive Species Database (GISD) Species Profile Equus asinus. Available online: http://www.iucngisd.org/gisd/species.php?sc=639 (accessed on 29 August 2023).
  5. Yılmaz, O.; Boztepe, S.; Ertuğrul, M. The Domesticated Donkey: I—Species Characteristics. Can. J. Appl. Sci. 2012, 4, 339–353. [Google Scholar]
  6. Gichure, M.; Onono, J.; Wahome, R.; Gathura, P. Assessment of Phenotypic Characteristics and Work Suitability for Working Donkeys in the Central Highlands in Kenya. Vet. Med. Int. 2020, 2020, 8816983. [Google Scholar] [CrossRef] [PubMed]
  7. Pearson, R.A.; Ouassat, M. Estimation of the Liveweight and Body Condition of Working Donkeys in Morocco. Vet. Rec. 1996, 138, 229–233. [Google Scholar] [CrossRef]
  8. Cappai, M.G.; Picciau, M.; Pinna, W. An Integrated Approach towards the Nutritional Assessment of the Sardinian Donkey: A Tool for Clinical Nutritionists. Ital. J. Anim. Sci. 2013, 12, e29. [Google Scholar] [CrossRef]
  9. Brinkmann, K.; Patzelt, A.; Dickhoefer, U.; Schlecht, E.; Buerkert, A. Vegetation Patterns and Diversity along an Altitudinal and a Grazing Gradient in the Jabal Al Akhdar Mountain Range of Northern Oman. J. Arid Environ. 2009, 73, 1035–1045. [Google Scholar] [CrossRef]
  10. Mohamed, W.F. Negative Effects of the Feral Donkey Equus asinus (Linnaeus, 1758) in Al-Ula Governorate, Western Saudi Arabia. Kuwait J. Sci. 2021, 49, 1–12. [Google Scholar] [CrossRef]
  11. FAO FAOSTAT. License: CC BY-NC-SA 3.0 IGO. Available online: https://www.fao.org/faostat/en/#data/QCL (accessed on 20 May 2022).
  12. Vetter, S. Investigating the Impacts of Donkeys on a Communal Range in Namaqualand: How Much Does a Donkey “Cost” in Goat Units? Bachelor’s Thesis, University of Cape Town, Cape Town, South Africa, 1996. [Google Scholar]
  13. Samuels, M.I.; Allsopp, N.; Hoffman, T. Mobility Patterns of Livestock Keepers in Semi-Arid Communal Rangelands of Namaqualand, South Africa. Nomad. People 2008, 12, 123–148. [Google Scholar] [CrossRef]
  14. Muller, H.; Bourne, A. Minimum Population Size and Potential Impact of Feral and Semi-Feral Donkeys and Horses in an Arid Rangeland. Afr. Zool. 2018, 53, 139–144. [Google Scholar] [CrossRef]
  15. Burrows, N.D. Feral Animals in the Semi-Arid and Arid Regions of Australia: Origins, Impacts and Control. In On the Ecology of Australia’s Arid Zone; Springer International Publishing: Cham, Switzerland, 2018; pp. 331–373. [Google Scholar]
  16. Perry, N.D.; Morey, P.; Miguel, G.S. Dominance of a Natural Water Source by Feral Horses. Southwest. Nat. 2015, 60, 390–393. [Google Scholar] [CrossRef]
  17. Hall, L.K.; Larsen, R.T.; Westover, M.D.; Day, C.C.; Knight, R.N.; McMillan, B.R. Influence of Exotic Horses on the Use of Water by Communities of Native Wildlife in a Semi-Arid Environment. J. Arid Environ. 2016, 127, 100–105. [Google Scholar] [CrossRef]
  18. Buckland, S.T.; Anderson, D.R.; Burnham, K.P.; Laake, J.L. Distance Sampling. In Encyclopedia of Biostatistics; John Wiley & Sons, Ltd.: Chichester, UK, 2005. [Google Scholar]
  19. Rudman, R. The Social Organisation of Feral Donkeys (Equus asinus) on a Small Caribbean Island (St. John, US Virgin Islands). Appl. Anim. Behav. Sci. 1998, 60, 211–228. [Google Scholar] [CrossRef]
  20. Hamrick, R.G.; Pirgalioglu, T.; Gunduz, S.; Carroll, J.P. Feral Donkey Equus asinus Populations on the Karpaz Peninsula, Cyprus. Eur. J. Wildl. Res. 2005, 51, 108–116. [Google Scholar] [CrossRef]
  21. Champion, G.; Downs, C.T. Status of the Nile Crocodile Population in Pongolapoort Dam after River Impoundment. Afr. Zool. 2017, 52, 55–63. [Google Scholar] [CrossRef]
  22. Hijmans, R.J. Raster: Geographic Data Analysis and Modeling. R Packag. Version 2.5-8. 2016. Available online: https://cran.r-project.org/package=raster (accessed on 20 January 2023).
  23. Bolker, B.M.; Brooks, M.E.; Clark, C.J.; Geange, S.W.; Poulsen, J.R.; Stevens, M.H.H.; White, J.-S.S. Generalized Linear Mixed Models: A Practical Guide for Ecology and Evolution. Trends Ecol. Evol. 2009, 24, 127–135. [Google Scholar] [CrossRef]
  24. Young, M.; Ierodiaconou, D.; Womersley, T. Forests of the Sea: Predictive Habitat Modelling to Assess the Abundance of Canopy Forming Kelp Forests on Temperate Reefs. Remote Sens. Environ. 2015, 170, 178–187. [Google Scholar] [CrossRef]
  25. Kaplan, R.H.; Rosamond, K.M.; Goded, S.; Soultan, A.; Glass, A.; Kim, D.H.; Arcilla, N. Bobolink (Dolichonyx oryzivorus) Declines Follow Bison (Bison bison) Reintroduction on Private Conservation Grasslands. Animals 2021, 11, 2661. [Google Scholar] [CrossRef]
  26. Brooks, M.E.; Kristensen, K.; van Benthem, K.J.; Magnusson, A.; Berg, C.W.; Nielsen, A.; Skaug, H.J.; Maechler, M.; Bolker, B.M. GlmmTMB Balances Speed and Flexibility among Packages for Zero-Inflated Generalized Linear Mixed Modeling. R J. 2017, 9, 378–400. [Google Scholar] [CrossRef]
  27. Maydeu-Olivares, A.; Garcia-Forero, C. Goodness-of-Fit Testing. Int. Encycl. Educ. 2010, 7, 190–196. [Google Scholar]
  28. Lu, G.Y.; Wong, D.W. An Adaptive Inverse-Distance Weighting Spatial Interpolation Technique. Comput. Geosci. 2008, 34, 1044–1055. [Google Scholar] [CrossRef]
  29. MacLaren, C.A. Climate Change Drives Decline of Juniperus seravschanica in Oman. J. Arid Environ. 2016, 128, 91–100. [Google Scholar] [CrossRef]
  30. Shaltout, K.H.; Eid, E.M.; Al-Sodany, Y.M.; Heneidy, S.Z.; Shaltout, S.K.; El-Masry, S.A. Effect of Protection of Mountainous Vegetation against Over-Grazing and Over-Cutting in South Sinai, Egypt. Diversity 2021, 13, 113. [Google Scholar] [CrossRef]
  31. Beschta, R.L.; Donahue, D.L.; DellaSala, D.A.; Rhodes, J.J.; Karr, J.R.; O’Brien, M.H.; Fleischner, T.L.; Deacon Williams, C. Adapting to Climate Change on Western Public Lands: Addressing the Ecological Effects of Domestic, Wild, and Feral Ungulates. Environ. Manag. 2013, 51, 474–491. [Google Scholar] [CrossRef]
  32. Danvir, R.E. Multiple-Use Management of Western U.S. Rangelands: Wild Horses, Wildlife, and Livestock. Hum. Wildl. Interact. 2018, 12, 4. [Google Scholar] [CrossRef]
  33. Beever, E.A.; Aldridge, C.L. Influences of Free-Roaming Equids on Sagebrush Ecosystems, with a Focus on Greater Sage-Grouse. In Greater Sage-Grouse: Ecology and Conservation of a Landscape Species and Its Habitats; University of California Press: Oakland, CA, USA, 2011; pp. 272–290. [Google Scholar]
  34. Peco, B.; Borghi, C.E.; Malo, J.E.; Acebes, P.; Almirón, M.; Campos, C.M. Effects of Bark Damage by Feral Herbivores on Columnar Cactus Echinopsis (=Trichocereus) terscheckii Reproductive Output. J. Arid Environ. 2011, 75, 981–985. [Google Scholar] [CrossRef]
  35. Marshal, J.P.; Bleich, V.C.; Andrew Marshal, N.G.; Bleich, J.P.; Andrew, V.C.; Bighorn, S.N.; Andrew, N.G. Evidence for Interspecific Competition between Feral Ass Equus asinus and Mountain Sheep Ovis canadensis in a Desert Environment. Wildl. Biol. 2008, 14, 228–236. [Google Scholar] [CrossRef]
  36. Attum, O.; El Noby, S.K.; Hassan, I.N. The Influence of Landscape Characteristics and Anthropogenic Factors on Waterhole Use by Vulnerable Nubian Ibex Capra nubiana. Oryx 2009, 43, 564. [Google Scholar] [CrossRef]
  37. Soultan, A.; Shah, M.S.; Almalki, A.M.; Browne, S.; Arcilla, N. Arabian Partridge Abundance and Limiting Factors at the Northern Boundary of Its Range. Birds 2023, 4, 179–189. [Google Scholar] [CrossRef]
  38. Grevemeyer, B.; Vandenplas, M.; Beigel, B.; Cho, E.; Willingham, A.; Verma, A. Detection of Leptospiral DNA in the Urine of Donkeys on the Caribbean Island of Saint Kitts. Vet. Sci. 2017, 4, 2. [Google Scholar] [CrossRef]
  39. Clancy, C.; Cooke, F.; Raw, Z. Entanglement, Autonomy and the Co-Production of Landscapes: Relational Geographies for Free-Roaming ‘Feral’ Donkeys (Equus asinus) in a Rapidly Changing World. Geoforum 2021, 123, 66–77. [Google Scholar] [CrossRef]
  40. Robin, L. Culling and Care: Ferals, Invasives and Conservation Icons in Australia. Aust. Zool. 2017, 39, 103–113. [Google Scholar] [CrossRef]
  41. Zeng, B. Camel Culling and Carbon Emissions in Rangelands in Central Australia. J. Environ. Plan. Manag. 2015, 58, 270–282. [Google Scholar] [CrossRef]
  42. Roberts, M.; Cresswell, W.; Hanley, N. Prioritising Invasive Species Control Actions: Evaluating Effectiveness, Costs, Willingness to Pay and Social Acceptance. Ecol. Econ. 2018, 152, 1–8. [Google Scholar] [CrossRef]
  43. Mallon, D.; Budd, K. Regional Red List Status of Carnivores in the Arabian Peninsula; IUCN, Gland & Environment and Protected Areas Authority: Sharjah, United Arab Emirates, 2011; pp. 1–49. [Google Scholar]
  44. Symes, A.; Taylor, J.; Mallon, D.; Porter, R.; Simms, C.; Budd, K. The Conservation Status and Distribution of the Breeding Birds of the Arabian Peninsula; IUCN, Gland & Environment and Protected Areas Authority: Sharjah, United Arab Emirates, 2015. [Google Scholar]
Figure 1. The study area (black area inside the circle in the right panel) and the location of the surveyed sites (black dots; n = 109) within the boundary (white polygon in the left panel) of AlUla region. The black polygons represent the nature reserves and the three areas in the region; abbreviations are defined as per the text.
Figure 1. The study area (black area inside the circle in the right panel) and the location of the surveyed sites (black dots; n = 109) within the boundary (white polygon in the left panel) of AlUla region. The black polygons represent the nature reserves and the three areas in the region; abbreviations are defined as per the text.
Animals 13 02885 g001
Figure 2. Association between the calculated minimum population number of feral donkeys and the landscape variables on the left (a1,a2) and the estimated spatial distribution of feral donkeys modeled using Inverse Distance Weighted interpolation (IDW) on the right (b).
Figure 2. Association between the calculated minimum population number of feral donkeys and the landscape variables on the left (a1,a2) and the estimated spatial distribution of feral donkeys modeled using Inverse Distance Weighted interpolation (IDW) on the right (b).
Animals 13 02885 g002
Table 1. The minimum population number (MPN) of feral donkeys in the study area.
Table 1. The minimum population number (MPN) of feral donkeys in the study area.
SiteMPN
Khaybar338
HUBR315
Jabal ElWard235
HZNR169
WNNR78
GNR3
Tayma0
Table 2. The model parameters for the best model.
Table 2. The model parameters for the best model.
EstimateSEZ Valuep-Value
(Intercept)−1.1610.419−2.7730.0001
Human presence−2.2950.488−4.7000.0001
Human presence (quadratic term)−0.6710.187−3.5810.0001
Vegetation cover2.8360.5465.1910.0001
Vegetation cover (quadratic term)−0.8610.198−4.3520.0001
Topographic roughness0.0850.0970.8670.3861
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Soultan, A.; Darwish, M.; Al-Johani, N.; Abdulkareem, A.; Alfaifi, Y.; Assaeed, A.M.; El-Bana, M.; Browne, S. Feral Donkey Distribution and Ecological Impacts in a Hyper-Arid Region. Animals 2023, 13, 2885. https://doi.org/10.3390/ani13182885

AMA Style

Soultan A, Darwish M, Al-Johani N, Abdulkareem A, Alfaifi Y, Assaeed AM, El-Bana M, Browne S. Feral Donkey Distribution and Ecological Impacts in a Hyper-Arid Region. Animals. 2023; 13(18):2885. https://doi.org/10.3390/ani13182885

Chicago/Turabian Style

Soultan, Alaaeldin, Mohammed Darwish, Nawaf Al-Johani, Ayman Abdulkareem, Yousef Alfaifi, Abdulaziz M. Assaeed, Magdy El-Bana, and Stephen Browne. 2023. "Feral Donkey Distribution and Ecological Impacts in a Hyper-Arid Region" Animals 13, no. 18: 2885. https://doi.org/10.3390/ani13182885

APA Style

Soultan, A., Darwish, M., Al-Johani, N., Abdulkareem, A., Alfaifi, Y., Assaeed, A. M., El-Bana, M., & Browne, S. (2023). Feral Donkey Distribution and Ecological Impacts in a Hyper-Arid Region. Animals, 13(18), 2885. https://doi.org/10.3390/ani13182885

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop