Mapping population change index in Southern Serbia (1961-2027) as a function of environmental factors

Niche analysis methods developed within the biogeography community are routinely used for species distribution modeling of wildlife and endangered species. So far, such techniques have not been used to explain distribution of people in an area, nor to assess spatio-temporal dynamics of human populations. In this paper, the MaxEnt approach to species distribution modeling and publicly available gridded predictors were used to analyze the population dynamics in Southern Serbia (South Pomoravlje Region) for the period 1961-2027. Population values from the census administrative units were first downscaled to 200 m grid using a detailed map of populated places and dasymetric interpolation. In the second step, a point pattern representing the whole population (468,500 inhabitants in 2002) was simulated using the R package spatstat. MaxEnt was then used to derive habitat suitability index (HSI) as a function of gridded predictors: distance to roads, elevation, slope, topographic wetness index, enhanced vegetation index and land cover classes. HSI and environmental predictors were further used to explain spatial patterns in the population change index (PCI) through regression modeling. The results show that inhabiting preference for year 1961 is mainly a function of topography (TWI, elevation). The HSI for year 2027 shows that large portions of remote areas are becoming less preferred for inhabiting. The results of cross-validation in MaxEnt show that distribution of population is distinctly controlled by environmental factors (AUC > 0.84). Population decrease is particularly significant in areas >25 km distant from the main road network. The results of regression analysis show that 40% of variability in the PCI values can be explained with these environmental maps, distance to roads and urban areas being the main drivers of migration process. This approach allows precise mapping of demographic patterns that otherwise would not be visible from the census data alone.

[1]  S. Milijić,et al.  Problems of spatial-functional organization of Južno Pomoravlje region's network of settlements , 2009 .

[2]  E. Farmer,et al.  A primer of GIS: fundamental geographic and cartographic concepts , 2009 .

[3]  Miroslav Dudík,et al.  Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation , 2008 .

[4]  Olaf Conrad,et al.  SAGA - Entwurf, Funktionsumfang und Anwendung eines Systems für Automatisierte Geowissenschaftliche Analysen , 2007 .

[5]  Michael N. DeMers,et al.  Fundamentals of Geographic Information Systems , 1996 .

[6]  A. Townsend Peterson,et al.  Novel methods improve prediction of species' distributions from occurrence data , 2006 .

[7]  John Bell,et al.  A review of methods for the assessment of prediction errors in conservation presence/absence models , 1997, Environmental Conservation.

[8]  J. E. Dobson,et al.  LandScan: A Global Population Database for Estimating Populations at Risk , 2000 .

[9]  E. Jaynes.  Information Theory and Statistical Mechanics , 1957 .

[10]  Mahbub ul Haq,et al.  Reflections on human development , 1995 .

[11]  Guangqing Chi,et al.  Applied Spatial Data Analysis with R , 2015 .

[12]  C. E. SHANNON,et al.  A mathematical theory of communication , 1948, MOCO.

[13]  Michael E. Bufalino,et al.  Street-Weighted Interpolation Techniques for Demographic Count Estimation in Incompatible Zone Systems , 2005 .

[14]  Robert P. Anderson,et al.  Maximum entropy modeling of species geographic distributions , 2006 .

[15]  Mitchel Langford,et al.  Obtaining population estimates in non-census reporting zones: An evaluation of the 3-class dasymetric method , 2006, Comput. Environ. Urban Syst..

[16]  Christopher R. Herrmann,et al.  Mapping Population Distribution in the Urban Environment: The Cadastral-based Expert Dasymetric System (CEDS) , 2007 .

[17]  Sunil Kumar,et al.  Maxent modeling for predicting suitable habitat for threatened and endangered tree Canacomyrica monticola in New Caledonia , 2009 .

[18]  C. Elvidge,et al.  Spatial analysis of global urban extent from DMSP-OLS night lights , 2005 .

[19]  P. Sutton.  Modeling population density with night-time satellite imagery and GIS , 1997 .

[20]  Jeremy Mennis,et al.  Intelligent Dasymetric Mapping and Its Application to Areal Interpolation , 2006 .

[21]  Brian C. O'Neill,et al.  A guide to global population projections. , 2001 .

[22]  Sisira Jayasuriya,et al.  Environment and Natural Resources , 2003 .

[23]  M. Araújo,et al.  An evaluation of methods for modelling species distributions , 2004 .

[24]  M. Araújo,et al.  Presence-absence versus presence-only modelling methods for predicting bird habitat suitability , 2004 .

[25]  B. Efron,et al.  A Leisurely Look at the Bootstrap, the Jackknife, and , 1983 .

[26]  Paul M. Torrens,et al.  Geosimulation: object-based modeling of urban phenomena , 2004, Comput. Environ. Urban Syst..

[27]  T. Oki,et al.  A SRES-based gridded global population dataset for 1990–2100 , 2007 .

[28]  Mitchel Langford,et al.  Rapid facilitation of dasymetric-based population interpolation by means of raster pixel maps , 2007, Comput. Environ. Urban Syst..

[29]  Heather Booth,et al.  Demographic forecasting: 1980 to 2005 in review , 2006 .

[30]  W. Tobler.  Smooth pycnophylactic interpolation for geographical regions. , 1979, Journal of the American Statistical Association.

[31]  Claude E. Shannon,et al.  A Mathematical Theory of Communications , 1948 .

[32]  Edzer J. Pebesma,et al.  Applied Spatial Data Analysis with R - Second Edition , 2008, Use R!.