Health and Well-Being in Protected Natural Areas—Visitors’ Satisfaction in Three Different Protected Natural Area Categories in Catalonia, Spain

Protected natural areas (PNAs) can be a source of health and well-being, but little research has been carried out regarding outcomes in terms of satisfaction (the difference between motivations and benefits). Inspired by previous research that examines the motivations and benefits perceived by visitors to various PNAs in Canada, and based on importance–performance analysis (IPA) and service quality gap (GAP) analysis theory, the aim of this study was to identify the outcomes generated by protected areas in terms of satisfaction, especially with regard to the PNAs’ individual protection categories. The study was based on survey data from visitors (n = 360) to three PNAs in Catalonia: one national park, one natural park and one periurban park. The results indicate that anticipated environmental, psychological, physical and social benefits were of major personal value in choosing to visit a PNA. The results indicate that, generally, visitors were satisfied with regard to the benefits anticipated. Differences between parks in this respect could be explained in part by sociodemographic factors and visitors’ behavior. The results are discussed in terms of their applicability and how they relate to the role of PNAs in the promotion of visitors’ health and well-being.


Introduction
At the global level, visits to protected natural areas (PNAs) have been increasing since the 1990s, reaching unprecedented numbers [1]. According to Europarc-Spain's statistics, 15 million people visited Spanish national parks in 2017, in contrast to only 8 million in 2007 [2]. In the case of Catalonia, depends on comparing visitors' expectations before visiting the PNA, and the experiences and images they come into contact with.
In the context of measuring the quality of services as perceived by visitors or visitor satisfaction in a tourist destination, including a number of PNAs, several approaches were developed. These approaches were based on two main types of analysis: importance-performance analysis (IPA) and gap analysis of service quality (GAP) [33][34][35][36][37][38][39]. The first approach combines a measure of performance with its associated importance in a two-dimensional grid. It provides a graphic representation of the performance of managers, suppliers or operators in providing a range of services. The second approach measures service performance as the difference between expected performance and an evaluation of the services as actually experienced. It is important to highlight the research carried out by Tonge and Moore [38] and Rice et al. [39], who reconceptualized IPA and GAP first to assess satisfaction [38], and more recently to examine motivation outcomes [39]. Tonge and Moore [40] replaced the performance element of IPA with satisfaction, and the GAP relied on the (statistical) means of importance and satisfaction. For Rice et al. [39], on the other hand, motivation is distinct from expectation-i.e., one's perceived likelihood of attaining the outcomes or benefits sought through a recreation experience. A negative statistically significant gap, where the importance/motivation mean is larger than the satisfaction/outcomes mean, suggests that management action is required. Conversely, a positive statistically significant gap, where the importance/motivation mean is lower than the satisfaction/outcomes mean, suggests no extra management measures are required.
In our framework, which takes these new approaches into account, the original IPA (importance-performance analysis) is conceptualized as Motivation and Benefits (see Figure 1). For instance, the Importance may be to relieve stress or to do physical activity. If successful, the benefit of the recreational experience would be the attainment of a relaxed state or to have an active visit. individuality of a particular PNA. Basically, for these authors, the evaluation of visitor satisfaction depends on comparing visitors' expectations before visiting the PNA, and the experiences and images they come into contact with.
In the context of measuring the quality of services as perceived by visitors or visitor satisfaction in a tourist destination, including a number of PNAs, several approaches were developed. These approaches were based on two main types of analysis: importance-performance analysis (IPA) and gap analysis of service quality (GAP) [33][34][35][36][37][38][39]. The first approach combines a measure of performance with its associated importance in a two-dimensional grid. It provides a graphic representation of the performance of managers, suppliers or operators in providing a range of services. The second approach measures service performance as the difference between expected performance and an evaluation of the services as actually experienced. It is important to highlight the research carried out by Tonge and Moore [38] and Rice et al. [39], who reconceptualized IPA and GAP first to assess satisfaction [38], and more recently to examine motivation outcomes [39]. Tonge and Moore [40] replaced the performance element of IPA with satisfaction, and the GAP relied on the (statistical) means of importance and satisfaction. For Rice et al. [39], on the other hand, motivation is distinct from expectation-i.e., one's perceived likelihood of attaining the outcomes or benefits sought through a recreation experience. A negative statistically significant gap, where the importance/motivation mean is larger than the satisfaction/outcomes mean, suggests that management action is required. Conversely, a positive statistically significant gap, where the importance/motivation mean is lower than the satisfaction/outcomes mean, suggests no extra management measures are required.
In our framework, which takes these new approaches into account, the original IPA (importance-performance analysis) is conceptualized as Motivation and Benefits (see Figure 1). For instance, the Importance may be to relieve stress or to do physical activity. If successful, the benefit of the recreational experience would be the attainment of a relaxed state or to have an active visit. The aim of the present study was to compare visitor satisfaction in terms of health and wellbeing motivations and perceived benefits at three different PNAs (one national park, one natural park and one periurban park), in order to identify any differences between different types of PNAs in their capacity to contribute to visitors' health and well-being. The key research questions addressed in this study follow: (Q1) To what extent did visitors to different categories of PNA differ in terms of their sociodemographics and visit behaviors?
(Q2) What were the most important health and well-being motivations and benefits (outcomes) identified by visitors?
(Q3) To what extent did motivations and benefits differ among PNAs? (Q4) To what extent did each individual PNA satisfy the health and well-being motivations pursued by visitors (i.e., GAP analysis)?
(Q5) For each PNA, were there any differences in the gaps between satisfaction and motivation that can be explained by sociodemographic characteristics and visit behaviors? The aim of the present study was to compare visitor satisfaction in terms of health and well-being motivations and perceived benefits at three different PNAs (one national park, one natural park and one periurban park), in order to identify any differences between different types of PNAs in their capacity to contribute to visitors' health and well-being. The key research questions addressed in this study follow: (Q1) To what extent did visitors to different categories of PNA differ in terms of their sociodemographics and visit behaviors?
(Q2) What were the most important health and well-being motivations and benefits (outcomes) identified by visitors?
(Q3) To what extent did motivations and benefits differ among PNAs? (Q4) To what extent did each individual PNA satisfy the health and well-being motivations pursued by visitors (i.e., GAP analysis)?
(Q5) For each PNA, were there any differences in the gaps between satisfaction and motivation that can be explained by sociodemographic characteristics and visit behaviors?

Study Areas
The study was carried out in three different PNAs in Catalonia: one national park, one natural park and one periurban park. Figure 2 shows their geographical locations: Aigüestortes i Estany de Sant Maurici National Park (A), Alt Pirineu Natural Park (P), and Serra de Collserola Periurban Park (C). The main reasons for choosing these three PNAs were that (1) all are Natura 2000 sites, (2) all are managed by public organizations, and (3) although they all provide similar opportunities in terms of recreational activities, they present different physical (i.e., topographical), social and managerial settings.

Study Areas
The study was carried out in three different PNAs in Catalonia: one national park, one natural park and one periurban park. Figure 2 shows their geographical locations: Aigüestortes i Estany de Sant Maurici National Park (A), Alt Pirineu Natural Park (P), and Serra de Collserola Periurban Park (C). The main reasons for choosing these three PNAs were that (1) all are Natura 2000 sites, (2) all are managed by public organizations, and (3) although they all provide similar opportunities in terms of recreational activities, they present different physical (i.e., topographical), social and managerial settings. According to the recreation opportunity spectrum proposed by Brown et al. [40], Aigüestortes i Estany de Sant Maurici National Park is semiprimitive nonmotorized (SPNM); Alt Pirineu Natural Park is semiprimitive motorized (SPM), and Serra de Collserola Periurban Park is urban (U).
Aigüestortes i Estany de Sant Maurici National Park, in the Pyrenees, was established in 1955 and extended in 1996, and is the only national park in Catalonia. Alt Pirineu Natural Park, which is also in the Pyrenees, is the largest natural park in Catalonia, with an area of over 79,300 ha. Serra de Collserola, although a natural park, is one of the largest periurban parks in Europe. It receives around 5,000,000 visits per year (Table 1). According to the recreation opportunity spectrum proposed by Brown et al. [40], Aigüestortes i Estany de Sant Maurici National Park is semiprimitive nonmotorized (SPNM); Alt Pirineu Natural Park is semiprimitive motorized (SPM), and Serra de Collserola Periurban Park is urban (U).
Aigüestortes i Estany de Sant Maurici National Park, in the Pyrenees, was established in 1955 and extended in 1996, and is the only national park in Catalonia. Alt Pirineu Natural Park, which is also in the Pyrenees, is the largest natural park in Catalonia, with an area of over 79,300 ha. Serra de Collserola, although a natural park, is one of the largest periurban parks in Europe. It receives around 5,000,000 visits per year (Table 1). One of the most important characteristics of Aigüestortes i Estany de Sant Maurici National Park is its relief-u-shaped valleys, glacial cirques, ponds, gorges and screes-a legacy of the erosive action of Quaternary glaciers, mostly on granite, but also slates, shales and other rocks. Alt Pirineu Natural Park, located near the Aigüestortes i Estany de Sant Maurici National Park, was established in 2003 and extended in 2018. It includes the highest peak of the Catalan Pyrenees and is one of the most popular hiking areas in the Pyrenees, suitable for both day-trippers and overnight hikers. Finally, Serra de Collserola Natural Park is located in the middle of one of the densest urban areas on the Mediterranean shore, Barcelona (1.6 million residents). It was established in 1987 by the Catalan government, covers an area of 8259 hectares (17 km long × 6 km wide) and is characterized by an extensive network of recreational trails that are used for multiple activities, including hiking, mountain biking and running.

Data Collection and Sampling Strategy
Data were collected using 360 structured questionnaires completed onsite, 120 in each PNA. Regarding sample strategy, the interviews were held from April to July 2016 and took place at popular locations in the parks (the main entrance, the main hut and the main landscape attraction point). The questionnaires were conducted by 21 trained research staff. Potential respondents were approached randomly at the main entrances on their way out, because most of the questions referred to the experience that the visitors had just had, during their visit. The average time spent on the interview was 8-10 min and participants did not receive any financial compensation for participating in the study. The response rate was 96%.
The final sample showed similar sociodemographic data to other visitors studied in previous research by Farías and Morera in Serra de Collserola Periurban Park [41] and Alt Pirineu Natural Park [42]. Sociodemographic outcomes such as gender and age among this and previous studies were below 7% of difference (1.5% and 2.7% for gender and age, respectively, for Serra de Collserola Periurban Park and 6% and 6.6% for Alt Pirineu Natural Park). Complete data about the results obtained are shown in Section 3.

Questionnaire
The survey was in three parts. First, visitors were asked about sociodemographic characteristics (e.g., age, place of residence). In the second part, they were asked about recreational behavior (e.g., frequency of visit, duration of their stay, main activity practiced during the visit). Finally, visitors were asked to rate the importance of each of six motivations and their satisfaction level or outcomes for each of six possible benefits, as suggested by Lemieux et al. [22].
The ratings for motivations used a 5-point Likert scale ranging from 1 (not at all important) to 5 (very important). The motivations were physical (the opportunity to do physical activity), psychological (the opportunity to relax, or to recover from stress or mental fatigue), social (the opportunity for increased social interaction and to meet new people), spiritual (the opportunity to connect with nature or look for inspiration), environmental (the opportunity to be outdoors, connect to nature or experience a sense of place) and intellectual (the opportunity to engage in creative and stimulating activities, or to acquire new knowledge).
Where benefits were concerned, visitors were asked to what extent they thought that the visit had increased their physical, psychological, social, spiritual, environmental and intellectual well-being (from 1 = did not improve at all, to 5 = improved greatly). The benefits were defined as follows: physical (perceived physical well-being), psychological (recovery from mental stress, improved mood), social (strengthened social relations, meeting new people), spiritual (finding inspiration, connecting with nature), environmental (enjoying being outdoors, experiencing a sense of place, fostering a harmonious human-nature relationship) and intellectual (learning something new, enjoying participating in creative activities).

Data Analysis
The data collected were analyzed using the Statistical Package for the Social Sciences 18.0 (SPSS, Chicago, IL, USA). First, preliminary analyses were carried out, including the identification of missing values and the distribution of the data, and descriptive statistics for each of the study's variables were defined. Secondly, we tested for possible differences between sociodemographic characteristics and visit behavior in the three PNAs. Chi-square goodness-of-fit at a 95% accuracy level was used to examine the differences between visitors with respect to categorical variables. Thirdly, since the data were not normally distributed, nonparametric Kruskal-Wallis and Mann-Whitney U tests were performed to explore the differences in motivations and benefits between visitors to the three different PNAs. Differences between the three parks were analyzed using Kruskal-Wallis tests, and differences between pairs of parks using Mann-Whitney tests. In the latter analysis, the level of significance was adjusted to p < 0.017 (0.05/3). Fourthly, in order to test the differences between visitors' motivations in visiting the PNA and their satisfaction after their visit, we calculated the gap value for each motivation and benefit. Gaps were estimated by subtracting the mean for motivation (importance) from the mean for benefits (satisfaction). Six gaps were calculated for each PNA, each of them referring to a different motivation/satisfaction pair (e.g., psychological). A negative gap (motivation exceeding benefits) suggests that satisfaction is negative and additional management action is required. Conversely, a positive gap (benefits exceeding motivation) suggests that satisfaction is positive. To test whether the gaps differed between PNAs, we also used Kruskal-Wallis and Mann-Whitney tests, following the same rationale as the one described above. Finally, for each PNA we tested whether characteristics of visit behavior were related to differences in the gaps between satisfaction and motivation. Kruskal-Wallis and Mann-Whitney tests were also used in this analysis.

Visitor Profiles and Recreational Behavior
As shown in Table 2, the most frequent visitors to the three areas were men (58%), middle-aged (40.4%), and held a university degree (57.5%). The educational level of the sample in the study was 10% higher than the average level for the Catalan population as a whole (last available census, 2017) [43]. Over two-thirds of the sample had visited the parks at least four times in the previous two years. People who were visiting the parks as part of a group most commonly came with a group of friends (43%); mostly spent about half a day in the park (78.3%); and during their visit, 60.5% undertook recreational hiking (slow walking) or hiking (brisk walking). Values are percentages of the number of individuals for each park. * Significant differences (Chi-square or Kruskal-Wallis test).
Regarding differences in sociodemographic characteristics and recreational or visit behavior between parks (Table 2), the results revealed significant differences in gender (p = 0.033), age (p = 0.021), frequency of visit (p < 0.001), group composition p < 0.001), duration of visit (p < 0.001) and recreational activity (p < 0.001). No significant differences were found concerning educational level (p = 0.094) and health perception (p = 0.286).
Specifically, visitors to Aigüestortes i Estany de Sant Maurici National Park were often women (50.8%) and older adults (41.6% were over 55), who were visiting the area for the first or second time (44.5%), in couples or family groups (57.5%). Almost 51% did slow walking (recreational hiking).
Visitors to the Serra de Collserola Periurban Park were characterized by being male (65.5%) and/or young (almost the 20% were aged [26][27][28][29][30][31][32][33][34][35], visit the area frequently (73.4%), came with a group of friends (40.8%), for half a day (99.2%). They carried out more intense physical activity than the visitors to the other two parks: more than 47% of them did trail running or mountain biking during their visit. Significantly, the second most common "group" composition was visitors who were there alone (26.7%), a much higher figure than for the other two PNAs.
Finally, visitors to the Alt Pirineu Natural Park, who had an intermediate profile, were mostly male (58.8%) and middle-aged (54.7%), and 50% accompanied by a group of friends. They also remained in the area longer than visitors in the other two parks (45%) and performed less physical activity. More than the 33% of the visitors to this park remained in or near the entrances, meaning that they did not practice any type of physical activity during their visit.

Health and Well-Being Motivation for Visiting Protected Natural Areas
As shown in Table 3

Motivation and Benefits Differences among Parks
As shown in Table 5, the results of Kruskal-Wallis tests revealed that the parks differed significantly for three of the six motivations relating to health and well-being. Interesting differences were observed for physical, social and intellectual motivations. Visitors to Serra de Collserola Periurban Park showed higher scores for physical motivation than those to Aigüestortes i Estany de Sant Maurici National Park and Alt Pirineu Natural Park (p < 0.001). Similarly, visitors to Aigüestortes i Estany de Sant Maurici National Park and Alt Pirineu Natural Park exhibited higher scores for social motivation compared to visitors to Serra de Collserola Periurban Park (p < 0.001). Furthermore, visitors to Aigüestortes i Estany Sant Maurici National Park presented higher scores for intellectual motivation than visitors to Alt Pirineu Natural Park and Serra de Collserola Periurban Park (p < 0.001). No statistically significant differences were observed for psychological, spiritual and environmental motivations, two of the three most important motivations reported in the set of the parks (environmental and psychological). Significant differences between parks were identified for four of the six underlying benefits: physical, social, environmental and intellectual (Table 6). Visitors to Aigüestortes i Estany de Sant Maurici National Park showed higher scores for environmental (p = 0.002) and intellectual (p < 0.001) benefits than those to Serra de Collserola Periurban Park and Alt Pirineu Natural Park. Likewise, Alt Pirineu Natural Park presented higher scores for social benefits (p < 0.001) and Serra de Collserola Periurban Park for physical benefits (p < 0.001).

Health and Well-Being Satisfaction Pursued by Visitors to Different PNAs
This subsection presents the differences (i.e., gaps) between PNAs for the six health and well-being dimensions. For the former, the results (see Table 7) revealed that the calculated gaps differed across parks only for the physical and spiritual dimensions. The two highest values of these dimensions were registered in Aigüestortes i Estany de Sant Maurici National Park, while the lowest values for these two dimensions were identified in Alt Pirineu Natural Park. However, a positive gap (higher mean for benefits than for motivation) was observed across all three parks for these same two dimensions, as well as for the psychological and physical dimensions. When further analyses were carried out, no negative gap was found in Aigüestortes i Estany de Sant Maurici National Park, while several were found in the other two PNAs-two in the case of Alt Pirineu Natural Park (environmental and intellectual), and one in Serra de Collserola Periurban Park (intellectual). Of these, only one, according to IPA theory-the intellectual dimension in Serra de Collserola Periurban Park (the value for which was outside acceptable margins)-would require intervention by management (Figure 3).

Influence of Sociodemographics and Recreation Behaviors on GAP Results
When examined in relation to sociodemographic and visit behavior, the gap values did not vary by age or educational level. However, several gap differences were obtained for the remaining independent variables (i.e., gender, health perception, frequency of visit, visitor group, duration of visit, and activity; see Table 8). Overall, eleven significant differences were found in Aigüestortes i Estany de Sant Maurici National Park (i.e., gender, frequency of visit, group composition, and activity); one in Alt Pirineu Natural Park (i.e., duration of visit), and four in Serra de Collserola Periurban Park (i.e., group composition, duration of visit, and activity). Subsequent post hoc analysis found six significant differences in Aigüestortes i Estany de Sant Maurici National Park, one in Alt Pirineu Natural Park, and two in Serra de Collserola Periurban Park, as the p value was corrected for multiple comparisons (see Appendix A). In the case of Aigüestortes i Estany de Sant Maurici National Park, differences were found according to gender

Influence of Sociodemographics and Recreation Behaviors on GAP Results
When examined in relation to sociodemographic and visit behavior, the gap values did not vary by age or educational level. However, several gap differences were obtained for the remaining independent variables (i.e., gender, health perception, frequency of visit, visitor group, duration of visit, and activity; see Table 8). Overall, eleven significant differences were found in Aigüestortes i Estany de Sant Maurici National Park (i.e., gender, frequency of visit, group composition, and activity); one in Alt Pirineu Natural Park (i.e., duration of visit), and four in Serra de Collserola Periurban Park (i.e., group composition, duration of visit, and activity). Subsequent post hoc analysis found six significant differences in Aigüestortes i Estany de Sant Maurici National Park, one in Alt Pirineu Natural Park, and two in Serra de Collserola Periurban Park, as the p value was corrected for multiple comparisons (see Appendix A). In the case of Aigüestortes i Estany de Sant Maurici National Park, differences were found according to gender

Discussion
Using three different PNAs as case studies, this research set out with the aim of answering five questions. The results obtained provide an interesting foundation for detailed discussion.

(Q1) To what extent did visitors to different categories of PNA differ in terms of their sociodemographics and visit behaviors?
The results obtained revealed certain patterns of characteristics and recreational behavior in common for all three parks, for all variables with the exception of educational level. However, each park presented significant differences in terms of visitor profiles. According to earlier studies (e.g., by Farías-Torbidoni; Baric, Anic and Bedoya; Múgica and Delucio or Delucio and Múgica; [44][45][46][47][48]), Aigüestortes i Estany de Sant Maurici National Park could be characterized as attracting visitors with a more general profile than the other two PNAs, and Serra de Collserola Periurban Park as being visited by people with a more sports-oriented profile. Here, it is important to highlight that Aigüestortes i Estany de Sant Maurici National Park and Serra de Collserola Periurban Park are situated at opposite ends of the recreational opportunity spectrum (referred to in Section 2): semiprimitive nonmotorized in the case of the former and urban in the case of the latter.
The results support the outcomes found previously regarding the influence of PNA characteristics (physical, social and managerial) not only on visitor profiles, but also on the main recreational activities carried out in the PNAs. As mentioned in Section 3 (Results), while visitors to Aigüestortes i Estany de Sant Maurici National Park were often first-or second-time visitors, in couples or families, there for half a day and doing recreational hiking (slow walking), visitors to Serra de Collserola Periurban Park tended to visit frequently, individually and to carry out higher levels of physical activity (moderate to vigorous).

(Q2) What were the most important health and well-being motivations and benefits (outcomes) identified by visitors?
The results obtained for main motivations and benefits were similar to those of earlier studies carried out in Canada [23,32], Finland [24] and, indeed, Spain [27]. It is worth underlining this coincidence with other studies regarding the four main motivations and benefits identified, although their order of significance may differ, for instance, physiological, social, physical and environmental aspects in the case of Canada, and environmental, psychological, physical and social in our case. Although the results are consistent, when we assess the differences between studies one limitation emerges, which is connected to the possible misunderstanding of some items. For example, while physical, psychological and social motivations and benefits could be clearly identified by users, this was not the case for environmental, spiritual and intellectual motivations or benefits. It is possible that respondents found it difficult to identify and interpret aspects such as connecting with nature or enjoying the outdoors (environmental), seeking inspiration (spiritual) or having the opportunity to participate in creative and stimulating activities in order to acquire new knowledge (intellectual). Thus, it is possible that not all the answers reflect the same interpretations.

(Q3) To what extent did motivations and benefits differ among PNAs?
The results obtained regarding this question not only support the existence of significant differences between parks, but also allow us to identify some specific well-being roles (related to the identity and character of each PNA) attributed by visitors to the parks. Although there is a perfect match between the three highest values for motivations and the perceived benefits (outcomes) of visiting each park, there are significant differences in the specific motivations and benefits associated with each park. For instance, environmental and intellectual benefits are clearly associated with Aigüestortes i Estany de Sant Maurici National Park, social benefits with Alt Pirineu Natural Park, and physical benefits with Serra de Collserola Periurban Park. However, these trends need to be carefully considered because the research did not go into depth with all possible associations; for instance, Serra de Collserola Periurban Park is visited weekly by visitors and their results show low scores for the intellectual dimension, which may be related with the high frequency of visits. These results are especially significant if we associate them with the visitor profiles identified by Q1. Furthermore, the results obtained in the case of Serra de Collserola Periurban Park (a more sports-oriented profile, with greater importance attached to physical motivation and benefits) again support the important role that this kind of area can play in the promotion of HEPA [14,28,31,49,50].

(Q4) To what extent did each individual PNA satisfy the health and well-being motivations pursued by visitors (i.e., Gap analysis)?
The data obtained show that generally, in this set of PNAs, there is a concordance between what visitors want (motivation) and what they get (benefits). In the case of these three parks, only three gap values emerged (i.e., where benefits derived by visitors did not meet with their motivations): the environmental and intellectual dimensions in Alt Pirineu Natural Park, and the intellectual dimension in Serra de Collserola Periurban Park. From IPA and GAP points of view, the intellectual dimension in Serra de Collserola Periurban Park is the only dimension in any of the three parks for which additional management action is required if we want to satisfy all visitors completely. However, if managerial actions aspire to increase the positive gap value (i.e., to exceed visitors' expectations), a more detailed study will be needed. As a starting point, we can mention that in our case, the highest positive gap value was found in the PNA with the highest protection category (Aigüestortes i Estany de Sant Maurici National Park; see Table 1), more pristine areas, and a low level of human intervention. If we were able to identify some specific features, beyond the natural characteristics of the PNA, that increased the positive gap, managers could be in a position to adapt their offers to make the parks more suitable and attractive to their target audiences where health motivations are concerned, thereby increasing health benefits.

(Q5) For each PNA, were there any differences in the gaps between satisfaction and motivation that can be explained by sociodemographic characteristics and visit behaviors?
In general, the differences found in our study can be explained by visitor behavior rather than by sociodemographic characteristics. For instance, only two of the sociodemographic variables (gender and health perception) differ in relation to level of satisfaction, whereas all five recreational behavior variables differ. While we do not have any specific results from other studies with which to compare our own directly, our results could relate to those of Lemieux et al. [22] or Romagosa [27], who found that women rather than men often perceived greater well-being benefits derived from visits. As it was expected according to classic studies targeting low-physical activity individuals and health outcomes [51,52], our results also show that visitors whose health perception was lower obtained indeed greater physical and social health benefits after doing physical activity in the PNA. Another important finding was the link between some kind of physical activity during the visit and higher levels of satisfaction in both the physical and the spiritual dimensions. However, more data are needed to establish any connections between sociodemographic characteristics or recreational behavior, and health and well-being satisfaction.

Study Limitations and Managerial Recommendations
The number of respondents is reliable and representative, as previously mentioned in data collection and sampling strategy section. One of our limitations was that data collection was performed during spring and early summer (April to July) and extended time of field work may identify different visitor profiles.
From our point of view, future studies should include more than one PNA according to their protection categories and even from different geographically and cultural settings, including, preferably, collecting data throughout the year.
Managerial recommendations specifically addressed to a particular PNA should first identify the assets of the PNA, to link the PNA features with specific health dimensions (e.g., spiritual health related with the PNA context) or public health policies and actions (e.g., www.eupap.org).
Finally, data obtained from such studies may be used for the development of a comprehensive monitoring, evaluation and reporting framework for the effectiveness of PNA management on health and well-being promotion.

Conclusions
Previous studies have analyzed and discussed the roles/benefits of PNAs in relation to human health and well-being. However, the real outcomes generated by these areas in terms of satisfaction (the difference between motivations and benefits) have not been examined empirically. This is the contribution of this study, which is based on the application of IPA and GAP analysis. The results obtained generally show that in our set of PNAs there is a concordance between what visitors want (motivation) and what visitors get (benefits), and consequently high visitor satisfaction levels. However, important differences between PNAs were found regarding human health and well-being motivations for visits and the expected benefits. Although these differences could be explained in part by sociodemographic characteristics and recreational behavior, looking at each area's identity or character as perceived by the visitor opens an interesting perspective. In this case, the specialization of each park with regard to some of the specific dimensions analyzed here could be a good way to optimize resources. Finally, the results obtained indicate, even if inconclusively, a link between PNA category and level of satisfaction, a link that would merit much more detailed study.      Note: For each score, M (SD) values with different superscripts indicate a significant difference using U Mann-Whitney tests as post hoc comparisons. P = Poor or fair, G = Good, V = Very good or excellent. The level of significance for the U Mann-Whitney tests was adjusted to p < 0.017 (p = 0.05/3).