A comparison of multiple indicator kriging and area-to-point Poisson kriging for mapping patterns of herbivore species abundance in Kruger National Park, South Africa

Kruger National Park (KNP), South Africa, provides protected habitats for the unique animals of the African savannah. For the past 40 years, annual aerial surveys of herbivores have been conducted to aid management decisions based on (1) the spatial distribution of species throughout the park and (2) total species populations in a year. The surveys are extremely time consuming and costly. For many years, the whole park was surveyed, but in 1998 a transect survey approach was adopted. This is cheaper and less time consuming but leaves gaps in the data spatially. Also the distance method currently employed by the park only gives estimates of total species populations but not their spatial distribution. We compare the ability of multiple indicator kriging and area-to-point Poisson kriging to accurately map species distribution in the park. A leave-one-out cross-validation approach indicates that multiple indicator kriging makes poor estimates of the number of animals, particularly the few large counts, as the indicator variograms for such high thresholds are pure nugget. Poisson kriging was applied to the prediction of two types of abundance data: spatial density and proportion of a given species. Both Poisson approaches had standardized mean absolute errors (St. MAEs) of animal counts at least an order of magnitude lower than multiple indicator kriging. The spatial density, Poisson approach (1), gave the lowest St. MAEs for the most abundant species and the proportion, Poisson approach (2), did for the least abundant species. Incorporating environmental data into Poisson approach (2) further reduced St. MAEs.

[1]  P. Monestieza,et al.  Geostatistical modelling of spatial distribution of Balaenoptera physalus in the Northwestern Mediterranean Sea from sparse count data and heterogeneous observation efforts , 2006 .

[2]  Sassan Saatchi,et al.  Woody Fractional Cover in Kruger National Park, South Africa: Remote Sensing–Based Maps and Ecological Insights , 2010 .

[3]  P. Kyriakidis.  A Geostatistical Framework for Area-to-Point Spatial Interpolation , 2004 .

[4]  Pierre Goovaerts,et al.  Comparative performance of indicator algorithms for modeling conditional probability distribution functions , 1994 .

[5]  I. Smit.  Resources driving landscape‐scale distribution patterns of grazers in an African savanna , 2011 .

[6]  Jaap van der Meer,et al.  Assessing the population size of the European storm petrel (Hydrobates pelagicus) using spatial autocorrelation between counts from segments of criss-cross ship transects , 1995 .

[7]  Pierre Goovaerts,et al.  Applying Geostatistical Analysis to Crime Data: Car-Related Thefts in the Baltic States. , 2010, Geographical analysis.

[8]  W. Getz,et al.  Biases in estimating population size from an aerial census: a case study in the Kruger National Park, South Africa , 2002 .

[9]  C. C. Grant,et al.  Managing surface-water in a large semi-arid savanna park: Effects on grazer distribution patterns , 2009 .

[10]  R. Kerry,et al.  Determining the effect of asymmetric data on the variogram. I. Underlying asymmetry , 2007, Comput. Geosci..

[11]  Chen-Wuing Liu,et al.  Evaluation of potential health risk of arsenic-affected groundwater using indicator kriging and dose response model. , 2007, The Science of the total environment.

[12]  Hulin Wu,et al.  Modelling the distribution of plant species using the autologistic regression model , 1997, Environmental and Ecological Statistics.

[13]  J.-P. Dubois,et al.  Assessing the risk of soil contamination in the Swiss Jura using indicator geostatistics , 1997, Environmental and Ecological Statistics.

[14]  Donald G. Bullock,et al.  A comparative study of interpolation methods for mapping soil properties , 1999 .

[15]  R. M. Lark,et al.  A comparison of some robust estimators of the variogram for use in soil survey , 2000 .

[16]  P. Beja,et al.  The use of sighting data to analyse Iberian lynx habitat and distribution , 1999 .

[17]  Pierre Goovaerts,et al.  Evaluating the probability of exceeding a site-specific soil cadmium contamination threshold , 2001 .

[18]  David J. Mulla,et al.  Geostatistical Tools for Modeling and Interpreting Ecological Spatial Dependence , 1992 .

[19]  Yu-Pin Lin,et al.  Factorial and indicator kriging methods using a geographic information system to delineate spatial variation and pollution sources of soil heavy metals , 2002 .

[20]  R. M. Lark,et al.  Two robust estimators of the cross‐variogram for multivariate geostatistical analysis of soil properties , 2003 .

[21]  B. Reilly,et al.  Application of distance sampling to estimate population densities of large herbivores in Kruger National Park , 2008 .

[22]  Marc Voltz,et al.  A comparison of kriging, cubic splines and classification for predicting soil properties from sample information , 1990 .

[23]  W. Gertenbach,et al.  Landscapes of the Kruger National Park , 1983 .

[24]  R. Kerry,et al.  Determining the effect of asymmetric data on the variogram. II. Outliers , 2007, Comput. Geosci..

[25]  N. Owen‐Smith,et al.  Manifold interactive influences on the population dynamics of a multispecies ungulate assemblage , 2006 .

[26]  S. McKenna,et al.  Delineating high-density areas in spatial Poisson fields from strip-transect sampling using indicator geostatistics: application to unexploded ordnance removal. , 2007, Journal of environmental management.

[27]  Stephen T. Buckland,et al.  Distance 4.1. Release 2: Research unit for wildlife population assessment [computer program] , 2004 .

[28]  Chung-Min Liao,et al.  Evaluation of arsenic contamination potential using indicator kriging in the Yun-Lin aquifer (Taiwan). , 2004, The Science of the total environment.

[29]  C. Gotway,et al.  Comparison of kriging and inverse-distance methods for mapping soil parameters , 1996 .

[30]  Pierre Goovaerts,et al.  Kriging and Semivariogram Deconvolution in the Presence of Irregular Geographical Units , 2008, Mathematical geology.

[31]  C. Cosner,et al.  Early detection of declining populations using floor and ceiling models , 2001 .

[32]  J. Ogutu,et al.  Efficiency of strip- and line-transect surveys of African savanna mammals , 2006 .

[33]  S. Buckland.  Introduction to distance sampling : estimating abundance of biological populations , 2001 .

[34]  Pierre Goovaerts,et al.  AUTO-IK: A 2D indicator kriging program for the automated non-parametric modeling of local uncertainty in earth sciences , 2009, Comput. Geosci..

[35]  Pierre Goovaerts,et al.  International Journal of Health Geographics Geostatistical Analysis of Disease Data: Estimation of Cancer Mortality Risk from Empirical Frequencies Using Poisson Kriging , 2022 .

[36]  Sam M. Ferreira,et al.  Management intervention affects river-bound spatial dynamics of elephants , 2010 .

[37]  N. Owen‐Smith,et al.  Establishing the causes of the roan antelope decline in the Kruger National Park, South Africa , 1999 .

[38]  P. Monestiez,et al.  Geostatistical Modelling of Wildlife Populations: A Non-stationary Hierarchical Model for Count Data , 2010 .

[39]  R. M. Lark,et al.  Robust estimation of the pseudo cross‐variogram for cokriging soil properties , 2002 .

[40]  J. Ogutu,et al.  ENSO, rainfall and temperature influences on extreme population declines among African savanna ungulates , 2003 .

[41]  Margaret A. Oliver,et al.  Variograms of Ancillary Data to Aid Sampling for Soil Surveys , 2003, Precision Agriculture.

[42]  Pierre Goovaerts,et al.  Accounting for measurement error in uncertainty modeling and decision‐making using indicator kriging and p‐field simulation: application to a dioxin contaminated site , 2002 .