The success of vaccination programs depends largely on the mechanisms used in vaccine delivery. National immunization programs offer childhood vaccines through fixed and outreach services within the health system and often, additional supplementary immunization activities (SIAs) are undertaken to fill gaps and boost coverage. Here, we map predicted coverage at 1 × 1 km spatial resolution in five low- and middle-income countries to identify areas that are under-vaccinated via each delivery method using Demographic and Health Surveys data. We compare estimates of the coverage of the third dose of diphtheria-tetanus-pertussis-containing vaccine (DTP3), which is typically delivered through routine immunization (RI), with those of measles-containing vaccine (MCV) for which SIAs are also undertaken. We find that SIAs have boosted MCV coverage in some places, but not in others, particularly where RI had been deficient, as depicted by DTP coverage. The modelling approaches outlined here can help to guide geographical prioritization and strategy design.

The expansion of childhood vaccination programs in low and middle income countries has been a substantial public health success story. Indicators of the performance of intervention programmes such as coverage levels and numbers covered are typically measured through national statistics or at the scale of large regions due to survey design, administrative convenience or operational limitations. These mask heterogeneities and 'coldspots' of low coverage that may allow diseases to persist, even if overall coverage is high. Hence, to decrease inequities and accelerate progress towards disease elimination goals, fine-scale variation in coverage should be better characterized. Using measles as an example, cluster-level Demographic and Health Surveys (DHS) data were used to map vaccination coverage at 1 km spatial resolution in Cambodia, Mozambique and Nigeria for varying age-group categories of children under five years, using Bayesian geostatistical techniques built on a suite of publicly available geospatial covariates and implemented via Markov Chain Monte Carlo (MCMC) methods. Measles vaccination coverage was found to be strongly predicted by just 4-5 covariates in geostatistical models, with remoteness consistently selected as a key variable. The output 1 × 1 km maps revealed significant heterogeneities within the three countries that were not captured using province-level summaries. Integration with population data showed that at the time of the surveys, few districts attained the 80% coverage, that is one component of the WHO Global Vaccine Action Plan 2020 targets. The elimination of vaccine-preventable diseases requires a strong evidence base to guide strategies and inform efficient use of limited resources. The approaches outlined here provide a route to moving beyond large area summaries of vaccination coverage that mask epidemiologically-important heterogeneities to detailed maps that capture subnational vulnerabilities. The output datasets are built on open data and methods, and in flexible format that can be aggregated to more operationally-relevant administrative unit levels.

Measles vaccination campaigns are conducted regularly in many low- and middle-income countries to boost measles control efforts and accelerate progress towards elimination. National and sometimes first-level administrative division campaign coverage may be estimated through post-campaign coverage surveys (PCCS). However, these large-area estimates mask significant geographic inequities in coverage at more granular levels. Here, we undertake a geospatial analysis of the Nigeria 2017-18 PCCS data to produce coverage estimates at 1 × 1 km resolution and the district level using binomial spatial regression models built on a suite of geospatial covariates and implemented in a Bayesian framework via the INLA-SPDE approach. We investigate the individual and combined performance of the campaign and routine immunization (RI) by mapping various indicators of coverage for children aged 9-59 months. Additionally, we compare estimated coverage before the campaign at 1 × 1 km and the district level with predicted coverage maps produced using other surveys conducted in 2013 and 2016-17. Coverage during the campaign was generally higher and more homogeneous than RI coverage but geospatial differences in the campaign's reach of previously unvaccinated children are shown. Persistent areas of low coverage highlight the need for improved RI performance. The results can help to guide the conduct of future campaigns, improve vaccination monitoring and measles elimination efforts. Moreover, the approaches used here can be readily extended to other countries.

Achieving equity in vaccination coverage has been a critical priority within the global health community. Despite increased efforts recently, certain populations still have a high proportion of un- and under-vaccinated children in many low- and middle-income countries (LMICs). These populations are often assumed to reside in remote-rural areas, urban slums and conflict-affected areas. Here, we investigate the effects of these key community-level factors, alongside a wide range of other individual, household and community level factors, on vaccination coverage. Using geospatial datasets, including cross-sectional data from the most recent Demographic and Health Surveys conducted between 2008 and 2018 in nine LMICs, we fitted Bayesian multi-level binary logistic regression models to determine key community-level and other factors significantly associated with non- and under-vaccination. We analyzed the odds of receipt of the first doses of diphtheria-tetanus-pertussis (DTP1) vaccine and measles-containing vaccine (MCV1), and receipt of all three recommended DTP doses (DTP3) independently, in children aged 12–23 months. In bivariate analyses, we found that remoteness increased the odds of non- and under-vaccination in nearly all the study countries. We also found evidence that living in conflict and urban slum areas reduced the odds of vaccination, but not in most cases as expected. However, the odds of vaccination were more likely to be lower in urban slums than formal urban areas. Our multivariate analyses revealed that the key community variables–remoteness, conflict and urban slum–were sometimes associated with non- and under-vaccination, but they were not frequently predictors of these outcomes after controlling for other factors. Individual and household factors such as maternal utilization of health services, maternal education and ethnicity, were more common predictors of vaccination. Reaching the Immunisation Agenda 2030 target of reducing the number of zero-dose children by 50% by 2030 will require country tailored analyses and strategies to identify and reach missed communities with reliable immunisation services.

After many decades of vaccination, measles epidemiology varies greatly between and within countries. National immunization programs are therefore encouraged to conduct regular situation analyses and to leverage models to adapt interventions to local needs. Here, we review applications of models to develop locally tailored interventions to support control and elimination efforts. In general, statistical and semi-mechanistic transmission models can be used to synthesize information from vaccination coverage, measles incidence, demographic, and/or serological data, offering a means to estimate the spatial and age-specific distribution of measles susceptibility. These estimates complete the picture provided by vaccination coverage alone, by accounting for natural immunity. Dynamic transmission models can then be used to evaluate the relative impact of candidate interventions for measles control and elimination and the expected future epidemiology. In most countries, models predict substantial numbers of susceptible individuals outside the age range of routine vaccination, which affects outbreak risk and necessitates additional intervention to achieve elimination. More effective use of models to inform both vaccination program planning and evaluation requires the development of training to enhance broader understanding of models and where feasible, building capacity for modelling in-country, pipelines for rapid evaluation of model predictions using surveillance data, and clear protocols for incorporating model results into decision-making.

Geographically precise identification and targeting of populations at risk of vaccine-preventable diseases has gained renewed attention within the global health community over the last few years. District level estimates of vaccination coverage and corresponding zero-dose prevalence constitute a potentially useful evidence base to evaluate the performance of vaccination strategies. These estimates are also valuable for identifying missed communities, hence enabling targeted interventions and better resource allocation. Here, we fit Bayesian geostatistical models to map the routine coverage of the first doses of diphtheria-tetanus-pertussis vaccine (DTP1) and measles-containing vaccine (MCV1) and corresponding zero-dose estimates in Nigeria at 1x1 km resolution and the district level using geospatial data sets. We also map MCV1 coverage before and after the 2019 measles vaccination campaign in the northern states to further explore variations in routine vaccine coverage and to evaluate the effectiveness of both routine immunization (RI) and campaigns in reaching zero-dose children. Additionally, we map the spatial distributions of reported measles cases during 2018 to 2020 and explore their relationships with MCV zero-dose prevalence to highlight the public health implications of varying performance of vaccination strategies across the country. Our analysis revealed strong similarities between the spatial distributions of DTP and MCV zero dose prevalence, with districts with the highest prevalence concentrated mostly in the northwest and the northeast, but also in other areas such as Lagos state and the Federal Capital Territory. Although the 2019 campaign reduced MCV zero-dose prevalence substantially in the north, pockets of vulnerabilities remained in areas that had among the highest prevalence prior to the campaign. Importantly, we found strong correlations between measles case counts and MCV RI zero-dose estimates, which provides a strong indication that measles incidence in the country is mostly affected by RI coverage. Our analyses reveal an urgent and highly significant need to strengthen the country's RI program as a longer-term measure for disease control, whilst ensuring effective campaigns in the short term.