Introduction
The COVID pandemic poses a significant challenge to mental health, with reports of increasing levels of anxiety, depressive symptoms, and psychological distress in the general population in multiple countries spanning across Asia, Europe, and the USA (2020; Pfefferbaum & North, Reference Pfefferbaum and North2020; Vindegaard & Benros, Reference Vindegaard and Benros2020). Several agencies and expert authorities have issued recommendations for promoting good mental health during the pandemic, including the Centers for Disease Control and World Health Organization (CDC, 2020; WHO, 2020). There are three central tenets in the recommendations. Firstly, individuals are encouraged to maintain social connections through virtual (phone, online) rather than physical contact (Brooks et al., Reference Brooks, Webster, Smith, Woodland, Wessely, Greenberg and Rubin2020; CDC, 2020; WHO, 2020). Secondly, limiting the use of substances has been encouraged during the pandemic (Brooks et al., Reference Brooks, Webster, Smith, Woodland, Wessely, Greenberg and Rubin2020; CDC, 2020), informed by evidence that isolation and stress are important risk factors for substance misuse, including alcohol use (Clay & Parker, Reference Clay and Parker2020). Finally, moderating exposure to COVID-related news has also been recommended (Amsalem, Dixon, & Neria, Reference Amsalem, Dixon and Neria2020; CDC, 2020). This recommendation is supported by reports of associations between news/media consumption about COVID-19 and poor mental health, including anxiety and depressive symptoms (Gao et al., Reference Gao, Zheng, Jia, Chen, Mao, Chen and Dai2020; Ni et al., Reference Ni, Yang, Leung, Li, Yao, Wang and Liao2020; Olagoke, Olagoke, & Hughes, Reference Olagoke, Olagoke and Hughes2020). The efficacy of these recommendations rests on the assumption that these behaviors impact subsequent mental health, rather than mental health impacting the behavior. In this study, we test three hypotheses that underlie these current public health recommendations. Specifically, in line with current recommendations, we hypothesized that better mental health will be promoted by (1) maintaining virtual contact with others; (2) limiting the use of substances (alcohol and illicit drugs); and (3) limiting the amount of news that is consumed. The three hypotheses are directional in that all predict that behaviors will impact on subsequent mental health, rather than the opposite direction of effects.
To test the hypotheses, we apply random intercept cross-lagged panel modeling (RI-CLPM) to daily ‘diary’ survey data acquired in June 2020, in order to parse the effects of behavior on mental health and vice versa. RI-CLPM disentangles dynamic, within-person processes relating behaviors to mental health from more stable individual differences (Hamaker, Kuiper, & Grasman, Reference Hamaker, Kuiper and Grasman2015). The dynamic processes relating behavior and mental health within persons are delineated as lagged effects, both autoregressive and cross-lagged. The within-person cross-lagged measures indicate the degree to which a daily deviation from usual levels of one variable is related to change the following day in another. In essence, the cross-lagged effect maps how ‘state-like’ features of behavior and mental health exert a mutual influence over time. The random intercepts that are included in these models capture the more ‘trait-like’, stable aspects of behavior and mental health. Correlations between the random intercepts thus represent associations between the more stable facets of behavioral and psychological measures. In summary, using daily diary survey data, we determine if daily variation in key pandemic-related behaviors – keeping social distance from others, maintaining virtual contact, and watching news/media – has an effect on the subsequent mental health; a direction of effects that would be consistent with current recommendations.
Methods
Individuals were recruited via Qualtrics' actively managed, double-opt-in market research panels. We aimed to oversample African American, Hispanic, and low-income populations, given evidence of a disparate impact of COVID on these groups (Fitzpatrick, Harris, & Drawve, Reference Fitzpatrick, Harris and Drawve2020; Sneed, Reference Sneed2020).
We thus created six categories to guide recruitment: African American and low income; African American and high income; Hispanic and low income, Hispanic and high income; White, non-Hispanic, and low income; and White, non-Hispanic, and high income. The sampling strategy thus aimed to recruit equal numbers from these six groups by targeting states that included counties with demographics that met this goal. As we were interested in the effects of social distancing, we further restricted recruitment to states with active ‘stay-at-home’ or ‘safer at home’ orders projected to be ongoing at the time daily diary recruitment aimed to be opened (May 15). This strategy led to the inclusion of the following states: Alabama, California, Oregon, Kentucky, Louisiana, Maine, Maryland, New Jersey, New York, Pennsylvania, Tennessee, Virginia, Washington, and Washington DC (online Supplementary Fig. S1a). After the study collection period ended, we evaluated infection rate trends during the daily diary data collection period (21 May through 16 June). Two trends emerged: states with rising and states with falling infections (online Supplementary Fig. S1a). This was used as a covariate in RICPLM models (see online Supplementary Table S10). Participants were compensated with $15 for completing the baseline survey and $1 for every daily survey completed.
Adults with online/smartphone access completed a baseline and daily survey, sent to their email/smartphone every evening. Written online consent was obtained from all participants. The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. The study was determined to be exempt from IRB review by the University of Pennsylvania.
We set a target recruitment of 1750. Beginning 5 May, Qualtrics invited approximately 94 000 individuals through their panels to participate in the study. Recruitment was closed on 26 May following successful recruitment of 1719 participants and to facilitate transition to the daily diary. Of these 1719 individuals, 1213 joined the daily diary (open to all who completed the baseline survey). Of those, 1148 completed daily surveys and met our inclusion criteria for entry in the cross-lagged analyses (66.8% of the total who completed the baseline survey). Those who completed the daily diary questionnaire (N = 1148) tended to be older, white/non-Hispanic, female, and have higher income than those who did not complete the daily survey (see online Supplementary Table S1).
The criteria for inclusion were (1) completion of at least one response in the first 3 days of the 15 days of data collection, (2) completion of at least three responses per week of responses, (3) completion of at least half of surveys within 1 day of receipt. The proportion of missing data among the 1148 who were included in the cross-lagged analyses was low. Our final data set included 23 028 observations out of a possible 24 108 (21 responses per 1148 participants), giving a missingness rate of 4.5%. Imputation of the missing values was performed using maximum likelihood estimation within the lavaan function for each pairwise model.
The daily diary questionnaire
The survey contained 11 questions on mental health completed by participants every evening – full survey in online Supplementary material. The 11 questions pertaining to general mental health were adapted from the Patient Health Questionnaire-9 (PHQ-9) (Kroenke, Spitzer, & Williams, Reference Kroenke, Spitzer and Williams2001) and General Anxiety Disorder-7 (GAD-7) (Spitzer, Kroenke, Williams, & Löwe, Reference Spitzer, Kroenke, Williams and Löwe2006). Items on COVID-related worrying were adaptations of the Disease Worry and Risk Perception scales, originally used in cancer research (Lerman et al., Reference Lerman, Trock, Rimer, Jepson, Brody and Boyce1991) – see online Supplementary Table S2 for breakdown of mental health questions. Longitudinal confirmatory factor analyses were run in lavaan, R, using the function measEq.syntax in the package semTools. Fitted models included a one-factor solution containing all 11 questions, a two-factor solution splitting mental health from COVID-related worry, a three-factor solution containing questions split by source (GAD-7, PHQ-9, DWRPC), and a three-factor solution split by positive valence mental health, negative valence mental health, and COVID-related worry. This final solution provided the best fit in longitudinal confirmatory factor analysis across all levels of invariance included in models [configural (CFI: 0.999, RMSEA: 0.025), threshold (CFI: 0.999, RMSEA: 0.025), metric (CFI: 0.999, RMSEA: 0.025), scalar (CFI: 0.999, RMSEA: 0.025), and strict (CFI: 0.999, RMSEA: 0.028)] – see online Supplementary Table S3 for full results. The extracted factors were defined as psychological struggling (levels of sadness, anxiety, loneliness, and stress), psychological striving (levels of pleasure, happiness, optimism), and COVID-related worry (concerns of infection of self or others, and its impact on daily functioning).
Daily news/media consumption was determined from three items: ‘How much time did you spend reading or listening to the news today?’, ‘Overall, how much time did you spend browsing social media today (e.g. Facebook, Twitter, TikTok, Instagram)?’, and ‘How much did you read or hear about the coronavirus pandemic today?’ This variable also showed a robust factor solution (configural CFI > 0.999, RMSEA 0.007; threshold CFI > 0.999, RMSEA < 0.0001; metric CFI > 0.999, RMSEA 0.01; scalar CFI > 0.999, RMSEA 0.025; strict CFI 0.999, RMSEA 0.035) when considering the repeated observations over time.
The daily amount of physical contact with others was determined for those within and outside of the home. Contact within the home was coded as the sum of people interacted with within the home environment (within 20 feet) and guests to the home environment (capping the maximum number of contacts at five). Contact outside the home was coded as the number of settings that involved close physical contact (within six feet) of others, with responses ranging from zero to seven settings. The settings were shopping/running errands, working, meeting a friend solo, meeting friends in a group setting, walking/leisurely activity alone, walking/leisurely activity with others, and attending organized meetings. Total physical contact was the sum of the within-home and outside-home contacts. Virtual contact was the summed amount of time participants reported interacting with family (zero to four, ranging from no contact to more than 2 h of contact) and friends (zero to four, ranging from no contact to more than 2 h of contact) each day.
Alcohol use coded as the number of daily alcoholic drinks consumed, from a minimum of zero to a maximum of five (for five or more alcoholic drinks). Substance use was coded as the number of substances used each day, including tobacco products, e.g. cigarettes and vapes, marijuana, opiates including heroin and other narcotics, and other drugs including amphetamines, cocaine, ecstasy, and hallucinogens. Responses ranged from zero (no substance use) to four (all four substances used).
Analysis
Random intercept cross-lagged panel analysis (implemented in Lavaan R) separated dynamic within-individual processes (measured as lagged effects between daily changes in behavior and mental health) from more stable individual differences in the variables (measured as the covariance between random intercepts – online Supplementary Fig. S1b, Supplementary Methods) (Hamaker et al., Reference Hamaker, Kuiper and Grasman2015; Mulder & Hamaker, Reference Mulder and Hamaker2020). The daily items are modeled as within-person latent factors and a latent intercept factor. The random intercepts are specified by creating a latent variable with the repeated measures as its indicators and fixing all the factor loadings to one. The covariance between the random intercepts indicates how stable individual differences in one measure are associated with another. The within-person components are specified as a latent variable for each measurement and constraining measurement error variances to zero. The model also specifies that within-person components at the first observation and the within-person residuals at all subsequent occasions are correlated within each occasion. For the between-person components, we specify that the random intercepts can be correlated. See online Supplementary Fig. S2 for a diagram of an example model.
The main analysis held the cross-lagged and auto-regressive effects constant over the 14 daily lags in the study. We used this approach as we did not expect relationships between mental health and behavior to vary substantially over this short time. The approach was also appropriate given the regular daily intervals between observations and favored on the grounds of parsimony as it estimates four parameters (two cross-lagged and two autoregresssive) compared with 56 parameters in the fully unconstrained model (four parameters for each of the 14 lags estimated in our 15-day study). We compared a constrained model against a fully unconstrained model (daily lag values allowed to vary). In line with recommendations, we used two indices of model fits: the difference between models in the comparative fit index (CFI), and differences in the root mean square error approximation (RMSEA) (Chen, Reference Chen2007; Cheung & Rensvold, Reference Cheung and Rensvold2002). We adopted the recommended critical values of ΔCFI > −0.01 and ΔRMSEA > 0.01 as indicating that the constrained model provided a worse fit to the data than the unconstrained model. This approach is preferred to the examination of differences in χ2 of model fits, which has several limitations including being highly sensitive to sample size (Brannick, Reference Brannick1995). As can be seen from online Supplementary Table S4, the constrained and unconstrained models showed consistently acceptable fit using the CFI and RMSEA metrics. When applying the constraint, the global fit indices did not change substantially, and thus we adopt the theoretically preferred, more parsimonious constrained model in the main report.
As noted in the introduction, our goal was to test the directional hypotheses that emerge from current public health recommendations: better mental health will be promoted by (1) maintaining virtual contact with others; (2) limiting the use of substances (alcohol and illicit drugs); and (3) limiting the amount of news that is consumed. By their nature, the cross-lagged analyses we used test other directional effects (both autoregressive and the effects of mental health on later behavior). For completeness, we also tested cross-lagged effects within domains (does one behavior impact another, etc.). In total, we thus tested the relationships between 28 pairs of variables, each with two cross-lagged and two auto-regressive effects, making 112 effects. Of those 112 effects, three pertained to the primary hypotheses, and the remainder parse all possible pairwise combinations of cross-lagged and autoregressive effects. Reported p values were adjusted using the false discovery rate procedure to account for multiple testing (Benjamini & Hochberg, Reference Benjamini and Hochberg1995).
Results
The final cohort comprised 1148 adults who met the criteria for adequate data completion (41 did not complete data satisfactorily – Methods). We also removed the 24 subjects who tested positive for SARS-CoV2 during the study (N = 24) – participant flow is shown in Fig. 1. We aimed to have equal representation of African American (non-Hispanic), White non-Hispanic and Hispanic groups, but the final cohort contained fewer Hispanic members than anticipated [N = 327 (28.5%); White, non-Hispanic N = 463 (40.3%); African American, non-Hispanic N = 358 (31.2%)]. The gender breakdown of the final cohort was N = 658 (57.2%) female, N = 485 (42.2%) male, and N = 5 (0.04%) non-binary (did not identify as male or female). The average age was 40.7 years (s.d. 11.5). We sampled equally from states that had increasing and decreasing rates of confirmed SARS-CoV2 infection during the study – online Supplementary Fig. S1a. In total, 760 (66.2%) subjects entered the study on 21 May 2020 and 388 (33.8%) entered 6 days later, and we confine these analyses to the overlapping 15 days.
Our primary focus was on the cross-lagged effects of behaviors (daily physical and virtual contact with others, substance use and consumption of news/media) on mental health. Daily means and standard deviations for these variables are located in online Supplementary Table S5. Variables were not centered. Using RI-CLPM, we found that daily fluctuations in news/media consumption showed a significant cross-lagged effect on subsequent COVID-related worrying [standardized estimate = 0.034 (95% CI 0.018–0.049), adjusted p = 0.00005] – see Fig. 2. Thus, for a given individual, a higher than usual amount of news consumption on a given day is followed the next day by an increase in COVID-related worries. The inverse cross-lagged effect of COVID-related worries on subsequent media consumption was also present [estimate 0.03 (0.012–0.048), adjusted p = 0.0017]. Full autoregressive and cross-lagged effects for all possible combinations are given in Table 1, with significant cross-lagged effects following FDR correction illustrated in Supplementary Fig. S3.
Effects are standardized with 95% confidence intervals and FDR-adjusted p values. Significant cross-lagged effects (following adjustment are shown in bold).
Daily fluctuations in news consumption also exacerbated psychological struggling [estimate 0.064 (0.03–0.098), adjusted p = 0.0005]. By contrast, striving had a negative cross-lagged effect on subsequent news consumption [−0.016 (−0.0027 to −0.03), adjusted p = 0.03], meaning above-average striving decreased the amount of news consumed the following day. A path from news to striving was not present.
Daily increases in substance use, but not alcohol consumption, were associated with subsequent decreases in psychological striving [−0.17 (−0.03 to −0.31), adjusted p = 0.03] and COVID-related worrying significantly decreased subsequent drug use [−0.0041 (−0.0015 to −0.0067), adjusted p = 0.003].
Neither physical nor virtual contact with others had a significant cross-lagged effect on mental health. Findings held when considering physical contact either within or outside the home: no significant lags with any mental health variable emerged when inside and outside contact were treated as distinct variables – online Supplementary Table S6. Additionally, all results were non-significant when virtual contact with work colleagues was added to virtual contact with close family and friends.
We found significant cross-lagged effects between the mental health items: for example, striving significantly decreased subsequent COVID-related worries [−0.02 (−0.033 to −0.007), adjusted p = 0.003]. By contrast, psychological struggling increased subsequent COVID-related worries [0.012 (0.0044–0.020), adjusted p = 0.003], and vice versa [0.071 (0.034–0.11), adjusted p = 0.0003]. We also found significant cross-lagged effects between the behavioral items. For example, daily increases in news intake were associated with subsequent increases in physical contact [0.021 (0.0062–0.036), adjusted p = 0.01].
Associations between the more stable components of the behavioral and psychological measures (i.e. the covariations between random intercepts) differed from cross-lagged effects (Table 2). For example, while individual differences in physical contact with others were associated with increased striving, alcohol and drug intake (i.e. random intercepts covaried), there were no significant cross-lagged, dynamic effects between the same variables.
The covariances between the random intercepts are given with 95% confidence intervals, with FDR-adjusted p values.
To examine whether race/ethnicity and income level significantly impact the main findings, we compared model fits in which lagged effects were allowed to vary by race/ethnicity against those that constrained the lags to be the same across race/ethnic groups. Fit indices in constrained and unconstrained models were nearly identical, and thus we conclude that moderation of the cross-lagged effects by racial/ethnic identity and income level did not emerge as a substantial effect (online Supplementary Table S7). Reports of a past diagnosis of either a depressive (N = 270) or anxiety disorder (N = 307) also did not emerge as significant moderators (online Supplementary Table S8).
We tested whether the impacts of COVID worry remained after treating worry as two factors: one pertaining to fear of becoming infected by SARS-CoV2 and the other relating to its impact on daily functioning. The main cross-lagged findings, wherein COVID-worry was significantly associated with news intake (in both directions) and substance use, held when COVID-related worry was treated as two factors (online Supplementary Table S9). Significant findings in the main analysis also held after covarying for gender, age, date of study entry, and state trends in confirmed infections (online Supplementary Table S10).
Discussion
We map a vicious cycle within individuals between news/media consumption and COVID-19-related worry. Watching more news than usual on a given day significantly increased COVID-19-related worry the next day, and vice versa. An adverse effect of media consumption was also detected on measures of psychological struggling. Notably, psychological striving reduced media consumption, though reduced media consumption did not impact subsequent striving.
We did not find cross-lagged effects between physical and virtual contact with others on subsequent change in mental health. This is surprising given previous evidence of social isolation having an adverse impact on mental health: a longitudinal study of older adults has shown that social disconnectedness and perceived isolation lead to increased depression and anxiety symptoms (although that study was not conducted during the current pandemic) (Santini et al., Reference Santini, Jose, York Cornwell, Koyanagi, Nielsen, Hinrichsen and Koushede2020). It is however noteworthy that the more stable measure of individual differences in overall levels of physical contact was associated with greater striving and less struggling (i.e. the respective random intercepts were correlated). It is possible that the strong associations between more stable individual differences in behavior and mood overshadowed the impact of daily fluctuations in mood and behavior on one another, as characterized by the cross-lagged effects. However, the direction of effects cannot be inferred from associations between the stable individual differences; such directionality can only be deduced from the cross-lagged effects.
Individual differences in virtual contact emerged as a double-edged sword, being associated with both positive and negative facets of mental health. The positive association may reflect the role of virtual contact in promoting a sense of social connectedness; the negative association may reflect the heavy presence of COVID-related information on social media. As expected, daily increases in the number of substances used had an adverse impact on mental health, decreasing subsequent levels of psychological striving. This adds to the evidence in favor of moderating substance use during the pandemic, particularly as individuals with substance use disorders may be at increased risk of being infected with COVID-19 (Ornell et al., Reference Ornell, Moura, Scherer, Pechansky, Kessler and von Diemen2020). However, this pattern does not account for fluctuations in the amount of other substances used, which was not assessed.
The analytic approach used (RI-CLPM) has several advantages over traditional cross-lagged analyses. Traditional cross-lag panel models are unable to fully account for variables that are at least partially time-invariant through just the autoregressive correlations (Hamaker et al., Reference Hamaker, Kuiper and Grasman2015). Additionally, by separating between-person and within-person differences, the RI-CLPM has also been found to provide more accurate mapping of the direction of cross-lagged effects than classic cross-lagged analytic approaches (Hamaker et al., Reference Hamaker, Kuiper and Grasman2015). Finally, RI-CPLM provides better fitting models than traditional cross-lagged analyses due to the capture of trait-like stability through the inclusion of random intercepts (Mulder & Hamaker, Reference Mulder and Hamaker2020). Nonetheless, while the RI-CLPM approach parses the direction of effects between behavior and mental health, it does not demonstrate causal connections, which requires other experimental designs.
The study has several limitations. Firstly, the survey was offered only in English, limiting the participation of some populations. The survey also required consistent Internet/smartphone access, possibly excluding some lower-income individuals. We limited survey items to not overburden participants, and thus coverage of mental health was narrowed. Twenty-four individuals tested positive for SARS-CoV2 during the study, a rate of approximately 1% per week, which is in line with general population rates. Given the small size of this group, we were unable to draw contrasts between those with and without confirmed infection. The study occurred at a phase in the pandemic when most states had implemented strict social distancing policies, and we are thus unable to assess the impact of the ongoing relaxation of these policies. Finally, our daily measures of substance use were limited, particularly for substances beyond alcohol, where we collected data only on the number not quantity of substances used. This approach might inflate the impact of relatively light use of several substances and underestimate the impact of heavy use of one substance.
As noted in the introduction, we were interested in potential disparate impacts on racial/ethnic minorities and those with limited income; by design, we oversampled these groups. This sampling strategy could limit the generalizability of our findings to the general population. However, given that these populations are particularly vulnerable to COVID and the economic risks associated with the pandemic, we felt that it was important to oversample these groups. It is also important to note that the cross-lagged effects between behavior and mental health measures did not differ by racial/ethnic groups and income, suggesting that the findings may hold for the more general population. Finally, we note that a recent study of online recruitment resources found that Qualtrics panels came closest to a national probability sample on a majority of demographic variables and types of media consumption, outperforming other platforms [Mechanical Turk (MTurk) and Facebook] (Boas, Christenson, & Glick, Reference Boas, Christenson and Glick2020).
We did not find lagged effects to vary substantially by race and ethnicity, given that the fits of models that held lags constant or allowed them to vary by race/ethnicity did not differ substantially. Despite this null finding, we note disparities in rates of many other COVID domains, including the risk of SARS-CoV2 infection, hospitalization, and death in African-American and Hispanic populations during this pandemic (Price-Haywood, Burton, Fort, & Seoane, Reference Price-Haywood, Burton, Fort and Seoane2020; Webb Hooper, Napoles, & Perez-Stable, Reference Webb Hooper, Napoles and Perez-Stable2020). Additionally, diagnosis of neither a past depressive nor past anxiety disorder emerged as a significant moderator of the lagged effects, a contrast with prior studies finding associations between a history of mental illness and worse mental health outcomes during the pandemic (Alonzi, La Torre, & Silverstein, Reference Alonzi, La Torre and Silverstein2020; Kwong et al., Reference Kwong, Pearson, Adams, Northstone, Tilling, Smith and Timpson2020; O'Connor et al., Reference O'Connor, Wetherall, Cleare, McClelland, Melson, Niedzwiedz and Robb2020).
In conclusion, we find that a relatively simple behavior change – watching less news/media – may reduce subsequent adverse mental health outcomes. These findings move beyond prior reports of associations between watching the news on highly distressing topics and poor mental health by demonstrating a dynamic impact of news on mental health that unfolds over days within the individual (Gao et al., Reference Gao, Zheng, Jia, Chen, Mao, Chen and Dai2020; Ni et al., Reference Ni, Yang, Leung, Li, Yao, Wang and Liao2020; Olagoke et al., Reference Olagoke, Olagoke and Hughes2020; Pfefferbaum & North, Reference Pfefferbaum and North2020). The benefits of media as an information source and means of engagement in public discourse may be balanced against its possible deleterious impact on mental health. Our finding that levels of psychological striving were associated with subsequent decreases in news intake provides another possible point of intervention. Thus, interventions that promote facets of striving such as happiness and optimism could help break the vicious cycle between news/media consumption and COVID-related worry that we delineate.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291721001896
Data
The de-identified data that support the findings of this study will be made available upon request to the corresponding author.
Conflict of interest
None.
Financial support
PS and LK received intramural NHGRI/NIMH awards from the National Human Genome Research Institute and the National Institute of Mental Health, as well as the National Institutes of Health's Office of Behavioral and Social Science Research.