Factors Associated with Spa Tourists’ Satisfaction

Most spas are surrounded by beautiful natural environments, so these facilities contribute to sustainable economic development and their guests are national and local tourists with a high average stay. This study has the objective of analyzing the factors associated with the level of satisfaction and the profile of spa tourists in Andalusia, southern Spain. The methodology used was a quantitative study based on a stratified cross-sectional survey. Each stratum was a spa, with a population of 53,231 users per year, whose final sample included 725 users. The results found were that the profile of Andalusian spa tourists is predominantly female, over 56 years old, and retired with average incomes below one thousand euros. The level of satisfaction of its users is above the average in the factors of accommodation, restaurants, spa, and public roads and natural environment and most of these clients would recommend it to others [92.4%] and express the intention to return [63.6%] which indicates that Andalusian spas show a very high level of customer loyalty. The study has the limitation of being cross-sectional quantitative research carried out in Andalusia that could be complemented in future research with a qualitative study with experts. The research presents contributions for companies since they could develop strategies to attract users with younger profiles with activities, prices, or spaces that are more appropriate to their needs and seek the conditions to increase the satisfaction of these segments. The study is original since this research is a quantitative study that shows the real profile of the Andalusian spa tourist and their high level of satisfaction with the spa and their loyalty. The surveys were conducted at each spa by a face-to-face researcher with a response rate of 98%.


Introduction
Thermal springs have long attracted flows of people primarily seeking to improve their health [1,2]. These mineral waters' benefits and properties have been combined with a new understanding of health tourism that requires higher-quality specialized offers that meet increasingly demanding consumers' needs [3][4][5][6].
Some authors confirm that tourism experiences connected to water correlate positively with the perception of the quality of life, satisfaction, and loyalty of people both towards the experience and towards the destination [7].
Internationally, according to the Global Wellness Institute [8], the hot springs market is an important business segment of the wellness tourism industry. More specifically, the spa sector in 2017 included more than 149,000 facilities, with $93.6 billion in revenue, and 2.5 million workers. The spa sector grew 9.9% annually between 2015 and 2017, and the forecast for 2022 is that earnings will reach $128 billion [8,9]. In addition, international wellness tourists spend an average of $1528 per trip, 53% more than the typical international tourist. Likewise, national wellness tourists spend $609 per trip, 178% more than the average national tourist [10]. In addition, health tourism is an expanding phenomenon [11] and one of the fastest growing forms of tourism at the moment [12] although according to the UNWTO (United Nations World Tourism Organization), due to the Covid-19 pandemic, tourism will take two or three years to recover those levels [13].
In this context, visits to thermal spring facilities and other derivatives such as spas are an increasingly common practice and there are studies that indicate that spa tourists seek a better quality of life and flee from mass tourism [14].
In the approach to the research topic, it was found that Andalusia had thermal stations that, precisely due to the characteristics of its waters in terms of temperature and minerals, act as a focus of attraction for tourists [15] and the tourism of spas in Andalusia was studied, which from the tourist point of view, is not as studied as sun and beach tourism.
Thus, the objective of this research was to analyze the factors associated with the level of satisfaction and the profile of spa tourists in Andalusia, southern Spain based on a stratified cross-sectional survey in which each stratum was a specific spa. The estimated study population was 53,231 spa users per year, and the final sample included 725 respondents.
This paper is structured as follows. Section 2 reviews the relevant literature. Section 3 covers the methodology, data collection procedures, and questionnaire reliability. The results, discussion, and conclusions are presented in Section 4, Section 5, and Section 6, respectively.

Literature Review
The global boom in health tourism [8] has led destinations and companies to specialize in health and wellness tourism as part of their corporate strategy. For some destinations, health tourism can be a way to combat seasonality and diversify coastal tourism [16]. Tourism has also contributed to the sustainable development of rural destinations.
Academic research has studied the factors that attract spa tourists. An example of this is the market segmentation study carried out with a survey of 2050 spa tourists in Poland [17] which has shown that demand is significantly dependent on consumers' choice of particular forms of therapeutic spa treatments based on multiple factors such as age, education, size of tour groups, past spa experiences, type of guest, length of stay, and expected benefits [18]. People use relaxation and tourism to enhance and manage their physical wellbeing in accordance with new concepts and dimensions of health, which the market has taken advantage of, significantly diversifying the offers and incorporating a new range of services. These take a holistic approach to the body that seek to take into account the quite different motivations of individual consumers [19,20].
Other studies have revealed that the feeling of enjoyment and better health is intensified by the wellbeing provided by the environment, the interaction with other users, and the staff who attend them [21].
In addition, other investigations add that in addition to the therapeutic qualities of the waters, the social aspects and the desire to know the geological elements of these places act as an attraction factor [22].
Spa companies offering a wide range of services are looking for new customer value propositions to attract new market segments as, until now, spas have tended to focus on their natural resources [23].
A study of spa tourism in Thailand found that the quality and satisfaction of tourists act as forces that promote loyalty to the destination [24].
Other studies have also shown that satisfaction is positively correlated with repeat purchases of services, loyalty, and intentions to recommend products to others [25].
Some research has concentrated on identifying which key destination attributes contribute to specific tourists' overall satisfaction, such as individuals who visit Mallorca in the low season, and analyzing these tourists' loyalty to the relevant destinations [26,27]. The existing literature focuses only on controllable attributes, so uncontrollable factors such as climate or landscape have not been taken into account because they are difficult to change when tourists are dissatisfied [28]. Conclusions drawn from these studies include that the attractions offered by destinations, facilities, activities, and services have a great impact on visitors' intention to return and recommend [26][27][28].
Tourists who repeatedly visit a destination are overall more satisfied than those who come for the first time. Younger tourists, in turn, are less satisfied than older tourists. Researchers have also concluded that destination recommendation's impact is more im-portant than intention to return since the latter is subject to multiple tourist motivations [26,27]. Many studies further link customer satisfaction with tourism establishments' location [29][30][31][32], so location is associated with these companies' survival [33][34][35].
The literature connecting tourist satisfaction and quality of service in tourism destinations is extensive [36][37][38][39]. For instance, customers can be satisfied with some services and not with others [40]. Other authors have focused on satisfaction as a construct that depends on perceived quality [18]. Given these findings, the present study's aims included verifying thermal spa tourism's contribution to the surrounding destinations' overall perceived quality as spa facilities can be important tourist attractions.
Relevant research has been carried out on quality in sun and sand tourism destinations specifically in the Costa del Sol [41]. The cited author designed a methodology and system of quality indicators that were then applied to three mature sun and beach tourism destinations-Torremolinos, Benalmádena, and Marbella. Otero Cordones [41] adapted previous methodologies used in studies evaluating the quality of services in this particular context. The cited research's cluster analysis revealed differences between tourists in terms of the weight that each tourist assigns to the factors under study, Thus, their investigations revealed that the factors with the highest average satisfaction in sun and beach tourism were accommodation and restaurants, while factors such as public infrastructure had a lower value than the average.
Studies have also used the chi-square automatic interaction detection [CHAID] algorithm to conduct market segmentation research on tourists. The results indicate that the most significant variable explaining intention to return is the number of previous visits to the destination in question, as a direct relationship exists between the two variables [17,[42][43][44][45]. When tourists visit a destination for the first time, the variable that most influences their intention to return is their weighted overall satisfaction. If the tourists have previously visited the destination, CHAID analysis highlights place of stay as the variable that most strongly affects their intention to return.
Researchers have also compared the SERVQUAL and SERVPERF scales' ability to measure the perceived quality of thermal spas [46]. The cited study was based on surveys of 12 thermal spring facilities in the Autonomous Community of Galicia. Some authors [46] found that most clients of these thermal spas are Spanish and, more specifically, Galician. In addition, the typical clients are individuals over 65 years old and married or with a life partner, without any formal education or with a primary education, as well as having a low income [i.e., an average family unit monthly income of less than 600 euros]. The results clearly confirm quality of service and customer satisfaction's influence on behavioral intentions in the tourism industry.
A market study conducted by Nextel in retail agencies throughout Spain found that thermal spa tourism has an extremely wide range of possible consumers [47]. To commercialize their products correctly, spas must take into account four very specific types of clients. The first segment is young weekend clients without specific needs who seek to enjoy a pleasant, relaxed weekend and who benefit from anti-stress treatments that are more fun than therapeutic. These consumers consider price decisive in their choice of destinations or products [47].
The second segment is clients between 35 and 55 years old with specific needs regarding therapeutic treatments, for whom price is a secondary factor compared to spas' specialization in specific medical conditions or facilities. The third segment is clients 55 years and older, who are also interested in specific establishments in which their ailments can be treated but whose choices are strongly affected by price. These consumers do not pay excessive attention to other factors such as establishments' category and additional services [47].
The last segment is company clients looking for new alternatives to traditional hotels in which to hold meetings involving internal departments, which can be given as gifts to their clients or distributors, and which offer the added value of thermal spa treatments. Price is not a decisive factor overall, although these consumers' selections must fit an existing advertising and/or promotion budget [48]. Basically, this final segment prefers 2-to 4-day stays with extremely simple treatments. Other studies focused especially on thermal spa therapies are in agreement with the cited researchers that this type of tourism mainly attracts older people [49].
The profile of the visitor to Spanish spas is mostly women over 60 years of age and 98% are of national origin. Further, 58% stay more than 8 days, 32% stay between 3 days and a week and respond to an age range of 30 to 45 years. Additionally, 40% of users receive an average of between 3 and 5 treatments, the most demanded being water services (68%), followed by manual techniques (12.7% of cases). The average cost of the client ranges between 251 and 500 euros (National Observatory of Thermalism). Health tourism in Andalusia is characterized by a high presence of national tourists (78.6% in 2016) of which 33.3% are from Andalusia itself, a greater presence of women (54%) than men (46%) [50].

Data Collection Procedures
Scientific research can follow two paths: the quantitative approach and the qualitative approach, both being equally valid given that each of these approaches has its own method. Each one establishes different procedures to generate observations, record them, analyze them, and apply the results according to their scope [51].
In market studies [52], both qualitative techniques (Delphi, expert panels, etc.) and quantitative techniques are commonly used, both descriptive or descriptive-analytical cross-sectional studies (as in this case), as well as essay or experiment type studies, or follow-up studies.
The following steps were taken to ensure an appropriate, adequate sample (see Table 1). The sample was designed based on information gathered from external sources, in particular, a telephone and/or email survey of individuals responsible for managing spas. The relevant data were used to calculate the number of users per year per spa, as well as their distribution by age groups and length of stay. The overall population size was estimated at 53,231 users per year with an age equal to or greater than 15 years old. These individuals were thus defined as this study's primary subjects. Stratified sampling was used, with each stratum coinciding with one of the participating spas. The sample size was estimated using Equation (1): in which n c is the size adjusted to ensure complete sampling (i.e., strata) and Deff is the design's effect or the relationship between the variance under stratified sampling and under simple random sampling. This variance was arbitrarily estimated a priori as 1.5, which was shown to be an overestimate. Any proportion calculated in this study had to meet the requirements of an accuracy of 5% and a 95% confidence level, with the worst possibility considered to be p = 0.5. This implied that n = 384 and n c = 576, which was set as the minimum. The sample size and/or stratum was the result of a proportional fixation (i.e., the largest spas received more than one visit during the survey). The final result was 725 valid questionnaires, which coincides with some authors' suggestion that sample size [53] be based on 10 subjects per analyzed variable. The sample size's adequacy was thus confirmed since 500 subjects would be a quite good sample size and 1000 or more subjects would be an excellent sample size.
The researchers of this study surveyed users in each spa in Andalusia, defined as a center with waters classified as medicinal mining [15] with prior authorization from those responsible. The same researchers asked the questionnaire questions face to face and asked all users on the days they visited each center, and only three users preferred not to participate. The choice of the tourist in the center was not a selection criterion at the origin.

Questionnaire Reliability
A cross-sectional survey design was followed, and the scale validated by [44] in a study assessing the quality of sun and sand tourism destinations on the Costa del Sol was adopted. The final questionnaire included some slight modifications from that used in the cited research. Once the scale was adapted, a diagnostic test was conducted to check test-retest reproducibility with 30 spa users between 15 and 28 years old, who did not belong to the main sample and who had visited 3 spas in the study. The intraclass correlation coefficient or kappa index was used depending on the type of variable (i.e., quantitative or qualitative), with values given based on Landis and Koch's scale [53].

Results
Regarding the values obtained for each dimension, the highest average values were given to items assessing accommodations (6.08 ± 1.15), spas (6.13 ± 1.08), restaurants (5.99 ± 1.29), and public roads and urban and natural environments (5.98 ± 1.22). The general satisfaction dimension stands out with an average of 6.17 ± 1.06, and the last item evaluating general satisfaction with quality and/or quality to price ratio also was assigned a high average value of 5.58 ± 1.29 (see Table 2).
At the population level, the average age of spa tourists in Andalusia is 56 years ± 0.8. Female users are a clear majority (61.8%). In terms of occupation, two large groups stand out among the tourists that visit Andalusian spas: retirees (53.1% ± 2.1) and employed visitors [34% ± 2.0].
Regarding monthly income, the sample includes three main income brackets. These are, first, those who have an average monthly income of less than 500 euros (29.8% ± 2.0); second, those who earn an average of 501 and 1000 euros per month (33.1% ± 2.0); and, third, those who have a monthly average income of between 1001 and 1500 euros (26.6% ± 1.9). Thus, a high percentage of the study population fall inside low income brackets.
Tourists who visit Andalusian spas come mostly from the Autonomous Community of Andalusia (82.9% ± 1.6). Users' provinces include those near the spa resorts such as Malaga (21.9% ± 1.7), Jaen (20.4% ± 1.6), and Granada (18.3% ± 1.6). The sample has no foreign users (see Table 3).  Slightly more than half (51.3% ± 2.2) of spa tourists in Andalusia travel as a couple, while 19.7% ± 1.7 are families. In terms of accommodations, 49.6% ± 1.7 choose three-star hotels. The highest percentage of overnight stays corresponds to those sponsored by the Institute for the Elderly and Social Services (Imserso) (48.4% ± 2.1), followed by weekend and long weekend tourists (28.1% ± 2.0). The average stay is one week. Close to half (46.8% ± 2.2) of these spa users repeat their visit, while nearly all book their trips without intermediaries and only use the spas' services (99.1% ± 0.4) (see Table 4).  (Spas) with SUDAAN. b : The "sd" (standard deviation) has been calculated in the case of the sampling distribution and the "se" (standard error) in the population distribution. c : In the sample, there was no user from abroad or from the Autonomous Communities of Aragon, As-turias, Balearic Islands, Galicia, or Navarra. (Spas) with SUDAAN. b : "sd" (standard deviation) has been calculated in the case of the sampling distribution and the "se" (standard error) in the population distribution. A full 92.4% ± 1.1 of the respondents reported that they intended to recommend the spas, while only 2.0% ± 0.6 had no intention of recommending them. Nearly two-thirds (63.6% ± 2.0) of the tourists surveyed said they intended to return the following year, and just 1.6% ± 0.5 had no intention of ever returning (see Table 5). Significant differences were found between the dimensions for the entire sample (p < 0.001). After comparing the 7 dimensions in pairs, 4 non-heterogeneous groups of dimensions were obtained. These are, from lowest to highest mean value: (1) sports facilities; (2) infrastructure, transportation, and other services; (3) leisure, culture, and shopping; and (4) accommodations, restaurants (i.e., meals), spas, and public roads and/or urban and natural environments (see Table 6). a Percentage distribution (without decimals for simplification) in 6 categories separated by hyphens, ranging from 1 ("Not at all satisfied") to 7 ("Very satisfied"): 1 a < 2; 2 a < 3; 3 a < 4; 4 a < 5; 5 a < 6; and 6 to 7; note that this is a weighted sample; b Arithmetic mean ± standard deviation, weighted. c standard error, ponderado y taking into account stratified sampling using SUDAAN software's DESCRIPT procedure; d 95% CI = confidence interval, calculated as mean ± 1.96 x se; e calculated using SUDAAN software's REGRESS procedure; numerical value of p given; f multiple pairwise comparisons using SUDAAN software's DESCRIPT procedure; means with different letters (w-z) are statistically different (p < 0.05, after Bonferroni correction for multiple comparisons).
Univariate analysis was conducted using simple linear regression models, and each spa facilities' statistical significance was obtained for the spa dimension. The results indicate that the spas under study are extremely heterogeneous (p < 0.05). More specifically, significant differences exist between facilities' visitors for the variables age, occupation, and income. Overall, older people, retirees, and those with incomes of less than or equal to 1000 euros are more satisfied.
In addition, significant results were obtained with reference to family versus other visitors traveling alone, in pairs, and with friends, with tourists who come with their family being less satisfied. Regarding accommodations in 3-and 4-star spa hotels, people staying in more luxurious hotels are more satisfied. The number of overnight stays is also important since satisfaction increases as guests spend more days in the spas (see Table 7). 0.02 ± 0.01 * 0.04 ± 0.01 * 0.02 ± 0.01 * −0.02 ± 0.01 0.08 ± 0.01 * 0.02 ± 0.01 * 0.03 ± 0.01 * 0.04 ± 0.01 * Number of previous visits in last 5 years 0.05 ± 0.02 * 0.01 ± 0.03 0.02 ± 0.02 −0.10 ± 0.04* −0.04 ± 0.03 −0.02 ± 0.02 -0.06 ± 0.03 * 0.00 ± 0.02 a : Models (one for each variable -column 1-and dimension/total satisfaction value) using the REGRESS procedure of SUDAAN (it takes into account stratified and weighted sampling). The coefficients βeta ± standard error are expressed in the Table against the category "ref." (reference category). The "*" indicates a category significantly different (p < 0.05) from the reference category. b : Dimensions. ACC: Accommodation, RES: Restoration, SPA: Spa, SPO: Sports facilities, LEI: Leisure/culture/shopping, PUB: Public roads/urban and natural environment, INF: Transport infrastructures and other services. c : The coefficient of this category for the variable "Accommodation" is meaningless, as they were NS/NC values that were substituted for the mean.
Multivariate analysis was conducted using linear regression models of the spa, sociodemographic, and dimension variables. The global scores overall reveal significant differences for four spas, indicating heterogeneity. Occupation is a significant factor, with retirees more satisfied than those in other occupations, and origin is as well, with tourists from Andalusia being more satisfied than those from other autonomous communities. When the scores were examined by dimensions for the age variable, significant values were found for spas and sports. Namely, older tourists reported greater satisfaction with spas, and younger visitors registered lower satisfaction with sports (see Table 8).  The "*" indicates significantly different category (p < 0.05) from the reference category. Using a "step-by-step backward" procedure based on statistical significance (p < 0.05 to include a variable **). b : Dimensions. ACC: Accommodation, RES: Restoration, SPA: Spa, SPO: Sports facilities, LEI: Leisure/culture/shopping, PUB: Public roads/urban and natural environment, INF: Transport infrastructures and other services. c : The coefficient of this category for the variable "Accommodation" is meaningless, since they were NS/NC values that were substituted by the mean.

Discussion
The spa tourists surveyed were, on average, 56 years old and predominantly female (61.8%). Regarding employment status, pensioners (53.1%) and employed individuals (34.0%) make up a majority of the study population. Similar correlations are detected in other studies [50].
Their average monthly income is below 1000 euros (62.9%), and they primarily come from the region of Andalusia (82.9%). These results confirm other studies' findings [53] that spa tourism is shaped by Spain's social security subsidies and language barriers.
However, internationally, wellness tourists spend an average of 1528 per trip and 53% more than the typical international tourist [8,9,23].
Overall, this profile differs from sun and sand tourism in the summer [41], whose tourists are mostly younger, actively employed individuals from abroad, with higher incomes. However, this description does not fit winter tourism, as 40% of these tourists are retirees older than 65 years. Spa tourism also shows similarities with Mallorca's winter tourism [26] except in terms of tourists' origin, since spas visitors are primarily from within Spain and nearby areas. [50].
Prior research has found evidence that elderly tourists have a positive effect on seasonality, which is also an important aspect of the current demand for spa tourism [44]. About half (51.3%) of spa tourists travel with their life partner and 83.9% choose to stay in hotels. These percentages are slightly lower than those found by studies of Mallorca [25] and sun and beach tourism on the Costa del Sol [44]. Notably, 5.8% of the tourists who visit spa facilities do not stay overnight, which is consistent with these individuals' place of origin-in many cases, quite local. The highest percentage of overnight stays (48.4%) is due to Imserso [54] subsidies that allow an average stay of one week [50].
The satisfaction result variable explains some factors such as age, origin, or occupation. The present results include significant differences associated with the variables of age, occupation, and income. On average, older spa users, retirees, and those with monthly incomes of less than or equal to 1000 euros are more satisfied, which matches. This could be explained by the benefits and properties of hot springs in improving health [1][2][3]6,26] by finding lower levels of satisfaction for younger spa tourists. This may be because it is a tourist profile not focused on therapeutic treatments, as evidenced by research on Polish spa tourists [17].
The current research's results for the study population reveal that 92.4% of the people who visit spas intend to recommend them and that 63.6% express an intention to return the following year that, together with a satisfaction above the average in accommodation, restoration, or Spa, confirms that experiences connected with water correlate positively with satisfaction [7] and correlate with other investigations that have shown that when tourists visit a destination the first time, the most significant variable is the weighted average satisfaction [42,45]. These data show a high degree of loyalty, which indicates the possibility of long-term policies, and the need to continue to maintain quality levels services, facilities, and surroundings, or it is necessary to take care of the image, maintain loyalty, and attract new clients.
The current findings also include a significant relationship (p < 0.001) between future intention variables and the scores given for items in the main questionnaire with regard to greater satisfaction. Various prior studies have found a significant positive connection between perceived quality, satisfaction, and loyalty, confirming that perceived quality affects tourists' satisfaction and loyalty and that this quality is determined by satisfaction with visit variables [41]. These data correlate with the world development of tourism [8] and with the sustainable development of rural destinations [16]. In addition, results indicate that satisfaction significantly affects intention to return. Other studies have additionally shown that satisfaction is positively correlated with repeat purchase of services, loyalty, and intention to recommend to others. The present study also confirmed this association since spa tourists who plan to repeat their visit are the most satisfied with their stay in destinations, and show, as in other studies, repeated purchases of services and products [25,54].
A careful analysis of the results shown in the above tables revealed that satisfaction levels are high in all dimensions, that is, above the scale's central value of four, so the study population includes no dissatisfied users but rather more or less satisfied spa tourists which may mean that the quality of spa destinations in Andalusia is high, as argued by the literature that has studied the relationship between service quality and tourist destinations [18,36,37,39,40].
The factors that show the highest average satisfaction values are 'accommodation', followed by 'spa', in third place is 'restaurants' and in fourth place with very little difference 'public roads/urban and natural environment'. Significant between dimensions, spa tourists have high satisfactions and are similar in accommodation, restaurants, spa, and public roads/urban and natural environment. For the lowest values, they are given in two groups, on the one hand, sports facilities and, on the other hand, 'transport infrastructures and other services'. It follows in this way that satisfaction is higher in the variables that depend directly on the spa. These issues correlate with those obtained in similar studies in Poland [Dryglas] and reinforce the individual motivations that demand these services [19,20,24,55] and companies strive to seek new value propositions to attract new market segments [24].
Other researchers [56] have also confirmed high levels of satisfaction among tourists supported by the National Health Insurance Fund. The items that were given values below the current scale's average level include such variables as soft sports facilities, tennis, other sports facilities, buses, bus lines, and tourism offices. These elements correlate with other studies' conclusions [41] that, in general, tourists are more satisfied with private services than with public services.

Conclusions
Thermal spring waters have been used to address health issues since ancient times, and, currently, thousands of people travel every year to thermal spring facilities in search of healthier vacations and improved health and wellbeing. Thermal spa tourism also has a high potential for attracting varied tourist segments. Spa tourism can further be an important way to diminish seasonality, which is one of the tourism industry's main problems.
The profile of Andalusian spa tourists is predominantly female, over 56 years old, and retired with average incomes below one thousand euros.
The level of satisfaction of its users is above the average in the factors of accommodation, restaurant, spa, and public roads and natural environment and most of these clients would recommend it to others and express the intention to return which indicates Andalusians show a very high level of user loyalty.
Spain's Andalusia has thermal spas that are major tourist attractions, especially as they are located in places of great natural beauty. These facilities are an important part of the sustainable development of inland towns' economies. These research's findings thus have important implications for tourism entrepreneurs, offering new strategies to encourage diversification and providing spa hotels with solutions that attract guests interested in enhancing their health and wellbeing.
In view of the results obtained, conclusions are drawn for management, such as companies could develop strategies to attract users with younger profiles with activities, prices, or spaces more appropriate to their needs and seek the conditions to increase the satisfaction of these segments.
This study presents as a limitation that does not compromise the results the lack of a specific questionnaire for spa tourism, but the authors found it interesting that Andalusia is one of the most important sun and beach destinations in the world, and the questionnaire used was validated for that population [Otero Cordones, 2003]. A few kilometers away, this tourism coexists with the spas, and it seemed a reference that would allow comparison.
This research was carried out in the Andalusia region, which presents a limitation to the potential external validity, but it will be up to the reader to interpret whether or not these results may be applicable to their environment, depending on the possible similarities. In our opinion, it would be extrapolated to any spa in Spain or countries with similar customs and weather conditions.
Our study is a cross-sectional study, which may be a limitation that could be accompanied by qualitative studies in future research lines.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.