Large-scale spatial population databases in infectious disease research

Population-based denominators matter: Bias in U.S. Virgin Islands COVID-19 vaccination coverage under changing population counts

PURPOSE

The U.S. Virgin Islands (USVI) receives an updated population count once every 10 years and used 2010 decennial census population counts to estimate COVID-19 vaccination coverage during the COVID-19 emergency response. We investigated whether using outdated (2010) or modeled (2020 international database [IDB]) population counts biased vaccination coverage estimates used to inform public health priorities during the 2020-2022 COVID-19 response.

METHODS

We estimated percentage of USVI residents with a completed primary COVID-19 vaccination series during December 16, 2020-September 20, 2022. Vaccination coverage was calculated as number of persons who completed the vaccination series divided by 2010 and 2020 decennial census population counts and 2020 IDB intercensal estimate.

RESULTS

COVID-19 vaccination coverage using the 2020 population count was 12 % points higher than coverage using 2010 denominator (2010 denominator: 51 %; 2020 denominator: 63 %). Vaccination coverage estimated using 2020 IDB was approximately equal with the 2010 decennial census estimate (52 %).

CONCLUSIONS

Using 2010 and modeled population counts underestimated 2020 USVI COVID-19 vaccination coverage given the 18 % population decline during 2010-2020, potentially limiting USVI's ability to assess vaccination progress. Identifying mechanisms for more reliable population enumeration or improved estimate modeling are essential for accurately guiding USVI public health decision-making.

A 31-year (1990-2020) global gridded population dataset generated by cluster analysis and statistical learning

Continuously monitoring global population spatial dynamics is crucial for implementing effective policies related to sustainable development, including epidemiology, urban planning, and global inequality. However, existing global gridded population data products lack consistent population estimates, making them unsuitable for time-series analysis. To address this issue, this study designed a data fusion framework based on cluster analysis and statistical learning approaches, which led to the generation of a continuous global gridded population dataset (GlobPOP). The GlobPOP dataset was evaluated through two-tier spatial and temporal validation to demonstrate its accuracy and applicability. The spatial validation results show that the GlobPOP dataset is highly accurate. The temporal validation results also reveal that the GlobPOP dataset performs consistently well across eight representative countries and cities despite their unique population dynamics. With the availability of GlobPOP datasets in both population count and population density formats, researchers and policymakers can leverage the new dataset to conduct time-series analysis of the population and explore the spatial patterns of population development at global, national, and city levels.

A data-driven framework to assess population dynamics during novel coronavirus outbreaks: A case study on Xiamen Island, China

The outbreak of the Coronavirus Disease 2019 (COVID-19) has profoundly influenced daily life, necessitating the understanding of the relationship between the epidemic's progression and population dynamics. In this study, we present a data-driven framework that integrates GIS-based data mining technology and a Susceptible, Exposed, Infected and Recovered (SEIR) model. This approach helps delineate population dynamics at the grid and community scales and analyze the impacts of government policies, urban functional areas, and intercity flows on population dynamics during the pandemic. Xiamen Island was selected as a case study to validate the effectiveness of the data-driven framework. The results of the high/low cluster analysis provide 99% certainty (P < 0.01) that the population distribution between January 23 and March 16, 2020, was not random, a phenomenon referred to as high-value clustering. The SEIR model predicts that a ten-day delay in implementing a lockdown policy during an epidemic can lead to a significant increase in the number of individuals infected by the virus. Throughout the epidemic prevention and control period (January 23 to February 21, 2020), residential and transportation areas housed more residents. After the resumption of regular activities, the population was mainly concentrated in residential, industrial, and transportation, as well as road facility areas. Notably, the migration patterns into and out of Xiamen were primarily centered on neighboring cities both before and after the outbreak. However, migration indices from cities outside the affected province drastically decreased and approached zero following the COVID-19 outbreak. Our findings offer new insights into the interplay between the epidemic's development and population dynamics, which enhances the prevention and control of the coronavirus epidemic.

Evaluating the Ability to Use Contextual Features Derived from Multi-Scale Satellite Imagery to Map Spatial Patterns of Urban Attributes and Population Distributions

With an increasing global population, accurate and timely population counts are essential for urban planning and disaster management. Previous research using contextual features, using mainly very-high-spatial-resolution imagery (<2 m spatial resolution) at subnational to city scales, has found strong correlations with population and poverty. Contextual features can be defined as the statistical quantification of edge patterns, pixel groups, gaps, textures, and the raw spectral signatures calculated over groups of pixels or neighborhoods. While they correlated with population and poverty, which components of the human-modified landscape were captured by the contextual features have not been investigated. Additionally, previous research has focused on more costly, less frequently acquired very-high-spatial-resolution imagery. Therefore, contextual features from both very-high-spatial-resolution imagery and lower-spatial-resolution Sentinel-2 (10 m pixels) imagery in Sri Lanka, Belize, and Accra, Ghana were calculated, and those outputs were correlated with OpenStreetMap building and road metrics. These relationships were compared to determine what components of the human-modified landscape the features capture, and how spatial resolution and location impact the predictive power of these relationships. The results suggest that contextual features can map urban attributes well, with out-of-sample R2 values up to 93%. Moreover, the degradation of spatial resolution did not significantly reduce the results, and for some urban attributes, the results actually improved. Based on these results, the ability of the lower resolution Sentinel-2 data to predict the population density of the smallest census units available was then assessed. The findings indicate that Sentinel-2 contextual features explained up to 84% of the out-of-sample variation for population density.

Temporal Geospatial Analysis of COVID-19 Pre-Infection Determinants of Risk in South Carolina

Disparities and their geospatial patterns exist in morbidity and mortality of COVID-19 patients. When it comes to the infection rate, there is a dearth of research with respect to the disparity structure, its geospatial characteristics, and the pre-infection determinants of risk (PIDRs). This work aimed to assess the temporal-geospatial associations between PIDRs and COVID-19 infection at the county level in South Carolina. We used the spatial error model (SEM), spatial lag model (SLM), and conditional autoregressive model (CAR) as global models and the geographically weighted regression model (GWR) as a local model. The data were retrieved from multiple sources including USAFacts, U.S. Census Bureau, and the Population Estimates Program. The percentage of males and the unemployed population were positively associated with geodistributions of COVID-19 infection (p values < 0.05) in global models throughout the time. The percentage of the white population and the obesity rate showed divergent spatial correlations at different times of the pandemic. GWR models fit better than global models, suggesting nonstationary correlations between a region and its neighbors. Characterized by temporal-geospatial patterns, disparities in COVID-19 infection rate and their PIDRs are different from the mortality and morbidity of COVID-19 patients. Our findings suggest the importance of prioritizing different populations and developing tailored interventions at different times of the pandemic.

COVID-19 Deaths in the United States: Shifts in Hot Spots over the Three Phases of the Pandemic and the Spatiotemporally Varying Impact of Pandemic Vulnerability

The geographic areas most impacted by COVID-19 may not remain static because public health measures/behaviors change dynamically, and the impacts of pandemic vulnerability also may vary geographically and temporally. The nature of the pandemic makes spatiotemporal methods essential to understanding the distribution of COVID-19 deaths and developing interventions. This study examines the spatiotemporal trends in COVID-19 death rates in the United States from March 2020 to May 2021 by performing an emerging hot spot analysis (EHSA). It then investigates the effects of the COVID-19 time-dependent and basic social vulnerability factors on COVID-19 death rates using geographically and temporally weighted regression (GTWR). The EHSA results demonstrate that over the three phases of the pandemic (first wave, second wave, and post-vaccine deployment), hot spots have shifted from densely populated cities and the states with a high percentage of socially vulnerable individuals to the states with relatively relaxed social distancing requirements, and then to the states with low vaccination rates. The GTWR results suggest that local infection and testing rates, social distancing interventions, and other social, environmental, and health risk factors show significant associations with COVID-19 death rates, but these associations vary over time and space. These findings can inform public health planning.

Impact of an accelerated melting of Greenland on malaria distribution over Africa

Studies about the impact of future climate change on diseases have mostly focused on standard Representative Concentration Pathway climate change scenarios. These scenarios do not account for the non-linear dynamics of the climate system. A rapid ice-sheet melting could occur, impacting climate and consequently societies. Here, we investigate the additional impact of a rapid ice-sheet melting of Greenland on climate and malaria transmission in Africa using several malaria models driven by Institute Pierre Simon Laplace climate simulations. Results reveal that our melting scenario could moderate the simulated increase in malaria risk over East Africa, due to cooling and drying effects, cause a largest decrease in malaria transmission risk over West Africa and drive malaria emergence in southern Africa associated with a significant southward shift of the African rain-belt. We argue that the effect of such ice-sheet melting should be investigated further in future public health and agriculture climate change risk assessments.

Understanding dynamics of population flood exposure in Canada with multiple high-resolution population datasets

In recent years, geospatial data (e.g. remote sensing imagery), and other relevant ancillary datasets (e.g. land use land cover, climate conditions) have been utilized through sophisticated algorithms to produce global population datasets. With a handful of such datasets, their performances and skill in flood exposure assessment have not been explored. This study proposes a comprehensive framework to understand the dynamics and differences in population flood exposure over Canada by employing four global population datasets alongside the census data from Statistics Canada as the reference. The flood exposure is quantified based on a set of floodplain maps (for 2015, 1 in 100-yr and 1 in 200-yr event) for Canada derived from the CaMa-Flood global flood model. To obtain further insights at the regional level, the methodology is implemented over six flood-prone River Basins in Canada. We find that about 9% (3.31 million) and 11% (3.90 million) of the Canadian population resides within 1 in 100-yr and 1 in 200-yr floodplains. We notice an excellent performance of WorldPop, and LandScan in most of the cases, which is unaffected by the representation of flood hazard, while Global Human Settlement and Gridded Population of the World showed large deviations. At last, we determined the long-term dynamics of population flood exposure and vulnerability from 2006 to 2019. Through this analysis, we also identify the regions that contain a significantly larger population exposed to floods. The relevant conclusions derived from the study highlight the need for careful selection of population datasets for preventing further amplification of uncertainties in flood risk. We recommend a detailed assessment of the severely exposed regions by including precise ground-level information. The results derived from this study may be useful not only for flood risk management but also contribute to understanding other disaster impacts on human-environment interrelationships.

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.

Generalized Vertical Components of built-up areas from global Digital Elevation Models by multi-scale linear regression modelling

The estimation of the vertical components of built-up areas from free Digital Elevation Model (DEM) global data filtered by multi-scale convolutional, morphological and textural transforms are generalized at the spatial resolution of 250 meters using linear least-squares regression techniques. Six test cases were selected: Hong Kong, London, New York, San Francisco, Sao Paulo, and Toronto. Five global DEM and two DEM composites are evaluated in terms of 60 combinations of linear, morphological and textural filtering and different generalization techniques. Four generalized vertical components estimates of built-up areas are introduced: the Average Gross Building Height (AGBH), the Average Net Building Height (ANBH), the Standard Deviation of Gross Building Height (SGBH), and the Standard Deviation of Net Building Height (SNBH). The study shows that the best estimation of the net GVC of built-up areas given by the ANBH and SNBH, always contains a greater error than their corresponding gross GVC estimation given by the AGBH and SGBH, both in terms of mean and standard deviation. Among the sources evaluated in this study, the best DEM source for estimating the GVC of built-up areas with univariate linear regression techniques is a composite of the 1-arcsec Shuttle Radar Topography Mission (SRTM30) and the Advanced Land Observing Satellite (ALOS) World 3D–30 m (AW3D30) using the union operator (CMP_SRTM30-AW3D30_U). A multivariate linear model was developed using 16 satellite features extracted from the CMP_SRTM30-AW3D30_U enriched by other land cover sources, to estimate the gross GVC. A RMSE of 2.40 m and 3.25 m was obtained for the AGBH and the SGBH, respectively. A similar multivariate linear model was developed to estimate the net GVC. A RMSE of 6.63 m and 4.38 m was obtained for the ANBH and the SNBH, respectively. The main limiting factors on the use of the available global DEMs for estimating the GVC of built-up areas are two. First, the horizontal resolution of these sources (circa 30 and 90 meters) corresponds to a sampling size that is larger than the expected average horizontal size of built-up structures as detected from nadir-angle Earth Observation (EO) data, producing more reliable estimates for gross vertical components than for net vertical component of built-up areas. Second, post-production processing targeting Digital Terrain Model specifications may purposely filter out the information on the vertical component of built-up areas that are contained in the global DEMs. Under the limitations of the study presented here, these results show a potential for using global DEM sources in order to derive statistically generalized parameters describing the vertical characteristics of built-up areas, at the scale of 250x250 meters. However, estimates need to be evaluated in terms of the specific requirements of target applications such as spatial population modelling, urban morphology, climate studies and so on.

Hotspots and correlates of soil-transmitted helminth infections in a Venezuelan rural community: Which are the "wormy" houses?

OBJECTIVES

To determine spatial clustering and risk factors for occurrence and intensity of infection for soil-transmitted helminthiasis (STH), namely Ascaris lumbricoides, Trichuris trichiura, Strongyloides stercoralis and hookworms in a Venezuelan rural community.

METHODS

MIF-fixed faecal samples were individually collected for STH testing. The Getis-Ord statistic was used to determine significant STH clustering within 25/50/100 m radiuses around houses. Individual- and house-level factors associated with STH occurrence and intensity of infection were determined using generalized estimating equations.

RESULTS

Significant clusters of "wormy" houses for one or multiple parasites were found at distances of 25-50 m around 13 houses. Risk factors differed between occurrence and intensity of infection. Overcrowding in the house increased occurrence of S. stercoralis, T. trichiura and hookworm infections, while poor housing conditions increased A. lumbricoides infection risk. Overcrowding, poor faecal disposal system, economic dependency and lack of basic services differentially influenced the STHs. The "wormy" houses were mainly those built with waste materials, under economic dependency and lacking indoor water supply.

CONCLUSIONS

STH distribution in a community is clustered, with significant hotspots of STH occurrence and intensity of infection and different associated risk factors. Targeting the "wormy" houses is expected to affect STH morbidity more efficiently.

Space-Time Conditional Autoregressive Modeling to Estimate Neighborhood-Level Risks for Dengue Fever in Cali, Colombia

Vector-borne diseases affect more than 1 billion people a year worldwide, causing more than 1 million deaths, and cost hundreds of billions of dollars in societal costs. Mosquitoes are the most common vectors responsible for transmitting a variety of arboviruses. Dengue fever (DENF) has been responsible for nearly 400 million infections annually. Dengue fever is primarily transmitted by female Aedes aegypti and Aedes albopictus mosquitoes. Because both Aedes species are peri-domestic and container-breeding mosquitoes, dengue surveillance should begin at the local level—where a variety of local factors may increase the risk of transmission. Dengue has been endemic in Colombia for decades and is notably hyperendemic in the city of Cali. For this study, we use weekly cases of DENF in Cali, Colombia, from 2015 to 2016 and develop space–time conditional autoregressive models to quantify how DENF risk is influenced by socioeconomic, environmental, and accessibility risk factors, and lagged weather variables. Our models identify high-risk neighborhoods for DENF throughout Cali. Statistical inference is drawn under Bayesian paradigm using Markov chain Monte Carlo techniques. The results provide detailed insight about the spatial heterogeneity of DENF risk and the associated risk factors (such as weather, proximity to Aedes habitats, and socioeconomic classification) at a fine level, informing public health officials to motivate at-risk neighborhoods to take an active role in vector surveillance and control, and improving educational and surveillance resources throughout the city of Cali.

Evaluating the Effects of Climate and Environmental Factors on Under-5 Children Malaria Spatial Distribution Using Generalized Additive Models (GAMs)

Although malaria burden has declined globally following scale up of intervention, the disease has remained a leading cause of hospitalization and deaths among children aged under-5 years in Nigeria. Malaria is known to be related to climate and environmental conditions. Previous research has usually studied the effects of these factors, neglecting possible correlation between them, high correlation among variables is a source of multicollinearity that induces overfitting in regression modelling. In this paper, a factor analysis was first introduced to circumvent the issue of multicollinearity and a Generalized Additive Model (GAM) was subsequently explored to identify the important risk factors that might influence the prevalence of childhood malaria in Nigeria. The GAM incorporated the complexity of the survey data, while simultaneously modelling the nonlinear and spatial random effects to allow a more precise identification of the major malaria risk factors that influence the geographical distribution of the disease. From our findings, the three latent factor components (constituted by humidity, precipitation, potential evapotranspiration, and wet days/maximum and minimum temperature/proximity to permanent waters, respectively) were significantly associated with malaria prevalence. Our analysis also detected statistically significant and nonlinear effect of altitude: the risk of malaria increased with lower values but declined sharply with higher values. A significant spatial variability in under-5 malaria prevalence across the survey clusters was also observed; malaria burden was higher in the northern part of Nigeria. Investigating the impact of important risk factors and geographical location on childhood malaria is of high relevance for the sustainable development goals (SDGs) 2015–2030 Agenda on malaria eradication, and we believe that the information obtained from this study and the generated risk maps can be useful to effectively target intervention efforts to high-risk areas based on climate and environmental context.

Information differences across spatial resolutions and scales for disease surveillance and analysis: The case of Visceral Leishmaniasis in Brazil

Nationwide disease surveillance at a high spatial resolution is desired for many infectious diseases, including Visceral Leishmaniasis. Statistical and mathematical models using data collected from surveillance activities often use a spatial resolution and scale either constrained by data availability or chosen arbitrarily. Sensitivity of model results to the choice of spatial resolution and scale is not, however, frequently evaluated. This study aims to determine if the choice of spatial resolution and scale are likely to impact statistical and mathematical analyses. Visceral Leishmaniasis in Brazil is used as a case study. Probabilistic characteristics of disease incidence, representing a likely outcome in a model, are compared across spatial resolutions and scales. Best fitting distributions were fit to annual incidence from 2004 to 2014 by municipality and by state. Best fits were defined as the distribution family and parameterization minimizing the sum of absolute error, evaluated through a simulated annealing algorithm. Gamma and Poisson distributions provided best fits for incidence, both among individual states and nationwide. Comparisons of distributions using Kullback-Leibler divergence shows that incidence by state and by municipality do not follow distributions that provide equivalent information. Few states with Gamma distributed incidence follow a distribution closely resembling that for national incidence. These results demonstrate empirically how choice of spatial resolution and scale can impact mathematical and statistical models.

An Introductory Framework for Choosing Spatiotemporal Analytical Tools in Population-Level Eco-Epidemiological Research

Spatiotemporal visualization and analytical tools (SATs) are increasingly being applied to risk-based surveillance/monitoring of adverse health events affecting humans, animals, and ecosystems. Different disciplines use diverse SATs to address similar research questions. The juxtaposition of these diverse techniques provides a list of options for researchers who are new to population-level spatial eco-epidemiology. Here, we are conducting a narrative review to provide an overview of the multiple available SATs, and introducing a framework for choosing among them when addressing common research questions across disciplines. The framework is comprised of three stages: (a) pre-hypothesis testing stage, in which hypotheses regarding the spatial dependence of events are generated; (b) primary hypothesis testing stage, in which the existence of spatial dependence and patterns are tested; and (c) secondary-hypothesis testing and spatial modeling stage, in which predictions and inferences were made based on the identified spatial dependences and associated covariates. In this step-wise process, six key research questions are formulated, and the answers to those questions should lead researchers to select one or more methods from four broad categories of SATs: (T1) visualization and descriptive analysis; (T2) spatial/spatiotemporal dependence and pattern recognition; (T3) spatial smoothing and interpolation; and (T4) geographic correlation studies (i.e., spatial modeling and regression). The SATs described here include both those used for decades and also other relatively new tools. Through this framework review, we intend to facilitate the choice among available SATs and promote their interdisciplinary use to support improving human, animal, and ecosystem health.

Routine data for malaria morbidity estimation in Africa: challenges and prospects

BACKGROUND

The burden of malaria in sub-Saharan Africa remains challenging to measure relying on epidemiological modelling to evaluate the impact of investments and providing an in-depth analysis of progress and trends in malaria response globally. In malaria-endemic countries of Africa, there is increasing use of routine surveillance data to define national strategic targets, estimate malaria case burdens and measure control progress to identify financing priorities. Existing research focuses mainly on the strengths of these data with less emphasis on existing challenges and opportunities presented.

CONCLUSION

Here we define the current imperfections common to routine malaria morbidity data at national levels and offer prospects into their future use to reflect changing disease burdens.

Geospatial Disaggregation of Population Data in Supporting SDG Assessments: A Case Study from Deqing County, China

Quantitative assessments and dynamic monitoring of indicators based on fine-scale population data are necessary to support the implementation of the United Nations (UN) 2030 Agenda and to comprehensively achieve its 17 Sustainable Development Goals (SDGs). However, most population data are collected by administrative units, and it is difficult to reflect true distribution and uniformity in space. To solve this problem, based on fine building information, a geospatial disaggregation method of population data for supporting SDG assessments is presented in this paper. First, Deqing County in China, which was divided into residential areas and nonresidential areas according to the idea of dasymetric mapping, was selected as the study area. Then, the town administrative areas were taken as control units, building area and number of floors were used as weighting factors to establish the disaggregation model, and population data with a resolution of 30 m in Deqing County in 2016 were obtained. After analyzing the statistical population of 160 villages and the disaggregation results, we found that the global average accuracy was 87.08%. Finally, by using the disaggregation population data, indicators 3.8.1, 4.a.1, and 9.1.1 were selected to conduct an accessibility analysis and a buffer analysis in a quantitative assessment of the SDGs. The results showed that the SDG measurement and assessment results based on the disaggregated population data were more accurate and effective than the results obtained using the traditional method.

An Improved Global Analysis of Population Distribution in Proximity to Active Volcanoes, 1975–2015

Better and more detailed analyses of global human exposure to hazards and associated disaster risk require improved geoinformation on population distribution and densities. In particular, issues of temporal and spatial resolution are important for determining the capacity for assessing changes in these distributions. We combine the best-available global population grids with latest data on volcanoes, to assess and characterize the worldwide distribution of population from 1975–2015 in relation to recent volcanism. Both Holocene volcanoes and those where there is evidence of significant eruptions are considered. A comparative analysis is conducted for the volcanic hot spots of Southeast Asia and Central America. Results indicate that more than 8% of the world’s 2015 population lived within 100 km of a volcano with at least one significant eruption, and more than 1 billion people (14.3%) lived within 100 km of a Holocene volcano, with human concentrations in this zone increasing since 1975 above the global population growth rate. While overall spatial patterns of population density have been relatively stable in time, their variation with distance is not monotonic, with a higher concentration of people between 10 and 20 km from volcanoes. We find that in last 40 years in Southeast Asia the highest population growth rates have occurred in close proximity to volcanoes (within 10 km), whereas in Central America these are observed farther away (beyond 50 km), especially after 1990 and for Holocene volcanoes.

How accurate are modelled birth and pregnancy estimates? Comparison of four models using high resolution maternal health census data in southern Mozambique

Background

Existence of inequalities in quality and access to healthcare services at subnational levels has been identified despite a decline in maternal and perinatal mortality rates at national levels, leading to the need to investigate such conditions using geographical analysis. The need to assess the accuracy of global demographic distribution datasets at all subnational levels arises from the current emphasis on subnational monitoring of maternal and perinatal health progress, by the new targets stated in the Sustainable Development Goals.

Methods

The analysis involved comparison of four models generated using Worldpop methods, incorporating region-specific input data, as measured through the Community Level Intervention for Pre-eclampsia (CLIP) project. Normalised root mean square error was used to determine and compare the models’ prediction errors at different administrative unit levels.

Results

The models’ prediction errors are lower at higher administrative unit levels. All datasets showed the same pattern for both the live birth and pregnancy estimates. The effect of improving spatial resolution and accuracy of input data was more prominent at higher administrative unit levels.

Conclusion

The validation successfully highlighted the impact of spatial resolution and accuracy of maternal and perinatal health data in modelling estimates of pregnancies and live births. There is a need for more data collection techniques that conduct comprehensive censuses like the CLIP project. It is also imperative for such projects to take advantage of the power of mapping tools at their disposal to fill the gaps in the availability of datasets for populated areas.

Population allocation at the housing unit level: estimates around underground natural gas storage wells in PA, OH, NY, WV, MI, and CA

BACKGROUND

Spatially accurate population data are critical for determining health impacts from many known risk factors. However, the utility of the increasing spatial resolution of disease mapping and environmental exposures is limited by the lack of receptor population data at similar sub-census block spatial scales.

METHODS

Here we apply an innovative method (Population Allocation by Occupied Domicile Estimation - ABODE) to disaggregate U.S. Census populations by allocating an average person per household to geospatially-identified residential housing units (RHU). We considered two possible sources of RHU location data: address point locations and building footprint centroids. We compared the performance of ABODE with the common proportional population allocation (PPA) method for estimating the nighttime residential populations within 200 m radii and setback areas (100 - 300 ft) around active underground natural gas storage (UGS) wells (n = 9834) in six U.S. states.

RESULTS

Address location data generally outperformed building footprint data in predicting total counts of census residential housing units, with correlations ranging from 0.67 to 0.81 at the census block level. Using residentially-sited addresses only, ABODE estimated upwards of 20,000 physical households with between 48,126 and 53,250 people living within 200 m of active UGS wells - likely encompassing the size of a proposed UGS Wellhead Safety Zone. Across the 9834 active wells assessed, ABODE estimated between 5074 and 10,198 more people living in these areas compare to PPA, and the difference was significant at the individual well level (p = < 0.0001). By either population estimation method, OH exhibits a substantial degree of hyperlocal land use conflict between populations and UGS wells - more so than other states assessed. In some rare cases, population estimates differed by more than 100 people for the small 200 m² well-areas. ABODE's explicit accounting of physical households confirmed over 50% of PPA predictions as false positives indicated by non-zero predictions in areas absent physical RHUs.

CONCLUSIONS

Compared to PPA - in allocating identical population data at sub-census block spatial scales -ABODE provides a more precise population at risk (PAR) estimate with higher confidence estimates of populations at greatest risk. 65% of UGS wells occupy residential urban and suburban areas indicating the unique land use conflicts presented by UGS systems that likely continue to experience population encroachment. Overall, ABODE confirms tens of thousands of homes and residents are likely located within the proposed UGS Wellhead Safety Zone - and in some cases within state's oil and gas well surface setback distances - of active UGS wells.

Gridded Population Maps Informed by Different Built Settlement Products

The spatial distribution of humans on the earth is critical knowledge that informs many disciplines and is available in a spatially explicit manner through gridded population techniques. While many approaches exist to produce specialized gridded population maps, little has been done to explore how remotely sensed, built-area datasets might be used to dasymetrically constrain these estimates. This study presents the effectiveness of three different high-resolution built area datasets for producing gridded population estimates through the dasymetric disaggregation of census counts in Haiti, Malawi, Madagascar, Nepal, Rwanda, and Thailand. Modeling techniques include a binary dasymetric redistribution, a random forest with a dasymetric component, and a hybrid of the previous two. The relative merits of these approaches and the data are discussed with regards to studying human populations and related spatially explicit phenomena. Results showed that the accuracy of random forest and hybrid models was comparable in five of six countries.

Unlocking pandemic potential: prevalence and spatial patterns of key substitutions in avian influenza H5N1 in Egyptian isolates

Background

Avian influenza H5N1 has a high human case fatality rate, but is not yet well-adapted to human hosts. Amino acid substitutions currently circulating in avian populations may enhance viral fitness in, and thus viral adaptation to, human hosts. Substitutions which could increase the risk of a human pandemic (through changes to host specificity, virulence, replication ability, transmissibility, or drug susceptibility) are termed key substitutions (KS). Egypt represents the epicenter of human H5N1 infections, with more confirmed cases than any other country. To date, however, there have not been any spatial analyses of KS in Egypt.

Methods

Using 925 viral samples of H5N1 from Egypt, we aligned protein sequences and scanned for KS. We geocoded isolates using dasymetric mapping, then carried out geospatial hot spot analyses to identify spatial clusters of high KS detection rates. KS prevalence and spatial clusters were evaluated for all detected KS, as well as when stratified by phenotypic consequence.

Results

A total of 39 distinct KS were detected in the wild, including 17 not previously reported in Egypt. KS were detected in 874 samples (94.5%). Detection rates varied by viral protein with most KS observed in the surface hemagglutinin (HA) and neuraminidase (NA) proteins, as well as the interior non-structural 1 (NS1) protein. The most frequently detected KS were associated with increased viral binding to mammalian cells and virulence. Samples with high overall detection rates of KS exhibited statistically significant spatial clustering in two governorates in the northwestern Nile delta, Alexandria and Beheira.

Conclusions

KS provide a possible mechanism by which avian influenza H5N1 could evolve into a pandemic candidate. With numerous KS circulating in Egypt, and non-random spatial clustering of KS detection rates, these findings suggest the need for increased surveillance in these areas.

Spatially disaggregated population estimates in the absence of national population and housing census data

Population numbers at local levels are fundamental data for many applications, including the delivery and planning of services, election preparation, and response to disasters. In resource-poor settings, recent and reliable demographic data at subnational scales can often be lacking. National population and housing census data can be outdated, inaccurate, or missing key groups or areas, while registry data are generally lacking or incomplete. Moreover, at local scales accurate boundary data are often limited, and high rates of migration and urban growth make existing data quickly outdated. Here we review past and ongoing work aimed at producing spatially disaggregated local-scale population estimates, and discuss how new technologies are now enabling robust and cost-effective solutions. Recent advances in the availability of detailed satellite imagery, geopositioning tools for field surveys, statistical methods, and computational power are enabling the development and application of approaches that can estimate population distributions at fine spatial scales across entire countries in the absence of census data. We outline the potential of such approaches as well as their limitations, emphasizing the political and operational hurdles for acceptance and sustainable implementation of new approaches, and the continued importance of traditional sources of national statistical data.

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.

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.

Plasmodium falciparum malaria importation from Africa to China and its mortality: an analysis of driving factors

Plasmodium falciparum malaria importation from Africa to China is rising with increasing Chinese overseas investment and international travel. Identifying networks and drivers of this phenomenon as well as the contributors to high case-fatality rate is a growing public health concern to enable efficient response. From 2011-2015, 8653 P. falciparum cases leading to 98 deaths (11.3 per 1000 cases) were imported from 41 sub-Saharan countries into China, with most cases (91.3%) occurring in labour-related Chinese travellers. Four strongly connected groupings of origin African countries with destination Chinese provinces were identified, and the number of imported cases was significantly associated with the volume of air passengers to China (P = 0.006), parasite prevalence in Africa (P < 0.001), and the amount of official development assistance from China (P < 0.001) with investment in resource extraction having the strongest relationship with parasite importation. Risk factors for deaths from imported cases were related to the capacity of malaria diagnosis and diverse socioeconomic factors. The spatial heterogeneity uncovered, principal drivers explored, and risk factors for mortality found in the rising rates of P. falciparum malaria importation to China can serve to refine malaria elimination strategies and the management of cases, and high risk groups and regions should be targeted.

Mapping global environmental suitability for Zika virus

Zika virus was discovered in Uganda in 1947 and is transmitted by Aedes mosquitoes, which also act as vectors for dengue and chikungunya viruses throughout much of the tropical world. In 2007, an outbreak in the Federated States of Micronesia sparked public health concern. In 2013, the virus began to spread across other parts of Oceania and in 2015, a large outbreak in Latin America began in Brazil. Possible associations with microcephaly and Guillain-Barré syndrome observed in this outbreak have raised concerns about continued global spread of Zika virus, prompting its declaration as a Public Health Emergency of International Concern by the World Health Organization. We conducted species distribution modelling to map environmental suitability for Zika. We show a large portion of tropical and sub-tropical regions globally have suitable environmental conditions with over 2.17 billion people inhabiting these areas.

DOI:http://dx.doi.org/10.7554/eLife.15272.001

Zika virus is transmitted between humans by mosquitoes. The majority of infections cause mild flu-like symptoms, but neurological complications in adults and infants have been found in recent outbreaks.

Although it was discovered in Uganda in 1947, Zika only caused sporadic infections in humans until 2007, when it caused a large outbreak in the Federated States of Micronesia. The virus later spread across Oceania, was first reported in Brazil in 2015 and has since rapidly spread across Latin America. This has led many people to question how far it will continue to spread. There was therefore a need to define the areas where the virus could be transmitted, including the human populations that might be risk in these areas.

Messina et al. have now mapped the areas that provide conditions that are highly suitable for the spread of the Zika virus. These areas occur in many tropical and sub-tropical regions around the globe. The largest areas of risk in the Americas lie in Brazil, Colombia and Venezuela. Although Zika has yet to be reported in the USA, a large portion of the southeast region from Texas through to Florida is highly suitable for transmission. Much of sub-Saharan Africa (where several sporadic cases have been reported since the 1950s) also presents an environment that is highly suitable for the Zika virus. While no cases have yet been reported in India, a large portion of the subcontinent is also suitable for Zika transmission.

Over 2 billion people live in Zika-suitable areas globally, and in the Americas alone, over 5.4 million births occurred in 2015 within such areas. It is important, however, to recognize that not all individuals living in suitable areas will necessarily be exposed to Zika.

We still lack a great deal of basic epidemiological information about Zika. More needs to be known about the species of mosquito that spreads the disease and how the Zika virus interacts with related viruses such as dengue. As such information becomes available and clinical cases become routinely diagnosed, the global evidence base will be strengthened, which will improve the accuracy of future maps.

DOI:http://dx.doi.org/10.7554/eLife.15272.002

Using Random Forest to Improve the Downscaling of Global Livestock Census Data

Large scale, high-resolution global data on farm animal distributions are essential for spatially explicit assessments of the epidemiological, environmental and socio-economic impacts of the livestock sector. This has been the major motivation behind the development of the Gridded Livestock of the World (GLW) database, which has been extensively used since its first publication in 2007. The database relies on a downscaling methodology whereby census counts of animals in sub-national administrative units are redistributed at the level of grid cells as a function of a series of spatial covariates. The recent upgrade of GLW1 to GLW2 involved automating the processing, improvement of input data, and downscaling at a spatial resolution of 1 km per cell (5 km per cell in the earlier version). The underlying statistical methodology, however, remained unchanged. In this paper, we evaluate new methods to downscale census data with a higher accuracy and increased processing efficiency. Two main factors were evaluated, based on sample census datasets of cattle in Africa and chickens in Asia. First, we implemented and evaluated Random Forest models (RF) instead of stratified regressions. Second, we investigated whether models that predicted the number of animals per rural person (per capita) could provide better downscaled estimates than the previous approach that predicted absolute densities (animals per km2). RF models consistently provided better predictions than the stratified regressions for both continents and species. The benefit of per capita over absolute density models varied according to the species and continent. In addition, different technical options were evaluated to reduce the processing time while maintaining their predictive power. Future GLW runs (GLW 3.0) will apply the new RF methodology with optimized modelling options. The potential benefit of per capita models will need to be further investigated with a better distinction between rural and agricultural populations.

Spatiotemporal patterns of population in mainland China, 1990 to 2010

According to UN forecasts, global population will increase to over 8 billion by 2025, with much of this anticipated population growth expected in urban areas. In China, the scale of urbanization has, and continues to be, unprecedented in terms of magnitude and rate of change. Since the late 1970s, the percentage of Chinese living in urban areas increased from ~18% to over 50%. To quantify these patterns spatially we use time-invariant or temporally-explicit data, including census data for 1990, 2000, and 2010 in an ensemble prediction model. Resulting multi-temporal, gridded population datasets are unique in terms of granularity and extent, providing fine-scale (~100 m) patterns of population distribution for mainland China. For consistency purposes, the Tibet Autonomous Region, Taiwan, and the islands in the South China Sea were excluded. The statistical model and considerations for temporally comparable maps are described, along with the resulting datasets. Final, mainland China population maps for 1990, 2000, and 2010 are freely available as products from 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.

Dynamic population mapping using mobile phone data

During the past few decades, technologies such as remote sensing, geographical information systems, and global positioning systems have transformed the way the distribution of human population is studied and modeled in space and time. However, the mapping of populations remains constrained by the logistics of censuses and surveys. Consequently, spatially detailed changes across scales of days, weeks, or months, or even year to year, are difficult to assess and limit the application of human population maps in situations in which timely information is required, such as disasters, conflicts, or epidemics. Mobile phones (MPs) now have an extremely high penetration rate across the globe, and analyzing the spatiotemporal distribution of MP calls geolocated to the tower level may overcome many limitations of census-based approaches, provided that the use of MP data is properly assessed and calibrated. Using datasets of more than 1 billion MP call records from Portugal and France, we show how spatially and temporarily explicit estimations of population densities can be produced at national scales, and how these estimates compare with outputs produced using alternative human population mapping methods. We also demonstrate how maps of human population changes can be produced over multiple timescales while preserving the anonymity of MP users. With similar data being collected every day by MP network providers across the world, the prospect of being able to map contemporary and changing human population distributions over relatively short intervals exists, paving the way for new applications and a near real-time understanding of patterns and processes in human geography.

Visualization and analytics tools for infectious disease epidemiology: a systematic review

BACKGROUND

A myriad of new tools and algorithms have been developed to help public health professionals analyze and visualize the complex data used in infectious disease control. To better understand approaches to meet these users' information needs, we conducted a systematic literature review focused on the landscape of infectious disease visualization tools for public health professionals, with a special emphasis on geographic information systems (GIS), molecular epidemiology, and social network analysis. The objectives of this review are to: (1) identify public health user needs and preferences for infectious disease information visualization tools; (2) identify existing infectious disease information visualization tools and characterize their architecture and features; (3) identify commonalities among approaches applied to different data types; and (4) describe tool usability evaluation efforts and barriers to the adoption of such tools.

METHODS

We identified articles published in English from January 1, 1980 to June 30, 2013 from five bibliographic databases. Articles with a primary focus on infectious disease visualization tools, needs of public health users, or usability of information visualizations were included in the review.

RESULTS

A total of 88 articles met our inclusion criteria. Users were found to have diverse needs, preferences and uses for infectious disease visualization tools, and the existing tools are correspondingly diverse. The architecture of the tools was inconsistently described, and few tools in the review discussed the incorporation of usability studies or plans for dissemination. Many studies identified concerns regarding data sharing, confidentiality and quality. Existing tools offer a range of features and functions that allow users to explore, analyze, and visualize their data, but the tools are often for siloed applications. Commonly cited barriers to widespread adoption included lack of organizational support, access issues, and misconceptions about tool use.

DISCUSSION AND CONCLUSION

As the volume and complexity of infectious disease data increases, public health professionals must synthesize highly disparate data to facilitate communication with the public and inform decisions regarding measures to protect the public's health. Our review identified several themes: consideration of users' needs, preferences, and computer literacy; integration of tools into routine workflow; complications associated with understanding and use of visualizations; and the role of user trust and organizational support in the adoption of these tools. Interoperability also emerged as a prominent theme, highlighting challenges associated with the increasingly collaborative and interdisciplinary nature of infectious disease control and prevention. Future work should address methods for representing uncertainty and missing data to avoid misleading users as well as strategies to minimize cognitive overload.

Zoom in at African country level: potential climate induced changes in areas of suitability for survival of malaria vectors

Background

Predicting anopheles vectors’ population densities and boundary shifts is crucial in preparing for malaria risks and unanticipated outbreaks. Although shifts in the distribution and boundaries of the major malaria vectors (Anopheles gambiae s.s. and An. arabiensis) across Africa have been predicted, quantified areas of absolute change in zone of suitability for their survival have not been defined. In this study, we have quantified areas of absolute change conducive for the establishment and survival of these vectors, per African country, under two climate change scenarios and based on our findings, highlight practical measures for effective malaria control in the face of changing climatic patterns.

Methods

We developed a model using CLIMEX simulation platform to estimate the potential geographical distribution and seasonal abundance of these malaria vectors in relation to climatic factors (temperature, rainfall and relative humidity). The model yielded an eco-climatic index (EI) describing the total favourable geographical locations for the species. The EI values were classified and exported to a GIS package. Using ArcGIS, the EI shape points were clipped to the extent of Africa and then converted to a raster layer using Inverse Distance Weighted (IDW) interpolation method. Generated maps were then transformed into polygon-based geo-referenced data set and their areas computed and expressed in square kilometers (km²).

Results

Five classes of EI were derived indicating the level of survivorship of these malaria vectors. The proportion of areas increasing or decreasing in level of survival of these malaria vectors will be more pronounced in eastern and southern African countries than those in western Africa. Angola, Ethiopia, Kenya, Mozambique, Tanzania, South Africa and Zambia appear most likely to be affected in terms of absolute change of malaria vectors suitability zones under the selected climate change scenarios.

Conclusion

The potential shifts of these malaria vectors have implications for human exposure to malaria, as recrudescence of the disease is likely to be recorded in several new areas and regions. Therefore, the need to develop, compile and share malaria preventive measures, which can be adapted to different climatic scenarios, remains crucial.

Unifying viral genetics and human transportation data to predict the global transmission dynamics of human influenza H3N2

Information on global human movement patterns is central to spatial epidemiological models used to predict the behavior of influenza and other infectious diseases. Yet it remains difficult to test which modes of dispersal drive pathogen spread at various geographic scales using standard epidemiological data alone. Evolutionary analyses of pathogen genome sequences increasingly provide insights into the spatial dynamics of influenza viruses, but to date they have largely neglected the wealth of information on human mobility, mainly because no statistical framework exists within which viral gene sequences and empirical data on host movement can be combined. Here, we address this problem by applying a phylogeographic approach to elucidate the global spread of human influenza subtype H3N2 and assess its ability to predict the spatial spread of human influenza A viruses worldwide. Using a framework that estimates the migration history of human influenza while simultaneously testing and quantifying a range of potential predictive variables of spatial spread, we show that the global dynamics of influenza H3N2 are driven by air passenger flows, whereas at more local scales spread is also determined by processes that correlate with geographic distance. Our analyses further confirm a central role for mainland China and Southeast Asia in maintaining a source population for global influenza diversity. By comparing model output with the known pandemic expansion of H1N1 during 2009, we demonstrate that predictions of influenza spatial spread are most accurate when data on human mobility and viral evolution are integrated. In conclusion, the global dynamics of influenza viruses are best explained by combining human mobility data with the spatial information inherent in sampled viral genomes. The integrated approach introduced here offers great potential for epidemiological surveillance through phylogeographic reconstructions and for improving predictive models of disease control.

What explains the geographic dispersal of emerging pathogens? Reconstructions of evolutionary history from pathogen gene sequences offer qualitative descriptions of spatial spread, but current approaches are poorly equipped to formally test and quantify the contribution of different potential explanatory factors, such as human mobility and demography. Here, we use a novel phylogeographic method to evaluate multiple potential predictors of viral spread in human influenza dynamics. We identify air travel as the predominant driver of global influenza migration, whilst also revealing the contribution of other mobility processes at more local scales. We demonstrate the power of our inter-disciplinary approach by using it to predict the global pandemic expansion of H1N1 influenza in 2009. Our study highlights the importance of integrating evolutionary and ecological information when studying the dynamics of infectious disease.

Mapping the capacities of fixed health facilities to cover people at risk of gambiense human African trypanosomiasis

Background

The emphasis placed on the activities of mobile teams in the detection of gambiense human African trypanosomiasis (HAT) can at times obscure the major role played by fixed health facilities in HAT control and surveillance. The lack of consistent and detailed data on the coverage of passive case-finding and treatment further constrains our ability to appreciate the full contribution of the health system to the control of HAT.

Methods

A survey was made of all fixed health facilities that are active in the control and surveillance of gambiense HAT. Information on their diagnostic and treatment capabilities was collected, reviewed and harmonized. Health facilities were geo-referenced. Time-cost distance analysis was conducted to estimate physical accessibility and the potential coverage of the population at-risk of gambiense HAT.

Results

Information provided by the National Sleeping Sickness Control Programmes revealed the existence of 632 fixed health facilities that are active in the control and surveillance of gambiense HAT in endemic countries having reported cases or having conducted active screening activities during the period 2000-2012. Different types of diagnosis (clinical, serological, parasitological and disease staging) are available from 622 facilities. Treatment with pentamidine for first-stage disease is provided by 495 health facilities, while for second-stage disease various types of treatment are available in 206 health facilities only. Over 80% of the population at-risk for gambiense HAT lives within 5-hour travel of a fixed health facility offering diagnosis and treatment for the disease.

Conclusions

Fixed health facilities have played a crucial role in the diagnosis, treatment and coverage of at-risk-population for gambiense HAT. As the number of reported cases continues to dwindle, their role will become increasingly important for the prospects of disease elimination. Future updates of the database here presented will regularly provide evidence to inform and monitor a rational deployment of control and surveillance efforts. Support to the development and, if successful, the implementation of new control tools (e.g. new diagnostics and new drugs) is crucial, both for strengthening and expanding the existing network of fixed health facilities by improving access to diagnosis and treatment and for securing a sustainable control and surveillance of gambiense HAT.

Mapping population and pathogen movements

For most of human history, populations have been relatively isolated from each other, and only recently has there been extensive contact between peoples, flora and fauna from both old and new worlds. The reach, volume and speed of modern travel are unprecedented, with human mobility increasing in high income countries by over 1000-fold since 1800. This growth is putting people at risk from the emergence of new strains of familiar diseases, and from completely new diseases, while ever more cases of the movement of both disease vectors and the diseases they carry are being seen. Pathogens and their vectors can now move further, faster and in greater numbers than ever before. Equally however, we now have access to the most detailed and comprehensive datasets on human mobility and pathogen distributions ever assembled, in order to combat these threats. This short review paper provides an overview of these datasets, with a particular focus on low income regions, and covers briefly approaches used to combine them to help us understand and control some of the negative effects of population and pathogen movements.

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.

Using global datasets to create environmental profiles for data-poor regions: a case from the Irrawaddy and Salween River Basins

The increasing availability of spatial data inspires the exploration of previously less-studied, yet regionally and nationally important areas, such as the Irrawaddy and Salween River Basins in Southeast Asia. This article documents our experience using global datasets to create environmental basin profiles in these two basins. Our approach draws on the concepts of freshwater vulnerability assessments that guided the selection of indicators. Data on land use, population distribution and fertilizer load were used. The unit of analysis was chosen to distinguish areas with similar bio-geographical characteristics, such as the critical delta areas. Results were further discussed for sub-areas that experience relatively the most pressure in terms of examined indicators within the studied area. The river mouths of both rivers had the most intensive land use and high population density. They are also home to important ecosystems and are sensitive to changes in upstream areas. Our study presents a concise and spatially distributed view of the environmental basin profiles of the Irrawaddy and Salween River Basins. The analysis also provides some interesting methodological insights about the potential of public macro-scale datasets for environmental assessment. The spatial approach allowed the analysis of different indicators, providing a platform for data integration as well as a visually powerful overview of the study area. Yet, the use of macro-scale datasets entails challenges. Despite improvements, the assessment process tends to be driven by the availability and quality of data, rather than by the actual research and management needs. The greatest utility of macro-scale datasets lies-at least in data-poor areas-in larger scale comparative analyses between the basins and their different sub-areas.

A geographic analysis of population density thresholds in the influenza pandemic of 1918-19

BACKGROUND

Geographic variables play an important role in the study of epidemics. The role of one such variable, population density, in the spread of influenza is controversial. Prior studies have tested for such a role using arbitrary thresholds for population density above or below which places are hypothesized to have higher or lower mortality. The results of such studies are mixed. The objective of this study is to estimate, rather than assume, a threshold level of population density that separates low-density regions from high-density regions on the basis of population loss during an influenza pandemic. We study the case of the influenza pandemic of 1918-19 in India, where over 15 million people died in the short span of less than one year.

METHODS

Using data from six censuses for 199 districts of India (n=1194), the country with the largest number of deaths from the influenza of 1918-19, we use a sample-splitting method embedded within a population growth model that explicitly quantifies population loss from the pandemic to estimate a threshold level of population density that separates low-density districts from high-density districts.

RESULTS

The results demonstrate a threshold level of population density of 175 people per square mile. A concurrent finding is that districts on the low side of the threshold experienced rates of population loss (3.72%) that were lower than districts on the high side of the threshold (4.69%).

CONCLUSIONS

This paper introduces a useful analytic tool to the health geographic literature. It illustrates an application of the tool to demonstrate that it can be useful for pandemic awareness and preparedness efforts. Specifically, it estimates a level of population density above which policies to socially distance, redistribute or quarantine populations are likely to be more effective than they are for areas with population densities that lie below the threshold.

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.

Global mapping of infectious disease

The primary aim of this review was to evaluate the state of knowledge of the geographical distribution of all infectious diseases of clinical significance to humans. A systematic review was conducted to enumerate cartographic progress, with respect to the data available for mapping and the methods currently applied. The results helped define the minimum information requirements for mapping infectious disease occurrence, and a quantitative framework for assessing the mapping opportunities for all infectious diseases. This revealed that of 355 infectious diseases identified, 174 (49%) have a strong rationale for mapping and of these only 7 (4%) had been comprehensively mapped. A variety of ambitions, such as the quantification of the global burden of infectious disease, international biosurveillance, assessing the likelihood of infectious disease outbreaks and exploring the propensity for infectious disease evolution and emergence, are limited by these omissions. An overview of the factors hindering progress in disease cartography is provided. It is argued that rapid improvement in the landscape of infectious diseases mapping can be made by embracing non-conventional data sources, automation of geo-positioning and mapping procedures enabled by machine learning and information technology, respectively, in addition to harnessing labour of the volunteer ‘cognitive surplus’ through crowdsourcing.

Mobile phones and malaria: modeling human and parasite travel

Human mobility plays an important role in the dissemination of malaria parasites between regions of variable transmission intensity. Asymptomatic individuals can unknowingly carry parasites to regions where mosquito vectors are available, for example, undermining control programs and contributing to transmission when they travel. Understanding how parasites are imported between regions in this way is therefore an important goal for elimination planning and the control of transmission, and would enable control programs to target the principal sources of malaria. Measuring human mobility has traditionally been difficult to do on a population scale, but the widespread adoption of mobile phones in low-income settings presents a unique opportunity to directly measure human movements that are relevant to the spread of malaria. Here, we discuss the opportunities for measuring human mobility using data from mobile phones, as well as some of the issues associated with combining mobility estimates with malaria infection risk maps to meaningfully estimate routes of parasite importation.

Uncertainties in measuring populations potentially impacted by sea level rise and coastal flooding

A better understanding of the impact of global climate change requires information on the locations and characteristics of populations affected. For instance, with global sea level predicted to rise and coastal flooding set to become more frequent and intense, high-resolution spatial population datasets are increasingly being used to estimate the size of vulnerable coastal populations. Many previous studies have undertaken this by quantifying the size of populations residing in low elevation coastal zones using one of two global spatial population datasets available - LandScan and the Global Rural Urban Mapping Project (GRUMP). This has been undertaken without consideration of the effects of this choice, which are a function of the quality of input datasets and differences in methods used to construct each spatial population dataset. Here we calculate estimated low elevation coastal zone resident population sizes from LandScan and GRUMP using previously adopted approaches, and quantify the absolute and relative differences achieved through switching datasets. Our findings suggest that the choice of one particular dataset over another can translate to a difference of more than 7.5 million vulnerable people for countries with extensive coastal populations, such as Indonesia and Japan. Our findings also show variations in estimates of proportions of national populations at risk range from <0.1% to 45% differences when switching between datasets, with large differences predominantly for countries where coarse and outdated input data were used in the construction of the spatial population datasets. The results highlight the need for the construction of spatial population datasets built on accurate, contemporary and detailed census data for use in climate change impact studies and the importance of acknowledging uncertainties inherent in existing spatial population datasets when estimating the demographic impacts of climate change.

Displaying R spatial statistics on Google dynamic maps with web applications created by Rwui

BACKGROUND

The R project includes a large variety of packages designed for spatial statistics. Google dynamic maps provide web based access to global maps and satellite imagery. We describe a method for displaying directly the spatial output from an R script on to a Google dynamic map.

METHODS

This is achieved by creating a Java based web application which runs the R script and then displays the results on the dynamic map. In order to make this method easy to implement by those unfamiliar with programming Java based web applications, we have added the method to the options available in the R Web User Interface (Rwui) application. Rwui is an established web application for creating web applications for running R scripts. A feature of Rwui is that all the code for the web application being created is generated automatically so that someone with no knowledge of web programming can make a fully functional web application for running an R script in a matter of minutes.

RESULTS

Rwui can now be used to create web applications that will display the results from an R script on a Google dynamic map. Results may be displayed as discrete markers and/or as continuous overlays. In addition, users of the web application may select regions of interest on the dynamic map with mouse clicks and the coordinates of the region of interest will automatically be made available for use by the R script.

CONCLUSIONS

This method of displaying R output on dynamic maps is designed to be of use in a number of areas. Firstly it allows statisticians, working in R and developing methods in spatial statistics, to easily visualise the results of applying their methods to real world data. Secondly, it allows researchers who are using R to study health geographics data, to display their results directly onto dynamic maps. Thirdly, by creating a web application for running an R script, a statistician can enable users entirely unfamiliar with R to run R coded statistical analyses of health geographics data. Fourthly, we envisage an educational role for such applications.

Mapping populations at risk: improving spatial demographic data for infectious disease modeling and metric derivation

The use of Global Positioning Systems (GPS) and Geographical Information Systems (GIS) in disease surveys and reporting is becoming increasingly routine, enabling a better understanding of spatial epidemiology and the improvement of surveillance and control strategies. In turn, the greater availability of spatially referenced epidemiological data is driving the rapid expansion of disease mapping and spatial modeling methods, which are becoming increasingly detailed and sophisticated, with rigorous handling of uncertainties. This expansion has, however, not been matched by advancements in the development of spatial datasets of human population distribution that accompany disease maps or spatial models.Where risks are heterogeneous across population groups or space or dependent on transmission between individuals, spatial data on human population distributions and demographic structures are required to estimate infectious disease risks, burdens, and dynamics. The disease impact in terms of morbidity, mortality, and speed of spread varies substantially with demographic profiles, so that identifying the most exposed or affected populations becomes a key aspect of planning and targeting interventions. Subnational breakdowns of population counts by age and sex are routinely collected during national censuses and maintained in finer detail within microcensus data. Moreover, demographic and health surveys continue to collect representative and contemporary samples from clusters of communities in low-income countries where census data may be less detailed and not collected regularly. Together, these freely available datasets form a rich resource for quantifying and understanding the spatial variations in the sizes and distributions of those most at risk of disease in low income regions, yet at present, they remain unconnected data scattered across national statistical offices and websites.In this paper we discuss the deficiencies of existing spatial population datasets and their limitations on epidemiological analyses. We review sources of detailed, contemporary, freely available and relevant spatial demographic data focusing on low income regions where such data are often sparse and highlight the value of incorporating these through a set of examples of their application in disease studies. Moreover, the importance of acknowledging, measuring, and accounting for uncertainty in spatial demographic datasets is outlined. Finally, a strategy for building an open-access database of spatial demographic data that is tailored to epidemiological applications is put forward.