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

Population Spatialization in Beijing City Based on Machine Learning and Multisource Remote Sensing Data

Remote sensing data have been widely used in research on population spatialization. Previous studies have generally divided study areas into several sub-areas with similar features by artificial or clustering algorithms and then developed models for these sub-areas separately using statistical methods. These approaches have drawbacks due to their subjectivity and uncertainty. In this paper, we present a study of population spatialization in Beijing City, China based on multisource remote sensing data and town-level population census data. Six predictive algorithms were compared for estimating population using the spatial variables derived from The National Polar-Orbiting Partnership/ Visible Infrared Imaging Radiometer Suite (NPP/VIIRS) night-time light and other remote sensing data. Random forest achieved the highest accuracy and therefore was employed for population spatialization. Feature selection was performed to determine the optimal variable combinations for population modeling by random forest. Cross-validation results indicated that the developed model achieved a mean absolute error (MAE) of 2129.52 people/km2 and a R2 of 0.63. The gridded population density in Beijing at a spatial resolution of 500 m produced by the random forest model was also adjusted to be consistent with the census population at the town scale. By comparison with Google Earth high-resolution images, the remotely-sensed population was qualitatively validated at the intra-town scale. Validation results indicated that remotely sensed results can effectively depict the spatial distribution of population within town-level districts. This study provides a valuable reference for urban planning, public health and disaster prevention in Beijing, and a reference for population mapping in other cities.

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.

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

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

GHS-POP Accuracy Assessment: Poland and Portugal Case Study

The Global Human Settlement Population Grid (GHS-POP) the latest released global gridded population dataset based on remotely sensed data and developed by the EU Joint Research Centre, depicts the distribution and density of the total population as the number of people per grid cell. This study aims to assess the GHS-POP data accuracy based on root mean square error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE) and the correlation coefficient. The study was conducted for Poland and Portugal, countries characterized by different population distribution as well as two spatial resolutions of 250 m and 1 km on the GHS-POP. The main findings show that as the size of administrative zones decreases (from NUTS (Nomenclature of Territorial Units for Statistics) to LAU (local administrative unit)) and the size of the GHS-POP increases, the difference between the population counts reported by the European Statistical Office and estimated by the GHS-POP algorithm becomes larger. At the national level, MAPE ranges from 1.8% to 4.5% for the 250 m and 1 km resolutions of GHS-POP data in Portugal and 1.5% to 1.6%, respectively in Poland. At the local level, however, the error rates range from 4.5% to 5.8% in Poland, for 250 m and 1 km, and 5.7% to 11.6% in Portugal, respectively. Moreover, the results show that for densely populated regions the GHS-POP underestimates the population number, while for thinly populated regions it overestimates. The conclusions of this study are expected to serve as a quality reference for potential users and producers of population density datasets.

Using gridded population and quadtree sampling units to support survey sample design in low-income settings

BACKGROUND

Household surveys are the main source of demographic, health and socio-economic data in low- and middle-income countries (LMICs). To conduct such a survey, census population information mapped into enumeration areas (EAs) typically serves a sampling frame from which to generate a random sample. However, the use of census information to generate this sample frame can be problematic as in many LMIC contexts, such data are often outdated or incomplete, potentially introducing coverage issues into the sample frame. Increasingly, where census data are outdated or unavailable, modelled population datasets in the gridded form are being used to create household survey sampling frames.

METHODS

Previously this process was done by either sampling from a set of the uniform grid cells (UGC) which are then manually subdivided to achieve the desired population size, or by sampling very small grid cells then aggregating cells into larger units to achieve a minimum population per survey cluster. The former approach is time and resource-intensive as well as results in substantial heterogeneity in the output sampling units, while the latter can complicate the calculation of unbiased sampling weights. Using the context of Somalia, which has not had a full census since 1987, we implemented a quadtree algorithm for the first time to create a population sampling frame. The approach uses gridded population estimates and it is based on the idea of a quadtree decomposition in which an area successively subdivided into four equal size quadrants, until the content of each quadrant is homogenous.

RESULTS

The quadtree approach used here produced much more homogeneous sampling units than the UGC (1 × 1 km and 3 × 3 km) approach. At the national and pre-war regional scale, the standard deviation and coefficient of variation, as indications of homogeneity, were calculated for the output sampling units using quadtree and UGC 1 × 1 km and 3 × 3 km approaches to create the sampling frame and the results showed outstanding performance for quadtree approach.

CONCLUSION

Our approach reduces the manual burden of manually subdividing UGC into highly populated areas, while allowing for correct calculation of sampling weights. The algorithm produces a relatively homogenous population counts within the sampling units, reducing the variation in the weights and improving the precision of the resulting estimates. Furthermore, a protocol of creating approximately equal-sized blocks and using tablets for randomized selection of a household in each block mitigated potential selection bias by enumerators. The approach shows labour, time and cost-saving and points to the potential use in wider contexts.

Scale and context dependency of deforestation drivers: Insights from spatial econometrics in the tropics

A better understanding of deforestation drivers across countries and spatial scales is a precondition for designing efficient international policies and coherent land use planning strategies such as REDD+. However, it is so far unclear if the well-studied drivers of tropical deforestation behave similarly across nested subnational jurisdictions, which is crucial for efficient policy implementation. We selected three countries in Africa, America and Asia, which present very different tropical contexts. Making use of spatial econometrics and a multi-level approach, we conducted a set of regressions comprising 3,035 administrative units from the three countries at micro-level, plus 361 and 49 at meso- and macro-level, respectively. We included forest cover as dependent variable and seven physio-geographic and socioeconomic indicators of well-known drivers of deforestation as explanatory variables. With this, we could provide a first set of highly significant econometric models of pantropical deforestation that consider subnational units. We identified recurrent drivers across countries and scales, namely population pressure and the natural condition of land suitability for crop production. The impacts of demography on forest cover were strikingly strong across contexts, suggesting clear limitations of sectoral policy. Our findings also revealed scale and context dependencies, such as an increased heterogeneity at local scopes, with a higher and more diverse number of significant determinants of forest cover. Additionally, we detected stronger spatial interactions at smaller levels, providing empirical evidence that certain deforestation forces occur independently of the existing de jure governance boundaries. We demonstrated that neglecting spatial dependencies in this type of studies can lead to several misinterpretations. We therefore advocate, that the design and enforcement of policy instruments-such as REDD+-should start from common international entry points that ensure for coherent agricultural and demographic policies. In order to achieve a long-term impact on the ground, these policies need to have enough flexibility to be modified and adapted to specific national, regional or local conditions.

Mitigation pathways towards national ambient air quality standards in India

Exposure to ambient particulate matter is a leading risk factor for environmental public health in India. While Indian authorities implemented several measures to reduce emissions from the power, industry and transportation sectors over the last years, such strategies appear to be insufficient to reduce the ambient fine particulate matter (PM_2.5) concentration below the Indian National Ambient Air Quality Standard (NAAQS) of 40 μg/m³ across the country. This study explores pathways towards achieving the NAAQS in India in the context of the dynamics of social and economic development. In addition, to inform action at the subnational levels in India, we estimate the exposure to ambient air pollution in the current legislations and alternative policy scenarios based on simulations with the GAINS integrated assessment model. The analysis reveals that in many of the Indian States emission sources that are outside of their immediate jurisdictions make the dominating contributions to (population-weighted) ambient pollution levels of PM_2.5. Consequently, most of the States cannot achieve significant improvements in their air quality and population exposure on their own without emission reductions in the surrounding regions, and any cost-effective strategy requires regionally coordinated approaches. Advanced technical emission control measures could provide NAAQS-compliant air quality for 60% of the Indian population. However, if combined with national sustainable development strategies, an additional 25% population will be provided with clean air, which appears to be a significant co-benefit on air quality (totaling 85%).

Modeling Population Density using a New Index Derived from Multi-Sensor Image Data

The detailed information about the spatial distribution of the population is crucial for analyzing economic growth, environmental change, and natural disaster damage. Using the nighttime light (NTL) imagery for population estimation has been a topic of interest in recent decades. However, the effectiveness of NTL data in population estimation has been impeded by some limitations such as the blooming effect and underestimation in rural regions. To overcome these limitations, we combine the NPP-VIIRS day/night band (DNB) data with normalized difference vegetation index (NDVI) and land surface temperature (LST) data derived from the moderate resolution imaging spectroradiometer (MODIS) onboard the Terra satellite, to create a new vegetation temperature light population index (VTLPI). A statistical model is developed to predict 250m grid-level population density based on the proposed VTLPI and the least square regression approach. After that, a case study is implemented using the data of Sichuan Province, China in 2015, and the results indicates that the VTLPI-estimated population density outperformed the results from other two methods based on nighttime light imagery or human settlement index, and the three publicized population products, LandScan, WorldPop, and GPW. When using the census data as reference, the mean relative error and median absolute relative error on a township level are 0.29 and 0.12, respectively, and the root-mean-square error is 212 persons/km2. The results show that our VTLPI-based model can achieve a better estimation of population density in rural areas and urban suburbs and characterize more spatial variations at 250m grid level both in both urban and rural areas. The resultant population density offers better population exposure data for assessing natural disaster risk and loss as well as other related applications.

Coastal proximity of populations in 22 Pacific Island Countries and Territories

The coastal zones of Small Island States are hotspots of human habitation and economic endeavour. In the Pacific region, as elsewhere, there are large gaps in understandings of the exposure and vulnerability of people in coastal zones. The 22 Pacific Countries and Territories (PICTs) are poorly represented in global analyses of vulnerability to seaward risks. We combine several data sources to estimate populations to zones 1, 5 and 10 km from the coastline in each of the PICTs. Regional patterns in the proximity of Pacific people to the coast are dominated by Papua New Guinea. Overall, ca. half the population of the Pacific resides within 10 km of the coast but this jumps to 97% when Papua New Guinea is excluded. A quarter of Pacific people live within 1 km of the coast, but without PNG this increases to slightly more than half. Excluding PNG, 90% of Pacific Islanders live within 5 km of the coast. All of the population in the coral atoll nations of Tokelau and Tuvalu live within a km of the ocean. Results using two global datasets, the SEDAC-CIESIN Gridded Population of the World v4 (GPWv4) and the Oak Ridge National Laboratory Landscan differed: Landscan under-dispersed population, overestimating numbers in urban centres and underestimating population in rural areas and GPWv4 over-dispersed the population. In addition to errors introduced by the allocation models of the two methods, errors were introduced as artefacts of allocating households to 1 km x 1 km grid cell data (30 arc-seconds) to polygons. The limited utility of LandScan and GPWv4 in advancing this analysis may be overcome with more spatially resolved census data and the inclusion of elevation above sea level as an important dimension of vulnerability.

Can More Accurate Night-Time Remote Sensing Data Simulate a More Detailed Population Distribution?

Aging, shrinking cities, urban agglomerations and other new key terms continue to emerge when describing the large-scale population changes in various cities in mainland China. It is important to simulate the distribution of residential populations at a coarse scale to manage cities as a whole, and at a fine scale for policy making in infrastructure development. This paper analyzes the relationship between the DN (Digital number, value assigned to a pixel in a digital image) value of NPP-VIIRS (the Suomi National Polar-orbiting Partnership satellite’s Visible Infrared Imaging Radiometer Suite) and LuoJia1-01 and the residential populations of urban areas at a district, sub-district, community and court level, to compare the influence of resolution of remote sensing data by taking urban land use to map out auxiliary data in which first-class (R1), second-class (R2) and third-class residential areas (R3) are distinguished by house price. The results show that LuoJia1-01 more accurately analyzes population distributions at a court level for second- and third-class residential areas, which account for over 85% of the total population. The accuracy of the LuoJia1-01 simulation data is higher than that of Landscan and GHS (European Commission Global Human Settlement) population. This can be used as an important tool for refining the simulation of residential population distributions. In the future, higher-resolution night-time light data could be used for research on accurate simulation analysis that scales down large-scale populations.

Contrasting the value of targeted versus area-wide mosquito control scenarios to limit arbovirus transmission with human mobility patterns based on different tropical urban population centers

Vector control is still our primary intervention for both prevention and mitigation of epidemics of many vector-borne diseases. Efficiently targeting control measures is important since control can involve substantial economic costs. Targeting is not always straightforward, as transmission of vector-borne diseases is affected by various types of host movement. Here we assess how taking daily commuting patterns into consideration can help improve vector control efforts. We examine three tropical urban centers (San Juan, Recife, and Jakarta) that have recently been exposed to Zika and/or dengue infections and consider whether the distribution of human populations and resulting commuting flows affects the optimal scale at which control interventions should be implemented. We developed a stochastic, spatial model and investigated four control scenarios. The scenarios differed in the spatial extent of their implementation and were: 1) a response at the level of an individual neighborhood; 2) a response targeted at a neighborhood in which infected humans were detected and the one with which it was most strongly connected by human movement; 3) a limited area-wide response where all neighborhoods within a certain radius of the focal area were included; and 4) a collective response where all participating neighborhoods implemented control. The relative effectiveness of the scenarios varied only slightly between different settings, with the number of infections averted over time increasing with the scale of implementation. This difference depended on the efficacy of control at the neighborhood level. At low levels of efficacy, the scenarios mirrored each other in infections averted. At high levels of efficacy, impact increased with the scale of the intervention. As a result, the choice between scenarios will not only be a function of the amount of effort decision-makers are willing to invest, but largely epend on the overall effectiveness of vector control approaches.

Control and prevention of Aedes-transmitted viruses, such as dengue, chikungunya, or Zika relies heavily on vector control approaches. Given the effort and cost involved in implementation of vector control, targeting of control measures is highly desirable. However, it is unclear to what extent the effectiveness of highly focal and reactive control measures depends on the commuting and movement patterns of humans. To investigate this question, we developed a model and four control scenarios that ranged from highly focal to area-wide larval control. The distribution of humans and their commuting patterns were modelled after three major tropical urban centers, San Juan, Recife, and Jakarta. We show that as implementation is applied across a wider area, a greater number of infections is averted. Critically, this only occurs if the efficacy of control at the neighborhood level is sufficiently high. A consistent outcome across the three settings was that the focal strategy was most likely to provide the best outcome at lower levels of effort, and when the efficacy of control was low. These outcomes suggest that optimal control strategies will likely have to be tailored to individual settings by decision makers and would benefit from localized cost-effectiveness modelling studies.

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

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

Decreasing measles burden by optimizing campaign timing

Measles vaccine is a highly effective healthcare intervention, but getting vaccine to those in need remains a major problem. Complicating the issue, high-burden countries typically have low-quality infrastructure, severely limiting the number of infections detected and therefore limiting our understanding of local epidemiology. Here we show that statistical disease models can be fitted to sparse case data from Pakistan using a fast linear regression approach. This method yields estimates of the effects of past interventions, the seasonal likelihood of measles transmission, and the magnitude of future outbreaks under different intervention policies. We use these models to understand in general when and where vaccine should be distributed, and these results were used to inform Pakistan’s 2018 vaccination campaign planning.

Measles remains a major contributor to preventable child mortality, and bridging gaps in measles immunity is a fundamental challenge to global health. In high-burden settings, mass vaccination campaigns are conducted to increase access to vaccine and address this issue. Ensuring that campaigns are optimally effective is a crucial step toward measles elimination; however, the relationship between campaign impact and disease dynamics is poorly understood. Here, we study measles in Pakistan, and we demonstrate that campaign timing can be tuned to optimally interact with local transmission seasonality and recent incidence history. We develop a mechanistic modeling approach to optimize timing in general high-burden settings, and we find that in Pakistan, hundreds of thousands of infections can be averted with no change in campaign cost.

A spatial assessment of Nipah virus transmission in Thailand pig farms using multi-criteria decision analysis

BACKGROUND

Thailand's Central Plain is identified as a contact zone between pigs and flying foxes, representing a potential zoonotic risk. Nipah virus (NiV) has been reported in flying foxes in Thailand, but it has never been found in pigs or humans. An assessment of the suitability of NiV transmission at the spatial and farm level would be useful for disease surveillance and prevention. Multi-criteria decision analysis (MCDA), a knowledge-driven model, was used to map contact zones between local epizootic risk factors as well as to quantify the suitability of NiV transmission at the pixel and farm level.

RESULTS

Spatial risk factors of NiV transmission in pigs were identified by experts as being of three types, including i) natural host factors (bat preferred areas and distance to the nearest bat colony), ii) intermediate host factors (pig population density), and iii) environmental factors (distance to the nearest forest, distance to the nearest orchard, distance to the nearest water body, and human population density). The resulting high suitable areas were concentrated around the bat colonies in three provinces in the East of Thailand, including Chacheongsao, Chonburi, and Nakhonnayok. The suitability of NiV transmission in pig farms in the study area was quantified as ranging from very low to medium suitability.

CONCLUSIONS

We believe that risk-based surveillance in the identified priority areas may increase the chances of finding out NiV and other bat-borne pathogens and thereby optimize the allocation of financial resources for disease surveillance. In the long run, improvements of biosecurity in those priority areas may also contribute to preventing the spread of potential emergence of NiV and other bat-borne pathogens.

The Impact of Portfolio Location Uncertainty on Probabilistic Seismic Risk Analysis

Probabilistic seismic risk analysis is a well-established method in the insurance industry for modeling portfolio losses from earthquake events. In this context, precise exposure locations are often unknown. However, so far, location uncertainty has not been in the focus of a large amount of research. In this article, we propose a novel framework for treatment of location uncertainty. As a case study, a large number of synthetic portfolios resembling typical real-world cases were created. We investigate the effect of portfolio characteristics such as value distribution, portfolio size, or proportion of risk items with unknown coordinates on the variability of loss frequency estimations. The results indicate that due to loss aggregation effects and spatial hazard variability, location uncertainty in isolation and in conjunction with ground motion uncertainty can induce significant variability to probabilistic loss results, especially for portfolios with a small number of risks. After quantifying its effect, we conclude that location uncertainty should not be neglected when assessing probabilistic seismic risk, but should be treated stochastically and the resulting variability should be visualized and interpreted carefully.

Spatial and temporal projections of the prevalence of active tuberculosis in Cambodia

Introduction

Cambodia is among the 30 highest burden of tuberculosis (TB) countries. Active TB prevalence has been estimated using nationally representative multistage sampling that represents urban, rural and remote parts of the country, but the prevalence in non-sampled communes remains unknown. This study uses geospatial Bayesian statistics to estimate point prevalence across Cambodia, and demographic modelling that accounts for secular trends in fertility, mortality, urbanisation and prevalence rates to project the future burden of active TB.

Methods

A Bayesian hierarchical model was developed for the 2011 National Tuberculosis Prevalence survey to estimate the differential effect of age, sex and geographic stratum on active TB prevalence; these estimates were then married with high-resolution geographic information system layers to project prevalence across Cambodia. Future TB projections under alternative scenarios were then derived by interfacing these estimates with an individual-based demographic model.

Results

Strong differences in risk by age and sex, together with geographically varying population structures, yielded the first estimated prevalence map at a 1 km scale. The projected number of active TB cases within the catchment area of each existing government healthcare facility was derived, together with projections to the year 2030 under three scenarios: no future improvement, continualreduction and GDPprojection.

Conclusion

Synthesis of health and geographic data allows likely disease rates to be mapped at a high resolution to facilitate resource planning, while demographic modelling allows scenarios to be projected, demonstrating the need for the acceleration of control efforts to achieve a substantive impact on the future burden of TB in Cambodia.

Dasymetric mapping of urban population in China based on radiance corrected DMSP-OLS nighttime light and land cover data

High spatial resolution urban population dataset is increasingly required for sustainable urban planning and management. Dasymetric mapping is an effective approach to create such dataset. However, the created gridded total population datasets usually have limitation for urban analysis in developing countries as they usually underestimate urban population because of the strong urban-rural difference. In this study, we aimed to create a dataset of gridded urban population with 1 km resolution in China in year 2000 and 2010. We proposed an index of urban nighttime light (UNTL) by integrating radiance corrected DMSP nighttime light (RcNTL) and urban land, which is then used as weight to disaggregate county-level urban population. The validation using township population in Beijing as references shows reasonable accuracy with a mean relative error of 38% and a R² of 68%. Using only two widely available datasets (RcNTL and urban land), the proposed method is simple and computing efficient compared with methods using multiple geospatial data (e.g., land use and land cover, distance to city center, slope) and that combined with remote sensing imagery. As the used two auxiliary datasets are accessible globally, the method has great potential to produce similar urban population dataset for other developing countries where fine scale census population datasets are scarce. The produced urban population dataset is valuable for enriching our understanding of the urbanization process and designing sustainable urban planning and management strategies in China.

Estimating the eco-toxicological risk of estrogens in China's rivers using a high-resolution contaminant fate model

The contamination of freshwater systems arises in many river basins due to industrialization and population growth, posing risks to ecosystems and human health. Despite these concerns, the fate and potential impact of many emerging pollutants are currently unknown, especially when the compounds are released into surface waters from populations distributed across large spatial scales. In order to address this shortcoming, a spatially-explicit contaminant fate model was developed as an extension of the global, vector-based river routing model HydroROUT. HydroROUT operates at very high spatial resolution (∼500 m), simulates river and stream chemical transport with in-stream removal, and contains links to a set of lakes and reservoirs, which act as a partial sink during the transport. The chemical fate model was applied to China and includes a consumption and release module based on county-level population demographics, considers point-source contributions from 2547 wastewater treatment plants, and accumulates contributions of rural and urban populations not connected to sewage treatment plants. As a case study, the sources and fates of the estrogens estrone (E₁), 17β-estradiol (E₂), estriol (E₃), as well as the synthetic estrogenic steroid hormone 17α-ethinylestradiol (EE₂) were modelled in Chinese surface water bodies. Preliminary validation of the results showed predictions to be within the ranges of concentrations reported in literature, with exception of EE₂. The total estrogenic mass in the entire river and lake system amounted to 370 tonnes of estrogens, with about 1.3 tonnes per year discharged to the ocean, neighboring countries or to inland sinks. Under a selected baseline scenario, eco-toxicological risk-i.e., contaminant concentrations in excess of the predicted no effect concentration (PNEC)-is found in 23.6% of all analyzed rivers of China with an average flow > 0.1 m³/s. Out of these, about 4% of rivers showed a high level of risk of 10 times or more above PNEC. Medium-to-large rivers are disproportionally affected; for example, 23.6%, 37.3%, 29.0% and 21.6% of river length are at risk in rivers of 1-10, 10-100, 100-1,000, and 1,000-10,000 m³/s of discharge, respectively, whereas no risk was predicted in the largest rivers (i.e., >10,000 m³/s) of China. Wastewater treatment plants process 22.5% of the total hormone load and thus play an important role in water quality control by reducing the risk in substantial portions of the river network, which would otherwise show elevated risk. Releases from untreated population dominate by far the overall contribution to risk.

Spatial distribution and populations at risk of A. lumbricoides and T. trichiura co-infections and infection intensity classes: an ecological study

BACKGROUND

Soil-transmitted helminth (STH) infections are highly prevalent in the Philippines. Mapping the prevalence and high-intensity of STH co-infections can help guide targeted intervention programmes to reduce morbidity, especially among vulnerable school-aged children. In this study, we aimed to predict the spatial distribution of the prevalence of Ascaris lumbricoides and Trichuris trichiura co-infection and infection intensity classes in the Philippines to identify populations most in need of interventions.

METHODS

Data on STH infections from 29,919 individuals during the nationwide parasitological survey in 2005 to 2007 were included in the analysis. To geographically predict the prevalence of A. lumbricoides and T. trichiura co-infections and infection intensity classes, Bayesian multinomial geostatistical models were built including age, sex, environmental variables and a geostatistical random effect. The number of individuals co-infected and belonging to each of the infection intensity classes in 2017 was forecast by combining our predictive prevalence maps with population density maps.

RESULTS

Our models showed that school-aged children (5-19 years) are most at risk of A. lumbricoides and T. trichiura co-infections and of moderate/high infection intensity compared to other age groups. We identified target provinces where the likelihood of STH-associated morbidity was highest: Luzon (Bulacan, Benguet, Cavite, Sorsogon, Metropolitan Manila, Pampanga and Rizal), the Visayas (Cebu, Iloilo, Leyte and Negros Occidental), and in Mindanao (Agusan Del Norte, Davao Del Sur, Davao Oriental, Lanao Del Sur, Maguindanao, Misamis Oriental, Sulu and Zamboanga Del Sur). Luzon had the highest estimated number of school-aged children with A. lumbricoides and T. trichiura co-infections (estimated total 89,400), followed by the Visayas (38,300) and Mindanao (20,200).

CONCLUSIONS

Our study provided epidemiological evidence to highlight national priority areas for controlling co-infections and high intensity infections in the Philippines. Our maps could assist more geographically targeted interventions to reduce the risk of STH-associated morbidity in the Philippines.

A spatial regression model for the disaggregation of areal unit based data to high-resolution grids with application to vaccination coverage mapping

The growing demand for spatially detailed data to advance the Sustainable Development Goals agenda of ‘leaving no one behind’ has resulted in a shift in focus from aggregate national and province-based metrics to small areas and high-resolution grids in the health and development arena. Vaccination coverage is customarily measured through aggregate-level statistics, which mask fine-scale heterogeneities and ‘coldspots’ of low coverage. This paper develops a methodology for high-resolution mapping of vaccination coverage using areal data in settings where point-referenced survey data are inaccessible. The proposed methodology is a binomial spatial regression model with a logit link and a combination of covariate data and random effects modelling two levels of spatial autocorrelation in the linear predictor. The principal aspect of the model is the melding of the misaligned areal data and the prediction grid points using the regression component and each of the conditional autoregressive and the Gaussian spatial process random effects. The Bayesian model is fitted using the INLA-SPDE approach. We demonstrate the predictive ability of the model using simulated data sets. The results obtained indicate a good predictive performance by the model, with correlations of between 0.66 and 0.98 obtained at the grid level between true and predicted values. The methodology is applied to predicting the coverage of measles and diphtheria-tetanus-pertussis vaccinations at 5 × 5 km² in Afghanistan and Pakistan using subnational Demographic and Health Surveys data. The predicted maps are used to highlight vaccination coldspots and assess progress towards coverage targets to facilitate the implementation of more geographically precise interventions. The proposed methodology can be readily applied to wider disaggregation problems in related contexts, including mapping other health and development indicators.

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

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

Rapid Population Growth throughout Asia's Earthquake-Prone Areas: A Multiscale Analysis

Assessing the changes of the population living throughout the most seismically hazardous area (MSHA) constitutes an important foundation for seismic risk assessment. However, the changes of the population living in the MSHA of Asia, which exhibits the highest number of earthquake related fatalities, were poorly understood. Therefore, this study analyzed the changes of the population in the MSHA between 2000 and 2015 at the continental, subcontinental, and national scales. We found that the population, especially the vulnerable population (i.e., children under or equal to the age of 14 and elderly people over or equal to the age of 65), in Asia’s MSHA increased rapidly between 2000 and 2015. The population in the MSHA increased by 185.88 million with a growth rate of 20.93%, which was 3.38% greater than that in the non-MSHA region. Meanwhile, the vulnerable population in the MSHA increased by 63.65 million with a growth rate of 19.73%. The increase of the vulnerable population in the MSHA was 19.93% greater than that in the non-MSHA region. We also found that urban population growth was a major factor impacting the increase in both the population and the vulnerable population throughout Asia’s MSHA. Therefore, attention should be paid to the changes of the population in Asia’s MSHA, whilst it is imperative to execute strict building codes and select the development location more carefully in the MSHA.

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

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

A spatio-temporal land use and land cover reconstruction for India from 1960-2010

In recent decades India has undergone substantial land use/land cover change as a result of population growth and economic development. Historical land use/land cover maps are necessary to quantify the impact of change at global and regional scales, improve predictions about the quantity and location of future change and support planning decisions. Here, a regional land use change model driven by district-level inventory data is used to generate an annual time series of high-resolution gridded land use/land cover maps for the Indian subcontinent between 1960-2010. The allocation procedure is based on statistical analysis of the relationship between contemporary land use/land cover and various spatially explicit covariates. A comparison of the simulated map for 1985 against remotely-sensed land use/land cover maps for 1985 and 2005 reveals considerable discrepancy between the simulated and remote sensing maps, much of which arises due to differences in the amount of land use/land cover change between the inventory data and the remote sensing maps.

Individually-tailored multifactorial intervention to reduce falls in the Malaysian Falls Assessment and Intervention Trial (MyFAIT): A randomized controlled trial

Objective

To determine the effectiveness of an individually-tailored multifactorial intervention in reducing falls among at risk older adult fallers in a multi-ethnic, middle-income nation in South-East Asia.

Design

Pragmatic, randomized-controlled trial.

Setting

Emergency room, medical outpatient and primary care clinic in a teaching hospital in Kuala Lumpur, Malaysia.

Participants

Individuals aged 65 years and above with two or more falls or one injurious fall in the past 12 months.

Intervention

Individually-tailored interventions, included a modified Otago exercise programme, HOMEFAST home hazards modification, visual intervention, cardiovascular intervention, medication review and falls education, was compared against a control group involving conventional treatment.

Primary and secondary outcome measures

The primary outcome was any fall recurrence at 12-month follow-up. Secondary outcomes were rate of fall and time to first fall.

Results

Two hundred and sixty-eight participants (mean age 75.3 ±7.2 SD years, 67% women) were randomized to multifactorial intervention (n = 134) or convention treatment (n = 134). All participants in the intervention group received medication review and falls education, 92 (68%) were prescribed Otago exercises, 86 (64%) visual intervention, 64 (47%) home hazards modification and 51 (38%) cardiovascular intervention. Fall recurrence did not differ between intervention and control groups at 12-months [Risk Ratio, RR = 1.037 (95% CI 0.613–1.753)]. Rate of fall [RR = 1.155 (95% CI 0.846–1.576], time to first fall [Hazard Ratio, HR = 0.948 (95% CI 0.782–1.522)] and mortality rate [RR = 0.896 (95% CI 0.335–2.400)] did not differ between groups.

Conclusion

Individually-tailored multifactorial intervention was ineffective as a strategy to reduce falls. Future research efforts are now required to develop culturally-appropriate and affordable methods of addressing this increasingly prominent public health issue in middle-income nations.

Trial registration

ISRCTN Registry no. ISRCTN11674947

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

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

Population sensitivity of acute flaccid paralysis and environmental surveillance for serotype 1 poliovirus in Pakistan: an observational study

Background

To support poliomyelitis eradication in Pakistan, environmental surveillance (ES) of wastewater has been expanded alongside surveillance for acute flaccid paralysis (AFP). ES is a relatively new method of surveillance, and the population sensitivity of detecting poliovirus within endemic settings requires estimation.

Methods

Data for wild serotype 1 poliovirus from AFP and ES from January 2011 to September 2015 from 14 districts in Pakistan were analysed using a multi-state model framework. This framework was used to estimate the sensitivity of poliovirus detection from each surveillance source and parameters such as the duration of infection within a community.

Results

The location and timing of poliomyelitis cases showed spatial and temporal variability. The sensitivity of AFP surveillance to detect serotype 1 poliovirus infection in a district and its neighbours per month was on average 30.0% (95% CI 24.8–35.8) and increased with the incidence of poliomyelitis cases. The average population sensitivity of a single environmental sample was 59.4% (95% CI 55.4–63.0), with significant variation in site-specific estimates (median varied from 33.3–79.2%). The combined population sensitivity of environmental and AFP surveillance in a given month was on average 98.1% (95% CI 97.2–98.7), assuming four samples per month for each site.

Conclusions

ES can be a highly sensitive supplement to AFP surveillance in areas with converging sewage systems. As ES for poliovirus is expanded, it will be important to identify factors associated with variation in site sensitivity, leading to improved site selection and surveillance system performance.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3070-4) contains supplementary material, which is available to authorized users.

Rat-bites of an epidemic proportion in Peshawar vale; a GIS based approach in risk assessment

Contemporary studies demonstrate that rodent bites do not occur frequently. However, a huge number of cases were reported from Peshawar vale, Pakistan during 2016. Two species, the local black rat Rattus rattus (Linnaeus, 1758) and the invasive brown rat Rattus norvegicus (Berkenhout, 1769) might be the suspected cause. Several studies indicated the invasion of brown rats into Pakistan presumably via port city of Karachi. In this study, we modeled geospatial distribution of rodent bites for risk assessment in the region. Bite cases reported to tertiary care lady reading hospital were monitored from January 1 to August 31, 2016. Among 1747 cases, statistically informative data (n = 1295) was used for analyses. MaxEnt algorithm was employed for geospatial modeling, taking into account various environmental variables (temperature, precipitation, humidity, and elevation) and anthropogenic factors (human population density, distance from roads, distance from water channels, and land use/land cover). MaxEnt results revealed that urban slums (84.5%) are at highest risk followed by croplands (10.9%) and shrublands (2.7%). Anthropogenic factors affecting incidence of rodent bites included host density (contribution: 34.7), distance from water channels (3.2), land use/land cover (2.8), and distance from roads (2). Most of the cases occurred within a radius of 0.3 km from roads and 5 km from water channels. Rodent bite incidence is currently at its peak in Peshawar vale. Factors significantly affecting rodents' bite activity and their distribution and dispersal include urbanization, distance from roads, and water channels. Further studies are needed to determine the impact of invasion by brown rat on bite incidence.

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.

Neglected Tropical Disease Control - The Case for Adaptive, Location-specific Solutions

The world is experiencing environmental and social change at an unprecedented rate, with the effects being felt at local, regional, and international scales. This phenomenon may disrupt interventions against neglected tropical diseases (NTDs) that operate on the basis of linear scaling and 'one-size-fits-all'. Here we argue that investment in field-based data collection and building modelling capacity is required; that it is important to consider unintended consequences of interventions; that inferences can be drawn from wildlife ecology; and that interventions should become more location-specific. Collectively, these ideas underpin the development of adaptive decision-support tools that are sufficiently flexible to address emerging issues within the Anthropocene.

A global map of travel time to cities to assess inequalities in accessibility in 2015

The economic and man-made resources that sustain human wellbeing are not distributed evenly across the world, but are instead heavily concentrated in cities. Poor access to opportunities and services offered by urban centres (a function of distance, transport infrastructure, and the spatial distribution of cities) is a major barrier to improved livelihoods and overall development. Advancing accessibility worldwide underpins the equity agenda of 'leaving no one behind' established by the Sustainable Development Goals of the United Nations. This has renewed international efforts to accurately measure accessibility and generate a metric that can inform the design and implementation of development policies. The only previous attempt to reliably map accessibility worldwide, which was published nearly a decade ago, predated the baseline for the Sustainable Development Goals and excluded the recent expansion in infrastructure networks, particularly in lower-resource settings. In parallel, new data sources provided by Open Street Map and Google now capture transportation networks with unprecedented detail and precision. Here we develop and validate a map that quantifies travel time to cities for 2015 at a spatial resolution of approximately one by one kilometre by integrating ten global-scale surfaces that characterize factors affecting human movement rates and 13,840 high-density urban centres within an established geospatial-modelling framework. Our results highlight disparities in accessibility relative to wealth as 50.9% of individuals living in low-income settings (concentrated in sub-Saharan Africa) reside within an hour of a city compared to 90.7% of individuals in high-income settings. By further triangulating this map against socioeconomic datasets, we demonstrate how access to urban centres stratifies the economic, educational, and health status of humanity.

Development of an Asset Value Map for Disaster Risk Assessment in China by Spatial Disaggregation Using Ancillary Remote Sensing Data

The extent of economic losses due to a natural hazard and disaster depends largely on the spatial distribution of asset values in relation to the hazard intensity distribution within the affected area. Given that statistical data on asset value are collected by administrative units in China, generating spatially explicit asset exposure maps remains a key challenge for rapid postdisaster economic loss assessment. The goal of this study is to introduce a top-down (or downscaling) approach to disaggregate administrative-unit level asset value to grid-cell level. To do so, finding the highly correlated "surrogate" indicators is the key. A combination of three data sets-nighttime light grid, LandScan population grid, and road density grid, is used as ancillary asset density distribution information for spatializing the asset value. As a result, a high spatial resolution asset value map of China for 2015 is generated. The spatial data set contains aggregated economic value at risk at 30 arc-second spatial resolution. Accuracy of the spatial disaggregation reflects redistribution errors introduced by the disaggregation process as well as errors from the original ancillary data sets. The overall accuracy of the results proves to be promising. The example of using the developed disaggregated asset value map in exposure assessment of watersheds demonstrates that the data set offers immense analytical flexibility for overlay analysis according to the hazard extent. This product will help current efforts to analyze spatial characteristics of exposure and to uncover the contributions of both physical and social drivers of natural hazard and disaster across space and time.

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.

Assessing ecosystem service provision under climate change to support conservation and development planning in Myanmar

Inclusion of ecosystem services (ES) information into national-scale development and climate adaptation planning has yet to become common practice, despite demand from decision makers. Identifying where ES originate and to whom the benefits flow-under current and future climate conditions-is especially critical in rapidly developing countries, where the risk of ES loss is high. Here, using Myanmar as a case study, we assess where and how ecosystems provide key benefits to the country's people and infrastructure. We model the supply of and demand for sediment retention, dry-season baseflows, flood risk reduction and coastal storm protection from multiple beneficiaries. We find that locations currently providing the greatest amount of services are likely to remain important under the range of climate conditions considered, demonstrating their importance in planning for climate resilience. Overlap between priority areas for ES provision and biodiversity conservation is higher than expected by chance overall, but the areas important for multiple ES are underrepresented in currently designated protected areas and Key Biodiversity Areas. Our results are contributing to development planning in Myanmar, and our approach could be extended to other contexts where there is demand for national-scale natural capital information to shape development plans and policies.

Environmental risk modelling and potential sand fly vectors of cutaneous leishmaniasis in Chitral district: a leishmanial focal point of mount Tirich Mir, Pakistan

OBJECTIVE

To provide baseline information about suspected vectors and the incidence, distribution and an active zone of transmission for cutaneous leishmaniasis (CL) in Chitral, Pakistan, using GIS tools; and to investigate the role of environmental factors in the disease dynamics.

METHOD

Two surveys in 2014 and 2016 as a basis for choropleth and environmental risk mapping.

RESULTS

A total of 769 captured specimens yielded 14 Phlebotomus and six Sergentomyia species including two potential vectors of CL, i.e. Phlebotomus papatasi and Phlebotomus sergenti. P. papatasi (71%) was dominant, followed by P. sergenti (18%). A choropleth map generated in Arcmap 10.1 based on 1560 CL case reports displayed maximum prevalence (0.92-2.5%) in Ayun, Broz, Charun, Chitral 1 and 2 and Darosh 1 and 2 union councils. An environmental risk map constructed by MaxEnt 3.3.3 defined an active zone of transmission based on leishmaniasis occurrence records (n = 315). The analysis of variable contribution in MaxEnt indicates significance of elevation (54.4%), population density (23.3%) and land use/land cover (6.6%) in CL disease dynamics.

CONCLUSION

The probability of CL increases (0.6-1 on logistic scale) in severely deforested areas, in lowland valleys and in regions with high-population density.

GridSample: an R package to generate household survey primary sampling units (PSUs) from gridded population data

BACKGROUND

Household survey data are collected by governments, international organizations, and companies to prioritize policies and allocate billions of dollars. Surveys are typically selected from recent census data; however, census data are often outdated or inaccurate. This paper describes how gridded population data might instead be used as a sample frame, and introduces the R GridSample algorithm for selecting primary sampling units (PSU) for complex household surveys with gridded population data. With a gridded population dataset and geographic boundary of the study area, GridSample allows a two-step process to sample "seed" cells with probability proportionate to estimated population size, then "grows" PSUs until a minimum population is achieved in each PSU. The algorithm permits stratification and oversampling of urban or rural areas. The approximately uniform size and shape of grid cells allows for spatial oversampling, not possible in typical surveys, possibly improving small area estimates with survey results.

RESULTS

We replicated the 2010 Rwanda Demographic and Health Survey (DHS) in GridSample by sampling the WorldPop 2010 UN-adjusted 100 m × 100 m gridded population dataset, stratifying by Rwanda's 30 districts, and oversampling in urban areas. The 2010 Rwanda DHS had 79 urban PSUs, 413 rural PSUs, with an average PSU population of 610 people. An equivalent sample in GridSample had 75 urban PSUs, 405 rural PSUs, and a median PSU population of 612 people. The number of PSUs differed because DHS added urban PSUs from specific districts while GridSample reallocated rural-to-urban PSUs across all districts.

CONCLUSIONS

Gridded population sampling is a promising alternative to typical census-based sampling when census data are moderately outdated or inaccurate. Four approaches to implementation have been tried: (1) using gridded PSU boundaries produced by GridSample, (2) manually segmenting gridded PSU using satellite imagery, (3) non-probability sampling (e.g. random-walk, "spin-the-pen"), and random sampling of households. Gridded population sampling is in its infancy, and further research is needed to assess the accuracy and feasibility of gridded population sampling. The GridSample R algorithm can be used to forward this research agenda.

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.

High resolution global gridded data for use in population studies

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

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

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

Global mapping of highly pathogenic avian influenza H5N1 and H5Nx clade 2.3.4.4 viruses with spatial cross-validation

Global disease suitability models are essential tools to inform surveillance systems and enable early detection. We present the first global suitability model of highly pathogenic avian influenza (HPAI) H5N1 and demonstrate that reliable predictions can be obtained at global scale. Best predictions are obtained using spatial predictor variables describing host distributions, rather than land use or eco-climatic spatial predictor variables, with a strong association with domestic duck and extensively raised chicken densities. Our results also support a more systematic use of spatial cross-validation in large-scale disease suitability modelling compared to standard random cross-validation that can lead to unreliable measure of extrapolation accuracy. A global suitability model of the H5 clade 2.3.4.4 viruses, a group of viruses that recently spread extensively in Asia and the US, shows in comparison a lower spatial extrapolation capacity than the HPAI H5N1 models, with a stronger association with intensively raised chicken densities and anthropogenic factors.

Faster Detection of Poliomyelitis Outbreaks to Support Polio Eradication

As the global eradication of poliomyelitis approaches the final stages, prompt detection of new outbreaks is critical to enable a fast and effective outbreak response. Surveillance relies on reporting of acute flaccid paralysis (AFP) cases and laboratory confirmation through isolation of poliovirus from stool. However, delayed sample collection and testing can delay outbreak detection. We investigated whether weekly testing for clusters of AFP by location and time, using the Kulldorff scan statistic, could provide an early warning for outbreaks in 20 countries. A mixed-effects regression model was used to predict background rates of nonpolio AFP at the district level. In Tajikistan and Congo, testing for AFP clusters would have resulted in an outbreak warning 39 and 11 days, respectively, before official confirmation of large outbreaks. This method has relatively high specificity and could be integrated into the current polio information system to support rapid outbreak response activities.

Spatial analysis and characteristics of pig farming in Thailand

Background

In Thailand, pig production intensified significantly during the last decade, with many economic, epidemiological and environmental implications. Strategies toward more sustainable future developments are currently investigated, and these could be informed by a detailed assessment of the main trends in the pig sector, and on how different production systems are geographically distributed. This study had two main objectives. First, we aimed to describe the main trends and geographic patterns of pig production systems in Thailand in terms of pig type (native, breeding, and fattening pigs), farm scales (smallholder and large-scale farming systems) and type of farming systems (farrow-to-finish, nursery, and finishing systems) based on a very detailed 2010 census. Second, we aimed to study the statistical spatial association between these different types of pig farming distribution and a set of spatial variables describing access to feed and markets.

Results

Over the last decades, pig population gradually increased, with a continuously increasing number of pigs per holder, suggesting a continuing intensification of the sector. The different pig-production systems showed very contrasted geographical distributions. The spatial distribution of large-scale pig farms corresponds with that of commercial pig breeds, and spatial analysis conducted using Random Forest distribution models indicated that these were concentrated in lowland urban or peri-urban areas, close to means of transportation, facilitating supply to major markets such as provincial capitals and the Bangkok Metropolitan region. Conversely the smallholders were distributed throughout the country, with higher densities located in highland, remote, and rural areas, where they supply local rural markets. A limitation of the study was that pig farming systems were defined from the number of animals per farm, resulting in their possible misclassification, but this should have a limited impact on the main patterns revealed by the analysis.

Conclusions

The very contrasted distribution of different pig production systems present opportunities for future regionalization of pig production. More specifically, the detailed geographical analysis of the different production systems will be used to spatially-inform planning decisions for pig farming accounting for the specific health, environment and economical implications of the different pig production systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-016-0849-7) contains supplementary material, which is available to authorized users.

Mapping internal connectivity through human migration in malaria endemic countries

Human mobility continues to increase in terms of volumes and reach, producing growing global connectivity. This connectivity hampers efforts to eliminate infectious diseases such as malaria through reintroductions of pathogens, and thus accounting for it becomes important in designing global, continental, regional, and national strategies. Recent works have shown that census-derived migration data provides a good proxy for internal connectivity, in terms of relative strengths of movement between administrative units, across temporal scales. To support global malaria eradication strategy efforts, here we describe the construction of an open access archive of estimated internal migration flows in endemic countries built through pooling of census microdata. These connectivity datasets, described here along with the approaches and methods used to create and validate them, are available both through the WorldPop website and the WorldPop Dataverse Repository.

Clade-level Spatial Modelling of HPAI H5N1 Dynamics in the Mekong Region Reveals New Patterns and Associations with Agro-Ecological Factors

The highly pathogenic avian influenza (HPAI) H5N1 virus has been circulating in Asia since 2003 and diversified into several genetic lineages, or clades. Although the spatial distribution of its outbreaks was extensively studied, differences in clades were never previously taken into account. We developed models to quantify associations over time and space between different HPAI H5N1 viruses from clade 1, 2.3.4 and 2.3.2 and agro-ecological factors. We found that the distribution of clades in the Mekong region from 2004 to 2013 was strongly regionalised, defining specific epidemiological zones, or epizones. Clade 1 became entrenched in the Mekong Delta and was not supplanted by newer clades, in association with a relatively higher presence of domestic ducks. In contrast, two new clades were introduced (2.3.4 and 2.3.2) in northern Viet Nam and were associated with higher chicken density and more intensive chicken production systems. We suggest that differences in poultry production systems in these different epizones may explain these associations, along with differences in introduction pressure from neighbouring countries. The different distribution patterns found at the clade level would not be otherwise apparent through analysis treating all outbreaks equally, which requires improved linking of disease outbreak records and genetic sequence data.

Population Estimation Using a 3D City Model: A Multi-Scale Country-Wide Study in the Netherlands

The remote estimation of a region’s population has for decades been a key application of geographic information science in demography. Most studies have used 2D data (maps, satellite imagery) to estimate population avoiding field surveys and questionnaires. As the availability of semantic 3D city models is constantly increasing, we investigate to what extent they can be used for the same purpose. Based on the assumption that housing space is a proxy for the number of its residents, we use two methods to estimate the population with 3D city models in two directions: (1) disaggregation (areal interpolation) to estimate the population of small administrative entities (e.g. neighbourhoods) from that of larger ones (e.g. municipalities); and (2) a statistical modelling approach to estimate the population of large entities from a sample composed of their smaller ones (e.g. one acquired by a government register). Starting from a complete Dutch census dataset at the neighbourhood level and a 3D model of all 9.9 million buildings in the Netherlands, we compare the population estimates obtained by both methods with the actual population as reported in the census, and use it to evaluate the quality that can be achieved by estimations at different administrative levels. We also analyse how the volume-based estimation enabled by 3D city models fares in comparison to 2D methods using building footprints and floor areas, as well as how it is affected by different levels of semantic detail in a 3D city model. We conclude that 3D city models are useful for estimations of large areas (e.g. for a country), and that the 3D approach has clear advantages over the 2D approach.