Stable Isotope Sourcing Using Sampling


Stable isotope mixing models are used to estimate proportional contributions of sources to a mixture, such as in the analysis of animal diets, plant nutrient use, geochemistry, pollution, and forensics. We describe an algorithm implemented as SISUS software for providing a user-specified number of probabilistic exact solutions derived quickly from the extended mixing model. Our method outperforms IsoSource [1], a deterministic algorithm for providing approximate solutions to represent the solution polytope. Our method is an approximate Bayesian large sample procedure. SISUS software is freely available at and as an R package at

Share and Cite:

Erhardt, E. , Wolf, B. , Ben-David, M. and Bedrick, E. (2014) Stable Isotope Sourcing Using Sampling. Open Journal of Ecology, 4, 289-298. doi: 10.4236/oje.2014.46027.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Phillips, D.L. and Gregg, J.W. (2003) Source Partitioning Using Stable Isotopes: Coping with Too Many Sources. Oecologia, 136, 261-269.
[2] Phillips, D.L. (2001) Mixing Models in Analyses of Diet Using Multiple Stable Isotopes: A Critique. Oecologia, 127, 166-170.
[3] Martinez del Rio, C. and Wolf, B.O. (2005) Mass-Balance Models for Animal Isotopic Ecology. In: Matthias Starck, J., and Wang, T., Eds., Physiological and Ecological Adaptations to Feeding in Vertebrates, Science Publishers, Inc., Enfield, Chapter 6, 141-174.
[4] Hobson, K.A. and Wassenaar, L.I. (Eds.) (2008) Tracking Animal Migration with Stable Isotopes. Academic Press.
[5] Ehleringer, J.R., Bowen, G.J., Chesson, L.A., West, A.G., Podlesak, D.W. and Cerling, T.E. (2008) Hydrogen and Oxygen Isotope Ratios in Human Hair Are Related to Geography. Proceedings of the National Academy of Sciences, 105, 2788.
[6] Kendall, C. and McDonnell, J.J. (Eds.) (1998) Isotope Tracers in Catchment Hydrology. Elsevier, Amsterdam.
[7] DeNiro, M. and Epstein, S. (1978) Influence of the Diet on the Distribution of Carbon Isotopes in Animals. Geochimicaet Cosmochimica Acta, 42, 495-506.
[8] Schwarcz, H.P. (1991) Some Theoretical Aspects of Isotope Paleodiet Studies. Journal of Archaeological Science, 18, 261-275.
[9] DeNiro, M. and Epstein, S. (1981) Influence of the Diet on the Distribution of Nitrogen Isotopes in Animals. Geochimicaet Cosmochimica Acta, 48, 341-351.
[10] Minagawa, M. and Wada, E. (1984) Stepwise Enrichment of 15N along Food Chains: Further Evidence and the Relation between N and Animal Age. Geochimicaet Cosmochimica Acta, 48, 1135-1140.
[11] Koch, P.L. and Phillips, D.L. (2002) Incorporating Concentration Dependence in Stable Isotope Mixing Models: A Reply to Robbins, Hilderbrand and Farley. Oecologia, 133, 14-18.
[12] Jacoby, M.E., Hilderbrand, G.V., Servheen, C., Schwartz, C.C., Arthor, S.M., Hanley, T.A., Robbins, C.T. and Michener, R. (1999) Trophic Relations of Brown and Black Bears in Several Western North American Ecosystems. Journal of Wildlife Management, 63, 921-929.
[13] Erhardt, E.B. and Bedrick, E.J. (2014) Inference for Stable Isotope Mixing Models: A Study of the Diet of Dunlin. Journal of the Royal Statistical Society: Series C, online.
[14] Erhardt, E.B. and Bedrick, E.J. (2013) A Bayesian Framework for Stable Isotope Mixing Models. Environmental and Ecological Statistics, 20, 377-397.
[15] Phillips, D.L. (2006) IsoSource (Version 1.3.1).
[16] Fukuda, K. and Prodon, A. (1996) Double Description Method Revisited. Lecture Notes in Computer Science, 1120, 91-111.
[17] Fukuda, K. (2005) C-Library Cddlib (Version 0.94f). home/cdd.html
[18] Boneh, A. and Golan, A. (1979) Constraints’ Redundancy and Feasible Region Boundedness by Random Feasible Point Generator (rfpg). 3rd European Congress on Operations Research (EURO III), Amsterdam, 9-11 April 1979.
[19] Smith, R.L. (1980) Monte Carlo Procedures for Generating Random Feasible Solutions to Mathematical Programs. ORSA/TIMS Conference, Washington DC.
[20] Smith, R.L. (1984) Efficient Monte Carlo Procedures for Generating Points Uniformly Distributed over Bounded Regions. Operations Research, 32, 1296-1308.
[21] R Development Core Team (2006) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna.
[22] Lazar, R. and Meeden, G. (2003) Isipta’03, Proceedings of the 3rd International Symposium on Imprecise Probabilities and Their Applications. In: Bernard, J.M., Seidenfeld, T. and Zaffalon, M., Eds., ISIPTA, Proceedings in Informatics, Vol. 18, 361-371.
[23] Meeden, G. and Lazar, R. (2006) Polyapost: Simulating from the Polya Posterior. R Package Version 1.1.
[24] Gelman, A., Carlin, J.B., Stern, H.S. and Rubin, D.B. (1995) Bayesian Data Analysis. Chapman and Hall, London.
[25] Ben-David, M., Hanley, T.A., Klein, D.R. and Schell, D.M. (1997) Seasonal Changes in Diets of Coastal and Riverine Mink: The Role of Spawning Pacific Salmon. Canadian Journal of Zoology, 75, 803-811.
[26] Moore, J.W. and Semmens, B.X. (2008) Incorporating Uncertainty and Prior Information into Stable Isotope Mixing Models. Ecology Letters, 11, 470-480.
[27] Parnell, A. and Jackson, A. (2008) Siar: Stable Isotope Analysis in R. R Package Version 3.3.
[28] Parnell, A.C., Inger, R., Bearhop, S. and Jackson, A.L. (2010) Source Partitioning Using Stable Isotopes: Coping with Too Much Variation. PLoS ONE, 5, Article ID: e9672.

Copyright © 2023 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.