Work Engagement among Rescue Workers: Psychometric Properties of the Portuguese UWES

Rescue workers have a stressful and risky occupation where being engaged is crucial to face physical and emotional risks in order to help other persons. This study aims to estimate work engagement levels of rescue workers (namely comparing nurses, firefighters, and police officers) and to assess the validity evidence related to the internal structure of the Portuguese versions of the UWES-17 and UWES-9, namely, dimensionality, measurement invariance between occupational groups, and reliability of the scores. To evaluate the dimensionality, we compared the fit of the three-factor model with the fit of a second-order model. A Portuguese version of the instrument was applied to a convenience sample of 3,887 rescue workers (50% nurses, 39% firefighters, and 11% police officers). Work engagement levels were moderate to high, with firefighters being the highest and nurses being the lowest engaged. Psychometric properties were evaluated in the three-factor original structure revealing acceptable fit to the data in the UWES-17, although the UWES-9 had better psychometric properties. Given the observed statistically significant correlations between the three original factors, we proposed a 2nd hierarchal structure that we named work engagement. The UWES-9 first-order model obtained full uniqueness measurement invariance, and the second-order model obtained partial (metric) second-order invariance.


INTRODUCTION
Rescue workers' main tasks are to bring people out of danger after accidents, attacks, disasters, etc., are the first on the scene of an emergency or provide care to victims who are suffering. Firefighters, police officers or health professionals working in urgency or pre-hospital departments can be considered as rescue workers, having a highly risky and stressful occupation (Stassen et al., 2013). However, they are very motivated and engaged with their tasks frequently facing physical and emotional risks to help other persons (Shakespeare-Finch, 2011). To better understand work engagement among rescue workers, the theoretical framework is presented beginning with the work engagement concept and its study among rescue workers, followed by some considerations about work engagement measurement using the Utrecht Work Engagement Scale (UWES).

Work Engagement Concept and Its Study Among Rescue Workers
Considering the frequent occurrence of disasters due to natural or human causes, particularly attention must be given to rescue workers dealing with critical incidents or disasters, since chronic emotional demands and coping with stressful situations harms their health and the services provided by their institution (Queirós et al., 2013;Rabjohn, 2013;International Labour Organization, 2016). Therefore, these professionals can develop negative emotional states, psychological disorders or mental diseases (e.g., anxiety, loneliness, panic attacks, depression, posttraumatic stress disorder, vicarious trauma), but particularly, they can experience chronic job stress and become vulnerable to burnout (Cieslak et al., 2014;Katsavouni et al., 2016;Krok, 2016). Several studies have emphasized the relationships between burnout and work engagement (Schaufeli et al., 2002a;Maslach and Leiter, 2008;Leon et al., 2015;Leiter and Maslach, 2017;Schaufeli and De Witte, 2017), suggesting that promoting work engagement can prevent burnout (Argentero and Setti, 2011), thus ameliorating rescue workers' mental health and helping their stress management (Nishi et al., 2016;O'Neill and Rothbard, 2017). In fact, work engagement seems to be an important psychological resource for rescue workers, since it may protect them from the risk of developing negative health effects (Setti and Argentero, 2014), depression, sleep disturbances, relational conflicts, burnout, compassion fatigue or vicarious trauma (Lenthall et al., 2009;Kumar, 2011;Ray et al., 2013), or presenteeism (Burton et al., 2017).
For several decades, Psychology has focused on pathology and psychological illness, which has led to a restricted view of human behavior that has consequently contributed to the depreciation of resources, qualities and skills of human beings as promoters of individual growth and psychological well-being (Diener et al., 1999). However, the recent emergence of Positive Psychology has counteracted this focus, valorizing the study of the positive part of human experience (Luthans, 2002), aiming to increase the quality of life, to promote individual qualities, resources, and characteristics that increase psychological wellbeing and growth (Seligman and Csikszentmihalyi, 2000). This approach within Positive Psychology explores building positive qualities, instead of focusing on repairing psychological illness. One of the theoretical concepts that arises within this perspective is work engagement (Schaufeli et al., 2002a), focused on the positive aspects of behaviors, the well-being and the optimal functioning of each person (Schaufeli and Salanova, 2007).
Although work engagement is a relatively new construct, it has been viewed as being comprised of slightly different aspects, such as physical, cognitive, emotional and mental expression by employees within the exercise of their roles (Kahn, 1990). Distinctions between organizational and job engagement have appeared regarding employees' different roles suggesting that engaged workers invest and dedicate themselves deeply into their tasks, which are viewed as motivated situations (Saks, 2006). Other authors Leiter, 1997, 2017;Maslach et al., 2001) have considered work engagement as a concept diametrically opposed to burnout, being defined as an energetic state of involvement in vital personal activities, but quite different from job satisfaction or organizational commitment (Hallberg and Schaufeli, 2006).
While the discussion about burnout and work engagement as different constructs or as opposite poles of the same construct continues (Schaufeli et al., 2002b;Macey et al., 2009;Cole et al., 2012;Leon et al., 2015;Leiter and Maslach, 2017;Schaufeli and De Witte, 2017), we adopt Schaufeli and colleagues' perspective on work engagement as a positive motivational, cognitiveaffective, and emotional persistent state (Schaufeli and Salanova, 2007;Schaufeli, 2017b), divided into three factors or dimensions: vigor, dedication and absorption. Vigor is related to persistence, energy and mental resilience during work, while dedication implies being involved in the work and living with enthusiasm, inspiration, and pride, and finally, absorption means to be concentrated and immersed in one's work, with the perception that time flies (Schaufeli et al., 2002b). Thus, workers with high levels of work engagement see themselves as full of energy when facing their professional activities (Albrecht, 2010).
Work engagement is also related with other important occupational health and organizational outcomes, including in-role and extra-role behavior (Bakker et al., 2004). It is positively correlated with performance (Cooper and Quick, 2017;Reijseger et al., 2017); more specifically, a direct and positive relation has been found with job performance and a direct and negative relation with turnover intention (Karatepe and Avci, 2017;Kim, 2017), proactive behaviors and employee intrapreneurship (Gawke et al., 2017), job crafting and personal resources (Xanthopoulou et al., 2009;Vogt et al., 2016), general job resources , structural and social resources (Harju et al., 2016), and job, career and life satisfaction (Karatepe and Karadas, 2015).

Work Engagement Measurement with the Utrecht Work Engagement Scale
In line with its three-dimensional definition, the work engagement measure was divided into three sub-scales, in a psychometric instrument called Utrecht Work Engagement Scale (UWES), a self-report instrument that includes the vigor, dedication and absorption factors (with six, five and six items respectively). Although the UWES presents good psychometric properties. it is common to find very strong correlations between the factors (Chughtai and Buckley, 2013;Agarwal, 2014), leading some studies to prefer the single factor measure (Shimazu et al., 2008;Seppälä et al., 2009;Klassen et al., 2012), while others (Shimazu et al., 2008;Alok, 2013) tested UWES as a second-order construct (Edwards, 2001). Kulikowski (2017) conducted a literature review with 21 studies using Confirmatory Factor Analysis (CFA) of UWES; six of these recognized the three-factor version as superior to the one-factor solution, six studies found the opposite, eight studies found both structures as equivalent, and one study failed to confirm both structures. However, usually, the CFAs present good results with the threefactor structure. Moreover, various studies found invariance between different occupational and racial groups and countries (Storm and Rothmann, 2003;Xanthopoulou et al., 2012), and relative longitudinal invariance for the short version (Seppälä et al., 2009). This shorter version with nine items is divided into the same three dimensions, with three items per dimension .
The UWES has been used in several studies with rescue workers and translated into many languages (Schaufeli and Bakker, 2003). Several studies have shown that the structure of the UWES remains unchanged among different occupational groups. Table 1 presents a synthesis of 30 studies between 2003 and 2017, conducted mainly in Europe but also in America, Asia, Africa and Oceania, with professionals with occupation statuses similar to the ones of our study, revealing that work engagement levels are moderate to high among these rescue workers.
Considering previous studies with UWES among these professionals, our research question was to investigate if UWES is an adequate instrument to assess rescue workers' engagement. Thus, this study aims to estimate work engagement levels of rescue workers (namely comparing nurses, firefighters, and police officers) and to assess the validity evidence related with the internal structure of the Portuguese versions of the UWES-17 and UWES-9, namely to evaluate their dimensionality (comparing the fit of the three-factor first-order model of the UWES-17 and UWES-9, and of a possible UWES-9 second-order model), measurement invariance between the three occupational groups, and reliability of the scores, according to the Standards for Educational and Psychological Testing framework (American Educational Research Association et al., 2014). According to these aims, we have as hypotheses: (H1) UWES-17 presents validity evidence that supports its usage among rescue workers; (H2) UWES-9 presents validity evidence that supports its usage among rescue workers; (H3) UWES-9 presents measurement invariance that allows for comparing different occupational groups; (H4) UWES-9 reveals different work engagement levels between the three different occupational groups of rescue workers.

Participants
A sample of 3,887 Portuguese rescue workers completed the UWES-17 (except 428 workers that only completed the UWES-9), with data collected from several research projects (seven studies) developed at the Psychosocial Rehabilitation Laboratory from FPCEUP, between 2010 and 2016. Data were from different Portuguese districts, but with all occupations contributing to the data during the same year. The average age was 33.35 years old (SD = 9.84), 50% male, 44% married (and 50% single), 62% at least had a college degree. Regarding occupations, 50% were nurses, 39% voluntary firefighters (according to the Portuguese firefighters' organization, they can be voluntary or professional firefighters; this study collected data from voluntary firefighters, a larger group than professional firefighters) and 11% police officers.

Measures
A sociodemographic questionnaire was administered, with items such as age, sex, marital status, school level and occupation.
Despite other questions specific to each rescuer occupation, these questions were asked to all groups. To assess work engagement, we used the UWES ( Table 2) developed by Schaufeli et al. (2002b) in its Portuguese version (Marques-Pinto and Picado, 2011).

Procedures
Participants were invited to voluntarily participate after formal authorization from each organization where they were working at and filled out printed versions of the UWES anonymously. The study was approved by the Ethics Committee of University of Porto, Portugal, and written informed consent from all research participants was obtained.

Data Analysis
To study the dimensionality of the UWES versions, a CFA was conducted to verify if the three-factor structure proposed by the UWES' authors presented an adequate fit for this study's sample. We used as the goodness-of-fit indices the TLI (Tucker Lewis Index), NFI (Normed Fit Index), χ 2/ df (ratio chi-square and degrees of freedom), CFI (comparative fit index) and the RMSEA (root mean square error of approximation). The fit of the model was considered good for χ 2/ df smaller than 5, CFI, NFI and TLI values above 0.95 and RMSEA values below 0.08 (Hoyle, 1995;Boomsma, 2000;McDonald and Ho, 2002;Marôco, 2014;Byrne, 2016).
To analyze convergent validity evidence, the average variance extracted (AVE) was estimated as described in Fornell and Larcker (1981). Values of AVE ≥ 0.5 were considered indicative of the constructs' convergent validity evidence, for both versions of the UWES per the proposal by Hair et al. (2009).
Discriminant validity indicates that the items that represent a dimension are not strongly correlated with other dimensions (Fornell and Larcker, 1981;Marôco, 2014). For two factors x and y, if AVEx and AVEy ≥ ρ 2 xy (squared correlation between the factors x and y) there is evidence of discriminant validity (Marôco, 2014).
The reliability of the UWES scores was investigated using internal consistency estimates, Composite Reliability (CR), the ordinal Cronbach's alpha coefficient (α) and the ordinal omega coefficient (ω) for each factor and for the total scale. The second-order factor reliability was also calculated with the omega coefficient (semTools Contributors, 2016).
The measurement invariance corresponds to the situation where links between the items and the construct studied are invariant for different contexts, such as countries, religions, times, or languages (Gregorich, 2006;Little, 2013;Brown, 2015;Hayduk, 2016). This should be tested on the construction or adaptation plan of any psychometric instrument (Edwards et al., 2017). This analysis will show if in different measurement conditions we are measuring the same attributes. Without evidence for this, the basis for drawing inferences will be severely compromised, and, as a result, the findings of possible differences between samples of different groups cannot be unambiguously interpreted (Horn and McArdle, 1992).
Establishing measurement invariance consists of the estimation of an increasingly constrained set of structural    , 2012). To detect whether the same three-factor model holds in each group of professionals, a group of nested models with indications of equivalence is needed (Marôco, 2014). For the three-factor first-order version of the UWES-9, the first model to be tested is one where only the factor loadings are equal across groups, but the thresholds are allowed to differ between groups; this is called metric invariance and assesses if the professionals in the different occupational groups attribute the same meaning to the three work engagement factors. Next, to assess scalar invariance, a model in which the thresholds and loadings are forced to be the same is tested; it indicates that the three work engagement latent variables (in other words, the loadings), and that the levels of the thresholds are equal for the three different kinds of professionals (van de Schoot et al., 2012). As a result, the occupational groups can be compared on their scores for the vigor, absorption and dedications variables. The next step tests a model in which the residual variances are also forced to be equal across the three groups of professionals, achieving the so-called full uniqueness Measurement Invariance; if obtained, the explained variance for all items is the same for the occupational groups. This means the three factors are measured identically across professional groups. Finally, a model in which the latent means are the same for the different groups is tested. This approach is based on the recommendations of Millsap and Yun-Tein (2004), who specifically discuss the measurement invariance for ordered-categorical measures using the theta-parameterization. Most researchers still use multi-group confirmatory maximum likelihood factor analysis of a Person covariance matrix irrespective of the nature of the measures (Koh and Zumbo, 2008). The approach used in this study tries to be more precise and adequate given the kind of measures adopted (Hirschfeld and von Brachel, 2014). For the second-order model, a different approach (with more models) is needed. A group of seven models were tested, based on the recommendations of Chen et al. (2005) to test the invariance in second-order models, and in the recommendations of Millsap and Yun-Tein (2004) due to the ordered-categorical nature of our measures. The following models with increasing constraints were compared: (a) configural invariance, (b) firstorder factor loadings, (c) second-order factor loadings, (d) thresholds of measured variables, (e) intercepts of first-order factors, (f) disturbances of first-order factors, and (g) residual variances of observed variables. Invariance is supported if the CFI < 0.01 between constrained and free models is presented (Cheung and Rensvold, 2002) and if χ 2 test comparing the fit of the constrained vs. free models is not statistically significant (Satorra and Bentler, 2001). Frontiers in Psychology | www.frontiersin.org All statistical analyses were performed with R (R Core Team, 2017) and Rstudio (RStudio Team, 2017). The descriptive statistics were obtained with the skimr package (Arino de la Rubia et al., 2017). The CFA analysis was conducted with the lavaan package (Rosseel, 2012) using the Weighted Least Squares Means and Variances (WLSMV) estimation method. The reliability estimates and measurement invariance (for the first-order model) were calculated with the semTools package (semTools Contributors, 2016), the measurement invariance for the second-order model was analyzed with lavaan (Rosseel, 2012), the Mardia Kurtosis (Mardia, 1970) was assessed with the psych package (Revelle, 2017), and the doBy package (Højsgaard and Halekoh, 2016) for the percentiles. The comparisons of the raw levels of the factors (estimated as the mean of the items for each factor) between the occupational groups were tested using one-way ANOVA with pairwise comparisons using the stats package (R Core Team, 2017). The effect sizes were calculated using Partial Eta Squared from the lsr package (Navarro, 2015).

RESULTS
Results are presented following the Standards for Educational and Psychological Testing framework. Table 3 (UWES-17) and Table 4 (UWES-9), there were no items with skewness and kurtosis values that were indicative of severe normality deviations (Kline, 2016), the Mardia's Multivariate Kurtosis for the UWES-17 items was 279.4, p < 0.001, and for the UWES-9 version it was 189.4, p < 0.001. No outliers were deleted. All possible answer values for each item were also present. The mean for most items was close to 4. Therefore, factor analysis with WLSMV, to account for the ordinal level of measurement of the items, can be done without concerns about the validity of its estimates.

Discriminant validity evidence
The discriminant validity was assessed by comparing the AVE of the factors with the squared correlation of the factors, as proposed by Fornell and Larcker (1981) and described in Marôco (2014). Discriminant validity evidence is obtained when the AVE for two scales is larger than the squared Pearson correlation between the two scales. For the UWES-17, AVE vigor = 0.65 and AVE dedication = 0.72 were smaller than r 2 VD = 0.86, AVE absorption = 0.52 and AVE dedication = 0.72 were smaller than r 2 AD = 0.85 and AVE vigor = 0.65 and AVE absorption = 0.52 were smaller than r 2 VA = 0.86. For the UWES-9, AVE vigor = 0.80 and AVE dedication = 0.78 were not larger than r 2 VD = 0.83, AVE absorption = 0.61 and AVE dedication = 0.78 were smaller than r 2 AD = 0.82, and AVE vigor = 0.80 and the AVE absorption = 0.61 are, respectively, larger and smaller than r 2 VA = 0.76.

Second-order construct
When there is a higher order latent variable (work engagement) that is hypothesized to consider the relations among the lower order factors (vigor, dedication, and absorption in this case), a second-order latent model is admissible (Chen et al., 2005). Another possibility to consider second-order latent models is when the first-order dimensions are substantially correlated with each other (Marôco, 2014), which is verified in the first-order three-factor UWES-9 model. For the UWES-9 model, a secondorder latent factor was proposed. A correlation between the residual variances of item 1 and item 4 of the Vigor dimension was added. Regarding the UWES-9 with a second-order latent factor, the fit was very good (Figure 2; χ 2 (24) = 498.849, p < 0.001, n = 3,623; CFI = 0.998; TLI = 0.997; NFI = 0.998; RMSEA = 0.074; P(rmsea ≤ 0.05) < 0.001; IC90 ]0.068; 0.080[).

Reliability of the Scores: Internal Consistency
To assess the internal consistency, first we used the Composite Reliability of the factors. Also standard Cronbach's alpha coefficient (α) and omega coefficient (ω) were used for each dimension proposed in the UWES-17 and in the UWES-9, and to assess the consistency of the total scales ( Table 5).
Reliability of the higher-order construct was very good. The proportion of the second-order factor explaining the total score (ω L1 ) was 0.93, the proportion of the second-order factor explaining the variance of the first-order factor level (ω L2 ) was 0.98, and the proportion of observed variance explained by the second-order factor after controlling for the uniqueness from the first-order factor (ω partial L1 ) was 0.94. The structural weights for the second-order factor model were: vigor (β = 0.99; p < 0.001); dedication (β = 0.98; p < 0.001); and absorption (β = 0.93; p < 0.001).

Measurement Invariance
To detect whether the same original three-factor (first-order model) holds in each occupation (Table 6), a group of nested models with indications of equivalence is needed (Marôco, 2014). The hypothesized scale invariance was supported by the Cheung and Rensvold (2002) CFI criterion, but not by the χ 2 criterion (Satorra and Bentler, 2001), although the second criterion is more restrictive. According to the Cheung and Rensvold (2002) criterion, the results supported the full uniqueness measurement invariance between the occupational groups (nurses, firefighters, and police officers). However, there were significant differences between the means of the UWES dimensions for the different occupational groups. It could be noticed that these comparisons are valid since the measures showed scalar invariance. The test of the measurement invariance for the second-order model for categorical variables is shown in Table 7. According to CFI criteria, only partial (metric) second-order invariance is supported in the second-order UWES-9 model for the three different occupational groups.

Dimension Comparisons
Considering rescue workers' occupations, Table 8 presents percentiles and mean levels of each work engagement dimension, and a comparison between nurses, firefighters, and police officers. For the hypothesized differences between occupational groups, some statistical significant differences in work engagement were found between the three groups (Table 9), being highest among firefighters and lowest among nurses. Percentile analyses reveal the exact value for each category and occupation, reflecting the tendency revealed by the mean analysis of each sample.

Hypotheses Findings
Considering the psychometric properties of data gathered with the UWES in the present sample, our hypotheses were confirmed since UWES-17 and UWES-9 presented validity evidence that supports its usage among rescue workers (respectively H1 and H2), and UWES-9 presented measurement invariance that allows for comparing occupational groups (H3). In fact, the data presented better psychometric properties for the UWES-9 and supported a second-order latent factor for the threefactor structure due to strong correlations between the first-order dimensions. This solution is supported by other international studies, presenting better goodness-of-fit indexes for the shorter version of the instrument (Nerstad et al., 2010;Hernandez-Vargas et al., 2016). The convergent validity evidence and reliability of the scores were good, although the discriminant validity was not acceptable, leading to the proposal of a secondorder latent factor. There have been several reports that discuss the dimensionality of work engagement as assessed by UWES (see Table 1). Our data suggests that work engagement can be defined as a second-order factor that is reflected in the workers' vigor, dedication, and absorption. Thus, we suggest that work engagement can be estimated as a single value in research contexts and in clinical settings facilitating the use of the UWES. A second-order latent factor has advantages over a single firstorder factor for the UWES as proposed by several researchers (see Table 1) since it allows for the estimation of a single value for work engagement while retaining the three first-order factors proposed (Schaufeli and Bakker, 2003). Regarding H4 (UWES-9 can evaluate different occupational groups of rescue workers), and considering rescue workers' occupation specificities, work engagement levels varied; despite being moderate to high, they were higher among firefighters and lower among nurses. This can be explained by the fact that Portuguese firefighters are traditionally volunteers (paid or not), being highly engaged in their tasks and frequently assuming physical risks to help in rescue situations (e.g., the tragedy of the Portuguese forest fire during June 2017 at Pedrógão Grande where a firefighter died, and several were severely burned; or the Portuguese forest fire during the summer of 2013 at Caramulo where eight firefighters died). According to Portuguese data on firefighters, in 2014 there were 42,000 firefighters, only 6,000 of whom were professionals and the remaining were volunteers, despite all of them being trained by the Portuguese Firefighters' National School. Moreover, some studies explain that aid workers, such as firefighters, can present a sense of coherence, a concept developed by Antonovsky (1987). Related to the salutogenic theory, it is defined as the "ability of the persons to understand what happens around them, to what extent they were able to manage the situation on their own or through significant others in their social network, and the ability to find meaning in the situation" (Eriksson and Lindström, 2005, p. 460). This can help someone to better cope with traumatic events life adversities (Veronese and Pepe, 2014;Veronese et al., 2016). The study of 464 firefighters by Dudek and Koniarek (2000) revealed that post-traumatic stress symptoms were associated with a low level of sense of coherence, especially if they already have post-traumatic stress disorder. Additionally, firefighters, like police officers, share some personality traits like sensation seeking (Perrott and Blenkarn, 2015), that increase their motivation to pursue certain occupations. Nurses presented the lowest work engagement values, maybe due to the external demands their occupation is facing, namely lack of human or material resources, conflicts with patients or being aggressed by them or their families, etc. (Lin et al., 2016;Marôco et al., 2016). Indeed, recent austerity measures adopted in countries such as Portugal have had significant negative effects on the availability of healthcare system resources, including reduced staffing and employment conditions (Dussault and Buchan, 2014). Regarding the work engagement dimensions, it is important to notice that they are differentiated according to each rescue occupation, since vigor is related to physical energy, used more by firefighters and police officers, while dedication is clearly higher among firefighters who voluntarily seek this activity. Finally, absorption must be interpreted according to each task, which was higher among firefighters and police officers, especially when they are in rescue operations where time "flies, " while nurses have more rigid protocols to obey, accomplishing schedules and some rigid tasks. Results are similar to other studies where rescue workers present moderate to high work engagement levels (see Table 1), with nurses having lower work engagement than police officers  , and a lack of studies about firefighter's work engagement (curiously the one found was with a Portuguese sample).

Limitations and Recommendations for Future Research
This study has some limitations, namely having a convenience sample collected using participants' voluntary participation. In the future, respecting voluntary participation, we would like to attract more participants, trying to have equal representation among the occupations. In fact, in this sample the percentages were not homogeneous, with nurses representing half of the sample, and police officers comprising only 11% of the sample. Firefighters were also included, but not other rescuers such as medical emergency staff, civil protection technicians, social workers or psychologists, who, for example, during the summer of 2017 in Portugal, participated in rescue tasks for the population due to the tragic forest fires. Furthermore, we did not have any source of evidence-based on relations to other variables, which can provide external validity evidence of UWES when compared with other concurrent or divergent measures (American Educational Research Association et al., 2014). This was a consequence of collecting data from each occupational group separately or during a different year, according to different research projects developed in the same University research unit, which were addressing different topics, having in common the measurement of work engagement with UWES. In the future, we would like to increase the sample, with replication efforts to collect data among more diverse rescue occupations. Also, we intend to have more homogeneity between rescue workers' samples sizes, having thus a stronger national representativeness, which can be obtained using web-based versions (Ilieva et al., 2002) of UWES (with other measures) through, for example, LimeSurvey (Limesurvey GmbH, 2017). The advantages of webbased surveys can also include increased delivery speed, control and monitoring, automatic data validation, and decreased data entry error (Falletta and Combs, 2002). Finally, we would like to combine work engagement measurement with other variables such as burnout, trauma or job satisfaction, as frequently studied in numerous research (e.g., Argentero and Setti, 2011;Ângelo and Chambel, 2015;Espinoza-Parra et al., 2015).

Theoretical Implications
The study has some theoretical implications, enabling a better understanding of work engagement dimensions related with each occupational group specificities (e.g., vigor is higher among firefighters and police officers who use more physical energy). Additionally, the recent Schaufeli (2017a) publication of the general engagement scale (UGES), with 9 items evaluating daily activities and not only work/job activities will broad the evaluation of engagement and its relationship with psychological variables such as positive and negative affect, satisfaction with life, personality traits (Schaufeli, 2017b), or sense of coherence and trauma (Dudek and Koniarek, 2000;Veronese and Pepe, 2017). Thus, work engagement seems to be an important concept to study among workers and organizations, considering that organizations are facing more competition, innovation, and constant changes, and also a decrease of resources (Graffigna, 2017). Engaged workers are thus an important element for organizations' productivity or to the care/services they provide to citizens.

Practical Implications
Regarding practical implications, this validation of a Portuguese version of UWES can contribute to the improvement of rescue workers' quality of life, since these results can alert individuals and institutions about the need to explore work engagement level to buffer against negative physiological health effects (Lenthall et al., 2009;Kumar, 2011;Ray et al., 2013;Setti and Argentero, 2014;Burton et al., 2017). Moreover, reference values for the Portuguese firefighters, nurses and police officers' populations make it possible to establish comparisons with other occupational groups of rescue workers (e.g., physicians, social workers). Taken together, the findings of the present study have important ecological implications for Portuguese policymakers and human resource managers who can structure rescue workers' work conditions, or develop intervention programs to promote their work engagement, and therefore help prevent burnout and other mental health conditions that affect these professionals' wellbeing.

CONCLUSIONS
Regarding the four research hypotheses, all were confirmed, since convergent validity evidence was found. Due to the lack of discriminant validity evidence, a second-order latent factor was proposed. Additionally, reliability evidence presented good values, and we tested measurement invariance for ordinal variables, using an innovative approach for categorical variables.
Moreover, UWES can discriminate work engagement levels for each occupational group of rescue workers. Thus, this UWES version can be used to collect more data for future comparative studies. Additionally, due to natural or human causes (e.g., forest fires, multi-victim accidents, terrorist attacks), rescue workers play an increasingly important role in modern society, and the work engagement concept will help to better understand rescue workers' motivation to continue performing their tasks even when they risk their own life and frequently face tragedies.

ETHICS STATEMENT
This study was carried out in accordance with the recommendations of guidelines of the Faculty of Psychology and Education Sciences of the University of Porto with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Ethics Committee of the Faculty of Psychology and Education Sciences of the University of Porto.