Global fertility in 204 countries and territories, 1950–2021, with forecasts to 2100: a comprehensive demographic analysis for the Global Burden of Disease Study 2021

BACKGROUND: Accurate assessments of current and future fertility—including overall trends and changing population age structures across countries and regions—are essential to help plan for the profound social, economic, environmental, and geopolitical challenges that these changes will bring. Estimates and projections of fertility are necessary to inform policies involving resource and health-care needs, labour supply, education, gender equality, and family planning and support. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 produced up-to-date and comprehensive demographic assessments of key fertility indicators at global, regional, and national levels from 1950 to 2021 and forecast fertility metrics to 2100 based on a reference scenario and key policy-dependent alternative scenarios. METHODS: To estimate fertility indicators from 1950 to 2021, mixed-effects regression models and spatiotemporal Gaussian process regression were used to synthesise data from 8709 country-years of vital and sample registrations, 1455 surveys and censuses, and 150 other sources, and to generate age-specific fertility rates (ASFRs) for 5-year age groups from age 10 years to 54 years. ASFRs were summed across age groups to produce estimates of total fertility rate (TFR). Livebirths were calculated by multiplying ASFR and age-specific female population, then summing across ages 10–54 years. To forecast future fertility up to 2100, our Institute for Health Metrics and Evaluation (IHME) forecasting model was based on projections of completed cohort fertility at age 50 years (CCF50; the average number of children born over time to females from a specified birth cohort), which yields more stable and accurate measures of fertility than directly modelling TFR. CCF50 was modelled using an ensemble approach in which three sub-models (with two, three, and four covariates variously consisting of female educational attainment, contraceptive met need, population density in habitable areas, and under-5 mortality) were given equal weights, and analyses were conducted utilising the MR-BRT (meta-regression—Bayesian, regularised, trimmed) tool. To capture time-series trends in CCF50 not explained by these covariates, we used a first-order autoregressive model on the residual term. CCF50 as a proportion of each 5-year ASFR was predicted using a linear mixed-effects model with fixed-effects covariates (female educational attainment and contraceptive met need) and random intercepts for geographical regions. Projected TFRs were then computed for each calendar year as the sum of single-year ASFRs across age groups. The reference forecast is our estimate of the most likely fertility future given the model, past fertility, forecasts of covariates, and historical relationships between covariates and fertility. We additionally produced forecasts for multiple alternative scenarios in each location: the UN Sustainable Development Goal (SDG) for education is achieved by 2030; the contraceptive met need SDG is achieved by 2030; pro-natal policies are enacted to create supportive environments for those who give birth; and the previous three scenarios combined. Uncertainty from past data inputs and model estimation was propagated throughout analyses by taking 1000 draws for past and present fertility estimates and 500 draws for future forecasts from the estimated distribution for each metric, with 95% uncertainty intervals (UIs) given as the 2·5 and 97·5 percentiles of the draws. To evaluate the forecasting performance of our model and others, we computed skill values—a metric assessing gain in forecasting accuracy—by comparing predicted versus observed ASFRs from the past 15 years (2007–21). A positive skill metric indicates that the model being evaluated performs better than the baseline model (here, a simplified model holding 2007 values constant in the future), and a negative metric indicates that the evaluated model performs worse than baseline. FINDINGS: During the period from 1950 to 2021, global TFR more than halved, from 4·84 (95% UI 4·63–5·06) to 2·23 (2·09–2·38). Global annual livebirths peaked in 2016 at 142 million (95% UI 137–147), declining to 129 million (121–138) in 2021. Fertility rates declined in all countries and territories since 1950, with TFR remaining above 2·1—canonically considered replacement-level fertility—in 94 (46·1%) countries and territories in 2021. This included 44 of 46 countries in sub-Saharan Africa, which was the super-region with the largest share of livebirths in 2021 (29·2% [28·7–29·6]). 47 countries and territories in which lowest estimated fertility between 1950 and 2021 was below replacement experienced one or more subsequent years with higher fertility; only three of these locations rebounded above replacement levels. Future fertility rates were projected to continue to decline worldwide, reaching a global TFR of 1·83 (1·59–2·08) in 2050 and 1·59 (1·25–1·96) in 2100 under the reference scenario. The number of countries and territories with fertility rates remaining above replacement was forecast to be 49 (24·0%) in 2050 and only six (2·9%) in 2100, with three of these six countries included in the 2021 World Bank-defined low-income group, all located in the GBD super-region of sub-Saharan Africa. The proportion of livebirths occurring in sub-Saharan Africa was forecast to increase to more than half of the world's livebirths in 2100, to 41·3% (39·6–43·1) in 2050 and 54·3% (47·1–59·5) in 2100. The share of livebirths was projected to decline between 2021 and 2100 in most of the six other super-regions—decreasing, for example, in south Asia from 24·8% (23·7–25·8) in 2021 to 16·7% (14·3–19·1) in 2050 and 7·1% (4·4–10·1) in 2100—but was forecast to increase modestly in the north Africa and Middle East and high-income super-regions. Forecast estimates for the alternative combined scenario suggest that meeting SDG targets for education and contraceptive met need, as well as implementing pro-natal policies, would result in global TFRs of 1·65 (1·40–1·92) in 2050 and 1·62 (1·35–1·95) in 2100. The forecasting skill metric values for the IHME model were positive across all age groups, indicating that the model is better than the constant prediction. INTERPRETATION: Fertility is declining globally, with rates in more than half of all countries and territories in 2021 below replacement level. Trends since 2000 show considerable heterogeneity in the steepness of declines, and only a small number of countries experienced even a slight fertility rebound after their lowest observed rate, with none reaching replacement level. Additionally, the distribution of livebirths across the globe is shifting, with a greater proportion occurring in the lowest-income countries. Future fertility rates will continue to decline worldwide and will remain low even under successful implementation of pro-natal policies. These changes will have far-reaching economic and societal consequences due to ageing populations and declining workforces in higher-income countries, combined with an increasing share of livebirths among the already poorest regions of the world. FUNDING: Bill & Melinda Gates Foundation

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

Global, regional, and national mortality among young people aged 10–24 years, 1950–2019: a systematic analysis for the Global Burden of Disease Study 2019

Summary: Background Documentation of patterns and long-term trends in mortality in young people, which reflect huge changes in demographic and social determinants of adolescent health, enables identification of global investment priorities for this age group. We aimed to analyse data on the number of deaths, years of life lost, and mortality rates by sex and age group in people aged 10–24 years in 204 countries and territories from 1950 to 2019 by use of estimates from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019. Methods We report trends in estimated total numbers of deaths and mortality rate per 100 000 population in young people aged 10–24 years by age group (10–14 years, 15–19 years, and 20–24 years) and sex in 204 countries and territories between 1950 and 2019 for all causes, and between 1980 and 2019 by cause of death. We analyse variation in outcomes by region, age group, and sex, and compare annual rate of change in mortality in young people aged 10–24 years with that in children aged 0–9 years from 1990 to 2019. We then analyse the association between mortality in people aged 10–24 years and socioeconomic development using the GBD Socio-demographic Index (SDI), a composite measure based on average national educational attainment in people older than 15 years, total fertility rate in people younger than 25 years, and income per capita. We assess the association between SDI and all-cause mortality in 2019, and analyse the ratio of observed to expected mortality by SDI using the most recent available data release (2017). Findings In 2019 there were 1·49 million deaths (95% uncertainty interval 1·39–1·59) worldwide in people aged 10–24 years, of which 61% occurred in males. 32·7% of all adolescent deaths were due to transport injuries, unintentional injuries, or interpersonal violence and conflict; 32·1% were due to communicable, nutritional, or maternal causes; 27·0% were due to non-communicable diseases; and 8·2% were due to self-harm. Since 1950, deaths in this age group decreased by 30·0% in females and 15·3% in males, and sex-based differences in mortality rate have widened in most regions of the world. Geographical variation has also increased, particularly in people aged 10–14 years. Since 1980, communicable and maternal causes of death have decreased sharply as a proportion of total deaths in most GBD super-regions, but remain some of the most common causes in sub-Saharan Africa and south Asia, where more than half of all adolescent deaths occur. Annual percentage decrease in all-cause mortality rate since 1990 in adolescents aged 15–19 years was 1·3% in males and 1·6% in females, almost half that of males aged 1–4 years (2·4%), and around a third less than in females aged 1–4 years (2·5%). The proportion of global deaths in people aged 0–24 years that occurred in people aged 10–24 years more than doubled between 1950 and 2019, from 9·5% to 21·6%. Interpretation Variation in adolescent mortality between countries and by sex is widening, driven by poor progress in reducing deaths in males and older adolescents. Improving global adolescent mortality will require action to address the specific vulnerabilities of this age group, which are being overlooked. Furthermore, indirect effects of the COVID-19 pandemic are likely to jeopardise efforts to improve health outcomes including mortality in young people aged 10–24 years. There is an urgent need to respond to the changing global burden of adolescent mortality, address inequities where they occur, and improve the availability and quality of primary mortality data in this age group

Mapping inequalities in exclusive breastfeeding in low- and middle-income countries, 2000–2018

Exclusive breastfeeding (EBF)-giving infants only breast-milk for the first 6 months of life-is a component of optimal breastfeeding practices effective in preventing child morbidity and mortality. EBF practices are known to vary by population and comparable subnational estimates of prevalence and progress across low- and middle-income countries (LMICs) are required for planning policy and interventions. Here we present a geospatial analysis of EBF prevalence estimates from 2000 to 2018 across 94 LMICs mapped to policy-relevant administrative units (for example, districts), quantify subnational inequalities and their changes over time, and estimate probabilities of meeting the World Health Organization's Global Nutrition Target (WHO GNT) of ≥70% EBF prevalence by 2030. While six LMICs are projected to meet the WHO GNT of ≥70% EBF prevalence at a national scale, only three are predicted to meet the target in all their district-level units by 2030.This work was primarily supported by grant no. OPP1132415 from the Bill & Melinda Gates Foundation. Co-authors used by the Bill & Melinda Gates Foundation (E.G.P. and R.R.3) provided feedback on initial maps and drafts of this manuscript. L.G.A. has received support from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brasil (CAPES), Código de Financiamento 001 and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (grant nos. 404710/2018-2 and 310797/2019-5). O.O.Adetokunboh acknowledges the National Research Foundation, Department of Science and Innovation and South African Centre for Epidemiological Modelling and Analysis. M.Ausloos, A.Pana and C.H. are partially supported by a grant from the Romanian National Authority for Scientific Research and Innovation, CNDS-UEFISCDI, project no. PN-III-P4-ID-PCCF-2016-0084. P.C.B. would like to acknowledge the support of F. Alam and A. Hussain. T.W.B. was supported by the Alexander von Humboldt Foundation through the Alexander von Humboldt Professor award, funded by the German Federal Ministry of Education and Research. K.Deribe is supported by the Wellcome Trust (grant no. 201900/Z/16/Z) as part of his international intermediate fellowship. C.H. and A.Pana are partially supported by a grant of the Romanian National Authority for Scientific Research and Innovation, CNDS-UEFISCDI, project no. PN-III-P2-2.1-SOL-2020-2-0351. B.Hwang is partially supported by China Medical University (CMU109-MF-63), Taichung, Taiwan. M.Khan acknowledges Jatiya Kabi Kazi Nazrul Islam University for their support. A.M.K. acknowledges the other collaborators and the corresponding author. Y.K. was supported by the Research Management Centre, Xiamen University Malaysia (grant no. XMUMRF/2020-C6/ITM/0004). K.Krishan is supported by a DST PURSE grant and UGC Centre of Advanced Study (CAS II) awarded to the Department of Anthropology, Panjab University, Chandigarh, India. M.Kumar would like to acknowledge FIC/NIH K43 TW010716-03. I.L. is a member of the Sistema Nacional de Investigación (SNI), which is supported by the Secretaría Nacional de Ciencia, Tecnología e Innovación (SENACYT), Panamá. M.L. was supported by China Medical University, Taiwan (CMU109-N-22 and CMU109-MF-118). W.M. is currently a programme analyst in Population and Development at the United Nations Population Fund (UNFPA) Country Office in Peru, which does not necessarily endorses this study. D.E.N. acknowledges Cochrane South Africa, South African Medical Research Council. G.C.P. is supported by an NHMRC research fellowship. P.Rathi acknowledges support from Kasturba Medical College, Mangalore, Manipal Academy of Higher Education, Manipal, India. Ramu Rawat acknowledges the support of the GBD Secretariat for supporting the reviewing and collaboration of this paper. B.R. acknowledges support from Manipal College of Health Professions, Manipal Academy of Higher Education, Manipal. A.Ribeiro was supported by National Funds through FCT, under the programme of ‘Stimulus of Scientific Employment—Individual Support’ within the contract no. info:eu-repo/grantAgreement/FCT/CEEC IND 2018/CEECIND/02386/2018/CP1538/CT0001/PT. S.Sajadi acknowledges colleagues at Global Burden of Diseases and Local Burden of Disease. A.M.S. acknowledges the support from the Egyptian Fulbright Mission Program. F.S. was supported by the Shenzhen Science and Technology Program (grant no. KQTD20190929172835662). A.Sheikh is supported by Health Data Research UK. B.K.S. acknowledges Kasturba Medical College, Mangalore, Manipal Academy of Higher Education, Manipal for all the academic support. B.U. acknowledges support from Manipal Academy of Higher Education, Manipal. C.S.W. is supported by the South African Medical Research Council. Y.Z. was supported by Science and Technology Research Project of Hubei Provincial Department of Education (grant no. Q20201104) and Outstanding Young and Middle-aged Technology Innovation Team Project of Hubei Provincial Department of Education (grant no. T2020003). The funders of the study had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication. All maps presented in this study are generated by the authors and no permissions are required to publish them