Examining the Spatial Variability of Geosite Assessment and Its Relevance in Geosite Management

Assessment is the initial step for experts on geoheritage and geotourism when designating geosites in a certain area. During this process, geologically interesting outcrops, formations and places are examined with the use of different criteria to see if they are suitable for geoconservation and geotourism purposes. A quantitative assessment method—Modified Geosite Assessment Model (M-GAM)—was applied in the study area, which is part of the Bakony–Balaton UNESCO Global Geopark in Hungary. M-GAM uses a weight factor (importance) that expresses the opinion of geotourists about 27 infrastructural, tourism and scientific indicators. This factor was examined by questionnaires at nine geosites in the area. At each site, we determined a unique importance value, which shows significant difference from site-to-site and reflects the opinions of visitors about the geosite. The M-GAM method is originally aimed at applying a common weight on each of the 27 criteria during the assessment of selected sites. While this approach is valid, we demonstrated that the method can be extended because the weights spatially vary and can be used to draw conclusions on geosite management. Practically, the evaluation of the factors obtained in this way offers an individual development plan for every site. The current state of improvement direction, the level of communication and the interpretability of the geo-objects can also be determined. In this way, we can get a more realistic development strategy for the geosites.


Introduction
Geotourism has multiple aims and purposes in the life of the community of an area. It presents geological objects as tourist attractions to the public and strengthens the identity of the local people. Its effect is notable on the area's economy as well as on job creation. As tourists without earth science background also visit these places, the sites should be well interpretable with the help of information boards, booklets or other kinds of infrastructure (Dowling 2011, Grant 2010. These places are called geosites (e.g., Brilha 2002;Joyce 2009). Geotourism has an important bond with geoheritage.
The spread of geosites is uneven on the Earth. There are places that have many geo-attractions, but some regions lack them-there are those kinds of objects that became tourist destinations, which would not stand out in geologically more varying regions. It is a difficult task to choose the sites that are worth presenting and need protection. Usually, experts who can determine the geoscientific, infrastructural and tourism value of a geosite do the assessment process. This explains why all infrastructure and tourism purposes stand on the basis of the scientific value of a geosite (Brilha 2016).
Quantitative models that include these three types of values are suitable to enhance the objectivity of the assessment and to make geosites comparable to each other (e.g., Reynard et al. 2007;Vujičić et al. 2011;Brilha 2016), while qualitative methods can pick those phenomena which are interesting in a geologically less diverse region or elements that are parts of a long-term strategy of establishing geoheritage promotion (Migoń 2018). There were several attempts to combine quantitative and qualitative methods to benefit from the numerical approach of the first one, and the recognition of local importance of the second one (e.g., Bruschi and Cendrero 2009;Pereira and Pereira 2010;Tomić and Božić 2014). Based on the outcomes of the evaluation, we can distinct those characteristics of a site that may generate geotourism. The result of an assessment refers to the actual state of the investigated sites whether it was created using quantitative or qualitative methods. Inventories of geosites need monitoring and reassessment from time to time due to anthropogenic or natural changes (e.g., infrastructural development, weathering and natural disasters), which can reduce or increase their geotourism potential (Brilha 2016). For example, the dinosaur footprints around Sintra, Portugal, are under massive erosion, while the newly built Paiva Walkways (near to Arouca, Portugal) experience higher visitor number owing to a significant infrastructural development (Brilha 2016).
If the geotourism potential of a geosite is high, it may become a subject of both tourism development and geoconservation. The assessment process has high responsibility, because only thorough and accurate evaluations may lead to appropriate results: just the truly important geosites should be under "geotourism protection". Such geositesdespite being frequented tourist attractions-can more successfully avoid harmful environmental impacts and investments compared to the unprotected ones (Leader European Observatory 1997).
In Hungary, the largest geopark with 3200 km 2 area covering parts of the Balaton Uplands National Park and the High Bakony Protected Landscape was established in 2012 (Bakony-Balaton Geopark 2019). Other geoparks in the country were also founded around this time, and several individual geosites in the country were surveyed as potential tourist destinations (Horváth and Csüllög 2011;Horváth and Lóczy 2015;Szepesi et al. 2017). The assessment methods in these works were often based on the methodology of the quantitative GAM-Geosite Assessment Model (Vujičić et al. 2011). However, the establishment of the Bakony-Balaton UNESCO Global Geopark (BBG) was preceded only by qualitative evaluation (Korbély 2011). This is mainly because the areal extent of the BBG is large and geo-attractions are usually concentrated as clusters, which were considered as single geosites, despite each of these clusters contain several individual attractions. For example, the basalt volcanoes of the Tapolca Basin are parts of one single "volcanic landscape" geosite (Korbély 2011;Harangi 2014).
The BBG hosts one of the most visited tourist destinations in Hungary: Lake Balaton. Although most of the tourists do not choose geosites as their primary destinations, the geotourism potential is definitely present (e.g., one of the largest geo-attractions, Lóczy Cave is visited by approximately 35 000 tourists per year). Quantitative assessment of the BBG geosites in this frequented region would be an advantage when doing strategic planning for infrastructural developments. In 2017, a detailed assessment process was started by  in the eastern part of the BBG using the Modified Geosite Assessment Model (M-GAM) of Tomić and Božić (2014). The assessment's results identified the most attractive geosites in the examined area .
M-GAM provided an opportunity for us to analyse the modificatory effect of the visitors' opinion on experts' evaluation. The model was derived from the GAM of Vujičić et al. (2011) to assess the geotourism potential of the Lazar Canyon area (Serbia). M-GAM differs from GAM by containing an Importance (Im) weight determined by interviewing geotourists at various geosites on their personal opinion of the importance of various indicators of the model. Then, the Im values were averaged to all geosites and used in the assessment as weights ranging between 0 and 1 for each indicator. The way M-GAM was applied by Tomić and Božić (2014) tells us which indicators are important for the visitors. However, by averaging the opinions, we lose information about the difference of these weights among the geosites.
When applying M-GAM in the assessment of the BBG geosites, our hypothesis was that the visitors' opinion of the criteria could give us a clue about their expectations concerning the geosites. Thus, the assessment of these sites does not only identify the most significant locations but also helps the geopark management to plan future infrastructural developments .
Such assessment is useful for management purposes, and the site-specific opinion (the Im-s) of the visitors of each criterion is crucial for every geosite. We suspected that the averaged Im values would distort the assessment if doing sitespecific assessment as they do not explicitly represent the geosite that is under evaluation. Based on the results of the previous assessment, we have clustered the most visited natural attractions of the study area to further examine how the visitors weight the GAM indicators' importance. This study treats the Im value as a variable that can be effectively used separately at every examined geosite. The spatial variability of the geosite assessment in the present study is expressed in the site-specific evaluation of the collected data: the Im values change from geosite to geosite. Our aim was to examine how visitors see the actual scientific content and tourism facilities of the selected geosites by applying the M-GAM assessment method in a new way.

Study Area: the Southeastern Part of the Bakony-Balaton UNESCO Global Geopark
One of the most visited tourist destinations in Hungary and Central Europe is Lake Balaton and its surroundings. Visitors with various interests can entertain themselves in the area: there are beaches on the shore, famous wine cellars and rich cultural heritage. Due to the variety of natural (geological, biological) values, this area is even more characteristic. The northern shore of the lake is part of the Transdanubian Range and is well-known for its rare geological formations (Korbély 2011). The climate has Mediterranean nature near the shores, but on higher areas, it is temperately cool and wet. The mountains and hills do not make the formation of rivers possible, although there are small streams called "séd" in the deep valleys. The "séd" valleys have cool microclimate even during hot summer days (Dövényi 2012;Futó 2013).
An extensive geological mapping was done here in the 1990s and a 1 : 50 000 scale map was published in 1999 (Budai et al. 1999). The first geological hiking map of the area was published in 2018 ). The rectangular shape of the examined territory ( Fig. 1) covers the extent of the hiking map as well, because the editing work was done in parallel with the assessment. Until now, this 1 : 30 000 map is the most detailed published geological map of the area.
The geological history dates back to the time when life was only present in seas. In the following half a billion years, volcanoes, deep and shallow seas, rivers crossing semi-deserts, tropical lagoons, lakes, streams and wind transformed the landscape. Today, humans are also remarkable surfacechanging factors . The oldest formations here are from the Silurian, which are more than 410 million years old. The Lovas Slate metasediment hides the small bodies of volcanic rocks called Alsóörs Metarhyolite. The outcrops and quarries of these formations are extremely rare Fig. 1 The extent of the Bakony-Balaton Geopark (red border) and the examined area (green border) and highly protected due to geoconservation aims for the preservation of their visible scientific values. The metamorphism of these sedimentary rocks took place about 325-310 million years ago during the Variscan orogeny (Haas 2013).
A specific rock type of the area was formed during the Permian in fluvial and limnic environment: the red Balatonfelvidék (Balaton Uplands) Sandstone, which is commonly used as a building stone. A large part of the study area's outcrops is from the Mesozoic sedimentary sequence, which mainly consists of shallow-marine carbonates of Triassic age. The area tectonically belongs to the East Alpine Nappe System, so such rocks as in the Alps can be found here from this era referring to the large-scale tectonism that has shifted these formations to their current place (Haas 2013). The thickness of the Triassic carbonate sequence in the Bakony-Balaton region can reach 2.5-3 km. Although Jurassic and Cretaceous sequences are not represented in the study area, they are found in the vicinity further to the north (Budai et al. 1999).
As the area was denuded during the Paleogene, only from the late Miocene can we find sediments next in connection with the development of the Pannonian Inland Sea (Korbély 2011). At this time, volcanic activity has also started forming phreatomagmatic explosions, basalt magmatism and postvolcanic geyser activities in the region (Németh et al. 2001). Only the latter is represented in the study area .
During the Pleistocene, climatic fluctuations created variable landforms and peculiar sediments like aeolian loess sheets and alluvial fans. The area was under periglacial conditions during the glacial periods, and wind erosion was dominant. The basin of Lake Balaton formed during this time due to the combined effects of gradual tectonic subsidence and deflation (Cserny and Nagy-Bodor 2000). In the more humid interglacial periods, fluvial erosion was dominant and the alluvial fans prograded into the Balaton Basin. The slow uplifting of the area during the Pleistocene helped the meteoric water to form caves in the carbonates and allowed the southflowing creeks to cut through ridges creating spectacular landforms on the surface (Korbély 2011).

The Applied Assessment Models
The assessment of the study area was carried out with the M-GAM method, which practically involves the quantitative method of the Geosite Assessment Model (GAM). GAM has already been successfully used in some countries, including Hungary (e.g., Szepesi et al. 2017). The geotourism potential of the evaluated geosite is described by the sum of two main variable groups: the main values (MV) and the additional values (AV).
The main values group consists of the scientific/educational (VSE), scenic/aesthetic (VSA) and protection (VPr) subgroups, while the additional values group is formed by the functional (VFn) and tourism (VTr) subgroups (Table 1). The MV groups contain 12, and the AV groups contain 15 parameters, which are called indicators altogether. Geosites are evaluated by scoring these indicators with 0 (the worst), 0.25, 0.5, 0.75 and 1 (the best) marks. The final GAM values are calculated using the following equations: These are the final results of a GAM assessment (Vujičić et al. 2011). The detailed process of scoring is found in Tables 1 and 2.
Although this is rather an objective way of assessment, a geosite usually holds some other kinds of values that are not considered in this method. There might be other important factors, such as psychological, historical, entertaining and religious characters that could be subjects of individual preference. In this sense, the geological "story" is superimposed by personal experiences, memories and aims. All these aspects are gathered in the Im weight of M-GAM.
The authors of M-GAM asked 96 geotourists at 3 different geosites to give scores for the 27 GAM indicators according to their opinion (on the scale of 0, 0.25, 0.5, 0.75 and 1). It was complemented in 2015, when a group of 197 visitors filled the questionnaires but not at the same sites (Božić and Tomić 2015). The visitor's votes were then averaged to get the Im value of a certain indicator. They calculated the M-GAM value of a site by multiplying the GAM scores by the Im values: With the use of this weighting method, the assessment includes the opinions of all interviewed geotourists strengthening the objectivity of the evaluation.
The authors of M-GAM used scatterplot matrices (Fig. 2) in order to express and visualise differences between the GAM and M-GAM values. The position of geosites in different matrix fields refers to their current status in geotourism. Values on the horizontal and vertical axes show the main and additional value scores, respectively; the higher the score, the more significant the geosite is. Although the Im values identified by Tomić and Božić (2014) helped to highlight the expectations and interests of tourists generally, it referred to the whole Lazar Canyon area, and to an even larger territory in their subsequent work (Božić and Tomić 2015). The present study used another approach of the M-GAM method by measuring the Im weights for all the geosites individually using on-site and online questionnaires. In this way, we were able to observe the spatial variability of the Im values, that is, the Im differences between each examined geosite.

Collecting and Processing Data from the Study Area
By asking geotourists to fill out questionnaires at every examined geosite, we have the opportunity to produce geositespecific data: the Im weights are determined independently of each other. A comparison of the results has been done to analyse the spatial variability of geosite assessment and to determine the unique character of each location. It must be noted that collecting site-specific answers requires either personnel or some infrastructure. Since the selected nine locations are not all frequently visited by tourists, it was decided to set some infrastructure at all sites by placing waterproof sheets on information boards calling for online questionnaire fill, and collect the answers mainly through them. The total number of fills gathered is 247 during a 14-month period.

Online Interviewing
To produce data comparable with other geosites in Hungary and especially with the originally published Lazar Canyon territory, our questionnaire was based on the original descriptions of GAM/M-GAM and Im. We asked the visitors how important they consider each indicator of the model. Tourists were not differentiated throughout the survey, although we collected demographic data for further research. We have designated the nine most important geosites of the area according to the preliminary GAM-M-GAM evaluation (Fig. 3). Most of them have relatively developed infrastructure (explanatory signs, information boards, benches, bins, parking places, etc.), but they might be utilised better with the help of this research.
To ensure that results come from only tourists or locals who definitely visited and know the geosite, we put laminated paper sheets next to the geosites (e.g., on the wooden frame of an information board). On each of these sheets, there was a short description, a QR-code and a link to a site-specific Google Form that made the on-site fill possible. On the webform, our research topic's importance was shortly described and all 27 indicator-related questions were asked. Visitors had to evaluate every indicator on a 0-to-1 scale: 0-not important, 0.25-rather not important, 0.5-moderately important, 0.75-rather important, 1-important. We have provided a box where geotourists, who are interested in this project, may type their e-mail address and we give them information about our work.

On-Site Interviewing
On-site interviewing took place at sites where the number of online fills did not rise in the expected pace (i.e., 2-4 fills per month). These were typically the ones with low level of infrastructure. The paper-based questionnaire was the same as the online version. This type of interviewing was not so successful as we expected: only~10% of the fills were collected in this way. Most of the visitors we asked were not confident and enthusiastic about scientific topics: they may have liked solitude when filling these types of questionnaires.

Results of the Evaluation
The interviewing took place between July 2018 and September 2019. During this time, 247 fills were collected from the designated geosites (Fig. 4). The tool for data processing, analysis and visualisation was MS Excel. The Im value that indicates the importance of a geosite was calculated in two ways: (1) the calculation of the arithmetic mean of fills per indicator per site was quasi-similar to the original M-GAM methodology, (2) each of the fills was considered as an individual assessment, and a site-wise statistical evaluation

The Local Im (Importance) Values
By summarising Im data from each geosite, we get unique geosite-specific values for each GAM indicator. If we insist on using the originally published methodology of the M-GAM evaluation, a further averaging of the Im values would be necessary for all the indicators. Instead, we compared our data with the available Im weights from elsewhere to analyse the differences.
According to the original description of M-GAM (Tomić and Božić 2014), the Im values reflect a general opinion about the importance of the model's indicators. They suggested that the identified Im weights can make the evaluation more objective by pointing out those values, which are of paramount importance for tourists. Several studies applied the determined values, including ourselves during the assessment around Csopak village in 2017 . This new method seemed to show more realistic results owing to the GAMmodifying Im factor. The analysis of the difference between the two models' results gave the opportunity to give hints about development strategies, but the question of the relevance of the averaged weights have emerged. Tomić and Božić published new importance values in 2015 (with 293 fills). These were produced by summing the 96 fills from the original survey from the Lazar Canyon and adding 197 from other sites. The difference between the two Im weights is surprisingly high, suggesting that the Im value is not a constant, but a variable that may have spatial dependency. However, these aspects were not mentioned in the publications.
The Im values of the two archive surveys and the Im-s of the present survey (Figs. 5 and 6) are greatly varying. There are two columns for the global factors (blue circles) and nine columns for the examined geosites in the Balaton Uplands (red circles).  Tomić and Božić (2014) and Božić and Tomić (2015), compared with the present results There is even larger variability between the global and local AV Im-s (Fig. 6) than in the case of MV Im-s (Fig. 5). The determined values expressively refer to the character of a geosite: for example, the difficulty level of accessibility (if easy, larger score, if hard, lower score) or tourism infrastructure (if exists, higher score, if does not exist, lower score).
The Im values for the study area represent the weight of the certain indicator, but by comparing the values calculated for  Tomić and Božić (2014) and Božić and Tomić (2015), compared with the present results the sites, the spatial variance of the indicators can also be given (Table 3). For this, we used the VAR.P function of Excel.

Site-Wise Statistical Evaluation
The individual Im values of an interviewed visitor produce a single M-GAM evaluation for the given site. Because of this, the number of questionnaire fills is equal to the number of individual assessments available at a certain geosite. By plotting these unique assessments into a scatterplot diagram of MVs and AVs (considered as independent variables), we can visualise how differently the interviewed visitors think about the geotourism potential of a geosite (Fig. 7). When plotting these values as x( MV ), y( AV ) functions, we get point clouds with various shapes in the scatterplot diagrams. These diagrams show how the visitors interpret the importance and significance of each geosite, and how these evaluations differ from each other. It is evident that the scientific value (MVs) and local infrastructure value (AVs) fit differently to the visitors' expectations and needs, but the shapes of the clouds statistically reflect to all the interviewed visitors' opinion about the examined sites. The point clouds can be differentiated based on their shape: LI = large isometric; SI = small isometric; EA = elongated, additional values dominated; EM = elongated, main values dominated.
The large and isometric (regular shaped) clouds refer to the diverse interests, expectations and scientific-infrastructural needs of geotourists. The variance (σ 2 ) is relatively high among these values (e.g., Lake Köcsi and Miske Cliff), because those who go there weight the indicators very differently, which means that the people who visited the site and filled the form had varying knowledge about the site. This could refer to that the geosite is not well communicated or promoted, and its features do not fit to the visitors' expectations. This situation could be treated by clear, well-designed commercial material, internet sources and interpretive panels.
If the cloud is small and isometric, the visitors have quite similar opinions (either good or bad) of the features of a certain geosite (e.g., P/T key section and Forrás Hill). In this case, each individual assessment gets nearly the same M-GAM value. This means that only those people who visited the site (filled out the form) were quite equally aware of the AV and MV state of the geosite.
When the shape of the point set is not regular, we can differentiate the visitors according to the MV and AV values. These elongated clouds could refer to different evaluation groups, which can be categorised by a linear function model fitted on the point cloud.
(1) The negative line-fit: lower values on the y (AV) axis but higher on the x (MV) axis refer to visitors who prefer the scientific values over the infrastructural ones, but there are others who consider geoscientific features not so important, so they gave lower Im scores for the x (MV) and weight the AVs as important indicators. If both types occur on the scatterplot, a line with negative steepness would fit on the points, meaning that people who expect very different things visit the site. Practically, the geosite is not well communicated (2) The positive line-fit: the lower values on the y (AV) axis have generally lower values on the x (MV) axis as well and occur together with assessments having higher values on both axes. A line with positive steepness would fit on the points meaning that there are people who are generally unsatisfied and people who are generally satisfied with both MVs and AVs. Normally, this situation is the most common. Finding the good management of these sites is quite difficult, because there is a need for infrastructural and scientific development too When analysing these non-isometric elongated clouds, we can use the steepness (S) of the linear functions plotted on the M-GAM score set. This steepness refers to the difference between the visitors' opinions. We can determine if there is higher demand for scientific or communicational-infrastructural development: when the steepness is greater than 1, additional values, if less than 1, main values are more important for the visitors in general.
The R 2 coefficient shows the goodness of the fit of a linear function on the personal assessment scores. If R 2 < 0.2, the fit is rather uncertain. If the R-square value is high, the cloud fits to the trend line and we can draw conclusions about further strategies. For example, the Koloska Valley geosite has nearly the highest coefficient value (R 2 = 0.353) and the steepness is The σ 2 is the variance of the dataset. This value shows how the plotted values differ from the statistical mean of the whole set of fills on each geosite. It also shows the extent of the difference of the visitors' personal evaluations.

Discussion
We have gathered 247 fills varying between 8 and 64 per site (Fig. 4). As it is not a statistical amount of data, we consider this study as an exploratory research. We found that interviewing at a geosite is a hard task because of the sporadic visits and a common indifference among tourists concerning such queries. The length of the questionnaire can also be caused for the relatively small number of answers: thinking over 27 questions and answering them requires patience.
However, the installation of QR codes directing to the online questionnaire seemed to be working as 90% of the collected data were gathered through this way.
In Table 4, results, consequences and proposed development strategies are summarised regarding the nine examined sites. The current state of the geosites is derived from the individual Im values, while the development strategies are summarised from the site-wise statistical evaluation.
The variances (σ 2 ) of the indicators show that geotourists consider each geosite important due to different marks, facilities and services (Table 3). This proves the uniqueness of the sites and highlights the importance of separate assessment factors. The sigma-square values in Table 3 are determined for the whole study area and, thus, they represent the spatial variance of the indicators. Among the MV indicators, MV12 (suitable number of visitors) shows the largest spatial variance (σ 2 = 0.012) and the MV9 (current state) the smallest one (σ 2 = 0.001). Considering the AV indicators, the AV13 (tour guide service) has the largest (σ 2 = 0.025) and the AV14 and 15 (hostelry and restaurant services) the smallest (σ 2 = 0.002) spatial variance, although the latter ones are considered equally unimportant by the answerers with Im values barely reaching 0.4 (Fig. 6). However, the AV10 (interpretation panels), showing also a small spatial variability (σ 2 = 0.004), was the most important indicator for the visitors with Im values around 0.9 at all the 9 sites ( Fig. 6 and Table 3). Our results also show that there is a large difference in the variance of the personal evaluations ranging from 1.108 (Vöröskő Nature Trail) to 4.116 (Ember Cliff) on the examined sites ( Fig. 7 and Table 4). These support our hypothesis about the spatial variability of the geosite assessment using the M-GAM method.
The collection of data continues on these nine and other selected sites on the adjacent areas of the BBG. The increasing number of fills does not only improve the statistical parameters of the survey but also if it lasts long enough and lives through several tourist seasons it can detect the change of the Im value in time. We suppose that besides the spatial variability, the modificatory effect of time can also be a property of geosite assessment because infrastructural developments and good promotional solutions can modify the calculated importance values. We hope that conclusions can be drawn about infrastructural and/or scientific questions in the following years (e.g., building of interpretive panels, visitors' centre and development of geoeducation). In the case of two examined sites, Lóczy Cave and Vöröskő Nature Trail, large infrastructural development is to be built, so the variability of Im values in time could be examined there. The site is not well-promoted (large isometric cloud). The tourism infrastructure is not maintained correctly.
Clear and well-designed commercial materials are needed to filter the diverse expectations of visitors. Vöröskő Nature Trail EM S:0.644 R 2 : 0.319 σ 2 :1.108 Good AV Im-s because of the neighbouring outdoor theatre. The scientific values are not balanced with it.
A visitors' centre is being built now. It will put more emphasis on geoheritage-MV Im-s and AV Im-s will be balanced. P/T key section SI S:0.007 R 2 : 0.001 σ 2 :1.819 The MV Im-s and AV Im-s are well-balanced. The location is not suitable for more tourism infrastructure. A car parking could make visiting more comfortable. Miske Cliff LI S:0.323 R 2 : 0.066 σ 2 :3.845 Good MV Im-s but poor AV Im-s because of the lack of tourism infrastructure.
Interpretive panels, good commercial materials and well-marked hiking routes are needed, it is hard to find the location from Felsőörs. Forrás Hill SI S:0.156 R 2 : 0.012 σ 2 :1.886 The site is not well-promoted (large isometric cloud). The tourism infrastructure is not maintained correctly.
The MV Im-s and AV Im-s seem to be in balance, but clear promotion is needed to reach a regular-shaped point cloud The site is visited not only for the geological heritage-it is a nice picnic area.
Better marketing could point out the large geotourism potential of the valley. Improvement of the infrastructure is needed.

Conclusion
Perfect assessment models do not exist. At least the subjectivity of the evaluating person always emerges while inconsistent and uncertain data and data sources may cause precarious and misleading results. We have used the M-GAM model during the analysis. The determination of the importance (Im) values should take the spatial variability of the evaluation into consideration. To count with it, we needed to collect visitors' opinions from each of the nine examined geosites in the research area in the Balaton Uplands. We demonstrated that the Im value of M-GAM is suitable for determining the visitors' opinion about the current state of the geosites. For this, the assessment was completed with the acquisition of on-site and online questionnaires at nine geosites of the BBG. The difference of the original and the newly determined Im values is explained by the difference between their meanings. While the previously published Im values are aimed at expressing a general weight for the GAM indicators for an examined area, the ones presented in this study are unique weights expressing the visitors' expectations for each examined geosite. It can be analysed to plan a better infrastructure or to put more emphasis on the interpretation of the scientific content. When interviewing, it is very important to ask wellcomposed and clear questions. However, this GAM-based method is quite complex with its 27 indicators and it should also be considered to work out, or apply a less complex model for on-site questionnaires.
Among the improvement of geosite assessment, our aim is to emphasise the role of geoheritage. Effective protection should always precede developments. Proper assessment using different criteria helps to compare geosites that are suitable for conservation, tourism or educational purposes. The nine examined natural values are all under international protection. They can only be the places to spend our free time actively if we properly take care of them.