Global, regional, and national age-sex-specific mortality and life expectancy, 1950–2017: a systematic analysis for the Global Burden of Disease Study 2017

BACKGROUND: Assessments of age-specific mortality and life expectancy have been done by the UN Population Division, Department of Economics and Social Affairs (UNPOP), the United States Census Bureau, WHO, and as part of previous iterations of the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD). Previous iterations of the GBD used population estimates from UNPOP, which were not derived in a way that was internally consistent with the estimates of the numbers of deaths in the GBD. The present iteration of the GBD, GBD 2017, improves on previous assessments and provides timely estimates of the mortality experience of populations globally. METHODS: The GBD uses all available data to produce estimates of mortality rates between 1950 and 2017 for 23 age groups, both sexes, and 918 locations, including 195 countries and territories and subnational locations for 16 countries. Data used include vital registration systems, sample registration systems, household surveys (complete birth histories, summary birth histories, sibling histories), censuses (summary birth histories, household deaths), and Demographic Surveillance Sites. In total, this analysis used 8259 data sources. Estimates of the probability of death between birth and the age of 5 years and between ages 15 and 60 years are generated and then input into a model life table system to produce complete life tables for all locations and years. Fatal discontinuities and mortality due to HIV/AIDS are analysed separately and then incorporated into the estimation. We analyse the relationship between age-specific mortality and development status using the Socio-demographic Index, a composite measure based on fertility under the age of 25 years, education, and income. There are four main methodological improvements in GBD 2017 compared with GBD 2016: 622 additional data sources have been incorporated; new estimates of population, generated by the GBD study, are used; statistical methods used in different components of the analysis have been further standardised and improved; and the analysis has been extended backwards in time by two decades to start in 1950. FINDINGS: Globally, 18·7% (95% uncertainty interval 18·4–19·0) of deaths were registered in 1950 and that proportion has been steadily increasing since, with 58·8% (58·2–59·3) of all deaths being registered in 2015. At the global level, between 1950 and 2017, life expectancy increased from 48·1 years (46·5–49·6) to 70·5 years (70·1–70·8) for men and from 52·9 years (51·7–54·0) to 75·6 years (75·3–75·9) for women. Despite this overall progress, there remains substantial variation in life expectancy at birth in 2017, which ranges from 49·1 years (46·5–51·7) for men in the Central African Republic to 87·6 years (86·9–88·1) among women in Singapore. The greatest progress across age groups was for children younger than 5 years; under-5 mortality dropped from 216·0 deaths (196·3–238·1) per 1000 livebirths in 1950 to 38·9 deaths (35·6–42·83) per 1000 livebirths in 2017, with huge reductions across countries. Nevertheless, there were still 5·4 million (5·2–5·6) deaths among children younger than 5 years in the world in 2017. Progress has been less pronounced and more variable for adults, especially for adult males, who had stagnant or increasing mortality rates in several countries. The gap between male and female life expectancy between 1950 and 2017, while relatively stable at the global level, shows distinctive patterns across super-regions and has consistently been the largest in central Europe, eastern Europe, and central Asia, and smallest in south Asia. Performance was also variable across countries and time in observed mortality rates compared with those expected on the basis of development. INTERPRETATION: This analysis of age-sex-specific mortality shows that there are remarkably complex patterns in population mortality across countries. The findings of this study highlight global successes, such as the large decline in under-5 mortality, which reflects significant local, national, and global commitment and investment over several decades. However, they also bring attention to mortality patterns that are a cause for concern, particularly among adult men and, to a lesser extent, women, whose mortality rates have stagnated in many countries over the time period of this study, and in some cases are increasing

Estimates of global, regional, and national incidence, prevalence, and mortality of HIV, 1980�2015: the Global Burden of Disease Study 2015

Background Timely assessment of the burden of HIV/AIDS is essential for policy setting and programme evaluation. In this report from the Global Burden of Disease Study 2015 (GBD 2015), we provide national estimates of levels and trends of HIV/AIDS incidence, prevalence, coverage of antiretroviral therapy (ART), and mortality for 195 countries and territories from 1980 to 2015. Methods For countries without high-quality vital registration data, we estimated prevalence and incidence with data from antenatal care clinics and population-based seroprevalence surveys, and with assumptions by age and sex on initial CD4 distribution at infection, CD4 progression rates (probability of progression from higher to lower CD4 cell-count category), on and off antiretroviral therapy (ART) mortality, and mortality from all other causes. Our estimation strategy links the GBD 2015 assessment of all-cause mortality and estimation of incidence and prevalence so that for each draw from the uncertainty distribution all assumptions used in each step are internally consistent. We estimated incidence, prevalence, and death with GBD versions of the Estimation and Projection Package (EPP) and Spectrum software originally developed by the Joint United Nations Programme on HIV/AIDS (UNAIDS). We used an open-source version of EPP and recoded Spectrum for speed, and used updated assumptions from systematic reviews of the literature and GBD demographic data. For countries with high-quality vital registration data, we developed the cohort incidence bias adjustment model to estimate HIV incidence and prevalence largely from the number of deaths caused by HIV recorded in cause-of-death statistics. We corrected these statistics for garbage coding and HIV misclassification. Findings Global HIV incidence reached its peak in 1997, at 3·3 million new infections (95 uncertainty interval UI 3·1�3·4 million). Annual incidence has stayed relatively constant at about 2·6 million per year (range 2·5�2·8 million) since 2005, after a period of fast decline between 1997 and 2005. The number of people living with HIV/AIDS has been steadily increasing and reached 38·8 million (95% UI 37·6�40·4 million) in 2015. At the same time, HIV/AIDS mortality has been declining at a steady pace, from a peak of 1·8 million deaths (95% UI 1·7�1·9 million) in 2005, to 1·2 million deaths (1·1�1·3 million) in 2015. We recorded substantial heterogeneity in the levels and trends of HIV/AIDS across countries. Although many countries have experienced decreases in HIV/AIDS mortality and in annual new infections, other countries have had slowdowns or increases in rates of change in annual new infections. Interpretation Scale-up of ART and prevention of mother-to-child transmission has been one of the great successes of global health in the past two decades. However, in the past decade, progress in reducing new infections has been slow, development assistance for health devoted to HIV has stagnated, and resources for health in low-income countries have grown slowly. Achievement of the new ambitious goals for HIV enshrined in Sustainable Development Goal 3 and the 90-90-90 UNAIDS targets will be challenging, and will need continued efforts from governments and international agencies in the next 15 years to end AIDS by 2030. Funding Bill & Melinda Gates Foundation, and National Institute of Mental Health and National Institute on Aging, National Institutes of Health. © 2016 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY licens

Global, regional, national, and selected subnational levels of stillbirths, neonatal, infant, and under-5 mortality, 1980�2015: a systematic analysis for the Global Burden of Disease Study 2015

Background Established in 2000, Millennium Development Goal 4 (MDG4) catalysed extraordinary political, financial, and social commitments to reduce under-5 mortality by two-thirds between 1990 and 2015. At the country level, the pace of progress in improving child survival has varied markedly, highlighting a crucial need to further examine potential drivers of accelerated or slowed decreases in child mortality. The Global Burden of Disease 2015 Study (GBD 2015) provides an analytical framework to comprehensively assess these trends for under-5 mortality, age-specific and cause-specific mortality among children under 5 years, and stillbirths by geography over time. Methods Drawing from analytical approaches developed and refined in previous iterations of the GBD study, we generated updated estimates of child mortality by age group (neonatal, post-neonatal, ages 1�4 years, and under 5) for 195 countries and territories and selected subnational geographies, from 1980�2015. We also estimated numbers and rates of stillbirths for these geographies and years. Gaussian process regression with data source adjustments for sampling and non-sampling bias was applied to synthesise input data for under-5 mortality for each geography. Age-specific mortality estimates were generated through a two-stage age�sex splitting process, and stillbirth estimates were produced with a mixed-effects model, which accounted for variable stillbirth definitions and data source-specific biases. For GBD 2015, we did a series of novel analyses to systematically quantify the drivers of trends in child mortality across geographies. First, we assessed observed and expected levels and annualised rates of decrease for under-5 mortality and stillbirths as they related to the Soci-demographic Index (SDI). Second, we examined the ratio of recorded and expected levels of child mortality, on the basis of SDI, across geographies, as well as differences in recorded and expected annualised rates of change for under-5 mortality. Third, we analysed levels and cause compositions of under-5 mortality, across time and geographies, as they related to rising SDI. Finally, we decomposed the changes in under-5 mortality to changes in SDI at the global level, as well as changes in leading causes of under-5 deaths for countries and territories. We documented each step of the GBD 2015 child mortality estimation process, as well as data sources, in accordance with the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER). Findings Globally, 5·8 million (95 uncertainty interval UI 5·7�6·0) children younger than 5 years died in 2015, representing a 52·0% (95% UI 50·7�53·3) decrease in the number of under-5 deaths since 1990. Neonatal deaths and stillbirths fell at a slower pace since 1990, decreasing by 42·4% (41·3�43·6) to 2·6 million (2·6�2·7) neonatal deaths and 47·0% (35·1�57·0) to 2·1 million (1·8-2·5) stillbirths in 2015. Between 1990 and 2015, global under-5 mortality decreased at an annualised rate of decrease of 3·0% (2·6�3·3), falling short of the 4·4% annualised rate of decrease required to achieve MDG4. During this time, 58 countries met or exceeded the pace of progress required to meet MDG4. Between 2000, the year MDG4 was formally enacted, and 2015, 28 additional countries that did not achieve the 4·4% rate of decrease from 1990 met the MDG4 pace of decrease. However, absolute levels of under-5 mortality remained high in many countries, with 11 countries still recording rates exceeding 100 per 1000 livebirths in 2015. Marked decreases in under-5 deaths due to a number of communicable diseases, including lower respiratory infections, diarrhoeal diseases, measles, and malaria, accounted for much of the progress in lowering overall under-5 mortality in low-income countries. Compared with gains achieved for infectious diseases and nutritional deficiencies, the persisting toll of neonatal conditions and congenital anomalies on child survival became evident, especially in low-income and low-middle-income countries. We found sizeable heterogeneities in comparing observed and expected rates of under-5 mortality, as well as differences in observed and expected rates of change for under-5 mortality. At the global level, we recorded a divergence in observed and expected levels of under-5 mortality starting in 2000, with the observed trend falling much faster than what was expected based on SDI through 2015. Between 2000 and 2015, the world recorded 10·3 million fewer under-5 deaths than expected on the basis of improving SDI alone. Interpretation Gains in child survival have been large, widespread, and in many places in the world, faster than what was anticipated based on improving levels of development. Yet some countries, particularly in sub-Saharan Africa, still had high rates of under-5 mortality in 2015. Unless these countries are able to accelerate reductions in child deaths at an extraordinary pace, their achievement of proposed SDG targets is unlikely. Improving the evidence base on drivers that might hasten the pace of progress for child survival, ranging from cost-effective intervention packages to innovative financing mechanisms, is vital to charting the pathways for ultimately ending preventable child deaths by 2030. Funding Bill & Melinda Gates Foundation. © 2016 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY license

Global, regional, and national disability-adjusted life-years (DALYs) for 315 diseases and injuries and healthy life expectancy (HALE), 1990�2015: a systematic analysis for the Global Burden of Disease Study 2015

Background Healthy life expectancy (HALE) and disability-adjusted life-years (DALYs) provide summary measures of health across geographies and time that can inform assessments of epidemiological patterns and health system performance, help to prioritise investments in research and development, and monitor progress toward the Sustainable Development Goals (SDGs). We aimed to provide updated HALE and DALYs for geographies worldwide and evaluate how disease burden changes with development. Methods We used results from the Global Burden of Diseases, Injuries, and Risk Factors Study 2015 (GBD 2015) for all-cause mortality, cause-specific mortality, and non-fatal disease burden to derive HALE and DALYs by sex for 195 countries and territories from 1990 to 2015. We calculated DALYs by summing years of life lost (YLLs) and years of life lived with disability (YLDs) for each geography, age group, sex, and year. We estimated HALE using the Sullivan method, which draws from age-specific death rates and YLDs per capita. We then assessed how observed levels of DALYs and HALE differed from expected trends calculated with the Socio-demographic Index (SDI), a composite indicator constructed from measures of income per capita, average years of schooling, and total fertility rate. Findings Total global DALYs remained largely unchanged from 1990 to 2015, with decreases in communicable, neonatal, maternal, and nutritional (Group 1) disease DALYs offset by increased DALYs due to non-communicable diseases (NCDs). Much of this epidemiological transition was caused by changes in population growth and ageing, but it was accelerated by widespread improvements in SDI that also correlated strongly with the increasing importance of NCDs. Both total DALYs and age-standardised DALY rates due to most Group 1 causes significantly decreased by 2015, and although total burden climbed for the majority of NCDs, age-standardised DALY rates due to NCDs declined. Nonetheless, age-standardised DALY rates due to several high-burden NCDs (including osteoarthritis, drug use disorders, depression, diabetes, congenital birth defects, and skin, oral, and sense organ diseases) either increased or remained unchanged, leading to increases in their relative ranking in many geographies. From 2005 to 2015, HALE at birth increased by an average of 2·9 years (95 uncertainty interval 2·9�3·0) for men and 3·5 years (3·4�3·7) for women, while HALE at age 65 years improved by 0·85 years (0·78�0·92) and 1·2 years (1·1�1·3), respectively. Rising SDI was associated with consistently higher HALE and a somewhat smaller proportion of life spent with functional health loss; however, rising SDI was related to increases in total disability. Many countries and territories in central America and eastern sub-Saharan Africa had increasingly lower rates of disease burden than expected given their SDI. At the same time, a subset of geographies recorded a growing gap between observed and expected levels of DALYs, a trend driven mainly by rising burden due to war, interpersonal violence, and various NCDs. Interpretation Health is improving globally, but this means more populations are spending more time with functional health loss, an absolute expansion of morbidity. The proportion of life spent in ill health decreases somewhat with increasing SDI, a relative compression of morbidity, which supports continued efforts to elevate personal income, improve education, and limit fertility. Our analysis of DALYs and HALE and their relationship to SDI represents a robust framework on which to benchmark geography-specific health performance and SDG progress. Country-specific drivers of disease burden, particularly for causes with higher-than-expected DALYs, should inform financial and research investments, prevention efforts, health policies, and health system improvement initiatives for all countries along the development continuum. Funding Bill & Melinda Gates Foundation. © 2016 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY licens

Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990�2015: a systematic analysis for the Global Burden of Disease Study 2015

Background The Global Burden of Diseases, Injuries, and Risk Factors Study 2015 provides an up-to-date synthesis of the evidence for risk factor exposure and the attributable burden of disease. By providing national and subnational assessments spanning the past 25 years, this study can inform debates on the importance of addressing risks in context. Methods We used the comparative risk assessment framework developed for previous iterations of the Global Burden of Disease Study to estimate attributable deaths, disability-adjusted life-years (DALYs), and trends in exposure by age group, sex, year, and geography for 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks from 1990 to 2015. This study included 388 risk-outcome pairs that met World Cancer Research Fund-defined criteria for convincing or probable evidence. We extracted relative risk and exposure estimates from randomised controlled trials, cohorts, pooled cohorts, household surveys, census data, satellite data, and other sources. We used statistical models to pool data, adjust for bias, and incorporate covariates. We developed a metric that allows comparisons of exposure across risk factors�the summary exposure value. Using the counterfactual scenario of theoretical minimum risk level, we estimated the portion of deaths and DALYs that could be attributed to a given risk. We decomposed trends in attributable burden into contributions from population growth, population age structure, risk exposure, and risk-deleted cause-specific DALY rates. We characterised risk exposure in relation to a Socio-demographic Index (SDI). Findings Between 1990 and 2015, global exposure to unsafe sanitation, household air pollution, childhood underweight, childhood stunting, and smoking each decreased by more than 25. Global exposure for several occupational risks, high body-mass index (BMI), and drug use increased by more than 25 over the same period. All risks jointly evaluated in 2015 accounted for 57·8 (95 CI 56·6�58·8) of global deaths and 41·2 (39·8�42·8) of DALYs. In 2015, the ten largest contributors to global DALYs among Level 3 risks were high systolic blood pressure (211·8 million 192·7 million to 231·1 million global DALYs), smoking (148·6 million 134·2 million to 163·1 million), high fasting plasma glucose (143·1 million 125·1 million to 163·5 million), high BMI (120·1 million 83·8 million to 158·4 million), childhood undernutrition (113·3 million 103·9 million to 123·4 million), ambient particulate matter (103·1 million 90·8 million to 115·1 million), high total cholesterol (88·7 million 74·6 million to 105·7 million), household air pollution (85·6 million 66·7 million to 106·1 million), alcohol use (85·0 million 77·2 million to 93·0 million), and diets high in sodium (83·0 million 49·3 million to 127·5 million). From 1990 to 2015, attributable DALYs declined for micronutrient deficiencies, childhood undernutrition, unsafe sanitation and water, and household air pollution; reductions in risk-deleted DALY rates rather than reductions in exposure drove these declines. Rising exposure contributed to notable increases in attributable DALYs from high BMI, high fasting plasma glucose, occupational carcinogens, and drug use. Environmental risks and childhood undernutrition declined steadily with SDI; low physical activity, high BMI, and high fasting plasma glucose increased with SDI. In 119 countries, metabolic risks, such as high BMI and fasting plasma glucose, contributed the most attributable DALYs in 2015. Regionally, smoking still ranked among the leading five risk factors for attributable DALYs in 109 countries; childhood underweight and unsafe sex remained primary drivers of early death and disability in much of sub-Saharan Africa. Interpretation Declines in some key environmental risks have contributed to declines in critical infectious diseases. Some risks appear to be invariant to SDI. Increasing risks, including high BMI, high fasting plasma glucose, drug use, and some occupational exposures, contribute to rising burden from some conditions, but also provide opportunities for intervention. Some highly preventable risks, such as smoking, remain major causes of attributable DALYs, even as exposure is declining. Public policy makers need to pay attention to the risks that are increasingly major contributors to global burden. Funding Bill & Melinda Gates Foundation. © 2016 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY licens

Success of post-fire plant recovery strategies varies with shifting fire seasonality

Wildfires are increasing in size and severity and fire seasons are lengthening, largely driven by climate and land-use change. Many plant species from fire-prone ecosystems are adapted to specific fire regimes corresponding to historical conditions and shifts beyond these bounds may have severe impacts on vegetation recovery and long-term species persistence. Here, we conduct a meta-analysis of field-based studies across different vegetation types and climate regions to investigate how post-fire plant recruitment, reproduction and survival are affected by fires that occur outside of the historical fire season. We find that fires outside of the historical fire season may lead to decreased post-fire recruitment, particularly in obligate seeding species. Conversely, we find a general increase in post-fire survival in resprouting species. Our results highlight the trade-offs that exist when considering the effects of changes in the seasonal timing of fire, an already present aspect of climate-related fire regime change

Mapping routine measles vaccination in low- and middle-income countries

The safe, highly effective measles vaccine has been recommended globally since 1974, yet in 2017 there were more than 17 million cases of measles and 83,400 deaths in children under 5 years old, and more than 99% of both occurred in low- and middle-income countries (LMICs)(1-4). Globally comparable, annual, local estimates of routine first-dose measles-containing vaccine (MCV1) coverage are critical for understanding geographically precise immunity patterns, progress towards the targets of the Global Vaccine Action Plan (GVAP), and high-risk areas amid disruptions to vaccination programmes caused by coronavirus disease 2019 (COVID-19)(5-8). Here we generated annual estimates of routine childhood MCV1 coverage at 5 x 5-km(2) pixel and second administrative levels from 2000 to 2019 in 101 LMICs, quantified geographical inequality and assessed vaccination status by geographical remoteness. After widespread MCV1 gains from 2000 to 2010, coverage regressed in more than half of the districts between 2010 and 2019, leaving many LMICs far from the GVAP goal of 80% coverage in all districts by 2019. MCV1 coverage was lower in rural than in urban locations, although a larger proportion of unvaccinated children overall lived in urban locations; strategies to provide essential vaccination services should address both geographical contexts. These results provide a tool for decision-makers to strengthen routine MCV1 immunization programmes and provide equitable disease protection for all children.Peer reviewe

Global injury morbidity and mortality from 1990 to 2017: Results from the global burden of disease study 2017

Background Past research in population health trends has shown that injuries form a substantial burden of population health loss. Regular updates to injury burden assessments are critical. We report Global Burden of Disease (GBD) 2017 Study estimates on morbidity and mortality for all injuries. methods We reviewed results for injuries from the GBD 2017 study. GBD 2017 measured injury-specific mortality and years of life lost (YLLs) using the Cause of Death Ensemble model. To measure non-fatal injuries, GBD 2017 modelled injury-specific incidence and converted this to prevalence and years lived with disability (YLDs). YLLs and YLDs were summed to calculate disability-adjusted life years (DALYs). Findings In 1990, there were 4 260 493 (4 085 700 to 4 396 138) injury deaths, which increased to 4 484 722 (4 332 010 to 4 585 554) deaths in 2017, while age-standardised mortality decreased from 1079 (1073 to 1086) to 738 (730 to 745) per 100 000. In 1990, there were 354 064 302 (95% uncertainty interval: 338 174 876 to 371 610 802) new cases of injury globally, which increased to 520 710 288 (493 430 247 to 547 988 635) new cases in 2017. During this time, age-standardised incidence decreased non-significantly from 6824 (6534 to 7147) to 6763 (6412 to 7118) per 100 000. Between 1990 and 2017, age-standardised DALYs decreased from 4947 (4655 to 5233) per 100 000 to 3267 (3058 to 3505). Interpretation Injuries are an important cause of health loss globally, though mortality has declined between 1990 and 2017. Future research in injury burden should focus on prevention in high-burden populations, improving data collection and ensuring access to medical care. © Author(s) (or their employer(s)) 2020. Re-use permitted under CC BY. Published by BMJ

Estimating global injuries morbidity and mortality: Methods and data used in the Global Burden of Disease 2017 study

Background: While there is a long history of measuring death and disability from injuries, modern research methods must account for the wide spectrum of disability that can occur in an injury, and must provide estimates with sufficient demographic, geographical and temporal detail to be useful for policy makers. The Global Burden of Disease (GBD) 2017 study used methods to provide highly detailed estimates of global injury burden that meet these criteria. Methods: In this study, we report and discuss the methods used in GBD 2017 for injury morbidity and mortality burden estimation. In summary, these methods included estimating cause-specific mortality for every cause of injury, and then estimating incidence for every cause of injury. Non-fatal disability for each cause is then calculated based on the probabilities of suffering from different types of bodily injury experienced. Results: GBD 2017 produced morbidity and mortality estimates for 38 causes of injury. Estimates were produced in terms of incidence, prevalence, years lived with disability, cause-specific mortality, years of life lost and disability-adjusted life-years for a 28-year period for 22 age groups, 195 countries and both sexes. Conclusions: GBD 2017 demonstrated a complex and sophisticated series of analytical steps using the largest known database of morbidity and mortality data on injuries. GBD 2017 results should be used to help inform injury prevention policy making and resource allocation. We also identify important avenues for improving injury burden estimation in the future. © Author(s) (or their employer(s)) 2020. Re-use permitted under CC B Y. Published by BMJ

Mapping routine measles vaccination in low- and middle-income countries

The safe, highly effective measles vaccine has been recommended globally since 1974, yet in 2017 there were more than 17 million cases of measles and 83,400 deaths in children under 5 years old, and more than 99% of both occurred in low- and middle-income countries (LMICs)1,2,3,4. Globally comparable, annual, local estimates of routine first-dose measles-containing vaccine (MCV1) coverage are critical for understanding geographically precise immunity patterns, progress towards the targets of the Global Vaccine Action Plan (GVAP), and high-risk areas amid disruptions to vaccination programmes caused by coronavirus disease 2019 (COVID-19)5,6,7,8. Here we generated annual estimates of routine childhood MCV1 coverage at 5 × 5-km2 pixel and second administrative levels from 2000 to 2019 in 101 LMICs, quantified geographical inequality and assessed vaccination status by geographical remoteness. After widespread MCV1 gains from 2000 to 2010, coverage regressed in more than half of the districts between 2010 and 2019, leaving many LMICs far from the GVAP goal of 80% coverage in all districts by 2019. MCV1 coverage was lower in rural than in urban locations, although a larger proportion of unvaccinated children overall lived in urban locations; strategies to provide essential vaccination services should address both geographical contexts. These results provide a tool for decision-makers to strengthen routine MCV1 immunization programmes and provide equitable disease protection for all children