International Journal of Ecology

International Journal of Ecology / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 491858 |

L. L. Eberhardt, J. M. Breiwick, "A Scheme for Evaluating Feral Horse Management Strategies", International Journal of Ecology, vol. 2012, Article ID 491858, 5 pages, 2012.

A Scheme for Evaluating Feral Horse Management Strategies

Academic Editor: Bradford Hawkins
Received26 Dec 2011
Revised26 Mar 2012
Accepted27 Mar 2012
Published29 May 2012


Context. Feral horses are an increasing problem in many countries and are popular with the public, making management difficult. Aims. To develop a scheme useful in planning management strategies. Methods. A model is developed and applied to four different feral horse herds, three of which have been quite accurately counted over the years. Key Results. The selected model has been tested on a variety of data sets, with emphasis on the four sets of feral horse data. An alternative, nonparametric model is used to check the selected parametric approach. Conclusions. A density-dependent response was observed in all 4 herds, even though only 8 observations were available in each case. Consistency in the model fits suggests that small starting herds can be used to test various management techniques. Implications. Management methods can be tested on actual, confined populations.

1. Introduction

Feral horses continue to be a problem in the Western United States [1, 2] and elsewhere [3]. Some recent studies provide the basis for a simple model that should be useful in developing management approaches, and these are assessed here, along with the basis for selecting a model.

We have located 4 data sets on feral horses that may serve as a basis for designing tests of management schemes. These data come from widely dispersed locations (Argentina, France, and the United States). The studies were conducted on horse herds confined to fenced areas. Two of the studies were conducted in relatively small areas, making replication quite feasible. The other two were based on sizable populations and would thus be expensive to replicate.

The Argentine study was conducted on an isolated (but sizable) herd that was enumerated very carefully. All 4 studies gave clear evidence of density dependence in a relatively short period (8 years). Our proposal is to emulate these studies (in fenced areas) to test various management schemes, such as reproductive inhibition and the removal of the young females that are readily adopted. A successful approach would be evidenced by a substantial reduction in population growth. Presumably control measures would be applied while the population is well below the inflection point in the growth curve and that the test populations could then be maintained long enough to ascertain the efforts required for continued control.

2. Material and Methods

2.1. Population Growth Curves

Perhaps the most popular model for population growth is the generalized logistic of Nelder [14], Pella and Tomlinson [15], best-known from Ayala et al. [16], and often termed the theta logistic. The model is The corresponding difference equation is Recently, this model has been shown to have some shortcomings [17, 18], so we here use a simplification in which the parameter controlling the inflection point () is set to 2.0: The basis for this result appears in Eberhardt and Breiwick [19]. The corresponding difference equation is (2) with. An alternative, nonparametric model is It can readily be fit as a linear model and used to derive confidence limits, following Draper and Smith [20].

If we set in (1), the well-known ordinary logistic curve results, tracing back to Verhulst in 1844: The proportion of the asymptotic population size present at the inflection point can be defined as [21]: With , this is . This result is to be expected in that it has been known for some time that the inflection point in population growth curves often is higher than the 50% level of the ordinary logistic curve [2225]. denotes the asymptotic value, and is the rate of increase, governs the inflection point of the curve, and denotes population size. The curves were fit with nonlinear least-squares [26], as implemented in the R-language [27].

Part of the data used here (Table 1) was obtained from the Global Population Dynamics Database (GPDD) maintained by the National Environmental Research Council at Sources for some additional sets and the corresponding references are listed in Table 1.

SpecieszStandard errorCoef. of Variation

Red deer 65853.891.060.272
Boar 94284.995.681.138
Gray fox 3193.781.540.407
Gray fox 3421.852.131.151
Coyote 2631.932.861.482
Gray seals Bonner (1975) [4]11.7610.370.882
Bison Fredin (1984) [5]0.750.470.627
Macaque 606816.2625.121.545
Beaver 1510.460.81.739
Beaver 1567.437.631.027
Elephant seals Stewart et al. (1994) [6]1.860.670.36
Seneca deer Hesselton et al. (1965) [7]0.690.450.652
Badger 781.764.082.29
Fur seals Kenyon et al. (1954) [8]
Blue Springs manatees Eberhardt and O’Shea (1995) [9]5.0810.712.108
Chinstrap penguins Fraser (pers. commun.)0.452.65.778
Gentoo penguins Fraser (pers. commun.)0.691.442.087
Polar bear 1436.136.911.127
Beaver 1567.437.631.02
Coyote 2476.4210.741.673
Gray fox 3281.882.431.293
Boar 94291.660.990.596
Marten 96951.22.271.892
Squirrel 96990.982.052.92
Gray fox 3340.453.482.48

2.2. Biases in Estimation

There is a large recent literature concerned with various error terms in population models ([18, 2830] provide examples). Much of the literature on biases is concerned with decreasing populations, rather than with the growth curves considered here, and it seems that there is little agreement yet on terms and methodology. We believe that the dominant sources are “observation error” resulting from inaccuracies in counting and “process error” which may result from various sources, such as fluctuations due to changing age structure as populations increase. Because our major proposal is concerned with devising methods to stop population growth, biases of the kind described above are largely neglected here. We have examined growth curves for some populations that have been counted very accurately, usually by corralling the populations, so that “observational error” is absent or negligible, and the fluctuations present are largely due to “process error.” Very likely, the two sources are independent (“observational error” can scarcely be influenced by “process error” unless the latter is very large indeed). Hence, we believe that it is reasonable to assume that the two sources are additive and might be treated, at least approximately, by assuming a single-error term in curve fitting. We thus assume that the calculations underlying Table 1 provide useful estimates of the parameter , governing the inflection point.

3. Results

Figure 1 shows fits of the modified logistic curve (3) to the four wild horse data sets in the left-hand column. Confidence limits (95%) for (4) are shown by crosses. The right-hand column contains the corresponding fits from the “Sibly equation” [3134] resulting from a rearrangement of (2), again with : Table 1 provides data on estimates of the parameter in (1) from 25 data sets and shows that the median value is close to the value (2) used in (3).

Table 2 gives the parameters from fitting (3) to the four data sets on wild horses. All four have fairly large values of , compared to many other birds and mammals, and relatively small standard errors, as might be expected due to the accurate counts. Two populations (Argentine and Beaty horses) started out with sizable initial populations. All four populations grew beyond the inflection point within 8 years. We fitted the generalized logistic curve (1), the logistic curve (5), and the Sibly curve (7) to the data and checked to see if the fitted curves passed through the confidence limits shown in Figure 1. The Beaty horse data were too variable for satisfactory convergence, and the Sibly model yielded values above the upper confidence limit for the Camargue (87) and Prezewalski (52) herd data. The remaining curves fell within the confidence interval.

HerdScientific nameKr Standard error of K Source

CamargueEquis cabella940.326101.08Grange et al. 2009 [10]
PrezewalskiEquis ferus530.26291.08Tatin et al. 2009 [11]
ArgentineEquis cabella2300.216915.12Scorolli and Lopez-Cazarola 2010 [12]
BeatyEquis cabella12020.2813628.23Eberhardt et al. 1982 [13]

4. Discussion

As indicated in the Introduction above, our purpose is to identify a simple model for the growth of wild horse populations that should be useful in evaluating management strategies. In spite of the small sample sizes (8 observations in each case), it appears that the curves provide useful fits to the data with moderate coefficients of variation for estimates of the asymptotic values (4% for the Camargue horses, 8% for the Prezewalski horses, 5% for the Argentine horses, and 17% for the Beaty horse example). The high CV for Beaty horses likely is a consequence of the fact that the maximal count fell some distance below the estimated asymptote. Rates of increase (Table 2) were in reasonable agreement with a value obtained from reproductive and survival data on feral horses [35]. It should be noted, however, that the rate of increase given by Eberhardt [35] is based on assuming an exponential model, which is not as suitable for feral horses as is (3) which accommodates the decline in growth rate as the population increases.

5. Management Implications

We have identified 4 cases (on 4 quite different locations) where feral horse data is usefully described by a simple model, the modified logistic curve of (3). All of the populations were either in areas remote from other feral horse populations (Argentine) or enclosed in fences and provided useful estimates of the asymptotic value and of the rates of increase. Actual applications for management purposes will depend on circumstances and the management method to be tested. We believe that tests should be conducted in fenced areas, with a known initial population. Very likely, the initial populations should be below those observed for the Argentine and Beaty herds. As noted in the Introduction, if control is achieved, the resulting population should be maintained for some time to determine how long a given method will maintain control and what efforts are required to maintain control. Costs of fencing and maintaining a site may be significant, and it may be possible to do some tests in areas isolated from other herds.


The authors acknowledge support from the National Marine Fisheries Agency and the US Marine Mammal Commission.


  1. B. Abbey, “Managing wild horses on public lands,” Wildlife Professional:10, 2010. View at: Google Scholar
  2. J. Jeffries and P. Rouch, “Lethal hoofbeats,” Wildlife Professional: 50-55, 2010. View at: Google Scholar
  3. M. Thomas, “Is shooting the answer downunder?” Wildlife Professional 57, 2010. View at: Google Scholar
  4. W. N. Bonner, “Population increase of grey seals at the Farne Islands. Biology of the seal,” in Rapports et Proces-Verbaux des Reunions, K. Ronald and A. W. Mansfield, Eds., vol. 169, pp. 366–370, Conseil International pour l’Exploration de la Mer, 1975. View at: Google Scholar
  5. R. A. Fredin, “Levels of maximum net productivity in populations of large terrestrial mammals. Reproduction in whales, dolphins, and porpoises,” in Reports of the International Whaling Commission, W. F. Perrin, R. L. Brownell, and D. P. DeMaster, Eds., no. 6, pp. 381–387, International Whaling Commission, Cambridge, UK, 1984. View at: Google Scholar
  6. B. S. Stewart, P. K. Yoachem, H. R. Huber et al., “History and present status of the northern elephant seal population,” in Elephant Seals, Population Ecology, Behavior, and Physiology, B. J. LeBouef and R. M. Laws, Eds., pp. 29–48, University of California Press, Berkeley, Calif, USA, 1994. View at: Google Scholar
  7. W. T. Hesselton, C. W. Sevringhaus, and J. E. Tanck, “Population dynamics of deer in the Seneca army depot,” New York Fish and Game Journal, vol. 12, pp. 17–30, 1965. View at: Google Scholar
  8. K. W. Kenyon, V. H. Scheffer, and D. G. Chapman, “A Population Study of the Alaska Fur Seal Herd,” Special Scientific Report, Wildlife No. 12, United States Department of the Interior, U. S. Fish and Wildlife Service, Washington, DC, USA, 1954. View at: Google Scholar
  9. L. L. Eberhardt and T. J. O'Shea, “Integrating manatee life history data and population model development,” in Population Biology of the Florida Manatee, T. J. O'Shea, B. B. Ackerman, and H. F. Percival, Eds., National Biological Service Technical Report Series, pp. 269–278, Information and Technology Report 1, Fort Collins, Colo, USA, 1995. View at: Google Scholar
  10. S. Grange, P. Duncan, and J. M. Gaillard, “Poor horse traders: large mammals trade survival for reproduction during the process of feralization,” Proceedings of the Royal Society B, vol. 276, no. 1663, pp. 1911–1919, 2009. View at: Publisher Site | Google Scholar
  11. L. Tatin, S. R. B. King, B. Munkhtuya, A. J. M. Hewison, and C. Feh, “Demography of a socially natural herd of Przewalski's horses: an example of a small, closed population,” Journal of Zoology, vol. 277, no. 2, pp. 134–140, 2009. View at: Publisher Site | Google Scholar
  12. Scorolli and A. L. A. C. Lopez-Cazarola, “Demography of feral horses (Equus caballus): a long-term study in Tompquist Park, Argentina,” Wildlife Research, vol. 37, no. 3, pp. 207–214, 2010. View at: Publisher Site | Google Scholar
  13. L. L. Eberhardt, A. K. Majorowicz, and J. A. Wilcox, “Apparent rates of increase for two feral horse herds,” The Journal of Wildlife Management, vol. 46, pp. 367–374, 1982. View at: Google Scholar
  14. J. A. Nelder, “The fitting of a generalization of the logistic curve,” Biometrics, vol. 17, pp. 89–110, 1961. View at: Google Scholar
  15. J. J. Pella and P. K. Tomlinson, “A generalized stock production model,” Bulletin of the Inter-American Tropical Tuna Commission, vol. 13, pp. 419–496, 1969. View at: Google Scholar
  16. F. J. Ayala, M. E. Gilpin, and J. G. Ehrenfeld, “Competition between species: theoretical models and experimental tests,” Theoretical Population Biology, vol. 4, no. 3, pp. 331–356, 1973. View at: Google Scholar
  17. L. Polansky, P. de Valpine, J. O. Lloyd-Smith, and W. M. Getz, “Likelihood ridges and multimodality in population growth rate models,” Ecology, vol. 90, no. 8, pp. 2313–2320, 2009. View at: Publisher Site | Google Scholar
  18. F. Clark, B. W. Brook, S. Delean, H. R. Akcakaya, and C. J. A. Bradshaw, “The theta-logistic is unreliable for modeling most census data,” Methods in Ecology and Evolution, vol. 1, pp. 253–262, 2010. View at: Google Scholar
  19. L. L. Eberhardt and J. M. Breiwick, “2012. Models for population growth curves,” ISRN Ecology, vol. 2012, Article ID 815016, 7 pages, 2012. View at: Publisher Site | Google Scholar
  20. N. R. Draper and H. Smith, Applied Regression Analysis, Wiley, 3rd edition, 1998.
  21. L. L. Eberhardt, “Population projections from simple models,” Journal of Applied Ecology, vol. 24, no. 1, pp. 103–118, 1987. View at: Google Scholar
  22. D. S. Butterworth and P. B. Best, “The origins and choice of 54% carrying capacity as the protection level for baleen whale stocks, and the implications thereof for management procedures,” Reports International Whaling Commission, vol. 44, pp. 491–497, 1994. View at: Google Scholar
  23. L. L. Eberhardt, “Optimal policies for conservation of large mammals, with special reference to marine ecosystems,” Environmental Conservation, vol. 4, no. 3, pp. 205–212, 1977. View at: Google Scholar
  24. L. L. Eberhardt and D. B. Siniff, “Population dynamics and marine mammal management policies,” Journal of the Fisheries Research Board of Canada, vol. 34, pp. 183–190, 1977. View at: Google Scholar
  25. C. W. Fowler, “Density dependence as related to life history strategy,” Ecology, vol. 62, pp. 602–610, 1981. View at: Google Scholar
  26. D. M. Bates and D. G. Watts, Nonlinear Regression Analysis and Its Applications, John Wiley & Sons, New York, NY, USA, 1988.
  27. R Development Core Team, “R: a language and environment for statistical computing,” R Foundation for Statistical Computing, Vienna, Austria, 2004, View at: Google Scholar
  28. J. P. Buonaccorsi and J. Staudenmayer, “Statistical methods to correct for observation error in a density-independent population model,” Ecological Monographs, vol. 79, no. 2, pp. 299–324, 2009. View at: Publisher Site | Google Scholar
  29. E. E. Holmes, J. L. Sabo, S. V. Viscido, and W. F. Fagan, “A statistical approach to quasi-extinction forecasting,” Ecology Letters, vol. 10, no. 12, pp. 1182–1198, 2007. View at: Publisher Site | Google Scholar
  30. D. F. Staples, M. L. Taper, and B. Dennis, “Estimating population trend and process variation for PVA in the presence of sampling error,” Ecology, vol. 85, no. 4, pp. 923–929, 2004. View at: Google Scholar
  31. R. M. Sibly, D. Barker, M. C. Denham, J. Hone, and M. Pagel, “Ecology: on the regulation of populations of mammals, birds, fish, and insects,” Science, vol. 309, no. 5734, pp. 607–610, 2005. View at: Publisher Site | Google Scholar
  32. B. W. Brook and C. J. A. Bradshaw, “Strength of evidence for density dependence in abundance time series of 1198 species,” Ecology, vol. 87, no. 6, pp. 1445–1451, 2006. View at: Google Scholar
  33. B. E. Sæther, V. Grotan, P. Tryjanowski, C. Barbraud, S. Engen, and M. Fulin, “Climate and spatio-temporal variation in the population dynamics of a long distance migrant, the white stork,” Journal of Animal Ecology, vol. 75, no. 1, pp. 80–90, 2006. View at: Publisher Site | Google Scholar
  34. A. Berryman and M. Lima, “Detecting the order of population dynamics from time series: nonlinearity causes spurious diagnosis,” Ecology, vol. 88, no. 8, pp. 2121–2123, 2007. View at: Publisher Site | Google Scholar
  35. L. L. Eberhardt, “A paradigm for population analysis of long-lived vertebrates,” Ecology, vol. 83, no. 10, pp. 2841–2854, 2002. View at: Google Scholar

Copyright © 2012 L. L. Eberhardt and J. M. Breiwick. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles