The Impact of the COVID-19 Pandemic on Mental Health, Occupational Functioning, and Professional Retention Among Health Care Workers and First Responders

Background The COVID-19 pandemic has greatly affected front-line health care workers (HCW) and first responders (FR). The specific components of COVID-19 related occupational stressors (CROS) associated with psychiatric symptoms and reduced occupational functioning or retention remain poorly understood. Objectives Examine the relationships between total and factored CROS, psychiatric symptoms, and occupational outcomes. Design Observational, self-report, single time-point online assessment. Participants A total of 510 US HCW (N = 301) and FR (N = 200) with occupational duties affected by the COVID-19 pandemic. Main Outcomes and Measures CROS were assessed using a custom 17-item questionnaire. Post-traumatic stress disorder (PTSD), depression, insomnia, and generalized anxiety symptoms were assessed using the PTSD Checklist-5 (PCL5), Patient Health Questionnaire-9 (PHQ9), Insomnia Severity Index (ISI), and General Anxiety Disorder-7 (GAD7). Respondents’ likelihood of leaving current field and occupational functioning were assessed with 2-item PROMIS subscales. Relationships were modeled using multivariable regression. Open-ended responses were coded using rapid template analysis. Results CROS total scores correlated significantly with all four psychiatric symptom domains (R’s = .42–.53), likelihood of leaving one’s current occupation (R = .18), and trouble doing usual work (R = .28), all p’s < .001. Half of HCW indicated a decreased likelihood of staying in their current occupation as a result of the pandemic. CROS were fit to a 3-factor model consisting of risk, demoralization, and volume factors. All CROS factors were associated with psychiatric symptom burden, but demoralization was most prominently associated with psychiatric symptoms and negative occupational outcomes. Among psychiatric symptoms, PTSD symptoms were most strongly associated with negative occupational outcomes. Open-ended statements emphasized lack of protection and support, increased occupational demands, and emotional impact of work duties. Conclusions and Relevance These results demonstrate potentially treatable psychiatric symptoms in HCW and FR experiencing CROS, impacting both wellbeing and the health care system. Mitigating CROS, particularly by addressing factors driving demoralization, may improve HCW and FR mental health, occupational functioning, and retention. Supplementary Information The online version contains supplementary material available at 10.1007/s11606-021-07252-z.


INTRODUCTION
Health care workers (HCW) and first responders (FR) working during a pandemic experience a variety of acute and sustained stressors, including fear for their own safety and that of their colleagues and family, exposure to death and suffering, separations from family, and prolonged periods of exhaustion and vigilance. They may also experience demoralization 1 related to inadequate support or seeing their contributions as ineffective. The impact of these types of experiences has been well documented in contexts such as military service, where both isolated and cumulative trauma exposure is strongly associated with high psychiatric symptom burden, as well as decreased occupational functioning and workplace participation [2][3][4] . Understanding how such occupational stressors affect HCW/FR is important to identify risks to HCW/FR and to our health care system, and develop strategies to reduce those risks 5 .
High rates of psychiatric symptoms have been documented in HCW working during the COVID-19 [6][7][8][9][10][11] and prior pandemics 12,13 . Pandemic-related stressors identified as potential risk factors include quarantine procedures, isolation from social supports, stigmatization, and infection risk 9,10,14,15 . Less is known about the experiences of FR, despite exposure to similar pandemic-related stressors, often in less controlled environments. Consistent with the potential for high risks for this group, a survey of HCW and emergency medical services (EMS) workers in Italy during the COVID-19 pandemic found high levels of distress in both groups, but increased anger and regret, increased intrusiveness related to trauma, and decreased perception of self-efficacy in EMS relative to HCW 16 .
Physical and emotional stressors among HCW/FR have also been associated with decreased professional longevity and poorer patient care 17,18 . The potential for significant professional attrition of HCW related to the COVID-19 pandemic has been emphasized in abundant media reports 19,20 as well as organizational surveys 21 , and a review of medical leave in firefighters and emergency medical services (EMS) workers in New York early in the COVID-19 pandemic found increased use of medical leave, leading to decreased workforce availability 22 . However, little published data are available regarding the factors driving or modulating the risk of increased attrition or decreased occupational functioning among HCW/FR during the COVID-19 pandemic.
To address these gaps, we conducted an observational, crosssectional study among HCW (clinical and support staff) and FR (EMS, fire, and law enforcement officers [LEO]) working in the USA during the COVID-19 pandemic. The goal of this study was to characterize and assess associations among COVID-19related occupational stressors (CROS), psychiatric symptoms (depression, anxiety, insomnia, and post-traumatic stress disorder [PTSD]), and self-reported functional impairment and likelihood of leaving one's current field. CROS were quantified as a total score and by factor analysis, to identify specific types of CROS and their potential for differential impact. It is hoped that a better understanding of these relationships will facilitate the development of targeted interventions to protect HCW and FR even during periods of increased risk and workload.

METHODS
The study was approved by the VA Puget Sound Health Care System Human Subjects Committee. Prior to enrollment, all participants were provided an information statement that detailed the purpose, risks, benefits, and alternatives to participation.

Participants
A convenience sample of 510 participants was recruited through targeted outreach and paid advertising on Facebook between September 15, 2020, and February 7, 2021. Targeted outreach included large and COVID-focused HCW Facebook groups and emails to professional organizations (e.g., unions) and list-serves. Participants were asked to self-attest that they were a HCW or FR who provided professional services affected by the COVID-19 pandemic . The sample for the present  study includes 301 HCW (60 physicians, 188 RN/LPN) and  200 FR (162 EMS, 54 firefighters, and 19 law enforcement  officers [LEO]). LEO and firefighters were merged for subgroup analyses below, and individuals reporting dual EMS and LEO/firefighting roles were categorized under LEO/fire. Responses spanned 47 states and 445 zip codes, with broad distribution across the rural-urban continuum 23 (Supplementary Fig. 1).

Procedures
Self-report assessments were collected using Qualtrics. To encourage broad and representative participation, participants were not required to provide their legal name and were able to skip questions. Email addresses were collected to allow longitudinal follow-up for up to 9 months. Compensation was not provided. Data presented represent an analysis of the baseline survey only, except for the free-text response analysis. In order to include more respondents' perspectives to this optional item, and a broader range of themes, free-text responses from any time point during the same calendar period were included in this analysis.

Measures
COVID-19 Related Occupational Stressors. Exposure to CROS was assessed using a measure designed for this study by a collaborative team of physicians working in New York City hospitals and emergency rooms in March and April 2020 and our research team (Appendix A). This instrument asks four yes/no questions about personal experiences of loss due to COVID-19, and 13 questions assessing the frequency over the past 2 weeks of caring for individuals with COVID-19, perceiving increased risk to self or family due to occupational exposure to COVID-19, experiences of patient suffering related to the impacts of COVID-19, inadequate support or protection related to COVID-19, and feeling unable to provide effective or adequate care due to COVID-19, with responses scored 0-3. The 13-Likert scale items were summed to provide a CROS total score, while the initial 4 yes/no questions were used to characterize the sample and as covariates in multivariable regression models. Following this structured assessment, participants were asked an optional open-ended item, "Is there anything else you would like to share with us that is increasing your stress from your occupational duties?".
Psychiatric Symptoms and Prior Trauma Exposure. Psychiatric symptoms were assessed using the Patient Health Questionnaire 9-item (PHQ9) 24 , the Generalized Anxiety Disorder 7-item (GAD7) 25 , the Insomnia Severity Index (ISI) 26 , and the PTSD Checklist for DSM-5 (PCL5) 27 . Prior trauma history was assessed using the Life Events Checklist 28 ; a numerical index of prior trauma exposure was calculated using items selected a priori as those most likely to meet criterion A of PTSD (items 7-8, 11, and 14-16), and constituted the sum of items a participant reported they had witnessed or experienced.
Occupational Outcomes. Functional impairment was assessed using the two work-related items from the PROMIS Short Form v2.0 Ability to Participate in Social Roles and Activities 8a measure 29 , modified to focus on occupational work (Supplemental Table 1). The likelihood of leaving one's current profession was assessed with two items: "How likely do you think it is that you will still be working in your current field in 5-10 years?" and "How have your experiences providing care during the COVID-19 pandemic affected your interest, willingness, or ability to continue working in your current field?".

Data Analysis
Data were analyzed using R and RStudio. Correlation coefficients were implemented using the stat_cor function from ggpubr 30 . Pearson's R was used for variables representing scale totals, while Spearman's correlation coefficients were used when one or more of the variables represented a single ordinal item or the sum of two ordinal items. Figures were created using ggplot2 31 and jtools 32 , tables using flextable 33 . For all analyses, we selected a p value of < 0.05 to determine statistical significance. Participants with missing data were excluded from analyses using items or measures with missing data. Missingness ranged from 2.8 to 30% (Supplementary  Table 2). Most missing data represented partial survey completers rather than missing items, so no data were imputed. Completers of the instruments that came later in the survey (e.g., GAD7, PHQ9) were more likely to be married, physicians, to have had COVID, and to report a higher number of prior trauma exposures.
Characterization and factor analysis of the CROS was implemented using psych and psychTools 34 . Cronbach's alpha was 0.88 (CI 0.86-0.9), indicating strong internal consistency. Following the scree plot (Supplemental Fig. 2), a 3-factor solution was chosen. Using varimax rotation produced a solution with mean item complexity of 1.5 and RMSR of 0.04. These data-driven factors emphasized 3 conceptually coherent, face-valid domains (Fig. 1). Items characterizing a total Fig. 1 Factor analysis of COVID-19-related occupational stressors (CROS). Results of a 3-factor analysis of COVID-19 exposure items into 3 factors, termed volume, demoralization, and risk, based on an interpretation of their most highly weighted items. Color shading is proportional to the numeric weight on each cell, and indicates the weight that item contributes to the factor below. See Appendix A for complete wording of CROS items volume of COVID-19-related care delivered were most heavily weighted by factor 1, which was termed the "volume" factor. Factor 2 most heavily weighted items addressing the ability to provide high quality and effective care to all patients, being asked to take unnecessary risks, or being unsupported by one's workplace; this factor was termed the "demoralization" factor. The third factor most heavily weighted one's personal risk or one's family's risk of contracting COVID-19, termed the "risk" factor. Multivariable linear regression models were implemented using nlme 35 , with either continuous CROS factor scores or continuous psychiatric symptom domain scores as the predictor variables. All multivariable models included age and gender as covariates. For models characterizing the relationship of CROS factors to psychiatric symptoms and occupational outcomes, an index of prior trauma exposure, personal history of COVID-19 infection, history of COVID-19 infection in a family member, and death of a family member or close colleague due to COVID-19 were additionally included as covariates.
Analysis of free-text responses was carried out using a rapid template qualitative analysis, an established rapid analysis method ideally suited for pragmatically describing textual data (e.g., the restricted range and fragmented text for this single open-ended questionnaire item). 36,37 A coding template summarizing broad "codes" was initially developed by the senior qualitative analyst from reading a random sample of responses. Initial codes were then applied independently by two coders. Codes were updated as needed when discrepancies or lack of clarity occurred. Discrepancies were resolved until an intercoder agreement was reached. Responses were excluded if they were not related to COVID (n = 1) or were redundant responses by the same respondent across assessments (n = 2).
Initial coding was done with the combined HCW and FR sample to identify relevant themes. We then performed stratified counts and identified issues most salient to HCW or FR; i.e., if an overrepresented proportion of that occupational category endorsed a theme, defined as > 5% higher than would be expected given the overall proportion of open-ended respondents from that occupational group. Table 1 presents sample characteristics, including demographics, CROS scores, psychiatric symptoms, and outcome assessments, for the entire sample and stratified by subgroups. In addition to comparing HCW and FR, sufficient respondents from the subgroup physicians, nurses, LEO plus fire, and EMS were present to analyze these groups independently. In general, responses were similar for HCW and FR, including similar scores on CROS measures. However, nurses scored significantly higher than physicians on both total exposure (CROS total, p < 0.01) and most psychiatric symptoms domains (PTSD, depression and insomnia p < 0.01, thoughts of suicide or self-harm p = 0.04). Similarly, CROS total was significantly higher for EMS versus LEO/fire (p < 0.01), as were symptoms of depression (p < 0.01). Fewer FR (40.6%) than HCW (55%) reported experiences working during the pandemic had decreased their likelihood of remaining in their current field (p = 0.004).

RESULTS
Bivariate Relationship of CROS Total Score to Psychiatric Symptoms and Occupational Outcomes Figure 2 presents bivariate associations between demographic variables, exposure scores, and psychiatric symptom scores. Of the psychiatric symptom domains, the CROS total was most strongly correlated with the PCL5 (R = 0.52) but was also significantly associated with the PHQ9 (R = 0.44), ISI (R = 0.41), and GAD7 (R = 0.43, all p < 1e-15). CROS total also was significantly correlated with thoughts of suicide or self-harm (R = 0.25, p < 1e-5), problems completing work tasks (R = 0.28, p < 1e-8), and likelihood of leaving field (R = 0.18, p < 0.001). For the subgroup that reported at least moderate insomnia symptoms and was asked about the impact of sleep problems on functioning at work (N = 367), CROS total was significantly related to reported impact of sleep problems on work performance (R = 0.37, p < 1e-11).

Multivariable Relationships of the 3 CROS Factors to Psychiatric Symptoms and Occupational Outcomes
We conducted a series of multivariable linear regression models. The relationships of CROS factors to psychiatric symptoms across all respondents are presented in Fig. 3A, and by subgroup in Supplemental Fig. 3. These models found that all 3 CROS factors were significantly and positively related to all 4 psychiatric symptom domains. Across all participants and for both HCW and FR, the demoralization factor was the strongest correlate of PTSD (HCW β = 0.37, FR β = 0.59), depression (HCW β = 0.30, FR β = 0.5) and anxiety symptoms (HCW β = 0.29, β = 0.42, all p < 0.0001). The relationships were similar across the four small subgroups (physicians, nurses, LEO + fire, and EMS) with the exception that for LEO + fire, the volume factor was as or more strongly associated with these outcomes than demoralization (β = 0.58 GAD7, β = 0.47 PTSD, β = 0.44, all p < 0.01). For insomnia symptoms, the demoralization factor was strongly associated with symptom intensity for physicians (β = 0.41, p < 0.05) and EMS (β = 0.35, p < 0.01), while the risk factor was associated with symptom intensity for nurses (β = 0.36, p < 0.01).
Across all respondents, thoughts of suicide or self-harm were significantly and positively related to both the volume  factor and the demoralization factor, as well as to prior trauma history ( Fig. 3B and Supplemental 3B). On subgroup analysis, this pattern was maintained for nurses and EMS, but not for physicians (no significant predictors) or LEO + fire (volume factor β = 0.95, p < 0.001).
In examining the relationship of CROS factors to occupational outcome measures across all respondents, only the demoralization factor was significantly related to the likelihood of leaving one's current field or problems completing work-related tasks. This pattern was preserved across subgroups, with the exception that for LEO + fire demoralization was significantly related to the likelihood of leaving one's field (β = 0.52, p < 0.01) but not to occupational functioning, while for EMS demoralization was significantly related to occupational functioning (β = 0.55, p < 0.0001) but not the likelihood of leaving one's field.

Relationship of Psychiatric Symptom Domains to Functional Outcome Measures and Suicidality
The relationship of psychiatric symptom domains (PTSD, anxiety, depression, and insomnia) to functional outcome measures and thoughts of suicide or self-harm were similarly characterized ( Fig. 3C and Supplemental 3C). Across all participants, PTSD symptom severity was significantly related to the likelihood of leaving one's current field, trouble completing work tasks, and thoughts of self-harm or suicide. In addition, older age had a significant positive relationship to the likelihood of leaving one's current field, and depression symptoms were significantly and positively related to thoughts of self-harm or suicide. However, based on subgroup analyses, the relationship of PTSD symptoms to the increased likelihood of leaving one's current field was driven most strongly by nurses (β = 0.41, p < 0.01) and LEO + fire (β = 0.74, p < 0.05), while the relationship of PTSD to occupational functioning was driven most strongly by physicians (β = 0.67, p < 0.05) and EMS (β = 0.43, p < 0.01). Interestingly, for LEO + fire, anxiety symptoms were strongly and positively associated with thoughts of suicide or self-harm (β = 1.2, p < 0.001) but strongly and negatively associated with thoughts of leaving one's field (β = 0.91, p < 0.01).

Qualitative Analysis of Free-Text Responses
Over one-third of participants (36% of total, 37.5% of HCW and 32.4% of FR) responded to the optional, open-ended question about sources of occupational stress during the COVID-19 pandemic, with an average response length of 31 words (range: 1-169). Responses vividly conveyed the Fig. 2 COVID-19-related occupational exposure (CROS total score) is strongly related to increased burden of psychiatric symptoms. A Bivariate relationships between demographic variables, exposure scores, and psychiatric symptom scores are represented by Spearman's correlation coefficients. Scatter plots are provided for the relationship of CROS total to PTSD symptoms (B), depression symptoms (C), insomnia symptoms (D), and anxiety symptoms (E). *p < .05, **p < .01, ***p < .001, ****p < .0001. "Close death" = death of a family member or close colleague from COVID-19, "Inc risk" = medical condition associated with increased risk from COVID-19 infection Fig. 3 Relationships between different factors of COVID-19-related occupational stressors (CROS factors), psychiatric symptom expression, and functional outcomes. Results of multivariable regression models relating CROS factors and covariates to psychiatric symptom clusters (A) and functional outcome measures along with thoughts of suicide or self-harm (B). C Results of independent multivariable regression models relating symptom clusters as measured by total scores on the PCL5, PHQ9, GAD7, and ISI, along with covariates of age and gender, to functional outcome measures challenges faced by HCW and FR and the magnitude of associated distress. A rapid template analytic approach identified 12 overarching codes within respondent data, which are described, along with illustrative quotes in Table 2.
Overall, responses were highly consistent with the quantitative analyses above. Among the 12 identified codes, "Lack of Protection and Support," "Increased Demands," and "Emotional Toll" were the most common. Respondents described their work during COVID-19 as involving "some of the saddest death stories" they had experienced and feeling "spread thin and exhausted." Notable for a relative lack of representation were comments addressing respondents' personal risk of COVID-19 infection, with most references to this type of risk taking the form of concern for the impact of this risk on others.
HCWs were overrepresented (> 5% deviation) in categories "staffing shortages," "increased demands," and "patient care." FRs were overrepresented in "fear of or enacted reprisal from leaders," "concern for the well-being of others," and "betrayal by colleagues." HCW and FR were represented in all categories.  "People begging for your help. I feel so evil and dirty having to place a BiPap on a patient begging me not to. They don't like it and cry and beg for me to let them die. I must put patients in restraints to keep them from pulling out their tubes. They cry for me to let them go. It's like a bad horror movie." Fear of or enacted reprisal from leaders (TOTAL RESPONSES = 7; HCW: 57%; FR: 43%) • Fear of or enacted retaliation for speaking out against workplace risks, hazards, and lack of safeguards; blackmail; threats to career and/or professional development for wanting to quit "I was exposed to a confirmed COVID patient in April who was in respiratory distress. I only had my N95 on, no face shield or gown had been provided. I wanted to quarantine; I was accused of borderline patient abandonment. They threatened to report me to the board. I quit, then I couldn't get unemployment because I quit." Betrayal by Colleagues (TOTAL RESPONSES = 10; HCW: 60%; FR: 40%) • Failure of colleagues to follow guidelines/safeguards "My coworkers are COVID deniers. I work in EMS and it makes it really hard." Betrayal by Society (TOTAL RESPONSES = 28; HCW: 71%; FR: 25%; OTHER: 4%) • Public health guidelines disregarded by the general public; COVID-19 skepticism and/or denial often fueled by political and public health leaders; lack of accommodation from society for needs such as childcare "The worst thing is dealing with incredible stress at work, and then realizing no one really cares… I separate from my kids at the first sign of symptoms because I'm heavily exposed at work, but then have to listen to people complaining about recommendations they don't have people over for Thanksgiving. It's very disheartening when the community doesn't do its part. I feel betrayed. We determined responses as being overrepresented when the proportions of respondents who endorsed the item deviated by > 5% from the occupational category's representation in the open-ended text responses (i.e., 67% HCW, 31% FR, and 2% other). For example, if 73% of a theme was endorsed by HCW, it was considered overrepresented, and thereby uniquely relevant to that occupational group

DISCUSSION
Our findings are consistent with previous demonstrations of high levels of psychiatric symptoms and distress in HCW working during the COVID-19 pandemic [6][7][8][9][10][11] . In the present study, a quantitative measure of CROS was associated with psychiatric distress. The relationships remained significant when HCW and FR were analyzed independently and when a quantitative measure of exposure to prior traumatic stressors was included as a covariate, suggesting the high rates of psychiatric symptoms are unlikely to represent baseline rates in this population independent of COVID-19-related occupational stressors.
The overall magnitude of most relationships identified was similar, or sometimes greater, in FR vs. HCW. This was particularly striking and concerning regarding thoughts of suicide or self-harm, reported by 19% of FR and 12% of HCW and strongly related to CROS. Given prior work demonstrating high trauma exposure and psychiatric symptoms in FR 38 , including high rates of both suicidal thoughts and behaviors, particularly in EMS 39 , results highlight the importance of addressing the impact of working during a pandemic on FR. Similarly, findings underscored the higher CROS exposure and psychiatric distress in nurses as compared with physicians, building on prior work finding increased distress, suicidal ideation, and suicidal behavior specifically in nurses 40 .
A substantial proportion of both HCW and FR reported their likelihood of staying in their current field had been somewhat or significantly decreased by their experiences working during the COVID-19 pandemic and that they at least sometimes had difficulty completing work-related tasks. These results are consistent with and build upon findings from previous pandemics 12,41 and are particularly worrisome given existing concerns about current and future shortages in the health care workforce 42,43 . The present study suggests healthcare staffing shortages are in and of themselves a COVID-19 stressor, so further shrinking the labor pool could have an exponential negative impact on HCW wellbeing and professional retention.
The elements in the CROS were well-described by three face-valid factors, termed "volume," "demoralization," and "risk." Demoralization showed the strongest relationship to both psychiatric symptom domains and functional outcome measures, while risk showed the smallest relationship. These findings suggest that while strategies such as vaccination that decrease the risk posed to HCW and their families by COVID-19 infection are important, they are unlikely to fully mitigate the impact of COVID-19-related occupational stressors on HCW and FR mental health or functional outcomes. In fact, the theme of betrayal by the community in the free-text responses raises the concern that high volumes of COVID-19related care driven by patients who have declined vaccination may be associated with increased psychiatric distress compared with similar volumes of care when vaccination was insufficient or simply not available.
Our results suggest a significant number of strategies that could help mitigate the effects of CROS on HCW and FR (Fig. 4), both by decreasing the volume of COVID-19 impacted care individuals are providing and the associated personal risk they experience while doing so (direct factors) and by changing the context in which this care is being provided (contextual factors). Although it is important to prospectively test the impact of such interventions, existing evidence supports the efficacy of systematic approaches to identifying and addressing factors that unnecessarily increase workload, while supporting health care workers' control and flexibility, meaning in work, and workplace community 5,[44][45][46] . Identifying and rectifying barriers to communication and transparency within a workplace, providing regular and effective feedback to leadership, and increasing the frontline worker input are often critical intervention components 5,[47][48][49] . Such interventions have not been found to come at the cost of other organizational goals or significant financial expense 48 . However, these approaches are most effectively implemented at the organizational level [47][48][49] .
Among psychiatric symptom domains, PTSD symptoms stood out for both the strength of their relationship with overall COVID-19-related occupational stressors, and as the symptom domain most related to adverse occupational outcomes. Along with depression symptoms, PTSD symptoms were also significantly related to thoughts of suicide or self-harm. This suggests that the identification and implementation of interventions to reduce the risk of PTSD [50][51][52][53] in HCW and FR should be a top priority. Moreover, PTSD symptoms may be particularly important to detect and treat in HCW and FR. Broadly implemented strategies to reduce risk and strengthen resilience while proactively identifying and providing accessible, confidential care for those who require more intensive services, maybe a particularly effective and cost-efficient approach 54,55 .
The current work has limitations. Our results highlight the significant variety of experiences of HCW and FR working during the COVID-19 pandemic and the importance of these differences for outcomes. However, our data do not include a detailed assessment of factors related to financial resources, family obligations, or position and influence within one's specific field or the health care system more broadly, each of which are likely to interact with many of the factors explored here 56 . Follow-up studies of the persistent impacts of working during previous pandemics underscore the need for longitudinal follow-up of the affected workforce 13,41 , and further work following the impact of CROS on psychiatric symptom burden and occupational functioning over time will be important. The work relies upon participants' self-report of exposure history, workplace experiences, and psychiatric symptoms, which may result in same source bias. The survey instrument was long, and some participants did not complete all instruments. In addition, the study used targeted outreach and paid advertising targeting regions with high rates of COVID-19 cases. The results obtained from the current sample may not reflect rates in all health care workers, nationally or internationally. At the same time, the broad geographic and rural/urban diversity of study respondents is a strength, particularly given the geographic variability in pandemic intensity and community response.
Finally, the relationships characterized are observational and cross-sectional, and cannot be assumed to represent unidirectional, causal relationships. For example, PTSD symptoms were found to be strongly associated with suicidality and negative occupational outcomes. This suggests that prevention and/or treatment of PTSD symptoms may decrease the risk of suicidal behavior and improve functioning and workforce retention. However, this hypothesis will need to be tested in a prospective trial.

Supplementary Information
The online version contains supplementary material available at https://doi.org/10.1007/s11606-021-07252-z. Acknowledgements: The authors would like to thank the many participants who generously shared their time and experiences with us. We would also like to acknowledge the significant number of staff members who contributed to this study often well beyond their usual work roles, including Carolyn LH Fort, Joseph D Clark, Kim Hart, Adam T McPartlin, Emma C Onstad-Hawes, Soleil S Groh, and Andrew Nicholls. We also thank George Sayre, Ashley Mog, and the Seattle-Denver COIN qualitative core group, who provided brief consultation on the qualitative methods, analysis, and interpretation. Author Contribution RCH designed the study, with input from RAS, BPC, CAM, and MAR. RCH supervised the implementation of the study, with RAS playing a significant role in all implementation steps. RCH carried out the analysis of all results with consultation by RGT, with the exception of the thematic analysis of free-text responses which was carried out by KDH and GRM. RCH wrote the initial draft of the manuscript, with contributions from RAS, KDH and MAR. All authors provided input into the interpretation of the analysis and the drafting of the manuscript.  Fig. 4 Potential schematic framework for considering direct and contextual factors contributing to occupational stress from the COVID-19 pandemic, and potential mitigation strategies. Stressors (left) and potential mitigation strategies (right) are divided into direct and contextual factors. Direct factors result primarily from the volume of COVID-19-related care being provided by an individual or their institution and resources available for the system to respond to these demands. Contextual factors can be addressed independent of the volume of COVID-19related care being provided, and include the responsiveness of the system to addressing and supporting HCW/FR's needs, and ensuring they are not put at unnecessary risk. Stressors represent a synthesis of factors identified from the quantitative and qualitative analyses; mitigation strategies represent concrete examples of ways in which the identified stressors could be modified, minimized, or mitigated. Mitigation or intervention approaches may vary depending on the most relevant occupational stressors for a specific group. For example, the strong relationship between demoralization and both psychiatric symptoms and adverse occupational outcomes in FR, along with the emphasis in freetext responses of fear of or enacted reprisal from leaders and betrayal by colleagues, suggest interventions focused on responsiveness and clear communication from leadership and protections of job and financial security may be particularly important for many FR. The high rates of PTSD symptoms and the relationship of these symptoms to a high likelihood of leaving one's current field for nurses may suggest that interventions focused on decreasing the risk of PTSD, and increasing the availability and utilization of treatment for PTSD, may be of particularly high priority for nurses. An alternative example of a conceptual framework for planning risk mitigation and interventions based on a literature review can be found in Schwartz et al 5