Article Text

Projecting sex imbalances at birth at global, regional and national levels from 2021 to 2100: scenario-based Bayesian probabilistic projections of the sex ratio at birth and missing female births based on 3.26 billion birth records
  1. Fengqing Chao1,
  2. Patrick Gerland2,
  3. Alex Richard Cook3,
  4. Christophe Z Guilmoto4,
  5. Leontine Alkema5
  1. 1 Statistics Program, Computer Electrical and Mathematical Science and Engineering Division, King Abdullah University of Science and Technology, Thuwal, Makkah, Saudi Arabia
  2. 2 Population Division, United Nations, New York, New York, USA
  3. 3 Saw Swee Hock School of Public Health, National University of Singapore, Singapore
  4. 4 CEPED/IRD, Centre de Sciences Humaines, New Delhi, India
  5. 5 School of Public Health and Health Sciences, University of Massachusetts Amherst, Amherst, Massachusetts, USA
  1. Correspondence to Dr Fengqing Chao; fengqing.chao{at}kaust.edu.sa

Abstract

Introduction Skewed levels of the sex ratio at birth (SRB) due to sex-selective abortions have been observed in several countries since the 1970s. They will lead to long-term sex imbalances in more than one-third of the world’s population with yet unknown social and economic impacts on affected countries. Understanding the potential evolution of sex imbalances at birth is therefore essential for anticipating and planning for changing sex structures across the world.

Methods We produced probabilistic SRB projections from 2021 to 2100 based on different scenarios of sex ratio transition and assessed their implications in terms of missing female births at global, regional and national levels. Based on a comprehensive SRB database with 3.26 billion birth records, we project the skewed SRB and missing female births with a Bayesian hierarchical time series mixture model. The SRB projections under reference scenario S1 assumed SRB transitions only for countries with strong statistical evidence of SRB inflation, and the more extreme scenario S2 assumed a sex ratio transition for countries at risk of SRB inflation but with no or limited evidence of ongoing inflation.

Results Under scenario S1, we projected 5.7 (95% uncertainty interval (1.2; 15.3)) million additional missing female births to occur by 2100. Countries affected will be those already affected in the past by imbalanced SRB, such as China and India. If all countries at risk of SRB inflation experience a sex ratio transition as in scenario S2, the projected missing female births increase to 22.1 (12.2; 39.8) million with a sizeable contribution of sub-Saharan Africa.

Conclusion The scenario-based projections provide important illustrations of the potential burden of future prenatal sex discrimination and the need to monitor SRBs in countries with son preference. Policy planning will be needed in the years to come to minimise future prenatal sex discrimination and its impact on social structures.

  • mathematical modelling
  • public health
  • medical demography

Data availability statement

Data are available in a public, open access repository. All data relevant to the study are included in the article or uploaded as online supplemental information. All data used in this paper are publicly available or included in the article.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Key questions

What is already known?

  • No prior study has constructed scenario-based projections for the sex ratio at birth and missing female births for all countries.

  • Approaches used for sex ratio at birth projections include the assumption of constant country-specific levels from 2017 to 2100 in the Global Burden of Disease study and projections based on expert opinion and linear extrapolation by the UN World Population Prospects and the US Census Bureau.

What are the new findings?

  • We provided projections of the sex ratio at birth and missing female births based on two scenarios for all countries till 2100.

  • Our projections were based on Bayesian hierarchical time series mixture models that were fitted to an extensive database on national sex ratio at birth. The hierarchical structure of the model allowed for sharing information on sex ratio transitions across countries and for projecting the potential future inflation process under different scenarios.

  • Our study provided insights into future changes in the sex ratio at birth as well as the uncertainty associated with these changes.

Key questions

What do the new findings imply?

  • The scenario-based projections of sex ratio at births showed that the sex ratio at birth is most likely to stabilise and decline within less than 20 years in countries currently affected by sex imbalances at birth.

  • However, there will be a deficit of more than 4.7 million female births between 2021 and 2030 under our conservative scenario.

  • A less optimistic scenario envisages a rise in the sex ratio in other countries with strong prevalence of son preference, including in populous countries such as Nigeria or Pakistan.

Introduction

Over the last 40 years, prenatal gender-biased sex selection has become the most visible consequence of son preference. Along with child marriage and female genital mutilation, sex selection is one of the key harmful practices defined by the United Nations (UN) and targeted under the Sustainable Development Goals (SDGs).1 Sex-selective abortions, the main mechanism behind sex selection, have been observed across a range of various countries from Southeast Europe to South and East Asia.2–11 They lead to a hike in the sex ratio at birth (SRB; namely the ratio of male to female live births) above its natural (biological) level and to the emergence of a surplus of male births. Along with excess female mortality, this SRB inflation is now a major contributor to the ‘missing women’, a concept first described by Sen to refer to a population with much fewer females than males.12 According to recent studies, prenatal sex selection accounts for about half of the recent deficit of females in the world during the previous decades.13 14 A male-biased sex structure in a society could lead to demographic issues such as marriage squeeze with lack of marriageable females.15 16 Fewer-than-expected females in a population could result in elevated levels of antisocial behaviour and violence, and may ultimately affect long-term stability and social sustainable development.17–20

Data limitations necessitate the use of models for estimating SRB levels and trends.13 Population statistics derived from birth registration, census figures and demographic surveys have contributed to estimates of the SRB in different countries and of its inflation due to prenatal sex selection.21 However, SRB series in many affected countries are spotty—for lack of reliable civil registration—and SRB rates must therefore often be re-estimated from a variety of sources. In prior work, we gathered a comprehensive database of SRB measurements and provided a new set of SRB estimates for all countries during 1950–2017 along with estimates of the natural SRB level in the absence of deliberate sex selection.13 The total number of female births missing over 1970–2017 due to prenatal sex selection was estimated at 45.0 million with a 95% uncertainty interval (34.4; 54.8) million. More than 95% of them were missing from China or India, countries combining severely skewed SRB levels and the largest numbers of annual births in the world.

The almost simultaneous rise of the SRB in various countries accounting for more than a third of the world’s population is an unprecedented phenomenon in demographic history. It will have long-lasting repercussions on demographic and social structures due to the mounting number of missing adult females, with long-term social consequences on affected societies.15 16 22 It is therefore crucial to be able to anticipate the future dynamics of SRB inflation across the world, in already affected countries and in those at risk. The main challenge is to understand whether birth masculinity will stay indefinitely skewed in countries affected by sex-selective abortions and whether new countries may be affected in the future.

No rigorous attempt has been made so far to outline and project the future evolution of the SRB imbalances in the world and to estimate its implications on the potential deficit of future female births. The most recent set of world population projections from the Global Burden of Disease Study assumes for instance birth masculinity to remain constant through 2100 at their 2017 levels, an unlikely feature in view of the rapid changes in SRB levels observed in countries with skewed SRBs.23 The population projections published every 2 years by the UN derives SRB trends for the next decades from expert opinion and linear extrapolation.24 The International Data Base developed by the US census Bureau similarly projects the SRB by linear extrapolation.25

The absence of evidence-based SRB projections stems mainly from the unique character of the human-made SRB imbalances observed over the last 40 years and the corresponding lack of historical references from which to infer its future evolution. A more technical factor relates to the apparently non-linear trends observed in several countries. Mortality and fertility tend to follow monotonic and near-linear patterns such as long-term declines that are easier to model and project. In contrast, the SRB trends observed so far appear to follow a more complex and non-monotonic cycle, characterised by a three-stage transformation: an increase, stabilisation and an ultimate decrease of the SRB, referred to as the sex ratio transition.26 The SRB transition process has been observed in several countries in the past decades.13 The curved bell-shaped pattern cannot be modelled by a simple linear regression model and requires additional hypotheses or covariates to recreate. However, the experience of countries already affected by SRB inflation can now help us delineate the expected patterns of SRB transformations in relation to its original drivers. We can use for this purpose the existing database of country-level SRB measurements from 1970 to 2017 and an extended modelling approach to develop a set of projection scenarios. This model will also make possible the estimation of the probability of a future SRB rise and subsequent decline in countries that have so far not been affected by SRB imbalances in spite of pronounced son preference.

Our paper provided probabilistic scenario-based projections of the SRB for 212 countries (‘countries’ or ‘areas’ as per the UN classification with population size greater than 90 000 in 2017) until 2100 based on the trends observed until 2020. We focused on 29 countries where son preference has been documented and for which we already have the SRB baseline levels used to assess the extent of the SRB inflation. They include 12 countries where SRB inflation has already happened in the past as well as 17 additional countries at risk of SRB inflation, that is, countries where the SRB may rise in the future. We produced Bayesian projections of the SRB for each country as well as the number of female births missing due to prenatal sex selection under different scenarios. The results are presented at global, regional, and national levels from 2021 to 2100 for two scenarios.

Data and methods

We summarised here the data and methods used. A more detailed description of the model specifications, statistical computation and model validation is available in a separate methodological paper.27 We focus in this paper on summarising model assumptions in less technical terms.

Data

We started from an extensive database of 3.26 billion birth records, corresponding to 12,341 SRB observations (ie, observed SRB for a certain country-year/period) from 204 countries, compiled based on information available in 2021. The database is publicly available.28 This database extends on a database used for prior work.29 The description of the construction of the original database and preprocessing steps are given elsewhere.13 The updated database was used for estimating the SRB from 1950 to 2020 and projecting the SRB till 2100. Compared with the previous database,29 we extended the database by adding in 1,506 observations in 166 countries from surveys, censuses and civil registration and vital statistics (CRVS). In particular, we added 653 CRVS data points since 2017 and 853 observations from censuses, Demographic and Health Surveys, and Multiple Indicator Cluster Surveys conducted since 2017. Online supplemental appendix provides detailed information on available data by country.

Supplemental material

We also used the latest UN demographic estimates and projections for fertility and birth trends from 1950 to 2100. The UN World Population Prospects (WPP) were released in 2019; we used the annual medium-variant projections.24 We used 1,000 total fertility rate (TFR) trajectories to account for uncertainty for each country-year.30 31

Sex ratio transition model

The findings rest on the Bayesian model used to project future SRB change based on past SRB measurements and fertility trends.27 The four guiding assumptions of the SRB transition model used here can be summarised as follows:

  1. Countries with past, ongoing or potential future SRB inflation can be identified using criteria related to son preference intensity and the timing of fertility decline.

  2. National baseline SRB estimates can be used as reference for the period under study.

  3. SRB inflation trajectories follow three successive stages: rise, plateau and decline.

  4. SRB trajectories in countries affected by SRB imbalances prior to 2021 can be used as a template for Bayesian projections of the future SRB levels.

The first assumption is based on the analysis of the three preconditions of the emergence of prenatal sex selection and the subsequent rise of the SRB: strong preference for male children; fertility decline constraining parents’ desired number of children; and access to modern sex-selection technology.26 The first precondition is used to identify the propensity to sex select based on various qualitative and quantitative indicators of son preference. It led to the identification of 29 countries with actual or potential prenatal sex selection.

The second assumption rests on the stability of natural (biological) SRB level and of the preexisting regional differentials.32–35 SRB series in industrialised countries with reliable registration system and no prenatal sex selection show that the SRB rarely change by more than 1% over a century. In addition, the presence of specific regional differentials—such as the lower natural SRB among populations of Sub-Saharan Africa or of African descent in the USA—is well established.21 36 Baseline estimates can therefore be used as a benchmark to model SRB inflation and missing female births. The SRB regional and national baselines are estimated based on a subset of the database: we exclude observations after 1970 from the 29 countries at risk of sex selection for the estimation of baseline levels such that baselines are not affected by data that may be subject to inflation.13 27 The SRB is modelled as the sum of baseline SRB level (with natural year-by-year fluctuations) and the product of an indicator to capture the presence or absence of inflation and a non-negative SRB inflation factor,

Embedded Image

The third assumption, that SRB trajectories successively rise, plateau and decline, is based on the non-linear patterns of SRB changes observed in countries affected by sex selection until today. The first stage is the initial increase of the SRB from its baseline level. This rise continues over several years and is followed by stabilisation at a higher SRB level, which constitutes the second stage of the cycle. The third stage begins after several years at a plateau level and corresponds to the final decline of the SRB back to its baseline natural level. Several countries (Hong Kong (SAR of China), Georgia and Republic of Korea) already illustrate these three completed phases of SRB transition. Other countries are currently going through the second stabilisation phase (Armenia, Azerbaijan, China, India and Vietnam) or the third declining phase (Albania, Montenegro, Taiwan (Province of China) and Tunisia). The parameterisation of the inflation factor follows therefore a trapezoid function to capture the three stages of the sex ratio transition: initial increase, subsequent plateau and final decrease. The model also captures regularities across countries in the patterns of the sex ratio transition and it makes use of the TFR as a covariate of the SRB inflation. It is used in particular for estimating the start year of the sex ratio transition process.

The fourth assumption, that SRB trajectories in countries affected by SRB imbalances prior to 2021 can be used as a template for Bayesian projections of the future SRB levels, is at the core of the Bayesian projection exercise. Conditioning on the third assumption, past observations are used to determine future trends in a probabilistic manner. The length of each transition stage and the maximum extent of the SRB inflation is estimated with uncertainty by using Bayesian hierarchical models. The hierarchical structure of the model allows for sharing of information about sex ratio transitions across countries. The implementation of the model together with its motivation and characteristics are fully described in the methodology paper.27

Scenario-based projections of SRB

Among the 29 countries at risk of sex selection, we identified 12 countries with a strong statistical evidence of SRB inflation when the posterior probability of presence of SRB inflation is more than 95%. For the 12 countries with strong evidence of existing SRB inflation, the sex ratio transition model is used to produce projections of future inflation.

For the remaining 17 countries at risk but without strong evidence of past/ongoing SRB inflation, we project SRB under scenarios with varying assumptions regarding the occurrence of a sex ratio transition (see online supplemental appendix for the explanation of scenario names and online supplemental appendix table 1 for the type of SRB projections for each country). The two scenarios are the following:

  • In the first scenario, S1, the 17 countries’ SRBs are projected without SRB inflation.

  • In the second scenario, S2, projections assume that the 17 countries will experience an SRB inflation with 100% probability.

S2 scenario accounts for the uncertainty of the TFR projections. The uncertainties in these projections become larger as we draw further from the current period. No other factors are used in projections under scenario S2.

Missing female births

To project the effect of SRB imbalances under different scenarios, we calculated the annual number of missing female births (AMFB) and the cumulative number of missing female births (CMFB) from 2021 to 2100. We define missing female births as female births prevented by sex selection, that is, to the number of sex-selective abortions, as defined in a recent study.37 The AMFB is computed as the difference between the numbers of female live births based on the baseline SRB and of those based on the projected SRB. The CMFB is the cumulated sum of the AMFB from 2021 till a given year. The number of live births is drawn from the 2019 UN WPP for 2021–2100.24

Patient and public involvement

Patients or the public were not involved in the design, or conduct, or reporting, or dissemination plans of our research.

Results

We produced scenario-based probabilistic projections of the SRB and number of missing female births on global, regional and national levels from 2021 to 2100. Using the reference year for the SDGs, results are presented for 2021–2030 and 2031–2100. Estimates for 1970–2020 and scenario-based projection results by country are in online supplemental appendix. Detailed annual estimates and scenario-based projections by country, region and world are also included in online supplemental file 3 and online supplemental file 4.

Supplemental material

Supplemental material

SRB projections for global and regional aggregates

Figure 1 illustrate the projected SRB for the world and by region for the two scenarios. At the global level, the SRB projections are similar. Under this scenario S1, in which sex imbalances are restricted to countries with strong statistical evidence of SRB inflation, the global SRB is projected to decline from 1.063 (1.054; 1.072) in 2020 to 1.052 (1.045; 1.061) by 2030 and to 1.045 (1.040; 1.050) by 2100 (online supplemental appendix table 3). The SRBs under this scenario are slightly lower than those under scenario S2 in which additional countries register a rise in their SRB. Under S2, the global SRB would decline to 1.055 (1.047; 1.064) by 2030 and 1.048 (1.041; 1.057) in 2100.

Figure 1

Global and regional SRB estimates and projections 1970–2100. Point estimates (solid lines) and 95% uncertainty intervals (shaded areas) for scenario-based projections. The green horizontal lines are regional baseline SRB median estimates. Selected regions contain at least one country at risk of sex imbalance. ENAN: the combination of countries in Europe, North America, Australia and New Zealand; SRB: sex ratio at birth.

The situation is more diverse across regions. The SRBs for S1 are projected to converge towards their aggregated regional baseline levels by 2100, ranging from 1.034 (1.026; 1.041) in sub-Saharan Africa to 1.074 (1.054; 1.095) in Oceania. Most of this convergence will take place around 2030, notably for the three largest contributors to world sex imbalances at birth (China, India and Vietnam). The decline is especially rapid in China, which displayed the highest SRB level in the world in 2010.

The projected SRBs for scenario S2 differ from those for S1 in regions where additional countries at risk of sex imbalances at birth are located: northern Africa, sub-Saharan Africa, Caucasus and central Asia, and southern Asia. In the Asian regions, under S2, the SRB is estimated to have peaked in southern Asia before 2020 and is projected to first decline in Caucasus and central Asia until 2040 and then rebound up to 1.070 (1.056; 1.093) in around 2060 before an ultimate return to natural level at the end of the century. The onset of the SRB inflation is expected to take place after 2030 in two countries of southern Asia (Afghanistan and Pakistan). In the African regions under scenario S2, the SRB is projected to peak in the mid-2040s in northern Africa at 1.071 (1.042; 1.132) and in the 2080s in sub-Saharan Africa at 1.039 (1.030; 1.059). This later rise of the SRB in sub-Saharan Africa corresponds to a later timing of fertility decline. This is notably the case for populated countries such as Nigeria and Tanzania where the SRB is expected to rise only after 2060 and to remain elevated till the end of the century in scenario S2 because the span of the projected SRB inflation process is 37 (15; 64) years.

By 2100, the differences in projected SRB under S1 and S2 are almost negligible in all regions, with the exception of sub-Saharan Africa where the SRB inflation is expected to take place mostly during the second half of the century. The projected rise in the SRB of sub-Saharan Africa from S1 to S2 in 2100 is projected to be modest but statistically different from zero at 0.004 (0.0004; 0.010).

Missing female births

The global burden of sex imbalances at birth during 2021–2100 is projected to be 5.7 (1.2; 15.3) million missing female births (CMFB) under S1. This number increases fourfold to 22.1 (12.2; 39.8) million under S2 (table 1, figures 2 and 3). Based on the first scenario, missing female births are projected to occur mostly before 2030, with the projected CMFB at 4.7 (1.1; 9.7) million during 2021–2030 and 1.0 (0.0; 6.0) million during 2031–2100. Under scenario S2, the projections for 2021–2030 are slightly higher than the first scenario at 5.9 (2.2; 11.4) million. However, the number of CMFB projected for 2031–2100 is much larger under S2 at 16.2 (7.8; 32.0) million.

Figure 2

Global and regional CMFB projections 2021–2100, by scenario. Regional CMFBs are stacked to obtain global CMFB. Regions are sorted by the timing of the imbalance (ie, the year in which the annual number of missing female births peaks). CMFB: cumulative number of missing female births; ENAN: the combination of countries in Europe, North America, Australia and New Zealand.

Figure 3

Projected SRB in 2100 and the CMFB during 2021–2100 for S1 and S2, by country. Countries are coloured by the levels of their SRB median estimates. radii of circles are proportional to CMFB for countries. CMFB: cumulative number of missing female births; SRB: sex ratio at birth.

Table 1

Global and regional CMFB for periods 1970–2020, 2021–2030, 2031–2100 and 2021–2100, by scenario

Under S2, we project that from 2031 till 2100, the cumulative deficit of female births in southern Asia and eastern Asia will be 5.1 (1.1; 12.0) million and 0.5 (0.0; 5.0) million CMFB, respectively, while sub-Saharan Africa may experience 8.6 (2.2; 22.1) million CMFB (table 1). In sub-Saharan Africa, sex imbalances at birth are projected to emerge only during the second half of the century, except for Uganda where the rise may take place earlier. As a result, almost all of the missing female births in sub-Saharan Africa until 2100 are projected to occur after 2050. The AMFB is projected to peak in the region during the 2080s at 137 (5; 662) thousand, declining to 103 (1; 624) thousand by 2100 (figure 2). This implies that sub-Saharan Africa would become a third locus of missing girls by 2100 (figure 3, Bottom), accounting for more than 38% of female births missing in the world during 2021–2100 under scenario S2.

Table 2 summarises the distribution of CMFB during 2021–2100 by scenario at country level. Among the 12 countries with strong evidence of existing SRB inflation, three countries are estimated to have completed sex ratio transitions by 2020: Republic of Korea in 2007 (95% uncertainty interval (2004; 2012)), Hong Kong (SAR of China) in 2013 (2013; 2014), and possibly Georgia in 2020 (2008; 2037). The end years for the remaining nine countries are projected to occur in 2020s (Albania, Armenia, Montenegro, Taiwan (Province of China) and Tunisia) and in 2030s (Azerbaijan, India, China and Vietnam). During 2021–2100, the projected CMFB in India and China are 2.5 (0.0; 8.0) million and 2.8 (0.0; 10.6) million, respectively. Under scenario S1, China and India are projected to contribute the majority of the CMFB around the world during 2021–2100, representing 49.5% (0.0%; 91.2%) and 43.8% (0.0%; 93.1%) of the total number, respectively.

Table 2

Scenario-based CMFB projections for period 2021–2100, by country

For the 17 countries at risk of inflation but without strong evidence today, only five countries in this group have already reached low fertility before 2021 (Bangladesh, Singapore, Morocco, Nepal and Turkey). Even though the start years for the five countries are estimated to be before 2021, the model also suggests low probability of inflation among them (Bangladesh at 77.5%, Morocco at 37.8%, Nepal at 81.9%, Singapore at 75.7% and Turkey at 49.1%). For the remaining 12 countries, the start years of inflation under scenario S2 are projected, in general, to be around 2030s for most Asian countries and around 2050 and later for all countries in sub-Saharan Africa. These projected start years correspond to the projected fertility decline profiles, with Asian countries expected to achieve low fertility before sub-Saharan Africa countries. Under S2, the CMFB distribution becomes more spread out and shift towards sub-Saharan African countries. The AMFB in Nigeria and Tanzania will be among the largest in the world during the second half of the century. In this scenario, China and India account in this scenario only for 13.8% (0.0%; 39.0%) and 12.2% (0.0%; 31.6%), respectively, of the global CMFB during 2021–2100. In contrast, Nigeria contributes 20.6% (0.0%; 52.4%) of the total, Pakistan 14.3% (0.0%; 38.1%) and Tanzania 7.2% (0.0%; 6.8%).

Discussion

It is crucial to be able to anticipate the potential course of sex imbalances at birth in order to strengthen policies against prenatal sex selection and to plan for the impact of future changes in sex structures across the world. To do so, we projected the SRB for all countries from 2021 to 2100 under two scenarios. The SRB model used here was based on a comprehensive database of 12,341 observations with 3.26 billion birth records from 204 countries, the experience of countries facing elevated SRB levels prior to 2021, and model assumptions regarding the sex ratio transition process. The Bayesian hierarchical time series mixture model used allows sharing information from data-rich country-years with those with limited or no data for estimating baseline levels and natural deviations. In addition, the model included an inflation probability and provides probabilistic SRB projections based on two scenarios. Based on the out-of-sample validation and simulation results, the model showed reasonable predictive performance.27

The strength of our projection procedures lies in the quality of the original database and the model we developed. The national SRB database is comprehensive and is based on rigorous assessment and quality checks.13 The Bayesian hierarchical time series mixture models enable (1) the estimation of SRB baselines on regional and national levels, (2) the use of the experience of countries affected over the last 40 years to model future change and (3) the incorporation of fertility declines as the covariate of future SRB change in countries with established son preference.

This study rests, however, on several model assumptions and is therefore subject to limitations.13 27 Limitations pertain to the estimation of baseline SRB and the estimation of sex-selective abortions. Regarding the estimation of baseline SRB, first, the only source of heterogeneity of baseline (natural) explicitly accounted for in the model is ethnicity groups approximated by region and country. There are numerous potential determinates of the baseline SRB (eg, age at conception, birth order, state of the sex chromosomes during pregnancy, maternal weight, socioeconomic conditions, maternal investment and the Trivers-Willard hypothesis).34 38–46 These factors are not taken into account here for lack of reliable and systematic information on these covariates globally. Second, factors related to conflict, crises (eg, a Chernobyl effect), and economic stress have been reported to cause short-term changes in SRB.43 45 47–51 Again, lacking specific information to include in our model, we do not capture these effects explicitly. Instead, the baseline SBR includes a general country-year-specific term to capture data-driven year-to-year fluctuations.27

Limitations pertaining to the estimation of sex-selective abortions include the absence of information on covariates to capture sex-selective abortion. Quality estimates of variables that might be relevant to sex-selective abortion, such as accessibility of modern reproductive technologies and to legal abortion, are not typically available. Instead, we used fertility decline as the only covariate in the model. In addition, no reliable measurement of son preference was available for 57 out of the 212 countries included in the study and accounting for 3.2% of the global births since 1970. We selected the countries at risk of SRB rise prior to model fitting, as opposed to incorporating the selection in the model, to resolve potential identifiability issues between natural fluctuations and inflation in the SRB.

There are also limitations specific to the scenario-based projections presented in this study. First, there is substantial uncertainty associated with the inflation predictions for countries prior to observing country-specific transition data due to the variability across countries that have started their SRB transition. Second, information that is available for extrapolation sex ratio imbalance at birth to countries without existing SRB inflation is limited, because the sex ratio transition has been observed in only 12 countries where most inflation processes are incomplete. Subsequently, long-term projections up to 2100 are subject to additional uncertainty not captured in our approach while we took into account uncertainty in future fertility trajectories. Lastly, we used the UN WPP birth projections for computing missing female births when alternative demographic projections incorporating our SRB scenarios may have given slightly different figures due to the declining female population and its impact on future births.

Our findings indicate that under the conservative scenario S1, where sex ratio transitions are included only for 12 countries with strong statistical evidence of past or ongoing SRB inflation, the cumulative global female births deficits are projected at 4.7 million (1.1; 9.7) during 2021–2030 compared with 1.0 (0.0; 6.0) during the rest of the century. This projected deficit is smaller than the figure for 1970–2020—estimated at 49.4 million missing female births—and stems from the projected gradual recovery from the skewed SRB in China and India in the coming years. This projected deficit is also significantly lower than previous estimates that put the future number of missing female births in 2021–2030 at 11.8 million.14 The use of our model signals a faster reduction of prenatal sex selection in these countries than previously anticipated.

The impact of the last phase in the sex ratio transition in the world remains nonetheless still large in absolute terms, with millions of predicted missing female births added to the already existing demographic deficit of women. The global figure is mostly driven by SRB trends in large countries, but SRB levels above 1.1 are also projected in other countries affecting a sizeable share of annual births. Even under scenario S1, immediate actions are needed in countries with ongoing sex ratio transitions—notably Albania, Armenia, Azerbaijan, China, India and Vietnam—to accelerate the return of the SRB back to normal levels and to mitigate the impact on these societies of the mounting surplus of men among adult populations and the resulting marriage squeeze.9 52

We provided an additional set of scenario-based projections that quantify the effect of SRB inflation in 17 countries at risk but without strong evidence of existing SRB inflation. Most of these 17 countries are characterised by higher fertility and son preference, and their sex ratio transitions would happen later in this century, when fertility approaches replacement level. Scenario S2 leads to a CMFB at 22.1 (12.2; 39.8) million by 2100. The assumptions made for scenario S2 remain hypothetical as they depend on the projected effect of fertility decline during the future decades in countries such as Afghanistan, Egypt, Nigeria, Pakistan or Tanzania. Yet, these projections offer useful simulations of the potential burden of future prenatal sex discrimination and highlight the need to anticipate the possible deterioration of the SRB in vulnerable countries if gender bias and son preference remain what they are today.

It should be added that projections in this paper are labelled as ‘scenarios’ because the experience of countries with past or ongoing sex ratio transitions, mostly in eastern Asia and eastern Europe, does not have to be predictive of the future evolution of countries belonging to different social and cultural environments with different cultures, for example, countries in sub-Saharan Africa. In this regard, we recommend using projections deriving from the conservative scenario S1—assuming no inflation for the 17 at-risk countries—for future population projections such as those conducted by the UN or the Global Burden of Disease study. While it is too early to project SRB inflation in countries where birth masculinity levels have so far remained constant, the scenario S2 helps visualise what might happen under the impact of prolonged fertility decline in these countries.

In view of the high level of heterogeneity in gender systems within large countries such as China, India and Viet Nam, national-level projections still represent an imperfect tool to simulate future imbalances. The well-documented subregional diversity of SRB levels within countries should encourage the development of subnational projections in such settings to provide additional insights into the geographic concentration of missing female births and future excess male populations.11 53–57

These findings underline the need to monitor SRBs in countries with son preference and to address the factors behind the persistence of gender bias in families and institutions. Measurement issues remain severe, notably in the absence of birth registration data, and regular efforts to estimate trends and differentials are required to monitor SRB dynamics in vulnerable countries. In addition, policies based on monitoring, advocacy campaigns, as well as direct and indirect measures to combat gender bias are required to slow down the rise of SRBs or to accelerate its decline.58 A broader objective relates to the need to influence gender norms which lie at the core of harmful practices such as prenatal sex selection. This calls for broader legal frameworks to ensure gender equality.1 59

Data availability statement

Data are available in a public, open access repository. All data relevant to the study are included in the article or uploaded as online supplemental information. All data used in this paper are publicly available or included in the article.

Ethics statements

Acknowledgments

We are grateful to Nicolas Todd for helpful comments and discussions.

References

Supplementary materials

Footnotes

  • Handling editor Sanni Yaya

  • Twitter @FengqingChao, @Patrick_Gerland, @krstoczg, @LeontineAlkema

  • Contributors FC and LA developed the Bayesian statistical model. FC carried out the analysis, drafted the initial manuscript, and prepared the appendix. FC and PG oversaw database construction, assessed, compiled and updated the database. CZG interpreted the results and its implications. All authors reviewed model results, edited the manuscript and agreed on the final version of the manuscript.

  • Funding Funding source was research grant R-608-000-125-646 from the National University of Singapore. The research was also supported by the University of Massachusetts Amherst.

  • Disclaimer The views expressed herein are those of the authors and do not necessarily reflect the views of the United Nations. Its contents have not been formally edited and cleared by the United Nations.

  • Map disclaimer The depiction of boundaries on the map(s) in this article does not imply the expression of any opinion whatsoever on the part of BMJ (or any member of its group) concerning the legal status of any country, territory, jurisdiction or area or of its authorities. The map(s) are provided without any warranty of any kind, either express or implied.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.