The COVID-19 pandemic emerged as a unique societal trauma that had a substantive impact on the psychological well-being of the global population. Reference Shevlin, Butter, McBride, Murphy, Gibson-Miller and Hartman1 Clinically high levels of anxiety-depression (21%) and post-traumatic stress (17%) were reported in the general population during the pandemic. Reference Shevlin, Butter, McBride, Murphy, Gibson-Miller and Hartman2 Recent studies estimate that, whereas there has been a drop-off in mental health difficulties from peak levels, post-pandemic rates remain elevated compared with pre-COVID-19. Reference Chen, Xiong, Jinfeng, Jun, Jiali and Xuan3,Reference Jamshaid, Bahadar, Jamshed, Rashid, Afzal and Tian4 Such cross-sectional measurements can be misleading, however, because the long-term course of mental health symptoms has been evidenced as non-uniform, with variation across specific groups within the wider population. Reference Solomon, Mikulincer, Ohry and Ginzburg5
Researchers conducting longitudinal studies have observed notable trends. López-Castro et al Reference López-Castro, Papini, Bauer, Swarbrick, Paul and Nizzi6 assessed post-traumatic stress disorder (PTSD) levels in the general public at four time points between April 2020 and July 2021. Four groups emerged with different symptom trajectories across the measurement points, specifically (a) a large subgroup that did not develop clinical levels of PTSD (73%, ‘resilient’ group); (b) a subgroup with clinical levels of PTSD at baseline that decreased to subclinical levels midway before increasing to clinical levels at final measurement (13.3%, ‘recurring’ group); (c) a subgroup with clinical levels of PTSD at baseline that decreased to subclinical levels at final measurement (8.3%, ‘recovering’ group); and (d) a subgroup with clinical levels of PTSD at baseline that were maintained at final measurement (5.5%, ‘chronic’ group).
Despite the explanatory utility of the findings, this study has several methodological limitations, including the fact that participants were not required to rate their PTSD symptoms using COVID-19 as the index trauma. Instead, participants could rate symptoms in relation to any potentially traumatic event they had experienced (e.g. physical assault). Moreover, the longitudinal measurements were taken during only the first year of the pandemic and did not extend to post-pandemic, which makes it questionable to draw conclusions about post-traumatic reactions following a societal trauma. Shevlin et al Reference Shevlin, Butter, McBride, Murphy, Gibson-Miller and Hartman1 corrected for the first of these issues, anchoring symptom ratings to COVID-19 in their study and measuring across five time points during the first year of the pandemic. That investigation found five classes, including three that were comparable to the groups of López-Castro et al: Reference López-Castro, Papini, Bauer, Swarbrick, Paul and Nizzi6 ‘resilient’, ‘recovering’ and ‘chronic’ (i.e. ‘resilient’ (67.7%), ‘adaptive’ (8.4%) and ‘chronic’ (4.9%), respectively), but with two additional classes. These were (a) a ‘moderate-stable’ subgroup that had consistently moderate clinical levels across time points; and (b) a ‘deteriorating’ subgroup with subclinical levels of PTSD at baseline that deteriorated across the pandemic to clinically high levels at final measurement. Reference Shevlin, Butter, McBride, Murphy, Gibson-Miller and Hartman1
Although the studies cited thus far provide key insights into the general long-term course of PTSD, as well as post-COVID-19 trauma responses, it is widely acknowledged that during the pandemic several subpopulations were at increased risk of post-traumatic reactions. Higher levels of clinically significant symptoms were reported in healthcare workers on anxiety (27%), depression (36%) and PTSD (32%) compared with the wider population. Explanations for these heightened symptoms include enhanced exposure to psychological and organisational stressors in healthcare staff job roles. Reference Jordan, Shannon, Browne, Carroll, Maguire and Kerrigan7,Reference Scott, Stevelink SA, Gafoor, Lamb, Carr and Bakolis8 Jordan et al Reference Jordan, Shannon, Browne, Carroll, Maguire and Kerrigan9 used longitudinal data from four peri-pandemic time points spanning November 2020 to August 2021, and identified two groups of healthcare workers with different PTSD symptom trajectories. Eighty-four per cent of staff had post-traumatic stress symptoms that fluctuated in line with COVID-19 in-patient levels but remained at subclinical levels, whereas a substantial minority (16%) had persistent clinical levels of PTSD symptoms across all phases of the pandemic.
Aims
The present study examined the longitudinal course of PTSD symptoms within healthcare staff both during and after the COVID-19 pandemic. Considering that healthcare staff were a vulnerable population during the outbreak and now play a pivotal role in societal recovery from COVID-19, it was deemed important to identify at-risk groups who may have developed debilitating PTSD trajectories. The investigation followed on from Jordan et al, Reference Jordan, Shannon, Browne, Carroll, Maguire and Kerrigan9 and used latent class growth analysis (LCGA) to examine the trajectories of post-traumatic stress symptoms reported by healthcare staff over a 2-year period (February 2021 to May 2023). In line with the three trajectories identified in studies reported by López-Castro et al Reference López-Castro, Papini, Bauer, Swarbrick, Paul and Nizzi6 and Shevlin et al, Reference Shevlin, Butter, McBride, Murphy, Gibson-Miller and Hartman1 it was predicted that the majority of healthcare staff would maintain subthreshold levels of post-traumatic stress symptoms peri- and post-pandemic (i.e. ‘resilient’ group), whereas a minority of healthcare staff would develop clinically significant levels of symptoms that would either resolve post-pandemic (i.e. ‘recovering’ group) or continue at these higher levels (i.e. ‘chronic’ group). It was also hypothesised that personal, demographic and organisational stressors would predict membership of post-traumatic stress trajectory classes (i.e. ‘recovering’ and ‘chronic’ groups).
Method
Participants
All health and social care staff in Northern Ireland were eligible to take part in the COVID-19 Staff Wellbeing Survey, which was administered at five time points spanning both the peri- and post-pandemic periods. The present study used data from 4 of the 5 survey time points: Time 2 (T2., 8–28 February 2021; n = 2898); T3 (10–30 May 2021; n = 2480); T4 (9–29 August 2021; n = 2119); and T5 (15 May–12 June 2023; n = 2047). T2–4 occurred during the pandemic, and T5 coincided with the declaration by the World Health Organization in May 2023 that the COVID-19 global health emergency was over. 10,Reference Wise11 In Northern Ireland, T5 also occurred simultaneously with the collection of many COVID-19 health statistics being discontinued following a sustained and substantive reduction in deaths caused by COVID-19. 12
The survey was advertised to staff via broadcast emails, health service social media, screensavers and laminated posters in staff areas. Data from the four time points were merged using work or personal email addresses, which was an optional response field in the survey. Analyses were conducted on data reported by a sample of 439 healthcare workers who participated and provided their email addresses on at least 3 of the 4 time points. The demographic profile of the longitudinal sample is shown in Table 1. The gender and age profile of the sample is consistent with Northern Ireland Health and Social Care Workforce census data. 13
Demographic characteristics of the sample at Time 2

NA, not available.
For the longitudinal sample, Time 3 data for gender and occupation were used if Time 2 data were missing.
Participants in occupational groups with desk-based roles, and hence with greater access to computers at work to facilitate completion of an online survey, tended to be over-represented in the sample (administrative and clerical Time 2 34% v. population data 19%; professional and technical Time 2 27% v. population data 15% 13 ). By contrast, occupations such as nursing and midwifery (T2 18% v. population data 34%) were underrepresented.
Measures
Post-traumatic stress symptoms
The 22-item Impact of Event Scale-Revised Reference Weiss, Marmar, Wilson and Keane14 was used to assess post-traumatic stress severity. The index trauma for ratings was anchored to COVID-19, with participants reporting symptoms relating to the COVID-19 outbreak within the past 7 days at T2–5. Participants were specifically instructed to ‘indicate how distressing each difficulty has been for you DURING THE PAST SEVEN DAYS with respect to the COVID-19 outbreak’. Responses to the items on this DSM-IV-compatible measure were assessed using 5-point Likert items (0 = not at all; 4 = extremely; scale range 0–88). A cut-off of ≥26 was used to indicate clinically significant moderate-to-severe symptoms of post-traumatic stress symptoms, in line with other international studies on healthcare staff well-being during COVID-19. Reference Lai, Ma, Wang, Cai, Hu and Wei15
Predictor variables
A range of demographic, occupational and support variables as measured at T2 were considered as potential predictors of post-traumatic stress class membership. Demographic variables included gender (male/female); age (years); a measure of exposure to COVID-19 (scale range 0–7, higher scores being indicative of greater exposure); and if they had at least 1 of 10 COVID-19 risk factors (coding 0 = no; 1 = yes) listed in the COVID-19 Pandemic Mental Health Questionnaire Reference Rek, Freeman, Reinhard, Bühner, Keeser and Padberg16 (e.g. diabetes, chronic liver disease).
Occupational information recorded included whether they had managed patients with COVID-19 (0 = no; 1 = yes); whether their employer had asked them to consider a redeployment opportunity during COVID-19 (0 = no; 1 = yes); and the perceived effectiveness of communication from their organisation on COVID-19-related matters (0 = not effective; 4 = very effective). The demographic and occupational variables are described in full in Jordan et al. Reference Jordan, Shannon, Browne, Carroll, Maguire and Kerrigan7
Participants were shown a range of team supports (e.g. Schwartz rounds, ‘buddy’ system) and asked how many were available to them during the last 3 months of the pandemic. Similarly, a range of staff supports (e.g. staff well-being helpline, drop-in centre) were listed, and participants indicated which of these they had used in these last 3 months. From this, 4 variables were constructed: whether at least 1 form of team supports was available (0 = no; 1 = yes); whether they had used at least 1 form of staff supports (0 = no; 1 = yes); total number of team supports available (0–8); and total number of staff supports used (0–9).
Procedure
Health- and social care workers were invited to take part in the survey online via the Survey Mechanics platform. Participants indicated consent at each time point by clicking to start the questionnaire once they had read the participant information sheet, which provided full details of the study and of their rights as 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 2013. All procedures involving human subjects/patients were approved by the West of Scotland Research Ethics Service (reference 20/WS/0122, 26 August 2020).
Statistical analysis
The data were prepared for analysis in IBM SPSS version 30 for Windows (IBM, Armonk, New York, USA; https://www.ibm.com/products/spss), then transferred to Mplus version 8.8 for Windows (Mplus, Los Angeles, California, USA; https://www.statmodel.com/programs.shtml) Reference Muthén and Muthén17 for inferential analyses. Missing data rates by time point for the longitudinal sample were as follows: T2 (12%), T3 (6%), T4 (15%) and T5 (44%). Logistic regression was used to assess whether demographic, occupational and pandemic-related T2 variables (e.g. gender; age; perceived effectiveness of communication regarding COVID-19; exposure to COVID-19; awareness of available supports) predicted drop-out at T5. Analysis revealed that being female and having lower awareness of available organisational staff well-being supports, as measured at T2, significantly predicted drop-out at T5. Missing values in the unconditional growth curve models and unconditional latent class growth analysis were dealt with using full-information, maximum-likelihood estimation. Initially, a series of unconditional growth curve models of varying shape (linear, quadratic and free time scores) were fitted to the T2–5 post-traumatic stress scores, and fit was evaluated with reference to acceptable range guidelines. Reference Hu and Bentler18 The freed loading model provided the best fit to the data (x 2(3) = 5.51, P = 0.138; comparative fit index 1.00; Tucker–Lewis index 0.99; root mean square error of approximation 0.04; standardised root mean residual 0.02). In the freed loading model, T2 was coded as 0 and T5 as 1, meaning that growth rate represented change between T2 and 5.
Missing data on T2 predictor variables were dealt with via multiple imputation using the Markov chain Monte Carlo method to generate 20 data-sets, under the assumption of missing at random. Significant predictors of missingness at T2 were included in the models, alongside the outcome variable and T3–5 versions of the T2 predictor variables, because these were strongly related. There was significant variation in both intercept and slope (Table 2) parameters, indicating the possible existence of subpopulations. Reference Kandauda, Lee, O’Neal and Lorenz19 A series of LCGA models with 2–5 classes specified were then fitted to the post-traumatic stress T2–5 data. Fit statistics and class trajectory plots were inspected to determine the number of classes that provided best fit to the data. Reference Feldman, Masyn and Conger20–Reference Tein, Coxe and Cham23 It was determined that a three-class model provided the optimal solution to the data (Table 3). Rates of missingness across time points varied among the classes ‘chronic’ (T2 = 13%, T3 = 5%, T4 = 15%, T5 = 42%), ‘recovering’ (T2 = 6%, T3 = 13%, T4 = 0%, T5 = 63%) and ‘resilient’ (T2 = 9%, T3 = 9%, T4 = 15%, T5 = 46%).
Unconditional growth curve model (freed loading) for post-traumatic stress

a Intercept with slope covariance.
Fit indices for unconditional growth mixture models with one to five classes

AIC, Akaike information criterion; BIC, Bayesian information criterion; SABIC, sample-size adjusted BIC; LMR-LRT, Lo–Mendell–Rubin likelihood ratio test; BS-LRT, bootstrapped likelihood ratio test.
Adopting the three-step manual approach, Reference Vermunt, Magidson, Hagenaars and McCutchen24 predictors were then incorporated into the models. Initially a series of univariate multinomial logistic regression models were run to assess whether demographic (gender, age, exposure to COVID-19, having at least one COVID-19 risk factor), organisational (managed patients with COVID-19, perceived effectiveness of communication regarding COVID-19, whether asked to consider being redeployed) and support variables (whether team supports were available, whether used staff well-being supports, number of team supports available, number of staff well-being supports used) predicted class membership. The significant predictors (exposure to COVID-19, having at least one COVID-19 risk factor, perceived effectiveness of communication regarding COVID-19) were then entered simultaneously into a multinomial logistic regression model.
Results
A 3-class model provided the best fit for post-traumatic stress symptom T2–5 scores (Table 3).
Post-traumatic stress trajectories
Three distinct classes were evident from the three-class LCGA model. Whereas Akaike information criterion, Bayesian information criterion (BIC) and sample-size adjusted BIC continued to decrease throughout classes 1–5, scree plots Reference Nylund-Gibson and Choi25 revealed that no further substantive decreases occurred after 3 classes. Entropy was highest for 3 classes, and the average latent class probabilities for classes 2 and 3 were all in the ‘ideal’ range (>0.90 Reference Weller, Bowen and Faubert26 ), in contrast to the class 4 and 5 models, which included multiple values in the ‘acceptable’ range. Lo–Mendell–Rubin likelihood ratio test suggested that a three-class solution provided a better fit than a two-class, but no significant improvement for a four-class solution compared with a three-class.
The means of the three classes weighted by estimated class probability are shown in Fig. 1, and intercept and slope (T2–5 change) parameter estimates are displayed in Table 4. The Impact of Event Scale-Revised cut-off score is also overlain, and indicates the threshold of clinically significant moderate–severe post-traumatic stress symptom severity, rather than a formal PTSD diagnosis. The largest class, ‘resilient’ (n = 324; 74%), is characterised by mean symptom scores consistently below the cut-off and a general decreasing trend over time, except for a small increase at T4 that corresponded with an increase in COVID-19 hospital activity in Northern Ireland relative to T3. The ‘recovering’ class (n = 99; 23%) also displayed mean levels that fluctuated in line with the peaks and troughs of the pandemic. By contrast, the ‘recovering’ class started with mean levels well in excess of the cut-off, then displayed symptom reduction over time to the point that symptom levels were actually below the cut-off by T5. The ‘chronic’ class (n = 16; 4%) was the smallest, starting with mean levels above the cut-off; however, their symptoms generally worsened over time and were at their highest by T5. In contrast to the other two classes, the ‘chronic’ class self-reported greater symptoms when COVID-19 hospital activity was at its lowest: specifically T3 and, in particular, T5.
Trajectories of the three-class, post-traumatic stress, latent class growth analysis model. T, time. Time is scaled to account for unequal intervals of measurement.

Mean scores for growth parameters in the three-class unconditional model

*P < 0.05, ***P < 0.001.
Prediction of class membership
Table 5 displays the logistic coefficients for predictors of post-traumatic stress symptom class membership.
Multinomial logistic regression of predictors of class membership

*P < 0.05, **P < 0.01, ***P < 0.001.
A higher level of exposure to COVID-19, having at least one COVID-19 risk factor and perceiving communication from their organisation to be less effective were all associated with a greater likelihood of being in the ‘recovering’ or ‘chronic’ classes than the ‘resilient’ class.
Discussion
In line with the study hypotheses, three distinct post-traumatic stress symptom trajectories were observed in healthcare workers during and after the COVID-19 pandemic, specifically (a) a ‘resilient’ group (74%), with symptoms that remained below cut-offs throughout the pandemic; (b) a ‘recovering’ group (23%), with initially moderate-to-severe post-traumatic stress symptoms that decreased to subthreshold levels post pandemic; and (c) a ‘chronic’ group (4%), with initially moderate-to-severe levels of symptoms that continued to increase post-pandemic. These trajectories highlight that the developmental course for the majority of healthcare staff (96%) was either not to develop clinically significant levels of post-traumatic stress symptoms or to recover by the end of the pandemic. Worryingly, however, for the small percentage of individuals (4%) who experienced chronic levels of clinically significant post-traumatic stress, their symptom trajectory implies a continued escalation of symptoms post-pandemic, making them an at-risk group. In contrast to wider anecdotal evidence, a number of factors concerned with healthcare organisational exposure to trauma and availability of staff supports (e.g. managing patients with COVID-19, usage of staff supports) were found not to be related to class membership. Instead, perception of poor communication within the healthcare organisation emerged as the only employment-related factor that predicted membership of both post-traumatic stress symptom classes. Other significant predictors were general risk variables such as increased exposure to COVID-19 in everyday life and personal health risk factors for COVID-19.
The present study found post-traumatic stress trajectories that were broadly similar to a number of those obtained in general population studies examining the longitudinal course of mental health symptoms across the pandemic, namely ‘resilient’, ‘recovering’ and ‘chronic’. Reference Shevlin, Butter, McBride, Murphy, Gibson-Miller and Hartman1,Reference López-Castro, Papini, Bauer, Swarbrick, Paul and Nizzi6 Although it is not possible to make direct comparisons among such studies, it nonetheless highlights the fact that the course of PTSD symptoms for healthcare workers seems to be somewhat comparable to that of the wider public, even though the proportions in each class may vary due to specific challenges faced by this population and methodological differences. This finding is perhaps unsurprising, given that healthcare staff differ primarily in regard to the severity of COVID-19 stressors and symptoms, whereas the fundamental patterns of change probably reflect universal traumatological processes that operate similarly across society, regardless of occupation.
The predominant role of poor communication in contributing to COVID-19-related PTSD among both the general public and healthcare staff has received support from several studies. Reference Jordan, Shannon, Browne, Carroll, Maguire and Kerrigan7,Reference Wang, Huang, Wang, Wang, Xu and Li27 The pandemic is widely acknowledged as having fostered a unique ‘internal communication crisis’ for workplaces. Reference Ruck and Men28 In such crisis circumstances, employees have both increased anxiety and a greater need for information from their host organisation, with provision of clear, frequent, transparent and useful information important for alleviating distress. Reference Ecklebe and Löffler29 However, the COVID-19 pandemic presented a number of obstacles that made such a response challenging, including the following: (a) healthcare managers were under considerable pressure, which reduced capacity for communication; (b) reliable and timely information was often unavailable due to the evolving nature of the pandemic; and (c) many established forms of communication were no longer viable (e.g. in-person team meetings), with a marked shift towards new methods (e.g. online meetings, newsletters) that proved to be a difficult adjustment. Reference Ecklebe and Löffler29,Reference Björk, Corin, Akerstrom, Jonsdottir, Innocenti and Wijk30
Although the precise role of communication in the trajectories observed requires further exploration, it is likely that challenges in this area, combined with staff members’ unmet need for information in order to alleviate distress, contributed to the trajectories of post-traumatic stress in the ‘chronic’ and ‘recovering’ groups. It is possible that a number of other psychological and employment variables played an important explanatory role, including moral distress and the availability of wider staff supports. Reference Dyer, Shannon, McCann, Mitchell, Kerrigan and McClements31,Reference Fattori, Comotti, Mazzaracca, Consonni, Bordini and Colombo32 For example, the present study identified that individual risk factors, such as increased exposure to COVID-19 in everyday life and personal health risks for COVID-19, were predictive of post-traumatic stress trajectories, suggesting that some individuals may have held perceptions of being more physically vulnerable when working during the pandemic, or of posing a heightened risk of infecting others. Such psychological appraisals would be in line with established PTSD risk factors, including vulnerability beliefs and moral distress. Reference Cai, Ding, Tang, Wu and Yang33,Reference Norman, Feingold, Kaye-Kauderer, Kaplan, Hurtado and Kachadourian34
Clinical implications
Long-term post-pandemic PTSD monitoring and support for the healthcare workforce may be warranted based on the symptom courses observed in the present study, with future studies repeating such symptom measurements. The ‘chronic’ group trajectory indicates an escalating trend in post-traumatic stress reactions at the final post-pandemic assessment point. It is therefore quite possible that these symptoms in the at-risk group may have continued or even become exacerbated as time passed, as documented in other longitudinal studies of chronic PTSD. Reference Solomon, Mikulincer, Ohry and Ginzburg5 Indeed, the 3 post-pandemic trajectories obtained in the present study evolved from 2 peri-pandemic trajectories (i.e. low-symptom class (84%), high-symptom class (16%)) obtained using the same population. Reference Jordan, Shannon, Browne, Carroll, Maguire and Kerrigan9 Considering the inherent variability over the course of PTSD, Reference Bonanno, Mancini, Horton, Powell, LeardMann and Boyko35 it could be the case that the other groups, most notably the ‘recovering’ group, had later deviations in their trajectory, potentially returning to clinically significant levels. A number of general PTSD studies have noted significant changes in symptom trajectories in longer-term, 7- to 29-year follow-ups. Reference Solomon, Mikulincer, Ohry and Ginzburg5,Reference Bonanno, Mancini, Horton, Powell, LeardMann and Boyko35 The risk of relapse is also potentially heightened following further exposure to similar traumatic experiences. Reference Solomon36 This paints a concerning picture for healthcare workers in the ‘chronic’ and ‘recovering’ groups, who are now likely to be working in healthcare systems under significant pressure following the pandemic (e.g. waiting list backlogs, severe staff shortages 37 ) and may be vulnerable to further aversive professional experiences and moral distress. Reference Jeffs, Heeney, Johnstone, Hunter, Loftus and Ginty38 Consequently, it is vital for the continued analysis of these trajectories, as well as both the empirical validation and development of PTSD interventions for healthcare staff, to offset mental health deterioration in these potentially reactivating circumstances.
The study also has implications for understanding PTSD following societal trauma, in providing support to the relatively consistent finding of ‘resilient’, ‘recovering’ and ‘chronic’ trajectories across populations and following large-scale traumatic events (e.g. warfare, pandemics). Ineffective internal and external communication has also been highlighted as a significant issue in regard to pandemic responses, particularly COVID-19. Reference Ruck and Men28 The present findings highlight the fact that internal communication can be particularly impactful on the well-being of healthcare staff during and after pandemic response, and therefore effective communications strategies should be researched and developed in the advent of future crises to ensure that the well-being of healthcare staff is preserved.
Limitations
The present study has a number of strengths and weaknesses. A key strength is the longitudinal design permitting both peri- and post-pandemic follow-up on the course of post-traumatic stress symptoms in a substantive sample of healthcare workers. However, despite the methodological rigour, a number of important analyses were not supportable due to sampling limitations: most notably, the small size of the ‘chronic’ group (4%) meant that statistical power was not available to examine predictors of ‘chronic’ versus ‘recovery’ classes, which would have been useful in understanding the mechanisms of trajectory development. However, because the effects of COVID-19 exposure, COVID-19 risk factors and perceived ineffective communication were consistently stronger for the ‘chronic’ compared with the ‘recovering’ group (when compared with the ‘resilient’ class), it is likely that there is some degree of a dose–response relationship between risk factors and likelihood of being in the ‘chronic’ class. Caution should be exercised when interpreting the predictors of class membership analysis, because statistical power is limited for the gender variable given that males comprise only a small proportion of the sample. As a non-probability sample was used, self-selection bias also cannot be ruled out and it is possible that those with mental health issues were more likely to have engaged in the study across and after the pandemic, Reference Eisma39 which limits sample representativeness and potentially inflates prevalence. It should also be acknowledged that, although missing data levels at T5 were high, evidence suggests that when multiple imputation models are carefully specified (e.g. auxiliary variables are factored in), as was done in the present study, estimation remains accurate. Reference Austin and van Buuren40 Longitudinal pre-pandemic data on the sample were not available and, in order to minimise survey fatigue, it was possible to carry out only one post-pandemic follow-up (i.e. T5). Therefore, the potential impact of prior mental health on likelihood of class membership and the exact shape of change in trajectories between August 2021 and June 2023 could not be examined.
Exposure to additional traumatic events during the peri- and post-pandemic periods was also not assessed in the present study, and thus their potential contribution to longer-term symptom trajectories cannot be entirely ruled out. However, the largely representative nature of the staff sample – and the fact that most demographic and occupational variables did not predict group membership – reduce the likelihood that trauma exposures were disproportionately or systematically concentrated within a single latent trajectory group (e.g. the ‘chronic’ group). Last, whereas brief self-report measures are appropriate for screening post-traumatic stress symptomatology, their use in determining clinical PTSD levels has been shown to inflate caseness relative to clinical interviews, due to their high sensitivity and low specificity. Reference Scott, Stevelink SA, Gafoor, Lamb, Carr and Bakolis8 Nevertheless, conservative interpretations of clinically significant symptom changes can still be made using these instruments, because repeated clinical interviews for PTSD symptoms across five time points would be extremely challenging within the current longitudinal design.
To conclude, three post-traumatic stress symptom trajectories were observed in healthcare workers during and after the COVID-19 pandemic, namely ‘resilient’, ‘recovering’ and ‘chronic’. Although these groups are similar to classes found in the wider population, the escalating trajectory in the ‘chronic’ group and the unique, demanding nature of working in healthcare settings highlight the need for monitoring staff in healthcare services for PTSD symptoms after the pandemic, and offering the workforce psychological support. Prevention research should focus on internal communication enhancement and strategies within healthcare systems, in order to develop robust protective processes in the advent of future health crises.
Data availability
The data that support the findings of this study are available from the corresponding author (K.F.W.D.), upon reasonable request. Demographic data will be aggregated in order to protect the identity of the participants.
Author contributions
K.F.W.D., J.-A.J. and C.S. formulated the research questions, secured funding for the study, designed and carried out the study, interpreted the findings and drafted the manuscript. The analyses were undertaken by J.-A.J. N.H. interpreted the findings and led in drafting sections of the manuscript. The final manuscript was approved by all authors prior to submission.
Funding
This work was supported by the Health and Social Care Research & Development Division, Public Health Agency (grant no. COM/5602/20).
Declaration of interest
None.






eLetters
No eLetters have been published for this article.