Mapping development and health effects of cooking with solid fuels in low-income and middle-income countries, 2000–18: a geospatial modelling study

Background More than 3 billion people do not have access to clean energy and primarily use solid fuels to cook. Use of solid fuels generates household air pollution, which was associated with more than 2 million deaths in 2019. Although local patterns in cooking vary systematically, subnational trends in use of solid fuels have yet to be comprehensively analysed. We estimated the prevalence of solid-fuel use with high spatial resolution to explore subnational inequalities, assess local progress, and assess the effects on health in low-income and middle-income countries (LMICs) without universal access to clean fuels. Methods We did a geospatial modelling study to map the prevalence of solid-fuel use for cooking at a 5 km × 5 km resolution in 98 LMICs based on 2·1 million household observations of the primary cooking fuel used from 663 population-based household surveys over the years 2000 to 2018. We use observed temporal patterns to forecast household air pollution in 2030 and to assess the probability of attaining the Sustainable Development Goal (SDG) target indicator for clean cooking. We aligned our estimates of household air pollution to geospatial estimates of ambient air pollution to establish the risk transition occurring in LMICs. Finally, we quantified the effect of residual primary solid-fuel use for cooking on child health by doing a counterfactual risk assessment to estimate the proportion of deaths from lower respiratory tract infections in children younger than 5 years that could be associated with household air pollution. Findings Although primary reliance on solid-fuel use for cooking has declined globally, it remains widespread. 593 million people live in districts where the prevalence of solid-fuel use for cooking exceeds 95%. 66% of people in LMICs live in districts that are not on track to meet the SDG target for universal access to clean energy by 2030. Household air pollution continues to be a major contributor to particulate exposure in LMICs, and rising ambient air pollution is undermining potential gains from reductions in the prevalence of solid-fuel use for cooking in many countries. We estimated that, in 2018, 205 000 (95% uncertainty interval 147 000–257 000) children younger than 5 years died from lower respiratory tract infections that could be attributed to household air pollution. Interpretation Efforts to accelerate the adoption of clean cooking fuels need to be substantially increased and recalibrated to account for subnational inequalities, because there are substantial opportunities to improve air quality and avert child mortality associated with household air pollution.This study was funded by the Bill & Melinda Gates Foundation. L G Abreu acknowledges support from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (Capes; finance Code 001), Conselho Nacional de Desenvolvimento Científico e Tecnológico, and Fundação de Amparo à Pesquisa do Estado de Minas Gerais. D A Bennett acknowledges support from the Oxford National Institute for Health Research (NIHR) Biomedical Research Centre (BRC). The views expressed are those of the author and not necessarily those of the NHS, the NIHR, or the UK Department of Health and Social Care. Z A Bhutta acknowledges support from the Institute for Global Health & Development at the Aga Khan University. F Carvalho acknowledges UID/MULTI/04378/2019 and UID/QUI/50006/2019 support with funding from FCT/MCTES through national funds. J-W De Neve is supported by the Alexander von Humboldt Foundation. S Dey acknowledges the support from the Centre of Excellence for Research on Clean Air, IIT Delhi. M Ausloos and C Herteliu are partly supported by a grant of the Romanian National Authority for Scientific Research and Innovation (project number PN-III-P4-ID-PCCF-2016-0084). C Herteliu is partly supported by a grant of the Romanian National Authority for Scientific Research and Innovation (project number PN-III-P2-2.1-SOL-2020-2-0351), the Romanian Ministry of Research Innovation and Digitalization (project number ID-585-CTR-42-PFE-2021), and the Romanian Ministry of Labour and Social Justice (30/PSCD/2018). M Jakovljevic acknowledges partial support through Grant OI 175 014 of the Ministry of Science Education and Technological Development of the Republic of Serbia. J S John acknowledges support from the Kunshan Government and China Center for Disease Control and Prevention. W Mendoza is a program analyst in population and development at the United Nations Population Fund country office in Peru, an institution that does not necessarily endorse this study. M N Khan acknowledges the support of Jatiya Kabi Kazi Nazrul Islam University, Bangladesh. K Krishan is supported by UGC Centre of Advanced Study (CAS II), awarded to the Department of Anthropology, Panjab University, Chandigarh, India. M Kumar acknowledges support (FIC/NIH funded K43 TW010716-04 study). B Lacey acknowledges support from UK Biobank, the NIHR Oxford Biomedical Research Centre, and the British Heart Foundation Oxford Centre of Research Excellence. B R Nascimento acknowledges support in part by CNPq (Bolsa de produtividade em pesquisa, 312382/2019-7), by the Edwards Lifesciences Foundation (Every Heartbeat Matters Program 2020) and by FAPEMIG (grant APQ-000627-20). A M Samy acknowledges the support from the Egyptian Fulbright Mission Program. M M Santric-Milicevic acknowledges the support of the Ministry of Education, Science and Technological Development of Serbia (contract 175087). A Sheikh acknowledges the support of Health Data Research UK. I N Soyiri acknowledges support from the University of Hull internal QR Global Challenges Research Fund. S B Zaman acknowledges receiving a scholarship from the Australian Government research training program in support of his academic career. Y Zhang was supported by Science and Technology Research Project of Hubei Provincial Department of Education (grant Q20201104) and Outstanding Young and Middle Aged Technology Innovation Team Project of Hubei Provincial Department of Education (grant T2020003).publishedVersio

Global, regional, and national sex differences in the global burden of tuberculosis by HIV status, 1990–2019: results from the Global Burden of Disease Study 2019

Tuberculosis is a major contributor to the global burden of disease, causing more than a million deaths annually. Given an emphasis on equity in access to diagnosis and treatment of tuberculosis in global health targets, evaluations of differences in tuberculosis burden by sex are crucial. We aimed to assess the levels and trends of the global burden of tuberculosis, with an emphasis on investigating differences in sex by HIV status for 204 countries and territories from 1990 to 2019.publishedVersio

Assessing performance of the Healthcare Access and Quality Index, overall and by select age groups, for 204 countries and territories, 1990–2019: a systematic analysis from the Global Burden of Disease Study 2019

Background Health-care needs change throughout the life course. It is thus crucial to assess whether health systems provide access to quality health care for all ages. Drawing from the Global Burden of Diseases, Injuries, and Risk Factors Study 2019 (GBD 2019), we measured the Healthcare Access and Quality (HAQ) Index overall and for select age groups in 204 locations from 1990 to 2019. Methods We distinguished the overall HAQ Index (ages 0–74 years) from scores for select age groups: the young (ages 0–14 years), working (ages 15–64 years), and post-working (ages 65–74 years) groups. For GBD 2019, HAQ Index construction methods were updated to use the arithmetic mean of scaled mortality-to-incidence ratios (MIRs) and risk-standardised death rates (RSDRs) for 32 causes of death that should not occur in the presence of timely, quality health care. Across locations and years, MIRs and RSDRs were scaled from 0 (worst) to 100 (best) separately, putting the HAQ Index on a different relative scale for each age group. We estimated absolute convergence for each group on the basis of whether the HAQ Index grew faster in absolute terms between 1990 and 2019 in countries with lower 1990 HAQ Index scores than countries with higher 1990 HAQ Index scores and by Socio-demographic Index (SDI) quintile. SDI is a summary metric of overall development. Findings Between 1990 and 2019, the HAQ Index increased overall (by 19·6 points, 95% uncertainty interval 17·9–21·3), as well as among the young (22·5, 19·9–24·7), working (17·2, 15·2–19·1), and post-working (15·1, 13·2–17·0) age groups. Large differences in HAQ Index scores were present across SDI levels in 2019, with the overall index ranging from 30·7 (28·6–33·0) on average in low-SDI countries to 83·4 (82·4–84·3) on average in highSDI countries. Similarly large ranges between low-SDI and high-SDI countries, respectively, were estimated in the HAQ Index for the young (40·4–89·0), working (33·8–82·8), and post-working (30·4–79·1) groups. Absolute convergence in HAQ Index was estimated in the young group only. In contrast, divergence was estimated among the working and post-working groups, driven by slow progress in low-SDI countries. Interpretation Although major gaps remain across levels of social and economic development, convergence in the young group is an encouraging sign of reduced disparities in health-care access and quality. However, divergence in the working and post-working groups indicates that health-care access and quality is lagging at lower levels of social and economic development. To meet the needs of ageing populations, health systems need to improve health-care access and quality for working-age adults and older populations while continuing to realise gains among the young.publishedVersio

Public health utility of cause of death data : applying empirical algorithms to improve data quality

Background: Accurate, comprehensive, cause-specific mortality estimates are crucial for informing public health decision making worldwide. Incorrectly or vaguely assigned deaths, defined as garbage-coded deaths, mask the true cause distribution. The Global Burden of Disease (GBD) study has developed methods to create comparable, timely, cause-specific mortality estimates; an impactful data processing method is the reallocation of garbage-coded deaths to a plausible underlying cause of death. We identify the pattern of garbage-coded deaths in the world and present the methods used to determine their redistribution to generate more plausible cause of death assignments. Methods: We describe the methods developed for the GBD 2019 study and subsequent iterations to redistribute garbage-coded deaths in vital registration data to plausible underlying causes. These methods include analysis of multiple cause data, negative correlation, impairment, and proportional redistribution. We classify garbage codes into classes according to the level of specificity of the reported cause of death (CoD) and capture trends in the global pattern of proportion of garbage-coded deaths, disaggregated by these classes, and the relationship between this proportion and the Socio-Demographic Index. We examine the relative importance of the top four garbage codes by age and sex and demonstrate the impact of redistribution on the annual GBD CoD rankings. Results: The proportion of least-specific (class 1 and 2) garbage-coded deaths ranged from 3.7% of all vital registration deaths to 67.3% in 2015, and the age-standardized proportion had an overall negative association with the Socio Demographic Index. When broken down by age and sex, the category for unspecified lower respiratory infections was responsible for nearly 30% of garbage-coded deaths in those under 1 year of age for both sexes, representing the largest proportion of garbage codes for that age group. We show how the cause distribution by number of deaths changes before and after redistribution for four countries: Brazil, the United States, Japan, and France, highlighting the necessity of accounting for garbage-coded deaths in the GBD

Oral rehydration therapies in Senegal, Mali, and Sierra Leone: a spatial analysis of changes over time and implications for policy.

BACKGROUND: Oral rehydration solution (ORS) is a simple intervention that can prevent childhood deaths from severe diarrhea and dehydration. In a previous study, we mapped the use of ORS treatment subnationally and found that ORS coverage increased over time, while the use of home-made alternatives or recommended home fluids (RHF) decreased, in many countries. These patterns were particularly striking within Senegal, Mali, and Sierra Leone. It was unclear, however, whether ORS replaced RHF in these locations or if children were left untreated, and if these patterns were associated with health policy changes. METHODS: We used a Bayesian geostatistical model and data from household surveys to map the percentage of children with diarrhea that received (1) any ORS, (2) only RHF, or (3) no oral rehydration treatment between 2000 and 2018. This approach allowed examination of whether RHF was replaced with ORS before and after interventions, policies, and external events that may have impacted healthcare access. RESULTS: We found that RHF was replaced with ORS in most Sierra Leone districts, except those most impacted by the Ebola outbreak. In addition, RHF was replaced in northern but not in southern Mali, and RHF was not replaced anywhere in Senegal. In Senegal, there was no statistical evidence that a national policy promoting ORS use was associated with increases in coverage. In Sierra Leone, ORS coverage increased following a national policy change that abolished health costs for children. CONCLUSIONS: Children in parts of Mali and Senegal have been left behind during ORS scale-up. Improved messaging on effective diarrhea treatment and/or increased ORS access such as through reducing treatment costs may be needed to prevent child deaths in these areas

Measuring the availability of human resources for health and its relationship to universal health coverage for 204 countries and territories from 1990 to 2019: a systematic analysis for the Global Burden of Disease Study 2019

Background: Human resources for health (HRH) include a range of occupations that aim to promote or improve human health. The UN Sustainable Development Goals (SDGs) and the WHO Health Workforce 2030 strategy have drawn attention to the importance of HRH for achieving policy priorities such as universal health coverage (UHC). Although previous research has found substantial global disparities in HRH, the absence of comparable cross-national estimates of existing workforces has hindered efforts to quantify workforce requirements to meet health system goals. We aimed to use comparable and standardised data sources to estimate HRH densities globally, and to examine the relationship between a subset of HRH cadres and UHC effective coverage performance. Methods: Through the International Labour Organization and Global Health Data Exchange databases, we identified 1404 country-years of data from labour force surveys and 69 country-years of census data, with detailed microdata on health-related employment. From the WHO National Health Workforce Accounts, we identified 2950 country-years of data. We mapped data from all occupational coding systems to the International Standard Classification of Occupations 1988 (ISCO-88), allowing for standardised estimation of densities for 16 categories of health workers across the full time series. Using data from 1990 to 2019 for 196 of 204 countries and territories, covering seven Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) super-regions and 21 regions, we applied spatiotemporal Gaussian process regression (ST-GPR) to model HRH densities from 1990 to 2019 for all countries and territories. We used stochastic frontier meta-regression to model the relationship between the UHC effective coverage index and densities for the four categories of health workers enumerated in SDG indicator 3.c.1 pertaining to HRH: physicians, nurses and midwives, dentistry personnel, and pharmaceutical personnel. We identified minimum workforce density thresholds required to meet a specified target of 80 out of 100 on the UHC effective coverage index, and quantified national shortages with respect to those minimum thresholds. Findings: We estimated that, in 2019, the world had 104·0 million (95% uncertainty interval 83·5–128·0) health workers, including 12·8 million (9·7–16·6) physicians, 29·8 million (23·3–37·7) nurses and midwives, 4·6 million (3·6–6·0) dentistry personnel, and 5·2 million (4·0–6·7) pharmaceutical personnel. We calculated a global physician density of 16·7 (12·6–21·6) per 10 000 population, and a nurse and midwife density of 38·6 (30·1–48·8) per 10 000 population. We found the GBD super-regions of sub-Saharan Africa, south Asia, and north Africa and the Middle East had the lowest HRH densities. To reach 80 out of 100 on the UHC effective coverage index, we estimated that, per 10 000 population, at least 20·7 physicians, 70·6 nurses and midwives, 8·2 dentistry personnel, and 9·4 pharmaceutical personnel would be needed. In total, the 2019 national health workforces fell short of these minimum thresholds by 6·4 million physicians, 30·6 million nurses and midwives, 3·3 million dentistry personnel, and 2·9 million pharmaceutical personnel. Interpretation: Considerable expansion of the world's health workforce is needed to achieve high levels of UHC effective coverage. The largest shortages are in low-income settings, highlighting the need for increased financing and coordination to train, employ, and retain human resources in the health sector. Actual HRH shortages might be larger than estimated because minimum thresholds for each cadre of health workers are benchmarked on health systems that most efficiently translate human resources into UHC attainment

Subnational mapping of HIV incidence and mortality among individuals aged 15–49 years in sub-Saharan Africa, 2000–18 : a modelling study

Background: High-resolution estimates of HIV burden across space and time provide an important tool for tracking and monitoring the progress of prevention and control efforts and assist with improving the precision and efficiency of targeting efforts. We aimed to assess HIV incidence and HIV mortality for all second-level administrative units across sub-Saharan Africa. Methods: In this modelling study, we developed a framework that used the geographically specific HIV prevalence data collected in seroprevalence surveys and antenatal care clinics to train a model that estimates HIV incidence and mortality among individuals aged 15–49 years. We used a model-based geostatistical framework to estimate HIV prevalence at the second administrative level in 44 countries in sub-Saharan Africa for 2000–18 and sought data on the number of individuals on antiretroviral therapy (ART) by second-level administrative unit. We then modified the Estimation and Projection Package (EPP) to use these HIV prevalence and treatment estimates to estimate HIV incidence and mortality by second-level administrative unit. Findings: The estimates suggest substantial variation in HIV incidence and mortality rates both between and within countries in sub-Saharan Africa, with 15 countries having a ten-times or greater difference in estimated HIV incidence between the second-level administrative units with the lowest and highest estimated incidence levels. Across all 44 countries in 2018, HIV incidence ranged from 2 ·8 (95% uncertainty interval 2·1–3·8) in Mauritania to 1585·9 (1369·4–1824·8) cases per 100 000 people in Lesotho and HIV mortality ranged from 0·8 (0·7–0·9) in Mauritania to 676· 5 (513· 6–888·0) deaths per 100 000 people in Lesotho. Variation in both incidence and mortality was substantially greater at the subnational level than at the national level and the highest estimated rates were accordingly higher. Among second-level administrative units, Guijá District, Gaza Province, Mozambique, had the highest estimated HIV incidence (4661·7 [2544·8–8120·3]) cases per 100000 people in 2018 and Inhassunge District, Zambezia Province, Mozambique, had the highest estimated HIV mortality rate (1163·0 [679·0–1866·8]) deaths per 100 000 people. Further, the rate of reduction in HIV incidence and mortality from 2000 to 2018, as well as the ratio of new infections to the number of people living with HIV was highly variable. Although most second-level administrative units had declines in the number of new cases (3316 [81· 1%] of 4087 units) and number of deaths (3325 [81·4%]), nearly all appeared well short of the targeted 75% reduction in new cases and deaths between 2010 and 2020. Interpretation: Our estimates suggest that most second-level administrative units in sub-Saharan Africa are falling short of the targeted 75% reduction in new cases and deaths by 2020, which is further compounded by substantial within-country variability. These estimates will help decision makers and programme implementers expand access to ART and better target health resources to higher burden subnational areas

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

Mapping geographical inequalities in access to drinking water and sanitation facilities in low-income and middle-income countries, 2000-17

Background: Universal access to safe drinking water and sanitation facilities is an essential human right, recognised in the Sustainable Development Goals as crucial for preventing disease and improving human wellbeing. Comprehensive, high-resolution estimates are important to inform progress towards achieving this goal. We aimed to produce high-resolution geospatial estimates of access to drinking water and sanitation facilities. Methods: We used a Bayesian geostatistical model and data from 600 sources across more than 88 low-income and middle-income countries (LMICs) to estimate access to drinking water and sanitation facilities on continuous continent-wide surfaces from 2000 to 2017, and aggregated results to policy-relevant administrative units. We estimated mutually exclusive and collectively exhaustive subcategories of facilities for drinking water (piped water on or off premises, other improved facilities, unimproved, and surface water) and sanitation facilities (septic or sewer sanitation, other improved, unimproved, and open defecation) with use of ordinal regression. We also estimated the number of diarrhoeal deaths in children younger than 5 years attributed to unsafe facilities and estimated deaths that were averted by increased access to safe facilities in 2017, and analysed geographical inequality in access within LMICs. Findings: Across LMICs, access to both piped water and improved water overall increased between 2000 and 2017, with progress varying spatially. For piped water, the safest water facility type, access increased from 40·0% (95% uncertainty interval [UI] 39·4–40·7) to 50·3% (50·0–50·5), but was lowest in sub-Saharan Africa, where access to piped water was mostly concentrated in urban centres. Access to both sewer or septic sanitation and improved sanitation overall also increased across all LMICs during the study period. For sewer or septic sanitation, access was 46·3% (95% UI 46·1–46·5) in 2017, compared with 28·7% (28·5–29·0) in 2000. Although some units improved access to the safest drinking water or sanitation facilities since 2000, a large absolute number of people continued to not have access in several units with high access to such facilities (>80%) in 2017. More than 253 000 people did not have access to sewer or septic sanitation facilities in the city of Harare, Zimbabwe, despite 88·6% (95% UI 87·2–89·7) access overall. Many units were able to transition from the least safe facilities in 2000 to safe facilities by 2017; for units in which populations primarily practised open defecation in 2000, 686 (95% UI 664–711) of the 1830 (1797–1863) units transitioned to the use of improved sanitation. Geographical disparities in access to improved water across units decreased in 76·1% (95% UI 71·6–80·7) of countries from 2000 to 2017, and in 53·9% (50·6–59·6) of countries for access to improved sanitation, but remained evident subnationally in most countries in 2017. Interpretation: Our estimates, combined with geospatial trends in diarrhoeal burden, identify where efforts to increase access to safe drinking water and sanitation facilities are most needed. By highlighting areas with successful approaches or in need of targeted interventions, our estimates can enable precision public health to effectively progress towards universal access to safe water and sanitation

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.Peer reviewe