Introduction
Job stress in police officers and its consequences is attracting the attention of the scientific community both from a descriptive and preventive perspective, as well as from the point of view of intervention (Baek et al., 2021; Bhowmick and Mulla, 2021; Foley and Massey, 2021; Luceño-Moreno et al., 2021; Márquez et al., 2021). Guaranteeing fundamental rights, optimising citizen security, and protecting private property are activities that require considerable psychological effort on the part of police officers. These activities, given the high degree of responsibility involved, tend to generate situations that can trigger stress responses in these professionals. For example, having to work urgently, not having sufficient human or material resources, having total availability for police service, having to work night or weekend shifts and the unpredictability of the work are aspects that have been identified as psychosocial factors that increase the level of stress in police officers (Cieślak, 2020; Díaz et al., 2020; Sherwood et al., 2019).
These working conditions, as well as involving aspects related to social interaction, have to do with job-related aspects. Among the factors related to job-related elements, we can highlight the content, the workload, the working hours, or the place of work. While the more social aspects include organisational culture, interpersonal relationships, the employee's role in the organisation or career development (Leka et al., 2008). These psychosocial factors can play both a positive and a negative role in job performance. When their effect is negative and can impact on the physical and psychological health of workers, they are called psychosocial risk factors (Cox and Griffiths, 2003). Some professions are more likely to be exposed to such factors. This is the case for public safety professionals such as police officers or healthcare workers (Teoh et al., 2020, 2022).
In police work, the presence of psychosocial risk factors has been found to increase the perception of job stress and have an impact on both job satisfaction and performance (Adams and Mastracci, 2018; Galbraith et al., 2021; Garbarino et al., 2013; Hickman et al., 2011; Hoven et al., 2015; Luceño-Moreno et al., 2016; Mushwana et al., 2019; Violanti et al., 2017). Recently, Talavera et al. (2021) added fear of severe acute respiratory syndrome (SARS-Cov-2) to the list of psychosocial risk factors in the police profession.
As previously mentioned, one of the consequences of exposure to psychosocial risk factors is work-related stress. The classic definition provided by Lazarus and Folkman (1984) states that stress is the result of an appraisal process in which the individual perceives that there is an imbalance between the demands of his or her environment and the resources that the individual considers he or she has to cope with these demands. The theory of job demands and resources (Bakker and Demerouti, 2007, 2013; Bakker and Vries, 2021) has recently been used as a model for understanding burnout as a function of job stress. Some work has shown that this model is useful for understanding fatigue, stress evolution over time, perceived organisational support, mental health, and organisational commitment in police officers from different parts of the world (Hu et al., 2017; Lan et al., 2020; Maran et al., 2020; Raper et al., 2020; Zeng et al., 2020; Zhou et al., 2020). From this theoretical framework it is possible to predict job characteristics (demands and resources) and the consequences of chronic job stress (e.g., absenteeism and job performance). For this model, excessive job demands are the predictors of negative consequences on workers' health, while resources are considered the antecedents of job satisfaction, motivation, and commitment. Thus, when a person perceives that demands are greater than the resources he or she must cope with a certain situation, an episode of stress occurs. When workers continuously experience work-related stress, it can disrupt the smooth functioning of organisations (Queirós et al., 2020), as well as having repercussions on physical and mental health (Kelley et al., 2019; Magnavita et al., 2018).
The concept and definition of burnout has evolved over time to adapt to the working conditions that have characterised different eras (Salanova and Llorens, 2008). According to Maslach and Leiter (2016), Burnout Syndrome "is a psychological syndrome emerging as a prolonged response to chronic interpersonal stressors on the job" (p. 103). This syndrome has been recognised by the World Health Organisation (WHO) as a factor affecting health (WHO, 2022) in the 11th revision of the International Classification of Diseases (ICD-11). In the Spanish version of the classification of diseases it is called Occupational Burnout Syndrome (code QD85) and it is the product of chronic exposure to occupational stress that has not been successfully managed. A recent international study has managed to reach a consensus definition of burnout considering the harmonised criteria of 50 experts from 29 countries. This work suggests that burnout is simply " an exhaustion due to prolonged exposure to work-related problems" (Guseva Canu et al., 2021, p. 104).
Various contributions in recent years (Maslach, 2017; Maslach and Leiter, 2016, 2017; Schaufeli, 2017) have confirmed the relationship between the characteristics of public service, or direct work with citizens, and burnout syndrome. In members of police forces, the syndrome manifests itself in the form of emotional exhaustion, depersonalisation, physical discomfort, insomnia, depressive symptoms, inappropriate consumption of psychotropic/antidepressant medication, or hospitalisation due to the presence of mental disorders (McCarty et al., 2019; Salvagioni et al., 2017). In this group of workers, psychosocial risk factors have the same harmful impact as work-related stress, generating high scores in emotional exhaustion, cynicism, depersonalisation, and low levels of professional efficacy (De la Fuente et al., 2013; McCarty et al., 2019; Padyab et al., 2016; Talavera et al., 2021). Burnout syndrome in police officers has also been studied in relation to certain sociodemographic variables such as marital status or gender (Aguayo et al., 2017; Galbraith et al., 2021; Mushwana et al., 2019; Roberts, 2013). In terms of marital status, several studies indicate that having a partner is a protective factor (Aguayo et al., 2017; Roberts et al., 2013). On the other hand, those in a situation of separation or without a partner showed lower professional efficacy, higher levels of emotional exhaustion, and greater distancing from both work and colleagues. In relation to gender, García et al. (2013) did not observe differences between men and women in levels of work-related stress; however, according to Acquadro et al. (2015), there are high levels of stress in women for organisational reasons and in dealing with violent situations compared to their colleagues.
The present research aims to identify those psychosocial factors that pose a higher risk for police officers. This can be considered as a first step to contribute to the development of appropriate interventions for the prevention of occupational risks derived from high levels of burnout (Galanis et al., 2021). For this purpose, networks of psychosocial risk factors in police officers will be generated. Network analysis is a contemporary application of classical graphical models (Ruiz-Ruano and Puga, 2020) and in the last decade it has received special attention from psychology as an appropriate method to study phenomena of a pathological nature (Borsboom et al., 2011; Cramer et al., 2010; Schmittmann et al., 2013). From this theoretical perspective, the mental disorder is the result of the causal interaction that occurs between the symptoms of a network (Fonseca-Pedrero, 2017, 2018). The aetiology of mental disorder is seen as the result of the increasing interconnectedness of symptoms in networks or statistical graphical models (Borsboom, 2017). In practical terms this means that loosely connected symptom networks highlight the absence of psychological disorders, while highly connected symptom networks would relate to psychopathological problems. Given that, as mentioned above, the WHO (2022) recognises burnout as a syndrome of psychological symptoms occurring in the work context, the relationships between psychosocial risk factors and the components of burnout will be modelled. In other words, networks of psychosocial risk factors will be estimated based on the characteristic profile of burnout (high emotional exhaustion, high cynicism, and low professional efficacy) in order to assess the relationship between these variables. Therefore, it is hypothesised that the psychosocial risk factors appear disconnected, or mostly disconnected, in police officers without a burnout profile. On the other hand, it is hypothesised that more widely connected symptom networks are generated when police officers show a characteristic burnout profile. The aim of this work is, therefore, to use networks of psychosocial risk factors to identify those elements that play a more significant role in the occurrence of burnout. In addition, this research aims to further study the relationship between sociodemographic variables and the risk of suffering burnout symptoms in local police officers. Although this police force has specific characteristics conditioned by its field of action, the factors studied in this research could be compared with those present in the National Police Force or other State Security Forces and Corps.
Method
Participants
The sample consisted of 323 participants (14.6% of a total of 2211), all of them members of the Local Police forces of the Region of Murcia. Eighty-seven percent of the sample were men, of whom 72.9% were married or living with a partner, 16% said they were single and 11% said they were widowed. Meanwhile, 76.2% of the women were married or cohabiting, 19% were single and 4.8% were widowed. Their ages ranged from 20 to 60 years (M = 42, SD = 8.32). Most of the participants had completed secondary education (48.3%), 41.2% of the sample had a university degree or professional specialisation, 3% of the participants had only completed primary education, and only 7.4% had a master's degree or doctorate. In terms of professional category, 242 were officers, 77 were sub-inspectors, three were inspectors and one was a commissioner. Participants' length of service as police officers ranged from two to 37 years, and the mean length of service of the sample was 16.42 years (SD = 7.79, Mdn = 13). Participants were accessed through a non-probability sampling procedure based on accessibility and therefore potential participants were not segmented according to hierarchical levels or organisational divisions.
Instruments
A form was used to collect information on socio-demographic variables (age, sex, marital status, level of academic studies, place of work, seniority, and hierarchical scale within the police force), and the questionnaires F-PSICO 4.0 for the assessment of psychosocial risk factors and the Maslach Burnout Inventory-General Survey for the assessment of burnout.
F-PSICO
For the assessment of psychosocial risk factors, the F-PSICO in its 4.0 version was used, which is edited by the National Institute for Safety and Health at Work (Instituto Nacional de Seguridad y Salud en el Trabajo, INSST, 2018). This instrument consists of 44 questions, some of them multiple, so that the number of items amounts to 89. It provides information on nine psychosocial factors: working time (WT): referring to the structure and order of working time as well as rest periods, autonomy (AU): understood as the responsibility, capacity and possibility of the worker to regulate his/her work activity as well as to make decisions, mental workload (W): referring to the amount of tasks that the worker has to face in his/her job, psychological demands (PD): understood as the cognitive demands derived from the task, variety/content (VC): referring to the perceived social usefulness of the job, participation/supervision (PS): referring to the degree of control to which the worker is subjected and his or her involvement in the work activity, interest in the worker/compensation (IWC): referring to the degree to which the company offers promotion or career development incentives to the worker, role performance (RP): referring to the clear definition of the functions required in the job, relationships and social support (RSS): referring to the aspects of social interaction circumscribed to the job.
For each of the factors, McDonald's ω coefficient was obtained to estimate the reliability of the instrument understood as internal consistency (Prieto and Muñiz, 2000). In this study the estimates have been considered satisfactory for the dimensions of this instrument, WT (ω = .62, 95% CI [.53, .69]), AU (ω = .91, 95% CI [.89, .92]), W (ω = .7, 95% CI [.63, .73]), PD (ω = .78, 95% CI [.78, .85]), VC (ω = .65, 95% CI [.51, .64]), PS (ω = .86, 95% CI [CI95% = .85, .9]), IWC (ω = .83, 95% CI [8, .86]), RP (ω = .88, 95% CI [.87, .9]) and RSS (ω = .7, 95% CI [.64, .9]).
MBI-GS (Maslach Burnout Inventory-General Survey)
In this study we used the Spanish adaptation of the scale developed by Salanova et al. (2000). The response scale used is a 7-point Likert-type scale with which the person has to indicate the degree of agreement with the items presented: (0) never; (1) almost never/few times a year; (2) sometimes/once a month or less; (3) regularly/few times a month; (4) quite often/once a week; (5) almost always/few times a week; and (6) always/every day. The fifteen items that make up the test are grouped into three dimensions that assess the degree to which the person is burned out at work: a) emotional exhaustion, which is focused on on physical and emotional fatigue, b) cynicism, from which information is collected about indifference and detachment from one's work, and c) professional efficacy, which is centred on the subjects' expectations of success.
For the definition of the burnout profile, the cut-off points specified in the scales proposed by Salanova et al. (2000) were considered. Based on these cut-off points, categories of low, medium, and high scores were created for the dimensions of emotional exhaustion, cynicism, and professional efficacy. A person was considered to show a profile compatible with burnout when he/she obtained high scores for emotional exhaustion and cynicism while presenting low scores for professional efficacy. The internal consistency estimates for emotional exhaustion (ω = .93, 95% CI [.92, .94]), cynicism (ω = .91, 95% CI [.9, .93]) and professional efficacy (ω = .89, 95% CI [.86, .9]) are considered acceptable for the purposes of this research.
Procedure
Data were collected through an online questionnaire that was administered between March and December 2019. Before answering the questionnaire, participants were informed that the study aimed to determine which psychosocial factors have the greatest impact on work-related stress in local police officers and the assessment of burnout. Participation was voluntary and responses to the electronic form were anonymous. The answers given in the forms were analysed in aggregate form to maintain confidentiality and their management was in accordance with the ethical guidelines of the UCAM Catholic University of Murcia. The data collection procedure, as well as its subsequent analysis, was governed by the ethical standards and good practice guidelines suggested by the American Psychological Association (2020). The only exclusion criterion was not belonging to any Local Police Force in the Region of Murcia. The responses were manually transferred to the computer application F-PSICO version 4.0 (INSST, 2018). Four participants were discarded due to the application of the exclusion criterion.
Data analysis
To estimate the network models, the EBICglasso algorithm (Epskamp et al., 2018; Fonseca-Pedrero, 2018) was used to apply a regularisation method to the models (Lever et al., 2016). The algorithm was run on the correlation matrices to generate networks with weighted and signed (positive or negative) links. Raw centrality measures were estimated and the λ (tuning parameter) of the algorithm was set to the default value of 0.5. The expected influence measures of each node in each network were estimated and interpreted as neat impact of each variable on the network. Expected influence measures are statistics derived from classical measures of centrality in graphical models (Borsboom and Cramer, 2013) and high values indicate higher impact of the variable on the network model. In practical terms these statistics can be used to identify target nodes or variables for psychological intervention. Statistical analyses were conducted with SPSS statistical software version 24 and JASP version 0.16.4 (JASP Team, 2022).
Results
36.2% of the participants showed a response profile compatible with burnout. This percentage of agents fulfils the criteria of risk of suffering burnout considering the questionnaire used (high scores in emotional exhaustion and cynicism, and low scores in professional efficacy). In addition, 53.6% of the sample showed high scores on emotional exhaustion and cynicism. No statistically significant differences were observed in the burnout profile as a function of age (rbp = .02, p = .75, bilateral contrast). There is also no relationship between gender and burnout profile (OR = 0.54, 95% CI [0.26, 1.13], p = .14). This suggests that there is no difference in the proportion of people with burnout when the sample is segmented by men and women. However, there are differences according to marital status. While 58.5% of single people show high levels of burnout (χ2 = 11.21, p <.001, OR = 1.41, 95% CI [1.21, 1.64]), only 30% of married participants show symptoms compatible with the syndrome (χ2 = 74.30, p <.001, OR = 0.43, 95% CI [0.36, 0.50]). Also, it should be noted that there is no significant relationship between job seniority (rbp = .02, p = .67) and burnout.
Regarding the relationship established between the rest of the variables studied and the sex of the participants, as can be seen in Table 1, statistically significant differences were obtained in several factors. Women scored significantly higher on the following variables: working time, workload, psychological demands and relationship and social support. However, no statistically significant differences were observed between men and women in the levels of emotional exhaustion, cynicism, and professional efficacy. When analysing the percentage of men and women who score high on the emotional exhaustion component of the scale, no statistically significant differences are observed (χ2 = .76, p = .38, OR = .75, 95% CI [0.37, 1.51]). Neither is being male or female significantly related to scoring low on the self-efficacy scale (χ2 = 2.22, p = .14, OR = 0.61, 95% CI [0.30, 1.24]). However, there was a statistically significant relationship between being female and scoring high on the cynicism scale (χ2 = 4.39, p = .28, OR = .47, 95% CI [0.21, 0.97]).
Table 1. Psychosocial risk factor scores and burnout dimensions as a function of gender.

Note. WT: working time, AU: autonomy, W: workload, PD: psychological demands, VC: variety/content, PS: participation/supervision, IWC: interest in work/compensation, RP: role performance, and RSS: relationships and social support.
*p < .05,
**p < .01, and
***p < .001; degrees of freedom: 321 in all comparisons.
As can be seen in Figure 1, as the levels of cynicism and emotional exhaustion increase in the police officers who participated in the study, the sparsity of the networks decreases and the number of links between psychosocial risk factors increases. That is, the higher the level of cynicism and emotional exhaustion in police officers, the higher the degree of association between psychosocial risk factors. Cynicism is the dimension most clearly affected by this phenomenon. While low levels of cynicism produce an unconnected network (zero links and sparsity equal to 1), intermediate levels produce a connected network (sparsity = .53, 17 links) and high levels of cynicism produce a network with the highest number of links (22) and the least sparsity (.39). Similarly, low levels of emotional exhaustion produce a sparser network (sparsity = .56) and containing fewer links (16). Whereas the degree of association between psychosocial factors increases at intermediate (sparsity = .5, 18 links) and high levels of emotional exhaustion (sparsity = .33, 24 links). In the case of professional efficacy this pattern is not so clear since both low and high levels of professional efficacy produce networks with the same sparsity (.47) and the same number of links (19).

Figure 1. Networks of psychosocial risk factors obtained for different levels of exhaustion, cynicism and professional efficacy.
Table 2 shows the coefficients of expected influence of each node in the psychosocial risk factor networks for the burnout dimensions. As can be seen, the three factors that are related to worrying levels of cynicism and professional efficacy are the same and in the same order of importance: autonomy, role performance and psychological demands. As can be seen in Figure 1, the expected higher influence of autonomy and role performance can be explained by the large number of connections that these factors have with the rest of the dimensions of the psychosocial factors scale. On the other hand, the high observed influence of psychological demands can be explained, despite its lower number of links with the rest of the nodes, by the very close association that this variable has with workload in both cases (see Figure 1). The three most influential factors for the emotional exhaustion network are role performance, autonomy and the dimension of relationships and social support. In this case, as can be seen in Figure 1, the greater degree of influence of these variables can be explained by the density of connections that these variables have with the rest of the nodes in the model.
Table 2. Expected influence (in absolute value) of each node on the estimated networks for high levels of emotional exhaustion, high cynicism and low professional efficacy.

Note: WT: working time, AU: autonomy, W: workload, PD: psychological demands, VC: variety/content, PS: participation/supervision, IWC: interest in work/compensation, RP: role performance and RSS: relationships and social support. The number of people in each group and their percentage of the total sample is shown in brackets.
Discussion
The results of this study suggest that high scores on emotional exhaustion and cynicism tend to generate more interconnected networks of psychosocial risk factors. This pattern of interconnectedness was not as clear for the professional efficacy dimension. Increased levels of interconnectedness in symptom networks have been considered as an indicator and explanatory phenomenon for the presence of psychological disorders (Borsboom, 2017; Borsboom et al., 2011; Cramer et al., 2010; Fonseca-Pedrero, 2017, 2018; Schmittmann et al., 2013). These results are consistent with the most recent version of the demands-resource theory (Bakker and Vries, 2021) as they show that certain psychosocial risk factors place high demands on police officers, and this may lead to an increased risk of occupational burnout.
On the other hand, the analysis of the expected influence of psychosocial risk factors in the network models for those who scored higher on the burnout dimensions can be used to design organisational/workplace policies aimed at minimising the risk of police officers experiencing problems arising from working conditions. The results suggest that autonomy, role performance and psychological demands are the variables most strongly related to high levels of cynicism and emotional exhaustion. Previous studies have argued in favour of these results. For example, the review by Violanti et al. (2017) suggested that there was a positive relationship between the psychosocial factors of working time, job autonomy and role performance with burnout. Other studies, albeit in a different social context, have also highlighted the relevant role of unclear police roles or reduced social support (Padyab et al., 2016) in generating high levels of stress and burnout in law enforcement agencies. Therefore, all this could be considered when designing and implementing measures aimed at optimising the psychological health of the people who form part of the Security Forces and Corps.
According to De la Fuente et al. (2013) 32.2% of the police officers showed a high level of burnout, 12.5% a medium level and 55.4% a low level. In this study, the results are very similar as 36.2% of the police officers obtained burnout scores considered worrying. The results of this study also point in the same direction as those reported by De la Fuente et al. (2013) on the relationship between psychosocial risk factors and burnout syndrome. However, in this study, a greater weight of the dimensions of emotional exhaustion and cynicism was observed.
About sociodemographic variables, the results suggest that marital status is related to the occurrence of burnout symptoms, as has been observed in previous research (Aguayo et al., 2017). Taken together, these results suggest that single participants, compared to participants who are in a relationship or married, exhibit higher levels of burnout symptoms such as emotional exhaustion and depersonalisation. This may be because that family has on police officers in terms of emotional support and sharing of family responsibilities. The data do not provide evidence of a positive relationship between length of service and high levels of burnout risk.
Regarding the difference between men and women, the results are in line with recent studies by Mushwana et al. (2019) and Galbraith (2021). In these studies, the dimensions of social support and workload were found to be important stressors for female police officers. However, while the work of Garcia et al. (2013) reported no significant differences in the incidence of stress between women and men, the results suggest such a positive relationship in several factors: time spent on daily work, the workload/volume of their shifts, excessive psychological demands, and lack of social support. The weight of these factors could be decisive in explaining the overall differences observed between men and women (Acquadro et al., 2015). Finally, this study also found a higher proportion of women than men with high levels of cynicism.
Considering the results of the study and the negative consequences of chronic stress on police officers (Adams and Mastracci, 2018; Galbraith et al., 2021; Luceño-Moreno et al., 2021; Mushwana et al., 2019; Violanti et al., 2017), some suggestions could be put forward so that the risks of law enforcement officers suffering from burnout problems could be minimised. The results of the study suggest that lack of autonomy, role-related ambiguities, excessive workload, high psychological demands, and low perceived social support from the organisation are linked to increased risk of occupational burnout in police officers. Therefore, actions to minimise the impact of these risk factors would improve the psychological health of police officers.
However, the present study suffers from certain limitations in relation to the participants in the study, both because of the difficulties of accessibility to the study population and because of the size of the sample. The influence of more socio-demographic variables (e.g., number of children, work quadrant, history of mood psychopathologies) that could be important in the impact of both psychosocial factors and the experience of stress and burnout have not been assessed and would be interesting to consider for future research. It should be borne in mind, as noted in the introduction, that the study participants belong to a specific type of police force, the local police force, which performs its duties in a specific place and under specific conditions. The results observed may not be generalisable to other police forces, such as the National Police, or to other State Security Forces and Corps. One of the key factors in this regard that has not been addressed in this study is the relationship of hierarchical dependence (e.g., the mayor's office or the State). Future work should shed light on these possible differences. In the sample used, men outnumber women, and this may have affected the results. The coping strategies that the agents use on a regular basis have not been evaluated either. Future work could explore this aspect further in order to identify efficient intervention strategies as suggested by Bakker and Vries (2021).
This research is pioneering in Spain with respect to the study of the relationship established between burnout and psychosocial risk factors that generate stress, considering the network analysis methodology used. Despite the characteristics of the sample used, it may be useful to try to implement preventive policies among security professionals to increase their well-being and productivity at work knowing the factors that condition the experience of stress and burnout (De la Fuente et al., 2013). The results of this study point to some of the factors that can be considered to minimise burnout in police officers. Therefore, human resources departments or personnel management sections of police forces can take advantage of these results to design healthier working environments for officers. In addition, it would be important that in cases where prevention is no longer possible, interventions to reduce both stress and burnout at work are implemented. Preventing and, secondly, mitigating burnout in these professionals will favour optimal levels of quality of life, which will result, on the other hand, in greater levels of well-being for society.










texto en 


