The Scenic Beauty of Geosites and Its Relation to Their Scientific Value and Geoscience Knowledge of Tourists: A Case Study from Southeastern Spain

Scenic beauty is one of the most-commonly used indicators in the inventory and assessment of geosites for geoconservation, geoheritage management and geotourism development. It is an important driver of tourists to visit natural areas and it also provides support for the protection of natural heritage. Previous studies on scenic beauty mainly focused on landscape preference and physical characteristics of geosites that affect scenic beauty appreciation. The relationships between the scenic beauty of geosites, their scientific value and the geoscience knowledge of tourists has not been empirically investigated in detail. Hence, this study investigates this relationship using 34 geosites from southeastern Spain. For this purpose, 29 respondents with a geoscience background and who all visited the 34 geosites, 43 respondents with a geoscience background but who did not visit the geosites, and 104 respondents with no geoscience background and who did not visit the geosites, participated in a survey. The first group rated the scenic beauty and the scientific value of the geosites based on a direct field visit during which the scientific background of these geosites was given. On the other hand, the latter two groups rated scenic beauty using representative photos of the geosites. A five-point Likert scale was used to rate the scenic beauty and the scientific value of the geosites. We found a significant relationship between the scenic beauty of geosites and their scientific value, and this relationship becomes more significant if the geoscientific knowledge of the respondents increases. One-way ANOVA results indicated that a geoscience background contributed to higher perceived scenic beauty, especially for those geosites that in general were considered as more scenic by all the respondent groups. It was also found that geosites with viewpoints received in general higher scenic beauty and scientific value ratings.


Introduction
Geodiversity, which has been recognized as a concept worth investigating from the 1990s onwards, is in recent times more frequently introduced into scientific nomenclature. It has aroused a strong interest of researchers from geology, geography, biology, spatial planning, general tourism as well as national geotourism and cultural heritage [1]. Geodiversity refers to the "natural range (diversity) of geological (rocks, minerals, fossils), geomorphological (landforms, topography, physical processes), soil (pedological) and hydrological features. It includes their assemblages, structures, systems and contributions to landscapes" ( [2], p. 14). Geodiversity is the "abiotic equivalent" or "natural twin" of  Table A1 in Appendix A; for their scenic beauty ratings, see Table  A4 in Appendix A; for their scientific value rating, see Table A5 in Appendix A; for their photos, see Questionnaire S1 in the Supplementary Materials.
The study area offers a unique opportunity for teaching field geology and geomorphology [70], and many student field trips are organized by European universities to the region [71,72].
Geomorphological and Geological Setting of the Study Area The study area is amongst the driest regions of Europe [73] which makes it an ideal  Table A1 in Appendix A; for their scenic beauty ratings, see Table A4 in Appendix A; for their scientific value rating, see Table A5 in Appendix A; for their photos, see Questionnaire S1 in the Supplementary Materials. The study area offers a unique opportunity for teaching field geology and geomorphology [70], and many student field trips are organized by European universities to the region [71,72].

Geomorphological and Geological Setting of the Study Area
The study area is amongst the driest regions of Europe [73] which makes it an ideal place to learn and enjoy about the geology and geomorphology as the vegetation cover is rather limited allowing many geomorphological and geological features to be observed easily [71]. The region is characterized by a series of mountain ranges or Sierras (Betic chain, resulting from the Alpine orogeny and mainly consisting of hard Paleozoic and Mesozoic rocks) and uplifted Cenozoic basins dominated by unconsolidated sediments [74][75][76][77]. From a lithological point of view, volcanic, sedimentary and metamorphic rock types outcrop in the area [78][79][80][81][82], resulting in diverse landform types. The region is also tectonically active [74,75,[82][83][84] which allows to observe active faults, horst and graben landforms as well as volcanic features [77,[83][84][85][86]. Erosional and depositional features include various types of mass movements, gullies, badlands, fans and landforms resulting from fluvial and coastal dynamics [87][88][89]. Several mineral deposits (e.g., lead, iron, gold amongst others) have been mined from prehistory into modern times [78,90,91], of which several traces are still preserved in the landscape. Next to mining, other traces of past human-environment interactions on the landscape can be observed such as old farms and agricultural terraces being abandoned in the mountains [92,93], intensive (greenhouse) agriculture in the coastal and alluvial plains [94][95][96][97] and tourism development in the coastal areas (littoralization) [98] with its large impacts on groundwater hydrology [96,[98][99][100].

Data Collection and Analysis
For this study, 34 geosites were selected from southeastern Spain (see Figure 1b for their location; Table A1 in Appendix A for their description, Table A4 in Appendix A for their scenic beauty rating; Questionnaire S1 in the Supplementary Materials for their photos), which were part of an educational 'Physical Geography' field excursion, held from 18-24 May 2019. These geosites were selected based on their scientific value in order to present a variety of topics related to the geomorphology, geology, pedology, hydrology, and archaeology of the region, and with a strong focus on human-environment interactions in a Mediterranean environment at various timescales. Most geosites (30) are related to geomorphology. In addition, 27 of them have viewpoints while 16 geosites can be directlylinked to human-environment interactions. Specifically, the geosites include volcanic cones and columns, horst and graben structures, faults, travertine dams, sand dunes, tafoni, alluvial fans, landslides, gullies, badlands, mining, and archaeological settlements. Whilst scenic beauty did not form a major criterion when selecting the sites visited during the field trip, several geosites have been selected as they provide spectacular views to surrounding landscapes and are, hence, ideal sites for a physical geography field trip.
A total of 176 respondents (actual and potential tourists) participated in this survey (see Table A3 in Appendix A for their socio-demographic background). Of these, 104 were persons with no geoscience background and who did not visit the Spanish geosites (hereafter called NGB-NV); 43 were persons with a geoscience background and who did not visit the selected geosites (hereafter called GB-NV); and 29 were persons with a geoscience background and who all visited the selected geosites during the educational 'physical geography' field excursion, held in 2019 (hereafter called GB-V). The NGB-NV group consisted of persons whose educational background is unrelated to geosciences. On the other hand, the GB-NV and GB-V group comprised persons who studied geography and/or geology, and whose education level was bachelor's degree and higher. The GB-NV group was purposefully included to control for the effects of direct experience to the geosites and expert information on the scenic beauty rating between the NGB-NV and GB-V groups. The NGB-NV and GB-NV groups rated scenic beauty in an online survey based on representative photos of the geosites, while the GB-V group made the assessment after visiting the geosites in 2019. The scenic beauty rating by these three groups was used to determine the contribution of geoscience knowledge for scenic beauty.
For the photo-based assessment, a total of 74 photos representing the 34 geosites were selected. These photos were presented in the order in which the GB-V group visited them in the field. Each geosite was represented with two to three photos to provide typical views of the geosites to the respondents. The first photo of each geosite showed an overview of the surrounding landscape, and the second (and third, if any) photo/s usually showed the geosite in more detail. From the photos used in the survey, 63 were photos taken during the field trip while the remaining 11 photos were taken from images in Google Earth and the internet due to the poor quality of some of the photos we had. While selecting the latter photos, care was taken to make them representative of what the GB-V group saw on-site and to illustrate in the best possible way the main geo-feature of the selected site.
The online survey had two sections. The first section consisted of photos and a fivepoint Likert scale (where 1 = not at all interesting, 2 = slightly interesting, 3 = moderately interesting, 4 = very interesting, and 5 = extremely interesting) for rating scenic beauty, while the second section comprised the socio-demographic background of the respondents (see Questionnaire S1, in the Supplementary Material).
The online survey was conducted from 18 February to 10 March 2021, and all authors of this study sent the online survey via email to people in their network. In some cases, people who were first contacted by the authors further distributed the survey to other people, and hence it was rather difficult to know the total respondents contacted. However, by counting those for which we had reliable data, it was estimated that the survey was distributed to over 550 people. A total of 154 completed surveys were received, of which, 7 were discarded because the respondents indicated that they had previously visited one or more of the geosites. Hence, the 147 completed responses (104 NGB-NV and 43 GB-NV) were used for further analysis of the scenic beauty of the geosites. The NGB-NV group comprised persons with educational background from ca. 25 disciplines such as archaeology, agriculture, biology, chemistry, engineering, economics, history, languages, management, medicine, psychology, sustainable development and tourism (see Figure A1 in Appendix A) while the GB-NV group consisted of 32 geographers and 11 geologists.
On the other hand, the GB-V group assessed the scenic beauty and scientific value of geosites during an educational field excursion to the Spanish study area. This group comprised a total of 32 students (from KU Leuven and Free University Brussels) who enrolled in the 3rd bachelor and 1st master in Geography. It also included two KU Leuven professors of physical geography (>20 years of experience in the region), who led the trip, and three field assistants with a master's degree and educational background related to geoscience. Hence, the scenic beauty and scientific value assessment questionnaire was distributed to the 37 participants of the field trip. A total of 29 persons completed the questionnaire, including the two professors and the three field assistants.
The GB-V group had learned about the geosites in southeastern Spain before and during the field excursion, which enabled them to evaluate the scientific value of each geosite. They collected and read research papers about the geosites and they also made short presentations in the classroom before the field excursion. In addition, they were also given on-site scientific information about the characteristics, genesis, importance for earth sciences as well as for human-environment interactions of each geosite by the two professors ( Figure 2). Furthermore, shortly after visiting all sites, the GB-V group also evaluated the scenic beauty of the 34 geosites.
The GB-V group was asked to rate the scenic beauty and scientific value of the geosites using a five-point Likert scale (where 1 = not at all interesting, 2 = slightly interesting, 3 = moderately interesting, 4 = very interesting, and 5 = extremely interesting), which is a similar scale as that provided to the other two groups of respondents (see Questionnaire S2 and Questionnaire S3 in the Supplementary Materials). In addition, the GB-V group was also asked to list interesting features of each geosite related to its scenic beauty and scientific value. The socio-demographic profile of the respondents was also collected. The GB-V group was asked to rate the scenic beauty and scientific value of the geosites using a five-point Likert scale (where 1 = not at all interesting, 2 = slightly interesting, 3 = moderately interesting, 4 = very interesting, and 5 = extremely interesting), which is a similar scale as that provided to the other two groups of respondents (see Questionnaire S2 and Questionnaire S3 in the Supplementary Materials). In addition, the GB-V group was also asked to list interesting features of each geosite related to its scenic beauty and scientific value. The socio-demographic profile of the respondents was also collected.
The GB-V group was briefed about the contents of the survey questionnaire in a classroom one day before departure to Spain. The questionnaire was distributed to them immediately before the start of the field excursion and collected at the end of the field excursion; this helped them to be familiar with the questions.
Correlation analysis and scatter plots and boxplots were made to test the relationship between the scenic beauty and scientific value of the geosites. In addition, to further investigate the relationship between scenic beauty and scientific value, a word cloud analysis was conducted. This revealed the most frequent keywords that the respondents reported in order to describe the interesting features of geosites reflecting their scientific value and scenic beauty. In the word cloud analysis, only keywords were selected, and co-occurring words were removed before the analysis. For example, if the respondent mentioned "view over sierra and sea", then the keyword 'view' was taken, and the words 'over', 'sierra' and 'sea' were dropped from the analysis. This helped to avoid unnecessary repetition of words.
In order to investigate the relationship between scenic beauty and geoscience knowledge, one-way ANOVA and post-hoc multiple comparisons were conducted for the mean scenic beauty ratings of the geosites among the three respondent groups. Additionally, to identify why some geosites received higher scenic beauty and scientific value ratings by the respondent groups, the 34 geosites were grouped into five clusters based on the presence of particular features of interest to tourists at each geosite. The features used as criteria were local geo-features (such as volcanic cones and columns, horst and graben structures, faults, travertine dams, sand dunes, tafoni, alluvial fans, landslides, gullies, The GB-V group was briefed about the contents of the survey questionnaire in a classroom one day before departure to Spain. The questionnaire was distributed to them immediately before the start of the field excursion and collected at the end of the field excursion; this helped them to be familiar with the questions.
Correlation analysis and scatter plots and boxplots were made to test the relationship between the scenic beauty and scientific value of the geosites. In addition, to further investigate the relationship between scenic beauty and scientific value, a word cloud analysis was conducted. This revealed the most frequent keywords that the respondents reported in order to describe the interesting features of geosites reflecting their scientific value and scenic beauty. In the word cloud analysis, only keywords were selected, and co-occurring words were removed before the analysis. For example, if the respondent mentioned "view over sierra and sea", then the keyword 'view' was taken, and the words 'over', 'sierra' and 'sea' were dropped from the analysis. This helped to avoid unnecessary repetition of words.
In order to investigate the relationship between scenic beauty and geoscience knowledge, one-way ANOVA and post-hoc multiple comparisons were conducted for the mean scenic beauty ratings of the geosites among the three respondent groups. Additionally, to identify why some geosites received higher scenic beauty and scientific value ratings by the respondent groups, the 34 geosites were grouped into five clusters based on the presence of particular features of interest to tourists at each geosite. The features used as criteria were local geo-features (such as volcanic cones and columns, horst and graben structures, faults, travertine dams, sand dunes, tafoni, alluvial fans, landslides, gullies, etc.), humanenvironment interaction (such as archaeological sites, agricultural fields-both currently in use and abandoned ones, dams and reservoirs) and viewpoints. The resulting clusters of geosites, in order of their numbers of geosites they have, are (1) HE = human-environment interaction features (2 geosites); (2) HE-LG-VP = human-environment interaction feature, local geo-feature and viewpoint (4 geosites); (3) LG = local geo-feature (5 geosites); (4) HE-VP = human-environment interaction feature and viewpoint (10 geosites); and Land 2021, 10, 460 8 of 27 (5) LG-VP = local geo-feature and viewpoint (13 geosites) (see Table A2 in Appendix A for the list of clustered geosites). Table 1 shows the mean scenic beauty ratings of the 34 geosites as a function of sociodemographic factors. It can be seen that, on average, women (mean = 3.35; sd = 0.51) rated scenic beauty relatively higher than men (mean = 3.25; sd = 0.52). Young people with age 18-29 (mean = 3.31; sd = 0.57) rated scenic beauty higher than other age groups. In terms of education, those with bachelor's degree (mean = 3.34; sd = 0.53) rated scenic beauty higher than those with master's and PhD degrees. Besides, those who live in Belgium (mean = 3.46; sd = 0.6) rated scenic beauty higher than those from other countries. Furthermore, those who did not travel outside their continent (mean = 3.42; sd = 0.47) rated the scenic beauty of the geosites higher than respondents who visited one or more other continents. It was also found that respondents with geoscience background and who visited the geosites (mean = 3.40; sd = 0.34) rated scenic beauty higher than those with no geoscience background and who did not visit the geosites. As to the respondents' general primary preferred types of attractions, those who prefer attractions related to biodiversity and geology-landscape (mean = 3.6; sd = 0.56) rated scenic beauty higher than other groups. Table 1. Mean scenic beauty ratings (scale range 1-5; standard deviations in parentheses) of the 34 geosites as a function of socio-demographic factors of the respondents (n is the number of persons). For information on the socio-demographic background of the respondents, see Table A3 in Appendix A. Does not include data of those who visited the geosites. b Includes temporary (e.g., students) and permanent residence. c Respondent group refers to respondents grouped following their geoscience background and field visits: NGB-NV = no geoscience background and no visit to geosites; GB-NV = geoscience background and no visit to geosites; GB-V = geoscience background and visit to geosites. d Significant at p < 0.1; scenic beauty rated on 1-5 scale.

Socio-Demographic Factor
However, one-way ANOVA results showed that gender, age, country of residence and number of continents visited did not significantly influence mean scenic beauty rating while the respondent group and primarily preferred attraction were found significant at p < 0.1 (Table 1, Figure 3). Post-hoc pairwise (multiple) comparisons showed that none of the socio-demographic factors, including the respondent group and primarily preferred attraction, revealed a significant difference (at p < 0.1) in the scenic beauty ratings of the 34 geosites.
Culture-history (n = 37) 3.09 (0.56) Culture-history and geology-landscape (n = 27) 3.37 (0.43) Geology-landscape (n = 51) 3.31 (0.45) 2.05 0.06 d a Does not include data of those who visited the geosites. b Includes temporary (e.g., students) and permanent residence. c Respondent group refers to respondents grouped following their geoscience background and field visits: NGB-NV = no geoscience background and no visit to geosites; GB-NV = geoscience background and no visit to geosites; GB-V = geoscience background and visit to geosites. d Significant at p < 0.1; scenic beauty rated on 1-5 scale.
To further investigate the effect of socio-demographic factors (excluding age as there was no data for the GB-V group in the "respondent group" factor) on scenic beauty rating, two-way ANOVA was conducted. Only the interaction of education level and country of residence were found significant (F value = 3.56, p < 0.05). However, further pairwise multiple comparisons indicated that none of the interaction effects were significant at p < 0.1. To further investigate the effect of socio-demographic factors (excluding age as there was no data for the GB-V group in the "respondent group" factor) on scenic beauty rating, two-way ANOVA was conducted. Only the interaction of education level and country of residence were found significant (F value = 3.56, p < 0.05). However, further pairwise multiple comparisons indicated that none of the interaction effects were significant at p < 0.1.

Correlation Analysis Results
Geosites with a higher scientific value were rated higher for their scenic beauty, and those with lower scientific value were also rated lower by all the respondent groups (Figures 4a-c and 5a-c). In addition, the mean scenic beauty rating by all the respondent groups is positively related to the scientific value of geosites ( Figure 5d).

Correlation Analysis Results
Geosites with a higher scientific value were rated higher for their scenic beauty, and those with lower scientific value were also rated lower by all the respondent groups (Figures 4a-c and 5a-c). In addition, the mean scenic beauty rating by all the respondent groups is positively related to the scientific value of geosites ( Figure 5d).

Figure 4(a-c).
Relationship between scenic beauty ratings of the 34 geosites (rated by the three groups of respondents, grouped following their geoscience background and field visits: NGB-NV = no geoscience background and no visit to geosites; GB-NV = geoscience background and no visit to geosites; GB-V = geoscience background and visit to geosites) and scientific value (rated by the GB-V group). Both scenic beauty and scientific value were rated on a 1-5 scale. r = Pearson's correlation coefficient between scenic beauty (as rated by each respondent group) and scientific value (rated by the GB-V group). See Figure 1b for the location of geosites (indicated by their number in the figure); Table A1 in Appendix A for their description, Table A4 in Appendix A for their mean scenic beauty ratings; Table A5 in Appendix A for their scientific value ratings; Questionnaire S1 in the Supplementary Materials for their photos. (a-c) Relationship between scenic beauty ratings of the 34 geosites (rated by the three groups of respondents, grouped following their geoscience background and field visits: NGB-NV = no geoscience background and no visit to geosites; GB-NV = geoscience background and no visit to geosites; GB-V = geoscience background and visit to geosites) and scientific value (rated by the GB-V group). Both scenic beauty and scientific value were rated on a 1-5 scale. r = Pearson's correlation coefficient between scenic beauty (as rated by each respondent group) and scientific value (rated by the GB-V group). See Figure 1b for the location of geosites (indicated by their number in the figure); Table A1 in Appendix A for their description, Table A4 in Appendix A for their mean scenic beauty ratings; Table A5 in Appendix A for their scientific value ratings; Questionnaire S1 in the Supplementary Materials for their photos.
Correlation analysis also indicated that there is a significant (at least at p < 0.05) positive relationship between the scenic beauty and scientific value of the geosites (Figures 4a-c and 5a-c). However, the strength of the relationship depends on the type of respondent group who rated the scenic beauty. A weak relationship (Figure 4a) was found between scenic beauty rating by the NGB-NV group and scientific value (rated by the GB-V group) while a moderate relationship (Figure 4b) was found between scenic beauty as rated by the GB-NV group and scientific value (as rated by the GB-V group). The relationship between the scenic beauty and scientific value of the geosites, both rated by the GB-V group, was strong ( Figure 4c). Furthermore, a positive relationship (r = 0.70) was found between the mean scenic beauty ratings of the geosites by all the three respondent groups and the scientific value ratings of the geosites. This scientific value rating of geosites was then grouped based on mean scientific value rating: Very high = 4.0 and above (6 geosites); High = 3.5-3.9 (11 geosites); Moderate = 3.0-3.5 (12 geosites); Low = below 3.0 (5 geosites). The scientific value rating was done by the GB-V group, 29 persons (see Table A5 in Appendix A). "Scenic beauty rating by all respondents" is the mean scenic beauty rating of the three respondent groups.
Correlation analysis also indicated that there is a significant (at least at p < 0.05) positive relationship between the scenic beauty and scientific value of the geosites (Figures 4a-c and 5a-c). However, the strength of the relationship depends on the type of respondent group who rated the scenic beauty. A weak relationship (Figure 4a) was found between scenic beauty rating by the NGB-NV group and scientific value (rated by the GB-V group) while a moderate relationship (Figure 4b) was found between scenic beauty as rated by the GB-NV group and scientific value (as rated by the GB-V group). The relationship between the scenic beauty and scientific value of the geosites, both rated by the GB-V group, was strong ( Figure 4c). Furthermore, a positive relationship (r = 0.70) was found between the mean scenic beauty ratings of the geosites by all the three respondent groups and the scientific value ratings of the geosites.

Interesting Features Explaining Scenic Beauty and Scientific Value Ratings of the Geosites
The GB-V group was asked to list the most interesting geo-features that make up the scenic beauty as well as the scientific value of the geosites ( Figure 6). The top five most This scientific value rating of geosites was then grouped based on mean scientific value rating: Very high = 4.0 and above (6 geosites); High = 3.5-3.9 (11 geosites); Moderate = 3.0-3.5 (12 geosites); Low = below 3.0 (5 geosites). The scientific value rating was done by the GB-V group, 29 persons (see Table A5 in Appendix A). "Scenic beauty rating by all respondents" is the mean scenic beauty rating of the three respondent groups.

Interesting Features Explaining Scenic Beauty and Scientific Value Ratings of the Geosites
The GB-V group was asked to list the most interesting geo-features that make up the scenic beauty as well as the scientific value of the geosites ( Figure 6). The top five most frequent reported words describing interesting scenic features (Figure 6a) were view = 146, landscape = 38, gully = 34, sea = 32, and travertine = 31.
The top five most frequent words which the respondents mentioned as interesting features of geosites contributing to their scientific value ( Figure 6b) were gully = 33, travertine = 31, terrace = 30, dam = 29, and badland = 21.
Among the most frequent words in the word cloud analysis, 22 are common to both scenic beauty and scientific value (which account for 45.8% of the words in Figure 6a and 57.9% in Figure 6b). These include, in alphabetical order, archaeology, atoll, badland, biodiversity, caldera, castle, dam, dune, flood, geology, gully, history, landslide, mining, rambla, rock, tafoni, terrace, travertine, valley, viewpoint and wind. These common words indicate that the GB-V respondent group appreciates scenic beauty and scientific value on many similar features, supporting our hypothesis that these two geosite values are interrelated.  Among the most frequent words in the word cloud analysis, 22 are common to both scenic beauty and scientific value (which account for 45.8% of the words in Figure 6a and 57.9% in Figure 6b). These include, in alphabetical order, archaeology, atoll, badland, biodiversity, caldera, castle, dam, dune, flood, geology, gully, history, landslide, mining, rambla, rock, tafoni, terrace, travertine, valley, viewpoint and wind. These common words indicate that the GB-V respondent group appreciates scenic beauty and scientific value on many similar features, supporting our hypothesis that these two geosite values are interrelated.

Contribution of Geoscience Knowledge to Perceived Scenic Beauty
As shown in Figure 7, there is a relative agreement among the three respondent groups in the scenic beauty rating. More specifically, from the three pairwise comparisons, there is a better agreement in scenic beauty ratings of the geosites by the NGB-NV and GB-NV groups as their ratings are close to the 1:1 line.
However, in absolute terms, the groups with geoscience background (GB-NV and GB-V) gave higher scenic beauty rating to more geosites (i.e., to 24 and 20 geosites, respectively) as compared to the NGB-NV group (Figure 7). On the other hand, the GB-NV group gave a higher scenic beauty rating to 19 geosites and vice versa for the remaining 15 geosites (Figure 7b). (a) scenic beauty (total words = 48, minimum and maximum word frequency = 5 and 146, respectively); (b) scientific value (total words = 38, minimum and maximum word frequency = 5 and 33, respectively). Note that the word size corresponds to the frequency of that word in its category, i.e., in scenic beauty or scientific value.

Contribution of Geoscience Knowledge to Perceived Scenic Beauty
As shown in Figure 7, there is a relative agreement among the three respondent groups in the scenic beauty rating. More specifically, from the three pairwise comparisons, there is a better agreement in scenic beauty ratings of the geosites by the NGB-NV and GB-NV groups as their ratings are close to the 1:1 line.
However, in absolute terms, the groups with geoscience background (GB-NV and GB-V) gave higher scenic beauty rating to more geosites (i.e., to 24 and 20 geosites, respectively) as compared to the NGB-NV group (Figure 7). On the other hand, the GB-NV group gave a higher scenic beauty rating to 19 geosites and vice versa for the remaining 15 geosites (Figure 7b).
The data also revealed that geosite clusters where viewpoints are present were among those that were rated higher for their scenic beauty as well as for their scientific value ( Figure 8). More specifically, all respondent groups rated scenic beauty higher for geosites that have both a local geo-feature and offer a viewpoint (Figure 8a-c). In addition, the scientific value was also rated higher by the GB-V group for this cluster of geosites (Figure 8d).
One-way ANOVA was used to compare the mean scenic beauty ratings of the 34 individual geosites by the three respondent groups (Figure 9). Significant differences were found for 22 individual geosites at different significance levels: p < 0.05, p < 0.01 and p < 0.001 (see Table A4 in Appendix A).
Post-hoc pairwise multiple comparisons showed significant mean scenic beauty rating differences between the three respondent groups ( Table 2). A significant difference in scenic beauty rating was found between GB-NV and NGB-NV groups for 7 out of 34 geosites. For all these geosites, the GB-NV group rated scenic beauty higher than the NGB-NV group. In addition, a significant difference in scenic beauty rating was found between NGB-NV and GB-V for 17 out of 34 geosites. For 11 of these geosites, the GB-V group rated scenic beauty higher than the NGB-NV group, and vice versa for the remaining six geosites. Furthermore, a significant difference in mean scenic beauty rating was found between GB-NV and GB-V groups for 14 out of 34 geosites. For eight of these geosites, the GB-V group rated scenic beauty higher than the GB-NV group and vice versa for the rest six geosites.  Table A1 in Appendix A for their description, Table A4 in Appendix A for their mean scenic beauty rating; Questionnaire S1 in the Supplementary Materials for their photos.
The data also revealed that geosite clusters where viewpoints are present were among those that were rated higher for their scenic beauty as well as for their scientific value (Figure 8). More specifically, all respondent groups rated scenic beauty higher for geosites that have both a local geo-feature and offer a viewpoint (Figure 8a-c). In addition, the scientific value was also rated higher by the GB-V group for this cluster of geosites (Figure 8d).  Table A1 in Appendix A for their description, Table A4 in Appendix A for their mean scenic beauty rating; Questionnaire S1 in the Supplementary Materials for their photos. LG-VP = local geo-feature and viewpoint, 13 geosites; HE-VP = human-environment interaction feature and viewpoint, 10 geosites; HE-LG-VP = human-environment interaction feature, local geo-feature and viewpoint, 4 geosites; LG = local geo-feature, 5 geosites; HE = human-environment interaction feature, 2 geosites). See Table A2 in Appendix A for the list of clustered geosites. The scientific value rating was done by the GB-V group (29 persons). Both scenic beauty and scientific value were rated on a 1-5 scale (see Tables A4 and A5 in Appendix A for their mean ratings).
One-way ANOVA was used to compare the mean scenic beauty ratings of the 34 individual geosites by the three respondent groups (Figure 9). Significant differences were found for 22 individual geosites at different significance levels: p < 0.05, p < 0.01 and p < 0.001 (see Table A4 in Appendix A).
Post-hoc pairwise multiple comparisons showed significant mean scenic beauty rating differences between the three respondent groups ( Table 2). A significant difference in scenic beauty rating was found between GB-NV and NGB-NV groups for 7 out of 34 geosites. For all these geosites, the GB-NV group rated scenic beauty higher than the NGB-NV group. In addition, a significant difference in scenic beauty rating was found between NGB-NV and GB-V for 17 out of 34 geosites. For 11 of these geosites, the GB-V group rated scenic beauty higher than the NGB-NV group, and vice versa for the remaining six geosites. Furthermore, a significant difference in mean scenic beauty rating was found between GB-NV and GB-V groups for 14 out of 34 geosites. For eight of these geosites, the GB-V group rated scenic beauty higher than the GB-NV group and vice versa for the rest six geosites. LG-VP = local geo-feature and viewpoint, 13 geosites; HE-VP = human-environment interaction feature and viewpoint, 10 geosites; HE-LG-VP = human-environment interaction feature, local geo-feature and viewpoint, 4 geosites; LG = local geo-feature, 5 geosites; HE = human-environment interaction feature, 2 geosites). See Table A2 in Appendix A for the list of clustered geosites. The scientific value rating was done by the GB-V group (29 persons). Both scenic beauty and scientific value were rated on a 1-5 scale (see Tables A4 and A5 in Appendix A for their mean ratings). Table 2. Pairwise multiple comparisons of mean scenic beauty ratings (standard deviations in parentheses) of the 22 individual geosites (where a significant difference in mean scenic beauty rating was found) by the three groups, grouped following their geoscience background and field visits: NGB-NV = no geoscience background and no visit to geosites; GB-NV = geoscience background and no visit to geosites; GB-V = geoscience background and visit to geosites. Scenic beauty rated on a 1-5 scale.

Geosite Number
Geosite Name   Table A1 in Appendix A for their description, Table A4 in Appendix A for their scenic beauty ratings; Questionnaire S1 in the Supplementary Materials for their photos.

Discussion and Conclusions
The objective of this research was to better understand the relationships between scenic beauty of geosites, their scientific value as well as geoscience knowledge of actual and potential tourists. The scenic beauty of 34 geosites in southeastern Spain was evaluated by 176 respondents, which were divided into three groups based on their geoscience background and visits to the geosites.
Socio-demographic factors such as gender, age, education level, country of residence and number of continents visited did not significantly affect the scenic beauty ratings of  Table A1 in Appendix A for their description, Table A4 in Appendix A for their scenic beauty ratings; Questionnaire S1 in the Supplementary Materials for their photos.

Discussion and Conclusions
The objective of this research was to better understand the relationships between scenic beauty of geosites, their scientific value as well as geoscience knowledge of actual and potential tourists. The scenic beauty of 34 geosites in southeastern Spain was evaluated by 176 respondents, which were divided into three groups based on their geoscience background and visits to the geosites.
Socio-demographic factors such as gender, age, education level, country of residence and number of continents visited did not significantly affect the scenic beauty ratings of the geosites by the respondents. Previous studies on the assessment of scenic beauty of landscapes [43,55,62,63] also found that age, gender and education did not significantly affect the scenic beauty ratings of landscapes. In addition, Frank et al. [42] also found that age, gender and personal qualification (i.e., layman, stakeholder and expert) did not significantly influence landscape scenic beauty assessment results. On the other hand, Skřivanová et al. [61] found that there is a significant difference in the scenic beauty rating of landscapes between women and men, the former generally rating it higher. Svobodova et al. [50] studied the visual preferences for physical attributes of mining and post-mining landscapes with 1050 persons in a web-based survey in the Czech Republic and found that visual preferences for landscapes significantly varied based on gender, education level and professional field or study focus. Women rated the scenic beauty of landscapes significantly higher than men. Those with education lower than university-level rated scenic beauty of landscapes significantly higher than those with a university degree. Moreover, respondents whose profession is not related to landscape management rated scenic beauty of landscapes higher than respondents whose education is related to landscape management (e.g., ecology, nature conservation, architecture). López-Martínez [55] pointed out that different findings among studies about the effect of socio-demographic factors on scenic beauty assessments of landscapes might be attributed to differences in study areas, and thereby landscapes being evaluated.
The correlation analysis in this study revealed that there is a positive relationship between scenic beauty and scientific value of geosites, and this relationship improved with an increase in geoscience knowledge of the tourists and with a field visit of these sites. The positive correlation between scenic beauty and scientific value was stronger for the GB-V group than the GB-NV group, and a possible source of difference would be that the former rated the scientific value, received scientific information on-site and visited the geosites. This positive correlation was stronger for both the GB-V and GB-NV groups in comparison with the NGB-NV group, which can be attributed to the geoscience education of both GB groups.
The word cloud analysis of keywords provided by the respondents to describe the most interesting features at each geosite for their scenic beauty and scientific value supports the finding that scenic beauty and scientific value are interrelated. It should, however, be noted in the word cloud analysis that a larger frequency of keywords does not necessarily imply that more respondents reported a particular word (feature); it might well be that few respondents repeatedly mentioned such a word. Moreover, a respondent might have been interested in more than one type of interesting feature for a given geosite, and hence the frequency of the words might be larger.
There is a general consensus among the respondents in their scenic beauty rating in that for most cases, geosites rated higher by one group were also rated higher by the other (see Figure 7). Kalivoda et al. [38] studied scenic beauty rating between experts and non-experts and found that the higher the scenic beauty ranking, the better the consensus between the two groups. They argued that such a consensus plays an important role for the legal protection of geosites.
Overall, the GB-NV and GB-V groups rated the scenic beauty of some geosites significantly higher than the NGB-NV group (Table 2), and this could be because their geoscience knowledge helped them to better appreciate geosites than the NGB-NV group. Reynard and Giusti [101] (p. 152) support this arguing that "perhaps the [scenic] beauty resides . . . less in the outburst of emotions than in the elements of understanding". In addition, the field visit to the geosites by the GB-V group could also have helped them to rate scenic beauty compared to the NGB-NV group.
The GB-V group rated some geosites significantly higher than the GB-NV group, which could be attributed to the expert information the former received while visiting the geosites. Obviously, the GB-V group has better geoscience knowledge about the geosites in this study than the GB-NV group due to the desk research the former had made about the geosites and the expert information it had received on-site. In addition, the field visit could also have had an impact on the scenic beauty rating, allowing the GB-V group to rate scenic beauty higher than the GB-NV group.
In addition to geoscience knowledge, the field visit and the relatively good weather conditions (dry and sunny) during the visit could also explain why the GB-V group rated some geosites significantly higher than the NGB-NV and GB-NV groups. The field visit allowed to see more detailed features of the geosites than the photos. Hull and Stewart [102] argued that photographs may not always capture all the landscape differences, allowing those who visit the sites to rate scenic beauty higher than the group who does not visit.
A good example where the on-site expert information and the field visit might have caused a significantly higher scenic beauty rating of the geosites by the GB-V group than the other two groups are the Gorafe and Alicun de las Torres geosites (geosite numbers 29 and 30, Table 2). Among the interesting features of Gorafe are the view over the Rio Gor valley cut into the Guadix-Baza sedimentary basin and the nearby Bronze-age megalithic park. The latter could not be clearly understood unless one travels to this geosite and also receives expert information about its archaeological history and significance. In addition, during the field excursion, the unique travertine features of Alicun de las Torres were shown and explained to the GB-V group and they were also able to appreciate the flowing spring water and landscape on-site.
On the other hand, the NGB-NV and GB-NV groups rated some geosites significantly higher than the GB-V group (for example, geosite numbers 1, 5, 22 and 28). These were among the lowest-rated geosites in terms of their scientific value by the GB-V group. As scientific value has a relationship with scenic beauty, it might also have contributed to a lower appreciation of these geosites. In other words, the on-site scientific interpretations provided to the GB-V group might have influenced them not to rate these geosites higher for their scenic beauty. For example, at Cerro Pistolas and El Puntal (geosites 5 and 11, respectively), the GB-V group received scientific information about the impacts of the greenhouses on the landscape and the related landscape degradation processes (land levelling, groundwater extraction and pollution by pesticides and degraded plastics), and this might have negatively affected their scenic beauty rating. There was also litter dumped at El Puntal which could have influenced the scenic beauty rating by the GB-V group. Stecker [67] argues that though scientific knowledge could enhance scenic beauty appreciation, it could also work the other way round, i.e., "knowledge prevents appreciation from being malfounded, from appreciating a part of nature for properties it does not have" (p. 400).
The findings in this study that scenic beauty and scientific value of geosites are interrelated and that geoscience knowledge contributes to higher perceived scenic beauty have important implications for geoconservation, geoheritage management and geotourism development. Erikstad [103] indicated that the necessity of geoconservation is not well developed and accepted in many countries around the world. The relationship between the two values could provide a strong support for the protection of the geosites. In addition, more geoscience education and geo-interpretation to people could help them to appreciate the scenic beauty of geosites. Beck et al. [104] indicated that one of the main objectives of interpretation is creating appreciation and deeper understanding of nature. It has also been argued that applying the concept of geotourism and geosites to particular landforms is the best way to transfer geoscience knowledge to society [17]. Thus, the geo-interpretations offered could help enhance visitor experience, thereby contributing to sustainable geotourism development. According to Newsome and Dowling [105] (p. 6), "visitors will always rate their experiences higher if they have also learned something about the landscape and geology they are visiting". People who appreciate the scenic beauty of geosites could in turn play their part for the conservation and management of these geosites. Education about the geosites also raises awareness for their protection [106]. The fact that geosites combining certain features are more interesting to all types of respondents, irrespective of their geoscience background, is important for geotourism development. The presence of a viewpoint was an important factor in the scenic beauty rating. Mikhailenko and Ruban [107] also indicated that the value of viewpoints (also called viewpoint geosites) in the western Caucasus (Russia), is strongly linked to their aesthetic properties. This study demonstrated that geosites combining a local geo-feature and a viewpoint were the most preferred by all groups of respondents. A study conducted in the Lake Tana region in Ethiopia also found that geosites combining a local geo-feature and a viewpoint were rated highest for their scenic beauty [16]. In addition, a survey of 582 visitors in the Albacete mountains (Spain), indicated that among the components that shape the character of the landscape (such as relief, water, vegetation, rural habitat, the combination of human and natural environments, climate, rural landscape, environmental quality and wildlife), the most important landscape component (as indicated by the rankings of these components by the respondents) was relief which includes mountains, gorges and valleys [108]. Such landscape components combine a local geo-feature and a viewpoint, and hence also support our finding. Geosites that combine a local geo-feature and a viewpoint can cater to the needs of different types of geotourists: from those that are purely interested in the geo-feature to those that just want to appreciate the scenic beauty of geosites.
Although there can be factors which might influence the perception of observers in a photo-based scenic beauty assessment such as the height of the horizon in the photograph and the shape of the photograph (e.g., square vs. wide angle), we believe that this had a minor impact on scenic beauty ratings of the geosites in our study. This is because the representative photographs of the geosites were carefully selected from many photos the first and second authors of this paper took during the field excursion in 2019. Where we believed that our photos were not representative, we selected some photos from previous excursions to the study area or from the internet. In addition, we showed two to three representative photos of each geosite to the respondents, in order to provide them with a typical view of the geosite. Moreover, the photo-based survey was conducted online using "Google Forms" and one photo was displayed per page and hence no downsizing of photos was made. In addition, we based our answers to hypothesis two mainly on the comparison of scenic beauty ratings between the NGB-NV and GB-NV groups, who rated scenic beauty based on photos of the geosites.
We acknowledge the limitation of our research in that the geosites selected were mainly related to geomorphology and human-environment interactions, and therefore lack more diversity. The number of participants (n = 176) in the survey was also relatively small and less diverse in their socio-demographic background. Moreover, the NGB-NV group rated scenic beauty based on photographs only. Future research could investigate scenic beauty rating with the following settings: (1) larger sample size of respondents and more diverse socio-demographic profiles; (2) more diverse geosites; (3) respondents with geoscience and no-geoscience background both groups visiting geosites in real life and one group offered scientific information and the other not, and repeating this for multiple groups; (4) scientific value of geosites rated by persons other than those who rate the scenic beauty. Although our study reveals some important relations, future research taking these recommendations into account will allow the drawing of more general conclusions about the relationship between scenic beauty and scientific value, as well as the contribution of geoscience knowledge to perceived scenic beauty rating.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The questionnaire and photos of geosites used in this study are available in the supplementary material here: https://zenodo.org/deposit/4656060#.

Acknowledgments:
The authors thank the Institutional University Cooperation-Bahir Dar University (IUC-BDU) VLIR-UOS project for providing a PhD. research fellowship to the first author. We express our gratitude to Marie Rose Poesen and Ashebir Sewale who helped in distributing the online questionnaire to potential respondents. We also thank the 176 respondents who provided valuable input via the questionnaire. The authors are also grateful to the two anonymous reviewers for their constructive comments.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Table A1. Description of the geosites in southeastern Spain. For their location, see Figure 1b; for their scenic beauty rating, see Table A3 in Appendix A; for their photos, see Questionnaire S1 in the Supplementary Materials.    El Hoyazo Local geo-feature, viewpoint Table A3. Socio-demographic profile of the three respondent groups (n = number of persons).