- Open Access
Gridded population survey sampling: a systematic scoping review of the field and strategic research agenda
International Journal of Health Geographics volume 19, Article number: 34 (2020)
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.
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.
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.
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.
Household surveys provide insight into the distribution of health, demographics, economics, and behaviours of populations, and are a primary resource for decision-making across low- and middle-income countries (LMICs). Household survey data are used to estimate more than a quarter of the Sustainable Development Goal (SDG) indicators, to generate small area estimates (SAEs) of indicators that support decision-making in decentralized health systems , and to inform the distribution of development funding to, and within, LMICs. Nevertheless, as the use of household surveys has increased over the last 40 years, data accuracy has likely decayed because survey methods have not changed while population characteristics and behaviours have—drastically.
Survey sampling methods have been mature for decades . The Demographic and Health Surveys (DHS) , Multiple Indicator Cluster Surveys (MICS) , and Living Standards Measurement Surveys (LSMS)  have collectively supported hundreds of multi-topic surveys in over 130 countries since 1980 using essentially the same methods. They follow a stratified two-stage cluster design in which first- or second-level administrative units (e.g. provinces) serve as strata. In stage one, census enumeration areas (EAs) are selected with probability proportionate to population size (PPS), and then a field-based mapping-listing activity is conducted in each selected cluster to fully list all households. In stage two, households are sampled from the full listing by an impartial central team, and interviewers return to selected households to administer questionnaires. Rapid needs assessments and public opinion surveys follow a similar design, but tend to use a faster, less-rigorous household selection protocol during stage two; rather than performing a full household listing, interviewers perform a random walk from a central point in the cluster and directly sample households in the field [6, 7]. This approach is considered less rigorous than a full listing because interviewers may consciously or sub-consciously avoid undesirable households; the protocol can result in a “main street” bias, and information needed to adjust for household sample probabilities and non-response are, generally, not collected .
The last 40 years have seen dramatic increases in mobility of LMIC populations, urbanisation, and socioeconomic disparities within cities . The urban poorest include climate and political refugees, seasonal migrants, and rural migrants, as well as multi-generation slum dwellers, street-sleepers, and marginalized minorities . Concurrently, availability of technologies (e.g., mobile phones) and new data (e.g., high-resolution satellite imagery) has rapidly increased, though few new technologies and datasets have been incorporated into standard survey practice. This mismatch has resulted in challenges to sample frame and field protocol accuracy [10, 11]. Furthermore, the SDGs have increased emphasis on disaggregated indicators , raising concerns about whether current survey designs are ideal for accurate SAEs, which we highlight below. To address these emerging issues, survey practitioners have begun to use modelled gridded population datasets as an alternative to census sample frames.
Gridded population datasets are estimates of the total population in small grid cells derived with a geo-statistical model using census or small area population counts and a number of other spatial datasets . The cells in gridded population estimates range in size from 30 × 30 m to 1 × 1 km, and many of these datasets are free and publicly available. In gridded population sampling, grid cells are often aggregated into clusters of a desired population size, and used in place of census EAs. To contextualise gridded population sampling, we provide further background on key reasons that teams have turned to gridded population sampling, and provide an overview of gridded population datasets.
The objectives of this paper are to provide a systematic scoping review of the datasets, tools, and methods used in existing gridded population surveys in LMICs, and outline a research agenda that would equip survey designers to decide when gridded population sampling can be viable and preferable to census-based sampling. We aim to encourage new research and practices that improve the accuracy of survey data and, ultimately, to improve accuracy of health and other household survey data to better target resources toward mobile and vulnerable populations.
Reasons for use of gridded population sampling
The main reason that survey practitioners have turned to gridded population sampling is lack of a current, accurate census sample frame. One in four LMICs has not had a census in the last 10 years . High rates of urban growth and mobility in LMICs mean that megacities in Asia, and soon Africa, grow by 1000 people per day . Since 2000, the average household survey sample frame in LMICs was 7 years old, with some surveys using 15 (Pakistan) and 30 (DR Congo) year old sample frames . Vulnerable populations are most likely to be excluded from surveys with an outdated sample frame because population growth is greater among lower-income households, and they are more likely to be undercounted in censuses .
The second reason for choosing gridded population sampling is that standard survey methods, largely developed for rural settings 40 years ago, struggle to sample mobile and vulnerable households accurately . Even if the census sample frame is complete and updated, a time gap between the household mapping-listing activity and interviews in DHS, MICS, LSMS, and similar surveys means that mobile and vulnerable households are more likely to be counted as non-responders or to be under-listed during survey fieldwork. Furthermore, the mappers-listers who are responsible for generating the final household sample in the DHS, MICS, LSMS, and similar surveys frame have short interactions (e.g. 5–15 min) with residents. With limited rapport, residents may be unwilling to describe informal households in the dwelling (living space, e.g. apartment), and/or the mapping-listing team assumes one household occupies each dwelling which is simply not the case in modern LMIC cities [16, 17]. In LMICs that do not have geocoded census EA boundaries, mapping-listing activities rely on hand-sketched paper maps and subjective descriptions of EA boundaries by local leaders, leading to further potential biases.
A third reason for choosing gridded population sampling is to produce improved small area estimates. In recent years, funders and decision-makers have pushed for important health outcomes to be measured at smaller administrative scales (e.g., district) for policy planning and evaluation [1, 12]. Increased availability of satellite imagery has enabled survey outcomes to be modelled at fine-scale using geostatistical SAE techniques . However, SAEs based on the stratified two-stage PPS design tend to have large uncertainty in sparsely-sampled rural areas and in heterogeneous urban settings [19,20,21]. Gridded population estimates can provide more up-to-date and detailed population counts than outdated census frames, permit new survey designs such as area-microcensus sampling to eliminate the time lag between mapping-listing and interviews, and facilitate spatial oversampling to improve survey-based SAEs.
Gridded population data
A number of gridded population datasets are available across LMICs (Table 1). “Top-down” datasets disaggregate census counts to grid cells, while “bottom-up” estimates are based on micro-census population counts . Currently, nine sources of “top-down” estimates are available for multiple LMICs, and two sources of “bottom-up” estimates are in production for multiple LMICs .
Top-down gridded population estimates
Nearly all gridded population datasets available at the time of this writing were derived from “top-down” models which disaggregate census or other full-coverage population counts into small grid cells. These models produce “pycnophylactic” estimates such that the cell-level counts re-aggregate to the counts of input administrative data . Generally input population counts are adjusted to UN population projections before modelling , however, this still means that countries with the greatest need for improved sample frames have the least accurate top-down gridded population datasets. Additional factors influence the accuracy of top-down modelled population estimates, namely the aggregation scale of the input census data, modelling approach, and area of the output grid cell.
Scale of input data
The most important factor for top-down gridded population accuracy is the aggregation scale of the model input population data (e.g., census) . This is intuitive; the more detailed and accurate the input dataset, the more precise and certain the output estimates will be in small grid squares.
The simplest top-down models assume that the population is spread evenly across grid cells within administrative units (e.g. GPWv4 [26, 27]) or are weighted by land cover types (GHS-POP [28, 29]; HRSL ; ESRI WPE ; WorldPop-Land Cover [32, 33]). These modelling techniques are more mechanical than statistical, and thus do not result in estimates of model error. These models produce reasonably accurate cell-level estimates if a highly accurate dataset of built-up areas is used to mask unpopulated areas, and the input population data is both disaggregated and recent , all of which are rare in LMICs.
Complex modelling techniques using multiple Earth Observation-, government-, and crowd-sourced spatial covariates (e.g., WorldPop-Random Forest [34, 35], WorldPop-Global [34, 35], LandScan-Global , Demobase ) are employed to produce substantially more accurate gridded population estimates. WorldPop-Random Forest and WorldPop-Global are 100 × 100 m datasets of the residential (night-time) population based on a regression tree machine-learning method, and are accompanied by prediction errors at the scale of the input population data . Neither WorldPop-Random Forest nor WorldPop-Global datasets mask built-up areas, thus they produce small, non-zero population predictions in deserts, savannahs, and forests (e.g., 0.0001 persons per cell). WorldPop-Global incorporates changes to urban extents over time, and is modelled from a reduced set of covariates that are available globally. Demobase is a free 100 × 100 m dataset of the residential (night-time) population in three countries based on semi-automated classification of high- and medium-resolution satellite imagery, with prediction errors at the scale of the input population data . LandScan-Global is an annual 1 × 1 km dataset of the “ambient” population; a 24-h average of daytime commuter population and night-time residential population . This dataset is derived with a smart interpolation approach and model error estimates are not provided .
A common issue across all top-down gridded population datasets is that they sometimes allocate population to airports, universities, factories, and government buildings, affecting cell-level accuracy in urban areas. This misallocation may be reduced by including covariates associated with variation in urban density (e.g. building footprints), and/or covariates that represent points of interest and infrastructure where people tend not to live.
Area of output grid cells
The geographic size of the output cells influences estimated population accuracy at the cell-level. Generally, estimates in smaller cells have greater uncertainty, and accuracy improves with cell size. For household survey sampling, however, cell-level accuracy must be balanced against feasibility of cell size for fieldwork; in dense urban contexts, a 100 × 100 m grid cell might contain 1000 s of people. Gridded population datasets with small cells are easy to aggregate into larger units, however, complex methods are required by users to disaggregate cells that are too populous for survey field work .
Bottom-up gridded population datasets
To generate gridded population estimates in countries without a recent or accurate census, “bottom-up” models are under development to estimate population counts based on recent micro-census samples rather than full censuses . These models draw on geo-statistical relationships between population density in a micro-census unit and settlement type, as well as other spatial covariates to predict population counts in un-sampled areas of the country. These census-independent gridded population estimates are produced by the GRID3 and LandScan-HD projects for multiple LMICs, and have the benefit of being constrained to settled areas [39, 40]. Other projects have resulted in a bottom-up gridded population estimate for a single country (e.g. Sierra Leone , Afghanistan ).
Gridded population sample frame attributes
Gridded population datasets are not provided with urban/rural classes, administrative unit names, or estimates of sub-populations because they are designed to be aggregated into any desired spatial unit. Publicly available datasets can be used to classify a gridded population dataset within a geographic information system (GIS) (e.g., ArcGIS, QGIS) or statistical program (e.g., R, Python). Urban/rural datasets include the Global Urban Footprint (GUF)  dataset of 85 × 85 m grid cells classified as built-up or not built-up, and the Global Human Settlement GHS-SMOD  dataset of 1 × 1 km grid cells classified as high-dense urban, low-dense urban, rural, and unsettled based on the GHS-POP population density and GHS-BUILT-UP datasets. Administrative boundaries are available as shapefiles through a number of initiatives including GADM , UN-SALB , and MapLibrary .
We conducted a systematic scoping review in Scopus using the terms: (“gridded” OR “landscan” OR “worldpop” OR “gpw” OR “ghs-pop” OR “hrsl” OR “wpe” OR “demobase”) AND (“population” OR “household”) AND “survey”. No limits were placed on the search (e.g. year or status of publication). Article abstracts were independently screened by co-authors DRT and DAR and retained if they referred to sampling of human populations. We additionally solicited reports, websites, and articles from colleagues. DRT performed a full-text review of all screened articles and reports, and retained those that described a method, tool, or survey based on gridded population data. Retained publications were reviewed for gridded population survey details including sample frame, sample design, sample size, target population, tools, and protocols used. This review followed PRISMA-ScR guidelines (see Additional file 1). A strategic gridded population survey research agenda was iteratively developed among co-authors with feedback from survey experts in a 2 day workshop and via email.
The review in Scopus identified 65 articles describing a gridded population survey, tool, or method. Solicitation of documents from colleagues resulted in seven additional publicly available resources, and awareness of five additional survey teams who described to us their unpublished gridded population surveys (Fig. 1). Although we did not restrict our search by geography, nearly all identified gridded population surveys were located in LMICs and were motivated by an outdated or unavailable census. This literature review resulted in 43 gridded population surveys across 29 LMICs: Bangladesh [16, 47, 48], Brazil , Burkina Faso, Cameroon, Colombia , Cote D’Ivoire, DR Congo [49, 50], Ghana , Guatemala , India , Indonesia, Iraq , Kenya , Mali, Mozambique , Myanmar , Nepal [16, 17, 47, 50], Niger, Nigeria , Rwanda , Somalia [53, 54], Tanzania, Thailand , Togo, Uganda , Uruguay, and Vietnam [16, 47] (Table 2). Additional gridded population surveys were conducted in Greece, Italy, and Slovenia, but excluded from this analysis (personal communication, S. Nichols, Gallup, 14 Jan 2020). Three resources described tools or methods for selecting gridded population survey clusters [38, 55, 56].
Most sample frames in early gridded population surveys were derived from LandScan-Global 1 × 1 km estimates [48, 49, 51, 52], while most recent surveys derived sample frames from WorldPop 100 × 100 m estimates (Table 2) [16, 17, 47, 50, 53]. The final selection of households followed two approaches. First, all eligible households in a cluster or segment were interviewed (called area-microcensus hereafter). Second, households were sampled within clusters or segments before interviewing (called two-stage hereafter). We note whether household sampling was conducted with a robust probability method (i.e., complete mapping-listing of households before sampling households), or a non-probability method (e.g., random-walk or spin-the-pen) .
Thirty-two of the 43 surveys (74%) had national coverage with 1000–4000 households each (Table 2). Nineteen surveys (47%) followed an area-microcensus design (Table 2) for one of four reasons. First, area-microcensus sampling saved time and costs by eliminating, or reducing, the mapping-listing activity [47, 50]. Second, it restricted fieldwork to one visit in insecure or hard-to-reach areas [50, 51]. Third, it provided a simple field protocol and required less training of interviewers which was assumed to ensure higher data quality . Fourth, in complex, dynamic urban environments, it removed the time lag between mapping-listing and interviewing, guarding against under-listing of mobile or vulnerable households, and placed responsibility for household identification with interviewers rather than mapper-listers [17, 47].
One survey compared area-microcensus and two-stage gridded population sampling in Kathmandu, Nepal, and found that when interviewers (area-microcensus) rather than the mapper-listers (two-stage) performed the household listing, non-family and single-adult households were more likely to be identified because interviewers spent substantially more time building rapport with residents in area-microcensus clusters during the interview process . This study also found lower design effects for socio-economic indicators in the area-microcensus design, suggesting better identification of heterogeneous “hidden” households, though household response rates were also lower in the area-microcensus sample .
Four tools and numerous ad-hoc geographic information system (GIS) approaches were described to select gridded population survey clusters (Table 3), and resulted in various forms of a gridded population sample frame, visualized in Fig. 2. The first gridded population sampling tool was the open-source GridSample R package, released by Thomson and colleagues in 2016  and used in six sub-national surveys [17, 47, 50]. The GridSample R algorithm treats the gridded population dataset as the sample frame and selects grid cells with PPS allowing for stratification, oversampling in urban/rural domains, and spatial oversampling . The algorithm runs on a personal computer and is limited by the computer’s memory. All datasets must be pre-processed and specified by the user, allowing use of any gridded population but also requiring GIS and/or R programming skills. The algorithm enables optional “growth” of clusters to a minimum population size or maximum area by randomly adding neighbouring cells after selection of “seed” cells with PPS. While this process results in clusters with roughly consistent population size for improved fieldwork, the population counts in the “grown” clusters do not reflect the population counts used for sample selection, and may skew sample weights . The output is a shapefile of cluster boundaries, with attributes of estimated population counts.
Second, the Geo-sampling survey tool was created by RTI and used in 14 national and sub-national surveys  (personal communication, J. Cajka, RTI, 9 Apr 2020). The Geo-sampling tool is designed for use with large grid cells (e.g. 1 × 1 km), and supports a multi-stage stratified sampling approach. Clients are provided with a shapefile of the final cluster boundaries and population counts. In 13 surveys conducted in 2014–15, administrative units were sampled with PPS, and then 1 × 1 km LandScan-Global cells were sampled with PPS. To improve fieldwork, 1 × 1 km cells with fewer than 250 persons were excluded, potentially biasing the sample toward higher-density populations. The sampled 1 × 1 km cells were partitioned into 150, 100 or 50 m grid cells depending on population density. Next, a deep-learning residential scene classification model was used to identify and exclude small cells without settlement, and disaggregate the 1 × 1 km population to remaining small cells. Finally, three of the small cells were selected at random for an area-microcensus sample . In a 2019 RTI survey, WorldPop-Random Forest estimates were aggregated to 400 × 400 m cells and used in place of 1 × 1 km cells, and a machine-learning building feature extraction algorithm was used to sample structures in the final stage of sampling (personal communication, J. Cajka, RTI, 9 Apr 2020).
Third, many gridded population surveys have developed ad-hoc approaches to sampling using GIS software, such as ArcGIS. Galway and colleagues sampled 1 × 1 km cells with PPS, then randomly selected one household in one building and performed a random walk . Thomson and colleagues converted 1 × 1 km population counts to random points, selected points at random, manually delineated clusters within cells around selected points, and performed an area-microcensus sample . Muñoz and Langeraar proposed an approach for 1 × 1 km cells, though it is unclear if a survey followed . In this approach, 1 × 1 km cells are aggregated to 3 × 3 km grid cells and sampled with PPS. Then 1 × 1 km grid cells are combined within selected 3 × 3 km cells to achieve a minimum population and sampled with PPS. Next, they select a 1 × 1 km (or larger) area and manually delineate segments of approximately 100 households each. One segment is randomly selected, households are listed via a mapping-listing activity, and finally a sample of households is selected . Sollom and colleagues joined 1 × 1 km gridded population estimates to rural village point locations and sampled points with PPS, and then used spin-the-pen to sample households in the field . Qader and colleagues used gridded population estimates to update census EA counts in urban areas where EA boundaries were available, and used a quadtree method to create different sized grid cells with approximately the same population each in rural areas . The combined frame was sampled with PPS before manually segmenting and randomly selecting one household per segment . Finally, Gallup polling teams aggregated 100 × 100 m WorldPop-Random Forest grid cells into larger units (e.g. 200 or 500 m cells) depending on local population density, sampled aggregated grid cells with PPS, and used satellite imagery to choose a central location from which to start a random walk (personal communication, S. Nichols, Gallup, 14 Jan 2020).
Fourth, GridSample.org is a free web-based tool released in late 2019 that runs the open-source GridSample2.0 algorithm developed by Flowminder Foundation. It provides a point-and-click interface, preloaded datasets, and guidance to enter parameters and select clusters for a gridded population survey. It also leverages gridEZ, a publicly-available algorithm, to group cells into clusters before sampling. Preloaded datasets include WorldPop-Global 100 × 100 m gridded population estimates, a bespoke version of WorldPop-Global 100 × 100 m estimates constrained to settled areas, GADM administrative boundaries, and GHS-SMOD urban/rural boundaries. All surveys are implicitly stratified by level of urbanicity; stratification and spatial oversampling are supported; and custom coverage, strata, or sample frame boundaries can be uploaded by users. GridSample.org is designed for low-bandwidth settings, running sample selection remotely on a super-computer. The user is emailed a shapefile of cluster boundaries, population estimates to calculate sample weights, and a report. The US Department of State (USDS) (personal communication, S. Staveteig Ford and M. Kirwin of USDS, 10 Apr 2020), ORB International (personal communication, C. Carter and Y. Dudaronak of ORB International, 9 Apr 2020), and the Nigerian Government (personal communication, R. Bhattarai of Flowminder Foundation and M. Imohi of the Nigeria National Bureau of Statistics, 10 Dec 2019) used GridSample.org to select national or state-level household surveys. USDS in Uruguay and the Nigerian Government used GridSample.org and WorldPop-Global 100 × 100 m grid cells to create a sample frame of “medium” gridEZ units (clusters) of approximately 500 people each, while ORB International used the tool to define “large” gridEZ units (clusters) of up to 1,200 people in a maximum area of 5 × 5 km in seven Sahel countries characterized by vast unsettled areas and low-density population. The USDS and ORB International survey used a random-walk method to sample households in the field, while the Nigerian Government performed a full listing in sampled clusters before sampling and interviewing households.
A range of simple-to-advanced tools have been used to implement gridded population surveys. Lower-tech field tools include use of paper maps displaying cluster boundaries over satellite imagery in Google Earth, and paper listing forms and questionnaires [49,50,51]. Higher-tech field tools include tablet-based applications for navigation [16, 48], paper field maps designed in GIS [16, 17, 50, 51, 53], and tablet-based household listing and/or questionnaires [7, 16, 17, 48, 50]. Satellite imagery was essential to all gridded populations surveys to manually segment along roads, rivers, and other features [47, 49, 56], and as a field map base layer [48,49,50,51, 53]. In some surveys, satellite imagery was used to digitize building footprints and roads in OpenStreetMap which was then displayed as a field map base layer [17, 47]. Many teams included points of interest from OpenStreetMap or GPS coordinates of recognizable intersections/structures on field maps to aid navigation [17, 47, 49, 53].
The successful implementation of more than 40 gridded population sample surveys across a variety of settings bodes well for this emerging field. Due to the use of English language search terms in this review, focus on academic literature, and use of gridded population sampling by some practitioners who do not publicly describe their survey methods, the number of gridded population surveys implemented is likely larger. The possibility that gridded population sampling might improve accuracy of data about vulnerable and mobile populations, especially in settings with outdated or inaccurate census data, is appealing to researchers and practitioners who work on health and social inequities in LMICs . However, a survey statistician considering whether to recommend an outdated census-based frame or a gridded population frame is faced with questions about sample frame accuracy, methods to form and select sample frame units, and optimal survey designs. Next, we outline a research agenda to equip survey designers to identify situations where gridded population sampling can be a feasible and trustworthy option. The agenda shows key stages of a gridded population survey and available options (Fig. 3).
Choose gridded population
Top-down gridded population datasets that restrict estimates to settled areas (e.g. LandScan-Global) are likely to underestimate rural, and overestimate urban, populations because small settlements are often undetected in the settlement layer. Conversely, datasets that estimate population in all landmasses (e.g. WorldPop-Global) likely overestimate rural, and underestimate urban, population because fractions of the population are allocated to unsettled cells. Factors that affect survey accuracy include the gridded population model accuracy, aggregation of the gridded population model input dataset, whether residential or ambient population is modelled, accuracy and type of covariates, and area of the cell in which population is estimated .
A major gap is that cell-level accuracy is not known for any top-down gridded population datasets. To assess accuracy, a recent census disaggregated to household locations would be needed, though this is rarely, if ever, available. The next best option is comparison of modelled gridded population estimates with micro-census counts from a sample of areas. Household listings from a recent geo-located household survey aggregated to cells might serve this purpose, but to our knowledge, data sharing agreements for such work have not been investigated or defined. Simulated household-level datasets are a third option .
Furthermore, survey designers will want to consider how uncertainty estimates might be used to improve sample designs or sample size calculations. Presently, some top-down datasets (e.g. WorldPop-Global, Demobase) include model prediction errors at the scale of the input population dataset based on internal validation, and new bottom-up datasets include cell-level uncertainty estimates. A clear understanding of cell-level accuracy is not only important to assess whether gridded population datasets are technically fit for purpose in practical applications that effects the public’s health and wellbeing ; this transparency is also a critical component of fostering political buy-in . DHS, MICS, LSMS, and other surveys are distributed via national statistical offices, and thus their sample frames are often mandated to come from official sources. Processes are needed for national statistical agencies to engage with gridded population dataset production so that official endorsements might be made .
Choose sample design
Area-microcensus sample designs in small clusters (e.g. 10–20 households) may prove to be faster and cheaper than two-stage designs in larger clusters (e.g. 100–300 households), and more accurately sample vulnerable urban populations; however, there can be a counter-balancing detriment of higher survey design effects due to variable numbers of respondents per cluster, greater within-cluster homogeneity, and lower response rates. For survey designers to assess these trade-offs and to select a sample size that will meet stakeholders’ goals for budget, timeline, and statistical precision, they need reliable projections of likely design effects in area-microcensus samples. The current limited evidence is mixed. A simulation study of a rural population in Namibia found that nearly twice as many area-microcensus clusters would be needed to achieve the same precision as a two-stage survey, holding constant the number of respondents per cluster . While a study in urban Nepal found higher design effects for demographic indicators and lower design effects for socio-economic indicators in an area-microcensus design versus a two-stage design .
Also, as urban settlement classification becomes increasingly possible , survey designers need to understand how within-urban stratification affects the various sample designs used in gridded population, and other, surveys. With no way to stratify urban populations, all surveys are at risk of under-sampling or omitting slums and other vulnerable populations [61, 62]. This threat to survey accuracy and social equity will only grow as LMIC urban population continue to expand in the coming decades. In addition, research is needed to balance survey designs that can support both precise design-based estimation of outcomes and precise SAEs of indicators at fine geographic scales to support local decision-making, SDGs, and other initiatives requiring spatially disaggregated data .
Create sample frame
Existing gridded population sample frame approaches result in squared-off, arbitrary cluster boundaries that are not recognizable on the ground. Improved methods are needed to use natural features such as rivers and roads to delineate cluster boundaries from gridded population data. To date, nearly all spatial feature datasets for LMICs have been produced by governments or volunteers (e.g. OpenStreetMap), neither of which are sufficiently detailed, complete, or spatially precise to support delineation of “natural” cluster boundaries across many LMICs, especially those with vast sparsely populated areas . However, this is rapidly changing with new availability of very high resolution imagery and supercomputing facilities (e.g. Maxar’s building footprint and road data in 51 African countries) which might lead to new approaches to delineating “natural” cluster boundaries for gridded population data . As the field continues to evolve, survey designers need to be confident that clusters will yield the right number of eligible respondents and have a geographic area that can be canvassed by a field team in the time budgeted for fieldwork.
Several gridded population sampling tools and approaches are available, and their feasibility is influenced by cost, transparency of the methods, clarity of documentation, and usability by survey design professionals in government agencies and organizations who may not have advanced programming and GIS skills. The GridSample R algorithm does not scale to large geographic areas nor does include an optimal method to create clusters from grid cells, and is thus not suitable for routine national surveys. GridSample.org is free, offers ease of use and clear documentation, but currently cannot be adapted for in-house (private) use by national statistical agencies without manipulation of the underlying GridSample2.0 algorithm. Use of the Geo-sampling tool requires the hiring and support of an external company, which prohibits widespread use.
The emerging field of gridded population survey sampling should recommend tools and protocols for both lower- and higher-tech settings. For example, a common protocol should be described to deal with arbitrary gridded population boundaries that intersect buildings (e.g. include buildings in north and east boundaries, exclude buildings on south and west boundaries). Uniquely, gridded population surveys rely on access to up-to-date high-resolution satellite imagery (0.5 m) for fieldwork. This is less of a challenge in urban areas worldwide thanks to Google Earth, Bing, and other free websites. However, imagery resolution in rural areas of LMICs is quite variable, with images sometimes being several years old. As a result, it would be difficult to implement gridded population surveys in areas of heavy forest or cloud cover. Furthermore, tools for implementing surveys (e.g. Survey123, OpenMapKit) tend to focus on questionnaires and often lack integration with satellite imagery, visualisation of cluster boundaries, and geo-location services in offline environments which means that multiple tools are often needed to conduct gridded survey field activities .
Organizations with skills in GIS and digital tools can successfully implement surveys with gridded population sample frames, which have the potential to yield samples that are more representative of mobile and vulnerable respondents than outdated census-based frames. However, census-based frames are likely to be considered a safe choice by many survey designers because censuses have long been the standard and their limitations are commonly accepted. To recommend a gridded population frame would involve risks and rewards that are currently difficult to quantify. New tools are needed to evaluate gridded population datasets and frames in specific country contexts, and to facilitate low-burden survey implementation. There are opportunities to develop tools for nearly every stage of survey planning and implementation, which ultimately will improve the accuracy of survey data.
Availability of data and materials
Demographic and Health Surveys
Geographic Information System
Global Positioning System
Low- or middle-income country
Living Standard Measurement Surveys
Multiple Indicator Cluster Surveys
Probability proportionate to (population) size
Small area estimate
Sustainable Development Goal
US Department of State
Victora CG, Black RE, Boerma JT, Bryce J. Measuring impact in the Millennium Development Goal era and beyond: a new approach to large-scale effectiveness evaluations. Lancet. 2011. https://doi.org/10.1016/S0140-6736(10)60810-0.
Hancioglu A, Arnold F. Measuring coverage in MNCH: tracking progress in health for women and children Using DHS and MICS household surveys. PLoS Med. 2013. https://doi.org/10.1371/journal.pmed.1001391.
ICF International. Demographic and Health Survey sampling and household listing manual. 2012. https://dhsprogram.com/pubs/pdf/DHSM4/DHS6_Sampling_Manual_Sept2012_DHSM4.pdf. Accessed 12 Apr 2020.
United Nations Children’s Fund. MICS 5 tools: manual for mapping and household listing. 2013. https://mics.unicef.org/tools?round=mics5. Accessed 12 Apr 2020.
Grosh MGE, Munoz J. A manual for planning and implementing the Living Standards Measurement Study survey. LSMS Working Paper No. 126. 1996. https://documents.worldbank.org/curated/en/363321467990016291/pdf/multi-page.pdf. Accessed 12 Apr 2020.
Pew Research Center. Methods: International Survey Research. 2020. https://www.pewresearch.org/methods/international-survey-research/. Accessed 12 Apr 2020.
Gallup. World Poll methodology. 2017 https://news.gallup.com/poll/105226/world-poll-methodology.aspx. Accessed 12 Apr 2020.
Grais RF, Rose AMC, Guthmann J. Don’t spin the pen: two alternative methods for second-stage sampling in urban cluster surveys. Emerg Themes Epidemiol. 2007. https://doi.org/10.1186/1742-7622-4-8.
United Nations Human Settlements Programme. Urbanization and development: emerging futures. 2016. https://unhabitat.org/world-cities-report. Accessed 12 Apr 2020.
Cobham A. Uncounted: power, inequalities and the post-2015 data revolution. Development. 2014. https://doi.org/10.1057/dev.2015.28.
Lucci P, Bhatkal T, Khan A. Are we underestimating urban poverty? World Dev. 2018. https://doi.org/10.1016/j.worlddev.2017.10.022.
Golding N, Burstein R, Longbottom J, Browne AJ, Fullman N, Osgood-Zimmerman A, et al. Mapping under-5 and neonatal mortality in Africa, 2000–15: a baseline analysis for the Sustainable Development Goals. Lancet. 2017. https://doi.org/10.1016/S0140-6736(17)31758-0.
Leyk S, Gaughan AE, Adamo SB, de Sherbinin A, Balk D, Freire S, et al. Allocating people to pixels: a review of large-scale gridded population data products and their fitness for use. Earth Syst Sci Data Discuss. 2019. https://doi.org/10.5194/essd-11-1385-2019.
United Nations Statistics Division. 2020 world population and housing census programme. Census dates for all countries. 2019. https://unstats.un.org/unsd/demographic-social/census/censusdates/. Accessed 12 Apr 2020.
United Nations Department of Economic and Social Affairs. World Urbanization Prospects: the 2018 revision. 2018. https://population.un.org/wup/DataQuery/. Accessed 12 Apr 2020.
Thomson DR, Bhattarai R, Khanal S, Manandhar S, Dhungel R, Gajurel S, et al. Addressing unintentional exclusion of vulnerable and mobile households in traditional surveys in Kathmandu, Dhaka and Hanoi: a mixed methods feasibility study. 2020. https://www.preprints.org/manuscript/201910.0320/v3. Accessed 21 Aug 2020.
Elsey H, Thomson D, Lin R, Maharjan U, Agarwal S, Newell J. Addressing inequities in urban health: Do decision-makers have the data they need? Report from the Urban Health Data special session at International Conference on Urban Health Dhaka 2015. J Urban Heal. 2016. https://doi.org/10.1007/s11524-016-0046-9.
Burgert-Brucker CR, Dontamsetti T, Gething PW. The DHS program’s modeled surfaces spatial datasets. Stud Fam Plann. 2018. https://doi.org/10.1111/sifp.12050.
Cutts FT, Dansereau E, Ferrari MJ, Hanson M, McCarthy KA, Metcalf CJE, et al. Using models to shape measles control and elimination strategies in low- and middle-income countries: a review of recent applications. Vaccine. 2020. https://doi.org/10.1016/j.vaccine.2019.11.020.
Mosser JF, Gagne-Maynard W, Rao PC, Osgood-Zimmerman A, Fullman N, Graetz N, et al. Mapping diphtheria-pertussis-tetanus vaccine coverage in Africa, 2000–2016: a spatial and temporal modelling study. Lancet. 2019. https://doi.org/10.1016/S0140-6736(19)30226-0.
Utazi CE, Thorley J, Alegana VA, Ferrari MJ, Takahashi S, Metcalf CJE, et al. High resolution age-structured mapping of childhood vaccination coverage in low and middle income countries. Vaccine. 2018. https://doi.org/10.1016/j.vaccine.2018.02.020.
Wardrop NA, Jochem WC, Bird TJ, Chamberlain HR, Clarke D, Kerr D, et al. Spatially disaggregated population estimates in the absence of national population and housing census data. Proc Natl Acad Sci. 2018. https://doi.org/10.1073/pnas.1715305115.
Tobler WR. Smooth Pycnophylactic Interpolation for Geographical Regions. J Am Stat Assoc. 1979. https://doi.org/10.2307/2286968.
United Nations Department of Economic and Social Affairs. World Population Prospects. 2019. https://population.un.org/wpp/Download/Standard/Population/. Accessed 12 Apr 2020.
Hay S, Noor A, Nelson A, Tatem A. The accuracy of human population maps for public health application. Trop Med Int Heal. 2005. https://doi.org/10.1111/j.1365-3156.2005.01487.x.
Center for International Earth Science Information Network—CIESIN—Columbia University. Gridded Population of the World, Version 4 (GPWv4), New York. 2016. https://sedac.ciesin.columbia.edu/data/collection/gpw-v4. Accessed 12 Apr 2020.
Doxsey-Whitfield E, MacManus K, Adamo SB, Pistolesi L, Squires J, Borkovska O, et al. Taking advantage of the improved availability of census data: a first look at the Gridded Population of the World, Version 4. Pap Appl Geogr. 2015. https://doi.org/10.1080/23754931.2015.1014272.
European Commission. Global human settlement city model (GHS-SMOD), Ispra. 2017. https://ghsl.jrc.ec.europa.eu/ghs_smod2019.php. Accessed 12 Apr 2020
Pesaresi M, Ehrlich D, Florczyk AJ, Freire S, Julea A, Kemper T, et al. Operating procedure for the production of the Global Human Settlement Layer from Landsat data of the epochs 1975, 1990, 2000, and 2014. 2016. https://publications.jrc.ec.europa.eu/repository/handle/JRC97705. Accessed 12 Apr 2020.
Facebook Connectivity Lab and Center for International Earth Science Information Network—CIESIN—Columbia University. High Resolution Settlement Layer (HRSL). 2016. https://ciesin.columbia.edu/data/hrsl/. Accessed 12 Apr 2020.
Frye C, Nordstrand E, Wright DJ, Terborgh C, Foust J. Using classified and unclassified land cover data to estimate the footprint of human settlement. Data Sci J. 2018. https://doi.org/10.5334/dsj-2018-020/.
Linard C, Gilbert M, Tatem AJ. Assessing the use of global land cover data for guiding large area population distribution modelling. GeoJournal. 2011. https://doi.org/10.1007/s10708-010-9364-8.
Gaughan AE, Stevens FR, Linard C, Jia P, Tatem AJ. High resolution population distribution maps for Southeast Asia in, 2010 and 2015. PLoS ONE. 2013. https://doi.org/10.1371/journal.pone.0055882.
Stevens FR, Gaughan AE, Linard C, Tatem AJ. Disaggregating census data for population mapping using random forests with remotely-sensed and ancillary data. PLoS ONE. 2015. https://doi.org/10.1371/journal.pone.0107042.
WorldPop. WorldPop-Global data 2000–2020, Southampton. 2019. www.worldpop.org/geodata/listing?id=29. Accessed 12 Apr 2020.
Dobson JE, Brlght EA, Coleman PR, Worley BA, Bright EA, Coleman PR, et al. LandScan: a global population database for estimating populations at risk. Photogramm Eng Remote Sensing. 2000;66(7):849–57.
Azar D, Engstrom R, Graesser J, Comenetz J. Generation of fine-scale population layers using multi-resolution satellite imagery and geospatial data. Remote Sens Environ. 2013. https://doi.org/10.1016/j.rse.2012.11.022.
Chew RF, Amer S, Jones K, Unangst J, Cajka J, Allpress J, et al. Residential scene classification for gridded population sampling in developing countries using deep convolutional neural networks on satellite imagery. Int J Health Geogr. 2018. https://doi.org/10.1186/s12942-018-0132-1.
Oak Ridge National Laboratories. LandScan HD. 2020. https://landscan.ornl.gov. Accessed 17 Aug 2020.
Geo-Referenced Infrastructure and Demographic Data for Development (GRID3). High-resolution population estimates. 2020. https://grid3.org/resources/data. Accessed 17 Aug 2020.
Hillson R, Alejandre JD, Jacobsen KH, Ansumana R, Bockarie AS, Bangura U, et al. Methods for determining the uncertainty of population estimates derived from satellite imagery and limited survey data: a case study of Bo City, Sierra Leone. PLoS ONE. 2014. https://doi.org/10.1371/journal.pone.0112241.
Juran S, Tatem A. New methodology: a hybrid census to generate spatially-disaggregated population estimates. 2018. www.unfpa.org/sites/default/files/resource-pdf/Hybrid_Census_Brief_v9.pdf. Accessed 12 Apr 2020.
DLR Earth Observation Center. Global Urban Footprint (GUF), Webling. 2017. www.dlr.de/eoc/en/desktopdefault.aspx/tabid-11725/20508_read-47944/. Accessed 12 Apr 2020.
GADM. Global administrative areas version 3.6. 2015. www.gadm.org. Accessed 12 Apr 2020.
United Nations Geospatial Information Section. United Nations Second Administrative Level Boundaries (UN-SALB), New York. 2019. www.unsalb.org/data. Accessed 12 Apr 2020.
Map Maker Ltd. Map Library, Argyll. 2007. www.maplibrary.org. Accessed 12 Apr 2020.
Elsey H, Poudel AN, Ensor T, Mirzoev T, Newell JN, Hicks JP, et al. Improving household surveys and use of data to address health inequities in three Asian cities: protocol for the Surveys for Urban Equity (SUE) mixed methods and feasibility study. BMJ Open. 2018. https://doi.org/10.1136/bmjopen-2018-024182.
Cajka J, Amer S, Ridenhour J, Allpress J. Geo-sampling in developing nations. Int J Soc Res Methodol. 2018. https://doi.org/10.1080/13645579.2018.1484989.
Thomson DR, Hadley MB, Greenough PG, Castro MC. Modelling strategic interventions in a population with a total fertility rate of 8.3: a cross-sectional study of Idjwi Island, DRC. BMC Public Health. 2012. https://doi.org/10.1186/1471-2458-12-959.
GridSample. Case Studies. 2019. https://gridsample.org/case-studies. Accessed 12 Apr 2020.
Galway LP, Bell N, Sae AS, Hagopian A, Burnham G, Flaxman A, et al. A two-stage cluster sampling method using gridded population data, a GIS, and Google Earth imagery in a population-based mortality survey in Iraq. Int J Health Geogr. 2012. https://doi.org/10.1186/1476-072X-11-12.
Sollom R, Richards AK, Parmar P, Mullany LC, Lian SB, Iacopino V, et al. Health and human rights in Chin State, Western Burma: a population-based assessment using multistaged household cluster sampling. PLoS Med. 2011. https://doi.org/10.1371/journal.pmed.1001007.
Qader SH, Lefebvre V, Ninneman A, Himelein K, Pape U, Bengtsson L, et al. A novel approach to the automatic designation of pre-defined census enumeration areas and population sampling frames: a case study in Somalia. Policy Research Working Paper No. 8972. 2019. https://documents.worldbank.org/curated/en/810771565268760366/pdf/A-Novel-Approach-to-the-Automatic-Designation-of-Predefined-Census-Enumeration-Areas-and-Population-Sampling-Frames-A-Case-Study-in-Somalia.pdf. Accessed 12 Apr 2020.
Pape U, Wollburg P. Estimation of poverty in Somalia using innovative methodologies. Policy Research Working Paper No. 8735. 2019. https://documents.worldbank.org/curated/en/509221549985694077/pdf/WPS8735.pdf. Accessed 12 Apr 2020.
Thomson DR, Stevens FR, Ruktanonchai NW, Tatem AJ, Castro MC. GridSample: an R package to generate household survey primary sampling units (PSUs) from gridded population data. Int J Health Geogr. 2017. https://doi.org/10.1186/s12942-017-0098-4.
Muñoz J, Langeraar W. A census-independent sampling strategy for a household survey in Myanmar. 2013. https://winegis.com/images/census-independent-GIS-based-sampling-strategy-for-household-surveys-plan-of-action%2520removed.pdf. Accessed 20 Apr 2020.
Thomson DR, Kools L, Jochem WC. Linking synthetic populations to household geolocations: a demonstration in Namibia. Data. 2018. https://doi.org/10.3390/data3030030.
Prakash M, Ramage S, Kavvada A, Goodman S. Open Earth observations for sustainable urban development. Remote Sens. 2020. https://doi.org/10.3390/rs12101646.
Kools L. One-stage versus two-stage cluster sampling, a simulation study. Essays on wealth, health, and data collection [Dissertation]. Leiden University; 2018. https://openaccess.leidenuniv.nl/handle/1887/67120. Accessed 12 Apr 2020.
Mahabir R, Croitoru A, Crooks A, Agouris P, Stefanidis A. A critical review of high and very high-resolution remote sensing approaches for detecting and mapping slums: trends, challenges and emerging opportunities. Urban Sci. 2018. https://doi.org/10.3390/urbansci2010008.
Lilford RJ, Oyebode O, Satterthwaite D, Chen Y, Mberu B, Watson SI, et al. Improving the health and welfare of people who live in slums. Lancet. 2017. https://doi.org/10.1016/S0140-6736(16)31848-7.
Thomson DR, Kuffer M, Boo G, Hati B, Grippa T, Elsey H, et al. Need for an integrated deprived area “slum” mapping system (IDEAMAPS) in low- and middle-income countries (LMICs). Soc Sci. 2020. https://doi.org/10.3390/socsci9050080.
Sinha P, Gaughan AE, Stevens FR, Nieves JJ, Sorichetta A, Tatem AJ. Assessing the spatial sensitivity of a random forest model: application in gridded population modeling. Comput Environ Urban Syst. 2019. https://doi.org/10.1016/j.compenvurbsys.2019.01.006.
Minaei M. Evolution, density and completeness of OpenStreetMap road networks in developing countries: the case of Iran. Appl Geogr. 2020. https://doi.org/10.1016/j.apgeog.2020.102246.
Many thanks to Madeline Fabric, Clara Burgert-Brucker, and Sarchil Qader for their review and valuable feedback on an early draft of this paper. Also thanks to the participants at the Gates Foundation 2018 Gridded Population Sampling meeting for discussions and synthesis of gridded population survey evidence: Vincent Seaman, Brian Blankespoor, Hil Lyons, Julien Nobili, Kathleen Wannemuehler, Laina Mercer, Laura A Dwyer-Lindgren, Linnea Zimmerman, Michael Wild, Benjamin L. Nygren, Sabhanaz Rashid Diya, Tove Ryman, Eric M Weber, Wenfeng Gong, Win Brown, Sophia Magalona, Fiona Lau, and Jon Wakefield.
Dana R Thomson received funding from the UK Economic and Social Research Council (ESRC) grant ID: ES/5500161/1. The funder had no role in the study design, execution, analyses, interpretation of the data, or decision to submit results.
Ethics approval and consent to participate
Consent for publication
DRT is the creator of the GridSample R and GridSample2.0 python algorithms, and formerly managed GridSample.org at Flowminder Foundation. AJT is the director of the WorldPop team at University of Southampton. DAR and MCC declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Thomson, D.R., Rhoda, D.A., Tatem, A.J. et al. Gridded population survey sampling: a systematic scoping review of the field and strategic research agenda. Int J Health Geogr 19, 34 (2020). https://doi.org/10.1186/s12942-020-00230-4
- Survey design
- Household survey