Next Article in Journal
Contraceptive Behavior in Appalachia: Exploring Use, Nonuse, and Contraceptive Attitudes
Next Article in Special Issue
Attention Deficits in Healthcare Workers with Non-Clinical Burnout: An Exploratory Investigation
Previous Article in Journal
Health Disparities in the Use of Primary Cesarean Delivery among Asian American Women
Previous Article in Special Issue
Mental Health Challenges during COVID-19 Pandemic: Experiences of Primary Healthcare Nurses in Durban, South Africa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trajectories of Adjustment Disorder and Well-Being in Austria and Croatia during 20 Months of the COVID-19 Pandemic

by
Irina Zrnić Novaković
1,*,
Alina Streicher
1,
Dean Ajduković
2,
Marina Ajduković
3,
Jana Kiralj Lacković
2,
Annett Lotzin
4,5 and
Brigitte Lueger-Schuster
1
1
Department of Clinical and Health Psychology, Faculty of Psychology, University of Vienna, 1010 Vienna, Austria
2
Department of Psychology, Faculty of Humanities and Social Sciences, University of Zagreb, 10000 Zagreb, Croatia
3
Department of Social Work, Faculty of Law, University of Zagreb, 10000 Zagreb, Croatia
4
Department of Psychiatry and Psychotherapy, University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany
5
Department of Psychology, MSH Medical School Hamburg, 20457 Hamburg, Germany
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(19), 6861; https://doi.org/10.3390/ijerph20196861
Submission received: 31 August 2023 / Revised: 26 September 2023 / Accepted: 27 September 2023 / Published: 29 September 2023

Abstract

:
The present study aimed to investigate the trajectories of adjustment disorder (AD) symptoms and well-being over 20 months of the COVID-19 pandemic in Austria and Croatia. Further objectives of this study were to examine whether sociodemographic characteristics and the symptoms of anxiety and depression could predict these trajectories. As part of the pan-European ESTSS ADJUST study, N = 1144 individuals were recruited using convenience sampling and assessed four times between June 2020 and January 2022 through an online survey. Latent growth curve modelling was applied to estimate the trajectories of AD symptoms and well-being. Over time, the prevalence of probable AD varied between 9.8% and 15.1%. The symptoms of AD tended to increase, whereas well-being tended to decrease. According to the majority of the models tested, women, participants from Austria and those with lower income had higher initial AD symptoms, whereas older participants and those from Croatia had higher initial well-being. In all models and at all timepoints, anxiety and depression significantly predicted AD and well-being scores. Overall, our study points to several predictors of AD and well-being and indicates high variability in people’s reactions to the pandemic. Psychosocial support for the general population is needed during pandemics and similar crises, with a special focus on vulnerable groups.

1. Introduction

The past three years have been marked by the COVID-19 pandemic and its repercussions. Compromised mental and physical health, disrupted routines and changed social dynamics, along with financial losses and work-related challenges, are just some of the stressors people have been repeatedly exposed to (for an overview of stressors related to the COVID-19 pandemic, see [1]). Due to its long-lasting and potentially traumatic impact [2], COVID-19 has been described as “continuous traumatic stress” [3]. According to the literature, adjustment disorder (AD) can best capture a maladaptive reaction to pandemic-related stress or trauma [4,5]. However, longitudinal studies on changes in AD symptoms over time are scarce, and potential cross-country differences amidst COVID-19 are still poorly understood. Acknowledging that AD symptoms are likely to develop in response to pandemics but also to other global challenges (e.g., the climate crisis, inflation, wars and refugee crises), the present study aims to better understand their trajectories and underlying factors in Austria and Croatia. Despite their geographical proximity, historical ties and cultural and economic exchanges [6], these two European Union countries differ in terms of their socioeconomic situations and degrees of life satisfaction [7,8], as well as their responses to the COVID-19 pandemic. For context, Croatia implemented stricter measures at the beginning of the pandemic, whereas government measures in Austria were stricter in the later stages of the pandemic [9].

1.1. Adjustment Disorder Amidst COVID-19

AD can be defined as a maladjustment to a stressful event or severe changes in life, such as job loss, financial strain or the illness of a loved one [10,11]. It is characterised by a preoccupation with the stressor and a failure to adapt which severely impacts one’s functioning. According to the World Health Organisation [12], the symptoms of AD usually resolve within six months. However, the symptoms might last longer if the stressor persists longer [13], which was the case for the COVID-19 pandemic and its consequences.
During the pandemic, international studies researching different samples reported prevalence rates of AD ranging from 8.2% to 61.3% [14,15]. In older people, diagnoses of AD were more often documented before (2017–2019) than during the pandemic (2020 [16]). Gender differences were also found, with 40.7% of women compared to 31% of men scoring above cut-off for probable AD in a large sample of Georgian participants [17].
Apart from the above-mentioned sociodemographic characteristics, adjustment difficulties have been linked to COVID-related stressors, such as social restrictions, governmental pandemic management, work and financial difficulties [18]. Notably, the mental health impact of these stressors varied between countries, with countries having had more experience with stress and trauma showing better outcomes [19]. Different pandemic loads in terms of the stringency index and incidence and death rates could also explain cross-country differences in some mental health outcomes, e.g., anxiety, depression and post-traumatic stress disorders [20,21]. However, associations between pandemic load and AD symptoms have not yet been investigated.
In general, some findings on factors associated with AD symptoms are inconsistent, whereas others are more robust. For instance, older participants were found to be at risk for developing AD in some studies [18], while a protective role of older age was proposed in other studies [22]. In contrast, a large body of research has shown that AD symptoms are associated with different emotional responses to the pandemic, including the fear of infection [23], COVID-19 concerns [17] and depression symptoms [24].

1.2. Trajectories of Adjustment Disorder Amidst COVID-19

To date, only a small number of studies have investigated changes in AD symptoms over the course of the COVID-19 pandemic. There is evidence that the symptoms of AD consolidated in autumn 2020 [25], although a German study comparing AD symptoms in the spring and autumn of 2020 found no significant differences [26].
A representative study conducted in Israel in the summer of 2020 (before the second lockdown) and in late autumn in 2020 (after the second lockdown) detected four trajectories of AD symptoms, suggesting different patterns of adjustment to the pandemic [27]. Remarkably, there is only limited knowledge of positive psychological adjustments to the pandemic over time since the indicators of positive mental health, such as well-being, have not yet been included in the research on AD amidst COVID-19. Thus far, scientific attention has mainly been directed toward vulnerable groups and risk factors for probable AD (e.g., [23,28]).
Shifting the focus toward indicators of positive mental health and the course of AD symptoms in the later phases of COVID-19 could yield a more differentiated picture of the psychological adjustment to the pandemic. Longitudinal studies over a longer period of time are particularly beneficial given that some individuals may experience adjustment difficulties only in the later phases of the pandemic [28]. Moreover, the changing context of the pandemic in terms of the strictness of COVID-19 restrictions and the incidence and death rates might foster an understanding of changes in both adjustment and well-being, yet these interrelations have not been explored.

1.3. Current Study

In summary, the existing research on AD amidst COVID-19 has some limitations. Firstly, given that the majority of studies were conducted in the first months of the pandemic and covered intervals of several weeks or a few months, little is known about the trajectories of AD symptoms over a longer period and in later phases of the pandemic. Secondly, the co-development of well-being and AD symptoms during the pandemic is still poorly understood. The inclusion of sociodemographic characteristics (e.g., income and education) and mental health indicators (e.g., anxiety symptoms) might increase knowledge about this co-development [29,30,31]. Finally, a significant amount of uncertainty still exists with respect to cross-country differences, and a closer examination of the course of the pandemic in different countries might help to explain the development and severity of adjustment difficulties.
To address these research gaps, the present study examined the trajectories of AD symptoms in Austria and Croatia over 20 months of the COVID-19 pandemic. Both countries are part of the pan-European ADJUST study on stressors, coping and the symptoms of AD [18,32]. In addition to the three waves of assessments conducted in all ADJUST member countries, a fourth wave was launched in Austria and Croatia in December 2021–January 2022. The aim was to assess psychological reactions to the pandemic one year after the large spike in COVID-19 cases and deaths over the winter months of 2020/2021 [33,34].
Several factors provide a basis for comparing Austrian and Croatian participants. Firstly, life satisfaction and quality of life have been higher in Austria than in Croatia in previous years, which could largely be explained by the higher gross domestic product (GDP) per capita among Austrians [35,36]. Secondly, the countries differ in terms of mental health, with Austrian citizens experiencing a higher mental health burden [37]. Thirdly, Austria and Croatia have been affected by COVID-19 to different extents, with, e.g., higher numbers of confirmed deaths being reported in Croatia in the second year of the pandemic [38]. Government policies for managing COVID-19 also differed between the countries [34]. For these reasons, a comparison of Croatian and Austrian data might shed light on the effects of a variety of socioeconomic and contextual factors on adjustment trajectories.

1.4. Aims and Research Questions

The present study set out to explore the trajectories of AD and well-being in the first two years of the COVID-19 pandemic using latent growth curve modelling. Our main aims were to answer the following research questions:
  • How do the trajectories of AD symptoms and well-being develop over time in Austria and Croatia?
  • To what extent can sociodemographic characteristics (country, age, gender, income and education) predict the trajectories of AD and well-being?
  • To what extent can the symptoms of anxiety and depression predict the trajectories of AD and well-being?
Acknowledging that people are capable of coping with stressful situations and adapting to changing situational demands [39,40], we hypothesised that the symptoms of AD would decrease, whereas well-being would increase over time. Based on the existing evidence, we further hypothesised that sociodemographic characteristics and the symptoms of anxiety and depression would predict the trajectories of AD and well-being.

2. Materials and Methods

The present study is part of the international ADJUST COVID-19 study involving 11 European countries [32]. For further information, please see the study’s website: https://estss.org/adjust/ (accessed on 26 September 2023) and the corresponding OSF registry where the study was preregistered: https://doi.org/10.17605/OSF.IO/8XHYG (accessed on 26 September 2023). The data used for the secondary analysis described in this paper were collected from the general populations of Austria and Croatia at four measurement timepoints. As described in the introduction, only these two countries launched a fourth assessment wave in the course of the ADJUST study. Though they are closely related, e.g., in terms of cultural and economic exchanges [6], Austria and Croatia differ in several aspects [35,38], thus forming an interesting basis for comparison. The periods of data collection and information about the COVID-19 situation in both countries at each timepoint are presented in Table 1.

2.1. Procedure and Participants

Over the course of the ADJUST study (i.e., the primary analysis), all assessment waves were conducted as online surveys using LimeSurvey, an online tool for creating anonymous surveys and questionnaires in line with the General Data Protection Regulation (GDPR; www.limesurvey.org, accessed on 31 August 2023). We used several strategies to recruit participants for the first assessment wave (T1). For instance, our study invitation was distributed among the staff of major companies and professional associations, as well as among members of numerous leisure and interest groups. We also shared the study link on social media and distributed flyers in public spaces. A detailed overview of the recruitment strategies can be found at www.ncbi.nlm.nih.gov/pmc/articles/PMC8725769/bin/ZEPT_A_1964197_SM8752.docx (see Supplement 2, accessed on 26 September 2023). The survey was accessible via a website link or a QR code. On the first page of the survey, potential participants could read important information about the study (the aims of the study, data management and the right to withdraw) and, if interested, provide their informed consent.
To participate in the present study (i.e., the secondary analysis), the following criteria were required: (1) the participants needed to be at least 18 years old, (2) have good command of the German or Croatian language, (3) live in Austria or Croatia, (4) agree to participate in the study, and (5) respond to at least two assessment waves. In Austria, ethical approval for the study was granted by the Ethics Committee of the University of Vienna (date of approval: 29 May 2020, reference number: 00554), whereas the Ethics Committee of the Department of Psychology in the Faculty of Humanities and Social Sciences at the University of Zagreb granted ethical approval in Croatia (date of approval: 21 May 2020). Prior to the initiation of the survey, informed consent was obtained online from all participants.
At baseline (T1), N = 1001 valid cases were obtained in Austria and N = 1787 were obtained in Croatia, out of which 52.5% Austrian participants (n = 525) and 57.8% Croatian participants (n = 1032) agreed to take part in the following waves. These participants were re-contacted after approximately 6, 12 and 18 months and asked to participate in further study waves (T2, T3 or T4, respectively).
In total, N = 1144 participants (Mage = 44.0, SDage = 13.5), who were predominantly female (73.5%) met the inclusion criteria and were included in the present analysis. Table 2 shows important characteristics of the study participants.

2.2. Measures

An overview of the measures included in the ADJUST study can be found in the respective study protocol [32]. For the present analysis, the following measures were used in the German and Croatian language, respectively.
The 8-item Adjustment Disorder New Module (ADNM-8 [41]) screens for life stressors and AD symptoms according to the 11th version of the International Classification of Diseases (ICD-11 [12]). It consists of four items on preoccupation and four items on adjustment difficulties; each are scored on a 4-point Likert scale ranging from 1 (never) to 4 (often). Higher scores indicate more intense AD symptoms, and a score of >22 represents a cut-off for probable AD. The ADNM-8 has shown good psychometric properties [41,42] and has been used to assess AD symptoms during the pandemic [28,43].
The 5-item World Health Organization Well-Being Index (WHO-5 [44]) measures subjective well-being on a 6-point Likert scale ranging from 0 (at no time) to 5 (all of the time). The total raw score, with possible values between 0 and 25, is multiplied by 4, resulting in a final score between 0 and 100, with the latter representing the highest degree of well-being. A score below 51 suggests poor well-being and might be indicative of depressive symptomatology [45,46]. The high validity and applicability of the WHO-5 have been demonstrated across different countries and study samples [46], and the scale has been widely used in the pandemic context [20,47].
The 4-item Patient Health Questionnaire (PHQ-4 [48]) measures the core symptoms of depression and anxiety on a 4-point Likert scale (ranging from 0 = not at all to 4 = more than half the days). It comprises a 2-item depression scale (the 2-item Patient Health Questionnaire (PHQ-2)) and a 2-item anxiety scale (the 2-item Generalized Anxiety Disorder Scale (GAD-2)). A score of >2 on one of these scales is used as cut-off point for possible depression or anxiety. All three scales (i.e., the PHQ-4, PHQ-2 and GAD-2) have shown good reliability and validity [49] and have been applied in various COVID-19 studies [50].
In the present study, each scale had a high level of internal consistency across all measurement timepoints, ranging from 0.81 to 0.92 (for details, please see the Supplementary Materials, Table S1). Depending on the measurement, Cronbach’s α or the Spearman–Brown coefficient was used, as appropriate [51]. In addition to the described scales, self-constructed items were used to assess the sociodemographic characteristics of the participants. Gender, age, education and income were included in the main analysis as these variables were found to be predictive of AD and/or well-being in previous studies [18,52,53,54] and might explain the differences between mental health outcomes in Austria and Croatia.

2.3. Data Analysis

The analyses were conducted using IBM SPSS Statistics (version 27 [55]) and RStudio (version 2022.12.0.353 [56]). After performing descriptive and missing data analyses (using the Mann–Whitney U test and the Welch t-test), measurement invariance was assessed. Cochran’s Q test, followed by a post hoc analysis using Dunn’s procedure, was conducted to explore whether the percentage of participants at risk for AD, depression and anxiety was different at the different timepoints. The complete dataset without missing values was used for this analysis (n = 531). In addition, we conducted a linear regression analysis to check whether the anxiety and depression scores attained at T1 could predict AD symptoms and well-being scores at T4. With this exploratory analysis, we wanted to indirectly account for the possible predictive value of pre-existing mental health symptoms since no reliable pre-pandemic data were available in the study.
Following these preparatory steps, latent growth curve models (LGCMs) of different levels of complexity were tested. General rules for evaluating model fit were used for all analyses [57,58].

2.3.1. Measurement Invariance

Measurement invariance considers the psychometric equivalence of a construct over time or across different groups [59]. The assessment of invariance is an important step prior to testing group differences (e.g., cross-country differences, as in this study) or comparing different measurement occasions (e.g., changes over time, as in this study). Therefore, we assessed invariance over populations by comparing the Austrian and Croatian samples, and we assessed longitudinal measurement invariance by comparing the data collected at each of the four measurement timepoints.
We tested measurement invariance in the structural equation modelling (SEM) framework using a confirmatory factor analysis (CFA [60]) for all relevant measures (i.e., the ADNM-8, WHO-5, PHQ-2 and GAD-2). Following the general guidelines regarding estimation methods [60,61], we used the weighted least square mean and variance adjusted (WLSMV) estimator for testing the invariance of measures with four or fewer answering options (i.e., the ADNM-8, PHQ-2 and GAD-2) and the robust maximum likelihood (MLR) as an estimator for testing the invariance of measures with more than four answering options (i.e., the WHO-5). To handle missing data, listwise deletion was used in combination with the WLSMV, and the Full Information Maximum Likelihood (FIML) in combination with the MLR [60]. Parcelling was applied for the analysis of the ADNM-8 as it was the only measure with more than five items.
Given that all measures used in the study are well-established, we assumed psychometric equivalence across groups and over time. Therefore, we applied the model-building approach (also known as the constrained baseline approach) to assess measurement invariance, using the most restricted model (i.e., strong invariance) as the baseline [60]. This analysis was performed using lavaan, an open source R package (version 0.6.12 [62]).

2.3.2. Latent Growth Modelling

Latent growth curve (LGC) modelling is a statistical method through which within-person changes (i.e., intra-individual patterns of change) and between-person differences (i.e., inter-individual variability) over time can be estimated [63,64,65]. Thus, LGC modelling allows one to examine individual differences in a longitudinal design (i.e., individual latent trajectories), which is not feasible using conventional methods, such as an ANOVA. Another advantage of LGC modelling is its ability to handle missing data effectively [66].
The essential elements of all LGCMs are fixed and random effects. Fixed effects are factors estimated from the data, namely the mean intercept (i.e., starting point) and the mean slope (i.e., the increment or rate of change). Jointly, these factors define the average latent trajectory for the entire sample. On the other hand, random effects capture the variance around the average trajectory, i.e., the differences in the starting point and the rate of change across individuals. Taken together, the fixed and random effects allow for the description of the general characteristics of the growth curve (i.e., trajectory) for the entire sample and the subgroups within this sample [63,64].
In the present study, LGC modelling was performed within the SEM framework [67], using lavaan (version 0.6.15 [62]). Given that a linear change was hypothesised and the four timepoints were approximately equally spaced, the following coding scheme was used to represent the slope: 0, 1, 2, 3. For the intercept, the factor loadings for four measurement timepoints were set to 1 (i.e., 1, 1, 1, 1). In line with the current recommendations [64,68], the MLR was the preferred estimator for all analyses. A state-of-the-art method of handling missing data, FIML, was used to deal with partially missing data [69].
A series of LGCMs were used to assess changes in AD symptoms (determined via the ADNM-8) and well-being (determined via the WHO-5) over time:
  • Univariate models: Firstly, univariate unconditional models were estimated to explore changes over time in AD and well-being separately. Next, the time-invariant covariates (TICs [64]) age, gender, country, education and income were included as predictors of the variability in the intercept and slope in two conditional models for AD symptoms and well-being, respectively. In a third step, depression (determined via the PHQ-2) and anxiety (determined via the GAD-2) scores were added as time-varying covariates (TVCs [64]).
  • Multivariate models: The same procedure was repeated for multivariate models (i.e., joint models for AD symptoms and well-being) to investigate how AD symptoms and well-being simultaneously unfold over time [65]. The final model included all TICs and TVCs to explore whether these can explain the residual variance in the joint model. Moreover, the covariance between the intercept and slope was estimated for each primary outcome. Residuals between the primary outcomes were allowed to correlate within one timepoint. Figure 1 shows the proposed multivariate model of AD symptoms and well-being with all predictors and covariances.
  • Multivariate models by country: In the last step, the final model was executed in the Austrian and Croatian samples separately in order to explore possible differences in the growth trajectory and predictors between the two countries. For this final analysis, the predictor “country” was obsolete; thus, it was excluded from the analysis. The other model specifications remained unchanged.

3. Results

3.1. Participant Flow and Missing Data

Participants who did not meet the inclusion criteria and those who did not provide answers for the primary outcome measures were excluded from the analysis, resulting in a total of N = 1144 participants (nAustria = 415; nCroatia = 729). Detailed information about participant flow is provided in the Supplementary Materials (Table S2).
An analysis of missing data was conducted by comparing the participants who responded to all the assessment waves (n = 531) with those who had missing data for one or more assessment waves (n = 613). The two groups did not differ in terms of gender, p = 0.493. However, the participants who responded to all four assessment waves were significantly older (M = 45.01; SD = 13.97) than those who responded to two or three waves (M = 43.20; SD = 13.05); t(1142) = -2.27 and p = 0.023. According to a Mann–Whitney U test, the participants who responded to all waves were slightly better educated (mean rank = 597.48) than the non-respondents (mean rank = 550.87); U = 176,013.50, z = 3.08 and p = 0.002. Moreover, the average monthly household income was higher in respondents (mean rank = 577.13) than in non-respondents (mean rank = 539.09); U = 163,964.50, z = 2.052 and p = 0.040. Please note that complete cases were used for the descriptive analyses, whereas the FIML was used for the LGCMs.

3.2. Cross-Country Differences and Measurement Invariance

Due to a violation of the homoscedasticity assumption, a Welch t-test was run to determine whether there were differences in the mean ages of the Austrian and Croatian participants. The Austrian participants were significantly older than the Croatian participants (M = 3.78, 95% CI [2.09, 5.47], t(755.11) = 4.40, p < 0.001). The difference in the proportions of male and female participants in the two countries was also significant (p = 0.007). Furthermore, the participants from Croatia reported higher levels of educational attainment (U = 184606.50, z = 8.04, p < 0.001), whereas the participants from Austria reported higher income (U = 124867.00, z = −3.05, p = 0.002).
Regarding measurement invariance, the model fit was satisfactory in the most restricted model for all measures tested (see the Supplementary Materials, Table S3), suggesting strong invariance for the ADNM-8, WHO-5, PHQ-2 and GAD-2. Thus, the results show that all the measures included in the analysis measured the same constructs over different timepoints and different groups (i.e., the Austrian and Croatian samples).

3.3. Mental Health Outcomes over Time

At T1, one in ten participants reported clinically relevant AD symptoms, and approximately 15% met the criteria for probable depression or an anxiety disorder. All mental health outcomes assessed in the study showed changes over time, with the highest AD and depression scores at T2 and the highest anxiety scores at baseline. At each timepoint, a considerable number of participants reported reduced well-being, with percentages ranging from 33.7% to 45.1%. Descriptive statistics of all relevant measures at each timepoint are shown in Table 3. The observed individual trajectories of the ANDM-8 and WHO-5 scores are provided in the Supplementary Materials (Figures S1–S4).
The percentages of self-reported AD differed significantly over time (χ2(3) = 22.71; p < 0.001). The highest percentage of participants at risk of developing AD was at T2 (15.4%; n = 82), which was significantly more than at T1 (9.2%; n = 49), T3 (10.4%; n = 55) and T4 (11.5%; n = 61), p < 0.001, p = 0.002 and p = 0.027, respectively.
A significant difference regarding probable depression was also found (χ2(3) = 8.80; p = 0.032). The percentage of participants at risk for probable depression was equally high at T2 and T4 (16.0%; n = 85), followed by T1 (13.9%; n = 74) and T3 (11.5%; n = 61), though Dunn’s procedure with a Bonferroni correction could not detect where the differences in proportions lay (p ≥ 0.065). The percentage of participants at risk for probable anxiety changed over time (χ2(3) = 11.38; p = 0.010), with the highest prevalence rates found at T1 and T4 (14.9%; n = 76). These estimated prevalence rates were significantly different from those found at T3 (9.6%; n = 51), p = 0.023, but not from those found at T2 (13.4%; n = 71), p > 0.050.
Linear regression analyses showed that both depression and anxiety scores at baseline could predict AD scores at the last timepoint (F(1, 672) = 110.72; p < 0.001 and F(1, 672) = 122.13; p < 0.001, respectively). The adjusted R2 was 14.0% for depression and 15.3% for anxiety, both indicating a moderate effect size according to Cohen [70].
The depression and anxiety scores at baseline were also predictive of the WHO-5 scores at the last timepoint (F(1, 672) = 165.20; p < 0.001 and F(1, 672) = 150.66; p < 0.001, respectively). The effect size was moderate, with R2 = 19.7% for depression and R2 = 18.3% for anxiety [70].

3.4. Univariate LGCMs

The model fit indices for the univariate LGCMs showed mixed results. While the Comparative Fit Index (CFI) and Tucker–Lewis Index (TLI) mostly indicated an acceptable fit, the Root Mean Square Error of Approximation (RMSEA) and Standardised Root Mean Square Residual (SRMR) occasionally indicated a poor fit. All fit indices and parameter estimates are presented in Table 4.
In all models, the intercept was significant, whereas the slope was significant only in the unconditional model and the model with TVCs. According to these two models, AD symptoms increased, whereas well-being scores decreased over time. The covariance between the slope and intercept did not reach statistical significance in any of the models tested. The same was true for the predictors of the slope, with only country being a nearly significant predictor of the change in well-being over time in the second model (Table 4).
On the other hand, several intercept predictors were found. Women, participants from Austria and those with lower income showed, on average, higher baseline scores on the ADNM-8, whereas older participants, participants from Croatia and those with higher income showed higher levels of well-being at baseline.
Regarding the time-varying predictors, both the depression and anxiety scores were predictive at all timepoints, with slightly higher estimates for anxiety in AD models and higher estimates for depression in well-being models. From the empirical–practical perspective, the partially poor fit indices and the significant variances were sufficient to justify testing additional models.
Figure 2 shows the estimated trajectories of AD and well-being based on the unconditional LGCM. For better clarity, individual trajectories are depicted only for a randomly selected subsample; the remaining individual trajectories are provided in the Supplementary Materials (Figures S5–S8).

3.5. Multivariate LGCMs

The multivariate models achieved a better fit with the data than the univariate models. The intercepts for both AD and well-being reached statistical significance in the unconditional and conditional model, as presented in Table 5. The covariance between the intercepts of both constructs was significant in both models, indicating that on average, participants with higher initial AD symptoms showed lower initial levels of well-being and vice versa. The significant covariance between the slopes of both constructs suggests a negative association between the rate of change in AD symptoms and the rate of change in well-being. However, this covariance was only significant in the unconditional model. The covariances of the residuals were significant at each timepoint (p < 0.001).
Several significant predictors were detected in the final conditional model with all TVCs and TICs (see Table 5 for details). Higher symptoms of AD at baseline were found in women, older participants, those coming from Austria and those with higher income, with the last one being only nearly significant. Initial well-being scores were found to be higher in older participants and those coming from Croatia. None of the slope estimates reached statistical significance. Regarding the TVCs, both anxiety and depression scores could predict AD and well-being scores at all timepoints. Again, the predictive power of the anxiety scores was higher in the case of AD, whereas the predictive power of the depression scores was higher in the case of well-being. The intercepts of both measures and the slope of AD showed significant variance in both multivariate models but were lower in the conditional model.

3.6. Multivariate LGCMs by Country

The proposed linear model showed no convergence issues in Austria. However, when estimating the model in the Croatian sample, a negative residual variance for the well-being slope was detected. This is not plausible, given that computation of residual variances requires squaring values which are always positive [71]. We addressed the problem by freeing two initially fixed well-being parameters as follows: 0, a, b, 1 (instead of: 0, 1, 2, 3).
The model fit was acceptable for the Austrian and the Croatian data, as shown in Table 6. The intercept for the AD scores was slightly higher in Austria, whereas the intercept for the well-being score was higher in Croatia. The slope for the AD score was positive in the Austrian model, indicating an increase over time, while the negative slope in the Croatian model indicates a decrease over time. A decrease in the level of well-being was found in both models. The intercepts of AD and well-being were significantly negatively correlated in the Austrian model, suggesting a negative association between the ANDM-8 and WHO-5 scores at baseline. None of the tested covariances were significant in the Croatian model.
While age and income were significant predictors of the AD intercept in the Austrian model, gender was the only intercept predictor in the Croatian model. None of the TICs were significant predictors of the change over time regarding neither AD nor well-being.
Again, depression and anxiety were shown to be predictive of AD and well-being scores. Regarding AD, higher estimates for depression were found at T2 and T4 and higher estimates for anxiety were found at T1 and T3 in the model with Austrian participants. The opposite was found in the model with Croatian participants. Regarding well-being, the depression scores had a higher predictive power in both countries at all timepoints.

4. Discussion

The present study is the first study to investigate the trajectories of AD symptoms and well-being during nearly two years of the COVID-19 pandemic. We found an association between the starting point of the two trajectories, sociodemographic characteristics and the symptoms of depression and anxiety. Furthermore, the comparison between Austrian and Croatian data allowed us to identify country-specific predictors and differences in mental health trajectories.

4.1. Mental Health Outcomes over Time

Approximately one in ten participants met the criteria for probable AD at T1, T3 and T4. The prevalence was higher (15.1%) at the second timepoint, which was characterised by a notable increase in the number of cases (compared to T1) and stricter government measures than at other timepoints (see Table 1 for details). The intensity of the pandemic and associated restrictions have already been linked to mental health problems [20,47] but not specifically to AD. Previous studies conducted during COVID-19 reported prevalence rates of AD ranging from 8.2% [15] and 18.8% [72] to 61.3% [14]. However, these rates are only partly comparable with our results since the data were collected in the initial stages of the pandemic (spring of 2020), whereas our study started in the summer of 2020.
The prevalence rates of self-reported anxiety and depression ranged between 11.9% and 16.2%. This is relatively low compared to the pooled prevalence rates found in a recent umbrella review [21], which might be explained by the timing of the study. Namely, as shown in the review, prevalence rates were generally higher in studies conducted in the earlier stages of the pandemic, and the present study covered later stages of the pandemic.
The WHO-5 scores, which point to poor well-being and the presence of depressive symptoms in more than one-third of the participants at all timepoints, are rather alarming. This might indicate subtle mental health difficulties related to COVID-19 which do not necessarily lead to a diagnosis of a mental disorder but might impact a person’s mood and challenge their ability to cope. Moreover, several studies found a link between well-being and COVID-19 [31,47].

4.2. Trajectories of Adjustment and Well-Being

In the first step of the analysis, univariate models were tested to estimate the trajectories of AD and well-being separately. Contrary to our expectations, three univariate models for AD (i.e., the independent trajectories of AD) showed a minor increase in symptoms over time. Only one model showed a decrease; however, this result was not significant. The current literature on adjustment during COVID-19 is scarce and limited to the initial stages of the pandemic. For instance, a German study found a minor, non-significant decrease in AD scores in the first year of the pandemic [26]. The literature on AD outside the context of COVID-19 suggests that AD symptoms can persist for longer than six months if the stressor persists longer [13]. It is thus plausible that the duration and the changing course of the pandemic, with the varying number of infections and different restrictions being imposed at different times, repeatedly challenged people’s adaptational capacities, hindering successful adjustment.
The results regarding the trajectories of well-being support this assumption. According to the negative slopes in all the well-being models, well-being decreased over time, emphasising the difficulties in maintaining mental health over almost two years of the pandemic. Acknowledging the human ability to cope with and adapt to stressful situations [39], we hypothesised an improvement in mental health over time. However, COVID-19 posed numerous challenges in different areas of life [1], recurrently forcing people to deploy additional resources to cope and further exacerbating their mental fragility. This might explain the deterioration of mental health found in our study. A study including a post-pandemic follow-up would help to better understand adjustment in the aftermath of the COVID-19 pandemic.
The results of the multivariate models (i.e., the joint trajectories of AD and well-being) were only partially in accordance with those of above-described univariate models. The simultaneous estimation of the trajectories of AD and well-being confirmed a decrease in well-being over time. However, the slope of AD did not reach statistical significance, suggesting a flat trajectory. This is contradictory to the rising trajectory identified in the univariate models. At this point, it must be mentioned that all the slope means in the univariate models were rather low, indicating that the estimated trajectories were only slightly steep. It is possible that the slope means in the multivariate models were not any more significant due to the increasing complexity of the models (e.g., due to the additional covariances included in the models). Future studies should also test models of different complexities but with some additional predictors to better understand the factors underlying the joint trajectories of AD and well-being.
One noteworthy result of the multivariate models is the significant negative relationship between the initial symptoms of AD and the initial well-being scores. Thus, people with adjustment difficulties were more likely to have poor well-being at baseline, echoing the results of previous studies [43]. In addition, the slopes of well-being and AD were correlated in the unconditional model (i.e., the model without predictors), indicating that an increase in AD is linked with a decrease in well-being. However, this relationship should be treated with caution because it was not significant in the conditional model.
Finally, the significant variances of the intercepts in all models tested suggest considerable individual variabilities in the initial levels of AD and well-being. This heterogeneity might be attributed to a combination of sociodemographic characteristics, different impacts of the pandemic at the individual level (e.g., losing a loved one or having had COVID-19) or prior mental health difficulties [73,74,75]. The fluctuations in AD symptoms also differed significantly across individuals. This is in line with a recent study which identified four distinct trajectories of AD [54], highlighting the high variability in people’s reactions to the pandemic, as demonstrated earlier [76,77].

4.3. Predictors of Adjustment and Well-Being Trajectories

In our univariate conditional models, several sociodemographic variables could significantly predict the starting points of the trajectories of AD and well-being. Higher initial AD symptoms were identified in female participants, participants from Austria and those with lower income. On the other hand, the initial well-being scores were higher in older participants, participants from Croatia and those with higher income. These results are in line with previous COVID-19 studies which reported higher risks for AD for women and people with financial difficulties [17,18]. Notably, both gender and income could predict only the starting point but not how the symptoms of AD and well-being would change over time. For this reason, it can be assumed that women and people with low income were at risk even before the pandemic, and this risk remained throughout the pandemic. This idea is supported by a recent meta-analysis which showed that female gender and low income are reliable predictors of AD in different contexts [78]. Thus, these factors cannot only be considered pandemic-specific but rather universal, pointing to the importance of providing regular, low-threshold psychosocial support for women and other vulnerable groups, such as people with low income.
In the present study, older participants tended to have higher levels of well-being but also higher AD scores at baseline. This is contrary to a Spanish study which found that the elderly adapted better to COVID-19-related restrictions but also had more worries [79]. This discrepancy could be attributed to the timing of the two studies: while the Spanish study investigated the immediate psychological impact of the first Spanish lockdown, our study started several weeks after the end of the first lockdown in Austria and Croatia. It is therefore conceivable that the role of age varies depending on the stage of the pandemic. Considering the inconsistencies in the COVID-19 literature regarding the effects of age, future studies should more closely examine age differences in adapting to new environments and possible mechanisms behind, such as emotion regulation [22].
The estimation of the joint trajectories of AD and well-being revealed results comparable to those of the independent trajectories. Age and country were related to the initial levels of both AD and well-being, while gender and income were related only to the initial AD score. Sociodemographic characteristics thus seemed to influence the mental health of the participants at baseline but not the trajectories of their mental health over time.
It should be noted that education level was not a significant predictor in any of the models tested. Previous research delivered mixed results on the effect of education on AD [78]. Regarding well-being, a Danish study found a greater decline in well-being in highly educated people [31], whereas a lower level of education was associated with greater psychological distress in a German study [52]. To disentangle the relationship between education and mental health, samples in forthcoming studies should accurately represent the levels of education in the studied population.
As expected, the included time-dependent covariates, namely, the anxiety and depression scores, could significantly predict both AD and well-being scores at all timepoints in the univariate and multivariate models. While the anxiety and depression scores were positively correlated with the symptoms of AD, they were negatively correlated with well-being. These relationships are not surprising given that the symptoms of anxiety and depression were common emotional responses during past epidemics and natural disasters [80] and were repeatedly associated with AD symptoms during COVID-19 [18,23,24]. Moreover, depression and anxiety symptoms were among the most common indicators of well-being in longitudinal COVID-19 studies [30].
Notably, when predicting AD symptoms, higher estimates were observed for anxiety, whereas higher estimates for depression were observed when predicting well-being. According to some sources, anxiety is a normal reaction to the novel and unpredictable pandemic situation and is thus expected during the adjustment process [26,75]. Therefore, it can be assumed that AD with anxiety was the most frequent subtype of AD during COVID-19, which might explain the higher estimates for anxiety when predicting AD. On the other hand, well-being and depression are highly correlated and sometimes even perceived as opposite poles of the same continuum [45], which was possibly reflected in the stronger association between well-being and depression scores. Further research during pandemics and similar crises should assess AD subtypes and investigate the relationship between well-being and depression in a systematic way.

4.4. Cross-Country Differences in Adjustment and Well-Being Trajectories

Our study demonstrated several cross-country differences. On average, the Austrian participants had higher initial levels of AD and lower well-being. This finding is consistent with a large-scale international study which showed a greater mental health burden in Austria than in Croatia [37]. In the pandemic context, there are several possible explanations for the observed differences between the two countries. Firstly, COVID-19 severely impacted the global economy and job market [81] and, for many, achieving financial or professional goals was more difficult than prior to the pandemic. This is relevant as financial and work-related problems are commonly associated with AD [10,11]. The reduced possibility to achieve goals might have impacted the mental health of the Austrian participants more strongly given that satisfaction with life achievements was found to be the strongest predictor of well-being among Austrians [82]. Secondly, the incidence of COVID-19 was higher in Austria than in Croatia during the first assessment timepoint, which might have had an effect on the participants’ well-being and their adjustment to the pandemic at the baseline assessment. Thirdly, a large pan-European study reported that Austria was among the countries which experienced significant growth in their vulnerability levels, suggesting a greater mental health impact of COVID-19 on their population compared to some other European countries [83]. Finally, there is evidence that countries with histories of stress and trauma showed better mental health outcomes during the pandemic [19]. It can therefore be assumed that the pandemic was perceived less adversely than previous crises in Croatia (e.g., the Homeland War, 1991–1995, and major earthquakes in 2020 [84,85,86]), which might have facilitated psychological adjustment to the pandemic situation.
When comparing the multivariate models by country, no significant change in AD and well-being over time could be found. This was not anticipated given the significant rates of change observed in previous models. This discrepancy can be attributed to the smaller sizes of the national subsamples in comparison to the overall sample which was used in the previous models. On a descriptive level, the level of well-being tended to decline over time in both countries, whereas the symptoms of AD tended to increase over time in Austria and decrease in Croatia. It is thus conceivable that not only were the Austrian participants at a higher risk for AD at baseline but they were also more likely to develop AD symptoms as the pandemic progressed. However, a note of caution is due here since the effect was not significant. Interestingly, a slightly greater variability in terms of the fluctuations in AD symptoms was observed among the Austrian participants. One plausible reason for this might lie in the stringency index [34], which varied more strongly in Austria (ranging between 38.59 and 77.58) than in Croatia (ranging between 35.19 and 58.57). It has already been suggested that the intensity of the pandemic might account for mental health differences over time [20,21].
In Austria, older age and lower income were related to higher initial AD symptoms, whereas female gender was related to higher initial AD symptoms in Croatia. This finding is important as it implies that the predictors of mental health amidst pandemics might be country-dependent. Another peculiar finding was that anxiety and depression had different predictive powers depending on the timepoint of assessment. Compared to depression, anxiety scores were a stronger predictor of AD at T1 and T3 in Austria and at T2 and T4 in Croatia, with the respective timepoints corresponding to higher COVID-19 incidence rates. By way of illustration, the average incidence was higher in Austria than in Croatia at T1, and at this timepoint, the estimated value for anxiety was higher than the estimated value for depression in Austria and also higher than the estimated value in Croatia. This might be due to the presence of a greater fear of infection at times of higher infection rates, resulting in higher anxiety scores. This idea is supported by several studies which demonstrated the relationship between the fear of infection and AD symptoms [18,23]. In conclusion, the observed differences between Austria and Croatia highlight the importance of national studies which would consider country-related factors, e.g., cultural values and the pandemic situation. Such studies might provide a more complete picture of the psychological adjustment to the COVID-19 pandemic, as well as to other similar crises.

4.5. Strengths and Limitations

To the best of our knowledge, this is the first study to address changes in AD symptoms and well-being which was conducted in the later stages of the pandemic and in two countries. The lack of longitudinal studies is considered to be an important gap in AD research [13] which was addressed in this paper. The inclusion of both time-invariant and time-varying predictors and the examination of the co-development of AD and well-being trajectories via sophisticated methodology, was also innovative. Moreover, the present study had a longer assessment period than most COVID-19 studies, thus providing novel insights into the mental health impact of COVID-19. Another strength of our work is the investigation of AD and well-being in two countries that have rarely been jointly investigated in the past. This cross-country comparison is particularly relevant in COVID- 19 research since it adds to a better understanding of different responses to the pandemic in different contexts. Finally, a total of N = 531 participants responded to all four assessment waves, which is surprisingly good considering the length of the questionnaire and the long period of data collection.
However, the generalisability of our results is subject to certain limitations. First, convenience sampling was used in both countries, and the assessment waves were conducted at slightly different times (see Table 1). Second, the Austrian and Croatian samples were not of the same size, and they differed in sociodemographic characteristics, which might have biased the results. Third, the model fit was only satisfactory for some of the models tested, with the best model fit being found for the unconditional multivariate model. Future studies should include different covariates to assess their predictive power and check whether this could improve the model fit. Fourth, the Heywood case, one common problem of LGCMs, occurred when testing the final model in the Croatian sample. We managed to resolve the problem by freeing some of the previously fixed parameters in the Croatian model only, which might have limited its comparability with the Austrian model. Furthermore, pre-pandemic data and the data collected during the very first stage of the pandemic would have provided additional insights into the development of AD symptoms, but these data were not available for the countries investigated. Since the exploration of different subgroups using a Latent Class Analysis or Latent Profile Analysis was beyond the scope of our study, this might be an interesting topic for future research. Lastly, timestamps for individual survey responses were not available in our surveys, which is why it was not possible to consider the values of the stringency index and the incidence and death rates for all the dates during the data assessment period. This approach could be used in forthcoming investigations.

5. Conclusions

The present study sought to examine the trajectories of AD and well-being across Austria and Croatia and the factors underlying these trajectories. By focusing on mental health and well-being, our study makes a modest but meaningful contribution to achieving the third goal of the 2030 Agenda for Sustainable Development [87]. We found gender, age and income, as well as anxiety and depression scores, to be predictive of AD and the level of well-being during the COVID-19 pandemic, with the predictive powers of these covariates varying in the two countries. According to our results, mental health did not improve over time, highlighting the importance of psychosocial support during pandemics and similar global crises. In times of social restrictions, therapist-assisted computerised interventions might be particularly helpful [88,89]. Importantly, our results reveal a great variability in people’s responses to the pandemic, supporting the context-sensitive perspective in mental health research and clinical practice. Special attention should be paid to vulnerable groups, such as women and people with financial difficulties, to avoid further loss spirals and the development of full-blown mental disorders.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijerph20196861/s1, Table S1: Internal consistency of all measures across four assessment waves; Table S2: Flow of participants across four measurement timepoints; Table S3: Fit indices in models using strong measurement invariance as baseline; Figure S1: Austrian ANDM-8 scores over all four measurement times (n = 415); Figure S2: Croatian ANDM-8 scores over all four measurement times (n = 729); Figure S3: Austrian WHO-5 scores over all four measurement times (n = 415). Figure S4: Croatian WHO-5 scores over all four measurement times (n = 729); Figure S5: Adjustment disorder trajectories (Subsample 2); Figure S6: Adjustment disorder trajectories (Subsample 3); Figure S7: Well-being trajectories (Subsample 2); Figure S8: Well-being trajectories (Subsample 3).

Author Contributions

Conceptualization, I.Z.N., B.L.-S., M.A. and D.A.; methodology, I.Z.N. and A.S.; validation, J.K.L. and I.Z.N.; formal analysis, A.S., I.Z.N. and J.K.L.; investigation, I.Z.N., A.L., B.L.-S. and D.A.; writing—original draft preparation, I.Z.N.; writing—review and editing, B.L.-S., M.A., D.A., A.S. and I.Z.N.; visualization, A.S. and I.Z.N.; supervision, B.L.-S., M.A. and D.A.; project administration, B.L.-S. All authors have read and agreed to the published version of the manuscript.

Funding

We received no funding for conducting this research. The article processing charge (APC) was funded by the University of Vienna.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the University of Vienna in Austria (date of approval: 29 May 2020; reference number: 00554) and the Ethics Committee of the Department of Psychology, Faculty of Humanities and Social Sciences, at the University of Zagreb in Croatia (date of approval: 21 May 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study and the R codes for the main analysis are available upon request from the corresponding author. The data are not publicly available since an additional assessment wave is planned. As long as data collection and analysis are in process, excerpts of the data can only be shared upon reasonable request and if all authors agree.

Acknowledgments

The authors thank Hewan Giorgio Ammaturo, Theresa Wagner, and Stephan Ruben Marin (Team Austria) for their support in the literature review and table preparation. We further thank Ines Rezo Bagarić (Team Croatia) for her contribution in the initial phase of planning this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Javakhishvili, J.D.; Ardino, V.; Bragesjö, M.; Gorniak, J.; Kazlauskas, E.; Schäfer, I.; Schäflein, E.; Williamson, V. ESTSS Recommendations on Mental Health and Psychosocial Care During Pandemics. 2020. Available online: https://estss.org/wp-content/uploads/2020/05/ESTSS.Covid-19-Recommendations-2020-26-May.pdf (accessed on 26 September 2023).
  2. Giusti, L.; Salza, A.; Mammarella, S.; Bianco, D.; Ussorio, D.; Casacchia, M.; Roncone, R. #Everything will be fine. Duration of home confinement and “all-or-nothing” cognitive thinking style as predictors of traumatic distress in young university students on a digital platform during the COVID-19 Italian lockdown. Front. Psychiatry 2020, 11, 574812. [Google Scholar] [CrossRef] [PubMed]
  3. Kira, I.A.; Shuwiekh, H.A.M.; Ashby, J.S.; Elwakeel, S.A.; Alhuwailah, A.; Sous, M.S.F.; Baali, S.B.A.; Azdaou, C.; Oliemat, E.M.; Jamil, H.J. The impact of COVID-19 traumatic stressors on mental health: Is COVID-19 a new trauma type. Int. J. Ment. Health Addict. 2023, 21, 51–70. [Google Scholar] [CrossRef] [PubMed]
  4. Kazlauskas, E.; Quero, S. Adjustment and coronavirus: How to prepare for COVID-19 pandemic-related adjustment disorder worldwide? Psychol. Trauma Theory Res. Pract. Policy 2020, 12, S22–S24. [Google Scholar] [CrossRef] [PubMed]
  5. Brunet, A.; Rivest-Beauregard, M.; Lonergan, M.; Cipolletta, S.; Rasmussen, A.; Meng, X.; Jaafari, N.; Romero, S.; Superka, J.; Brown, A.D.; et al. PTSD is not the emblematic disorder of the COVID-19 pandemic; adjustment disorder is. BMC Psychiatry 2022, 22, 300. [Google Scholar] [CrossRef]
  6. Federal Ministry for European and International Affairs. Austria’s Neighbours. Available online: https://www.bmeia.gv.at/en/european-foreign-policy/foreign-policy/europe/austrias-neighbours/ (accessed on 26 September 2023).
  7. Helliwell, J.F.; Layard, R.; Sachs, J.D.; De Neve, J.-E.; Aknin, L.B.; Wang, S. World Happiness Report 2023; Sustainable Development Solutions Network: New York, NY, USA, 2023; ISBN 978-1-7348080-5-6. [Google Scholar]
  8. The World Bank. World Development Indicators. Available online: https://databank.worldbank.org/source/world-development-indicators (accessed on 31 August 2023).
  9. Hale, T.; Angrist, N.; Goldszmidt, R.; Kira, B.; Petherick, A.; Phillips, T.; Webster, S.; Cameron-Blake, E.; Hallas, L.; Majumdar, S.; et al. A global panel database of pandemic policies (Oxford COVID-19 Government Response Tracker). Nat. Hum. Behav. 2021, 5, 529–538. [Google Scholar] [CrossRef]
  10. Zelviene, P.; Kazlauskas, E. Adjustment disorder: Current perspectives. Neuropsychiatr. Dis. Treat. 2018, 14, 375–381. [Google Scholar] [CrossRef]
  11. O’Donnell, M.L.; Agathos, J.A.; Metcalf, O.; Gibson, K.; Lau, W. Adjustment disorder: Current developments and future directions. Int. J. Environ. Res. Public Health 2019, 16, 2537. [Google Scholar] [CrossRef]
  12. World Health Organization. The 11th Revision of the International Classification of Diseases. Available online: https://icd.who.int/browse11/l-m/en (accessed on 31 August 2023).
  13. O’Donnell, M.L.; Alkemade, N.; Creamer, M.; McFarlane, A.C.; Silove, D.; Bryant, R.A.; Felmingham, K.; Steel, Z.; Forbes, D. A longitudinal study of adjustment disorder after trauma exposure. Am. J. Psychiatry 2016, 173, 1231–1238. [Google Scholar] [CrossRef]
  14. Vancappel, A.; Jansen, E.; Bachem, R.; Bray, A.; Egreteau, L.; Réveillère, C.; Maercker, A.; El-Hage, W. Validation of the French ADNM-20 in the assessment of emotional difficulties resulting from COVID-19 quarantine and outbreak. BMC Psychol. 2021, 9, 3–9. [Google Scholar] [CrossRef]
  15. Jannini, T.B.; Rossi, R.; Socci, V.; Reda, F.; Pacitti, F.; Di Lorenzo, G. Psychometric and factorial validity of the International Adjustment Disorder Questionnaire (IADQ) in an Italian sample: A validation and prevalence estimate study. Clin. Psychol Psychother. 2022, 30, 436–445. [Google Scholar] [CrossRef]
  16. Bohlken, J.; Kostev, K.; Riedel-Heller, S.; Hoffmann, W.; Michalowsky, B. Effect of the COVID-19 pandemic on stress, anxiety, and depressive disorders in German primary care: A cross-sectional study. J. Psychiatr. Res. 2021, 143, 43–49. [Google Scholar] [CrossRef] [PubMed]
  17. Makhashvili, N.; Javakhishvili, J.D.; Sturua, L.; Pilauri, K.; Fuhr, D.C.; Roberts, B. The influence of concern about COVID-19 on mental health in the Republic of Georgia: A cross-sectional study. Glob. Health 2020, 16, 111. [Google Scholar] [CrossRef] [PubMed]
  18. Lotzin, A.; Krause, L.; Acquarini, E.; Ajdukovic, D.; Ardino, V.; Arnberg, F.; Böttche, M.; Bragesjö, M.; Dragan, M.; Figueiredo-Braga, M.; et al. Risk and protective factors, stressors, and symptoms of adjustment disorder during the COVID-19 pandemic—First results of the ESTSS COVID-19 pan-European ADJUST Study. Eur. J. Psychotraumatol. 2021, 12, 1964197. [Google Scholar] [CrossRef] [PubMed]
  19. Coiro, M.J.; Asraf, K.; Tzischinsky, O.; Hadar-Shoval, D.; Tannous-Haddad, L.; Wolfson, A.R. Sleep quality and COVID-19-related stress in relation to mental health symptoms among Israeli and U.S. adults. Sleep Health 2021, 7, 127–133. [Google Scholar] [CrossRef]
  20. Fioravanti, G.; Bocci Benucci, S.; Prostamo, A.; Banchi, V.; Casale, S. Effects of the COVID-19 pandemic on psychological health in a sample of Italian adults: A three-wave longitudinal study. Psychiatry Res. 2022, 315, 114705. [Google Scholar] [CrossRef] [PubMed]
  21. Witteveen, A.B.; Young, S.Y.; Cuijpers, P.; Ayuso-Mateos, J.L.; Barbui, C.; Bertolini, F.; Cabello, M.; Cadorin, C.; Downes, N.; Franzoi, D.; et al. COVID-19 and common mental health symptoms in the early phase of the pandemic: An umbrella review of the evidence. PLoS Med. 2023, 20, e1004206. [Google Scholar] [CrossRef]
  22. Dworakowski, O.; Huber, Z.M.; Meier, T.; Boyd, R.L.; Horn, A.B. Emotion regulation across the lifespan: Age differences in intrapersonal and interpersonal strategies for the adjustment to the COVID-19 pandemic in four countries. Aging Ment. Health 2022, 26, 2048–2053. [Google Scholar] [CrossRef]
  23. Mahat-Shamir, M.; Pitcho-Prelorentzos, S.; Kagan, M.; Kestler-Peleg, M.; Lavenda, O. Adjustment disorder in the face of COVID-19 outbreak: The impact of death anxiety, media exposure, fear of contagion and hypochondriasis symptoms. Omega 2021, 87, 1189–1206. [Google Scholar] [CrossRef]
  24. Hamama-Raz, Y.; Goodwin, R.; Leshem, E.; Ben-Ezra, M. The toll of a second lockdown: A longitudinal study. J. Affect. Disord. 2021, 294, 60–62. [Google Scholar] [CrossRef]
  25. Levin, Y.; Bachem, R.; Goodwin, R.; Hamama-Raz, Y.; Leshem, E.; Ben-Ezra, M. Relationship between adjustment disorder symptoms and probable diagnosis before and after second lockdown in Israel: Longitudinal symptom network analysis. BJPsych Open 2022, 8, e186. [Google Scholar] [CrossRef]
  26. Bartels, C.; Hessmann, P.; Schmidt, U.; Vogelgsang, J.; Ruhleder, M.; Kratzenberg, A.; Treptow, M.; Reh-Bergen, T.; Abdel-Hamid, M.; Heß, L.; et al. Medium-term and peri-lockdown course of psychosocial burden during the ongoing COVID-19 pandemic: A longitudinal study on patients with pre-existing mental disorders. Eur. Arch. Psychiatry Clin. Neurosci. 2022, 272, 757–771. [Google Scholar] [CrossRef] [PubMed]
  27. Ben-Ezra, M.; Hamama-Raz, Y.; Goodwin, R.; Leshem, E.; Levin, Y. Association between mental health trajectories and somatic symptoms following a second lockdown in Israel: A longitudinal study. BMJ Open 2021, 11, e050480. [Google Scholar] [CrossRef] [PubMed]
  28. Kazlauskas, E.; Gelezelyte, O.; Nomeikaite, A.; Zelviene, P. Posttraumatic stress disorder and adjustment disorder in Lithuanian healthcare in 2018–2020: A nation-wide cohort study of the effects of COVID-19 pandemic. Healthcare 2021, 9, 1422. [Google Scholar] [CrossRef] [PubMed]
  29. Salanti, G.; Peter, N.; Tonia, T.; Holloway, A.; White, I.R.; Darwish, L.; Low, N.; Egger, M.; Haas, A.D.; Fazel, S.; et al. The impact of the COVID-19 pandemic and associated control measures on the mental health of the general population: A systematic review and dose-response meta-analysis. Ann. Intern. Med. 2022, 175, 1560–1571. [Google Scholar] [CrossRef] [PubMed]
  30. Farrell, A.H.; Vitoroulis, I.; Eriksson, M.; Vaillancourt, T. Loneliness and well-being in children and adolescents during the COVID-19 pandemic: A systematic review. Children 2023, 10, 279. [Google Scholar] [CrossRef]
  31. Thygesen, L.C.; Møller, S.P.; Ersbøll, A.K.; Santini, Z.I.; Nielsen, M.B.D.; Grønbæk, M.K.; Ekholm, O. Decreasing mental well-being during the COVID-19 pandemic: A longitudinal study among Danes before and during the pandemic. J. Psychiatr. Res. 2021, 144, 151–157. [Google Scholar] [CrossRef] [PubMed]
  32. Lotzin, A.; Acquarini, E.; Ajdukovic, D.; Ardino, V.; Böttche, M.; Bondjers, K.; Bragesjö, M.; Dragan, M.; Grajewski, P.; Figueiredo-Braga, M.; et al. Stressors, coping and symptoms of adjustment disorder in the course of the COVID-19 pandemic—Study protocol of the European Society for Traumatic Stress Studies (ESTSS) pan-European study. Eur. J. Psychotraumatol. 2020, 11, 1780832. [Google Scholar] [CrossRef]
  33. Looi, M.-K. Covid-19: Is a second wave hitting Europe? BMJ 2020, 371, m4113. [Google Scholar] [CrossRef]
  34. Mathieu, E.; Ritchie, H.; Rodés-Guirao, L.; Appel, C.; Giattino, C.; Hasell, J.; Macdonald, B.; Dattani, S.; Beltekian, D.; Ortiz-Ospina, E.; et al. Coronavirus Pandemic (COVID-19). Available online: https://ourworldindata.org/coronavirus (accessed on 23 March 2023).
  35. Helliwell, J.F.; Layard, R.; Sachs, J.D.; De Neve, J.-E.; Aknin, L.B.; Wang, S. World Happiness Report 2022; Sustainable Development Solutions Network: New York, NY, USA, 2022; ISBN 978-1-7348080-2-5. [Google Scholar]
  36. Ivanović, V.; Šimac, B.; Trako Poljak, T. Subjective wellbeing in rural and urban central Europe: Evidence from the European Social Survey (2008 and 2018). Int. J. Sociol. 2022, 52, 397–419. [Google Scholar] [CrossRef]
  37. GBD 2019 Mental Disorders Collaborators. Global, regional, and national burden of 12 mental disorders in 204 countries and territories, 1990–2019: A systematic analysis for the Global Burden of Disease Study 2019. Lancet Psychiatry 2022, 9, 137–150. [Google Scholar] [CrossRef]
  38. World Health Organization. WHO Coronavirus (COVID-19) Dashboard. Available online: https://covid19.who.int/ (accessed on 30 August 2023).
  39. Bonanno, G.A.; Westphal, M.; Mancini, A.D. Resilience to loss and potential trauma. Annu. Rev. Clin. Psychol. 2011, 7, 511–535. [Google Scholar] [CrossRef] [PubMed]
  40. Tao, T.J.; Liang, L.; Liu, H.; Hobfoll, S.E.; Hou, W.K.; Bonanno, G.A. The interrelations between psychological outcome trajectories and resource changes amid large-scale disasters: A growth mixture modeling analysis. Transl. Psychiatry 2023, 13, 57. [Google Scholar] [CrossRef] [PubMed]
  41. Kazlauskas, E.; Gegieckaite, G.; Eimontas, J.; Zelviene, P.; Maercker, A. A brief measure of the International Classification of Diseases-11 adjustment disorder: Investigation of psychometric properties in an adult help-seeking sample. Psychopathology 2018, 51, 10–15. [Google Scholar] [CrossRef] [PubMed]
  42. Ben-Ezra, M.; Mahat-Shamir, M.; Lorenz, L.; Lavenda, O.; Maercker, A. Screening of adjustment disorder: Scale based on the ICD-11 and the Adjustment Disorder New Module. J. Psychiatr. Res. 2018, 103, 91–96. [Google Scholar] [CrossRef]
  43. Truskauskaite-Kuneviciene, I.; Kazlauskas, E.; Ostreikaite-Jurevice, R.; Brailovskaia, J.; Margraf, J. Positive mental health and adjustment following life-stressors among young adults. Curr. Psychol. 2022, 41, 1951–1956. [Google Scholar] [CrossRef]
  44. WHO. WHO (Five) Well-Being Index; WHO Regional office for Europe: Hillerod, Denmark, 1998. [Google Scholar]
  45. Krieger, T.; Zimmermann, J.; Huffziger, S.; Ubl, B.; Diener, C.; Kuehner, C.; Grosse Holtforth, M. Measuring depression with a well-being index: Further evidence for the validity of the WHO Well-Being Index (WHO-5) as a measure of the severity of depression. J. Affect. Disord. 2014, 156, 240–244. [Google Scholar] [CrossRef]
  46. Topp, C.W.; Østergaard, S.D.; Søndergaard, S.; Bech, P. The WHO-5 Well-Being Index: A systematic review of the literature. Psychother. Psychosom. 2015, 84, 167–176. [Google Scholar] [CrossRef]
  47. Toffolutti, V.; Plach, S.; Maksimovic, T.; Piccitto, G.; Mascherini, M.; Mencarini, L.; Aassve, A. The association between COVID-19 policy responses and mental well-being: Evidence from 28 European countries. Soc. Sci. Med. 2022, 301, 114906. [Google Scholar] [CrossRef]
  48. Kroenke, K.; Spitzer, R.L.; Williams, J.B.W.; Löwe, B. An ultra-brief screening scale for anxiety and depression: The PHQ-4. Psychosomatics 2009, 50, 613–621. [Google Scholar] [CrossRef]
  49. Löwe, B.; Wahl, I.; Rose, M.; Spitzer, C.; Glaesmer, H.; Wingenfeld, K.; Schneider, A.; Brähler, E. A 4-item measure of depression and anxiety: Validation and standardization of the Patient Health Questionnaire-4 (PHQ-4) in the general population. J. Affect. Disord. 2010, 122, 86–95. [Google Scholar] [CrossRef]
  50. Marvaldi, M.; Mallet, J.; Dubertret, C.; Moro, M.R.; Guessoum, S.B. Anxiety, depression, trauma-related, and sleep disorders among healthcare workers during the COVID-19 pandemic: A systematic review and meta-analysis. Neurosci. Biobehav. Rev. 2021, 126, 252–264. [Google Scholar] [CrossRef] [PubMed]
  51. Eisinga, R.; te Grotenhuis, M.; Pelzer, B. The reliability of a two-item scale: Pearson, Cronbach, or Spearman-Brown? Int. J. Public Health 2013, 58, 637–642. [Google Scholar] [CrossRef] [PubMed]
  52. Liu, S.; Haucke, M.N.; Heinzel, S.; Heinz, A. Long-Term Impact of Economic Downturn and Loneliness on Psychological Distress: Triple Crises of COVID-19 Pandemic. J. Clin. Med. 2021, 10, 4596. [Google Scholar] [CrossRef] [PubMed]
  53. Zhou, M.; Kan, M.-Y. The varying impacts of COVID-19 and its related measures in the UK: A year in review. PLoS ONE 2021, 16, e0257286. [Google Scholar] [CrossRef]
  54. Kazlauskas, E.; Elklit, A.; Truskauskaite, I. A Longitudinal Course of ICD-11 Adjustment Disorder Symptom Profiles: A 12-Month Follow-Up Study. Psychopathology 2022, 55, 373–381. [Google Scholar] [CrossRef]
  55. IBM Corp. IBM SPSS Statistics for Windows, Version 27.0; IBM Corp.: Armonk, NY, USA, 2020. [Google Scholar]
  56. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022. [Google Scholar]
  57. West, S.G.; Taylor, A.B.; Wu, W. Model fit and model selection in structural equation modeling. In Handbook of Structural Equation Modeling; Hoyle, R.H., Ed.; Guilford Press: New York, NY, USA, 2012; pp. 209–231. ISBN 9781606230770. [Google Scholar]
  58. Little, T.D.; Card, N.A. Longitudinal Structural Equation Modeling; The Guilford Press: New York, NY, USA; London, UK, 2013; ISBN 9781462510160. [Google Scholar]
  59. Putnick, D.L.; Bornstein, M.H. Measurement invariance conventions and reporting: The state of the art and future directions for psychological research. Dev. Rev. 2016, 41, 71–90. [Google Scholar] [CrossRef]
  60. Kline, R.B. Principles and Practice of Structural Equation Modeling, 4th ed.; The Guilford Press: New York, NY, USA; London, UK, 2016; ISBN 9781462523344. [Google Scholar]
  61. Li, C.-H. Confirmatory factor analysis with ordinal data: Comparing robust maximum likelihood and diagonally weighted least squares. Behav. Res. Methods 2016, 48, 936–949. [Google Scholar] [CrossRef]
  62. Rosseel, Y. lavaan: An R package for structural equation modeling. J. Stat. Softw. 2012, 48, 1–36. [Google Scholar] [CrossRef]
  63. Curran, P.J.; Obeidat, K.; Losardo, D. Twelve frequently asked questions about growth curve modeling. J. Cogn. Dev. 2010, 11, 121–136. [Google Scholar] [CrossRef]
  64. Newsom, J.T. Longitudinal Structural Equation Modeling; Routledge: New York, NY, USA; East Sussex, UK, 2015; ISBN 9781317975359. [Google Scholar]
  65. Bollen, K.A.; Curran, P.J. Latent Curve Models: A Structural Equation Perspective; Wiley-Interscience: Hoboken, NJ, USA, 2006; ISBN 978-0-471-45592-9. [Google Scholar]
  66. Preacher, K.; Wichman, A.; MacCallum, R.; Briggs, N. Latent Growth Curve Modeling; SAGE Publications, Inc.: Thousand Oaks, CA, USA, 2008; ISBN 9781412939553. [Google Scholar]
  67. Gana, K. Structural Equation with Lavaan; John Wiley & Sons Incorporated: Newark, NJ, USA, 2018; ISBN 978-1-786-30369-1. [Google Scholar]
  68. Rosseel, Y. The Lavaan Tutorial; Department of Data Analysis, Ghent University: Ghent, Belgium, 2023; Available online: https://lavaan.ugent.be/tutorial/tutorial.pdf (accessed on 31 July 2023).
  69. Little, T.D.; Jorgensen, T.D.; Lang, K.M.; Moore, E.W.G. On the joys of missing data. J. Pediatr. Psychol. 2014, 39, 151–162. [Google Scholar] [CrossRef]
  70. Cohen, J. Statistical Power Analysis for the Behavioral Sciences, 2nd ed.; Lawrence Erlbaum Associates: Hillsdale, NJ, USA, 1988; ISBN 0805802835. [Google Scholar]
  71. Berlin, K.S.; Parra, G.R.; Williams, N.A. An introduction to latent variable mixture modeling (part 2): Longitudinal latent class growth analysis and growth mixture models. J. Pediatr. Psychol. 2014, 39, 188–203. [Google Scholar] [CrossRef] [PubMed]
  72. Levin, Y.; Bachem, R.; Hyland, P.; Karatzias, T.; Shevlin, M.; Ben-Ezra, M.; Maercker, A. Validation of the International Adjustment Disorder Questionnaire in Israel and Switzerland. Clin. Psychol. Psychother. 2022, 29, 1321–1330. [Google Scholar] [CrossRef] [PubMed]
  73. Fineberg, N.A.; Pellegrini, L.; Wellsted, D.; Hall, N.; Corazza, O.; Giorgetti, V.; Cicconcelli, D.; Theofanous, E.; Sireau, N.; Adam, D.; et al. Facing the “new normal”: How adjusting to the easing of COVID-19 lockdown restrictions exposes mental health inequalities. J. Psychiatr. Res. 2021, 141, 276–286. [Google Scholar] [CrossRef]
  74. Gersons, B.P.R.; Smid, G.E.; Smit, A.S.; Kazlauskas, E.; McFarlane, A.C. Can a ‘second disaster’ during and after the COVID-19 pandemic be mitigated? Eur. J. Psychotraumatol. 2020, 11, 1815283. [Google Scholar] [CrossRef] [PubMed]
  75. Chen, S.; Bi, K.; Sun, P.; Bonanno, G.A. Psychopathology and resilience following strict COVID-19 lockdowns in Hubei, China: Examining person- and context-level predictors for longitudinal trajectories. Am. Psychol. 2022, 77, 262–275. [Google Scholar] [CrossRef]
  76. Sandbakken, E.M.; Moss, S.M. “Now we are all in the same boat. At the same time, we are not.” Meaning-making and coping under COVID-19 lockdown in Norway. Hum. Arenas 2021, 6, 201–225. [Google Scholar] [CrossRef]
  77. Zrnić Novaković, I.; Ajduković, D.; Bakić, H.; Borges, C.; Figueiredo-Braga, M.; Lotzin, A.; Anastassiou-Hadjicharalambous, X.; Lioupi, C.; Javakhishvili, J.D.; Tsiskarishvili, L.; et al. Shaped by the COVID-19 pandemic: Psychological responses from a subjective perspective-A longitudinal mixed-methods study across five European countries. PLoS ONE 2023, 18, e0285078. [Google Scholar] [CrossRef]
  78. Spanovic Kelber, M.; Morgan, M.A.; Beech, E.H.; Smolenski, D.J.; Bellanti, D.; Galloway, L.; Ojha, S.; Otto, J.L.; Wilson, A.L.G.; Bush, N.; et al. Systematic review and meta-analysis of predictors of adjustment disorders in adults. J. Affect. Disord. 2022, 304, 43–58. [Google Scholar] [CrossRef]
  79. Morales-Vives, F.; Dueñas, J.-M.; Vigil-Colet, A.; Camarero-Figuerola, M. Psychological variables related to adaptation to the COVID-19 lockdown in Spain. Front. Psychol. 2020, 11, 565634. [Google Scholar] [CrossRef]
  80. Esterwood, E.; Saeed, S.A. Past epidemics, natural disasters, COVID19, and mental health: Learning from history as we deal with the present and prepare for the future. Psychiatr. Q. 2020, 91, 1121–1133. [Google Scholar] [CrossRef]
  81. Nicola, M.; Alsafi, Z.; Sohrabi, C.; Kerwan, A.; Al-Jabir, A.; Iosifidis, C.; Agha, M.; Agha, R. The socio-economic implications of the coronavirus pandemic (COVID-19): A review. Int. J. Surg. 2020, 78, 185–193. [Google Scholar] [CrossRef] [PubMed]
  82. Jovanović, V. Evaluation of domain satisfaction measure in cross-national perspective: Evidence from Austria and four countries of the former Yugoslavia. Soc. Indic. Res. 2018, 141, 1369–1385. [Google Scholar] [CrossRef]
  83. Grané, A.; Albarrán, I.; Merchán, D.E. Impact of Pandemic on European Well-Being: Visualizing Scenarios from the SHARE Database. Int. J. Environ. Res. Public Health 2021, 18, 4620. [Google Scholar] [CrossRef] [PubMed]
  84. Gizdic, A.; Baxter, T.; Barrantes-Vidal, N.; Park, S. Loneliness and psychosocial predictors of psychosis-proneness during COVID-19: Preliminary findings from Croatia. Psychiatry Res. 2022, 317, 114900. [Google Scholar] [CrossRef] [PubMed]
  85. Romic, I.; Silovski, H.; Mance, M.; Pavlek, G.; Petrovic, I.; Figl, J.; Grbavac, D.; Moric, T.; Romic, R.; Bakula, B.; et al. Psychological effects of “double crisis” (COVID-19 pandemic and earthquakes) on Croatian medical students. Psychiatr. Danub. 2021, 33, 120–125. [Google Scholar]
  86. Matić, I.; Takšić, I.; Božičević, M. The sense of coherence and subjective well-being as resources of resilience in the time of stressful situations: COVID-19 outbreak and earthquakes. Psychiatr. Danub. 2021, 33, 639–645. [Google Scholar] [CrossRef]
  87. UN General Assembly. Transforming Our World: The 2030 Agenda for Sustainable Development; A/RES/70/1; 2015. Available online: https://undocs.org/en/A/RES/70/1 (accessed on 18 September 2023).
  88. Ehlers, A.; Wild, J.; Warnock-Parkes, E.; Grey, N.; Murray, H.; Kerr, A.; Rozental, A.; Thew, G.; Janecka, M.; Beierl, E.T.; et al. Therapist-assisted online psychological therapies differing in trauma focus for post-traumatic stress disorder (STOP-PTSD): A UK-based, single-blind, randomised controlled trial. Lancet Psychiatry 2023, 10, 608–622. [Google Scholar] [CrossRef]
  89. Salza, A.; Giusti, L.; Ussorio, D.; Casacchia, M.; Roncone, R. Cognitive behavioral therapy (CBT) anxiety management and reasoning bias modification in young adults with anxiety disorders: A real-world study of a therapist-assisted computerized (TACCBT) program vs. “person-to-person” group CBT. Internet Interv. 2020, 19, 100305. [Google Scholar] [CrossRef]
Figure 1. Conceptional model of AD symptoms and well-being with covariances, time-invariant covariates (TICs; yellow lines) and time-varying covariates (TVCs; blue lines). The solid grey lines represent the factor loadings, and the dashed grey lines represent the covariances. Note: ADNM-8 = 8-item Adjustment Disorder New Module; AD = adjustment disorder; WB = well-being; AS = anxiety scores (PHQ-2 = 2-item Patient Health Questionnaire); DS = depression score (GAD-2 = 2-item Generalized Anxiety Disorder Scale).
Figure 1. Conceptional model of AD symptoms and well-being with covariances, time-invariant covariates (TICs; yellow lines) and time-varying covariates (TVCs; blue lines). The solid grey lines represent the factor loadings, and the dashed grey lines represent the covariances. Note: ADNM-8 = 8-item Adjustment Disorder New Module; AD = adjustment disorder; WB = well-being; AS = anxiety scores (PHQ-2 = 2-item Patient Health Questionnaire); DS = depression score (GAD-2 = 2-item Generalized Anxiety Disorder Scale).
Ijerph 20 06861 g001
Figure 2. The estimated growth trajectories based on the unconditional univariate models. (a) Adjustment disorder trajectory. (b) Well-being trajectory.
Figure 2. The estimated growth trajectories based on the unconditional univariate models. (a) Adjustment disorder trajectory. (b) Well-being trajectory.
Ijerph 20 06861 g002
Table 1. The COVID-19 situation in Austria and Croatia at four data collection timepoints.
Table 1. The COVID-19 situation in Austria and Croatia at four data collection timepoints.
T1T2T3T4
AustriaCroatiaAustriaCroatiaAustriaCroatiaAustriaCroatia
Recruitment period27 June 2020–
22 September 2020
15 June 2020–
16 August 2020
14 January 2021–
29 March 2021
30 November 2020–
7 January 2021
13 July 2021–
8 October 2021
21 June 2021–
26 July 2021
26 November 2021–
13 December 2021 a
8 December 2021–
11 January 2022
Duration of data collection88 days ≈
13 weeks
63 days ≈
9 weeks
75 days ≈
11 weeks
39 days ≈
6 weeks
88 days ≈ 13 weeks36 days ≈
5 weeks
18 days ≈
3 weeks
35 days
≈ 5 weeks
Stringency index
M (SD)38.5944.9877.5858.5750.1135.1951.7036.97
Range36.11–50.0035.19–54.6373.15–82.4147.22–67.5946.23–55.0530.49–41.6749.30–52.0633.26–38.08
Incidence
M (SD)24.73 (20.42)14.55 (7.20)214.75 (67.75)639.83 (243.16)126.58 (72.32)22.81 (5.40)413.46 (418.22)613.41 (350.14)
Range4.25–83.340.18–31.09142.82–351.49276.26–905.0217.45–232.9816.91–36.370.00–1257.420.00–981.76
Deaths
M (SD)0.19 (0.06)0.22 (0.19)4.57 (1.38)16.35 (1.92)0.59 (0.43)0.59 (0.29)7.50 (0.34)12.15 (1.88)
Range0.06–0.370.00–0.603.02–7.3412.73–19.570.16–1.370.25–1.496.76–7.869.39–15.24
Note: The stringency index is a composite measure of government policies, with higher scores indicating a stricter response (i.e., 100 = strictest response [9]). Data on incidence are defined as new confirmed cases of COVID-19 per one million people. Data on death rates are defined as new deaths attributed to COVID-19 per one million people. The values reported in the table are the averages of the daily values recorded during the respective recruitment periods. All data are based on the dataset provided by Our World in Data (last accessed on 23 March 2023 [34]). a This data collection period was shorter than the others in order to reflect the duration of the lockdown implemented in Austria in the winter of 2021.
Table 2. Participants’ characteristics at T1.
Table 2. Participants’ characteristics at T1.
Total Sample
N = 1144
Austria
n = 415
Croatia
n = 729
AgeM = 44.0
SD = 13.5
M = 46.5
SD = 14.7
M = 42.7
SD = 12.6
n (%)n (%)n (%)
Gender a
Male300 (26.2)128 (30.8)172 (23.6)
Female841 (73.6)285 (68.7)556 (76.4)
Other2 (0.2)2 (0.5)0 (0.0)
Education
Low6 (0.5)5 (1.2)1 (0.1)
Middle306 (26.7)166 (40.0)140 (19.2)
High832 (72.7)244 (58.8)588 (80.7)
Income b
Very low126 (11.3)43 (11.2)83 (11.4)
Low346 (31.1)140 (36.5)206 (28.3)
Medium344 (30.9)31 (8.1)313 (43.0)
High296 (26.6)170 (44.3)126 (17.3)
Relationship status
Single269 (23.5)102 (24.6)167 (22.9)
In a relationship875 (76.5)313 (75.4)562 (77.1)
Employment status c
Training/Study124 (25.5)52 (12.5)72 (9.9)
Employed part-time154 (34.5)122 (29.4)32 (4.4)
Employed full-time742 (77.0)193 (46.5)549 (75.3)
Self-employed68 (15.2)36 (8.7)32 (4.4)
Retired84 (19.0)57 (13.7)27 (3.7)
Job-seeking57 (12.3)8 (1.9)49 (6.7)
Other26 (6.0)10 (2.4)16 (2.2)
Diagnosis of a mental disorder
Yes167 (14.6)96 (23.1)71 (9.7)
No977 (85.4)319 (76.9)658 (90.3)
Note: a N = 1143 (nAustria = 415; nCroatia = 728). b N = 1112 (nAustria = 384; nCroatia = 728). c Percentages sum to more than 100 as multiple responses were possible.
Table 3. Average scores of the assessed mental health outcomes and self-reported prevalence rates (i.e., the percentage of participants above the cut-off) of mental health disorders at four timepoints.
Table 3. Average scores of the assessed mental health outcomes and self-reported prevalence rates (i.e., the percentage of participants above the cut-off) of mental health disorders at four timepoints.
TimepointADNM-8PHQ-2GAD-2WHO-5 a
M (SD)% (n)M (SD)% (n)M (SD)% (n)M (SD)% (n)
T1
n = 1144
14.3 (5.3)9.8 (112)1.3 (1.5)15.4 (176)1.4 (1.3)15.8 (181)58.0 (21.4)33.7 (386)
T2
n = 993
15.8 (5.7)15.1 (150)1.5 (1.5)16.2 (161)1.3 (1.5)15.5 (154)51.5 (22.1)45.1 (448)
T3
n = 837
14.3 (5.4)9.9 (83)1.2 (1.4)12.9 (108)1.1 (1.4)11.6 (97)57.3 (21.6)34.5 (289)
T4
n = 674
15.0 (5.6)11.7 (79)1.3 (1.6)16.2 (109)1.3 (1.5)14.5 (98)53.8 (22.7)40.8 (275)
Note: ADNM-8 = 8-item Adjustment Disorder New Module; WHO-5 = 5-item World Health Organization Well-Being Index; PHQ-2 = 2-item Patient Health Questionnaire; GAD-2 = 2-item Generalized Anxiety Disorder Scale. The prevalence rates were calculated based on the following established cut-off scores: ADNM-8 > 22, PHQ-2 > 2. GAD-2 > 2 and WHO-5 ≤ 50. a The percentages suggest suboptimal well-being and likely depression [46].
Table 4. Comparison of fit indices and parameter estimates in univariate models.
Table 4. Comparison of fit indices and parameter estimates in univariate models.
Univariate
Unconditional Model
Univariate Conditional Model with TICsUnivariate Conditional Model with TVCsUnivariate Model with TICs and TVCs
ADWBADWBADWBADWB
Fit indices
Χ2 (df)105.89 *** (5)136.80 *** (5)131.84 *** (15)158.36 *** (15)148.70 *** (29)140.68 *** (29)169.82 *** (39)147.65 *** (39)
CFI0.9330.9250.9290.9240.9250.9400.9190.941
TLI0.9190.9100.8760.8690.9010.9210.8790.912
RMSEA
[90% CI]
0.165
[0.136, 0.196]
0.182 [0.153, 0.213]0.101
[0.084, 0.119]
0.107
[0.091, 0.125]
0.097
[0.082, 0.113]
0.091
[0.076, 0.107]
0.089
[0.075, 0.103]
0.078
[0.065, 0.092]
SRMR0.0630.0680.0340.0370.1190.1060.0850.070
Parameter estimates (SE)
Intercept14.57 *** (0.16)56.77 *** (0.64)16.34 *** (1.28)39.88 *** (4.70)11.73 *** (0.22)69.32 *** (0.85)10.78 *** (1.34)61.49 *** (5.17)
Slope0.15 * (0.06)−1.10 *** (0.22)0.17 (0.47)−2.76 (1.67)0.24 ** (0.09)−1.34 *** (0.30)−0.33 (0.56)−0.76 (1.84)
Covariances (SE)
i ~~ s−0.78 (0.51)8.31 (6.71)−0.97 (0.50)5.98 (6.86)0.17 (0.37)4.04 (5.05)0.07 (0.37)2.85 (5.10)
Variances (SE)
Intercept18.87 *** (1.38)286.64 *** (19.74)18.31 *** (1.36)274.99 *** (20.02)7.13 ***
(1.11)
114.14 *** (14.58)6.52 ***
(1.16)
115.11 ***
(14.53)
Slope1.24 *** (0.30)3.66 (3.80)1.30 *** (0.29)3.78 (3.82)0.59* (0.23)0.83 (2.88)0.58 * (0.23)1.47 (2.94)
TICs (SE)
Gender--i:
s:
1.25 ** (0.37)
0.04 (0.13)
−1.55 (1.48)
−0.52 (0.49)
--1.10 ** (0.39)
0.04 (0.15)
−0.88 (1.60)
−0.09 (0.51)
Age--i:
s:
−0.00 (0.01)
0.01 (0.00)
0.22 *** (0.05)
−0.00 (0.02)
--0.04 ** (0.01)
0.01 (0.01)
0.11 * (0.05)
−0.01 (0.02)
Country--i:
s:
−1.12 ** (0.35)
−0.09 (0.14)
4.04 ** (1.36)
0.93 (0.47)
--−0.81 * (0.38)
−0.00 (0.15)
2.75 (1.40)
−0.42 (0.50)
Education--i:
s:
−0.64 (0.39)
−0.08 (0.16)
−0.01 (1.52)
0.57 (0.56)
--−0.42 (0.44)
0.05 (0.18)
0.51 (1.76)
0.34 (0.58)
Income--i:
s:
−0.51 ** (0.17)
−0.01 (0.07)
1.86 ** (0.67)
0.08 (0.23)
--−0.36 (0.19)
0.03 (0.08)
−0.57 (0.70)
0.14 (0.25)
TVCs (SE)
T1Depression
Anxiety
----1.00 *** (0.17)
1.01 *** (0.17)
−5.28 *** (0.52)
−3.25 *** (0.48)
0.98 *** (0.17)
1.00 *** (0.18)
−4.99 *** (0.53)
−3.42 *** (0.48)
T2Depression
Anxiety
----1.29 *** (0.17)
1.31 *** (0.18)
−6.31 *** (0.54)
−4.54 *** (0.53)
1.26 *** (0.18)
1.32 *** (0.18)
−6.30 *** (0.55)
−4.46 *** (0.55)
T3Depression
Anxiety
----0.85 *** (0.17)
1.27 *** (0.17)
−5.63 *** (0.62)
−2.97 *** (0.67)
0.94 *** (0.17)
1.16 *** (0.17)
−5.39 *** (0.64)
−3.23 *** (0.68)
T4Depression
Anxiety
----0.94 *** (0.14)
1.01 *** (0.15)
−6.03 *** (0.58)
−3.23 *** (0.59)
0.96 *** (0.15)
1.04 *** (0.16)
−5.92 *** (0.62)
−3.36 *** (0.62)
Note: TICs = time-invariant covariates; TVCs = time-varying covariates; AD = adjustment disorder; WB = well-being; CFI = Comparative Fit Index; TLI = Tucker–Lewis Index; RMSEA = Root Mean Square Error of Approximation; CI = Confidence Interval; SRMR = Standardised Root Mean Square Residual; SE = Standard Error; i = intercept. s = slope. The robust CFI, TLI and RMSEA, along with the scaled SRMR, are reported. The literature suggests considering a lower CI limit CI when interpreting the RMSEA [60,65]; conventionally, values of 0.05–0.08 indicate an acceptable fit, and values ≤ 0.05 indicate a very good fit [58]. p = 0.05. * p < 0.05. ** p < 0.01. *** p < 0.001.
Table 5. Comparison of fit indices and parameter estimates in the multivariate models.
Table 5. Comparison of fit indices and parameter estimates in the multivariate models.
Multivariate Unconditional Model aMultivariate Conditional Model
(All Predictors)
Fit indices
Χ2 (df)188.49 *** (21)310.69 *** (89)
CFI0.9780.937
TLI0.9710.906
RMSEA
[90% CI]
0.081
[0.071, 0.092]
0.075
[0.066, 0.084]
SRMR0.0600.103
Parameter estimates   (SE)
Interceptadnm_i: 14.58 *** (0.22)well_i: 57.33 *** (0.86)adnm_i: 10.86 *** (1.36)well_i: 61.39 *** (5.21)
Slopeadnm_s: 0.10 (0.07)well_s: −1.03 *** (0.25)adnm_s: −0.35 (0.56)well_s: −0.92 (1.83)
Covariances   (SE)
adnm_i ~~ adnm_s−0.85 (0.55)0.09 (0.37)
well_i ~~ well_s9.08 (6.69)3.02 (5.04)
adnm_i ~~ well_i−52.34 *** (5.45)−7.41 * (3.10)
adnm_s ~~ well_s−1.19 * (0.55)0.04 (0.38)
Variances   (SE)
Interceptadnm: 20.18 *** (1.94)well: 277.56 *** (24.23)adnm: 6.65 *** (1.18)well: 116.78 *** (14.73)
Slopeadnm: 1.33 *** (0.35)well: 4.66 (3.99)adnm: 0.52 * (0.23)well: 0.98 (2.88)
Time-invariant covariates   (SE)
Gender--adnm_i: 1.13 ** (0.40)
adnm_s: 0.03 (0.15)
well_i: −0.94 (1.60)
well_s: −0.07 (0.52)
Age--adnm_i: 0.04 ** (0.01)
adnm_s: 0.01 (0.01)
well_i: 0.11 * (0.05)
well_s: −0.01 (0.02)
Country--adnm_i: −0.82 * (0.38)
adnm_s: 0.01 (0.15)
well_i: 2.74 (1.40)
well_s: −0.39 (0.50)
Education--adnm_i: −0.43 (0.44)
adnm_s: 0.05 (0.18)
well_i: 0.54 (1.76)
well_s: 0.32 (0.58)
Income--adnm_i: −0.37 (0.19)
adnm_s: 0.03 (0.08)
well_i: −0.56 (0.71)
well_s: 0.16 (0.25)
Time-varying covariates   (SE) Adjustment disorderWell-being
T1Depression
Anxiety
0.96 *** (0.18)
0.98 *** (0.18)
−4.95 *** (0.53)
−3.37 *** (0.49)
T2Depression
Anxiety
1.24 *** (0.18)
1.30 *** (0.18)
−6.26 *** (0.56)
−4.38 *** (0.56)
T3Depression
Anxiety
0.93 *** (0.17)
1.15 *** (0.17)
−5.31 *** (0.63)
−3.15 *** (0.68)
T4Depression
Anxiety
0.94 *** (0.15)
1.04 *** (0.16)
−5.86 *** (0.61)
−3.22 *** (0.61)
Note: CFI = Comparative Fit Index; TLI = Tucker–Lewis Index; RMSEA = Root Mean Square Error of Approximation; CI = Confidence Interval; SRMR = Standardised Root Mean Square Residual; SE = Standard Error; i = intercept; s = slope. The robust CFI, TLI and RMSEA, along with the scaled SRMR, are reported. The literature suggests considering a lower CI limit CI when interpreting the RMSEA [60,65]; conventionally, values of 0.05–0.08 indicate an acceptable fit, and values ≤ 0.05 indicate a very good fit [58]. a Since some estimated variances were negative in the model using the MLR as an estimator, the WLSMV was used as an alternative estimator. p = 0.05. * p < 0.05. ** p < 0.01. *** p < 0.001.
Table 6. Comparison of multivariate models by country.
Table 6. Comparison of multivariate models by country.
AustriaCroatia a
Fit indices
  Χ2 (df)185.19 *** (85)229.42 *** (84)
  CFI0.9320.938
  TLI0.9010.908
  RMSEA
[90% CI]
0.081
[0.064, 0.097]
0.076
[0.064, 0.088]
  SRMR0.1090.112
Parameter
estimates   (SE)
  Interceptadnm_i: 9.93 *** (2.00)well_i: 61.37 *** (7.59)adnm_i: 9.44 *** (1.50)well_i: 76.33 *** (6.85)
  Slopeadnm_s: 0.12 (0.91)well_s: −3.28 (2.49)adnm_s: −0.84 (0.57)well_s: −7.05 (4.81)
Covariances   (SE)
  adnm_i ~~ adnm_s−0.13 (0.71)0.10 (0.39)
  well_i ~~ well_s−2.82 (10.05)−3.90 (29.67)
  adnm_i ~~ well_i−13.05 * (6.59)−2.72 (3.10)
  adnm_s ~~ well_s0.40 (0.77)−0.72 (1.53)
Variances   (SE)
  Interceptadnm: 6.93 ** (2.17)well: 156.40 *** (28.84)adnm: 6.46 *** (1.27)well: 100.84 ** (34.12)
  Slopeadnm: 0.85 (0.45)well: 3.37 (5.56)adnm: 0.36 (0.21)well: 16.97 (34.22)
Time-invariant
covariates   (SE)
  Genderadnm_i: 0.98 (0.66)
adnm_s: −0.14 (0.27)
well_i: 2.27 (2.49)
well_s: 1.25 (0.83)
adnm_i: 1.09 * (0.48)
adnm_s: 0.23 (0.16)
well_i: −4.37 (2.38)
well_s: −1.08 (1.87)
  Ageadnm_i: 0.06 ** (0.02)
adnm_s: 0.00 (0.01)
well_i: 0.08 (0.08)
well_s: −0.02 (0.02)
adnm_i: 0.03 (0.02)
adnm_s: 0.01 (0.01)
well_i: 0.07 (0.07)
well_s: 0.04 (0.05)
  Educationadnm_i: −0.62 (0.63)
adnm_s: −0.01 (0.30)
well_i: 0.82 (2.43)
well_s: 0.13 (0.80)
adnm_i: −0.35 (0.62)
adnm_s: 0.14 (0.21)
well_i: 0.90 (2.62)
well_s: −0.12 (1.92)
  Incomeadnm_i: −0.64 * (0.28)
adnm_s: 0.08 (0.12)
well_i: −0.14 (1.02)
well_s: 0.07 (0.33)
adnm_i: −0.08 (0.25)
adnm_s: −0.05 (0.09)
well_i: −1.31 (1.06)
well_s: 1.25 (0.91)
Time-varying
covariates   (SE)
Adjustment disorderWell-beingAdjustment disorderWell-being
T1Depression
Anxiety
0.90 ** (0.31)
1.54 *** (0.31)
−5.36 *** (0.83)
−5.09 *** (0.85)
0.93 *** (0.21)
0.75 *** (0.19)
−5.00 *** (0.68)
−3.42 *** (0.61)
T2Depression
Anxiety
1.61 *** (0.26)
1.18 *** (0.26)
−7.72 *** (0.85)
−4.33 *** (0.84)
0.92 *** (0.23)
1.45 *** (0.25)
−4.24 *** (0.70)
−3.55 *** (0.74)
T3Depression
Anxiety
0.75 * (0.29)
1.70 *** (0.25)
−5.21 *** (0.82)
−4.57 *** (1.01)
1.10 *** (0.20)
0.69 ** (0.21)
−5.55 *** (0.88)
−2.07 ** (0.94)
T4Depression
Anxiety
1.37 *** (0.26)
1.11 *** (0.27)
−6.59 *** (1.02)
−3.27 *** (0.94)
0.53 ** (0.19)
1.03 *** (0.20)
−5.24 *** (0.73)
−3.59 *** (0.78)
Note: CFI = Comparative Fit Index; TLI = Tucker–Lewis Index; RMSEA = Root Mean Square Error of Approximation; CI = Confidence Interval; SRMR = Standardised Root Mean Square Residual; SE = Standard Error; i = intercept; s = slope. The robust CFI, TLI and RMSEA, along with the scaled SRMR, are reported. The literature suggests considering a lower CI limit CI when interpreting the RMSEA [60,65]; conventionally, values of 0.05–0.08 indicate an acceptable fit, and values ≤ 0.05 indicate a very good fit [58]. Since some estimated variances were negative in the model utilizing the MLR as an estimator, the WLSMV was used as an alternative estimator. a A free-loading model was tested, with the following loadings for the slope: 0, a, b, 1. p = 0.06. * p < 0.05. ** p < 0.01. *** p < 0.001.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zrnić Novaković, I.; Streicher, A.; Ajduković, D.; Ajduković, M.; Kiralj Lacković, J.; Lotzin, A.; Lueger-Schuster, B. Trajectories of Adjustment Disorder and Well-Being in Austria and Croatia during 20 Months of the COVID-19 Pandemic. Int. J. Environ. Res. Public Health 2023, 20, 6861. https://doi.org/10.3390/ijerph20196861

AMA Style

Zrnić Novaković I, Streicher A, Ajduković D, Ajduković M, Kiralj Lacković J, Lotzin A, Lueger-Schuster B. Trajectories of Adjustment Disorder and Well-Being in Austria and Croatia during 20 Months of the COVID-19 Pandemic. International Journal of Environmental Research and Public Health. 2023; 20(19):6861. https://doi.org/10.3390/ijerph20196861

Chicago/Turabian Style

Zrnić Novaković, Irina, Alina Streicher, Dean Ajduković, Marina Ajduković, Jana Kiralj Lacković, Annett Lotzin, and Brigitte Lueger-Schuster. 2023. "Trajectories of Adjustment Disorder and Well-Being in Austria and Croatia during 20 Months of the COVID-19 Pandemic" International Journal of Environmental Research and Public Health 20, no. 19: 6861. https://doi.org/10.3390/ijerph20196861

APA Style

Zrnić Novaković, I., Streicher, A., Ajduković, D., Ajduković, M., Kiralj Lacković, J., Lotzin, A., & Lueger-Schuster, B. (2023). Trajectories of Adjustment Disorder and Well-Being in Austria and Croatia during 20 Months of the COVID-19 Pandemic. International Journal of Environmental Research and Public Health, 20(19), 6861. https://doi.org/10.3390/ijerph20196861

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop