Abstract
Importance Deaths of parents and grandparent caregivers linked to social and health crises threaten child wellbeing due to losses of nurturance, financial support, physical safety, family stability, and care. Little is known about the full burden of all-causes and leading cause-specific orphanhood and caregiver death beyond estimates from select causes.
Objective To estimate 2000-2021 prevalence and incidence trends of all-cause orphanhood and caregiver death among children <18, by cause, age, race/ethnicity, and state.
Data Sources National Center for Health Statistics (NCHS) birth, death, race/ethnicity, and population data to estimate fertility rates and identify causes of death; 1983-1998 ICD-9 causes-of-death harmonized to ICD-10 classifications; 1999-2021 ICD-10 causes-of-death; CDC WONDER for state-specific estimates; and American Community Survey for grandparent population estimates.
Data extraction and synthesis We extracted U.S. population-level death, birth, population size, race, and ethnicity data from NCHS and attributed to each deceased individual the average number of children left behind according to subgroup-specific fertility rates in the previous 0-17 years. We examined prevalence and incidence of orphanhood by leading causes-of-death, including COVID-19, the leading 5 causes-of-death for 1983-2021, and additional leading causes for ages 15-44. We extended these to obtain state-level outcome estimates.
Main outcome measures National incidence and prevalence of orphanhood and caregiver death from 2000-2021, with orphanhood by year, parental cause-of-death and sex, child age, race/ethnicity, and state.
Results From 2000-2021, orphanhood and custodial/co-residing grandparent caregiver loss annual incidence and prevalence trends increased 49.2% and 8.3%, respectively. By 2021, 2.9 million children (4% of all children) had experienced prevalent orphanhood and caregiver death. Populations disproportionately affected by orphanhood included 5.0% of all adolescents; 6.5%, 4.8%, and 3.9% respectively of non-Hispanic American Indian/Alaska Native, non-Hispanic Black, and non-Hispanic White children; and children in New Mexico and Southern and Eastern States. Parental death due to drug overdose during 2020-2021 surpassed COVID-19 as the leading cause of incident and prevalent orphanhood during the COVID-19 pandemic.
Conclusions and Relevance Policies, programs, and practices aimed at orphanhood prevention, identification, and linkage to services and support of nearly 3 million bereaved children are needed, foremost prioritizing rapidly increasing overdose-linked orphanhood.
Key Point Question: What are U.S. trends in all-cause and cause-specific orphanhood and caregiver death among children <18?
Findings: From 2000-2021, orphanhood and caregiver loss incidence and prevalence increased 49.2% and 8.3%, respectively. By 2021, 2.9 million children (4% of all children) were affected. Populations disproportionately affected by orphanhood included 1.7 million adolescents ages 10-17; 6.5%, 4.8%, and 3.9% respectively of non-Hispanic American Indian/Alaska Native, non-Hispanic Black, and non-Hispanic White children; and children in New Mexico, Southern and Eastern States. Drug overdose was the leading cause of orphanhood during the COVID-19 pandemic.
Meaning: Evidence-based programs and policies are needed to prevent orphanhood and support these bereaved children.
Introduction
Deaths of parents and co-residing/custodial grandparent caregivers, in the context of intersecting social and health crises, threaten child wellbeing1 due to losses of nurturance, financial support, housing, and care.2–4 Orphanhood (death of one or both parents)5,6 and caregiver loss (death of co-resident/custodial grandparents) are considered Adverse Childhood Experiences (ACEs),7,8 and may have enduring consequences that persist well into adulthood.9–12 ACE’s are associated with increased chronic stress and lifelong risks of mental health illness, suicide, post-traumatic stress disorder, violence, insecure housing,13 and chronic and infectious diseases.14,15 These impacts often lead to ongoing needs for mental health, parenting, educational, and economic support for affected children, and for foster care or adoption services for children bereft of care.15
In the United States, mortality data informs public health prevention and response policies according to leading causes-of-death.16 However, one disproportionately affected population often remains largely invisible – those children left behind by adult deaths.2,6,17 One exception are children experiencing COVID-19-associated orphanhood, where real-time incidence estimates18–20 led to policies to support children, including recommendations in the National COVID-19 Pandemic Preparedness Plan21 and investments in some states,22 such as financial, bereavement, and mental health support for bereaved children and surviving caregivers. These established frameworks raise the possibility that new standards of care can be extended for any child experiencing orphanhood, regardless of cause.23
To inform evidence-based, compassionate prevention and response strategies for affected children and families, understanding the numbers, time-trends, causes, locations and disparities in all-cause orphanhood is essential. Reports show COVID-19 excess deaths caused a surge in orphanhood and caregiver death, affecting over 10.5 million children globally.24 In the U.S., the pandemic was linked to increased challenges including substance use, economic crises, mental health distress, and may have been linked to deaths due to overdose, suicide, and excessive alcohol use.25,26 Though reports have described the U.S. orphanhood burden related to HIV,27 maternal cancers,23 COVID-19,15 and child bereavement,28 up-to-date estimates on all-cause orphanhood do not exist.
Therefore, we estimate prevalence, incidence, and trends from 2000-2021 in all-cause orphanhood and caregiver death among children in the U.S. We further characterize leading causes of orphanhood and caregiver loss and identify whether the compound crises of the COVID-19 pandemic and drug overdose epidemic in 2020 and 2021 were associated with escalating orphanhood. We aim to identify populations disproportionately affected, by age, race/ethnicity, and state, to strengthen evidence-based response strategies. Our findings can contribute to identifying children at risk among population subgroups and geographic areas, and to targeting effective prevention and protection programming for affected children and parents/caregivers.
Methods
Overview
We adapted methods from Hillis et al19 to estimate the magnitude of orphanhood and caregiver loss among children ages 0-17 from all causes of deaths, according to Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER). We extracted subgroup-specific individual live birth and death data from National Center for Health Statistics (NCHS), along with population size estimates. We then attributed to each deceased individual the average number of children orphaned using subgroup-specific fertility rates in the previous 0-17 years. Our outcomes were orphanhood incidence and prevalence from 2000 through 2021, disaggregated by year, child age, and sex, race/ethnicity, cause-of-death of the deceased parent (numbers and percent of child population). We extended these estimates to include loss of co-residing/custodial grandparent caregivers, and state-level estimates of incidence and prevalence of orphanhood and grandparent caregiver loss (Supplementary Text (ST) S1, Fig. 2). We calculated 95% uncertainty intervals around estimates using bootstrap resampling (ST S1).
National mortality data
We extracted mortality data of U.S. residents by underlying cause-of-death from NCHS Vital Statistics for 1983-2021.29 The World Health Organization (WHO) defines underlying cause-of-death as “the disease or injury which initiated the train of events leading directly to death, or the circumstances of the accident or violence which produced the fatal injury”.30 Mortality records29 were mapped and aggregated to one of 53 caregiver loss causes-of-death that are closely related to the 52 rankable causes-of-death in the NCHS 113 cause-of-death list31, or ‘Other’ caregiver loss cause-of-death (ST S1). Cause-of-death classifications changed with the Tenth Revision of the International Classification of Disease (ICD-10), and we harmonized 1983-1998 cause-of-death data to ICD-10 classifications to avoid discontinuities in orphanhood estimates from 1998-1999 (ST S1).
Our cause-specific analyses for orphanhood estimates focused on leading causes of death in 2021: COVID-19, drug overdose, the remaining top five causes-of-death overall (heart disease, cancer, unintentional injuries, cerebrovascular diseases, chronic lung disease), and any additional causes of death in the top five for men and women ages 15-44 (suicide and homicide), as these ages include those more likely to be parents. We used a hierarachical approach for classifying cause of death; therefore, we consistently excluded drug overdose from each of these three categories: suicide, homicide, and unintentional injuries. Race and ethnicity in death certificates are typically reported by next of kin. We converted data that were incrementally reported by multiple races to standardised race categories (Table S7).
National orphanhood estimates
To estimate the number of children each death leaves behind, we obtained age, race, and ethnicity data of fathers and mothers from NCHS live birth records.29,31 We used corresponding 1990-2021 population sizes32 to calculate annualized age- and sex-specific national-level fertility rates by standardised race/ethnicity categories (ST S1). Before 1990, we assumed fertility rates of 1990 because population sizes by race/ethnicity were unavailable before 1990, and report corresponding sensitivity analyses (ST S2).
Incidence
We estimated annual incidence of orphanhood for a given year among children by multiplying the number of cause-specific deaths among adults stratified by sex, race/ethnicity, and five-year age bands in each year 1983-2021, with the expected number of children aged 0-17 born to each adult, with the corresponding characteristics in the same year. Calculations were performed by race/ethnicity groups and then aggregated to U.S. totals due to correlations between mortality and fertility rates by race/ethnicity (ST S1). Age of child was attributed by considering fertility rates in each of the current and preceding 17 years, after adjustment for pediatric survival probabilities (Table S4).
Prevalence
We estimated annual prevalence of orphanhood among children in any given year in 2000-2021 by cumulatively summing incidence estimates by one-year age of child in the current and previous 17 years, excluding children who would have since turned 18; calculations of orphanhood prevalence for 2000-2015 required including annual orphanhood incidence estimates for 1983-1999, in order to include the cumulative prior 17-year risk for experiencing prevalent parental death. We removed duplicate orphanhood counts by adjusting for those who had lost both parents (ST 1, Table S5). Incidence and prevalence rates were calculated relative to child population sizes in the corresponding year.
National grandparent caregiver loss
We estimated annual incidence of children losing grandparent caregivers in 1983-2021 by multiplying the number of cause-specific deaths in each year in adults aged 30 years and above stratified by sex and race/ethnicity with the corresponding proportions of adults co-residing with their grandchild(ren) and providing some or all basic needs for grandchild(ren), using U.S. Census American Community Surveys (ACS) 2010-2021 data (ST S1).33 We assumed only one child was affected by each grandparent caregiver death. We also developed estimates for the age of children affected (Fig S21).
State-level orphanhood and grandparent caregiver death estimates
We generated state-level estimates of cause-specific incidence and prevalence of orphanhood and caregiver death for 2021, using data from CDC WONDER.34,35 (ST S1). State-level orphanhood and caregiver death were calculated as above, with exceptions for male fertility rates for 2005-2015 (ST S1). We developed correction factors to adjust for bias resulting from data not stratified by race and ethnicity (ST S1).
Results
All-cause orphanhood and caregiver loss
In 2021, we estimate 498,249 (95% CI, 496,938-499,639) childr experienced incident orphanhood or caregiver death from any cause (0.72% of children, Table 1, Fig. 2). Most children lost a parent (81.6%), while the remainder lost a grandparent caregiver (18.4%). Before the COVID-19 pandemic, incidence of orphanhood, and of orphanhood/caregiver loss (combined) increased rapidly from 2000 to 2019 by 11.4% and 11.7%, respectively (Table 1). Through the pandemic, incidence of orphanhood, and of orphanhood/caregiver loss (combined) increased further from 2019 to 2021 by 39.5% and 33.6%, respectively. Prevalence of orphanhood, and of orphanhood/caregiver loss decreased slightly by 2.4% and 0.9% from 2000 to 2019, then increased by 10.1% and 9.3% during the period spanning the pandemic years (2019 to 2021). In 2021, we estimate 2,966,008 (95%CI, 2,963,030-2,968,931) children experienced prevalent orphanhood or caregiver death, representing 4.0% of all children.
Cause-specific trends
From 2000-2021, orphanhood caused by drug overdose increased markedly, as did orphanhood due to COVID-19 between 2019-2021. Orphanhood due to suicide decreased until 2010 and subsequently increased. Orphanhood due to unintentional injuries and homicide decreased until 2019, then increased until 2021; orphanhood due to diseases of the heart remained relatively unchanged until 2019, then increased. Only orphanhood due to malignant neoplasms consistently decreased from 2000-2021. By 2021, both incidence and prevalence of orphanhood due to fatal injuries – including drug overdose, suicide, homicide, and unintentional injuries (e.g. motor vehicle crashes) – exceeded orphanhood due to chronic diseases (Fig 1A, B and Table S1). Orphanhood caused by drug overdose rose substantially from 2012-2019, then jumped sharply during 2020-2021 (Fig 1C). From 2000-2019, orphanhood incidence due to drug overdose increased from an estimated 0.015% of children in 2000, to 0.06% in 2019, 0.085% in 2020 after the pandemic onset, and 0.09% in 2021. Drug overdose was the leading cause of orphanhood incidence and prevalence in 2020 and 2021, surpassing both chronic diseases and COVID-19 (Fig 1C, D). The disproportionate impact of fatal injuries on orphanhood was evident in 2021, with drug overdose contributing 17.5% to orphanhood incidence, yet only 3.1% to adult mortality (Fig 1E); similarly, unintentional injuries, suicide, and homicide contributed 8.3%, 5.0%, and 3.7% to orphanhood, versus 3.6%, 1.2% and 0.7%, respectively, to adult deaths.
(A) Estimated number of U.S. children newly experiencing orphanhood by any cause and over time. (B) Estimated number of U.S. children experiencing orphanhood in their lifetime (over the past 18 years) by any cause and over time. (C) Incidence rates of orphanhood among US children. (D) Prevalence rates of orphanhood among US children. (E) Main contributors to children experiencing orphanhood (parental death) versus the main contributors to adult deaths in 2021.
Disparities by age and race/ethnicity
From 2000-2021, children ages 10-17 years were 5.4 times more likely than children ages 0-4 to experience prevalent orphanhood (Fig 2A, Table S2). By 2021, orphanhood prevalence among ages 10-17 was 5.0%, affecting 1,704,671 children. Substantial disparities across race and ethnicity occurred, accentuated by the COVID-19 pandemic. By 2021, orphanhood affected 6.5% of non-Hispanic American Indian or Alaska Native children, 4.8% of non-Hispanic Black children, 3.9% of non-Hispanic White children, 1.7% of Hispanic children, and 1.9% of non-Hispanic Asian children (Fig 2C).
(A) Prevalence rates of orphanhood among US children aged 0-4 years, 5-9 years, and 10-17 years. (B) Prevalence rates of orphanhood among children by race and ethnicity of the deceased parent. (C) Prevalence rate of all-cause orphanhood by age and standardized race & ethnicity in 2021. (D) Time trends in sex and cause-specific incidence rates of orphanhood according to race among US children. On the y-axis are shown 2021 incidence rates of orphanhood among US children, by sex and cause of death of the deceased parent. On the x-axis are shown the differences in incidence rates in 2021 minus those in 2000, with positive differences indicating increases in incidence rates and negative differences indicating decreases in incidence rates. Numbers show the cause of orphanhood. Shapes show the sex of the deceased parent. The size of points indicates the contribution of each cause to new cases of orphanhood respectively among US children in 2021.
Causes underpinning these disparities differed by parental race/ethnicity and sex (Fig 2D, Table S3). Among non-Hispanic American Indian or Alaska Native children, leading causes of maternal orphanhood in 2021 were, in order, COVID-19, drug overdose, and chronic liver disease and cirrhosis; while leading causes of paternal orphanhood were chronic liver disease and cirrhosis, unintentional injuries, and diseases of the heart – all of which increased as causes of orphanhood since 2000. Among non-Hispanic Black children, leading causes of maternal orphanhood were COVID-19, diseases of heart, and drug overdose; while paternal orphanhood was primariy caused by diseases of heart, drug overdose, and homicide. In non-Hispanic White children, drug overdose deaths among fathers increased substantially over the past two decades; in Hispanic children, COVID-19 deaths among fathers predominated. In non-Hispanic Asian children, primary causes of orphanhood were heart disease, malignant neoplasms, and COVID-19.
State-level Estimates
In 2021, 33 states had >3% of all children experiencing prevalent orphanhood; Southeastern, Northeastern, Southern Border, and Midwest states were most affected (Table 2, Fig. 3, and Tables S4-S6). California, Texas, and Florida had the highest numbers of children experiencing orphanhood incidence and prevalence in 2021. West Virginia and New Mexico had the highest incidence (0.8-0.9%) and prevalence rates (4.5-5.0%) of orphanhood in 2021. Injury-associated parental deaths – including overdose, suicide, and/or unintentional injury – were among the top two causes of orphanhood prevalence in 47 states (Table 2). Drug overdose was the leading cause of orphanhood prevalence in 30 states, particularly affecting Southeastern, Southwestern, Southern Border, and Midwest states, and Southwestern States with large American Indian and Alaska Native populations.
Discussion
Over the past two decades, prevalence of all-cause orphanhood and caregiver death decreased slightly until 2019 and then increased markedly during 2020-2021 with intersecting crises of the overdose epidemic and COVID-19 pandemic.36,37 Parental death due to drug overdose was the leading cause of orphanhood incidence and prevalence during the pandemic. By 2021, over 2.9 million children – 4.0% of all children in the U.S. – had lost a parent or caregiver in their lifetime. Populations disproportionately affected by orphanhood included over 1.7 million adolescents ages 10-17 years (1 of every 20 adolescents); and children of non-Hispanic American Indian or Alaska Native, non-Hispanic Black, and non-Hispanic White race/ethnicities (approximately 1 of 15, 1 of 20, and 1 of 25 children, respectively). Orphanhood prevalence was highest in West Virginia, New Mexico, Mississippi, Louisiana, and Kentucky, affecting 1 of every 25 children. Compared to previous cause-specific reports (an estimated 97,376 children orphaned by HIV/AIDS (1998),27 157,183 by maternal cancers (2020),23 and 218,800 by COVID-19 in 2020-202238), data on all-cause orphanhood and co-residing grandparent caregiver death are over ten times greater, quantifying the full burden.
By 2021, orphanhood incidence and prevalence due to fatal injuries (drug overdose, suicide, homicide, and unintentional injuries), exceeded those linked to chronic diseases. Together, injuries and chronic diseases accounted for over half of all children orphaned. Fatalities due to injuries and chronic diseases among minoritized populations, are more common among younger adults still caring for children. For 47 states, drug overdose, suicide, and/or unintentional injuries were among the top two causes of orphanhood. Although drug overdose deaths were increasing before 2020, pandemic-linked increases of these deaths exceeded those forecasted; and the pandemic appeared to increase drug overdose risk.39,40 The evolving nature of drug supply, including the proliferation of illegally made fentanyl and resurgence of stimulants such as methamphetamine, has affected almost every state.41,42 By 2021, our data show drug overdose contributed 17.5% to orphanhood compared to 3.1% to mortality; establishing comprehensive care programs for children orphaned due to drug overdose could better mitigate the long-term negative effects that these children are disproportionately facing.
Our findings reveal social and health disparities in orphanhood rates and leading caregiver causes-of-death between age, sex, and racial/ethnic groups, and across geographies. These results can inform and contextualize prevention and response programming for children affected by orphanhood. For non-Hispanic American Indian or Alaska Native persons, for example, alcohol-related deaths were a leading cause, while non-Hispanic Black populations showed homicide among fathers as a leading cause. Many of these disparities are related to historical racism and derived practices that foster inequitable access to housing, education, and employment. Populations most affected in different states may have been exposed to such inequalities.43 In parts of the northeastern and midwestern U.S., rising unemployment, limited access to health and social services during the pandemic, increased access to a more lethal supply of illicit drugs or availability of other substances such as alcohol particularly among affected white men, may be linked to drug overdose emerging as the top cause of incident orphanhood in these regions.44 Circumstances of parental death influence grief-related psychopathology in surviving children,45 and evidence-based responses can improve short- and long-term outcomes.9,11,12 It is paramount to keep children in their families, whenever possible. This requires ensuring bereaved families receive support, and those needing kinship or foster care are served rapidly.46 Child resilience after parental loss can be strengthened by programs and policies that promote safe, stable, and nurturing relationships and environments, and that address childhood adversity, including preventing violence and abuse.46
Despite serious risks for bereaved children,47 timely, multifactored responses using an age- and life-stage approach may restore hope and build resilience.48 As grief is an immediate reaction that can extend years after parental loss, healing and support are urgent needs for children, who may be less equipped to recover than adults, particularly when the loss is sudden and unexpected.11,27 The governor of Utah49 and others50 have called for the addition to the death certificate of a checkbox to identify children living in the home of the deceased, so that bereaved children can be linked to support and services. Such responses typically provide parenting, economic, and education support, and are consistent with the WHO INSPIRE package for ending violence against children.51,52
Our study has several limitations. First, cause-specific estimates of children experiencing orphanhood and caregiver deaths are derived from mortality statistics on underlying causes-of-death, and may be underestimates53 for causes associated with erroneous or incomplete reporting, uncertainty in the chain of events preceding death, or coding limitations - such as for COVID-19, drug overdose, or suicide.54 Second, as we attributed only one child per grandparent caregiver death, our estimates of grandparent caregiver loss are minimum estimates. Third, we assumed current mortality is unrelated to historic fertility for years before 1990, and sensitivity analyses suggest this may lead to inaccuracies in orphanhood prevalence estimates until 2007 (ST S2). Fourth, publicly available state-specific data were partly suppressed due to small counts; therefore, we cannot exclude bias in state-specific orphanhood estimates. Finally, to characterise orphanhood prevalence in 2000-2021, we had to consider vital statistics since 1983 and consolidate changes in cause-of-death and race/ethnicity coding. Our sensitivity analyses, including detailed comparisons of NCHS death data to CDC WONDER vital statistics for each of the 21 years’ calculations (ST S2), suggest our incidence and prevalence estimates are robust, and should be interpreted as minimum estimates of orphanhood and caregiver death.
In conclusion, evidence-based policies and programs that provide healing and support for nearly 3 million children and adolescents in the U.S. who have experienced orphanhood and caregiver loss are urgently needed to reduce long term negative affects of this adverse childhood experience. This is especially relevant given unprecedented rates of drug overdoses, increasing the possibility that children who are living in a household where a parent or caregiver is negatively affected by substance use, may also experience the loss of a parent.52,55 Evidence highlights three essential components of orphanhood prevention and response that effectively promote their recovery and resilience: (1) prevent death of parents and caregivers through timely prevention and treatment of leading causes of death and ensured access to affordable compassionate health and mental health care for all; (2) prepare families to provide safe and nurturing alternative care; and (3) protect children affected by orphanhood and vulnerabilities, through grief and mental health counseling, and parenting, economic, and educational support that can be contextualized and delivered at scale.19,56
Ethical approval
All data used for this study are publicly available through U.S. Government websites.
Data sharing statement
All data used to calculate estimates of mortality and orphanhood are publicly available.
Code availability
Code to reproduce all analyses is freely available on Github version 1.1.2 under the GNU General Public License version 3.0 at the repository (https://github.com/MLGlobalHealth/orphanhood-caregiver-death-in-US-from-all-causes-of-mortality).
Competing interests
OR reports grants from the Bill & Melinda Gates Foundation, the EPSRC, the AIDSFonds, and the NIH during the conduct of this study.
Copyright statement
For the purpose of open access, the authors have applied a ‘Creative Commons Attribution’ (CC BY) license to any Author Accepted Manuscript version arising.
Author contributions
AV, SH and OR designed the study; AV, YC, ST, AB, LG, VKM, DW, SH and OR performed the analysis; LG, VKM, DW, JTU and SF checked analysis; AV, LC, LS, FA, GM, JL, LR, RN, CAN, SF, JTU, SH and OR oversaw data interpretation; AV, YC, ST, SH and OR wrote the first draft; All authors discussed the first draft and contributed to the final manuscript.
Supplementary Materials
S0 Supplementary Figures and Tables
S1 Extended methods
This Supplementary Text describes the data sources and statistical methods used to estimate cause-specific incidence and prevalence of orphanhood and grandparent caregiver loss among U.S. residents between 2000 and 2021. Our methods are overall similar to Hillis et al. [2021b,a], and extended here to estimate the magnitude of children aged 0-17 years experiencing the death of one or both parents from all causes or the death of a grandparent caregiver from all causes.
The United Nations International Children’s Emergency Fund (UNICEF) defines orphanhood as children experiencing the death of one or both parents [UNAIDS, 2009]. As in Hillis et al. [2021a], we considered mothers of age 15-66 years and fathers of age 15-94 years so the maximum ages of parents at birth of a child were respectively 49 and 77 years. Mortality data were recorded among U.S. residents, and for this reason, orphanhood estimates are restricted to children of U.S. residents. Grandparents play indispensable roles as caregivers for children [United States Census bureau, 2021, Pew Research Center, 2013, US Census Bureau, 2023b, Sowden et al., 2021, Ellis and Simmons, 2014]. We further considered to provide minimum numbers of children who lost a grandparent caregiver, defined as the death of a grandparent caregiver aged older than 30 years living with children in the same household. Mortality data were recorded among U.S. residents, and so grandparent caregiver death estimates are also restricted to children of U.S. resident grandparent caregivers. In Hillis et al. [2021a], we previously stratified grandparent caregivers based on the care they provided into primary grandparent caregivers who were either the guardians of children or co-resided and provided for basic needs of children, as well as secondary grandparent caregivers who co-resided and provided for some needs of children, following the characterizations of Ellis and Simmons [2014]. Here, we report throughout the sum of primary and secondary grandparent caregiver loss (Supplementary Fig. S1).
S1.1 National-level NCHS mortality data by rankable causes of death, 1983-2021
We obtained line-list mortality data on U.S. residents from the NCHS Vital Statistics portal for each year from 1983 to 2021 [National Center for Health Statistics, 2023]. Data were collected from 1983 onwards because the corresponding children who lost a caregiver in 1983 at age 0 were of age 17 in 2000, and so entered our estimation of orphanhood prevalence in 2000. For each mortality record we retained year of death, the corresponding codes of the underlying cause-of-death (ICD-9 code and 282 cause recode before 1999; ICD-10 code and 113 cause recode after 1999), and demographic data of the decedent including sex, age at death, and information on race and Hispanic origin [National Center for Health Statistics, 2022].
Information about the race and Hispanic origin of decedents in death certificates is typically reported by the surviving next of kin, or on the basis of observation in the absence of an informant [Rosenberg, 1999]. Race and Hispanic origin were reported in different formats across the study period [Ingram et al., 2003], which we mapped to standardized race and Hispanic origin categories as described in Supplementary Table S7. Specifically, we grouped individuals of Hispanic origin and all individuals of non-Hispanic origin by their race, i.e. ‘Hispanic’, ‘Non-Hispanic American Indian or Alaska Native’, ‘Non-Hispanic Asian or Pacific Islander’, ‘Non-Hispanic Black’, ‘Non-Hispanic White’, and we refer to the resulting categories as “standardized race & ethnicity” for simplicity. Data on Hispanic origin were not available for 1983 (see below). Individuals of other races in 1984-1991 and individuals of more than one race in 2021 were not coded consistently and not included in this study. We note this approach to harmonising race reporting over the study period did not necessarily make the primary data fully comparable. Earlier research shows inaccuracies are limited [Heron, 2021], which indicates that the incremental implementation of multiple race reporting in the U.S. is unlikely to have introduced notable bias in our orphanhood estimates.
The underlying cause-of-death is defined by the World Health Organization (WHO) as ‘the disease or injury which initiated the train of events leading directly to death, or the circumstances of the accident or violence which produced the fatal injury’ [World Health Organization, 2016]. For 1983-1998, underlying causes-of-death were in the data classified with the Ninth Revision of the International Classification of Diseases (ICD-9) [World Health Organization, 1977], and grouped further into 282 selected causes of death, termed ‘282 cause recode’ [National Center for Health Statistics, 1979]. For 1999-2021, underlying causes-of-death were in the data classified with the Tenth Revision of the International Classification of Diseases (ICD-10), and grouped further into ‘113 Selected Causes of Death’ [National Center for Health Statistics, 2020].
Due to our interest in orphanhood and caregiver loss, we defined 53 non-overlapping rankable underlying causes-of-death that we termed ‘caregiver loss causes-of-death’ and that apart from small modifications are identical to the 52 rankable underlying causes-of-death from the NCHS 113 Selected Causes of Death list [National Center for Health Statistics, 2020, Anderson et al., 2001]. The modifications are listed in Supplementary Table S8. Specifically, we re-categorized drug-induced causes-of-death into ‘Drug overdose’, joining the ICD-9/ICD-10 causes-of-death ‘Intentional self-poisoning by and exposure to nonopioid analgesics, antipyretics and antirheumatics’, ‘Intentional self-poisoning by and exposure to antiepileptic, sedative-hypnotic, antiparkinsonism and psychotropic drugs, not elsewhere classified’, ‘Intentional self-poisoning by and exposure to narcotics and psychodysleptics [hallucinogens], not elsewhere classified’, ‘Intentional self-poisoning by and exposure to other drugs acting on the autonomic nervous system’, ‘Intentional self-poisoning by and exposure to other and unspecified drugs, medicaments and biological substances’ (E950.0-E950.5; X60-X64); ‘Assault by drugs, medicaments and biological substances’ (E962.0; X85); ‘Accidental poisoning by and exposure to nonopioid analgesics, antipyretics and antirheumatics’, ‘Accidental poisoning by and exposure to antiepileptic, sedative-hypnotic, antiparkinsonism and psychotropic drugs, not elsewhere classified’, ‘Accidental poisoning by and exposure to narcotics and psychodysleptics [hallucinogens], not elsewhere classified’, ‘Accidental poisoning by and exposure to other drugs acting on the autonomic nervous system’, ‘Accidental poisoning by and exposure to other and unspecified drugs, medicaments and biological substances’ (E850-E858; X40-X44); and ‘Poisoning by and exposure to nonopioid analgesics, antipyretics and antirheumatics, undetermined intent’, ‘Poisoning by and exposure to antiepileptic, sedative-hypnotic, antiparkinsonism and psychotropic drugs, not elsewhere classified, undetermined intent’, ‘Poisoning by and exposure to narcotics and psychodysleptics [hallucinogens], not elsewhere classified, undetermined intent’, ‘Poisoning by and exposure to other drugs acting on the autonomic nervous system, undetermined intent’, ‘Poisoning by and exposure to other and unspecified drugs, medicaments and biological substances, undetermined intent’ (E980.0-E980.5; Y10-Y14).
We retained these four drug overdose causes-of-death sub-categories as separate drug overdose sub-groups for the purpose of data harmonization (see below). Correspondingly, we removed the related three drug overdose causes-of-death sub-categories respectively from ‘Intentional self-harm’, ‘Assault’ and ‘Accidents’. We renamed these resulting three categories as ‘Suicide excluding drug overdose’, ‘Homicide excluding drug overdose’ and ‘Unintentional injuries excluding drug overdose’ (see Supplementary Table S8). Finally, to simplify presentation as in Figure 1 in the main text, we primarily considered COVID-19 and the leading five causes-of-death between 1983 and 2021 due to our interest in all caregiver loss. Additionally, we included any cause in the top five for adult ages 15-44 years, as these ages specifically include younger adults who are more likely to be parents of younger children. The only one addition to our causes when taking this approach was ‘Homicide excluding drug overdose’. Supplementary Table S9 summarises our aggregation of the 53 rankable causes-of-death into these seven leading caregiver loss cause-of-death groups and ‘Other’ caregiver loss causes-of-death. The ‘Other’ caregiver loss causes-of-death comprise the remaining 46 rankable caregiver loss causes-of-death and any other causes-of-death that are not included in the NCHS 52 rankable causes-of-death.
We then mapped and aggregated line-list mortality records to the 53 rankable caregiver loss causes-of-death. This was done by mapping the 1983-1998 line list data using the 282 recodes to the 52 NCHS rankable causes-of-death as described in [National Center for Health Statistics, 1988] on page 282ff and page 290ff. For drug-induced causes, we used Table 2 in [National Center for Injury Prevention and Control (U.S.). Division of Unintentional Injury Prevention. Health Systems and Trauma Systems Branch. Prescription Drug Overdose Team, 2013]. Supplementary Table S9 summarises our approach.
For 1999-2021, we mapped ICD-10 113 cause recodes to NCHS 52 rankable causes-of-death based on Table A in [Curtin SC, 2023], and subsequently mapped the NCHS 52 rankable causes-of-death to the 53 rankable caregiver loss causes-of-death based on Supplementary Table S9.
Next, we aggregated line-list death records to annualised death counts by the 53 rankable caregiver loss causes-of-death for each year in 1983-2021, as well as by sex, age band (15-19, 20-24, 25-29, 30-34, 35-39, 40-44, 45-49, 50-54, 55-59, 60-64, 65-69, 70-74, 75-79, 80-84 and 85+ years) and standardized race & ethnicity of the decedent. Data on race and ethnicity were not available for 1983, and for this year we attributed mortality counts to standardized race categories according to the age-, sex- and cause-of-death-specific standardized race compositions in 1984. To harmonise cause-of-death data from 1983-1999 to the ICD-10 cause-of-death classifications and avoid discontinuities in our orphanhood estimates from 1998 to 1999, we used where available the comparability ratios in Table 1 of [Anderson et al., 2001]. Thirteen comparability ratios of rankable causes of death were not provided in [Anderson et al., 2001] when underlying estimations were considered imprecise, and in these cases we set the comparability ratios to 1. Supplementary Fig. S4 illustrates the aggregated annual mortality data among U.S. residents by standardized race categories, and Supplementary Fig. S5 by leading causes-of-death as relevant for caregiver loss and orphanhood.
S1.2 National-level live births data, 1969-2021
To calculate fertility rates and attribute children experiencing orphanhood to decendents between 1983 to 2021, we required for our approach to estimating orphanhood live births data from before 1983 because the children who lost a caregiver in 1983 at age 1-17 years were born between 1966-1982. However due to limitations in publicly available population size data, we only considered in the central analysis live births data from 1990 and assumed constant fertility rates before 1990. We investigated the sensitivity of our orphanhood estimates to this assumption using live birth records since 1980 together with corresponding population size estimates, and found that orphanhood prevalence estimates since 2008 were not affected by our assumptions on historic fertility, while in the sensitivity analysis orphanhood prevalence estimates for 2000-2007 were slightly lower due to overall lower fertility rates in the 1980s (Supplementary Text S2). Line-list live births with demographic data on both mothers and fathers were available and downloaded from the NCHS Vital Statistics portal for each year between 1969 to 2021 [National Center for Health Statistics, 2023]. We considered only live birth records to U.S. resident mothers for consistency with the mortality data available. Information on residency status of fathers was not available, and we assumed fathers were also U.S. residents in the retained birth records associated with U.S. resident mothers. For each live birth, we retained year of birth, the age of mothers and fathers, and information on race and Hispanic ethnicity [National Center for Health Statistics, 2023].
Age of mothers and fathers was reported by single year of age. Following Hillis et al. [2021a], we considered live births to women aged 15-49 years and men aged 15-77 years, to match the mortality data of women aged 15-66 years and men aged 15-94 years. In total, less than 0.01% of line-list records were removed from further analysis because parent were outside of these age ranges. Age information was mapped to 5-year age bands 15-19 years, …, 45-49 years for mothers and age bands 15-19 years, …, 50-54 years, and 55 years and older (55+ years) for fathers.
Information about the race and Hispanic origin of mothers and fathers was between 1969-2021 reported in different standards due to the revision of the U.S. certificates of live births on race in 1989 and 2003 [National Center for Health Statistics, 2006]. Information on the race of mothers and fathers is publicly available since 1969, and information on their ethnicity since 1978. For 1978-2021, we mapped available information on race and ethnicity to the same standardized race & ethnicity categories used to stratify the mortality data, according to the mapping described in Supplementary Table S10.
Next, we aggregated line-list live birth records to annualized live birth counts to mothers (in age bands 15-19, …, 45-49 years) and fathers (in age bands 15-19, …, 55+ years) for each year in 1969-1977, and stratified further by standardized race & ethnicity categories for each year in 1978-2021. Before 1985, data reporting was incomplete for some U.S. states and we used NCHS sampling weights reported in each year (see Appendix A of National Center for Health Statistics [1978]) to extrapolate reported live birth counts to state populations. Supplementary Fig. S6 illustrates the aggregated annual live birth counts by mothers and fathers and standardized race & ethnicity categories since 1978.
S1.3 National-level population size data, 1990-2021
To calculate fertility rates, we further required population size data of U.S. residents in the age-, sex- and standardized race & ethnicity categories of the live births data. We obtained from the CDC WONDER Vintage bridged-race postcensal & ethnicity population size estimates for each year in 1990-1999 and 2000-2020 [CDC, 2021a], and single race population size estimates for 2021 [CDC, 2021b]. Data were extracted by 5-year age bands (15-19 years, …, 70-74 years) and single years of age from 75 to 77 years, and data for individuals aged 55-77 years were summed into a single age band as required for the purposes of our analyses. Race-specific population size estimates were aggregated to the standardized race & ethnicity categories in Supplementary Table S11.
Bridged-race & ethnicity population size data were not publicly available prior to 1990. For sensitivity analyses, we also obtained national-level population size estimates by 5-year age bands, sex, and U.S. states without race & ethnicity stratification for each year in 1969-1989 from the U.S. National Cancer Institute Surveillance, Epidemiology, and End Results (SEER) Program [National Cancer Institute, 2023].
S1.4 Estimates of national-level fertility rates, 1990-2021
We calculated age-, sex- and standardized race & ethnicity-specific fertility rates in each year in 1990-2021 for women in one of the age bands a ∈ {15 − 19, …, 45 − 49} years and men in one of the age bands a ∈ {15 − 19,, 55 − 77} years according to
where the number of live births and population sizes in each strata are denoted by By,a,s,rand Py,a,s,r respectively. Calculations were done for each of the five standardized race categories ‘Hispanic’, ‘Non-Hispanic American Indian or Alaska Native’, ‘Non-Hispanic Asian or Pacific Islander’, ‘Non-Hispanic Black’ and ‘Non-Hispanic White’. In the central analysis, we assumed the same fertility rates in each year in 1966-1989 as in 1990, and considered alternative assumptions in several sensitivity analyses (see below). Supplementary Fig. S7 shows the calculated U.S. fertility rates by standardized race & ethnicity.
Mortality rates were calculated in analogy to Eq. (1). We found correlations by standardized race & ethnicity between fertility and mortality rates that changed primarily by age of mothers and fathers, and less so over calendar years (Supplementary Fig. S8 for women); findings were similar in men. These correlations prompted us to estimate national-level incidence and prevalence of orphanhood by standardized race & ethnicity and then sum the standardized race & ethnicity-specific estimates to obtain national-level estimates.
S1.5 Estimates of national-level orphanhood, 1983-2021
We estimated the number of children who newly experienced orphanhood in year y, y = 1983, …, 2021, from the population-level mortality records of U.S. residents in year y and the number of children each decedent was expected to leave behind. Our approach follows the methods described in [Hillis et al., 2021a, Unwin et al., 2022], with modifications to account for causes-of-death.
We obtained the expected number of children per one U.S. resident of age a years, sex s, and standardized race & ethnicity r in year y who are of age b = 0, 1, …, 17 years (denoted with Cy,a,g,r,b) by multiplying the corresponding fertility rates of Eq. (1) with pediatric survival probabilities of children born in year y − b and surviving until age b + 1 (denoted with psurvive). Specifically,
where y ranges from 1983 to 2021, the single year of age of mothers ranges from a ∈ {15, …, 66} years, the single year of age of fathers ranges from a ∈ {15, …, 94} years, the standardized race & ethnicity categories are as described before, and b = 0, 1, …, 17 years. We obtained the pediatric survival probabilities from the child mortality data in [United Nations Population Division, 2020] and use in Eq. (2) the fertility rates of the age band that includes age a − b. We then estimated the number of children aged b who newly experienced in year y the death of a parent s of age specified in 5-year age bands a′ ∈ A = {15 − 19, · · ·, 80 − 84, 85 +years}, and standardized race & ethnicity r who died of caregiver loss cause-of-death c by
Due to the correlations between standardized race & ethnicity specific fertility and mortality rates, we calculated Eq. (3) for each standardized race & ethnicity. In Eq. (3), we assumed that population-level fertility rates are not correlated with population-level mortality rates, which may lead to upwards or downwards bias in orphanhood estimates that we explored in several sensitivity analyses (see below). We also assumed that parents and their children have the same standardized race & ethnicity.
Since orphanhood considers children who experienced the death of their mother, father or both, we are interested in the sum of Eq. (3) for both mothers and fathers, but need to subtract children who lost their other parent in the previous b − 1 years, or who lost their other parent in the same year y. Throughout, we assumed that the other parent 1 − s is in the same age band a′and of the same standardized race & ethnicity as parent s. The probability that the other biological parent of the children in Eq. (3) died in the same year y is based on standard life table calculations [Shryock et al., 1975, Lee and Wang, 2003, Arias et al., 2022], specifically the probability that an individual of sex 1 − s and standardized race & ethnicity r died in year y between (continuous) age a and a + n conditional on survival up to age a, where n = 5 corresponds to the width of the 5-year age bands considered. This mortality hazard is approximated using midpoints x = (a + a + 5)/2 in each age interval, through
where for ease of readability we have suppressed y, s and r, and nPa and nDaare respectively the estimated population sizes and observed death counts by the end of the corresponding calendar year in each 5-year age band a′ = [a, a + 5). The intermediate quantities in (4) are the approximated (unknown) mortality probability density function nf (x) at age midpoint x and (unknown) survival function S(a) up to age a, that can be expressed in terms of the age-specific mortality rate nqa, the proportion of individuals alive at age a and who die before reaching age a + n in the corresponding calendar year. It is standard to estimate nqa with (nDa)/(nPa + 1 nDa), from which Eq. (4) follows [Shryock et al., 1975, Lee and Wang, 2003, Arias et al., 2022]. Using Eq. (4), we estimated the number of children aged b who newly experienced in year y the death of one parent due to cause-of-death c and the death of the other parent due to any cause with
where x is the midpoint age in age band a′. In Eq. (5), we assumed that deaths among parents occurred independently of each other and ignored correlations of deaths among parents by the same cause-of-death such as COVID-19 [Hillis et al., 2021a], as well as correlations of deaths among parents who died of different causes-of-death. Similarly, we estimated that the number of children aged b who newly experienced in year y the death of one parent s due to cause-of-death c and the death of their other parent 1 − s due to any cause in any of the previous i = 1, …, b − 1 years with
With these considerations, we estimated the number of children aged b who newly experienced orphanhood in year y = 1983, …, 2021, by the death of one or both parents of age a′ and standardized race & ethnicity r who died of cause-of-death c with
Eq. (7b) subtracts the children who lost in year y a parent due to cause c and the other parent due to any cause, which are counted twice in Eq. (7a). Without line-list family data and working from individual-level live birth and death statistics, the two terms Onew-double, Onew-double are not identical and for this reason we subtracted the average of both. Eq. (7c) subtracts the children who already lost the other parent in previous years. In our previous studies on COVID-19 associated orphanhood, we did not consider the possibility of reinfection with COVID-19 and for this reason did not subtract children who already lost the other parent due to COVID-19 in previous years in these studies [Hillis et al., 2021a]. Further arguments show that across ages, standardized race & ethnic groups, and caregiver loss causes-of-death, the values in Eq. 7b-7c are approximately equal to orphanhood prevalence divided by four, and so do not exceed 1% of the average values in Eq. 7a.
To estimate the prevalence of orphanhood in year y = 2000, …, 2021, we accrued the number of children who newly experienced orphanhood in the previous 17 years and current year y, while accounting for ageing. Specifically, for each calendar year since 2000, we estimated the total number of children aged b = 0, …, 17 years in calendar year y and of race & ethnicity r who experienced orphanhood in their lifetime by cause-of-death c in one or both parents with
which sums over children who newly experienced orphanhood at younger ages in previous years.
Following Eqs. (7-8), we derived additional key quantities such as the number of children aged b in calendar year y = 1983, …, 2021 and standardized race & ethnicity r who newly experienced orphanhood in year y by parental sex s and cause-of-death c,
We chose the specific form of Eq. (9) over simpler approaches to ensure that the sum of Onew and equals Eq. 7. Further, we derived the number of children aged b in calendar year y = 2000, …, 2021 and standardized race & ethnicity r who experienced orphanhood in their lifetime by parent sex s and cause-of-death c by
and the number of children of standardized race & ethnicity r who experienced orphanhood in calendar year y = 2000, …, 2021 in their lifetime by cause-of-death c in one or both parents by
All other total numbers reported in this manuscript are aggregations of Eqs. (7-11).
S1.6 Estimates of national-level grandparent caregiver loss, 1983-2021
We estimated the number of children who newly experienced grandparent caregiver death in year y, y = 1983, …, 2021 from the population-level mortality records in year y of U.S. residents aged 30 years and above (30+), and assuming that each decedent who is estimated to have been a grandparent caregiver leaves a minimum of one grandchild behind. Our approach follows the methods described in [Hillis et al., 2021a], with modifications to estimate the age of children who experience the death of a grandparent caregiver.
Starting from 2010, the American Community Survey (ACS) [Ellis and Simmons, 2014] collects data on the total number of adults aged 30+ years living with grandchildren of age 17 or under in the U.S., and proportions of these by sex, and separately by bridged-race and Hispanic origin [United States Census bureau, 2021].
For each year y = 2010, …, 2021, we estimated the number of adults aged 30+ years living with grandchildren of age 17 or under for each sex and standardized race & ethnicity by multiplying totals and proportions made available by the ACS. We then divided these numbers with the respective population size estimates to estimate sex- and standardized race & ethnicity-specific proportions of adults aged 30+ years who live with their grandchildren of age 17 or under, which we denoted with γy,s,r. For y = 1983, …, 2009, we assumed that γy,s,r is the same as in 2010. We further assumed that each grandparent caregiver leaves upon death a minimum of one child behind (corresponding to Eq. (2)) and estimated the minimum number of grandchildren who newly experienced grandparent caregiver death in year y = 1983, …, 2021, with a U.S. resident grandparent caregiver aged 30+ years, sex s and standardized race & ethnicity category r who died of leading cause c by
The ACS surveys did not collect data on whether both grandparents are alive and live with their grandchildren of age 17 or under, and for this reason we did not adjust Eq. (12) further for loss of other grandparents. We investigated in Sensitivity analyses our assumption that γy,s,r was approximately constant from 1983 to 2010 using longitudinal United Nations Population Division data on Households and Living Arrangments of Older Persons for the U.S., which suggested that the proportion of older persons who live with children or who are the primary caregivers of children has in the U.S. remained fairly constant since 1990 (Supplementary Fig. S23).
To estimate the minimum number of children who experienced grandparent caregiver death in their lifetime, we additionally need disaggregations of Eq. (12) by single year of age b = 0, …, 17 years. For the central analysis, we assumed that the age composition of Gnew is the same as the age composition of children who lost parents older than 30 years, are obtained from Eq. (7) by summing over the years y = 2000, …, 2021, and c is one of the leading caregiver loss causes-of-death. Supplementary Fig. S9 shows that the age compositions of children who lost parents older than 30 years
differ across causes-of-death while they are relatively stable across standardized race & ethnicity. It is plausible that the age composition of Gnew may differ from theage composition of children who lost parents older than 30 years, and we explored alternative approaches to Eq. (13) in sensitivity analyses; none of these had a considerable impact on our overall estimates.
To obtain minimum estimates of caregiver loss, we then summed Eq. (9) and Eq. (13).
S1.7 Uncertainty quantification in national estimates
This section describes how we derived uncertainty intervals in the national-level estimates of orphanhood incidence and prevalence (Eqs. (7-11)) and grandparent caregiver loss (Eq. (13)). To capture uncertainty, we followed the same phenomenological approach as in [Hillis et al., 2021a], and added Poisson noise around mortality, natality, and population size data for all year, sex, age, standardized race & ethnicity, and cause-of-death strata. We generated 10,000 Poisson noise data sets, and in terms of the uncertainty in mortality data, additionally considered uncertainty related to the harmonization of causes-of-death counts from 1983-1999 to the ICD-10 classification to avoid discontinuities in orphanhood estimates between 1998 and 1999. Specifically, we resampled the comparability ratios 10,000 times, assuming these were normally distributed with the standard deviations reported in Table 1 of [Anderson et al., 2001] where available, and multiplied these ratios with the Poisson noise mortality data. Then, we repeated orphanhood incidence and prevalence estimates on each data set, and calculated 2.5% and 97.5% quantiles to generate 95% uncertainty intervals around central estimates.
In the ACS data, we accounted for uncertainty in the estimated number of adults aged 30+ years living with grandchildren and in the attribution to sex, and standardized race & ethnicity groups using published uncertainty ranges. The ACS provides for each of their published estimates margins of error that correspond to 90% confidence intervals. We converted the 90% margins of error into standard deviations and bootstrap re-sampled the available totals and separate proportions 1,000 times assuming normal distributions around their central estimates and the corresponding standard deviations, multiplied the bootstrap resampled totals and proportions, and then divided by the corresponding population sizes. Across years, we found these sources of uncertainty amounted on average to 95% bootstrap intervals on the order of ±0.79% of the central estimate of the number of Hispanic women aged 30+ years who live with grandchildren, ±8.23% among Non-Hispanic American Indian or Alaska Native women, ±1.93% among Non-Hispanic Asian women, ±0.97% among Non-Hispanic Black women and ±0.69% among Non-Hispanic White women; and similarly for men. We similarly generated 10,000 resampled data sets, repeated grandparent caregiver loss incidence and prevalence estimates on each data set, and then calculated 2.5% and 97.5% quantiles to generate 95% uncertainty intervals around central estimates.
S1.8 Estimates of state-level orphanhood, 2021
To guide policy at the state level, we further estimated orphanhood and grandparent caregiver loss in each of the 50 U.S. states and District of Columbia. We focussed on estimating the number of children who newly experienced orphanhood (incidence) and the number of children who experienced orphanhood in their lifetime (prevalence) in each state in 2021, by the leading caregiver loss causes-of-death (see Supplementary Table S9). To obtain state-level orphanhood prevalence estimates for 2021, as before we had to accrue and age state-level orphanhood incidence estimates from 2004, and to obtain state-level fertility rates from 1987. Overall, we proceeded analogously as for the national estimations, but obtained data from CDC WONDER [CDC, 2023] as state-level mortality data were not publicly available from NCHS from 2005 onwards [National Center for Health Statistics, 2023]. Counts below 10 were suppressed in CDC WONDER and for this reason we performed all estimations without stratification by standardized race & ethnicity. We adjusted for data suppression, and known biases in this approach that are due to correlations in mortality and fertility rates across standardized race & ethnicity. We note these are modelled adjustments and we cannot exclude that this resulted in bias in state-specific orphanhood estimates.
We extracted annual death counts by state, sex, age bands (15 − 19 years, …, 95 − 99 years, > 100 years), and cause-of-death (ICD-10 113 Selected Causes of Death) from the CDC WONDER mortality portal [CDC, 2023] from 2005 to 2021; and combined these data with the previously described, cause-specific NCHS mortality data sets from 2004 for which information on the state of residence of the decedent was publicly available. Throughout, we focused only on the eight leading caregiver-loss causes of death chronic liver disease and cirrhosis, COVID-19, diseases of the heart, drug overdose, homicide excluding drug overdose, malignant neoplasms, suicide excluding drug overdose, and unintentional injuries excluding drug overdose. For the leading caregiver-loss causes of death, the average discrepancy in the state-, sex-, age-, and cause-of-death-specific data from CDC WONDER was, when aggregated to national-level and compared to the corresponding NCHS mortality data, 22.9% among women and 13.5% among men. As a first step, we imputed suppressed entries with a value of 2, which resulted in average discrepancies of 6.80% among women and 2.53% among men. To account for these discrepancies, as a second step, we adjusted the state-level mortality counts according to
where y in 2005 to 2021, a′ denotes the age band of the parent, s the sex of the parent and c the leading caregiver loss cause-of-death of the parent, and Dy,a′,s,c are the national-level deaths derived from the line-list NCHS data in Eq. (3), and Dy,a′,s,l,c are the state-specific deaths derived from CDC WONDER for each of the 50 U.S. states and District of Columbia, indexed by l. The suppression adjusted, cause-specific, state-level mortality counts matched the cause-specific, national-level mortality counts well, except when all deaths Dy,a′,s,c,l were entirely suppressed across all states, which was the case for deaths due to homicide in women aged 55 years or older (Supplementary Fig. S10, compare pink versus blue). Overall, these suppression adjustments did not noticeably change the contribution of causes-of-death to mortality, because the majority of death counts were not suppressed when the data were aggregated to leading caregiver loss causes-of-death (Supplementary Fig. S11).
We next extracted live birth counts for mothers by state and age bands (15 − 19, …, 44 − 49 years) from the CDC WONDER natality portal [CDC WONDER, 2022] from 2005 to 2021; and combined these data with the previously described NCHS live birth data sets from 1987 to 2004 for which information on the state of residence of mothers was available. We proceeded analogously for fathers using the age bands 15 − 19, …, 44 − 49, 50 − 54, 55+ years from 2016 to 2021, as natality records were not available by demographic characteristics of fathers from 2005 to 2015. Data suppression was not an issue for the stratifications required to estimate orphanhood: for the years 2005-2021, the average discrepancy in the CDC WONDER and NCHS natality data sets was 0.026% for mothers and 0.032% for fathers (as compared to male NCHS live birth data with ages up to 77 years). To compute fertility rates at the state level, we further extracted age- and sex-stratified annual population size estimates for each state from 1987 to 1989 from the NIH National Cancer Institute Surveillance, Epidemiology, and End Results (SEER) database [National Cancer Institute, 2023]; and CDC WONDER from 1990 to 2021 [CDC, 2021c, 2022]. State-level fertility rates were calculated as in Eq. (1) when live births counts were available, but without stratification by race & ethnicity. For mothers, live birth counts were not publicly available for women aged 44-49 years in several years and 8 states (Alaska, Delaware, Montana, North Dakota, South Dakota, Vermont, West Virginia and Wyoming), and in these cases we interpolated state-specific fertility rates with locally estimated scatterplot smoothing (LOESS) as implemented in the R stat package version 4.2.3 with span argument 0.85. For Wyoming, no data were publicly available after 2019 and so interpolation was not possible and we used 2019 values. For fathers, live births counts were not publicly available from 2005 to 2015, and for these years we estimated state-specific fertility rates with LOESS with span argument 0.85, and using NCHS data from 2000 to 2004 and CDC WONDER data from 2016 to 2020 (Supplementary Fig. S12).
We then estimated state-specific incidence and prevalence of orphanhood from the suppression-adjusted state-level, cause-specific mortality counts and partially imputed fertility rates as outlined in Eqs. (7-11), with one modification. As described above the state-level mortality and fertility data were not disaggregated by race & ethnicity. To account for potential bias arising from correlations in mortality and fertility rates across race & ethnicity,
Overall, the correction factors νy,s,c tended to be close to one, except for Homicide excluding drug overdose in women and ‘Other’ caregiver loss causes-of-death (Supplementary Fig. S13). We then adjusted the state-specific estimates as follows,
we compared the national-level estimates of the number of children who newly experienced orphanhood in (9) to the state-specific estimates with the correction factors
and used Eq. (16) to calculate the state-specific analogues to Eqs. (8-11).
We estimated the number of children who newly experienced grandparent caregiver death in year y, y = 2004, …, 2021 from the state-level mortality records in year y of U.S. residents aged 30 years and above (30+), and assuming that each decedent who is estimated to have been a grandparent caregiver in each state leaves a minimum of one grandchild behind. Overall, we proceeded as for the national-level estimation of grandparent caregiver loss as ACS also published state-specific data on adults aged 30+ years who live with grandchildren of age 17 or under [United States Census bureau, 2021]. For each year y = 2010, …, 2021, we calculated the expected number of adults aged 30+ years who live with their grandchildren of age 17 or under by multiplying the corresponding total co-resident numbers with the sex-specific proportions for each state; and then divided these with corresponding population sizes to obtain the corresponding proportions γy,s,l. Overall, these proportions were considerably more uncertain than in the national-level analysis due to smaller ACS sample sizes (see below). For y = 2004, …, 2009, we assumed that γy,s,lis the same as in 2010. We further assumed that each grandparent caregiver leaves upon death a minimum of one child behind (corresponding to Eq. (2)) and estimated the minimum number of grandchildren who newly experienced grandparent caregiver death in year y = 2004, …, 2021, with a U.S. resident grandparent caregiver aged 30+ years, sex s and state category l who died of leading caregiver loss cause-of-death c by
To estimate the minimum number of children who experienced grandparent caregiver death in their lifetime in each state, we additionally need disaggregations of Eq. (17) by single year of age b = 0, …, 17 years. For the central analysis, we assumed that the age composition of Gnew is the same as the age composition of children who lost parents older than 30 years in the same state,
are obtained from Eq. (16) by summing over the years y = 2004, …, 2021, and c is one of the leading caregiver loss causes-of-death. Supplementary Fig. S14 shows that the age compositions of children who lost parents older than 30 years differed across causes-of-death while they are relatively stable across states.
S1.9 Uncertainty quantification in state-level estimates
This section describes how we derived uncertainty intervals in the state-level estimates of orphanhood incidence and prevalence (Eqs. (7-11)) and grandparent caregiver loss (Eq. (13)). Overall, the approach here is similar to the national-level uncertainty analysis.
We added Poisson noise around state-specific mortality, natality, and population size data for all year, sex, age, and cause-of-death strata. We generated 10,000 Poisson noise data sets, repeated orphanhood incidence and prevalence estimates on each data set, and then calculated 2.5% and 97.5% quantiles to generate 95% uncertainty intervals around central estimates.
To estimate grandparent caregiver loss, we accounted for uncertainty in the estimated number of adults aged 30+ years living with grandchildren by U.S. state using published uncertainty ranges. The ACS provides for each of their state-specific estimates margins of error that correspond to 90% confidence intervals. We converted the 90% margins of error into standard deviations and bootstrap re-sampled the available totals and separate proportions 10,000 times assuming normal distributions around their central estimates and the corresponding standard deviations, and then multiplied the sampled totals and proportions. As expected given available ACS sample sizes, across years we found that uncertainty in the proportions γy,s,l was in some states considerably larger than for the national race & ethnicity-specific analysis, with 95% bootstrap intervals of up to ±13.5% of the central estimates (Supplementary Fig. S15). We similarly generated 10,000 resampled data sets, repeated grandparent caregiver loss incidence and prevalence estimates on each data set, and then calculated 2.5% and 97.5% quantiles to generate 95% uncertainty intervals around central estimates.
S1.10 Estimates of racial & ethnic groups impacted by leading causes of caregiver loss in U.S. states, 2021
Characterizing the racial & ethnic groups of children that are impacted by orphanhood and grandparent caregiver loss at state-level is challenging due to limits in publicly available data. We focused on estimating and comparing 2021 orphanhood and grandparent caregiver loss prevalence rates by standardized race & ethnicity in the 10 U.S. states with highest orphanhood prevalence rates and for the primary (first-ranked) caregiver-loss cause-of-death only, West Virginia (‘Drug overdose’), New Mexico (‘Drug overdose’), Mississippi (‘Unintentional injuries’), Louisiana (‘Drug overdose’), Kentucky (‘Drug overdose’), Tennessee (‘Drug overdose’), Alabama (‘Diseases of heart’), Oklahoma (‘Unintentional injuries excluding drug overdose’), Ohio (‘Drug overdose’) and Florida (‘Drug overdose’). Throughout, we obtained publicly available data from CDC WONDER [CDC, 2023] from 2005 onwards [National Center for Health Statistics, 2023] and adjusted for data suppression, but note these are modelled adjustments and we cannot exclude that this resulted in bias in state-specific orphanhood estimates.
Annualized live birth counts were extracted as described for the state-level analysis for 2005 to 2019, but now also stratified by bridged-race and Hispanic origins, and then combined with NCHS live births data from 1987 to 2004. For 2020-2021, live birth counts to women were not available by bridged-race at state-level. Analogously, annual death counts for the leading caregiver loss causes-of-death were compiled as described above for 2005 to 2021, and stratified into the standardized race & ethnicity categories in Supplementary Table S7, and then combined with the NCHS mortality data sets from 2004. Suppressed values were imputed by 1 and adjusted further for suppression as in Eq. (14). We followed previous criteria on data reliability and excluded from consideration for each state those race & ethnicity groups that had more than two age bands with fewer than 20 live birth counts Mathews and Hamilton [2019]. Supplementary Table S6 lists the corresponding strata as ‘small populations’, and also describes the average completeness of live birth records for the years 1995-2004 when directly comparable line-list live birth records were also available from NCHS. For the remaining year-state-sex-race & ethnicity strata, we considered the corresponding mortality data for the leading caregiver-loss cause-of-death. Again, we excluded from consideration those standardized race & ethnicity groups that had more than two age bands with fewer than 20 death counts Mathews and Hamilton [2019]. Supplementary Table S6 lists the corresponding strata as ‘small death counts’, and also describes the average completeness of death records for the years 1999-2004 when directly comparable line-list death records were also available from NCHS. Supplementary Fig. S16 illustrates the state and standardized race & ethnicity specific live birth data from CDC WONDER, the line-list live birth data from NCHS and the state and standardized race & ethnicity specific live birth data for the three states with highest orphanhood prevalence rates in 2021. Similarly, Supplementary Fig. S17 illustrates the state and standardized race & ethnicity specific mortality data from CDC WONDER for the primary caregiver loss cause-of-death, the line-list live death from NCHS and the adjusted state and standardized race & ethnicity specific mortality data for the three states with highest orphanhood prevalence rates in 2021. Population-size estimates for each state were obtained from the CDC WONDER population size portal for 1990 to 2021 [CDC, 2021c, 2022] by bridged-race, and converted to the standardized race & ethnicity categories described in Supplementary Table S11.
We estimated race & ethnicity- and state-specific incidence and prevalence of orphanhood as outlined in Eqs. (7-11), except that female race & ethnicity- and state-specific fertility rates in 2020-2021 were as assumed to be as in 2019 due to limitations in publicly available data. To investigate estimation accuracy, we compared the resulting sum of the race & ethnicity- and state-specific orphanhood incidence estimates to the previous state-specific orphanhood incidence estimates (Supplementary Fig. S18). For Alaska and Oklahoma, the sum of the race & ethnicity- and state-specific orphanhood incidence estimates was more than 20% below the state-specific orphanhood incidence estimates attributable to the leading caregiver loss cause-of-death. Supplementary Table S6 lists these states as ‘large discrepancy in estimates’. The table then reports, for all remaining states, 2021 race & ethnicity- and state-specific orphanhood prevalence rate estimates that are attributable to the leading caregiver loss cause-of-death.
S2 Sensitivity analyses
Sensitivity in mortality data
In the central analysis, we derived mortality counts by year, sex, age band, standardized race & ethnicity and caregiver loss causes-of-death from NCHS line-list mortality records from 1983 to 2021. Aggregate mortality counts are also available from CDC WONDER CDC [2023] from 1999 to 2021. We extracted data from the NCHS Vital Statistics portal because this allowed us to use the same data source across all years, bypass data suppression, and incorporate uncertainties in standardized race & ethnicity reporting. The CDC WONDER mortality counts allowed us to check our in-house data aggregations, if we obtain data from CDC WONDER at coarser population strata. For this purpose, we extracted annual death counts by sex, and age bands (15-19 years, …, 95-99 years, >100 years) from the CDC WONDER mortality portal CDC [2023] from 2000 to 2021. The overall mortality counts that we aggregated from the NCHS line list data were identical with the CDC WONDER mortality data without cause-of-death stratification. Secondly, we extracted annual death counts by sex, age bands (15-19 years, …, 95-99 years, >100 years) and cause-of-death (ICD-10 113 Selected Causes of Death) from the CDC WONDER mortality portal CDC [2023] from 2000 to 2021, without stratification by race & ethnicity. We then compared the mortality counts from the two data sources for each year and the leading caregiver loss causes-of-death (Supplementary Fig. S19). We found that across years, the maximum discrepancy in the CDC WONDER counts relative to the NCHS mortality counts was 0.0068% in women and 0.0057% in men, which reflected data suppression especially due to homicide excluding drug overdose in women, and overall indicated consistency of our data with that from CDC WONDER.
Sensitivity in live births data
In the central analysis, we derived live birth counts by year, sex, age band, and standardized race & ethnicity from NCHS line-list natality records from 1969 to 2021. Aggregate live birth counts are also available from CDC WONDER CDC WONDER [2022] from 1995 to 2021. We chose to extract data from the NCHS Vital Statistics portal because this allowed us to use the same data source across all years. The CDC WONDER live birth counts allowed us to check our in-house data aggregations. For this purpose, we extracted annual live birth counts from the CDC WONDER natality portal for women by age bands 15-19 years, …, 45-49 years from 2000 to 2021, and separately for men by age bands 15-19 years, …, 55+ years from 2016 to 2021 without stratification by race & ethnicity, which avoids data suppression. We then compared the live birth counts from the two data sources for each year and found that both data sets matched across all years.
Sensitivity in national-level orphanhood estimates to assumptions on historic fertility rates
In the central analysis, we assumed that male and female fertility rates by age and standardized race & ethnicity were in 1966-1989 constant and as in 1990 (Supplementary Fig. S7). This assumption was made because population size data stratified by detailed race categories (more than three categories) were not publicly available before 1990. Considering that age- and sex-specific live births data were available by race & ethnicity since 1978, we estimated in this sensitivity analysis the historic composition of population sizes by race & ethnicity in 1980-1989, and then updated the corresponding fertility rates and national-level orphanhood incidence and prevalence estimates. More specifically, we considered time trends in the proportion of each race & ethnicity in each 5-year age band and both sexes in the U.S. intercensal population size estimates in 1990 to 1995 (Supplementary Fig. S20A), and extrapolated these with linear models backwards in time to 1980-1989. We do not think that this estimation approach resulted in accurate estimates of population sizes, but rather that this approach conveys possible sensitivities in our orphanhood estimates. We identified minor sensitivities to national orphanhood prevalence estimates up to 2007, with no impact on incidence and prevalence estimates for recent years (Supplementary Fig. S20).
Sensitivity in national-level orphanhood estimates to potentially correlated fertility rates
In the central analysis, we assumed that population-level fertility rates were not correlated with population-level mortality rates. We considered in sensitivity analyses possible deviations from this assumption that might lead to upwards bias in our central estimates. For example, a person who died in year y may have experienced poor health in the preceding years y − 1, y − 2, … and in this case may also have been less likely to mother or father a child in years y, y − 1, y − 2, We modelled this possibility through four scenarios of dampened fertility rates close to individual death events. Specifically, we considered cumulative logistic adjustment factors that dampened fertility rates to either zero or half of the corresponding population-level, year-, age-, sex-, and standardized race & ethnicity-specific fertility rates in the year of death y. We also considered sensitivity analyses so that the onset of lower fertility rates preceded the year of death y by 1 or 3 years as shown in Supplementary Fig. S21A. We found that orphanhood incidence estimates were up to 8.4% lower than in the central analysis, and orphanhood prevalence estimates were up to 15% lower than in the central analysis (Supplementary Fig. S21B). It is also possible that our central orphanhood estimates might be biased downward. For example, women experiencing premature death are more likely economically disadvantaged and in turn are more likely to have had reduced access to contraceptives and corresponding higher fertility rates [Price et al., 2013], or earlier sexual debut resulting in higher cumulative fertility until death. It is also possible that overall, at population-level, there is no strong correlation between fertility and mortality rates especially as the lag between birth and death events is typically relatively large, and in the absence of data we have opted for this middle approach in the central analysis.
Sensitivity in national-level grandparent caregiver loss estimates to assumptions on the age of children experiencing loss of a grandparent caregiver
In the central analysis, we considered in Eq. (13) as a proxy to the age composition of children experiencing grandparent caregiver loss the age composition of children experiencing orphanhood. Calculations were done independently for each leading caregiver loss cause-of-death. We performed two analyses to characterise the sensitivity of this approach to national-level grandparent caregiver loss estimates. First, we repeated calculations using as proxy the age composition of children experiencing orphanhood across all causes of caregiver loss, i.e.
Second, we used data from the ‘Topical survey’ of the National Survey of Children’s Health (NSCH) [US Census Bureau, 2023a], filtered to respondents reporting to be grandparents and living with at least one child in the household. Respondents were asked to report on demographic characteristics including the age and race & ethnicity of one randomly selected child. We pooled these data across six years from 2016 to 2021 to characterise the age composition of children who live with grandparents due to small sample sizes in each survey after filtering (approximately 2, 500 grandparents per survey round). Supplementary Fig. S22A compares the age compositions used in the central analysis to those in the two sensitivity analyses. The comparison of the grandparent caregivers prevalence estimates are shown between the central estimate and two sensitivity analyses by age of child in Supplementary Fig. S22B. We found that grandparent caregiver incidence estimates deviated by up to ±0.028% and ±0.032% of the central estimate in the two sensitivity analyses due to the rounding numbers, respectively. Additionally, in the two sensitivity analyses, grandparent caregiver prevalence estimates were respectively 12.0% fewer and 13.6% higher than in the central analysis.
Sensitivity in national-level caregiver loss estimates to assumptions on historic numbers of grandparent caregivers
In the central analysis, we assumed that the proportion of adults aged 30+ years who live with their grandchildren of age 17 or under (denoted γy,s,r) was the same in the years y = 1983, …, 2009 as in 2010. This assumption was made because data were not available prior to 2010. To investigate the implications of this assumption, we considered longitudinal United Nations Population Division data on Households and Living Arrangments of Older Persons for the U.S. since 1960 United [2022]. These data indicate that the proportion of older persons who live with children or who are the primary caregivers of children has in the U.S. remained fairly constant since 1980 (Supplementary Fig. S23), and for this reason suggests that our assumptions on the historic number of grandparent caregivers is unlikely to have a substantive impact on caregiver loss estimates.
Acknowledgements
We thank the Global Reference Group for Children In Crisis, reviewers at the CDC and NCHS especially Dr. Robert Anderson for his helpful suggestions on interpreting and classifying disease groups and race groups using existing NCHS data. We also thank Prof. Chris Desmond for his comments on early versions of this work. We thank the Imperial College Research Computing Service (https://doi.org/10.14469/hpc/2232) for providing the computational resources to perform this study; and Zulip for sponsoring team communications through the Zulip Cloud Standard chat app. This study was supported by the Oak Foundation (to LC, LS); the Moderna Charitable Foundation (to OR); the World Health Organisation (to SF); the Engineering and Physical Sciences Research Council (EPSRC) through the EPSRC Centre for Doctoral Training in Modern Statistics and Statistical Machine Learning at Imperial College London and Oxford University (EP/S023151/1 to A. Gandy); the Imperial College London President’s PhD Scholarship fund to YC; Imperial College London Undergraduate Research Bursaries to LG and VKM; and London Mathematical Society Undergraduate Research Bursary to DW (URB-2023-86). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. The findings and conclusions in this report are those of the author(s) and do not necessarily represent the official position of the Centers for Disease Control and Prevention.
Footnotes
↵* Joint first authors
↵§ Joint senior authors
Disclaimer: The findings in this report are those of the authors. They do not necessarily represent the official position of the US Centers for Disease Control and Prevention. For the purpose of open access, the authors have applied a ‘Creative Commons Attribution’ (CC-BY) license to any Author Accepted Manuscript version arising.
“The death of a parent is one of the last taboos in public health. The acute and chronic pain of parental loss and the resulting tragedy for the children … left behind are almost too distressing to contemplate.” The Editorial, The Lancet Public Health, August 2022.