Introduction

The loss of a loved one is a common and often highly stressful human experience. Although most people ultimately adapt to their loss [1], there is a sizeable minority of bereaved individuals who experience an enduring and disabling grief response [2]. This occurrence of maladaptive grief responses has been recognized in the most recent version of the International Classification of Diseases [ICD-11; 3] via the introduction of Prolonged Grief Disorder (PGD). Essential characteristics of PGD include (i) a history of bereavement, (ii) persistent and pervasive longing for, or preoccupation with, the deceased associated with intense emotional pain (e.g., sadness, guilt, denial), (iii) functional impairment, and (iv) symptoms persisting for an atypically long time relative to cultural norms [3].

Recently, the International Prolonged Grief Disorder Scale (IPGDS; 4) was developed as the first self-report measure of ICD-11 PGD symptoms. The IPGDS is a pan-cultural measure designed for the assessment of the ICD-11 PGD diagnosis. Despite being a novel instrument, emerging research provides support for the psychometric properties of the scale [4, 5]. This scale can generate scores to reflect PGD severity, but there is a lack of consensus regarding how best to identify levels of PGD that may require clinical help. Two approaches to identifying probable caseness have been proposed for the IPGDS, the ‘moderate’ and ‘strict’ approaches (hereafter referred to as ‘moderate PGD’ and ‘strict PGD’ respectively) [4, 5]. The required frequency of symptoms varies between the two approaches; for moderate PGD, symptoms must appear “sometimes”, “often”, or “always”, and for strict PGD, symptoms must appear “frequently”, or “always”. There has been limited research on the levels of ‘moderate PGD’ and ‘strict PGD’ as assessed using the IPGDS in samples of bereaved adults, with prevalence estimates ranging from 19.8 to 33.5% for moderate, and 6.9 to 37.8% for strict (in adult samples from Germany, China, Switzerland, and Portugal [4, 5]. To date, no study has examined the prevalence or correlates of moderate and strict PGD as measured by the IPGDS in a representative general population sample from the United Kingdom (UK).

It should be emphasized that various instruments have been developed for the assessment of ICD-11 PGD symptoms, including the extended version of the Prolonged Grief Disorder-13 (PG13 + 9; 6) and the Traumatic Grief Inventory-Self Report Plus (TGI-SR+; 7). The “strict” approach—also referred to as ‘liberal’ approach—is used by both the PG13 + 9 and the TGI-SR + , while the TGI-SR + also offers the option of a ‘conservative’ approach, which increases the number of symptoms necessary for diagnosis. Other studies have used proxy items derived from established measures to measure ‘strict’ and ‘conservative’ PGD [e.g., 8, 9]. Prevalence rates of probable ‘strict PGD’ have ranged from 1.5% in the German general population to 34% in a community sample of bereaved Danish adults [6, 7, 9]. Establishing prevalence rates of both probable ‘strict’ and ‘moderate’ PGD using different measures including the IPGDS is necessary to determine the diagnostic concordance between these measures. This is especially important considering how much emphasis the grief literature has placed on determining the best method for diagnosing PGD. There are currently five different criteria sets proposed within the literature, each of which generates diverging PGD prevalence rates [10]. Hence, more research is needed to determine the best way to implement the ICD formulation of PGD to ensure consistency in prevalence estimates across studies. Finally, it is important to consider the potential influence of the COVID-19 pandemic on prevalence rates. Emerging research indicates elevated levels of PGD among individuals bereaved during the COVID-19 pandemic [e.g., 11, 12], likely attributable to more unexpected deaths, changes to customary mourning and grief rituals, and lack of social support due to government restrictions [13].

This study aimed to report the levels of PGD symptoms, and estimate rates of moderate/strict PGD, in the UK general population and identify significant correlates. It was anticipated that correlates of PGD would be similar to those identified in other countries including sociodemographic factors such as gender [14], unemployment, low income, presence of children in the household [15], and being in a committed relationship [17]. It was anticipated that less time since passing [9, 14] and younger age of the deceased [16, 17] would be significant correlates of PGD. Finally, it was anticipated that mental-health factors would be significantly associated with PGD including depression [15,16,17], and loneliness [e.g., 18]. The present study also sought to look at the association between prior mental health treatment and PGD, however, no a priori hypotheses were developed due to the paucity of research in this area.

Methods

Participants

This study used data from wave five of the COVID-19 Psychological Research Consortium Study (C19PRC-UK), a longitudinal internet-based survey assessing the psychological and socioeconomic impact of the COVID-19 pandemic [19]. The C19PRC data have been discussed extensively elsewhere [20]. In brief, data for wave one of the C19PRC study was collected between the 23rd and 28th of March 2020 through the survey company Qualtrics. Quota sampling methods, a non-probability based sampling method whereby quotas of participants with certain characteristics are obtained, were used to ensure national representativity of the sample regarding age, sex, and household income. The final survey sample comprised of 2025 adults aged 18 years and above. Data collection for wave five of the C19PRC-UK study began on the 22nd of March 2021, approximately one year following the completion of the baseline survey [20]. Data collection was conducted in two phases: in phase 1 (24th March to 20th April 2021) adults who participated in Wave 4 of the survey (n = 3867) were contacted by Qualtrics either by email, SMS, or by in-app notifications and invited to participate further in the survey. In phase 2 (8th April to 20th April 2021) participants who had completed any other wave (i.e., waves one to three) were recontacted and invited to participate in wave five. Phase one of this fieldwork resulted in a total of 2377 participants from wave four completing wave five (61.5% recontact rate), while phase 2 resulted in an additional 143 participants from wave one to wave three re-entering the survey [19]. In total, there were 2520 participants for Wave 5. For the purposes of the present study, only those who were bereaved were eligible for inclusion. Specifically, respondents were asked “At any time in your life, has someone close to you died (e.g., a partner, parent, child, friend)?” and “When did the death occur?”. If the respondent indicated ‘yes’ to the former and more than six months ago to the latter, they were screened in to complete the IPGDS. Hence, the current study has a final sample size of 2025 participants. Ethical approval for the study was granted by the University of Sheffield (Ethical approval ref no. 033759). The wave five data used in the current study is available at: https://osf.io/ducgs.

Measures

ICD-11 Prolonged Grief Disorder

The International Prolonged Grief Disorder Scale (IPGDS; 4) is a fourteen-item self-report measure of ICD-11 PGD symptoms. The IPGDS contains two items measuring core symptoms (i.e., “longing or yearning” and “preoccupied with thoughts”) and an additional ten items measuring associated emotional symptoms including anger, avoidance, sorrow, and emotional numbness (see Table 1; items 3–12). One item measures functional impairment, while another measures the extent to which symptoms deviate from community or cultural norms. Participants rate symptom frequency during the past week using a five-point Likert scale ranging from ‘Not at all’ (1) to ‘Always’ (5). The reliability of the twelve symptom items of the IPGDS in the current study was α = 0.94.

Table 1 Endorsement rates for IPGDS items

Correlates

Bereavement timeframe

The IPGDS includes an additional item which asks participants to indicate when their loss occurred. The following response options are provided: (a) less than 6 months ago, (b) 6–12 months ago, (c) 1–5 years ago, (d) 5–10 years ago, (e) 10–20 years ago, and (f) more than 20 years ago. Because the current study included only those bereaved longer than six months (due to the requirement for symptoms to be present for ≥ 6 months in the ICD-11) the first category was redundant. The final two categories were collapsed into a single category because of low numbers of participants in each, leaving a total of four categories for the time since bereavement variable (1 = 6–12 months, 2 = 1–5 years, 3 = 5–10 years, 4 =  > 10 years).

Age of deceased

Participants were also asked “how old was the person when they died?”, with responses reported in years. A multi-categorical variable was created to represent the age group of the deceased (i.e., 1 = 0–17 years, 2 = 18–34 years, 3 = 35–59 years, 4 = 60–79 years, 5 = 80 years and above).

Demographics

Participants indicated their sex (0 = male, 1 = female) and age which was collapsed into four categories (18–34 years; 35–44 years; 45–54 years; 65 years and over).

Living location

Area (city, suburb, town, rural) and country (England, Scotland, Wales, Northern Ireland) of residence were assessed.

Relationship status

Participants were asked to indicate their current relationship status using six categories (single - not currently in a committed relationship; in a committed relationship but not living together; cohabiting; married; in a civil partnership; single- never been in a committed relationship). Relationship status was recoded into a binary variable for the purposes of the current study (0 = not in a committed relationship, 1 = in a committed relationship).

Religious identify

Participants were asked to indicate their religious identity and were provided with eleven response options: Atheist, Agnostic, Catholic, Protestant, Sunni, Shia, Jewish, Buddhist, Sikh, Hindu, Other. The data was recoded into a binary variable representing religious status [0 = not religious (i.e., Atheist, Agnostic), 1 = religious (i.e., all other categories)].

Parental status

Participants were asked to identify their parental status; categories included: (1) I do not have any children; (2) I have a child/children under 18 years of age, and he/she/they primarily live with me in my household; (3) I have a child/children under 18 years of age, but he/she/they primarily live elsewhere; (4) I have a child/children aged 18 years or over, and he/she/they primarily live with me in my household; (5) I have a child/children aged 18 years or over, but he/she/they primarily live elsewhere, and (6) someone else’s child/children under 18 years of age live with me in my household. The data were recoded into a binary variable to represent parental status (0 = child/children 18 years or below in the household (i.e., 1, 3, 4, 5), 1 = child/children 18 years or above in the household (i.e., items 2, 6).

Finances and employment

Participants were asked to indicate their self-estimated gross annual income for 2019 (1 = £0–15,490; 2 = £15, 491–£25, 340; 3 = £25,341–£38,740; 4 = £38,741–£57,930; 5 = £57,931 or more) and if they were currently in receipt of any government benefits excluding child support and state pension (0 = not in receipt, 1 = in receipt). Current employment status was initially measured using 10 categories and was recoded into a binary variable representing economic activity (0 = not economically active (unemployed, retired, full-time student), 1 = economically active (employed full or part-time, self-employed full or part-time, furloughed due to COVID-19 Pandemic).

Mental health

Depression

Symptoms of major depressive disorder (MDD) were assessed using the PHQ-9 [21]. Participants were asked to indicate how often, over the last two weeks, they had been bothered by each of the nine symptoms contained within the measure. Responses ranged from ‘not at all’ (0) to ‘nearly every day’ (3). Possible PHQ-9 scores ranged from 0 to 27, with higher scores indicative of higher levels of depression symptomology. Scores of ≥ 5, ≥ 10, and ≥ 15 indicate mild, moderate, and severe levels of depression, respectively (13), with the recommended cut-off score of ≥ 10 used in the current study to indicate meeting caseness criteria [22]. This cut-off score has been shown to have adequate sensitivity (0.85) and specificity (0.89) for detecting cases of MDD [23]. The psychometric properties of the PHQ-9 scores have been widely supported [23] and the internal consistency of the scale in this study was α = 0.92.

Anxiety

The GAD-7 [23] was used to measure symptoms of generalized anxiety disorder (GAD). Participants were asked to indicate, using a 4-point Likert scale ranging from ‘not at all’ (0) to ‘nearly every day’ (3), the extent to which they were bothered by each of the seven anxiety symptoms contained within the GAD-7 over the past week. Possible GAD-7 scores range from 0 to 21, with higher scores indicating higher levels of anxiety symptomology. For the current study, the recommended cut-off score of ≥ 10 was used to indicate clinically significant anxiety symptoms as this cut-off score has been shown to have adequate sensitivity (0.89) and specificity (0.82) for detecting cases of GAD [24]. The validity and the reliability of the GAD-7 as a measure of anxiety has been supported for use in the general population [24]. The internal consistency of the GAD-7 in the present study was α = 0.95.

Mental health treatment seeking

Participants were presented with the following statement: “mental health difficulties are very common. It will help us understand our survey results if you would tell us whether you currently or have in the past received treatment (medication or talking therapies) for these kinds of difficulties.” Participants were required to tick all options that applied to them which included: (1) I have never received treatment for mental health problems, (2) I have received treatment for mental health problems in the past, (3) I am currently receiving treatment for mental health problems, (4) I am currently receiving treatment for mental health problems but it has been cancelled temporarily due to the lockdown, (5) I am currently on a waiting list to receive treatment for a mental health problem, and (6) Prefer not to answer”. A binary variable was created to represent current mental health help seeking [0 = no (1, 2, 6), 1 = yes (3, 4, 5)].

Loneliness

Loneliness was measured using the Loneliness Scale [25], a three-item measure specifically designed for use in large-scaled population surveys. Respondents were asked how often they felt (no time frame specified): (1) that they lacked companionship, (2) left out; and (3) isolated from others. Responses were scored on a three-point Likert scale ranging from ‘hardly ever’ (1) to ‘often’ (3). Possible scores ranged from 3 to 9, with higher scores indicating higher levels of loneliness. A cut-off score of 6 or more was used to indicate loneliness in the present study, as has been done in prior research [e.g., 26]. The internal consistency of the Loneliness Scale in the current study was α = 0.90.

Analytic procedure

Diagnostic algorithms

Although the ICD-11 lists essential requirements for PGD it does not specify the number or severity of symptoms required to meet the clinical criteria. Two diagnostic algorithms have been proposed to date [4, 5]. PGD requires the (1) endorsement of one or more core symptom, (2) one or more emotional symptom, (3) significant functional impairment, and (4) the grief response persisting beyond cultural or community norms. For the IPGDS the ‘strict’ requirement is met when the items are rated ≥ 4 (i.e., ‘often’, ‘always’) on the Likert scale, and ≥ 3 (i.e., ‘sometimes’, ‘often’, ‘always’) for ‘moderate’.

Statistical analysis

The frequencies of PGD symptom-endorsement and the proportions of participants meeting requirements for ICD-11 PGD according to the ‘strict’ and ‘moderate’ algorithms were calculated. Next, a tripartite variable was created representing whether a participant met probable caseness using the (1) strict or (2) moderate criterion, or (3) did not meet probable caseness for either. Chi-square tests of association were conducted to examine the association between each correlate and this grouping variable, with the strength of these associations quantified using Cramer’s V (≤ 0.2 = weak, 0.2–0.6 = moderate, > 0.6 strong). Adjusted standardised residuals ≥ 1.96 for cells in the cross-tabulations were used to assess significant differences between the observed and expected counts. Following this, all multi-categorical correlates were dummy coded. The bivariate associations between each of the correlates variables and the three-level dependent variable were calculated using multinomial logistic regressions, with the ‘no criteria’ group serving as the reference category. Results are reported as odds ratios (ORs) with 95% confidence intervals. Finally, all correlates were then entered simultaneously in the logistic regression to examine the effects of each correlate while controlling for all other potentially influential correlates, with results reported as adjusted ORs (AOR) with 95% confidence intervals; this was the multivariate model.

Results

Symptoms and levels of PGD

Of the 2,025 participants, 1771 (70.5%) had reported having lost someone more than 6 months ago. The frequencies of PGD symptom-endorsement are presented in Table 1.The most commonly endorsed PGD symptoms at both the moderate and strict level were “longing or yearning for the deceased” (moderate = 46.8%; strict = 12.9%), “intense feelings of sorrow, related to the deceased” (moderate = 42.1%; strict = 15.5%), “lost a part of myself” (moderate = 40.3%; strict = 16.5%), and “preoccupied with thoughts about the deceased or circumstances of the death” (moderate = 29.4%; strict = 8.4%). The least commonly endorsed PGD symptoms at both the moderate and strict level were “grief significantly interferes with my ability to work, socialize or function in everyday life” (moderate = 13.9%; strict = 4.5%), “my grief would be considered worse than for others from my community or culture” (moderate = 14.6%, strict = 5.5%) andI blame others or the circumstances for the death” (moderate = 14.0%, strict = 4.7%). In this sample, 7.9% (95% CI 6.6, 9.1: n = 140) met probable caseness using the moderate criteria, and 2.4% (95% CI 1.7, 3.1: n = 43) met probable caseness using the ‘strict’ criteria for ICD-11 PGD.

Chi square tests of association results

Chi-square tests of association were conducted to examine the bivariate relationship between the correlates and diagnostic status for those participants who were bereaved for longer than six months (see Table 2). All correlates were significantly associated with diagnostic status except for gender, religious status, country, and area of residence. For the sociodemographic variables the largest effect sizes were for being in receipt of benefits (V = 0.14) and Age (V = 0.13: significantly more 65 + years in the moderate and strict criteria groups), and these effects were small. The effects sizes for the psychosocial variables were larger with participants screening positive for anxiety (V = 0.27), depression (V = 0.32), loneliness (V = 0.23), and mental health help-seeking (V = 0.16) all more likely to be in the moderate or strict PGD groups. The effect sizes were generally moderate.

Table 2 Crosstabulation of diagnostic status and sociodemographic, loss-related and mental health factors

Regression results

The bivariate association between each of the correlates and group membership were examined (see Table 3). The reference group for all analyses was ‘no criteria met’ group. When the unadjusted ORs were calculated, significant positive correlates of moderate PGD included lower income, receiving government benefits, being economically active, living in a city or town, living in Northern Ireland, younger age, having a child/children under 18 years living in the household, younger age of deceased, and all mental health variables. Regarding strict PGD membership, significant positive correlates in the unadjusted analyses included earning between £0–£15,490 per annum, receiving government benefits, younger age, younger age of deceased, and all mental health variables. Significant negative correlates of strict PGD included being in a committed relationship.

Table 3 Correlates of ICD-11 PGD diagnostic status (unadjusted and adjusted) for bereaved for greater than six months sample

The multivariate association between the correlates and group membership were then examined (see Table 3). The multinomial logistic regression model of diagnostic status was statistically significant (χ2 (64) = 296.91, p < 0.001). Positive correlates of moderate PGD included low income (two lowest income categories compared to highest), living in Northern Ireland (compared to living in England), living in a town (compared to rural), being aged 18–34 years (compared to ≥ 65 years), identifying as religious (compared to not religious), age of deceased being 0–17 years (compared to ≥ 80 years), less time since bereavement (compared to 10 + years since bereavement), depression (compared to not meeting caseness criteria), and loneliness (compared to not meeting caseness criteria). Correlates of strict PGD included religiosity (compared to not religious), younger age of the deceased (i.e., 0–17 years, 18–34 years, 35–59 years compared to ≥ 80 years), less time since the bereavement (i.e., 6–12 months, 1–2 years compared to ≥ 10 years), anxiety (compared to not meeting caseness criteria), and depression (compared to not meeting caseness criteria).

Discussion

The current study had two primary aims: (1) report the levels of PGD symptoms, and estimate rates of moderate/strict PGD, and (2) identify sociodemographic, loss, and mental health related correlates of ICD-11 PGD (moderate and strict) using data from a large, nationally representative sample of UK adults.

Levels of PGD symptoms and rates of ICD-11 PGD

Results demonstrated that common features of PGD in the UK population (at the moderate and strict level) include longing and preoccupation with the deceased, intense sadness, and loss of identity. A substantially smaller proportion of bereaved adults reported functional impairment associated with their PGD symptoms or persistence of PGD symptoms for an atypically long time relative to cultural norms, both of what are essential features of ICD-11 PGD. These findings are not surprising given that the symptoms comprising ICD-11 PGD are not considered inherently pathological but rather it is the persistence of symptoms and their ability to interfere with daily functioning which separates ICD-11 PGD from a “normal” grief response [27]. These findings are similar to prior studies where a relatively small albeit significant proportion of bereaved individuals meet criteria for probable diagnosis of PGD [e.g., 2, 5, 6, 28]. Hence, high endorsement of functional impairment and the cultural criterion would not be expected given that only a small proportion of bereaved individuals are postulated to experience an enduring and disabling grief response [2]. Findings suggested that the prevalence of ICD-11 PGD was between 2.4 and 8%. The prevalence of moderate PGD (i.e., 8%) fell within the lower range of strict PGD rates observed in prior studies [4, 5]. This may have been related to the small samples used in other studies and/or the sampling methods (which included the recruitment of participants from bereavement support groups [4, 5]). The proportion of participants who met probable caseness for PGD using the ‘strict approach’ was slightly higher than that observed in a general population sample of German adults [6], suggesting that different measures may produce divergent prevalence rates. That being said, the probable diagnostic rates observed in the present study align with a recent systematic review suggesting that 6.8 to 14.0% experience symptoms consistent with PGD [2]. However, the variability in prevalence rates across studies underscores the importance of determining where a threshold for ICD-11 PGD caseness should be set. In ICD-11, PGD belongs to the diagnostic category of stress-related disorders, alongside Posttraumatic Stress Disorder (PTSD), Complex PTSD (CPTSD), and Adjustment Disorder (AjD) [3] and each require an external stressor for diagnosis (a traumatic event in the case of PTSD and CPTSD, a psychosocial stressor or multiple stressors in the case of AjD, and the loss of a loved one in the case of PGD). The thresholds used in the self-report measures developed to assess these other stress-related disorders—the International Trauma Questionnaire [ITQ; 28] for ICD-11 PTSD and CPTSD and the International Adjustment Disorder Questionnaire [IADQ; 29] for ICD-11 AdJ—require that symptoms be present at least sometimes (i.e., a score of ≥ 3 on Likert scale of symptom frequency). These thresholds align to the ‘moderate’ algorithm for the IPGDS. Hence, it could be argued that the diagnostic algorithm employed to determine caseness of ICD-11 PGD should mirror these other stress-related disorders. Given that almost half (i.e., 40.6%) of the present sample reported their loss as occurring more than five years ago and that research has shown more recent losses to be associated with higher rates of PGD [e.g., 9] further research is required to understand the extent of PGD among recently bereaved adults.

Sociodemographic, loss-related, and mental health correlates of PGD (moderate and strict)

Findings from the bivariate analyses highlighted several sociodemographic correlates of both moderate and strict PGD including low income, receiving government benefits, younger age of bereaved, and not being in a committed relationship. After adjusting for all other correlates, only low income and younger age of the bereaved remained significant correlates of moderate but not strict PGD. Specifically, participants in the two lowest income categories were over twice as likely to meet criteria for PGD than those in the highest income category. These effect sizes were much larger than those observed in prior studies assessing the association between low income and CMDs where ORs ranged from 1.18, for those earning between £100–£200 pounds per week (equivalent to our lowest income category), to 0.96 for those earning between £400–£500 pounds per week (equivalent to our second lowest income category) [30]. Moreover, the greater likelihood of moderate PGD among young adults aligns with the broader psychological literature evidencing a higher prevalence of mental disorders among young adults [31]. Contradicting the well-established gender gap in mental health [e.g., 32], no gender effects were observed in the current study. The absence of a gender effect may be considered in the context of a recent study investigating gender differences in trajectories of grief which found that although males and females were equally as likely to experience prolonged grief symptoms, males typically experienced acute PGD symptoms which abated with time, while females experienced worsening symptoms over time [33]. Hence, the commonly observed gender effects in mental health may not extend to PGD however, the disorder trajectory may differ in males and females. Other correlates of moderate PGD from the bivariate and multivariate analyses included living in NI as compared to England and living in a town as compared to rural areas. The former finding may be linked to “the Troubles” in NI, a thirty-year long period of political violence in NI which resulted in more than 3600 deaths [34], while the latter finding is consistent with literature indicating a lower prevalence of CMDs among rural inhabitants [35]. Overall, many of the sociodemographic correlates investigated were associated with increased risk of moderate but not strict PGD, suggesting that sociodemographic factors may be less predictive of more severe grief responses when considered in the context of other potentially more relevant factors such as loss-related and mental health factors.

Indeed, both bivariate and multivariate results demonstrated that loss-related correlates and psychological wellbeing were risk factors for moderate and strict PGD, with these effects being greatest for the latter group. Consistent with prior PGD research, bereavement recency was a risk factor for both moderate and strict PGD [9, 14]. Research has shown that time since trauma represents an important risk factor for PTSD and CPTSD [36], indicating that recency of the external stressor may represent a common risk factor among the stress-related disorders. Alternatively, data for this study was collected during the COVID-19 pandemic, and it is likely that some participants had experienced bereavements during the pandemic. There is a growing body of research illustrating a high prevalence of PGD among individuals bereaved due to COVID-19 [e.g., 11,12,13, 37], and hence further research is required to examine whether deaths that have occurred during the COVID-19 pandemic (both COVID-related and non-COVID related) have been associated with elevated prevalence of ICD-11 PGD. Consistent with prior research [e.g., 16, 17], younger age of the deceased was identified as a risk factor for moderate and strict PGD, and these effects were greatest for those reporting the death of a child. Similar to the death of a child, the death of a young or middle-aged adult is often unexpected and untimely which may lead to more maladaptive grief responses. It was interesting that religiosity (as compared to non-religiosity) was associated with increased risk of moderate and strict PGD, especially given that some prior research has shown those with strong spiritual beliefs demonstrate more adaptive mental health responses to grief [e.g., 38]. Further research is required to unpack the association between religiosity and ICD-11 PGD.

Finally, in terms of mental health correlates, depression, anxiety, loneliness, and mental health help-seeking were associated with both moderate and strict PGD in the unadjusted analyses. After controlling for other correlates, depression and loneliness remained significant correlates of moderate PGD while depression and anxiety remained significant correlates of strict PGD. Previous factor analytic work has shown that PGD, depression, and anxiety represent distinct yet highly correlated constructs [e.g., 39], and that the co-occurrence of these disorders is not uncommon. Similarly, there are several commonalities among the symptoms captured by the Loneliness Scale (i.e., lacking companionship, feeling left out, isolation from others) and the IPGDS (e.g., “I feel that I lost a part of myself”, “I have difficulties engaging in activities”, “grief significantly interferes with my ability to work, socialize, or function in everyday life”). Interestingly, anxiety uniquely predicted strict PGD while loneliness uniquely predicted moderate PGD in the multivariate analyses. This suggested that anxiety may be more prevalent among those with more severe grief symptoms while loneliness may be more relevant to those with moderate symptoms of grief. Future research may benefit from unpacking these associations, particularly given that understanding comorbidities associated with PGD remains an important research priority [40]. It is well-established that anxiety is a common response to bereavement due to the separation from a significant other, confrontation with one’s mortality, and exposure to stressors including financial adversity [41]. Indeed, research has shown how high levels of comorbidity between PGD and adult separation anxiety disorder (ASAD) [42] and generalized anxiety disorder [43]. It should be noted that intense loneliness is acknowledged as a PGD symptom in the text version of the DSM-5 [44] and as such, may be regarded a feature rather than a correlate of PGD. Whether intense loneliness should be listed as a symptom of PGD in the ICD-11 is a matter which warrants investigation.

Strengths and limitations

Findings from this study should be considered in terms of several limitations. First, the cross-sectional nature of the data does not allow for conclusions to be drawn regarding causality. Second, the use of self-report measures are subject to several biases including misinterpretations of items and response categories [45]. Nevertheless, it is important to highlight that a clinician administered tool is not yet available. Third, although we aimed to include a wide range of sociodemographic, loss-related, and mental health correlates, it is likely that other potentially relevant indicators were excluded. Because the C19PRC assesses a broad range of social, economic, and mental health variables in the context of the COVID-19 pandemic, a thorough examination of grief-related variables was not possible. For instance, relationship to the deceased and cause of death has been shown to be important correlates of PGD [e.g., 14,15,16], which were not assessed within the C19PRC survey. Fourth, it should be noted that the current study did not investigate the psychometric properties of the IPGDS and hence, the validity of this measure within this particular sample is uncertain. Further research is now necessary to examine the psychometric properties of the IPGDS within this sample, especially given the paucity of such studies. Finally, only a small proportion of participants met probable caseness for 'strict PGD' and several categories of the correlates examined were relatively small (e.g., proportion of participants living in NI).This should be taken into account when evaluating the results as it is possible that it led to unstable estimates.

Conclusion

In conclusion, this study is the first to examine the prevalence and correlates of ICD-11 PGD as measured via the IPGDS in a representative general-population sample of bereaved adults. This study highlights that approximately one in ten bereaved UK adults exhibit moderate levels of prolonged grief symptoms (i.e., moderate PGD) while one in forty exhibit high levels of prolonged grief symptoms (i.e., strict PGD). Considering the high prevalence of PGD in the general population, we recommend routine screening for PGD following a bereavement. Further exploration of the role of income, religiosity, area and country of residence, age of deceased, and mental health factors in predicting maladaptive grief responses is necessary. The variability in prevalence rates and risk factors when using the moderate and strict diagnostic algorithms emphasizes the need for an agreed upon expert consensus as to which algorithm to carry forward. This is especially imperative from a clinical perspective where prognosis and selection of clinical interventions are determined by diagnosis [46]. Future studies may benefit from determining the clinical meaningfulness of moderate PGD symptoms as compared to strict PGD symptoms to ensure that the algorithm employed to determine PGD caseness does not discount bereaved individuals in need of clinical intervention.