Abstract
The goals of this study were to identify patterns of polysubstance use and their associations with stressful life events among U.S. late middle-aged and older adults and examine whether gender moderates these associations. Adults aged 50 and older (N = 14,738) from the National Epidemiological Survey on Alcohol and Related Conditions-III were included. Latent class analysis was conducted to identify patterns of polysubstance use. Weighted multinomial logistic regression was estimated with a generalized structural equation model. Three different polysubstance use patterns (non-users/low substance users; cannabis and excessive alcohol users; painkiller and sedative/tranquilizer misusers) were identified. Higher levels of stressful life events were associated with patterns of polysubstance use. Gender moderated the association between stressful life events and co-misusing painkillers and sedatives/tranquilizers (p < 0.05). Substance use prevention efforts should consider aging adults’ patterns of polysubstance use and associated stressful life events when designing and implementing gender-specific polysubstance use prevention interventions.
Keywords
• This is the first study that used latent class analysis to identify patterns of polysubstance use among U.S. late middle-aged and older adults on a national scale. • Stressful life events were associated with an increased likelihood of polysubstance use. • Gender moderated the association between stressful life events and polysubstance use (i.e., co-misusing painkillers and sedatives/tranquilizers).
• Given that the aging population is growing rapidly, substance use prevention efforts should consider aging adults’ patterns of polysubstance use. • Substance use prevention efforts for preventing polysubstance use, such as early screening and health education, should target aging adults who experience stressful life events as an at-risk population. • Intervention efforts to reduce polysubstance use should address stressful life events, implement gender-specific interventions, and customize approaches based on the patterns of polysubstance use.What this paper adds
Applications of study findings
Introduction
Numerous studies conducted over recent decades have shown that one of the greatest public health challenges in the United States (U.S.) is substance use. The increase of substance use among the aging population has been a notable public health concern. For example, one study projected that there will be approximately 72 million older adults misusing substances by 2050 (Youdin, 2019). Substance use has been shown to be associated with drug overdose, attempted suicide, and increased mortality (Schulte & Hser, 2013). Moreover, substance use causes an additional and serious economic burden to families, communities, and governments (Florence et al., 2016). Recently, studies have noted that polysubstance use, referring to using two or more substances concurrently, is an emerging public health concern in the U.S. (Ogbu et al., 2015). For instance, one study found that opioid use commonly overlapped with other substances such as alcohol and sedatives, which has been shown to cause adverse events including overdose deaths (Compton et al., 2021). The additive effects vary by different patterns of polysubstance use. For example, one study reviewed laboratory studies of alcohol and tobacco interactions, suggesting that co-use of alcohol and tobacco was associated with craving and thus increased alcohol and tobacco consumption (Verplaetse & McKee, 2017). In addition, one review study concluded that co-use of opioids and sedatives could result in respiratory depression, which in turn leads to overdose deaths (Hassamal et al., 2016). Given that different patterns of polysubstance use pertain to different additive effects of substances and are associated with various adverse health effects, a better understanding of patterns of polysubstance use is critical for designing tailored prevention efforts.
Polysubstance use has been assessed by summing up multiple substance use indicators such as alcohol and cannabis use statuses (Fendrich et al., 2021). However, such measurement strategies for polysubstance use neither distinguish polysubstance use patterns nor consider the effect of interactions between substances. Moreover, previous polysubstance studies often focused on a single pattern of polysubstance use such as co-misuse of opioids and tranquilizers or co-use of illicit substances and cannabis (Jones et al., 2012; Olthuis et al., 2013). In fact, Timko et al. (2018) conducted a latent class analysis to investigate the patterns of polysubstance use among stimulant users. Yet, the sample characteristic—the mean age of the sample was 32.5 years—limits the study’s generalizability to the aging population. Additionally, the way in which sociodemographics, physical characteristics, and psychological risk factors are associated with each pattern of polysubstance use are not fully explored. A more thorough understanding of polysubstance use patterns and their associations with sociodemographic characteristics and physical and psychological risk factors can better inform interventions to prevent and reduce polysubstance use. It is particularly important to study these relations among older adults given the demographic shifts in our population resulting in increasing numbers of individual in older age groups.
Aging is associated with social, emotional, and physical changes. Age-related changes as well as associated increasing chronic illness may increase vulnerability to substance use and misuse (Moos et al., 2004). In particular, older adults are at great risk of the adverse effects from substance use because of polypharmacy (i.e., simultaneous use of multiple medications) and pharmacokinetic changes (e.g., decreased renal elimination of drugs, decreased hepatic metabolism) (Kalapatapu & Sullivan, 2010). Further, aging adults may be less likely to perceive substance use as a problematic behavior (Wu & Blazer, 2011). In the U.S., the sizeable aging boomer cohort started turning 65 years old in 2011. By 2050, over 85 million, approximately 22%, Americans will be over age 65 (United States Census Bureau, 2018). Given that the projected increase in the aging population and their substance use and misuse, a better understanding of substance use behaviors and associated risk factors among the aging cohort in the U.S. is imperative. However, previous studies have focused primarily on single substance use among the aging cohort or on polysubstance use among the general adult population. One study suggested that using a lifespan perspective approach to investigate substance use prevention is needed (Schulte & Hser, 2013), implying that it is crucial to break down the study population into different age cohorts to investigate (poly)substance use among middle-aged and older adults. Thus, with a better understanding of polysubstance use among aging cohorts, substance use prevention and intervention programs can be better tailored by addressing multiple substance use behaviors.
A stressful life event is defined as an event that results in psychological and physiological stress responses for an individual (Myers et al., 2014). A growing body of evidence has shown that stressful life events are connected to substance use (Verplaetse et al., 2018) and problem drinking behaviors (Brennan & Moos, 1996). Supported by empirical evidence, general strain theory (Agnew, 2001) posits that individuals who experience strains (e.g., stressful life events) may be upset or produce negative emotions, and substance use serves as a way to manage negative emotions (Akins et al., 2010). Yet, an important limitation of existing studies that explored stressful life events and associated substance use is that they focused on single substance use rather than polysubstance use. In addition, previous studies suggest that gender differences exist in the associations between stressful life events and substance use (Sacco et al., 2014). For instance, compared with females, males are more likely to cope with strains through delinquent behaviors such as substance use (Brezina, 2017). Therefore, building on these previous findings, it is important to further investigate whether gender moderates the associations between stressful life events and polysubstance use. To the best of our knowledge, no studies have investigated the moderating role of gender on the associations between stressful life events and polysubstance use patterns, in particularly among middle-aged and older adults. Understanding gender differences can inform gender-specific intervention design in the context of stressful life events and polysubstance use for the aging population.
Taken together, this study aimed to fill literature gaps by investigating polysubstance use among adults aged 50 and older. The purpose of this study was three-fold: (a) to identify patterns of polysubstance users; and (b) to test the associations between past-year stressful life events and polysubstance use patterns; and (c) to investigate the moderating effects of gender on the relationship between stressful life events and polysubstance use patterns in a nationally representative sample of U.S. adults aged 50 and older. This study hypothesized that: (a) adults aged 50 and older who experienced stressful life events would report more polysubstance use; and (b) male adults aged 50 and older who experienced high levels of stressful life events would have more polysubstance use than their female counterparts.
Methods
Data and Study Sample
Data were extracted from Wave 3 of the National Epidemiological Survey on Alcohol and Related Conditions (NESARC-III), a cross-sectional survey of a nationally representative sample of U.S. non-institutionalized civilians. It was conducted by the National Institute on Alcohol Abuse and Alcoholism (NIAAA) in 2012-2013. NESARC-III used multistage probability sampling to randomly select U.S. adults aged 18 and older (Grant et al., 2015). Of the original sample from NESARC-III, a total of 14,738 adults aged 50 and older were included in the analyses. This study was verified as a non-human subjects study by the Institutional Review Board of the author’s institution due to the deidentified data of NESARC-III (protocol number: 1703680184).
Measurement
Outcome Variable
Polysubstance use pattern was constructed as the outcome variable using latent class analysis (LCA). This study included six, binary (yes/no), indicators of substance use in the past year. Substances that have a euphoriant effect and are commonly used/misused were included, including excessive alcohol use (i.e., binge drinking), cannabis use, cocaine use, and prescriptions drug misuse (painkillers, sedatives/tranquilizers, stimulants) (Dackis & O'Brien, 2001; Mukherjee et al., 2008; Weaver, 2015). Literature shows that excessive alcohol use, illicit drug use (such as cocaine), and prescriptions drug misuse increase the risk for injury and overdose death (Seth et al., 2018). Also, this study included cannabis use as more and more states have implemented medical and recreational cannabis legalization, which has been linked to increased cannabis use among older adults (Han & Palamar, 2020). NESARC-III respondents were asked about their past-year experiences of using these substances. For example, during the interview, they were asked “Have you ever used sedatives or tranquilizers non-medically? For example, barbs, downers, Ambien, Lunesta, phenobarbital, pentobarbital, Halcion, Tuinal, Nembutal, Seconal, Librium, Valium, Xanax, benzodiazepines, tranks, Ativa.” Excessive alcohol use was measured as having 5 or more drinks (for men) or 4 or more drinks (for women) in a period of two hours or less. Notably, this study did not include tobacco use because recent literature suggests that nicotine produces little euphoriant effect (Dar et al., 2007; Wei et al., 2018).
LCA was conducted to identify the patterns of polysubstance use in the study sample. LCA can be used to identify distinct subgroups among study population based on observed categorical variables, classifying individuals with similar behaviors or characteristics into latent classes (Lanza & Rhoades, 2013). A key assumption of LCA is local independence, requiring individuals’ responses to items to be independent. We first estimated the 1-class model, serving as a comparison class with the other classes. We then increased the number of classes until the model no longer converged or had meaningful applications (Lanza et al., 2020). Optimal model selection was based on fit statistics including the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC), as well as model stability and interpretability (Nylund et al., 2007). A lower value (closer to zero) of AIC and BIC indicates a better balance model fit and model parsimony. We finally determined the number of classes and the membership of each class among all participants. This variable was a nominal variable.
Main Predictor
Number of stressful life events was the main predictor, which was assessed based on sixteen questions regarding past-year stressful experiences, including moving or someone new moving in, being laid off, being unemployed and looking for new job, having interpersonal conflict at workplace, starting a new job, getting divorce, experiencing death of a family member or close friend, having financial problems (unable to repay debt, declared bankruptcy), having interpersonal conflict with friends or relatives, being a victim of a crime (property has been stolen or damaged), having a family member or close friend who was the victim of a crime (being physically assaulted or attacked), having trouble with the police or the law, having family member or close friend who had trouble with the police or the law, and being homeless (Verplaetse et al., 2018). Each question was coded as “1” if responded yes. A summed total number of experienced stressful events (range 0 to 16) was used as continuous predictor variable. This scale has been widely used in previous studies and has been shown to be reliable (Lehavot et al., 2018; Vásquez et al., 2019). In this study, Cronbach’s alpha was .82, indicating high internal consistency of this scale.
Moderator
This study tested the moderating effect of gender (male = 0; female = 1) on the associations between stressful life events and polysubstance use patterns.
Covariates
Sociodemographic variables were included as covariates, including age (50–64, 65 and above), gender (male or female), race/ethnicity (non-Hispanic White, non-Hispanic Black, Asian/others, Hispanic), education (less than high school, high school, some college, college, or graduate school), marital status (married or not), unemployment (yes or no), and urbanicity (urban or rural).
Analytical Approach
Descriptive Analysis
This study first computed descriptive statistics for each variable. Unweighted frequencies and weighted percentages were reported for all categorical variables. Weighted means and standard errors were reported for continuous variables.
Generalized Structural Equation Modeling
Generalized structural equation modeling (GSEM) was used to conduct latent class analysis to identify subgroups of participants with distinct polysubstance use behaviors. After identifying the best-fitting model and the number of classes (as described in the measurement section), a multinomial logistic regression was estimated using GSEM to examine the associations between stressful life events and polysubstance use, controlling for all covariates. In the multinomial logistic regression, Class 2 (non-users/low substance users) was used as a reference group. Additionally, an interaction term between gender and stressful life events was included in the analyses to test the moderating effect of gender, which was further validated by an invariance test as a post-hoc test under GSEM. To produce nationally generalizable estimates, all analyses were weighed using the NESARC-III survey weights and sampling scheme. All statistical analyses were conducted using Stata 15.1. GSEM was fitted using the svy: gsem command in Stata.
Results
Descriptive Statistics of the Study Sample.
Noted: Weighted N = 102,465,952; Unweighted N = 14,738.
Data source: National Epidemiological Survey on Alcohol and Related Conditions (NESARC-III).
SE = standard error.
Weighted mean and standard error are computed for continuous variables.
Latent Class Model Fit Indices.
Noted: AIC = Akaike information criterion; BIC = Bayesian information criterion.
Bolded AIC and BIC indicate the best model fit.

Latent class analysis of polysubstance use patterns.
Weighted Multinomial Logistic Regression Estimates of Polysubstance Use by Generalized Structural Equation Modeling.
Noted: Weighted N = 102,465,952; Unweighted N = 14,738.
Data source: National Epidemiological Survey on Alcohol and Related Conditions-III (NESARC-III).
AOR: adjusted odds ratio; CI: confidence interval.
*p < .05, **p < .01, ***p < .001.
Regarding the moderating effect of gender on the associations between stressful life events and polysubstance use patterns in this sample of older adults, the results indicate that when experiencing high levels of stressful life events, males have an even higher likelihood to co-misuse painkillers and sedatives/tranquilizers than females (AOR = .84 [female], 95% CI: .72-.99, p < .05). The invariance test further confirmed the significant difference between male and female participants on this association (χ2 = 99.31, d.f. = 58, p < .05). However, the moderating effect of gender was not observed on the association between stressful life events and co-use cannabis and excessive alcohol.
Discussion
This study contributes to the literature by identifying patterns of polysubstance use and investigating associations between stressful life events and polysubstance use patterns among U.S. adults aged 50 and older, using data from a nationally representative study. Unlike previous studies, this study considered multiple substance use behaviors concurrently to identify patterns of polysubstance users (i.e., cannabis/excessive alcohol, painkillers/sedatives or tranquilizers). Moreover, this study extended previous studies of stressful life events and their associations with substance use. We found that past-year stressful life events were associated with co-using cannabis and excessive alcohol as well as co-misusing painkillers/sedatives or tranquilizers, and that gender moderated the later association.
Our results showed that 3.8% of U.S. adults aged 50 and older reported polysubstance use in the past year, including 2.4% co-using cannabis and excessive alcohol and 1.4% co-misusing painkillers and sedatives/tranquilizers. The most common polysubstance use pattern among adults aged 50 and older was co-using cannabis and excessive alcohol. Such a pattern is consistent with previous finding that middle-aged and older cannabis users were more likely to report excessive alcohol use (Salas-Wright et al., 2017). Previous studies found a higher blood tetrahydrocannabinol (THC) level in individuals co-using cannabis and alcohol compared to those who used only cannabis, indicating an increased risk for risky behaviors such as impaired driving among those co-using cannabis and alcohol (Hartman et al., 2015). In addition, middle-aged and older adults may experience age-related cognitive decline as their age increases (Salthouse, 2009), thus, co-use of cannabis and alcohol might exacerbate the problem resulting in increased incidence of adverse outcomes including dementia, being scammed, and traffic-related injuries and deaths. Substance use prevention programs and other interventions should enhance awareness of the consequences on co-using cannabis and alcohol.
This study identified co-misusing prescription drugs (i.e., painkillers and sedatives/tranquilizers) as the second common polysubstance use pattern among middle-aged and older adults. This is also consistent with prior findings that painkiller misusers are more likely to misuse sedatives/tranquilizers (Jones et al., 2012). Previous findings suggest that the main reason for misusing painkillers is to manage chronic pain (Lipari et al., 2017) as chronic pain is highly prevalent among older adults and often occurs among those who have comorbidities (Reid et al., 2015). In addition, previous studies indicated that two in three adults with chronic pain suffered from insomnia, resulting in higher likelihood of co-misusing painkillers and sedatives/tranquilizers among middle-aged and older adult population (Miller et al., 2018; Rubio et al., 2016). In fact, to mitigate the opioid epidemic and reduce overdose deaths due to concurrent use of opioids and sedatives, the Centers for Disease Control and Prevention (CDC) published guidelines for prescribing opioid analgesics to treat chronic pain, cautioning that concurrent use of opioids and sedatives should be avoided. The CDC also suggested that non-pharmacological therapy combined with non-opioid pharmacological therapy are preferred for adults to manage chronic pain (Dowell et al., 2016). In summary, although only 3.8% of the study sample were past-year polysubstance users, whereas 96.2% were non-users or used only a small amount of substances in past year, one study suggested that 1 in 3 overdose deaths was associated with polysubstance use (Kariisa et al., 2019). Thus, incorporating polysubstance use aspects into substance use intervention and control programs may reduce drug overdose and subsequent mortality. Moreover, because aging is associated with physical changes and increased odds of multi-morbidity, interventions to reduce polysubstance use should be better tailored to the characteristics of middle-aged and older adults.
A key finding from our study was that higher levels of stressful life events were associated with an increased likelihood of co-using cannabis and excessive alcohol as well as co-misusing painkillers and sedatives/tranquilizers among adults aged 50 and older, even after controlling for multiple sociodemographic variables. These findings suggest that substance use prevention efforts for late middle-aged and older adults, such as early screening and health education, should target those who experience stressful life events as an at-risk population. Specifically, around one-third of our study participants (35.24%) experienced the death of family members or a close friend in the past year. A previous study suggested that the death of a family member is associated with a higher likelihood of substance use (Glass et al., 1995). Thus, special attention to substance use prevention efforts should be given to this group. In addition, evidence-based programs such as the Stress Management and Resiliency Training (SMART) Program have been shown to be effective in managing psychological stress of late-life spousal loss and in enhancing resilience (Bui et al., 2018; Chesak et al., 2019). Such programs may have a secondary effect on preventing polysubstance use.
Our study results also showed that gender moderated the association between stressful life events and co-misusing painkillers and sedatives/tranquilizers. In the current study sample, males who experienced high levels of stressful life events were more likely to co-misuse painkillers and sedatives/tranquilizers than females experiencing the same level of stressful life events. This result is in line with general strain theory and previous findings that males are more likely to cope with strains by using a substance (Rubio et al., 2016). Therefore, prevention and intervention efforts that pay special attention to males who experience stressful life events to enhance their coping skills should be considered. Moreover, general strain theory suggests that coping skills mediate the associations between strains and substance use. Thus, future studies should incorporate a moderated mediation approach to examine coping as a pathway through which stressful life events and polysubstance use are connected, and how gender may moderate this pathway.
The current study also has secondary findings. Adults aged 65 and older were less likely to engage in polysubstance use compared with adults aged between 50 and 64. In particular, the majority of our study participants aged 50 to 64 were baby boomers who had greater exposure to substances in their youth stage than earlier generations (Specht et al., 2021). The result implies that targeting baby boomers in substance use prevention and reduction programs should be emphasized. This study found that married participants were at a lower risk for polysubstance use, which is consistent with previous studies (e.g., Salvatore et al., 2020). This study also, aligned with prior studies (e.g., Moody et al., 2016), found that unemployment is a risk factor of polysubstance use. However, the causal relationship between unemployment and polysubstance use is still unclear. Future study is needed before drawing conclusions of causality such as whether some people consider polysubstance use as a strategy to cope with unemployment.
There are limitations of the present study. First, due to the cross-sectional nature of NESARC-III, this study cannot determine the temporal sequence between stressful life events and polysubstance use. Future studies should employ a longitudinal study design to examine causal relationships between stressful life events over the life course and polysubstance use. Second, the information about stressful life events and past-year substance use/misuse was self-reported by respondents, thus, the information is subject to recall and response bias. Further, as described earlier, the general strain theory suggests that coping may serve as a critical mediator of the associations between stressful life events and substance use, which implies that coping skills may be potential pathways through which stressful life events are connected with polysubstance use. However, coping was not assessed as part of NESARC-III and the data are cross-sectional in nature, preventing consideration of mediation.
Despite these limitations, to our knowledge this is the first study that used LCA to identify the patterns of polysubstance use among adults aged 50 and older in the U.S. and that examined the associations between past-year stressful life events and polysubstance use. The LCA reveals the three most common polysubstance use patterns among aging adults on a national scale, making a unique contribution to present knowledge of polysubstance use among middle-aged and older adults. Given that the aging population is growing rapidly, our study findings could inform substance use intervention programs to incorporate aging adults’ polysubstance use behaviors. Substance intervention efforts to reduce polysubstance use should address stressful life events, implement gender-specific interventions, and customize approaches based on the patterns of polysubstance use. Through well-designed, tailored programs targeting at-risk subgroups, polysubstance use among middle-aged and older adults in the U.S. could be reduced.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Ethical Approval
This study was verified as a non-human subjects study by the Institutional Review Board of Indiana University due to the deidentified data of NESARC-III (Protocol Number: 1703680184).
