Assessing the use of global land cover data for guiding large area population distribution modelling

High-resolution gridded population datasets for Latin America and the Caribbean using official statistics

"Leaving no one behind" is the fundamental objective of the 2030 Agenda for Sustainable Development. Latin America and the Caribbean is marked by social inequalities, whilst its total population is projected to increase to almost 760 million by 2050. In this context, contemporary and spatially detailed datasets that accurately capture the distribution of residential population are critical to appropriately inform and support environmental, health, and developmental applications at subnational levels. Existing datasets are under-utilised by governments due to the non-alignment with their own statistics. Therefore, official statistics at the finest level of administrative units available have been implemented to construct an open-access repository of high-resolution gridded population datasets for 40 countries in Latin American and the Caribbean. These datasets are detailed here, alongside the 'top-down' approach and methods to generate and validate them. Population distribution datasets for each country were created at a resolution of 3 arc-seconds (approximately 100 m at the equator), and are all available from the WorldPop Data Repository.

Exploring the Relationship between the Spatial Distribution of Different Age Populations and Points of Interest (POI) in China

Population spatialization data is crucial to conducting scientific studies of coupled human–environment systems. Although significant progress has been made in population spatialization, the spatialization of different age populations is still weak. POI data with rich information have great potential to simulate the spatial distribution of different age populations, but the relationship between spatial distributions of POI and different age populations is still unclear, and whether it can be used as an auxiliary variable for the different age population spatialization remains to be explored. Therefore, this study collected and sorted out the number of different age populations and POIs in 2846 county-level administrative units of the Chinese mainland in 2010, divided the research data by region and city size, and explored the relationship between the different age populations and POIs. We found that there is a complex relationship between POI and different age populations. Firstly, there are positive, moderate-to-strong linear correlations between POI and population indicators. Secondly, POI has a different explanatory power for different age populations, and it has a higher explanatory power for the young and middle-aged population than the child and old population. Thirdly, the explanatory power of POI to different age populations is positively correlated with the urban economic development level. Finally, a small number of a certain kinds of POIs can be used to effectively simulate the spatial distributions of different age populations, which can improve the efficiency of obtaining spatialization data of different age populations and greatly save on costs. The study can provide data support for the precise spatialization of different age populations and inspire the spatialization of the other population attributes by POI in the future.

Evaluating the Accuracy of Gridded Population Estimates in Slums: A Case Study in Nigeria and Kenya

Low- and middle-income country cities face unprecedented urbanization and growth in slums. Gridded population data (e.g., ~100 × 100 m) derived from demographic and spatial data are a promising source of population estimates, but face limitations in slums due to the dynamic nature of this population as well as modelling assumptions. In this study, we compared field-referenced boundaries and population counts from Slum Dwellers International in Lagos (Nigeria), Port Harcourt (Nigeria), and Nairobi (Kenya) with nine gridded population datasets to assess their statistical accuracy in slums. We found that all gridded population estimates vastly underestimated population in slums (RMSE: 4958 to 14,422, Bias: −2853 to −7638), with the most accurate dataset (HRSL) estimating just 39 per cent of slum residents. Using a modelled map of all slums in Lagos to compare gridded population datasets in terms of SDG 11.1.1 (percent of population living in deprived areas), all gridded population datasets estimated this indicator at just 1–3 per cent compared to 56 per cent using UN-Habitat’s approach. We outline steps that might improve that accuracy of each gridded population dataset in deprived urban areas. While gridded population estimates are not yet sufficiently accurate to estimate SDG 11.1.1, we are optimistic that some could be used in the future following updates to their modelling approaches.

Gridded population mapping for Germany based on building density, height and type from Earth Observation data using census disaggregation and bottom-up estimates

Gridded population data is widely used to map fine scale population patterns and dynamics to understand associated human-environmental processes for global change research, disaster risk assessment and other domains. This study mapped gridded population across Germany using weighting layers from building density, building height (both from previous studies) and building type datasets, all created from freely available, temporally and globally consistent Copernicus Sentinel-1 and Sentinel-2 data. We first produced and validated a nation-wide dataset of predominant residential and non-residential building types. We then examined the impact of different weighting layers from density, type and height on top-down dasymetric mapping quality across scales. We finally performed a nation-wide bottom-up population estimate based on the three datasets. We found that integrating building types into dasymetric mapping is helpful at fine scale, as population is not redistributed to non-residential areas. Building density improved the overall quality of population estimates at all scales compared to using a binary building layer. Most importantly, we found that the combined use of density and height, i.e. volume, considerably increased mapping quality in general and with regard to regional discrepancy by largely eliminating systematic underestimation in dense agglomerations and overestimation in rural areas. We also found that building density, type and volume, together with living floor area per capita, are suitable to produce accurate large-area bottom-up population estimates.

Intraday Variation Mapping of Population Age Structure via Urban-Functional-Region-Based Scaling

The spatial distribution of the population is uneven for various reasons, such as urban-rural differences and geographical conditions differences. As the basic element of the natural structure of the population, the age structure composition of populations also varies considerably across the world. Obtaining accurate and spatiotemporal population age structure maps is crucial for calculating population size at risk, analyzing populations mobility patterns, or calculating health and development indicators. During the past decades, many population maps in the form of administrative units and grids have been produced. However, these population maps are limited by the lack of information on the change of population distribution within a day and the age structure of the population. Urban functional regions (UFRs) are closely related to population mobility patterns, which can provide information about population variation intraday. Focusing on the area within the Beijing Fifth Ring Road, the political and economic center of Beijing, we showed how to use the temporal scaling factors obtained by analyzing the population survey sampling data and population dasymetric maps in different categories of UFRs to realize the intraday variation mapping of elderly individuals and children. The population dasymetric maps were generated on the basis of covariates related to population. In this article, 50 covariates were calculated from remote sensing data and geospatial data. However, not all covariates are associate with population distribution. In order to improve the accuracy of dasymetric maps and reduce the cost of mapping, it is necessary to select the optimal subset for the dasymetric model of elderly and children. The random forest recursive feature elimination (RF-RFE) algorithm was introduced to obtain the optimal subset of different age groups of people and generate the population dasymetric model in this article, as well as to screen out the optimal subset with 38 covariates and 26 covariates for the dasymetric models of the elderly and children, respectively. An accurate UFR identification method combining point of interest (POI) data and OpenStreetMap (OSM) road network data is also introduced in this article. The overall accuracy of the identification results of UFRs was 70.97%, which is quite accurate. The intraday variation maps of population age structure on weekdays and weekends were made within the Beijing Fifth Ring Road. Accuracy evaluation based on sampling data found that the overall accuracy was relatively high—R2 for each time period was higher than 0.5 and root mean square error (RMSE) was less than 0.05. On weekdays in particular, R2 for each time period was higher than 0.61 and RMSE was less than 0.02.

Gridded population survey sampling: a systematic scoping review of the field and strategic research agenda

INTRODUCTION

In low- and middle-income countries (LMICs), household survey data are a main source of information for planning, evaluation, and decision-making. Standard surveys are based on censuses, however, for many LMICs it has been more than 10 years since their last census and they face high urban growth rates. Over the last decade, survey designers have begun to use modelled gridded population estimates as sample frames. We summarize the state of the emerging field of gridded population survey sampling, focussing on LMICs.

METHODS

We performed a systematic scoping review in Scopus of specific gridded population datasets and "population" or "household" "survey" reports, and solicited additional published and unpublished sources from colleagues.

RESULTS

We identified 43 national and sub-national gridded population-based household surveys implemented across 29 LMICs. Gridded population surveys used automated and manual approaches to derive clusters from WorldPop and LandScan gridded population estimates. After sampling, some survey teams interviewed all households in each cluster or segment, and others sampled households from larger clusters. Tools to select gridded population survey clusters include the GridSample R package, Geo-sampling tool, and GridSample.org. In the field, gridded population surveys generally relied on geographically accurate maps based on satellite imagery or OpenStreetMap, and a tablet or GPS technology for navigation.

CONCLUSIONS

For gridded population survey sampling to be adopted more widely, several strategic questions need answering regarding cell-level accuracy and uncertainty of gridded population estimates, the methods used to group/split cells into sample frame units, design effects of new sample designs, and feasibility of tools and methods to implement surveys across diverse settings.

Local Population Mapping Using a Random Forest Model Based on Remote and Social Sensing Data: A Case Study in Zhengzhou, China

High-resolution gridded population data are important for understanding and responding to many socioeconomic and environmental problems. Local estimates of the population allow officials and researchers to make a better local planning (e.g., optimizing public services and facilities). This study used a random forest algorithm, on the basis of remote sensing (i.e., satellite imagery) and social sensing data (i.e., point-of-interest and building footprint), to disaggregate census population data for the five municipal districts of Zhengzhou city, China, onto 100 × 100 m grid cells. We used a statistical tool to detect areas with an abnormal population density; e.g., areas containing many empty houses or houses rented by more people than allowed, and conducted field work to validate our findings. Results showed that some categories of points-of-interest, such as residential communities, parking lots, banks, and government buildings were the most important contributing elements in modeling the spatial distribution of the residential population in Zhengzhou City. The exclusion of areas with an abnormal population density from model training and dasymetric mapping increased the accuracy of population estimates in other areas with a more common population density. We compared our product with three widely used gridded population products: Worldpop, the Gridded Population of the World, and the 1-km Grid Population Dataset of China. The relative accuracy of our modeling approach was higher than that of those three products in the five municipal districts of Zhengzhou. This study demonstrated potential for the combination of remote and social sensing data to more accurately estimate the population density in urban areas, with minimum disturbance from the abnormal population density.

Analysis of accessibility to emergency rooms by dynamic population from mobile phone data: Geography of social inequity in South Korea

Accessibility of emergency medical care is one of the crucial factors in evaluating national primary medical care systems. While many studies have focused on this issue, there was a fundamental limit to the measurement of accessibility of emergency rooms, because the commonly used census-based population data are difficult to provide realistic information in terms of time and space. In this study, we evaluated the geographical accessibility of emergency rooms in South Korea by using dynamic population counts from mobile phone data. Such population counts were more accurate and up-to-date because they are obtained by aggregating the number of mobile phone users in a 50-by-50 m grid of a locational field, weighted by stay time. Considering both supply and demand of emergency rooms, the 2-step floating catchment analysis was implemented. As a result, urban areas, including the capital city Seoul, showed lower accessibility to emergency rooms, whereas rural areas recorded higher accessibility. This result was contrary to the results analyzed by us based on census-based population data: higher accessibility in urban areas and lower in rural. This implies that using solely census data for accessibility analysis could lead to certain errors, and adopting mobile-based population data would represent the real-world situations for solving problems of social inequity in primary medical care.

Downscaling livestock census data using multivariate predictive models: Sensitivity to modifiable areal unit problem

The analysis of census data aggregated by administrative units introduces a statistical bias known as the modifiable areal unit problem (MAUP). Previous researches have mostly assessed the effect of MAUP on upscaling models. The present study contributes to clarify the effects of MAUP on the downscaling methodologies, highlighting how a priori choices of scales and shapes could influence the results. We aggregated chicken and duck fine-resolution census in Thailand, using three administrative census levels in regular and irregular shapes. We then disaggregated the data within the Gridded Livestock of the World analytical framework, sampling predictors in two different ways. A sensitivity analysis on Pearson's r correlation statistics and RMSE was carried out to understand how size and shapes of the response variables affect the goodness-of-fit and downscaling performances. We showed that scale, rather than shapes and sampling methods, affected downscaling precision, suggesting that training the model using the finest administrative level available is preferable. Moreover, datasets showing non-homogeneous distribution but instead spatial clustering seemed less affected by MAUP, yielding higher Pearson's r values and lower RMSE compared to a more spatially homogenous dataset. Implementing aggregation sensitivity analysis in spatial studies could help to interpret complex results and disseminate robust products.

Evaluation of the Equity of Urban Park Green Space Based on Population Data Spatialization: A Case Study of a Central Area of Wuhan, China

To measure the equity of urban park green space, spatial matching between service supply and user group demand should be taken into consideration. However, if the demographic data, with the administrative division as the basic unit, are directly applied to characterize the spatial distribution of a user group, it may introduce inevitable deviation into the evaluation results due to the low-resolution nature and modifiable areal unit problem of such data. Taking the central area of Wuhan as an example, the population data spatialization method based on land use modeling was used to build a geographically weighted regression (GWR) model of land cover type and demographic data, and the spatial distribution of the population of the 150 m grid was obtained by inversion. Then, the equity of park green space in Wuhan central city was evaluated by population spatial data and network accessibility. The results showed that (1) the range of park green space in the central urban area of Wuhan was within a walking distance of 15 min, accounting for 25.8% of the total study area and covering 54.2% of the population in the study area; (2) the equity of park green space in Hongshan District was the worst; (3) and the use of population spatial data can measure equity on a more precise scale.

Global spatio-temporally harmonised datasets for producing high-resolution gridded population distribution datasets

Multi-temporal, globally consistent, high-resolution human population datasets provide consistent and comparable population distributions in support of mapping sub-national heterogeneities in health, wealth, and resource access, and monitoring change in these over time. The production of more reliable and spatially detailed population datasets is increasingly necessary due to the importance of improving metrics at sub-national and multi-temporal scales. This is in support of measurement and monitoring of UN Sustainable Development Goals and related agendas. In response to these agendas, a method has been developed to assemble and harmonise a unique, open access, archive of geospatial datasets. Datasets are provided as global, annual time series, where pertinent at the timescale of population analyses and where data is available, for use in the construction of population distribution layers. The archive includes sub-national census-based population estimates, matched to a geospatial layer denoting administrative unit boundaries, and a number of co-registered gridded geospatial factors that correlate strongly with population presence and density. Here, we describe these harmonised datasets and their limitations, along with the production workflow. Further, we demonstrate applications of the archive by producing multi-temporal gridded population outputs for Africa and using these to derive health and development metrics. The geospatial archive is available at https://doi.org/10.5258/SOTON/WP00650.

Using Information on Settlement Patterns to Improve the Spatial Distribution of Population in Coastal Impact Assessments

Broad-scale impact and vulnerability assessments are essential for informing decisions on long-term adaptation planning at the national, regional, or global level. These assessments rely on population data for quantifying exposure to different types of hazards. Existing population datasets covering the entire globe at resolutions of 2.5 degrees to 30 arc-seconds are based on information available at administrative-unit level and implicitly assume uniform population densities within these units. This assumption can lead to errors in impact assessments and particularly in coastal areas that are densely populated. This study proposes and compares simple approaches to regionalize population within administrative units in the German Baltic Sea region using solely information on urban extent from the Global Urban Footprint (GUF). Our results show that approaches using GUF can reduce the error in predicting population totals of municipalities by factor 2 to 3. When assessing exposed population, we find that the assumption of uniform population densities leads to an overestimation of 120% to 140%. Using GUF to regionalise population within administrative units reduce these errors by up to 50%. Our results suggest that the proposed simple modeling approaches can result in significantly improved distribution of population within administrative units and substantially improve the results of exposure analyses.

Data-enriched Interpolation for Temporally Consistent Population Compositions

This research evaluates the performance of areal interpolation coupled with dasymetric refinement to estimate different demographic attributes, namely population sub-groups based on race, age structure and urban residence, within consistent census tract boundaries from 1990 to 2010 in Massachusetts. The creation of such consistent estimates facilitates the study of the nuanced micro-scale evolution of different aspects of population, which is impossible using temporally incompatible small-area census geographies from different points in time. Various unexplored ancillary variables, including the Global Human Settlement Layer (GHSL), the National Land-Cover Database (NLCD), parcels, building footprints and the proprietary ZTRAX^® dataset are utilized for dasymetric refinement prior to areal interpolation to examine their effectiveness in improving the accuracy of multi-temporal population estimates. Different areal interpolation methods including Areal Weighting (AW), Target Density Weighting (TDW), Expectation Maximization (EM) and its data-extended approach are coupled with different dasymetric refinement scenarios based on these ancillary variables. The resulting consistent small area estimates of white and black subpopulations, people of age 18-65 and urban population show that dasymetrically refined areal interpolation is particularly effective when the analysis spans a longer time period (1990-2010 instead of 2000-2010) and the enumerated population is sufficiently large (e.g., counts of white vs. black). The results also demonstrate that current census-defined urban areas overestimate the spatial distribution of urban population and dasymetrically refined areal interpolation improves estimates of urban population. Refined TDW using building footprints or the ZTRAX^® dataset outperforms all other methods. The implementation of areal interpolation enriched by dasymetric refinement represents a promising strategy to create more reliable multi-temporal and consistent estimates of different population subgroups and thus demographic compositions. This methodological foundation has the potential to advance micro-scale modeling of various subpopulations, particularly urban population to inform studies of urbanization and population change over time as well as future population projections.

Examining the correlates and drivers of human population distributions across low- and middle-income countries

Geographical factors have influenced the distributions and densities of global human population distributions for centuries. Climatic regimes have made some regions more habitable than others, harsh topography has discouraged human settlement, and transport links have encouraged population growth. A better understanding of these types of relationships enables both improved mapping of population distributions today and modelling of future scenarios. However, few comprehensive studies of the relationships between population spatial distributions and the range of drivers and correlates that exist have been undertaken at all, much less at high spatial resolutions, and particularly across the low- and middle-income countries. Here, we quantify the relative importance of multiple types of drivers and covariates in explaining observed population densities across 32 low- and middle-income countries over four continents using machine-learning approaches. We find that, while relationships between population densities and geographical factors show some variation between regions, they are generally remarkably consistent, pointing to universal drivers of human population distribution. Here, we find that a set of geographical features relating to the built environment, ecology and topography consistently explain the majority of variability in population distributions at fine spatial scales across the low- and middle-income regions of the world.

High-resolution reconstruction of the United States human population distribution, 1790 to 2010

Where do people live, and how has this changed over timescales of centuries? High-resolution spatial information on historical human population distribution is of great significance to understand human-environment interactions and their temporal dynamics. However, the complex relationship between population distribution and various influencing factors coupled with limited data availability make it a challenge to reconstruct human population distribution over timescales of centuries. This study generated 1-km decadal population maps for the conterminous US from 1790 to 2010 using parsimonious models based on natural suitability, socioeconomic desirability, and inhabitability. Five models of increasing complexity were evaluated. The models were validated with census tract and county subdivision population data in 2000 and were applied to generate five sets of 22 historical population maps from 1790-2010. Separating urban and rural areas and excluding non-inhabitable areas were the most important factors for improving the overall accuracy. The generated gridded population datasets and the production and validation methods are described here.

A large-area, spatially continuous assessment of land cover map error and its impact on downstream analyses

Land cover maps increasingly underlie research into socioeconomic and environmental patterns and processes, including global change. It is known that map errors impact our understanding of these phenomena, but quantifying these impacts is difficult because many areas lack adequate reference data. We used a highly accurate, high-resolution map of South African cropland to assess (1) the magnitude of error in several current generation land cover maps, and (2) how these errors propagate in downstream studies. We first quantified pixel-wise errors in the cropland classes of four widely used land cover maps at resolutions ranging from 1 to 100 km, and then calculated errors in several representative "downstream" (map-based) analyses, including assessments of vegetative carbon stocks, evapotranspiration, crop production, and household food security. We also evaluated maps' spatial accuracy based on how precisely they could be used to locate specific landscape features. We found that cropland maps can have substantial biases and poor accuracy at all resolutions (e.g., at 1 km resolution, up to ∼45% underestimates of cropland (bias) and nearly 50% mean absolute error (MAE, describing accuracy); at 100 km, up to 15% underestimates and nearly 20% MAE). National-scale maps derived from higher-resolution imagery were most accurate, followed by multi-map fusion products. Constraining mapped values to match survey statistics may be effective at minimizing bias (provided the statistics are accurate). Errors in downstream analyses could be substantially amplified or muted, depending on the values ascribed to cropland-adjacent covers (e.g., with forest as adjacent cover, carbon map error was 200%-500% greater than in input cropland maps, but ∼40% less for sparse cover types). The average locational error was 6 km (600%). These findings provide deeper insight into the causes and potential consequences of land cover map error, and suggest several recommendations for land cover map users.

Modelling changing population distributions: an example of the Kenyan Coast, 1979-2009

Large-scale gridded population datasets are usually produced for the year of input census data using a top-down approach and projected backward and forward in time using national growth rates. Such temporal projections do not include any subnational variation in population distribution trends and ignore changes in geographical covariates such as urban land cover changes. Improved predictions of population distribution changes over time require the use of a limited number of covariates that are time-invariant or temporally explicit. Here we make use of recently released multi-temporal high-resolution global settlement layers, historical census data and latest developments in population distribution modelling methods to reconstruct population distribution changes over 30 years across the Kenyan Coast. We explore the methodological challenges associated with the production of gridded population distribution time-series in data-scarce countries and show that trade-offs have to be found between spatial and temporal resolutions when selecting the best modelling approach. Strategies used to fill data gaps may vary according to the local context and the objective of the study. This work will hopefully serve as a benchmark for future developments of population distribution time-series that are increasingly required for population-at-risk estimations and spatial modelling in various fields.

Spatial accessibility to basic public health services in South Sudan

At independence in 2011, South Sudan's health sector was almost non-existent. The first national health strategic plan aimed to achieve an integrated health facility network that would mean that 70% of the population were within 5 km of a health service provider. Publically available data on functioning and closed health facilities, population distribution, road networks, land use and elevation were used to compute the fraction of the population within 1 hour walking distance of the nearest public health facility offering curative services. This metric was summarised for each of the 78 counties in South Sudan and compared with simpler metrics of the proportion of the population within 5 km of a health facility. In 2016, it is estimated that there were 1747 public health facilities, out of which 294 were non-functional in part due to the on-going civil conflict. Access to a service provider was poor with only 25.7% of the population living within one-hour walking time to a facility and 28.6% of the population within 5 km. These metrics, when applied sub-nationally, identified the same high priority, most vulnerable counties. Simple metrics based upon population distribution and location of facilities might be as valuable as more complex models of health access, where attribute data on travel routes are imperfect or incomplete and sparse. Disparities exist in South Sudan among counties and those with the poorest health access should be targeted for priority expansion of clinical services.

Improving Large Area Population Mapping Using Geotweet Densities

Many different methods are used to disaggregate census data and predict population densities to construct finer scale, gridded population data sets. These methods often involve a range of high resolution geospatial covariate datasets on aspects such as urban areas, infrastructure, land cover and topography; such covariates, however, are not directly indicative of the presence of people. Here we tested the potential of geo-located tweets from the social media application, Twitter, as a covariate in the production of population maps. The density of geo-located tweets in 1x1 km grid cells over a 2-month period across Indonesia, a country with one of the highest Twitter usage rates in the world, was input as a covariate into a previously published random forests-based census disaggregation method. Comparison of internal measures of accuracy and external assessments between models built with and without the geotweets showed that increases in population mapping accuracy could be obtained using the geotweet densities as a covariate layer. The work highlights the potential for such social media-derived data in improving our understanding of population distributions and offers promise for more dynamic mapping with such data being continually produced and freely available.

Comprehensive framework for visualizing and analyzing spatio-temporal dynamics of racial diversity in the entire United States

The United States is increasingly becoming a multi-racial society. To understand multiple consequences of this overall trend to our neighborhoods we need a methodology capable of spatio-temporal analysis of racial diversity at the local level but also across the entire U.S. Furthermore, such methodology should be accessible to stakeholders ranging from analysts to decision makers. In this paper we present a comprehensive framework for visualizing and analyzing diversity data that fulfills such requirements. The first component of our framework is a U.S.-wide, multi-year database of race sub-population grids which is freely available for download. These 30 m resolution grids have being developed using dasymetric modeling and are available for 1990-2000-2010. We summarize numerous advantages of gridded population data over commonly used Census tract-aggregated data. Using these grids frees analysts from constructing their own and allows them to focus on diversity analysis. The second component of our framework is a set of U.S.-wide, multi-year diversity maps at 30 m resolution. A diversity map is our product that classifies the gridded population into 39 communities based on their degrees of diversity, dominant race, and population density. It provides spatial information on diversity in a single, easy-to-understand map that can be utilized by analysts and end users alike. Maps based on subsequent Censuses provide information about spatio-temporal dynamics of diversity. Diversity maps are accessible through the GeoWeb application SocScape (http://sil.uc.edu/webapps/socscape_usa/) for an immediate online exploration. The third component of our framework is a proposal to quantitatively analyze diversity maps using a set of landscape metrics. Because of its form, a grid-based diversity map could be thought of as a diversity “landscape” and analyzed quantitatively using landscape metrics. We give a brief summary of most pertinent metrics and demonstrate how they can be applied to diversity maps.

High resolution global gridded data for use in population studies

Recent years have seen substantial growth in openly available satellite and other geospatial data layers, which represent a range of metrics relevant to global human population mapping at fine spatial scales. The specifications of such data differ widely and therefore the harmonisation of data layers is a prerequisite to constructing detailed and contemporary spatial datasets which accurately describe population distributions. Such datasets are vital to measure impacts of population growth, monitor change, and plan interventions. To this end the WorldPop Project has produced an open access archive of 3 and 30 arc-second resolution gridded data. Four tiled raster datasets form the basis of the archive: (i) Viewfinder Panoramas topography clipped to Global ADMinistrative area (GADM) coastlines; (ii) a matching ISO 3166 country identification grid; (iii) country area; (iv) and slope layer. Further layers include transport networks, landcover, nightlights, precipitation, travel time to major cities, and waterways. Datasets and production methodology are here described. The archive can be downloaded both from the WorldPop Dataverse Repository and the WorldPop Project website.

Research on Grid Size Suitability of Gridded Population Distribution in Urban Area: A Case Study in Urban Area of Xuanzhou District, China

The research on the grid size suitability is important to provide improvement in accuracies of gridded population distribution. It contributes to reveal the actual spatial distribution of population. However, currently little research has been done in this area. Many well-modeled gridded population dataset are basically built at a single grid scale. If the grid cell size is not appropriate, it will result in spatial information loss or data redundancy. Therefore, in order to capture the desired spatial variation of population within the area of interest, it is necessary to conduct research on grid size suitability. This study summarized three expressed levels to analyze grid size suitability, which include location expressed level, numeric information expressed level, and spatial relationship expressed level. This study elaborated the reasons for choosing the five indexes to explore expression suitability. These five indexes are consistency measure, shape index rate, standard deviation of population density, patches diversity index, and the average local variance. The suitable grid size was determined by constructing grid size-indicator value curves and suitable grid size scheme. Results revealed that the three expressed levels on 10m grid scale are satisfying. And the population distribution raster data with 10m grid size provide excellent accuracy without loss. The 10m grid size is recommended as the appropriate scale for generating a high-quality gridded population distribution in our study area. Based on this preliminary study, it indicates the five indexes are coordinated with each other and reasonable and effective to assess grid size suitability. We also suggest choosing these five indexes in three perspectives of expressed level to carry out the research on grid size suitability of gridded population distribution.

Fine resolution mapping of population age-structures for health and development applications

The age-group composition of populations varies considerably across the world, and obtaining accurate, spatially detailed estimates of numbers of children under 5 years is important in designing vaccination strategies, educational planning or maternal healthcare delivery. Traditionally, such estimates are derived from population censuses, but these can often be unreliable, outdated and of coarse resolution for resource-poor settings. Focusing on Nigeria, we use nationally representative household surveys and their cluster locations to predict the proportion of the under-five population in 1 × 1 km using a Bayesian hierarchical spatio-temporal model. Results showed that land cover, travel time to major settlements, night-time lights and vegetation index were good predictors and that accounting for fine-scale variation, rather than assuming a uniform proportion of under 5 year olds can result in significant differences in health metrics. The largest gaps in estimated bednet and vaccination coverage were in Kano, Katsina and Jigawa. Geolocated household surveys are a valuable resource for providing detailed, contemporary and regularly updated population age-structure data in the absence of recent census data. By combining these with covariate layers, age-structure maps of unprecedented detail can be produced to guide the targeting of interventions in resource-poor settings.

High-resolution gridded population datasets for Latin America and the Caribbean in 2010, 2015, and 2020

The Latin America and the Caribbean region is one of the most urbanized regions in the world, with a total population of around 630 million that is expected to increase by 25% by 2050. In this context, detailed and contemporary datasets accurately describing the distribution of residential population in the region are required for measuring the impacts of population growth, monitoring changes, supporting environmental and health applications, and planning interventions. To support these needs, an open access archive of high-resolution gridded population datasets was created through disaggregation of the most recent official population count data available for 28 countries located in the region. These datasets are described here along with the approach and methods used to create and validate them. For each country, population distribution datasets, having a resolution of 3 arc seconds (approximately 100 m at the equator), were produced for the population count year, as well as for 2010, 2015, and 2020. All these products are available both through the WorldPop Project website and the WorldPop Dataverse Repository.

Disaggregating census data for population mapping using random forests with remotely-sensed and ancillary data

High resolution, contemporary data on human population distributions are vital for measuring impacts of population growth, monitoring human-environment interactions and for planning and policy development. Many methods are used to disaggregate census data and predict population densities for finer scale, gridded population data sets. We present a new semi-automated dasymetric modeling approach that incorporates detailed census and ancillary data in a flexible, "Random Forest" estimation technique. We outline the combination of widely available, remotely-sensed and geospatial data that contribute to the modeled dasymetric weights and then use the Random Forest model to generate a gridded prediction of population density at ~100 m spatial resolution. This prediction layer is then used as the weighting surface to perform dasymetric redistribution of the census counts at a country level. As a case study we compare the new algorithm and its products for three countries (Vietnam, Cambodia, and Kenya) with other common gridded population data production methodologies. We discuss the advantages of the new method and increases over the accuracy and flexibility of those previous approaches. Finally, we outline how this algorithm will be extended to provide freely-available gridded population data sets for Africa, Asia and Latin America.

Mapping for maternal and newborn health: the distributions of women of childbearing age, pregnancies and births

BACKGROUND

The health and survival of women and their new-born babies in low income countries has been a key priority in public health since the 1990s. However, basic planning data, such as numbers of pregnancies and births, remain difficult to obtain and information is also lacking on geographic access to key services, such as facilities with skilled health workers. For maternal and newborn health and survival, planning for safer births and healthier newborns could be improved by more accurate estimations of the distributions of women of childbearing age. Moreover, subnational estimates of projected future numbers of pregnancies are needed for more effective strategies on human resources and infrastructure, while there is a need to link information on pregnancies to better information on health facilities in districts and regions so that coverage of services can be assessed.

METHODS

This paper outlines demographic mapping methods based on freely available data for the production of high resolution datasets depicting estimates of numbers of people, women of childbearing age, live births and pregnancies, and distribution of comprehensive EmONC facilities in four large high burden countries: Afghanistan, Bangladesh, Ethiopia and Tanzania. Satellite derived maps of settlements and land cover were constructed and used to redistribute areal census counts to produce detailed maps of the distributions of women of childbearing age. Household survey data, UN statistics and other sources on growth rates, age specific fertility rates, live births, stillbirths and abortions were then integrated to convert the population distribution datasets to gridded estimates of births and pregnancies.

RESULTS AND CONCLUSIONS

These estimates, which can be produced for current, past or future years based on standard demographic projections, can provide the basis for strategic intelligence, planning services, and provide denominators for subnational indicators to track progress. The datasets produced are part of national midwifery workforce assessments conducted in collaboration with the respective Ministries of Health and the United Nations Population Fund (UNFPA) to identify disparities between population needs, health infrastructure and workforce supply. The datasets are available to the respective Ministries as part of the UNFPA programme to inform midwifery workforce planning and also publicly available through the WorldPop population mapping project.

Modelling spatial patterns of urban growth in Africa

The population of Africa is predicted to double over the next 40 years, driving exceptionally high urban expansion rates that will induce significant socio-economic, environmental and health changes. In order to prepare for these changes, it is important to better understand urban growth dynamics in Africa and better predict the spatial pattern of rural-urban conversions. Previous work on urban expansion has been carried out at the city level or at the global level with a relatively coarse 5-10 km resolution. The main objective of the present paper was to develop a modelling approach at an intermediate scale in order to identify factors that influence spatial patterns of urban expansion in Africa. Boosted Regression Tree models were developed to predict the spatial pattern of rural-urban conversions in every large African city. Urban change data between circa 1990 and circa 2000 available for 20 large cities across Africa were used as training data. Results showed that the urban land in a 1 km neighbourhood and the accessibility to the city centre were the most influential variables. Results obtained were generally more accurate than results obtained using a distance-based urban expansion model and showed that the spatial pattern of small, compact and fast growing cities were easier to simulate than cities with lower population densities and a lower growth rate. The simulation method developed here will allow the production of spatially detailed urban expansion forecasts for 2020 and 2025 for Africa, data that are increasingly required by global change modellers.

Millennium development health metrics: where do Africa's children and women of childbearing age live?

The Millennium Development Goals (MDGs) have prompted an expansion in approaches to deriving health metrics to measure progress toward their achievement. Accurate measurements should take into account the high degrees of spatial heterogeneity in health risks across countries, and this has prompted the development of sophisticated cartographic techniques for mapping and modeling risks. Conversion of these risks to relevant population-based metrics requires equally detailed information on the spatial distribution and attributes of the denominator populations. However, spatial information on age and sex composition over large areas is lacking, prompting many influential studies that have rigorously accounted for health risk heterogeneities to overlook the substantial demographic variations that exist subnationally and merely apply national-level adjustments.

Here we outline the development of high resolution age- and sex-structured spatial population datasets for Africa in 2000-2015 built from over a million measurements from more than 20,000 subnational units, increasing input data detail from previous studies by over 400-fold. We analyze the large spatial variations seen within countries and across the continent for key MDG indicator groups, focusing on children under 5 and women of childbearing age, and find that substantial differences in health and development indicators can result through using only national level statistics, compared to accounting for subnational variation.

Progress toward meeting the MDGs will be measured through national-level indicators that mask substantial inequalities and heterogeneities across nations. Cartographic approaches are providing opportunities for quantitative assessments of these inequalities and the targeting of interventions, but demographic spatial datasets to support such efforts remain reliant on coarse and outdated input data for accurately locating risk groups. We have shown here that sufficient data exist to map the distribution of key vulnerable groups, and that doing so has substantial impacts on derived metrics through accounting for spatial demographic heterogeneities that exist within nations across Africa.

High resolution population distribution maps for Southeast Asia in 2010 and 2015

Spatially accurate, contemporary data on human population distributions are vitally important to many applied and theoretical researchers. The Southeast Asia region has undergone rapid urbanization and population growth over the past decade, yet existing spatial population distribution datasets covering the region are based principally on population count data from censuses circa 2000, with often insufficient spatial resolution or input data to map settlements precisely. Here we outline approaches to construct a database of GIS-linked circa 2010 census data and methods used to construct fine-scale (∼100 meters spatial resolution) population distribution datasets for each country in the Southeast Asia region. Landsat-derived settlement maps and land cover information were combined with ancillary datasets on infrastructure to model population distributions for 2010 and 2015. These products were compared with those from two other methods used to construct commonly used global population datasets. Results indicate mapping accuracies are consistently higher when incorporating land cover and settlement information into the AsiaPop modelling process. Using existing data, it is possible to produce detailed, contemporary and easily updatable population distribution datasets for Southeast Asia. The 2010 and 2015 datasets produced are freely available as a product of the AsiaPop Project and can be downloaded from: www.asiapop.org.

Genetics: A New Landscape for Medical Geography

The emergence and re-emergence of human pathogens resistant to medical treatment will present a challenge to the international public health community in the coming decades. Geography is uniquely positioned to examine the progressive evolution of pathogens across space and through time, and to link molecular change to interactions between population and environmental drivers. Landscape as an organizing principle for the integration of natural and cultural forces has a long history in geography, and, more specifically, in medical geography. Here, we explore the role of landscape in medical geography, the emergent field of landscape genetics, and the great potential that exists in the combination of these two disciplines. We argue that landscape genetics can enhance medical geographic studies of local-level disease environments with quantitative tests of how human-environment interactions influence pathogenic characteristics. In turn, such analyses can expand theories of disease diffusion to the molecular scale and distinguish the important factors in ecologies of disease that drive genetic change of pathogens.

Large-scale spatial population databases in infectious disease research

Modelling studies on the spatial distribution and spread of infectious diseases are becoming increasingly detailed and sophisticated, with global risk mapping and epidemic modelling studies now popular. Yet, in deriving populations at risk of disease estimates, these spatial models must rely on existing global and regional datasets on population distribution, which are often based on outdated and coarse resolution data. Moreover, a variety of different methods have been used to model population distribution at large spatial scales. In this review we describe the main global gridded population datasets that are freely available for health researchers and compare their construction methods, and highlight the uncertainties inherent in these population datasets. We review their application in past studies on disease risk and dynamics, and discuss how the choice of dataset can affect results. Moreover, we highlight how the lack of contemporary, detailed and reliable data on human population distribution in low income countries is proving a barrier to obtaining accurate large-scale estimates of population at risk and constructing reliable models of disease spread, and suggest research directions required to further reduce these barriers.

Population distribution, settlement patterns and accessibility across Africa in 2010

The spatial distribution of populations and settlements across a country and their interconnectivity and accessibility from urban areas are important for delivering healthcare, distributing resources and economic development. However, existing spatially explicit population data across Africa are generally based on outdated, low resolution input demographic data, and provide insufficient detail to quantify rural settlement patterns and, thus, accurately measure population concentration and accessibility. Here we outline approaches to developing a new high resolution population distribution dataset for Africa and analyse rural accessibility to population centers. Contemporary population count data were combined with detailed satellite-derived settlement extents to map population distributions across Africa at a finer spatial resolution than ever before. Substantial heterogeneity in settlement patterns, population concentration and spatial accessibility to major population centres is exhibited across the continent. In Africa, 90% of the population is concentrated in less than 21% of the land surface and the average per-person travel time to settlements of more than 50,000 inhabitants is around 3.5 hours, with Central and East Africa displaying the longest average travel times. The analyses highlight large inequities in access, the isolation of many rural populations and the challenges that exist between countries and regions in providing access to services. The datasets presented are freely available as part of the AfriPop project, providing an evidence base for guiding strategic decisions.

Spatial modelling of healthcare utilisation for treatment of fever in Namibia

Background

Health care utilization is affected by several factors including geographic accessibility. Empirical data on utilization of health facilities is important to understanding geographic accessibility and defining health facility catchments at a national level. Accurately defining catchment population improves the analysis of gaps in access, commodity needs and interpretation of disease incidence. Here, empirical household survey data on treatment seeking for fever were used to model the utilisation of public health facilities and define their catchment areas and populations in northern Namibia.

Method

This study uses data from the Malaria Indicator Survey (MIS) of 2009 on treatment seeking for fever among children under the age of five years to characterize facility utilisation. Probability of attendance of public health facilities for fever treatment was modelled against a theoretical surface of travel times using a three parameter logistic model. The fitted model was then applied to a population surface to predict the number of children likely to use a public health facility during an episode of fever in northern Namibia.

Results

Overall, from the MIS survey, the prevalence of fever among children was 17.6% CI [16.0-19.1] (401 of 2,283 children) while public health facility attendance for fever was 51.1%, [95%CI: 46.2-56.0]. The coefficients of the logistic model of travel time against fever treatment at public health facilities were all significant (p < 0.001). From this model, probability of facility attendance remained relatively high up to 180 minutes (3 hours) and thereafter decreased steadily. Total public health facility catchment population of children under the age five was estimated to be 162,286 in northern Namibia with an estimated fever burden of 24,830 children. Of the estimated fevers, 8,021 (32.3%) were within 30 minutes of travel time to the nearest health facility while 14,902 (60.0%) were within 1 hour.

Conclusion

This study demonstrates the potential of routine household surveys to empirically model health care utilisation for the treatment of childhood fever and define catchment populations enhancing the possibilities of accurate commodity needs assessment and calculation of disease incidence. These methods could be extended to other African countries where detailed mapping of health facilities exists.

The effects of spatial population dataset choice on estimates of population at risk of disease

Background

The spatial modeling of infectious disease distributions and dynamics is increasingly being undertaken for health services planning and disease control monitoring, implementation, and evaluation. Where risks are heterogeneous in space or dependent on person-to-person transmission, spatial data on human population distributions are required to estimate infectious disease risks, burdens, and dynamics. Several different modeled human population distribution datasets are available and widely used, but the disparities among them and the implications for enumerating disease burdens and populations at risk have not been considered systematically. Here, we quantify some of these effects using global estimates of populations at risk (PAR) of P. falciparum malaria as an example.

Methods

The recent construction of a global map of P. falciparum malaria endemicity enabled the testing of different gridded population datasets for providing estimates of PAR by endemicity class. The estimated population numbers within each class were calculated for each country using four different global gridded human population datasets: GRUMP (~1 km spatial resolution), LandScan (~1 km), UNEP Global Population Databases (~5 km), and GPW3 (~5 km). More detailed assessments of PAR variation and accuracy were conducted for three African countries where census data were available at a higher administrative-unit level than used by any of the four gridded population datasets.

Results

The estimates of PAR based on the datasets varied by more than 10 million people for some countries, even accounting for the fact that estimates of population totals made by different agencies are used to correct national totals in these datasets and can vary by more than 5% for many low-income countries. In many cases, these variations in PAR estimates comprised more than 10% of the total national population. The detailed country-level assessments suggested that none of the datasets was consistently more accurate than the others in estimating PAR. The sizes of such differences among modeled human populations were related to variations in the methods, input resolution, and date of the census data underlying each dataset. Data quality varied from country to country within the spatial population datasets.

Conclusions

Detailed, highly spatially resolved human population data are an essential resource for planning health service delivery for disease control, for the spatial modeling of epidemics, and for decision-making processes related to public health. However, our results highlight that for the low-income regions of the world where disease burden is greatest, existing datasets display substantial variations in estimated population distributions, resulting in uncertainty in disease assessments that utilize them. Increased efforts are required to gather contemporary and spatially detailed demographic data to reduce this uncertainty, particularly in Africa, and to develop population distribution modeling methods that match the rigor, sophistication, and ability to handle uncertainty of contemporary disease mapping and spread modeling. In the meantime, studies that utilize a particular spatial population dataset need to acknowledge the uncertainties inherent within them and consider how the methods and data that comprise each will affect conclusions.

A high resolution spatial population database of Somalia for disease risk mapping

BACKGROUND

Millions of Somali have been deprived of basic health services due to the unstable political situation of their country. Attempts are being made to reconstruct the health sector, in particular to estimate the extent of infectious disease burden. However, any approach that requires the use of modelled disease rates requires reasonable information on population distribution. In a low-income country such as Somalia, population data are lacking, are of poor quality, or become outdated rapidly. Modelling methods are therefore needed for the production of contemporary and spatially detailed population data.

RESULTS

Here land cover information derived from satellite imagery and existing settlement point datasets were used for the spatial reallocation of populations within census units. We used simple and semi-automated methods that can be implemented with free image processing software to produce an easily updatable gridded population dataset at 100 × 100 meters spatial resolution. The 2010 population dataset was matched to administrative population totals projected by the UN. Comparison tests between the new dataset and existing population datasets revealed important differences in population size distributions, and in population at risk of malaria estimates. These differences are particularly important in more densely populated areas and strongly depend on the settlement data used in the modelling approach.

CONCLUSIONS

The results show that it is possible to produce detailed, contemporary and easily updatable settlement and population distribution datasets of Somalia using existing data. The 2010 population dataset produced is freely available as a product of the AfriPop Project and can be downloaded from: http://www.afripop.org.