1. Introduction
Nowadays, there is no doubt that global climate change is a reality, and different factors affect its occurrence, such as greenhouse gases mainly due to emissions from human activities, melting and loss of glaciers, ocean warming, sea level rise, and other extreme weather events (droughts, floods and intense tropical cyclones, heat waves or forest fires) [
1].
The World Health Organization (WHO) [
2] and The Lancet [
3] have stated that climate change will adversely affect some of the fundamental determinants of health: clean air, clean water, sufficient food, and safe housing.
Although the health effects of climate change occur in all people, children are the population group most vulnerable to exposure to environmental risk factors due to their immaturity, psychosocial dependence, and also the lack of positive stimuli and communication in the first three years of life that will negatively determine the development of physical, mental and social habits [
4,
5,
6]. Over two-thirds of children's health is threatened by environmental risks in their homes, where they learn and play [
7]. Children under the age of five are particularly vulnerable, as evidenced by the fact that this age group accounts for more than 88% of the diseases caused by climate change [
8]. Future generations of children will experience a disproportionate burden of the adverse effects of climate change compared to previous generations. The impact of climate change on the health of this population is wide-ranging, from physical aspects such as increased vulnerability to disease (asthma, allergies, or vector-borne diseases) to other effects that damage mental health (depression or climate anxiety) [
9]. The WHO estimates that approximately one-third of the disease burden in developing countries is attributable to modifiable environmental factors, such as indoor and outdoor air pollution, unsafe water, inadequate sanitation and hygiene [
10].
The 2030 Sustainable Development Goals (SDGs) clearly state the need for Education for Sustainable Development (ESD) to sustain planetary health ([
11]). Therefore, it is necessary to include environmental issues in undergraduate nursing education and tools to assess their effectiveness in increasing environmental competencies, attitudes, knowledge, and skills in children's environmental health, as future nurses will be on the frontline of addressing the health consequences of climate change. This is why it is important to determine environmental competencies in nursing students, and previous research carried out in the Spanish context identified through validated instruments that the attitudes, knowledge, and skills related to children's environmental health in undergraduate nursing students were not adequate [
12].
There are different instruments to determine attitudes, knowledge, and skills on children's environmental health: the Sustainability Attitudes in Nursing Survey 2 (SANS_2), the Children's Environmental Health Knowledge Questionnaire (ChEHK-Q) and the Children's Environmental Health Skills Questionnaire (ChEHS-Q), which were developed and used in the Spanish context ([
13,
14]) and which have been adapted in some countries such as Korea [
15], China [
16] and United Kingdom [
17]; however, there is currently no measurement instrument validated in Portuguese. Therefore, the main objective of this study is to validate in the Portuguese context three scales to assess attitudes, knowledge and skills on children's environmental health. The second objective is to determine the knowledge, attitudes, and skills of the Portuguese nursing students in the sample studied.
2. Materials and Methods
A cross-sectional observational study was developed to translate, adapt, and validate the questionnaire consisting of the following three scales:
Attitude Scale (SANS_2). This scale was designed to assess nursing students' attitudes toward sustainability and climate change. It consists of five items, with response options ranging from one to seven on a Likert scale, with a maximum score of 35 points. The score range calculated by the study by Álvarez-García et al. [
12] is: Excellent attitude (> 90% of perceived attitudes), very good attitude (70-89%), good attitude (50-69%), insufficient attitude (30-49%) and poor attitude (< 29%). Reliability analysis of this tool previously revealed a Cronbach's alpha of 0.82, and the five items loaded on a single factor explained 58% of the variance [
13].
Knowledge scale (ChEHK-Q). This scale measures nursing students' knowledge of children's environmental health. It consists of 26 items, with true, false and I don't know response options. The maximum score on the scale is 26 points. The range of scores is: Excellent knowledge (> 90% correct answers), very good knowledge (89-80% correct answers), good knowledge (79-60% correct answers), insufficient knowledge (59-40% correct answers), and poor knowledge (< 39% correct answers). This scale shows a good fit and a reliability of 0.98 for items and 0.70 for people based on Rash's model in the Spanish population [
14].
Skills scale (ChEHS-Q). This scale measures the skills of nursing students in children's environmental health. It consists of 12 items with response options ranging from one to five on a Likert scale. The maximum scale score is 60 points. The range of scores is: Excellent skills (> 90% of perceived skills), very good skills (89-80% of perceived skills), good skills (80-70% of perceived skills), insufficient skills (79-50% of perceived skills), and poor skills (< 49% of perceived skills). This scale shows a good fit and a reliability of 0.87 for items and 0.76 for people according to Andrich's model in the Spanish population [
14].
In this study, the guidelines of the International Test Commission [
18] were followed for the translation, adaptation and validation of the questionnaire. The questionnaires are available in
Appendix A. This was carried out in two phases:
2.1. Phase 1: Translation and Adaptation Process
The questionnaire was translated into Portuguese using a centered or asymmetric method [
19], as the questionnaire was intended to be translated without affecting the wording of the original scales. It was decided to do this because the scales were used in the first part of the study conducted in Spain [
14] with a similar population of nursing students and have distinctive and desirable characteristics that should remain consistent.
To achieve semantic equivalence, a translation was carried out by two Portuguese native speakers in the research team. This was followed by a back-translation by two other native Spanish speakers. A matrix was constructed to visualize the original item, the translated items, and the back-translated items. The members of the research team then compared the original and back-translated version, and the necessary adjustments were made (
Appendix B). In order to achieve conceptual equivalence, the translators were members of the research team who are experts in the study construct and faculty members of the institution where the questionnaire will be used, so they are familiar with the target population of the questionnaire.
The translated version of the questionnaire was pilot-tested with 10 nursing students to assess their opinions on the readability and comprehension of the items [
19].
2.2. Phase 2: Validation Process
The population used for the psychometric assessment was undergraduate nursing students from a university in northern Portugal. The sample was reached by convenience, obtaining a sample of 326 students following the methodological recommendations of Polit and Beck [
19]. The questionnaire was sent by email, the institutional platform, or other digital applications during the academic year 2023-2024.
The psychometric analysis of the scales was intended to be conducted using one-parameter models based on item response theory, which is based on three assumptions [
20]. First, the unidimensionality of the scales was tested as suggested by Jiménez-Alfaro and Montero-Rojas [
21]. This was done, through an exploratory factor analysis of each scale and then Horn's parallel analysis to compare the variance values explained with those obtained in the simulation [
22]. Furthermore, at least 50 % of the variance needed to be explained, and each extracted factor needed to explain at least 5 % of the variance to consider the multidimensionality of the relevant constructs [
23]. There also had to be theoretical concordance with the previous known dimensions and the analysis results. Secondly, local independence was also tested using Yen's Q3 statistic. The recommended values are set at ± 0.3 [
24]. Thirdly, the Rasch Model analyzed subject (θ) and item (β) parameters on the attitude, knowledge, and skill variables, which will be shown in an item map. The infit and outfit statistics were used to fit the data to the models. Optimal fit values range from 0.7 to 1.3 [
25].
Finally, because the attitude scale (SANS_2) did not meet the above assumptions, it was decided to validate according to the classical measurement theory (since it was initially developed under the assumptions of this theory), testing construct validity with an exploratory factor analysis followed by a confirmatory factor analysis and reliability measures such as Cronbach's α and inter-item correlation with Pearson's coefficient.
The knowledge scale (ChEHK-Q) was analyzed using the Rasch model. The difficulty of the items was calculated, with 0 being the average level of difficulty. The separation index was also calculated, ideally having values greater than 2 [
26].
The ability scale (ChEHS-Q) was analyzed using Andrich's rating scale model, which is based on the same assumptions and parameters as the Rasch model. In addition to the (mis)fit of the items, this model analyzed the ordering of the thresholds by inspecting the locations of the estimated thresholds along the latent trait, as disordered thresholds may indicate inadequate functioning of the response categories.
For reliability analysis, an internal consistency value of 0.70 or higher was considered acceptable [
27].
These analyses were carried out with the jMetrik and JASP software.
2.3. Descriptive Analysis
Sociodemographic data of the participants were described. Descriptive values were calculated for the questionnaire scores, the number of hits, the ignorance index for the knowledge items, and the means for the attitude and skills items. The number of participants with excellent, very good, good, insufficient, or poor attitudes, knowledge, and skills was determined following the guidelines of the instruments. Finally, after checking the normality of the data, mean difference analyses (ANOVA with Bonferroni post-hoc tests) were carried out to compare the results obtained according to the sociodemographic variables collected. The data was analyzed in JASP.
3. Results
3.1. Cognitive Piloting
Cognitive piloting involved 10 students, 8 females (80%) and 2 males (20%); mainly from third year (80%, n = 8) or second year (20 %, n = 2). Most of them had not attended sessions on sustainability and nursing (60%, n = 6), although 3 students did (30%) and 1 student did but more than 3 months ago (10%).
Regarding the attitude items, the responses were distributed in the established range of 1 to 7 points, with the mean at a high level, reaching a value of 5.87 ± 0.74.
The knowledge items had a mean score of 12.22 ± 3.15 in a range of scores between 9-19. The most correct items were 1 and 24 referring to children's special vulnerability to environmental risks, which were correct for all participants. While the least correct items were 8 and 19 referring to smoking and pesticides, respectively, with a correct score of 22.23%. The rest of the items were distributed in a uniform range of correct answers.
Finally, the skills items obtained a mean score of 3.02 ± 0.34 evenly distributed on the scale of 1 to 5 proposed for these items.
Thus, although there were two knowledge items correct for all participants, no irregular response patterns were detected, as there should always be items in a range from very easy to very difficult to adequately classify respondents; as well as items in the different intervals of the Likert scales in the case of the attitude and skills items.
Post-questionnaire questions showed that the instructions made it clear how to complete the questionnaire, that the meaning of the items was understandable in all cases and that the format made the questionnaire easy to complete. The completion time ranged from 10-12 minutes.
3.2. Psychometric Study
The total number of students who completed the questionnaire was 326. The mean age was 21.87 ± 6.19 years, ranging between 17 and 56 years. The majority of participants were females 83.44% (n = 272) compared to 16.26% (n = 53) males and 0.31 % transgender (n = 1). The highest grade in which the students were enrolled was mostly third-grade with 108 participants (33.13%), followed by second grade with 81 participants (24.85%), 69 first-grade participants (21.17%), and 68 in fourth grade (20.86%). The majority had attended a session on sustainability and nursing, but more than three months ago, 67.49% (n = 220), 17.79% (n = 58) had not attended any session, and 14.72% (n = 48) had attended.
3.2.1. Attitude Scale (SANS_2)
First, an exploratory factor analysis was carried out to check construct validity. Factor analysis was found to be possible, as the Kaiser-Meyer-Olkin score was 0.794, and Bartlett's test of sphericity was statistically significant (
p < 0.001). One factor was found to explain 56.10% of the variance, confirming the original unidimensional structure of the scale (
Table 1).
Secondly, a confirmatory factor analysis was performed. The factor loadings ranged from 0.29 to 0.93 (
Figure 1). The model fit indices were determined as follows: χ
2 = 33.66, df = 4,
p < 0.001, RMSEA = 0.137, SRMR = 0.04, GFI = 0.99, CFI = 0.99.
The reliability measured by Cronbach's α coefficient was 0.84. The inter-item correlation values were lower than 0.20 only for item pairs 2 and 5 (0.19), 4 and 5 (0.17;
Appendix C).
3.2.2. Rasch Model of the Knowledge Scale (ChEHK-Q)
First, the unidimensionality of the knowledge scale was demonstrated. For this purpose, it was found that factor analysis was possible for the knowledge scale, as the Kaiser-Meyer-Olkin score was 0.762, and Bartlett's test of sphericity was statistically significant (
p < 0.001). In none of the analyses did we obtain good values of variance explained, nor did the items loading on the different components agree with our theoretical assumptions (
Table 2), therefore, the scale should be considered unidimensional as we did not find an internal structure that would justify the contrary.
To test the local independence of the items, Yen's Q3 correlations were calculated, and the values of the statistics did not exceed the established cut-off value (
Appendix D).
Regarding model fit, all 26 items showed central values of Infit and Outfit at 1, confirming the unidimensionality of the construct. The mean Infit was 0.98 (0.85 to 1.23), standardized Infit -0.41 (-3.01 to 4.13), Outfit 0.97 (0.71 to 1.23) and standardized Outfit -0.37 (-2.95 to 2.79). The difficulty values ranged from -2.08 for item 1 (The pediatric population is more susceptible to environmental threats due to their biological immaturity), the easiest item, to 2.78 for item 5 (Nitrogen oxide from fossil fuels in the home and tobacco smoke causes redness and burns on the skin), the most difficult. The item separation index was 7.05, and the person separation index was 1.75. Item separation was robust, but person separation showed some restrictions. The reliability for the item set was 0.98 and for the person 0.75. Again, these values reflect strong item quality but some restrictions in the range of persons in the sample, although both reliability values exceed the established limits. The individual item fit values can be examined in
Table 3. To appreciate the performance of the items, they are plotted in
Figure 2. In this figure, it can be seen that most of the items are concentrated at intermediate levels of ability, with some items being easier and some more difficult to measure adequately for the higher and lower-ability students in the sample.
3.2.3. Rating Scaling Model of the Ability Scale (ChEHS-Q)
The ability scale also proved to be amenable to factor analysis, as the Kaiser-Meyer-Olkin score was 0.885, and Bartlett's test of sphericity was statistically significant (
p < 0.001). This was followed by an exploratory factor analysis, using principal components as the extraction method. Using Horn's parallel analysis comparing the explained variance values with those obtained in the simulation, two factors were found (
Table 4). However, no justification was found for the two factors extracted after the oblique promax rotation, as they only reflected the items with positive and the items with negative sense (
Table 5). Therefore, the skills scale should also be considered as unidimensional.
To test the items' local independence, Yen correlations were calculated. The values of the Q3 statistic show deviations between items 3 and 6, 2 and 7, 6 and 8, 7 and 11, and 7 and 12. The highest correlation value is shown between items 3 and 6, -0.42. These correlations are detailed in
Appendix E.
In the skills scale model, all 12 items showed central values at 1 for Infit and Outfit, confirming the unidimensionality of the construct. The mean Infit was 1.00 (0.76 to 1.47), standardized Infit -0.12 (-2.76 to 4.39), Outfit 0.99 (0.73 to 1.44) and standardized Outfit -0.17 (-3.06 to 3.96). The difficulty values ranged from -0.41 for item 2 (I am NOT able to identify the environmental risks that can cause respiratory diseases in a child), the easiest item, to 0.35 for item 3 (I am able to identify the environmental risks that can cause neoplastic diseases in a child), the most difficult. The item separation index was 2.15, and the person separation index was 2.70. Item and sample separation were robust. The reliability for the set of items was 0.82, and for the person 0.88. In this case the reliability values are equally adequate in both cases.
Examining the individual item values (
Table 6), item 2 slightly exceeded the maximum values set for Infit and Outfit, without reaching the 1.50 value that would degrade the measure. In order to appreciate the performance of the items, they were represented in
Figure 3, which shows that the items are concentrated at intermediate and low levels of ability. The high ability of the students in the sample can be seen.
The score limits on the Likert scale range from -3.05 to 4.21, increasing in value progressively along the scale, with no setbacks, as shown in
Table 7. In addition, the values of Infit and Outfit are in line with the established values, with only the outfit of item 2 showing an under-adjustment.
3.3. Descriptive Analysis of the Questionnaire
Table 8,
Table 9 and
Table 10 below show the responses to the items on attitudes, knowledge and skills on children's environmental health.
The final mean score of the attitude scale was 28.15 (5-35) ± 4.61. 94 students (30.62%) had an excellent attitude, 152 students had a very good attitude (49.51%), 55 students had a good attitude (17.92%), 5 students had an insufficient attitude (1.63%), and 1 student had a poor attitude (0.33%).
On the knowledge scale, the mean final score was 14.92 (0-26) ± 4.51. There were 6 students with excellent knowledge (2.18%), 33 students with very good knowledge (12%), 124 students with good knowledge (45.09%), 84 students with insufficient knowledge (30.55%), and 28 students with poor knowledge (10.18%). The ignorance index was higher than 50% only for item 8 (Passive smoking is associated with the development of acute leukemias in children).
The final mean score on the skills scale was 42.51 (24-60) ± 6.41. Twelve students (4.53%) obtained excellent skills scores, 49 students had very good skills (18.49%), 95 students had good skills (35.85%), 101 students had insufficient skills (38.11%), and 8 students had poor skills (3.02%).
3.4. Bivariate analysis of the questionnaire
In a gendered analysis, male students obtained better environmental health skills (16.41 ± 0.88) than female students (14.93 ± 0.47); F = 4.601, p = 0.033, η2= 0.04.
Regarding the highest grade in which the student was enrolled, the third-grade student had better attitudes (29.50 ± 0.85, p = 0.004) and knowledge (16.10 ± 0.77, p = 0.001) compared to the first-grade student with lower attitudes (26.81 ± 1.17) and knowledge (13.14 ± 1.05); F = 4.420, p = 0.005, η2= 0.04. and F = 4.889, p = 0.003, η2= 0.05, respectively.
Taking into account the sessions on sustainability and nursing attended by the students, those who had attended more than three months ago had better knowledge (16.51 ± 0.66, p = 0.046) but lower skills (41.22 ± 1.02, p = 0.018) than the students who had attended recent sessions, with lower knowledge (15.44 ± 0.82) but higher skills (44.32 ± 1.26); F = 4.843, p = 0.009, η2= 0.03 and F = 4.843, p = 0.009, η2= 0.018.
4. Discussion
The main objective of the research was to translate, adapt, and validate in the Portuguese context three scales to assess knowledge, skills, and attitudes on children's environmental health in nursing students. The Portuguese versions of ChEHK-Q, ChEHS-Q and SANS_2 present adequate psychometric properties, providing tools that effectively measure knowledge, skills, and attitude on children's environmental health in Portuguese undergraduate nursing students.
The translation, adaptation and validation of instruments in different languages facilitates safe comparisons if the resulting scales are valid and reliable as well as semantically equivalent [
28]. Prior to data collection, a pilot test was conducted on a group of participants that was representative of the target population in order to identify problems related to item clarity [
29]. In this sense, a rigorous process has been followed to obtain instruments that allow for adequate and comparable results over time and across organizations.
In the psychometric analysis of the scales, two types of psychometric theories were used: classical measurement theory and item response theory (Rasch model and Andrich's model). The knowledge and application of these theories according to the fulfillment of each of the criteria they assume ensures the robustness of the instruments obtained [
30].
The use of item response theory techniques converts the raw data to a linear scale that facilitates analysis. The representation of the latent trait of the respondents with the difficulty of the items in the item maps shows the subject-to-subject and item-to-item relationship facilitating the understanding of the results. The model fit parameters for the ChEHK-Q and ChEHS-Q questionnaires in the studied population are considered optimal, except for item 2 of the skills questionnaire, which slightly exceeds the threshold value. Although the analysis of the mapping of the items of the ChEHS-Q questionnaire shows that they are concentrated at intermediate and low levels of ability, the separation index of both items and individuals indicates that the items are adequately distributed in terms of the range of difficulty and the questionnaire is able to detect differences between students with different abilities. Even so, the separation values found are below those found in studies conducted in China [
16]. In the case of the knowledge questionnaire, the separation index of the students is very similar to the values found in research carried out in Spanish [
14] and United Kingdom [
17] students, with values of 1.53 and 1.75, respectively, below the index considered adequate. The item and person reliability values obtained in the knowledge and skills questionnaire are considered acceptable when reaching the set values (equal to or higher than 0.70) [
27], as found in previous research [
12,
14], with higher reliability in the sample of Chinese students [
16].
The values for knowledge (14.92 ± 4.51) and skills (42.51 ± 6.41) of Portuguese students are similar to those found in the study carried out in Spain (knowledge: 15.19 ± 4.08 and skills: 41.10 ± 7.61) [
14] and higher than those found in Chinese nursing professionals (knowledge: 13.04 ± 6.27 and skills: 35.28 ± 8.16) [
16] and English students (knowledge: 11.59 ± 5.09 and skills: 38.36 ± 7.86) [
12]. According to the established cut-off points, and grouping the values into adequate (including the values of excellent, very good, or good) or not adequate (including the values of poor and insufficient), in the Portuguese sample, 40.73% and 41.14% of the students had neither adequate knowledge nor adequate skills, compared to 52.27% and 55.19% of the Spanish students in knowledge and skills respectively, but with better results than those found in English students where only 22.41% and 33.62% had adequate knowledge and skills or the data found in Chinese nursing professionals with values of 19.56% and 19.93% in terms of representation of professionals with good knowledge and skills. The values suggest that nursing professionals possess adequate knowledge and skills. Efforts are therefore required to improve aspects related to environmental health and to comply with professional standards established and defined by some bodies such as the International Council of Nurses [
31]. As shown in our study, more general issues, such as those related to vulnerable populations like children or air pollution, are known but more specific pesticides that may cause neoplasms are unknown.
In relation to the SANS_2 questionnaire, both reliability and validity have been assessed in the population following classical measurement theory. In order to provide a more valid and reliable factor solution, the exploratory factor analysis took into account two methods for extracting the factors, the choice of those with eigenvalues ˃1 (Kaiser's criterion), and Horn's parallel analysis [
32]. The results obtained from the factor loadings in terms of construct validity in the confirmatory analysis, although lower in one of them than in other studies conducted in nursing professionals [
15] generally indicate an adequate structure for the proposed measurement model. The model fit indices show data that do not reach values within the accepted ranges, such is the case of RMSEA and SRMR. This may be due to the inclusion in the sample of students from different educational levels, including first-year students whose training is limited to basic subjects, as it may lead to diverse responses that do not adequately reflect the fit in the proposed model. The reliability measured with Cronbach's alpha indicates acceptable values reaching a value of 0.84, similar to those obtained in other studies of nursing students and professionals [
13,
15], higher than the value of 0.77 obtained in other research carried out in China [
33] and below the results found (α = 0.926) in a sample of 159 American nursing students and teachers [
34].
Attitude is defined as ‘a disposition to respond favorably or unfavorably to an object, person, institution or event’ [
35]. According to the theory of planned behavior, if attitudes affect behavior, nursing students who have a more positive attitude are more likely to engage in sustainable behaviors [
36]. The mean score of the attitude scale in Portuguese nursing students was 28.15 ± 4.61, with 80.13% of the surveyed students scoring excellent and good attitude. Several factors appear to influence differences in attitudes, for example, Cruz et al. [
37] found significant differences between attitude scale scores in students from Saudi Arabia, Egypt, Iraq, and the Palestinian Territories, and found that country of residence, rural/urban location, and climate change education were variables that could contribute to such differences. Other variables, such as political affiliation, could influence attitude scores [
38]. Some research, like ours, indicates the analysis of attitude scores between first-year students and students in more advanced grades, with differences between the two results [
34]. This could be because first-year students would not have received training and education in sustainability and environmental issues, which students in higher grades would see included in more advanced and specific subjects. It also seems that the male gender has a higher competence, according to the results of previous studies [
17] and the present one, although the samples of male nursing students are very low and this could affect the statistical results found, so this aspect should be further studied, since in most of the studies reviewed there is no analysis with a gender perspective.
The complexity that arises as to how knowledge, attitudes, or skills might influence how people might behave toward issues related to environmental sustainability has led some authors to put forward theories in this regard. The value-belief-norm theory has offered the best explanation for supporting the environmental movement. In this theory, attitudes, environmental knowledge, and values are the central elements in pro-environmental behavior [
39]. This behavior will be driven by an individual's awareness of environmental issues as well as environmental values. Educational interventions that teach about environmental management and sustainability are required so that the users who receive them are not only able to understand but are also able to initiate change [
40]. This supports the results of our study, which shows the transformation of learning from environmental sessions into skills and the integration of knowledge in the long term.
In the current literature, there is a growing body of research using educational interventions that significantly increase attitudes toward climate change and sustainability [
41,
42,
43,
44,
45]; however, sometimes, although this training involves an improvement in attitudes, it is found that students have difficulties in applying sustainability and addressing unsustainable practices in future work environments [
41]. Therefore, this education should reinforce the empowerment and leadership of these future nursing professionals as promoters of change.
In spite of the described strengths of our study, it has limitations because, although the study sample is large, the student body is limited to a single training center.
Author Contributions
Conceptualization, Cristina Álvarez-García, Beatriz Edra, Goreti Marques, Catarina Simões and Mª Dolores López-Franco; Data curation, Cristina Álvarez-García and Mª Dolores López-Franco; Formal analysis, Cristina Álvarez-García and Mª Dolores López-Franco; Funding acquisition, Cristina Álvarez-García and Mª Dolores López-Franco; Investigation, Cristina Álvarez-García, Beatriz Edra, Goreti Marques, Catarina Simões and Mª Dolores López-Franco; Methodology, Cristina Álvarez-García, Beatriz Edra, Goreti Marques and Mª Dolores López-Franco; Project administration, Cristina Álvarez-García, Beatriz Edra, Goreti Marques, Catarina Simões and Mª Dolores López-Franco; Resources, Beatriz Edra, Goreti Marques and Catarina Simões; Software, Cristina Álvarez-García and Mª Dolores López-Franco; Supervision, Cristina Álvarez-García, Beatriz Edra, Goreti Marques, Catarina Simões and Mª Dolores López-Franco; Validation, Cristina Álvarez-García, Beatriz Edra, Goreti Marques, Catarina Simões and Mª Dolores López-Franco; Visualization, Cristina Álvarez-García, Beatriz Edra, Goreti Marques, Catarina Simões and Mª Dolores López-Franco; Writing – original draft, Cristina Álvarez-García and Mª Dolores López-Franco; Writing – review & editing, Cristina Álvarez-García, Beatriz Edra, Goreti Marques, Catarina Simões and Mª Dolores López-Franco.