Next Article in Journal
Individuals with Psychopathic Traits and Poor Attitudes towards Animals Can Recognise Infant Features But Give Them Reduced Attentional Priority
Previous Article in Journal
Use of Probiotics in Intravaginal Sponges in Sheep: A Pilot Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Potential for New Donkey Farming Systems to Supply the Growing Demand for Hides

School of Agriculture, Policy and Development, University of Reading, Reading RG6 6AR, UK
*
Author to whom correspondence should be addressed.
Animals 2020, 10(4), 718; https://doi.org/10.3390/ani10040718
Submission received: 18 February 2020 / Revised: 24 March 2020 / Accepted: 9 April 2020 / Published: 20 April 2020
(This article belongs to the Section Equids)

Abstract

:

Simple Summary

The demand for donkey hides for gelatin used in Traditional Chinese Medicine, called ejiao, has increased greatly over recent years, resulting in high and rapidly increasing prices. This has put pressure on donkey populations globally and led to theft and illegal trade in donkeys, resulting in concerns for donkey welfare and the livelihoods of those who rely on donkeys. New donkey farming systems have been set up in China in response. We model how quickly these systems may be able to meet the demand for donkey hides. Results show that it will take at least 10–15 years or more for new donkey farming systems to be able to supply the demand for hides. This means that prices of donkeys and donkey hides will continue to increase with continued thefts and illegal trade. This will have further negative impacts on donkey welfare and on the livelihoods of those who rely on donkeys. Chinese ejiao producers will continue to try to source donkey hides from around the world, putting even greater pressure on donkey populations globally.

Abstract

The demand for donkey hides for ejiao, a Traditional Chinese Medicine, has resulted in rapidly increasing prices for donkey hides and donkeys. This has put pressure on donkey populations globally and has implications for donkey welfare and the livelihoods of those who rely on donkeys as working animals. The aim of the research was to explore the feasibility of setting up new donkey farming systems to supply the rising demand for ejiao using a system dynamics model of donkey production. Results show that the size of the initial female breeding herd, reproductive performance, age of reproduction, percentage of female births and average breeding life of donkeys are key variables affecting the time to build up the donkey population to supply the demand for hides, which will be at least ten to fifteen years. The implications of this are: (i) prices for donkey hides will continue to increase, (ii) companies producing ejiao will use other ingredients, (iii) China will continue to source donkey hides from around the world, and (iv) there will be continued theft and illegal trade of donkeys and concerns for rural households reliant on donkeys for their livelihoods and adverse impacts on donkey welfare.

1. Introduction

Globally, donkey products, such as milk or meat, have generally been secondary to the use of donkeys for traction and transport. Nevertheless, since antiquity, donkey products have been attributed medicinal, rejuvenating and beautifying properties. In China, medicinal and rejuvenating effects have been attributed to a gelatin, called ejiao, produced from donkey hides. Today, ejiao is considered one of the most valuable products within Traditional Chinese Medicine (TCM) and is a very popular health tonic medicine in China [1,2] and for those using TCM in other countries, with rapidly increasing demand and prices over recent years. This is primarily due to increased disposable incomes and promotion of TCM (and ejiao) within China. Clinical uses of ejiao within TCM are to enrich the blood and stop bleeding (it is commonly taken by some women at/around the time of menstruation), improve the immune system and to treat insomnia and dizziness [3,4,5]. Ejiao is also ascribed anti-aging and rejuvenating effects [5,6]. Although there are numerous companies producing ejiao products, the previously state-owned Dong’e Ejiao company in Shandong Province is the market leader and premium market brand for ejiao. This company only uses donkey hides sourced in China and has set up its own large-scale donkey farms to facilitate this. It refers to donkey farming for ejiao as a new source of wealth of benefit to farmers and to the health and well-being of consumers.
Some have attributed the decline in the Chinese donkey population over the last three decades mainly to the increased production of ejiao [7,8,9], although increased industrialization and mechanization have certainly played their part because communities which previously relied on donkey traction now use motorized vehicles and machinery.
The most recent estimate for the size of the global donkey population is 45.8 million donkeys [10]. The main use of donkeys has traditionally been, and still largely is, for traction and transport of goods and people. For this reason, economic development with the mechanization of transport and agriculture has generally had dramatic impacts on donkey populations [11,12,13,14]. In Europe, the donkey population has declined by more than 80 percent since the Second World War. Countries with large donkey populations in the past, such as Greece and Italy, have seen the most dramatic declines of more than 95 percent [11]. In China, between 1990 and 2018, when the Chinese economy grew at an extraordinary pace, the donkey population decreased by 77 percent from 11.1 million to 2.5 million donkeys [10,15]. Given the large decline in donkey populations in other countries that have seen strong economic development, it does not seem justified to lay the primary cause of the decline in the China donkey population on the ejiao industry.
It is undeniable though that the strong and increasing demand for ejiao is putting pressure on donkey populations worldwide through increased prices for donkeys [16]. Prices of ejiao, donkey hides and donkeys have increased significantly over the last decade or so in particular. Reliable price information is lacking but all evidence suggests strongly increasing prices. Figure 1 shows prices of donkey hides between 2000 and 2015.
Anecdotal evidence suggests that prices have continued to increase. Hancock and Xueqiao (2018) put the value of a donkey hide in 2018 at 3000 yuan (USD 473) [17], while others report an increase from 20 yuan in the year 2000 to 3000 yuan by the end of 2017 [18]. Bloomberg News [19] report prices of donkey hides of up to 8000 yuan (USD 1160).
This rise in prices has been accompanied by increased efforts by companies in China to source supplies of donkey hides from around the world including from Africa and South America and exploration of the use of feral donkey populations in countries such as Australia. A number of governments have been sufficiently concerned about the impact of rising prices on their donkey populations that they have restricted exports of donkeys or donkey hides to protect them. For example, Burkina Faso, Uganda, Tanzania, Botswana, Niger, Mali and Senegal have put restrictions on the trade of donkeys and/or donkey hides as increased prices for hides have led to reports of widespread thefts of donkeys [18,20]. This loss of donkeys for hides has implications for the livelihoods of the rural poor in Africa and other countries given the importance of donkeys as working animals used for a variety of purposes from cultivation of land to transport of materials and people [21]. In addition, there are concerns for the welfare of donkeys, with many documented stories of the abuse and suffering of donkeys as a result of the trade in donkey hides [20,22]. Thus, although the major demand for ejiao comes from within China, it has far reaching global implications for donkey populations.
In this paper, we focus on the potential to supply donkey hides (not meat, milk or other donkey derived products) from donkey farming systems. New donkey farms for the supply of hides have been set up in China in recent years [23]. It is reported that, in recent years, 15 provinces and 22 cities have announced their intention to provide subsidies to stimulate donkey breeding [17]. Donkeys are relatively costly to breed because of their long gestation period, which is a disincentive to farmers to breed them. A representative from one of the market leaders of ejiao production is said to have indicated that, by 2020, his company will be able to meet their basic needs for hides through their donkey farms [17]. These donkey farming systems to produce hides are unlike any seen before globally in terms of their scale, intensity, technology and level of investment.
The main aim of this paper is twofold. Firstly, we bring together information on donkey reproduction, donkey farming and the ejiao market. Secondly, we use this information in a system dynamics model to assess the potential for donkey hide production from a farming system in China to meet demand in the short and longer term.
System dynamics models help to clarify, understand and simulate the behaviour of complex systems, which is often highly counterintuitive [24]. Livestock systems are complex, with feedback loops and often counterintuitive dynamic impacts. System dynamics modelling has been used in a number of papers to better understand herd dynamics and manage farming systems [25,26,27,28,29,30,31,32]. However, to date, to our knowledge, system dynamics modelling has not been applied to the evaluation of the dynamic behaviour of a donkey farming system.
The primary focus of system dynamics models is not on making precise predictions but on making imprecise, conditional projections of dynamic behaviour [33]. In the case of the donkey farming system, this approach is particularly useful as the dynamic behaviour of the system depends on a number of reproduction and farming parameters that are currently uncertain. In this situation of uncertainty, conditional “what if” projections under different scenarios can be particularly insightful.

2. Materials and Methods

2.1. Approach

The focus of this paper is on the potential of a donkey farming system to meet the demand for donkey hides (not including meat, milk or other donkey-derived products) from the ejiao market. Information on ejiao production and demand is scarce and mainly comes from the grey literature and media reports [16]. In recent years, information on donkey reproduction has improved due to an increased interest in donkey breeding, mainly in Europe and the Americas. Nevertheless, compared to horses and other domesticated animals, knowledge about the reproductive physiology of donkeys and donkey breeding management is still very limited [34,35,36,37,38,39,40]. In addition, intensive donkey farming is a relatively recent development, mainly for the production of donkey dairy products in Europe and for donkey hides but also meat and dairy in China [23,34,41,42,43,44,45]. To date, intensive donkey farming is still under-researched and knowledge on donkey farming systems is very limited. System dynamics (SD) modelling with its focus on conditional projections of dynamic behaviour rather than on precise predictions is therefore particularly well suited to the study of a donkey farming system to meet the demand for donkey hides for the ejiao industry. See Turner et al. for an overview of SD modelling applied to agricultural and natural resource issues [46].
Two SD models were built and run using Stella® Architect 1.3.1 (isee systems 2017). Monthly time steps were chosen to make the interpretation of changes in parameters such as gestation length or age at slaughter more intuitive (e.g., 13 months rather than 1.08 years). Parameters such as mortality, which are generally defined as annual rates, can be entered into the model as annual rates and are converted into monthly rates within the model. Interfaces are included in the model to facilitate communication of the results to interested stakeholders.
In a first step, a SD model is applied to assess how gestation length, conception rate, abortion rate and the percentage of twins affect the reproductive performance of female donkeys (jennies). Reproductive performance is defined as the average number of foals per jenny per year in a stable population. The model is run over 240 months (i.e., 20 years) to ensure that a steady state is reached. The results of the last 12 month period are presented and discussed.
In a second step, the output from the reproductive performance model is included in a SD model of a donkey farming system. The objective of the model is to explore the impact of reproductive performance, farming characteristics and hide demand on (a) the time it takes to build up a farming system to meet donkey hide demand for ejiao production, (b) the size of the required breeding population and (c) the total herd size. The model is run over 600 months (i.e., 50 years).

2.2. Reproductive Performance Model

In the first step, reproductive performance of a stable population of female donkeys is modelled based on gestation length, conception rate, abortion rate and the percentage of twins. Here, reproductive performance is the average number of foals per jenny per year (Appendix A for the equations of the reproductive performance model and Appendix Figure A1 for the user interface of the model). Figure 2 shows the stock and flow diagram of the model.
The population is constant at 1000 jennies (chosen as a convenient reasonably substantial round figure from which outcomes for larger or smaller starting populations can easily be extrapolated and which reflects the size of some enterprises already set up in China). The stock of jennies is divided into pregnant jennies (JP) and non-pregnant jennies (JNP). The system is initialised with 900 pregnant jennies and 100 non-pregnant jennies. Conception (C) is the flow of non-pregnant jennies to pregnant jennies and is determined by the conception rate. The stock of pregnant jennies is a conveyor stock. The time it takes jennies to move through the pregnant jennies conveyor stock is set by the gestation parameter. Abortions (AB) are a leakage from the conveyor stock and are determined by the abortion rate. The outflow of the pregnant jennies conveyor stock is giving birth (GB).
J P t = J P t d t + C t G B t A B t
The abortion outflow and the giving birth outflow make up the inflow into the non-pregnant jennies stock in the next period. Annual reproductive performance (RP) is defined as the total number of births (TB) over 12 months divided by the total number of jennies (TJ) in the steady state (i.e., R P = ( i = 1 12 T B i ) / T J ).
Knowledge about the reproductive physiology of donkeys and donkey breeding management is incomplete. The majority of the published research has been carried out in Europe [11]. Donkey farming in China mainly uses Chinese breeds with the Dezhou breed the preferred, but not the only, breed used by donkey farms that supply the hide market [43,45,47]. Where information is available for Chinese breeds, more weight is given to the values reported for Chinese breeds than to information on other donkey breeds. For each of the four main inputs into the reproduction performance model, low, middle and high parameter values are defined based on the literature, wherever possible.
Gestation lengths have been studied for a number of donkey breeds in recent years [12,34,35,38,48,49,50] but not for Chinese breeds. The reported average gestation length ranges from 353 days [12] to 371 days [48,49,50]. The shortest and longest gestation are reported by Galisteo and Perez-Marin [35] as 331 and 421 days, respectively. Simulations are run with average gestation lengths of 11.5, 12 and 12.5 months, equivalent to 350, 365 and 380 days, respectively.
Reported conception rates within and between different breeding methods (natural breeding, hand breeding, artificial insemination (AI) with fresh, cooled or frozen-thawed semen) vary widely. In a large-scale farming system, the use of AI is likely to play a major role. Outside China, AI for donkeys with frozen-thawed semen has only recently received attention [37,40,51,52,53]. In a review of the use of AI in equids in China between 1935 and 2012, Deng et al. [43] report conception rates between 52.3 percent and 73.2 percent, whilst Chang et al. [54] report conception rates between 76.0 percent and 70.8 percent for timed AI with frozen-thawed semen, which is high compared to studies from outside China. China has a long history of using AI in donkeys going back to the 1940s and AI use has been widespread since the 1950s [43]. Wu [55] reports that using oestrus synchronisation and other techniques, the conception rates for Dezhou donkeys in the National Black Donkey Breeding Center in Shandong, China, have increased from 45 percent in 2013 to 92 percent in 2016.
Based on the literature on Dezhou donkeys, we use a middle value for the conception rate of 0.70, which is at the lower end of the conception rates reported for Dezhou donkeys. We use a conservative value for two main reasons. In the model, after birth of a foal, jennies are immediately subject to insemination again, i.e., during foal heat. Foal heat, the first heat after giving birth, usually occurs in the first two weeks after a jenny has given birth to a foal [51]. In donkeys, foal heat can be used successfully but with slightly lower conception rates than in the following cycles [38,49]. Secondly, the conception rates reported in the literature are mainly based on controlled research studies. It is likely that conception rates in a (less controlled) large-scale farming system will be lower. On the other hand, in the literature, conception rates are often defined by the oestrous cycle, which is approximately 25 days for most donkey breeds studied [39]. The model defines the conception rate on a monthly basis. A low value of 0.60 foals per jenny per year and a high value of 0.80 are used in our simulations.
Information on the abortion rate is scarce. Most studies on donkey reproduction do not report the number of abortions, with a few notable exceptions [12,34,37,49]. In Carluccio et al. [49], one out of 73 pregnancies, i.e., less than two percent, is reported to have ended in abortion. The abortion occurred before day 45 of the pregnancy. Nervo et al. [34] report eight abortions out of 38 pregnancies, i.e., over 20 percent. Abortions occurred between days 236 and 342. Wang et al. [45] report an unusually high number of abortions on an intensive donkey farm in China. Over a period of several months, 61 of approximately 500 pregnant jennies, i.e., approximately twelve percent, aborted their foals between the 200th and 300th day of pregnancy. This was reported as an unusual event, suggesting that abortion rates on the farm are usually substantially lower. Based on the literature, abortion rates are likely to be somewhere between 2 percent and 20 percent. We use 10 percent as the middle value, 5 percent as the low value and 15 percent as the high value for abortion rates.
Little is known about the potential for donkeys to carry viable twins. Multiple ovulations are common in donkeys. The percentage of multiple ovulations reported are between 5 and 70 percent [39]. Yang et al. [56] find that depending on the season, 16 to 34 percent of Dezhou black donkey jennies experience double ovulation. Often double ovulations do not result in viable twin pregnancies [34]. A widespread practice is that if twin embryos are detected, one of them is crushed [49]. By contrast, in the study carried out by Galisteo and Perez-Marin [35], 3 out of the 58 pregnancies were twin pregnancies that were carried to term and resulted in viable foals. It would seem possible that the twin percentage could be increased in an intensive farming system. For sheep in 1990, 101 lambs were weaned per 100 ewes presented for breeding. By 2015 the number of lambs weaned per ewe had increased to 125 [25]. We base the middle value on Galisteo and Perez-Marin [35] assuming 5 out of 100 pregnancies are twins. The low value assumes no twins, i.e., the assumption that twin embryos are crushed. A high value of 10 percent of pregnancies being twin pregnancies is used. This assumption would require advances in viable multiple births similar to those seen in sheep.

2.3. The Donkey Farming System Model

In the second step, a model of a farming system is applied to shed light on the likely time frame for developing a large-scale donkey farming system as well as on the herd characteristics of such a system. Figure 3 shows the interface of the model (Appendix C for the equations of the farming systems model and Appendix Figure A2 for the model diagram).
In the farming system model, stocks of male and female donkeys are modelled as different sub-models. Each sub-model has four equivalent stocks: (1) a conveyor stock from birth to slaughter maturity (BS), (2) a reservoir stock at slaughter age (SA), (3) a conveyor stock from slaughter age to reproduction age (RA) and (4) a reservoir stock of breeding donkeys (BD).
These stocks are determined by the following relationships.
B S t = B S t d t + B t S M t M t
where B are births, SM are donkeys reaching slaughter maturity and M is mortality.
S A t = S A t d t + S M t R t S t
where R are replacements for breeding and S are animals slaughtered.
R A t = R A t d t + R t M R t M t
where MR are mature replacements
B D t = B D t d t + M R t M t C t
where C are animals culled at the end of their breeding life.
The model minimises the time it takes to build up a breeding donkey population that can supply a target number of hides on an annual basis. In the build-up phase of the farming systems, all female offspring enter the breeding herd, while most of the male offspring are slaughtered for their hides. The model approximates the jenny requirement for the target hide production and compares this to the replacements already in the pipeline plus the current number of breeding jennies, adjusted for future culls. When the required number of breeding jennies is reached, only jennies required to keep the breeding herd stable enter the breeding herd. The hides of donkeys that are culled from the breeding herd are also used for ejiao production. Donkey hides in the mortality flows are assumed not to count towards hide production.
The variables included in the interface (see Figure 3) are presented in this section.
The target hide variable determines the number of hides per year that the farming system is built to achieve. Production of ejiao is estimated at between five and six thousand tonnes annually [4,18,57]. Information on how much ejiao can be produced from one donkey hide varies widely. Based on these different sources, we estimate that on average one donkey hide produces between 1.25 kg [18] and 2.65 kg (report on Dong’e Ejiao (DEEJ), a large producer of top-quality ejiao on a Chinese platform (in Chinese) http://www.sohu.com/a/129809208_135357) of ejiao. These estimates, together with an estimated annual demand for ejiao production of between five and six thousand tonnes, translate into a donkey hide demand of somewhere between two and five million hides per annum. We run our model with an annual target hide production of 2.0, 3.5 and 5.0 million hides a year.
The farming system starts with an initial female breeding stock. Intensive farming of donkeys is already happening in China [23,44,47]. The latest donkey population estimate is 2.5 million donkeys in China in 2018 [15]. Of these, probably approximately one million are adult female donkeys. It is unknown how many of these are kept as pack and draught animals and how many are part of a donkey farming system. We run simulations with initial female breeding stocks of 50,000, 100,000 and 200,000 jennies, i.e., approximately five, ten and twenty percent of the total adult female donkeys in China.
The reproductive performance values are based on our analysis in step 1. Donkeys enter puberty between 1 and 2 years of age [36,51]. In many research studies, the youngest breeding animals included are two or three years of age [34,35,38,39,53,54,55,56]. The oldest animals in the non-Chinese studies are often between fifteen and twenty years but, for the two Chinese research studies, the oldest jennies were seven and ten years [54,56]. Information on farming systems is scarce but Fanelli et al. [40] studied a donkey dairy farm with jennies aged between three and twenty years old. On the donkey dairy farms studied by Dai at al. [41], the oldest animal was 29 years. We simulate the model with donkeys entering the breeding herd at 24, 30 and 36 months. The literature suggests that donkeys can enjoy a long breeding life, with animals between the ages of fifteen and twenty years regularly used for breeding, but the two Chinese studies do not use animals over ten years of age. Based on this information, the average breeding life from entering the breeding herd until culling is set at 10, 12 and 15 years.
In the build-up phase of a donkey farming system, having more female offspring is desirable. Sperm sexing technology exists and has been applied to and assessed for dairy and beef cattle but more recently also to polo horses [58,59]. In a recent study on sperm sexing in donkeys, Domínguez et al. [59] report that after the sexing procedure, 90 percent of sperm has the X chromosome. Sexing sperm can negatively impact conception rates although the new sexing technique had no impact on conception rates in horses. The authors report that research is under way to assess the impact on conception rates in donkeys. Since the technology is not widely used yet, we simulate the model with the percentage of female foals of 50, 55 and 65 percent.
Those donkeys not retained for breeding are slaughtered for their hides. Wu [55] reports that, in traditional systems, it takes 24 months to raise a donkey for its hide but that the period has been shortened to 16 to 18 months. We simulate the model with slaughter ages of 17, 22 and 24 months.
The model includes two mortality rates. Foal mortality is the percentage of donkeys that die before they reach slaughter maturity. Mortality of breeding donkeys is the percentage of breeding animals that die per year. Mortality of breeding donkeys is separate to and in addition to the cull rate of breeding animals. The cull rate is derived from the average breeding life variable. An average breeding life of ten years translates into an annual cull rate of 10 percent.
Data on mortality are very scarce. Few research papers report mortalities. It is not clear whether the lack of reporting of mortality is due to the fact that none of the donkeys died during the study or whether animals that died were excluded from the studies. Quaresma et al. [13] carried out a survey of donkey owners whose donkeys were registered in a studbook. The authors found that 725 of the 760 donkeys born between 1982 and 2012 and registered in the studbook were alive at the end of 2012. This means that less than five percent had died and suggests an annual mortality rate of substantially below five percent. Data on the donkey population and deaths are available for Ethiopia from the 2017/2018 agricultural survey [60]. The data suggest an annual mortality rate of approximately five percent. The deaths are likely to mainly be of older donkeys. In the farming system, these would be equivalent to the culled animals. Mortality in the farming system model includes deaths that occur before and in addition to the animals culled at the end of their breeding life. The Ethiopian figures therefore also suggest a mortality rate for breeding donkeys well below five percent. We use foal mortality figures of 10, 7 and 5 percent and additionally adult donkey mortality figures of 4, 3 and 2 percent per year.
The model also includes the jenny to jack ratio. This ratio will depend on the breeding method used, with more jacks needed if natural breeding is predominant. We use jenny to jack ratios of 10, 20 and 30 to one in our simulations.

3. Results

3.1. Reproductive Performance

Table 1 summarises the parameter values and the results for the three reproductive performance scenarios.
The middle scenario results in a reproductive performance of 0.95, or 95 foals being born per 100 jennies in the breeding herd. The reproductive performance of the less favourable scenario is 83 foals per 100 jennies in the breeding herd, i.e., a 13 percent lower reproductive performance than in the middle scenario. For the more favourable scenario, the model indicates that 109 foals are born per 100 jennies in the breeding herd—14 percent higher than in the middle scenario.
Table 2 shows the results of sensitivity analysis. We vary each parameter at a time and set all the other parameters at the values of the middle scenario.
The results show that the twin percentage has the largest impact followed by gestation length. The reproductive performance rates from three main scenarios are used in the farming system model described in the following section.

3.2. The Donkey Farming System

We simulate the model for nine different scenarios. We combine the three different target hide settings with the less favourable, middle and more favourable scenarios. In the middle scenario, all parameter values are set to their middle values as reported in Section 2. The middle scenario is most likely to represent what is currently achievable. However, it seems appropriate to re-emphasize here that the aim of this study is not to make precise predictions but to produce conditional projections of dynamic behaviour under “what if” scenarios—what if a donkey farming system could achieve the reproductive and farming performances reported in the literature.
Table 3 summarizes the parameters used in the simulations and the main results.
One important question is how long it might take to set up a farming system that can meet the demand for hides. We report the time in years, not months, to reflect the uncertainty present. We also report the total number of donkeys in the farming system and the number of breeding jennies.
Figure 4 shows that with regard to the time period it takes to reach the target hides, the different farming system performance scenarios have a greater impact than the target hide scenarios.
The target hide output is reached in 20 to 25 years in the middle scenario for all three target hide settings. It takes 38 to 46 years (over 80 percent longer) in the less favourable scenario and it takes 11 to 14 years (between 45 and 50 percent shorter) in the more favourable simulation setting to achieve the three target hide outputs. Comparing the 5.0 million hide scenario with the corresponding 3.5 million hide scenario extends the time taken to produce the target output by between 7 and 11 percent. Comparing the respective 2.0 million target hide and 3.5 million target hide scenarios, the time taken reduces by between 11 and 15 percent.
Figure 5 shows that in terms of the number of animals required in the farming system to produce the target hides, the number of animals is affected more by the variation in target hides than by the variation in farming system performance.
In the middle scenario for farming system performance, farming systems with 6.5 million donkeys, 11.4 million donkeys and 16.4 million donkeys are required for annual outputs of 2.0 million, 3.5 million and 5.0 donkey hides, respectively. Under the less (more) favourable scenario, the number of donkeys varies from 8.1 (5.8) million for 2.0 million hides to 20.2 (14.6) million for 5.0 million hides.
Table 4 shows the results of sensitivity analysis for the individual parameters. The impact of varying individual parameters on the time it takes to reach the target hides and the number of donkeys in the farming system are presented in percentage terms. The time taken to reach the target hide output is most sensitive to changes in the initial breeding herd and reproductive performance. The donkey numbers are most sensitive to the target hide values and the reproductive performance, with the latter mainly impacting on the number of breeding jennies as would be expected. The reproductive age, the lengths of the average breeding life and the percentage of female births have a substantial impact on the time it takes to reach the target hides but little impact on total donkey number and little or no impact on the number of breeding jennies. The contributions of other variables to changes in the time or donkey numbers are small.

4. Discussion

Our simulations suggest that it is likely to take at least ten to fifteen years, if not much longer, to build up a Chinese farming system to supply the current demand for donkey hides from the ejiao industry. Our model results indicate that, for an annual production of 2.0 million donkey hides, the number of donkeys in the donkey farming system needs to be 5.8 million under favourable reproduction and farming conditions, approximately double the current total donkey population in China. For an annual production of 5.0 million donkey hides, our model suggests that 14.6 million donkeys would be required under the favourable farming system performance scenario—or almost six times the current donkey population in China. Of course, in 10–15 years, demand for hides could be significantly higher than currently (e.g., due to increased population and disposable incomes in China), providing an even greater challenge to donkey farming systems to supply the market. Indeed, the demand for the use of TCM (including ejiao) has received a likely boost from the WHO inclusion of a chapter on TCM in its highly influential International Classification of Diseases (11th edition) which helps set the medical agenda globally and may contribute to the acceptance of TCM alongside conventional (western) medicine [9].
The assumptions underlying the more favourable estimates are rather optimistic. Although there are opportunities for better control and more efficient management on large-scale intensive donkey farms, large farms pose management challenges and disease risks [23,47]. For example, Yang et al. [61] report on an equine influenza outbreak on a 300 head donkey farm in Shandong which resulted in 25 percent mortality. Wang et al. [45] study the increase in abortions due to Salmonella abortus equi on a farm with over 1000 jennies. These examples show that it may take time to achieve reproductive and other performance levels found in studies carried out on smaller, less-intensive farms. From the above analysis, we can conclude that it is likely to take much longer than a decade for China to meet the demand for donkey hides from a domestic donkey farming system. This has negative implications for donkey welfare in that the economic incentive will continue for theft, illegal trade and trade (including transport and slaughter of donkeys) which gives little or no consideration to the welfare of the donkeys involved (see The Donkey Sanctuary, 2019 for some documented case studies). In the longer term, new large-scale donkey farming systems will have some economic incentive to ensure reasonable levels of donkey welfare and there is an important role for organizations such as the International Office for Animal Health and the Food and Agricultural Organization to agree international health and welfare standards for these systems. There is also an important role for non-governmental animal welfare charities to advise companies setting up such systems in relation to welfare standards and the commercial benefits of good donkey welfare (and there is already evidence of this happening in China).
The system dynamics model developed in this study does not include any economic considerations. The main reason is the lack of information on the costs of donkey production. Further, little reliable information is available on the price that donkey farmers receive for a hide. There are likely to be economic and business challenges to building up donkey farming systems. It takes approximately three years for any revenue from donkey farming (approximately one year of gestation and two years for the foal to reach slaughter weight). As a result, substantial initial investment is needed, with donkey farming only becoming profitable after several years. On this basis, calls for subsidies have been made by ejiao industry representatives as far back as 2010 [7]. The economic incentive for donkey farmers to increase breeding activities could be provided by subsidies or by further increases in the price of donkey hides. The latter would make investment in donkey farming more attractive but would not address the challenge of the long time period required before newly-established farms become profitable.
As already noted, the donkey hide targets in the simulations are subject to some uncertainty. Firstly, reliable figures on ejiao production are not available. Secondly, different sources on how many donkey hides are required to produce one tonne of ejiao give quite different answers. Thirdly, the potential of imports of donkey hides is also uncertain. According to the latest estimates from the Food and Agriculture Organization, there are just over 40 million donkeys outside China [10]. In most countries, the hides of donkeys do not seem to be used upon their death. The most detailed figures on population and deaths are published by the Central Statistical Agency of Ethiopia [60]. In 2017/2018, the total Ethiopian donkey population was estimated at 8.8 million and deaths were estimated at 481 thousand, i.e., approximately 5 percent of the donkeys died in 2017/2018. Using the same ratio for the total population outside China, there would be a natural annual supply of donkey hides of 2.2 million. However, there are significant economic, cultural and governance barriers to the export of donkey hides. Nevertheless, there is some potential for donkey hide trade to contribute to the supply in the short and longer term, but trade is likely to remain relatively small compared to overall demand. Finally, the current use of donkey hides also depends on the size of the counterfeit market. Some sources suggest that up to 40 percent of ejiao is counterfeit [57]. Counterfeit production is a well-known problem and research on identifying genuine product is going on in China [2,3,4,62,63]. Our results assume that the totality of the production of ejiao is genuine, i.e., made from donkey hides only.
There are a number of possible extensions of this study. As noted above, economic factors could negatively impact on the development of a donkey farming system. The model can be extended to include economic information. A second possible extension is to include milk and/or meat production in the model. At least some of the farms in China produce hides, milk and meat [45].

5. Conclusions

The system dynamics model of donkey farming for hides has helped to explore the potential for new systems of donkey farming in China (and elsewhere) to supply the growing demand for hides for ejiao production. Reliable information on the likely productive parameters for donkeys bred in a large, intensive production system using the latest technologies is lacking, with some uncertainty surrounding key variables such as reproductive performance, mortality and breeding life. The model has enabled such uncertainties and their effects on the speed with which the donkey population and output of hides can be increased to be explored and assessed. The study has found that it is likely to take at least 10 to 15 years, if not longer, to build up a Chinese farming system to supply the (current) demand for donkey hides for the ejiao industry, depending on the extent to which the productive performance of donkey breeding herds can be maximized. The implications of this are that (i) prices for donkey hides, and for donkeys, are likely to continue to increase over the short and medium term; (ii) there will be a strong incentive for companies producing ejiao products to diversify their products and ingredients to reduce the amount of ejiao and, for example, include gelatin from other species and other ingredients to ‘stretch’ the amount of ejiao available across products; (iii) China will continue to explore non-domestic sources of donkey hides from around the world until China’s domestic supply of donkeys can be increased; and (iv) there will be a continued financial incentive for the theft and illegal trade of donkeys. Concern will also remain about the possible impact of this pressure on donkey populations for low-income rural households reliant on donkeys for their livelihoods and the potential adverse impact on the welfare of donkeys themselves.
There is a need to continue to monitor the market for ejiao and donkey hides and its impact on donkey populations globally. There is also a need to monitor the development of new donkey farming systems that produce hides and to ensure that these systems have satisfactory standards in relation to biosecurity, donkey welfare and environmental impacts.

Author Contributions

R.B. was Principal Investigator on the research and S.P. took the lead in the dynamic systems modelling of donkey populations. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Donkey Sanctuary, grant number H5414100.

Acknowledgments

The authors would like to thank Mofei Jia for assistance with translations from Chinese.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Reproduction model: Stella® model equations
JenniesNotPregnant(t) = JenniesNotPregnant(t - dt) + (NewCycle - conception) * dt
 INIT JenniesNotPregnant = 100
 INFLOWS:
  NewCycle = nextCycle+Abortions
 OUTFLOWS:
  conception = JenniesNotPregnant*conception_rate
JenniesPregnant(t) = JenniesPregnant(t - dt) + (conception - GivingBirth - Abortions) * dt
 INIT JenniesPregnant = 900
  TRANSIT TIME = gestation-1
  CAPACITY = INF
  INFLOW LIMIT = INF
 INFLOWS:
  conception = JenniesNotPregnant*conception_rate
 OUTFLOWS:
  GivingBirth = CONVEYOR OUTFLOW
  Abortions = LEAKAGE OUTFLOW
   LEAKAGE FRACTION = AbortionRate
AbortionRate = 0.2
conception_rate = 0.4
gestation = 12
nextCycle = GivingBirth
ReproductivePerformance = (TotalBirths*12)/(JenniesNotPregnant+JenniesPregnant)
TotalBirths = GivingBirth*(TwinPrecentage/100+1)
TwinPrecentage = 5
{ The model has 13 (13) variables (array expansion in parens).
 In root model and 0 additional modules with 0 sectors.
 Stocks: 2 (2) Flows: 4 (4) Converters: 7 (7)
 Constants: 5 (5) Equations: 6 (6) Graphicals: 0 (0)

Appendix B

Figure A1. Reproduction model interface.
Figure A1. Reproduction model interface.
Animals 10 00718 g0a1

Appendix C

Farming system model: Stella® model equations
BreedingHerdDynamics:
BreedingJacks(t) = BreedingJacks(t - dt) + (ReplacementJacks - MaleEndOfBreedingCull - JackMortality) * dt
 INIT BreedingJacks = InitialJennies/FemaleMaleBreedingRatio
 INFLOWS:
  ReplacementJacks = CONVEYOR OUTFLOW
 OUTFLOWS:
  MaleEndOfBreedingCull = IF TIME < AverageBreedingLengthM/2 THEN 0 ELSE BreedingJacks/AverageBreedingLengthM
  JackMortality = BreedingJacks*JackMortalityM
BreedingJennies(t) = BreedingJennies(t - dt) + (ReplacementJennies - FemalesEndOfBreedingCull - JennieMortality) * dt
 INIT BreedingJennies = InitialJennies
 INFLOWS:
  ReplacementJennies = CONVEYOR OUTFLOW
 OUTFLOWS:
 FemalesEndOfBreedingCull = IF TIME < AverageBreedingLengthM/2 THEN 0 ELSE BreedingJennies/AverageBreedingLengthM
 JennieMortality = JennieMortalityM*BreedingJennies
FemaletoReproductionAge(t) = FemaletoReproductionAge(t - dt) + (FemaleforReplacement - ReplacementJennies - FemaleMortalityMaturing) * dt
 INIT FemaletoReproductionAge = 0
  TRANSIT TIME = FemaleMaturingTimeReproduction
  CAPACITY = INF
  INFLOW LIMIT = INF
 INFLOWS:
  FemaleforReplacement = IF (ApproximateJennieRequirement > (BreedingJennies*((AverageBreedingLengthM-FemaleMaturingTimeReproduction)/AverageBreedingLengthM)+FemaletoReproductionAge))/(1-JennieMortalityM) THEN JenniesSlaughterAge ELSE (IF TargetHides> HistoricalHides THEN (BreedingJennies*((AverageBreedingLengthM-FemaleMaturingTimeReproduction)/AverageBreedingLengthM)+FemaletoReproductionAge) * (1/AverageBreedingLengthM+JennieMortalityM) ELSE BreedingJennies*(1/AverageBreedingLengthM+JennieMortalityM))
 OUTFLOWS:
  ReplacementJennies = CONVEYOR OUTFLOW
  FemaleMortalityMaturing = LEAKAGE OUTFLOW
  LEAKAGE FRACTION = 0
FemaleToSlaughterMaturity(t) = FemaleToSlaughterMaturity(t - dt) + (FemaleBeingBorn - FemaleMaturingToSlaughter - FemaleFoalMortality) * dt
 INIT FemaleToSlaughterMaturity = 0
  TRANSIT TIME = FemaleSlaughterAge-1
  CAPACITY = INF
  INFLOW LIMIT = INF
 INFLOWS:
  FemaleBeingBorn = BreedingJennies*ReproductivePerformance*PercentageFemaleBirths/100
OUTFLOWS:
  FemaleMaturingToSlaughter = CONVEYOR OUTFLOW
  FemaleFoalMortality = LEAKAGE OUTFLOW
   LEAKAGE FRACTION = FemaleNonBreedingMortality
JacksSlaughterAge(t) = JacksSlaughterAge(t - dt) + (MaleMaturingToSlaughter - MaleSlaughter - MaleforReplacement) * dt
 INIT JacksSlaughterAge = 0
 INFLOWS:
  MaleMaturingToSlaughter = CONVEYOR OUTFLOW
 OUTFLOWS:
  MaleSlaughter = JacksSlaughterAge-MaleforReplacement
  MaleforReplacement = BreedingJacks*1/AverageBreedingLengthM - (BreedingJacks-1/FemaleMaleBreedingRatio*BreedingJennies)
JenniesSlaughterAge(t) = JenniesSlaughterAge(t - dt) + (FemaleMaturingToSlaughter - FemaleforReplacement - FemaleSlaughter) * dt
 INIT JenniesSlaughterAge = 0
 INFLOWS:
  FemaleMaturingToSlaughter = CONVEYOR OUTFLOW
 OUTFLOWS:
  FemaleforReplacement = IF (ApproximateJennieRequirement > (BreedingJennies*((AverageBreedingLengthM-FemaleMaturingTimeReproduction)/AverageBreedingLengthM)+FemaletoReproductionAge))/(1-JennieMortalityM) THEN JenniesSlaughterAge ELSE (IF TargetHides> HistoricalHides THEN (BreedingJennies*((AverageBreedingLengthM-FemaleMaturingTimeReproduction)/AverageBreedingLengthM)+FemaletoReproductionAge) * (1/AverageBreedingLengthM+JennieMortalityM) ELSE BreedingJennies*(1/AverageBreedingLengthM+JennieMortalityM))
 FemaleSlaughter = JenniesSlaughterAge-FemaleforReplacement
MaleToReproductionAge(t) = MaleToReproductionAge(t - dt) + (MaleforReplacement - ReplacementJacks - MaleMortalityMaturing) * dt
 INIT MaleToReproductionAge = 0
  TRANSIT TIME = 1
  CAPACITY = INF
  INFLOW LIMIT = INF
 INFLOWS:
  MaleforReplacement = BreedingJacks*1/AverageBreedingLengthM - (BreedingJacks-1/FemaleMaleBreedingRatio*BreedingJennies)
 OUTFLOWS:
  ReplacementJacks = CONVEYOR OUTFLOW
  MaleMortalityMaturing = LEAKAGE OUTFLOW
   LEAKAGE FRACTION = 0
MaleToSlaughterMaturity(t) = MaleToSlaughterMaturity(t - dt) + (MaleBeingBorn - MaleMaturingToSlaughter - MaleFoalMortality) * dt
 INIT MaleToSlaughterMaturity = 0
  TRANSIT TIME = MaleSlaughterAge
  CAPACITY = INF
  INFLOW LIMIT = INF
 INFLOWS:
  MaleBeingBorn = BreedingJennies*ReproductivePerformance*(1-PercentageFemaleBirths/100)
 OUTFLOWS:
  MaleMaturingToSlaughter = CONVEYOR OUTFLOW
  MaleFoalMortality = LEAKAGE OUTFLOW
   LEAKAGE FRACTION = MaleNonBreedingMortality
“12monthshides” = TotalHides + DELAY(TotalHides,1, 0)+ DELAY(TotalHides,2, 0) + DELAY(TotalHides,3, 0)+ DELAY(TotalHides,4, 0) + DELAY(TotalHides,5, 0) + DELAY(TotalHides,6, 0) + DELAY(TotalHides,7, 0) + DELAY(TotalHides,8, 0) + DELAY(TotalHides,9, 0) + DELAY(TotalHides,10, 0) + DELAY(TotalHides,11, 0)
AdultDonkeyMortalityY = 0.1
ApproximateJennieRequirement = TargetHides/(ReproductivePerformance*(1-FemaleNonBreedingMortality)-JennieMortalityM)*(1+(1/FemaleMaleBreedingRatio)/2)
AverageBreedingLengthM = AverageBreedingLengthY*12
AverageBreedingLengthY = 10
CleanHides = FemaleSlaughter + MaleSlaughter
CulHides = FemalesEndOfBreedingCull + MaleEndOfBreedingCull
DeadDonkeys = FemaleFoalMortality + FemaleMortalityMaturing + JackMortality + JennieMortality + MaleFoalMortality + MaleMortalityMaturing
FemaleMaleBreedingRatio = 10
FemaleMaturingTimeReproduction = FemaleStartReproduction-FemaleSlaughterAge
FemaleNonBreedingMortality = 0.01
FemaleSlaughterAge = 20
FemaleStartReproduction = 26
HistoricalHides = DELAY(TotalHides,1, 0)
InitialJennies = 100000
JackMortalityM = -((1-JackMortalityY)^(1/12)-1)
JackMortalityY = AdultDonkeyMortalityY
JennieMortalityM = -((1-JennieMortalityY)^(1/12)-1)
JennieMortalityY = AdultDonkeyMortalityY
MaleMaturingTimeReproduction = MaleStartReproduction-MaleSlaughterAge
MaleNonBreedingMortality = FemaleNonBreedingMortality
MaleSlaughterAge = 15
MaleStartReproduction = 24
Overallmortality = DeadDonkeys*12/TotalDonkeys
PercentageFemaleBirths = 50
ReproductivePerformance = ReproductivePerformanceY/12
ReproductivePerformanceY = 0.8
TargetHides = TargetHidesY/12
TargetHidesY = 4200000
TotalDonkeys = BreedingJacks + BreedingJennies + FemaleToSlaughterMaturity + FemaletoReproductionAge + JacksSlaughterAge + JenniesSlaughterAge + MaleToReproductionAge + MaleToSlaughterMaturity
TotalHides = CleanHides + CulHides
{ The model has 57 (57) variables (array expansion in parens).
 In this module and 0 additional modules with 0 sectors.
 Stocks: 8 (8) Flows: 18 (18) Converters: 31 (31)
 Constants: 18 (18) Equations: 31 (31) Graphicals: 0 (0)
 }

Appendix D

Figure A2. Farming system model: Stock and flow diagram.
Figure A2. Farming system model: Stock and flow diagram.
Animals 10 00718 g0a2

References

  1. Ge, S.; He, T.T.; Hu, H. Popularity and customer preferences for over-the-counter Chinese medicines perceived by community pharmacists in Shanghai and Guangzhou: A questionnaire survey study. Chin. Med. (U. K.) 2014, 9, 1–8. [Google Scholar] [CrossRef] [Green Version]
  2. Zuo, H.-L.; Zhao, J.; Wang, Y.-T.; Xia, Z.-N.; Hu, Y.-J.; Yang, F.-Z. Identification of the Adulterated Asini Corii Colla with Cytochrome c Oxidase Subunit I Gene-based Polymerase Chain Reaction. Pharmacogn. Res. 2017, 9, 313–318. [Google Scholar]
  3. Zhang, W.; Cui, S.; Cheng, X.L.; Wei, F.; Ma, S. An optimized TaqMan real-time PCR method for authentication of ASINI CORII COLLA (donkey-hide gelatin). J. Pharm. Biomed. Anal. 2019, 170, 196–203. [Google Scholar] [CrossRef]
  4. Li, X.; Shi, F.; Gong, L.; Hang, B.; Li, D.; Chi, L. Species-specific identification of collagen components in Colla corii asini using a nano-liquid chromatography tandem mass spectrometry proteomics approach. Int. J. Nanomed. 2017, 12, 4443–4454. [Google Scholar] [CrossRef] [Green Version]
  5. Yang, F.; Song, J.; Liang, M.; Ma, F.; Mao, X.; Ma, C.W.; Zhang, W.; Huang, Z. Overview of Beverages with Anti-Aging Functions in Chinese Market. Rejuvenation Res. 2014, 17, 1–4. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, D.; Liu, M.; Cao, J.; Cheng, Y.; Zhuo, C.; Xu, H.; Tian, S.; Zhang, Y.; Zhang, J.; Wang, F. Effect of Colla corii asini (Ejiao) on D-Galactose Induced Aging Mice. Biol. Pharm. Bull. 2012, 35, 2128–2132. [Google Scholar] [CrossRef] [Green Version]
  7. Swift, A. Ejiao: Donkeys Are Paying the Price for China’s Rapid Progress and Lost Traditions. Available online: https://donkeytime.org/2016/12/18/ejiao-donkeys-are-paying-the-price-for-rapid-progress-and-lost-traditions/ (accessed on 3 March 2019).
  8. Cox, L. Dramatic decline in donkey populations blamed on consumer demand for “ejiao”. Vet. Rec. 2017, 181, 610. [Google Scholar]
  9. Nature Editorial Break with tradition. The World Health Organization’s decision to include traditional Chinese medicine in its global diagnostic compendium could backfire. Nature 2019, 570, 5. [Google Scholar] [CrossRef] [Green Version]
  10. Food and Agriculture Organization of the United Nations (FAO). Statistics Division FAOStat. Available online: http://www.fao.org/faostat/en/#data/QA (accessed on 2 August 2019).
  11. Camillo, F.; Rota, A.; Biagini, L.; Tesi, M.; Fanelli, D.; Panzani, D. The Current Situation and Trend of Donkey Industry in Europe. J. Equine Vet. Sci. 2018, 65, 44–49. [Google Scholar] [CrossRef]
  12. Crisci, A.; Rota, A.; Panzani, D.; Sgorbini, M.; Ousey, J.C.; Camillo, F. Clinical, ultrasonographic, and endocrinological studies on donkey pregnancy. Theriogenology 2014, 81, 275–283. [Google Scholar] [CrossRef]
  13. Quaresma, M.; Martins, A.M.F.; Rodrigues, J.B.; Colaço, J.; Payan-Carreira, R. Pedigree and herd characterization of a donkey breed vulnerable to extinction. Animal 2014, 8, 354–359. [Google Scholar] [CrossRef] [Green Version]
  14. Arsenos, G.; Gelasakis, A.I.; Papadopoulos, E.I. The status of Donkeys (Equus asinus) in Greece. J. Hell. Vet. Med. Soc. 2017, 61, 212–219. [Google Scholar] [CrossRef] [Green Version]
  15. National Bureau of Statistics of China. Available online: http://data.stats.gov.cn/english/easyquery.htm?cn=C01 (accessed on 2 August 2019).
  16. Bennett, R.; Pfuderer, S. Demand for Donkey Hides and Implications for Global Donkey Populations. Available online: https://econpapers.repec.org/paper/agsaesc19/289683.htm (accessed on 8 September 2019).
  17. Hancock, T.; Xueqiao, W. China Bets on Donkey Breeding to Curb Africa Imports. Available online: https://www.ft.com/content/bee2b852-0190-11e8-9650-9c0ad2d7c5b5 (accessed on 12 August 2018).
  18. Xinhua News Agency. China Cuts Tariff on Donkey Hides. CGTN. Business 21:49, 31 December 2017, 1. Available online: https://news.cgtn.com/news/3067444d35637a6333566d54/index.html (accessed on 15 August 2019).
  19. Bloomberg News Have a Donkey to Trade? China Has an Exchange for That. Available online: https://www.bloomberg.com/news/articles/2017-03-29/china-has-exchange-traded-donkeys (accessed on 12 August 2018).
  20. The Donkey Sanctuary. Under the Skin: The Emerging Trade in Donkey Skins and its Implications for Donkey Welfare and Livelihoods. Available online: https://www.thedonkeysanctuary.org.uk/sites/uk/files/2017-11/under_the_skin_report.pdf (accessed on 5 July 2018).
  21. The Brook Invisible Workers: The Economic Contributions of Working Donkeys, Horses, and Mules to Livelihoods. Available online: https://www.thebrooke.org/sites/default/files/Advocacy-and-policy/Invisible-workers-report.pdf (accessed on 8 February 2020).
  22. The Donkey Sanctuary. Under the Skin: Update on the Global Crisis for Donkeys and the People Who Depend on Them. Available online: https://www.thedonkeysanctuary.org.uk/sites/uk/files/2019-11/under-the-skin-report-revised-2019.pdf (accessed on 29 November 2019).
  23. Dong, J.; Gao, N.; Waller, A.S.; Cook, F.R.; Fan, S.; Yuan, D.; Du, Y.; Li, F.; Norimine, J.; Zhu, W. An outbreak of strangles associated with a novel genotype of Streptococcus equi subspecies equi in donkeys in China during 2018. Equine Vet. J. 2019, 51, 743–748. [Google Scholar] [CrossRef]
  24. Sterman, J.D. All models are wrong: Reflections on becoming a systems scientist. Syst. Dyn. Rev. 2002, 18, 501–531. [Google Scholar] [CrossRef]
  25. Farrell, L.J.; Tozer, P.R.; Kenyon, P.R.; Ramilan, T.; Cranston, L.M. The effect of ewe wastage in New Zealand sheep and beef farms on flock productivity and farm profitability. Agric. Syst. 2019, 174, 125–132. [Google Scholar] [CrossRef]
  26. Shane, D.D.; Larson, R.L.; Sanderson, M.W.; Miesner, M.; White, B.J. A deterministic, dynamic systems model of cow–calf production: The effects of breeding replacement heifers before mature cows over a 10-year horizon. J. Anim. Sci. 2017, 95, 4533–4542. [Google Scholar] [CrossRef]
  27. Lie, H.; Rich, K.M.; van der Hoek, R.; Dizyee, K. An empirical evaluation of policy options for inclusive dairy value chain development in Nicaragua: A system dynamics approach. Agric. Syst. 2018, 164, 193–222. [Google Scholar] [CrossRef]
  28. Mumba, C.; Skjerve, E.; Rich, M.; Rich, K.M. Application of system dynamics and participatory spatial group model building in animal health: A case study of East Coast Fever interventions in Lundazi and Monze districts of Zambia. PLoS ONE 2017, 12, 1–21. [Google Scholar] [CrossRef] [Green Version]
  29. Turner, B.L.; Rhoades, R.D.; Tedeschi, L.O.; Hanagriff, R.D.; McCuistion, K.C.; Dunn, B.H. Analyzing ranch profitability from varying cow sales and heifer replacement rates for beef cow-calf production using system dynamics. Agric. Syst. 2013, 114, 6–14. [Google Scholar] [CrossRef]
  30. Shane, D.D.; Larson, R.L.; Sanderson, M.W.; Miesner, M.; White, B.J. A deterministic, dynamic systems model of cow-calf production: The effects of the duration of postpartum anestrus on production parameters over a 10-year horizon. J. Anim. Sci. 2017, 95, 1680–1695. [Google Scholar] [CrossRef]
  31. Gebre, K.T.; Wurzinger, M.; Gizaw, S.; Haile, A.; Rischkowsky, B.; Sölkner, J. The use of system dynamics modelling methodologies in sheep breeding programs and management systems. Presented at the Conference on International Research on Food Security, Natural Resource Management and Rural Development, University of Natural Resources and Life Sciences, Vienna, Austria, 18–21 September 2016; pp. 1–4. Available online: https://repo.mel.cgiar.org/handle/20.500.11766/6180 (accessed on 15 December 2019).
  32. Diaz-Solis, H.; Kothmann, M.M.; Grant, W.E.; De Luna-Villarreal, R. Use of irrigated pastures in semi-arid grazinglands: A dynamic model for stocking rate decisions. Agric. Syst. 2006, 88, 316–331. [Google Scholar] [CrossRef]
  33. Duggan, J. Systems Dynamics Modeling with R.; Springer International Publishing: Cham, Switzerland, 2016. [Google Scholar]
  34. Nervo, T.; Bertero, A.; Poletto, M.; Pregel, P.; Leone, R.; Toffoli, V.; Vincenti, L. Field ultrasound evaluation of some gestational parameters in jennies. Theriogenology 2019, 126, 95–105. [Google Scholar] [CrossRef] [PubMed]
  35. Galisteo, J.; Perez-Marin, C.C. Factors affecting gestation length and estrus cycle characteristics in Spanish donkey breeds reared in southern Spain. Theriogenology 2010, 74, 443–450. [Google Scholar] [CrossRef] [PubMed]
  36. Rota, A.; Puddu, B.; Sabatini, C.; Panzani, D.; Lainé, A.L.; Camillo, F. Reproductive parameters of donkey jacks undergoing puberty. Anim. Reprod. Sci. 2018, 192, 119–125. [Google Scholar] [CrossRef]
  37. Saragusty, J.; Lemma, A.; Hildebrandt, T.B.; Göritz, F. Follicular size predicts success in artificial insemination with frozen-Thawed sperm in donkeys. PLoS ONE 2017, 12, 1–16. [Google Scholar] [CrossRef]
  38. Tosi, U.; Bernabò, N.; Verni, F.; Valbonetti, L.; Muttini, A.; Mattioli, M.; Barboni, B. Postpartum reproductive activities and gestation length in Martina Franca jennies, an endangered Italian donkey breed. Theriogenology 2013, 80, 120–124. [Google Scholar] [CrossRef]
  39. Quaresma, M.; Payan-Carreira, R. Characterization of the estrous cycle of Asinina de Miranda jennies (Equus asinus). Theriogenology 2015, 83, 616–624. [Google Scholar] [CrossRef]
  40. Fanelli, D.; Tesi, M.; Bagnato, G.; Salari, F.; Martini, M.; Panzani, D.; Camillo, F.; Rota, A. New simplified protocols for timed artificial insemination (TAI) in milk-producing donkeys. Theriogenology 2019, 139, 126–131. [Google Scholar] [CrossRef]
  41. Dai, F.; Segati, G.; Brscic, M.; Chincarini, M.; Costa, E.D.; Ferrari, L.; Burden, F.; Judge, A.; Minero, M. Effects of management practices on the welfare of dairy donkeys and risk factors associated with signs of hoof neglect. J. Dairy Res. 2018, 85, 30–38. [Google Scholar] [CrossRef] [Green Version]
  42. Dai, F.; Segati, G.; Dalla Costa, E.; Burden, F.; Judge, A.; Minero, M. Management practices and milk production in dairy donkey farms distributed over the Italian territory. Maced. Vet. Rev. 2017, 40, 131–136. [Google Scholar] [CrossRef] [Green Version]
  43. Deng, L.; Duan, H.; Zhang, X.; Zeng, S.; Wu, C.; Han, G. Advances in the research and application of artificial insemination to equids in China: 1935-2012. J. Equine Vet. Sci. 2014, 34, 351–359. [Google Scholar] [CrossRef]
  44. Meng, Q.-F.; Li, D.; Yao, G.-Z.; Zou, Y.; Cong, W.; Shan, X.-F. Seroprevalence of Toxoplasma gondii infection and variables associated with seropositivity in donkeys in eastern China. Parasite 2018, 25, 66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wang, H.; Liu, K.J.; Sun, Y.H.; Cui, L.Y.; Meng, X.; Jiang, G.M.; Zhao, F.W.; Li, J.J. Abortion in donkeys associated with Salmonella abortus equi infection. Equine Vet. J. 2019, 51, 756–759. [Google Scholar] [CrossRef] [PubMed]
  46. Turner, B.L.; Menendez, H.M.; Gates, R.; Tedeschi, L.O.; Atzori, A.S. System dynamics modeling for agricultural and natural resource management issues: Review of some past cases and forecasting future roles. Resources 2016, 5, 40. [Google Scholar] [CrossRef] [Green Version]
  47. Cong, W.; Nie, L.-B.; Qin, S.-Y.; Wang, W.-L.; Qian, A.-D.; Meng, Q.-F. Prevalence of Neospora spp. in donkeys in China. Parasite 2018, 25, 16. [Google Scholar] [CrossRef] [Green Version]
  48. Carluccio, A.; Gloria, A.; Veronesi, M.C.; De Amicis, I.; Noto, F.; Contri, A. Factors affecting pregnancy length and phases of parturition in Martina Franca jennies. Theriogenology 2015, 84, 650–655. [Google Scholar] [CrossRef]
  49. Carluccio, A.; Gloria, A.; Robbe, D.; Veronesi, M.C.; De Amicis, I.; Cairoli, F.; Contri, A. Reproductive characteristics of foal heat in female donkeys. Animal 2017, 11, 461–465. [Google Scholar] [CrossRef]
  50. Veronesi, M.C.; Villani, M.; Wilsher, S.; Contri, A.; Carluccio, A. A comparative stereological study of the term placenta in the donkey, pony and Thoroughbred. Theriogenology 2010, 74, 627–631. [Google Scholar] [CrossRef]
  51. Miragaya, M.H.; Neild, D.M.; Alonso, A.E. A Review of Reproductive Biology and Biotechnologies in Donkeys. J. Equine Vet. Sci. 2018, 65, 55–61. [Google Scholar] [CrossRef]
  52. De Oliveira, J.V.; de Oliveira, P.V.L.F.; Melo e Oña, C.M.; Guasti, P.N.; Monteiro, G.A.; da Sancler Silva, Y.F.R.; de Papa, P.M.; Alvarenga, M.A.; Dell’Aqua Junior, J.A.; Papa, F.O. Strategies to improve the fertility of fresh and frozen donkey semen. Theriogenology 2016, 85, 1267–1273. [Google Scholar] [CrossRef] [Green Version]
  53. Rota, A.; Panzani, D.; Sabatini, C.; Camillo, F. Donkey jack (Equus asinus) semen cryopreservation: Studies of seminal parameters, post breeding inflammatory response, and fertility in donkey jennies. Theriogenology 2012, 78, 1846–1854. [Google Scholar] [CrossRef] [PubMed]
  54. Chang, Z.L.; Li, B.X.; Liu, B.; Yao, L.; Yu, J.; Jiang, G.M.; Tan, J.H. Effects of FSH and the weather during induced ovulation and timed artificial insemination to increase jenny conception rates. Sci. Rep. 2019, 9, 1–9. [Google Scholar] [CrossRef] [PubMed]
  55. Wu, C. Status of Chinese donkey industry. In Proceedings of the 2017 1st International Symposium on Donkey Science, Dong’E, Shandong, China, 14–17 August 2017. [Google Scholar]
  56. Yang, F.; Wu, S.; Zhou, J.; Yang, W.; Liu, B.; Zhang, R.; Ma, J.; Yu, J.; Ji, C.; Zeng, S. The Proportion of Multiple-ovulation Changes Among Different Months in Dezhou Black Jennies. J. Equine Vet. Sci. 2018, 66, 155–156. [Google Scholar] [CrossRef]
  57. Xinhua News Agency. Shortage of Donkey Skin Breeds TCM Fakes. Available online: http://www.chinadaily.com.cn/china/2016-01/28/content_23287841.htm (accessed on 5 September 2017).
  58. Hessle, A.; Kumm, K.I.; Bertilsson, J.; Stenberg, B.; Sonesson, U. Combining environmentally and economically sustainable dairy and beef production in Sweden. Agric. Syst. 2017, 156, 105–114. [Google Scholar] [CrossRef]
  59. Domínguez, E.; Moreno-Irusta, A.; Castex, H.R.; Bragulat, A.F.; Ugaz, C.; Clemente, H.; Giojalas, L.; Losinno, L. Sperm Sexing Mediated by Magnetic Nanoparticles in Donkeys, a Preliminary In Vitro Study. J. Equine Vet. Sci. 2018, 65, 123–127. [Google Scholar] [CrossRef]
  60. Central Statistical Agency of the Federal Democratic Republic of Ethiopia. Report on Livestock and Livestock Characteristics, Agricultural Sample Survey 2017/18; Statistical Bulletin; Central Statistical Agency of the Federal Democratic Republic of Ethiopia: Addis Ababa, Ethiopia, 2018; Volume II.
  61. Yang, H.; Xiao, Y.; Meng, F.; Sun, F.; Chen, M.; Cheng, Z.; Chen, Y.; Liu, S.; Chen, H. Emergence of H3N8 equine influenza virus in donkeys in China in 2017. Vet. Microbiol. 2018, 214, 1–6. [Google Scholar] [CrossRef] [PubMed]
  62. Kumeta, Y.; Maruyama, T.; Asama, H.; Yamamoto, Y.; Hakamatsuka, T.; Goda, Y. Species identification of Asini Corii Collas (donkey glue) by PCR amplification of cytochrome b gene. J. Nat. Med. 2014, 68, 181–185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Wang, H.; Wang, K.; Wang, B.; Lv, Y.; Tao, X.; Zhang, X.; Tan, M. Integration of Artificial Neural Network Modeling and Hyperspectral Data Preprocessing for Discrimination of Colla Corii Asini Adulteration. J. Food Qual. 2018, 2018, 3487985. [Google Scholar] [CrossRef]
Figure 1. The price of donkey hides in China between 2000 and 2015 (in yuan). (Source: http://www.bestchinanews.com/Domestic/11652.html.).
Figure 1. The price of donkey hides in China between 2000 and 2015 (in yuan). (Source: http://www.bestchinanews.com/Domestic/11652.html.).
Animals 10 00718 g001
Figure 2. Stock and flow diagram of the reproductive performance model.
Figure 2. Stock and flow diagram of the reproductive performance model.
Animals 10 00718 g002
Figure 3. User interface of the donkey farming model.
Figure 3. User interface of the donkey farming model.
Animals 10 00718 g003
Figure 4. Number of years until the target hides of 2.0, 3.5 and 5.0 million per annum are reached under the three scenarios.
Figure 4. Number of years until the target hides of 2.0, 3.5 and 5.0 million per annum are reached under the three scenarios.
Animals 10 00718 g004
Figure 5. Total donkey herd required for the production of 2.0, 3.5 and 5.0 million hides per annum under the three scenarios.
Figure 5. Total donkey herd required for the production of 2.0, 3.5 and 5.0 million hides per annum under the three scenarios.
Animals 10 00718 g005
Table 1. Summary of the reproductive performance parameters and results.
Table 1. Summary of the reproductive performance parameters and results.
VariableLess Favourable ScenarioMiddle ScenarioMore Favourable Scenario
Conception rate0.60.70.8
Gestation length12.512.011.5
Abortion rate0.150.100.05
Twin percentage0510
Results
Reproductive performance0.830.951.09
Table 2. Percentage change in reproductive performance compared to middle scenario using less and more favourable values for individual variables.
Table 2. Percentage change in reproductive performance compared to middle scenario using less and more favourable values for individual variables.
VariableLess Favourable Value% Change in Reproductive PerformanceMore Favourable Value% Change in Reproductive Performance
Conception rate0.6−1.9%0.81.6%
Gestation length12.5−3.8%11.54.2%
Abortion rate0.15−3.3%0.053.2%
Twin percentage0−4.8%104.8%
Table 3. Summary of farming model parameter values and results for hide production of 2.0, 3.5 and 5.0 million hides.
Table 3. Summary of farming model parameter values and results for hide production of 2.0, 3.5 and 5.0 million hides.
VariableLess FavourableMiddle ScenarioMore Favourable
Initial female breeding herd50,000100,000200,000
Reproductive performance0.830.951.09
Reproduction age36 months30 months24 months
Average breeding life (annual cull rate)10 years (10.0%)12 years (8.3%)15 years (6.7%)
Percentage female births50 percent55 percent65 percent
Slaughter age24 months22 months17 months
Mortality—foals10 percent7 percent5 percent
Mortality—breeding animals (annual)4 percent3 percent2 percent
Jennie to jack ratio102030
Results for 2.0 million hides
Hides target reached years38 years20 years11 years
Breeding jennies2.9 million2.4 million2.1 million
Donkey herd8.1 million6.5 million5.8 million
Results for 3.5 million hides
Hides target reached years43 years23 years12 years
Breeding jennies5.1 million4.1 million3.6 million
Donkey herd14.1 million11.4 million10.2 million
Results for 5.0 million hides
Hides target reached years46 years25 years14 years
Breeding jennies7.3 million5.9 million5.1 million
Donkey herd20.2 million16.4 million14.6 million
Table 4. Percentage change in time to reach target hides, total donkeys and breeding jennies compared to middle scenario using less and more favourable values for individual variables.
Table 4. Percentage change in time to reach target hides, total donkeys and breeding jennies compared to middle scenario using less and more favourable values for individual variables.
Variables (LF, M, MF a)% Change in Time to Target Hides% Change in Total Donkeys% Change in Breeding Jennies
LFMFLFMFLFMF
Target hides
(5, 3.5, 2 million)
9−1343−4343−43
Initial female breeding herd
(50, 100, 200 thousand)
16−160000
Reproductive performance
(0.83, 0.95, 1.09)
15−127−615−13
Reproduction age
(36, 30, 24 months)
9−103−111
Average breeding life
(10, 12, 15 years)
8−41−100
Percentage female births
(50, 55, 65)
8−110−100
Slaughter age
(24, 22, 17 months)
1−12−600
Mortality—foals
(10, 7, 5 percent)
4−22−23−2
Mortality—breeding animals (annual)
(4, 3, 2 percent)
3−32−12−1
Jennie to jack ratio
(10,20,30)
004−120
a LF = less favourable, M = middle, and MF = more favourable.

Share and Cite

MDPI and ACS Style

Bennett, R.; Pfuderer, S. The Potential for New Donkey Farming Systems to Supply the Growing Demand for Hides. Animals 2020, 10, 718. https://doi.org/10.3390/ani10040718

AMA Style

Bennett R, Pfuderer S. The Potential for New Donkey Farming Systems to Supply the Growing Demand for Hides. Animals. 2020; 10(4):718. https://doi.org/10.3390/ani10040718

Chicago/Turabian Style

Bennett, Richard, and Simone Pfuderer. 2020. "The Potential for New Donkey Farming Systems to Supply the Growing Demand for Hides" Animals 10, no. 4: 718. https://doi.org/10.3390/ani10040718

APA Style

Bennett, R., & Pfuderer, S. (2020). The Potential for New Donkey Farming Systems to Supply the Growing Demand for Hides. Animals, 10(4), 718. https://doi.org/10.3390/ani10040718

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