Next Article in Journal
An In-Depth Analysis of Parking Patterns in a Typical Chinese Danwei via Customized Data Collection App
Previous Article in Journal
Defining a Model for Integrating Indoor and Outdoor Network Data to Support Seamless Navigation Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Refinement Proposals for Geodiversity Assessment—A Case Study in the Bakony–Balaton UNESCO Global Geopark, Hungary

1
Doctoral School of Earth Sciences, Faculty of Science, ELTE Eötvös Loránd University, Pázmány P. sétány 1/C, H-1117 Budapest, Hungary
2
Institute of Cartography and Geoinformatics, Faculty of Informatics, ELTE Eötvös Loránd University, Pázmány P. sétány 1/A, H-1117 Budapest, Hungary
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2021, 10(8), 566; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10080566
Submission received: 5 July 2021 / Revised: 16 August 2021 / Accepted: 20 August 2021 / Published: 22 August 2021

Abstract

:
Geodiversity is the variety of natural elements that are excluded from biodiversity, such as: geological, geomorphological, and soil features including their properties, systems, and relationships. Geodiversity assessment measures these features, emphasising the characteristics and physical fragility of the examined areas. In this study, a quantitative methodology has been applied in Bakony–Balaton UGGp, Hungary. The Geopark’s area was divided into 2 × 2 km cells in which geodiversity indices were calculated using various data: maps, spatial databases, and elevation models. However, data sources differ significantly in each country: thematic information may not be entirely public or does not have the appropriate scale and complexity. We proposed to use universal data—geomorphons and a watercourse network—derived from Digital Elevation Models (DEMs) to calculate geomorphological diversity. Making a balance between the base materials was also an aim of this research. As sources with different data densities are used, some abiotic elements may be overrepresented, while others seem to have less significance. The normalisation of thematic layers solves this problem: it gives a proportion to each sub-element and creates a balanced index. By applying worldwide accessible digital base data and statistical standardization methods, abiotic nature quantification may open new perspectives in geoconservation.

1. Introduction

The quantitative assessment of biotic and abiotic nature elements is a relatively new research field in geosciences. Biology was the first to determine the concept ‘biological diversity—biodiversity’ at Rio Earth Summit in 1992 [1,2]. It is the variability of the living organisms of the Earth [1], or more precisely: ‘the diversity within species, between species and of ecosystems’ [2]. The introduction of the term ‘geodiversity’ came shortly after the Rio Earth Summit. It first appeared in Sharples’s work [3], and evolved in the following years [4,5]. According to Gray [6], geodiversity is ‘the natural range (diversity) of geological (rocks, minerals, fossils), geomorphological (landform, processes) and soil features. It includes their assemblages, relationships, properties, interpretation and systems.’ Gray [7] analyses this concept and states that geodiversity is ‘a theoretical framework of ideas’ and ‘a set of assumptions, concepts, values and practices that constitutes a way of viewing reality’.
Finally, Gray [8] extends the preceding definition. It is considered as the backbone of geoheritage, valuing abiotic nature, and geoconservation. The publication reveals a connection between modern society and geodiversity: humans have always used and earned benefits from abiotic features of the Earth for various purposes throughout history. It is also stated that geoconservation programmes and geoheritage are often based on the calculated spatial evaluation of geodiversity.
This definition also solves the contradiction between the triad of terms geodiversity, geological heritage and geoconservation by defining a practical approach for each of them. While geodiversity refers to all abiotic elements of nature, geological heritage (or geoheritage) is only a subset of geodiversity elements (with geological, educational or tourism importance). Geoconservation can be considered as a toolset: it is responsible for the identification, evaluation, conservation, and tourism promotion of geological heritage [9,10].
Assessment models to measure geodiversity evolved in parallel with the definition itself [10,11,12,13,14,15]. Most of them are based on cartographic and spatial database materials [16,17,18]. Processing geological, geomorphological and pedological features within an area for assessment tasks is done in a GIS (geographic information system) environment. Studies use various types of national data sources all over the world—even in tropical and developing countries [19,20,21].
In Hungary, no geodiversity assessment has been made before. However, many areas are promoted as having geoheritage and geotourism values and are under geoconservational protection. There are two UNESCO Global Geoparks in Hungary (one common with Slovakia) and other protected areas, as well as National Parks that focus on geological and earth scientific heritage. This research aims to set a basis for further assessment works by following these guidelines:
  • The methodology should be used anywhere.
  • Data sources (maps, databases) should be available to researchers.
  • The method should reach a balance between different data sources—none of them should be overrepresented—this is statistical normalisation in this case.
  • Try to build universal variables (data sources) into the model (instead of country-specific data) without loss in scientific value to establish the potential of international comparison. These are digital elevation models and derived data from them (geomorphons and watercourse network).
The geodiversity assessment of Bakony–Balaton UGGp is the first attempt to follow these pillars to determine the geodiversity index. Apart from theoretical benefits, this study can open new perspectives to the Geopark by highlighting areas with high geodiversity.

2. Geography and Diversity of the Area

The sample area—Bakony–Balaton UNESCO Global Geopark—is located in the western part of Hungary, next to Lake Balaton, the largest lake in Central Europe. The 200 km long and 30–50 km wide Transdanubian Mountains cross through the Geopark, with an average elevation of 400–600 m. It is mostly built up by Mesozoic limestone and dolomite, with volcanic areas and karst horsts. Between the hills, basins are covered with Tertiary/Quaternary deposits. When strong tectonic forces were present, horsts, grabens and basins were formed by deflation that exhumed the basalt buttes of the area [22]. The highest parts of the Geopark are in the Bakony Mts. The amount of precipitation is in tight connection with the NE-SW direction of the mountains: the annual rainfall is around 600 mm in the NW parts and 500–600 mm in the SE parts. In a geomorphological context, the area is divided by distinct faults into smaller units with various attributions: lifted highlands (Keszthely Mts), horsts and grabens (Balaton Uplands, Northern Bakony) and carstic–basaltic undulated relief (Southern Bakony) also appear. Due to the accentuated surface, the vegetation cover shows great variability. In the sub-Atlantic Bakony Mts, beech woods are the most common due to the relatively large amount of rainfall. As the elevation decreases, oak and hornbeam woods become common. The sub-Mediterranean climate of the surrounding of Lake Balaton is good for chestnut and almond trees too [23,24].
The Bakony–Balaton UNESCO Global Geopark (Figure 1) is founded and managed by Balaton Uplands National Park. It has been a member of the Global Geoparks Network (GGN) since 2012. It has an area of 3244 km2 and encompasses 151 settlements. Its area is mainly located in the Bakony Mountains. Since the Geopark has a geological history of 440 million years, its earth scientific and cultural heritage is internationally outstanding. Many brochures, educational trails and visitors’ centres help geotourists in gathering scientific and cultural information [25,26,27].
The eastern part of Lake Balaton is bordered by the 15-thousand-year-old Kenese High Shores that are built up by the marine sediments of the Upper Miocene Pannonian Sea [26]. This area is also home to a botanical rarity, Crambe tataria (Brassicaceae), which is an endemic species of the Pontic-Pannonic region [29]. The most common building stone of the Balaton Uplands is the Permian red sandstone. There are educational trails (e.g., Alsóörs) and rehabilitated quarries (e.g., Balatonalmádi) with explanatory panels introducing the scientific background of this geological formation. Lóczy Cave is the most famous geosite in the eastern part of the Geopark, where the Upper Triassic Füred Limestone layers (with thermal water effect on them) can be observed [30]. The area is also famous for its wines.
The Tihany Peninsula is well-known for its Benedictine monastery and picturesque landscape. Both biotic and abiotic features are worth mentioning here. Lavender is the most famous domesticated plant of the peninsula—it is grown here for centuries [31]. Even the Geopark’s visitors’ centre in Tihany is called Lavender House. The first phase of Balaton Uplands volcanism took place here 8 million years ago. The stone walls and cliffs of “Monks’ Houses” (Barátlakások) and the geyserite cones of, and around the “Golden House” (Aranyház) geosite also show the immense power of volcanic activity [32].
Famous vineyards and karren-surfaces (karst microforms) make the Pécsely Basin a tourist destination. From the quartz-cemented late Miocene sandstones, block fields (“seas of stones”) with spectacular pseudokarren features had formed on the borders of the Káli Basin that have a cultural value too: millstones have been carved from these blocks throughout centuries [33]. On the eastern part of the basin, an abandoned basalt quarry was converted into the Hegyestű Geological Visitor Site, where all the nearby rock types are exhibited among the basalt ‘organ pipes’ of the 7.9 Ma old volcano [34].
The Tapolca Basin and its residual hills (e.g., Badacsony, Csobánc) made of basalt form probably the most well-known landscape in Hungary. In this area, phreatic volcanism (producing scoria cones, lava fountains and lakes) was followed by erosional processes, and the soft pyroclastic rocks eroded, exposing the columnar jointed structure of the lava lake rocks as a result [35]. The uplift of the area resulted in stronger erosion: wind and water formed the current landscape of the Basin. Cultural heritage is also relevant here because these basalt mountains were important through history: castles were built on top of them (e.g., Szigliget, Csobánc) for defensive purposes. Many endangered plants grow on the karst terrains of the area [36]. Tapolca Lake Cave Visitors’ Centre is located under the city of Tapolca, presenting spectacular karst morphology. Csodabogyós Cave is the most important geosite of Keszthely Mountains being famous for its stalactites and stalagmites [37].
The educational trail of the “Ancient Karst of Úrkút” (paleokarst of Úrkút village) shows the tradition of manganese mining in the Southern Bakony. The first outdoor mining museum of Hungary can be found in Ajka—the surroundings of the city was famous for bauxite mining and alum earth production. The archaeological findings of Szentgáli Rockhole can be discovered during cave tours.
Zirc is often called the “capital” of the Bakony Mountains (located in the Northern Bakony). The Natural History Museum of the Bakony Mountains is located here, presenting the Geopark’s rich earth scientific and cultural history. The famous dinosaurs of Iharkút are exhibited in the Hungarian Natural History Museum (Budapest), but there are open days during fieldworks in the summer [26,36].

3. Assessment Methodology and Process

The quantitative methodology of this study is based on the Portuguese school [10]. There have been other modified assessment techniques rooting from this [15,37], but this original framework suits best to our standardising propositions. The process is built up by various numerical sub-indices that represent geodiversity elements: geology, geomorphology (relief and hydrology), pedology (soil science), palaeontology and mineral occurrences (precious stones, minerals, geological energy sources, mineral waters, and springs). The outcome of the evaluation is the geodiversity index—it is the summed value of the preliminarily determined sub-indices. Its purpose is to express all mentioned aspects in the most balanced way without emphasising any of them.
This study has applied only open-access software (QGIS 3.16.0, SAGA GIS 2.3.2 and GRASS GIS 7.8.4) to ensure reproducibility. The initial step of the assessment was the division of the sample area. The 3244 km2 of the area was sliced up into a grid of 2 × 2 km2 squares (Figure 2). Each square represented a homogenous unit of the Geopark. All sub-indices and the final geodiversity index were defined in these areas. The reason for choosing this resolution was in connection with the extent of the Geopark and the scale of the available data sources to produce the most realistic index values in the cells. If the grid is too dense and we have a small-scale thematic dataset, the relations between the thematic features are not correctly detected. If the grid cells are big and the data large-scale, each grid cell gets large index values, but it does not correctly represent the diversity. Otherwise, if the cells are too small or the data too small-scale, the value of the indices will be too low, and again we will not get a meaningful picture of diversity. In brief, grid resolution depends on the desired complexity of the assessment and the scale of the available data sources.

3.1. The Normalisation of Subindices

Each sub-index requires a dataset with specific thematic data. We used a 1:100,000 scale geological map [22], a 1:100,000 pedological map [38], a 90 m-resolution DEM [28], and explanatory books of geological formations of Hungary [39,40,41]. Usually, it is hard to find and use layers of the same scale and complexity. Scale differences may result in the overrepresentation of more complex thematic layers. This problem can be improved by the normalisation of sub-indices:
(1/SNDXmax) ∗ SNDXcell = SNDXnorm,cell
The multiplicative inverse (or reciprocal) of the maximum value of a certain sub-index multiplied by the sub-index value of a specified cell gives the normalised value for the cell. Therefore, the maximum of a sub-index can be 1, and the other values range from 0 to 1 proportionally. If this process is repeated in case of every thematic sub-index, a normalised geodiversity index is calculated by summing them up. The normalisation reduces the distortion of scale and complexity differences from the assessment process, and thus, a more balanced geodiversity index can be produced.

3.2. Geological Sub-Index

The geological sub-index was calculated using the 1:100,000-scale geological map of Hungary, cropped for the Geopark [22]. The sub-index was determined by counting the number of different lithostratigraphic units in each grid cell (Figure 3). For this, the geological vector layer was rasterised and the ‘Zonal Statistics’ geoprocessing tool’s ‘Variety’ attribute gave the different raster values (geological types) of each cell. After this, the normalisation formula was applied to get balanced sub-index values.

3.3. Pedological Sub-Index

This sub-index was calculated by using ‘AGROTOPO 100’ [38] that is a 1:100,000-scale soil map and digital database containing genetic soil classes. Using the different soil class occurrences in each cell, the sub-index was calculated in a way similar to the geological sub-index. After rasterization and running ‘Zonal Statistics’, the normalisation formula was applied (Figure 4).

3.4. Geomorphological Sub-Index

The geomorphological sub-index is the composition of relief and hydrological indices. In preceding publications [10,15], the relief index was calculated by using geomorphological maps, and the method was similar to the geological and pedological indices we described: the number of different geomorphological units determined the score for each cell (some models counted in the borders between morphostructural units too). In the study area, a medium scale geomorphological map was not available, so we used the MERIT DEM [28] digital elevation model for this task, which is an improved version of the SRTM 3 surface model. The DEM was categorised into morphological types called geomorphons, based on the morphometric derivatives. The theory of geomorphons [42] is a relatively new concept of terrain analysis. It uses 8-tuple patterns of the visibility neighbourhood to detect pixel values and determine terrain forms in the eight principal directions to a surface point. Two parameters can be manually changed: search radius and flatness threshold—they affect the complexity of form recognition. By analysing DEMs this way, we can delineate 10 different morphological types (Figure 5). The production of relief sub-index was produced similarly to geology and pedology from this raster data.
The hydrological index counts with the effect of hydrological elements on geomorphology. As in the case of the relief index, we tried to find a universal way to construct a watercourse network for the sample area. To do so, the hydrological modelling tools of SAGA GIS were used. After generating a depressionless elevation layer (‘Fill sinks’ tool), channel network can be extracted (‘Channel network’ tool). We have also used a ‘Flow Direction’ raster (also generated with SAGA) that enhanced the results. However, current hydrography does not always reflect the phenomena expected from geomorphology. This is because the climate may have changed since the valley was formed, or the stream may flow in a different direction due to erosion or tilt, or it may disappear underground because it is in karstic rock. In any case, the valley is there. Due to this, the generated channel network was corrected with the use of OpenStreetMap and large-scale topographic maps to avoid that the extracted network contains valleys where water is not present. With proper parametrisation and resolution, it is not too time-consuming—mostly limestone and karstic formations need moderate supervision. The resulting layer can be categorised by applying the tool ‘Strahler order’ [43]. According to this system, the lowest hierarchical group (with a value of 1) encompasses streams while the highest value (5) is assigned to large rivers, lakes, and coastal areas. The value of the index in each cell is the highest hierarchical level divided by 2, rounded up to the nearest integer (Figure 6). For instance: a score of 3 can be assigned by having a level ‘5′ watercourse (5/2 rounded to 3) in a certain cell. Zero value is determined if there is no hydrological element in the grid cell.
Finally, the geomorphological sub-index was calculated by summing the relief and hydrological indices and then performing normalization on the summed cell values.

3.5. Palaeontological Sub-Index

Unlike the previous cases, this sub-index was not determined with the use of geodatabase sources since freely available medium-scale fossil map was not available. However, the method based on the geological map we used for calculating the geological sub-index. The original dataset of the geological map does not contain information about the fossils, but the explanatory books of geological formations of Hungary [39,40,41] do. A new attribute was created for all the lithostratigraphic units using their descriptions. Fossiliferous formations were noted (based on the formation descriptions), and the number of different fossils per formation was added to each geological polygon. This value in every cell was considered as the palaeontological sub-index. The maximum number of different fossils in a cell was 9 and the minimum was 0. The average number of fossils in the cells was 2.31. Normalisation was also applied in this case.

3.6. Mineral Occurrences Sub-Index

This sub-index was also determined with the help of the geological map and the formation description booklets [39,40,41], as there were no medium-scale thematic maps available. The following elements were searched and written in the database in case of each geological formation polygon: (1) precious stones—e.g., diamond, gold, silver; (2) metallic minerals—e.g., lead, copper, manganese; (3) industrial minerals—e.g., quartz sand, calcite, bauxite, gypsum; (4) geological energy sources—e.g., lignite, oil shale, natural gas; (5) mineral waters and springs. In each cell, we counted all the different elements (enumerated before) and considered their number as the sub-index. The maximum number of different mineral elements in a cell was 5 and the minimum was 0. The average number of mineral elements in the cells was 0.68. The values were also normalised.

4. Results

Each sub-index value set was visualised using the 2 × 2 km grid in QGIS. Five cluster groups were formed by ‘Natural Breaks’ method [44]—it minimises normal deviation between the elements of groups and maximises the difference between the groups. The distribution of geodiversity elements can be observed in the produced maps (Figure 7).
The visualization of the geological sub-index shows that Balaton Uplands (south-eastern and central part—for geographic names see Figure 1) is geologically the most diverse area in the Geopark. There are some large-score cells in the Bakony Mountains (northern part), Tapolca Basin and the Keszthely Mountains (western part). The shore of Lake Balaton and the valley of stream Séd around Veszprém are the least diverse areas in the Geopark in this aspect (Figure 7a). Before normalisation, the maximum geological sub-index value was 21 and the minimum 0.
The highest pedological sub-index values are in the eastern part of the Geopark (around the city of Várpalota), Tapolca Basin, and the southern shore of Lake Balaton around Fonyód. The lowest scores occur around Veszprém and in the Bakony Mountains (Figure 7b). The maximum pedological subindex value before normalisation was 5, the minimum 0.
The distribution of geodiversity sub-index is rather even. The balanced relief index resulting from geomorphon analysis is supplemented with the hydrological index. Due to this, the shores of Lake Balaton and watercourses (e.g., Séd stream) can be spotted and tracked in the grid map (Figure 7c). Deep valleys and relatively high altitudes resulted large scores in some cells in the Bakony Mountains, Keszthely Mountains and Balaton Uplands. Before normalisation, the relief index value was between 1 and 10, and the hydrographic index value was between 0 and 3.
Palaeontology slightly correlates with the geological sub-index. Some sedimentary rocks (e.g., various formations of Triassic/Jurassic limestone and dolomite) of Balaton Uplands and even more in the Bakony Mountains have many more kinds of fossils than quaternary sediments or volcanic rocks (Figure 7d). The maximum value of palaeontological sub-index before normalisation was 12, and the minimum 0.
The sub-index of mineral occurrences, mineral waters and springs shows the opposite trend that of palaeontology. Geopark areas with metamorphic and volcanic features show a greater variety of these geodiversity elements: Tapolca Basin, the western part of the Bakony Mountains and the Keszthely Mountains are outstanding in this aspect (Figure 7e). The maximum value of this sub-index before normalisation was 5, and the minimum 0.

Geodiversity Index Map

By summarising normalised sub-indices in each cell, geodiversity index was produced. It was also visualised on a grid map (Figure 7f). The highest geodiversity values occur in the Balaton Uplands, Tapolca Basin, Káli Basin and Bakony Mountains (Figure 8). Less diverse are the valley of Séd stream and in the southern and eastern surroundings of the shores of Lake Balaton.

5. Discussion—The Use of Geomorphons and the Comparison of Raw and Normalised Geodiversity Index Values

In this paper, we presented two new methodological approaches: the use of geomorphons and the normalisation of geodiversity values.
The significance of the geomorphon concept in the production of a relief sub-index is that the use of DEMs eliminates the country- and scale-dependent factor, making the results more easily comparable. The algorithm used to generate the geomorphon parameters ensure that the method produces results that can be applied to any base data scale and grid resolution. This way, the relief sub-index can be produced easier where DEMs are available, thus simplifying the geospatial processing. However, for methods where other geomorphological parameters (e.g., age of landforms) are used to determine this sub-index in addition to the classification of landform types, this method may not be appropriate.
The summing up of the five different thematic indices without normalisation results in a geodiversity value in which some of the sub-indices are overrepresented (Table 1). Such assessment may cause a large inequality when taking sub-indices into account. Two different percentage values were counted. ‘Max % of geodiversity’ means the proportion of each sub-index maximum value in the maximum geodiversity index, while ‘Median % of geodiversity’ means the proportion of each sub-index median value in the median geodiversity index. If looking at the median values, it occurs that the geological (33.3% of points) and geomorphological values (42.8% of points) are largely overrepresented in the assessment.
This unbalanced situation is mostly caused by the various scales and different types (vector, raster, DEM, formation descriptions) of base data. Accordingly, this kind of error can occur anywhere in the world where datasets from various sources are used for geodiversity assessment. The difference between the two methods is visible in the geodiversity maps generated from the two types (raw and normalised) of indices (Figure 9).
The normalisation reduced the highlighting of both diverse and less diverse cells. Areas of large geodiversity are confirmed by the database of visitors’ centres and visitors’ sites of the Geopark [45]: they are mostly located in the most geologically diverse parts of the Geopark. Figure 10 supports this correlation. The location of the 77 most important geosites can be examined in it: most of them are in ‘highlighted’, diverse areas.
Undoubtedly, the best practice would be to use consistent data from one single source with similar complexity. However, even in countries with well-developed science and policy, this is quite difficult, but researchers working in more difficult areas have a much harder task.

6. Conclusions

In this study, the refined geodiversity assessment of Bakony–Balaton UNESCO Global Geopark was presented. This was the first application of geodiversity quantification models in Hungary. The five thematic sub-indices (geological, pedological, geomorphological, palaeontological, and mineral occurrences) were determined by using digital maps and data sources (Figure 7). Only open-source software solutions were involved: QGIS, GRASS GIS and SAGA GIS.
This research applied a widely used model of the Portuguese school [10], but the calculation of the geomorphological sub-index is different in our case. In the case of geomorphological sub-index, a digital elevation model was used to conduct both relief (Figure 5) and hydrological (Figure 6) evaluation. Although the remaining sub-indices require more country-specific data, the choice of worldwide available data sources makes it easier to compare geodiversity values in different areas.
Our modified model uses the statistical method of normalisation that has two important impacts on the assessment data:
  • It eliminates the overrepresentation of sub-indices, therefore they become equally weighted. This results in more realistic modelling of geodiversity and separating diverse and less diverse areas (Figure 9).
  • Normalisation also helps researchers working with datasets from various sources and with different scales. Larger geodiversity value is generated from more complex layers. This method smoothens the differences between sub-indices without causing any harm to the scientific substance of a certain index.
This kind of quantification holds the possibility to discover highly diverse cells with unexplored values—especially in new areas to be protected. Nevertheless, it is also useful for regions with active nature conservation since new elements and aspects can be found that enrich scientific message and stimulate tourism activities.
Geodiversity assessment is the initial step before establishing any infrastructure of geoheritage and geoconservation. Although this has been done before in the study area using other methods, no numerical survey has yet been carried out prior to this study. It was confirmed that tourist destinations with tourism and nature protection potential related to earth sciences are mainly located in cells with high diversity value. Through the methods presented, this paper helps researchers to conduct a thematically balanced assessment. Yet, the real contribution of the presented model to conservation and sustainable development needs to be tested in other areas.

Author Contributions

Conceptualization, Márton Pál and Gáspár Albert; methodology, Márton Pál, Gáspár Albert; software, Márton Pál; validation, Márton Pál, Gáspár Albert; formal analysis, Márton Pál, Gáspár Albert; investigation, Márton Pál, Gáspár Albert; resources, MÁRTON PÁL; data curation, Márton Pál; writing—original draft preparation, Márton Pál; writing—review and editing, Gáspár Albert, Márton Pál; visualization, Márton Pál, Gáspár Albert; supervision, Gáspár Albert; project administration, Gáspár Albert; funding acquisition, Gáspár Albert All authors have read and agreed to the published version of the manuscript.

Funding

The study was partly supported by the National Research Development and Innovation (NRDI) Fund of Hungary, financed under the Thematic Excellence Programme TKP2020-NKA-06 (National Challenges Subprogramme) funding scheme.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The resulting geodiversity database (including all subindex values) can be downloaded from http://mercator.elte.hu/~marchello/bbalaton/geodiv.gpkg.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wilson, E.O. The Diversity of Life; Belknap Press: Cambridge, MA, USA, 1992. [Google Scholar]
  2. United Nations. Rio 92. Conferentia de las Naciones Unidas Sobre el Medio Ambiente y Desarrollo; Ministerio de Obras Püblicas y Urbanismo: Madrid, Spain, 1993.
  3. Sharples, C. A Methodology for the Identification of Significant Landforms and Geological Sites for Geoconservation Purposes; Forestry Commission: Hobart, Australia, 1993.
  4. Kiernan, K. Conserving Geodiversity and Geoheritage: The Conservation of Glacial Landforms; Report to the Australian Heritage Commission; Forest Practices Unit: Hobart, Australia, 1996.
  5. Kiernan, K. The Conservation of Landforms of Coastal Origin: Conserving Tasmania’s Geodiversity and Geoheritage; Forest Practices Unit: Hobart, Australia, 1997; p. 273.
  6. Gray, M. Geodiversity: Valuing and Conserving Abiotic Nature; J. Wiley: New York, NY, USA, 2004; ISBN 0470848952. [Google Scholar]
  7. Gray, M. Geodiversity: The origin and evolution of a paradigm. In The History of Geoconservation; Burek, C.D., Prosser, C.D., Eds.; Special Publication 300; The Geological Society: London, UK, 2008; pp. 31–36. [Google Scholar]
  8. Gray, M. Geodiversity: The backbone of geoheritage and geoconservation. In Geoheritage, 1st ed.; Reynard, E., Brilha, J., Eds.; Assessment, Protection and Management; Elsevier: Amsterdam, The Netherlands, 2018; pp. 13–25. [Google Scholar] [CrossRef]
  9. Henriques, M.H.; Pena dos Reis, R.; Brilha, J.; Mota, T.S. Geoconservation as an emerging geoscience. Geoheritage 2011, 3, 117–128. [Google Scholar] [CrossRef] [Green Version]
  10. Pereira, D.I.; Pereira, P.; Brilha, J.; Santos, L. Geodiversity Assessment of Parana’ State (Brazil): An Innovative Approach. Environ. Manag. 2013, 52, 541–552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Cendrero, A. Propuesta sobre criterios para la clasificación y catalogación dei Patrimonio Geológico. In El Patrimonio Geolögico. Bases Para su Valoración, Protection, Conservation y Utilization; Monografias dei Ministerio de Obras Públicas, Transportes y Medio Ambiente: Madrid, Spain, 1996; pp. 29–38. [Google Scholar]
  12. Durán, J.J.; Brusi, D.; Palli, L.; Lopez, J.; Palacio, J.; Vallejo, M. Geologja ecológica, geodiversidad, geoconservation y patrimonio geológico: La Declaración de Girona. In Proceedings of the Comunicaciones de la IV Reunión Nacional de Patrimonio Geológico, Madrid, Spain, June–July 1998; Durán, J.J., Vallejo, M., Eds.; Sociedad Geolögjca de Espana: Madrid, Spain, 1998; pp. 69–72. [Google Scholar]
  13. Kozlowski, S. Geodiversity: The concept and scope of geodiversity. Prz. Geol. 2004, 52, 833–837. [Google Scholar]
  14. Serrano, E.; Ruiz-Flaño, P. Gedoiversity: A theoretical and applied concept. Geogr. Helv. 2007, 62, 140–147. [Google Scholar] [CrossRef]
  15. De Paula Silva, J.; Rodrigues, C.; Pereira, D.I. Mapping and Analysis of Geodiversity Indices in the Xingu River Basin, Amazonia, Brazil. Geoheritage 2015, 7, 337–350. [Google Scholar] [CrossRef] [Green Version]
  16. Pérez-Umaña, D.; Quesada-Román, A.; Tefogoum, G.Z. Geomorphological heritage inventory of Irazú Volcano, Costa Rica. Int. J. Geoherit. Parks 2020, 8, 31–47. [Google Scholar] [CrossRef]
  17. Tefogoum, G.Z.; Quesada-Román, A.; Pérez-Umaña, D. Geomorphosites inventory in the Eboga Volcano (Cameroon): Contribution for geotourism promotion. Géomorphologie: Relief Process. Environ. 2020, 26, 19–33. [Google Scholar] [CrossRef] [Green Version]
  18. Quesada-Román, A.; Pérez-Umaña, D. Tropical paleoglacial geoheritage inventory for geotourism management of Chirripó National Park, Costa Rica. Geoheritage 2020, 12, 58. [Google Scholar] [CrossRef]
  19. Pérez-Umaña, D.; Quesada-Román, A.; De Jesús-Rojas, J.C.; Zamorano-Orozco, J.J.; Dóniz-Páez, J.; Becerra-Ramírez, R. Comparative analysis of geomorphosites in volcanoes of Costa Rica, Mexico, and Spain. Geoheritage 2018, 11, 545–559. [Google Scholar] [CrossRef]
  20. Quesada-Román, A.; Tefogoum, G.Z.; Pérez-Umaña, D. Geomorphosites comparative analysis in Costa Rica and Cameroon volcanoes. Geoheritage 2020, 12, 90. [Google Scholar] [CrossRef]
  21. Quesada-Román, A.; Pérez-Umaña, D. State of the Art of Geodiversity, Geoconservation, and Geotourism in Costa Rica. Geosciences 2020, 10, 211. [Google Scholar] [CrossRef]
  22. Gyalog, L.; Síkhegyi, F. (Eds.) Geological Map of Hungary (1:100,000); Geological Institute of Hungary: Budapest, Hungary, 2005. Available online: https://map.mbfsz.gov.hu/fdt100/ (accessed on 30 January 2021).
  23. Dövényi, Z. (Ed.) Geography of the Carpathian Basin (A Kárpát-Medence Földrajza); Academic Press: Budapest, Hungary, 2012. [Google Scholar]
  24. Károly, K. (Ed.) National Atlas of Hungary: Natural Environment; Hungarian Academy of Sciences, Research Centre for Astronomy and Earth Sciences: Budapest, Hungary, 2018. [Google Scholar]
  25. Kocsis, K.; Schweitzer, F. (Eds.) Hungary in Maps; Geographical Research Institute, Hungarian Academy of Sciences: Budapest, Hungary, 2009; pp. 1–211. [Google Scholar]
  26. Bakony–Balaton Geopark. Geological Heritage of the Bakony-Balaton Geopark. 2012. Available online: http://www.geopark.hu/en/home/bakony-balaton-geopark/geological-heritage (accessed on 16 November 2020).
  27. Horváth, G.; Lóczy, D. Geoheritage, geoconservation, geomorphosites in Hungary. In Landscapes and Landforms of Hungary; Lóczy, D., Ed.; Springer: New York, NY, USA, 2015; pp. 281–288. [Google Scholar] [CrossRef]
  28. Yamazaki, D.; Ikeshima, D.; Tawatari, R.; Yamaguchi, T.; O’Loughlin, F.; Neal, J.C.; Sampson, C.C.; Kanae, S.; Bates, P.D. A high accuracy map of global terrain elevations. Geophys. Res. Lett. 2017, 44, 5844–5853. [Google Scholar] [CrossRef] [Green Version]
  29. Lóczy, D. Geomorphological regions. In Landscapes and Landforms of Hungary; Lóczy, D., Ed.; Springer: New York, NY, USA, 2015; pp. 281–288. [Google Scholar] [CrossRef]
  30. Prina, A.O. Taxonomic review of the genus “Crambe” sect. “Crambe” (Brassicaceae, Brassiceae). In Anales del Jardín Botánico de Madrid; Real Jardín Botánico: Madrid, Spain, 2009; Volume 66, pp. 7–24. [Google Scholar] [CrossRef] [Green Version]
  31. Budai, T.; Császár, G.; Csillag, G.; Dudko, A.; Koloszár, L.; Majoros, G. A Balaton-Felvidék Földtana—Magyarázó a Balaton-Felvidék Földtani Térképéhez, 1:50 000; Geological Institute of Hungary: Budapest, Hungary, 1999. [Google Scholar]
  32. Barczi, A.; Grónás, V.; Szalai, D.Á.; Nagy, V. Assessing landscape changes on Tihany peninsula to promote sustainable land use. In XVI. Nemzetközi Tudományos Napok: A Tudományos Napok Publikációi: Fenntarthatósági Kihívások és Válaszok; EKE Líceum Press: Eger, Hungary, 2018; pp. 271–280. [Google Scholar]
  33. Harangi, S. Volcanic heritage of the Carpathian—Pannonian region in eastern-Central Europe. In Volcanic Tourist Destinations; Erfurt-Cooper, P., Ed.; Springer: New York, NY, USA, 2014; pp. 103–123. [Google Scholar] [CrossRef]
  34. Veress, M.; Tóth, G. Kál Basin: Pseudokarren on Sandstone. In Landscapes and Landforms of Hungary; Lóczy, D., Ed.; Springer: New York, NY, USA, 2015; pp. 79–87. [Google Scholar] [CrossRef]
  35. Móga, J.; Strat, D.; Iván, V.; Mari, L.; Kiss, K.; Szabó, M.; Csüllög, G. Changes of the karst landscape and epikarst system in the area of the Tapolca karst terrains, North-West Balaton Highlands, Hungary. In Forum Geografic; University of Craiova, Department of Geography: Craiova, Romania, 2017; Volume 16, p. 12. [Google Scholar]
  36. Balaton Uplands National Park. 2020. Available online: https://www.bfnp.hu/en/latogatohelyek (accessed on 16 November 2020).
  37. Araujo, A.M.; Pereira, D.I. A New Methodological Contribution for the Geodiversity Assessment: Applicability to Ceará State (Brazil). Geoheritage 2018, 10, 591–605. [Google Scholar] [CrossRef]
  38. ATK. TAKI (Institute for Soil Sciences and Agricultural Chemistry): Az Agrotopográfiai Adatbázisról. 2020. Available online: https://www.mta-taki.hu/hu/osztalyok/kornyezetinformatikai-osztaly/agrotopo (accessed on 17 November 2020).
  39. Gyalog, L. (Ed.) A Földtani Térképek Jelkulcsa és a Rétegtani Egységek Rövid Leírása; Geological Institute of Hungary: Budapest, Hungary, 1996. [Google Scholar]
  40. Császár, G. Basic Lithostratigraphic Units of Hungary; Charts and short descriptions; Geological Institute of Hungary: Budapest, Hungary, 1997. [Google Scholar]
  41. Gyalog, L. (Ed.) Magyarázó Magyarország Fedett Földtani Térképéhez (az Egységek Rövid Leírása); Geological Institute of Hungary: Budapest, Hungary, 2005. [Google Scholar]
  42. Jasiewicz, J.; Stepinski, T. Geomorphons—A pattern recognition approach to classification and mapping of landforms. Geomorphology 2013, 182, 147–156. [Google Scholar] [CrossRef]
  43. Strahler, A.N. Quantitative analysis of watershed geomorphology. Eos Trans. Am. Geophys. Union 1957, 38, 913–920. [Google Scholar] [CrossRef] [Green Version]
  44. Jenks, G.F. The Data Model Concept in Statistical Mapping. Int. Yearb. Cartogr. 1967, 7, 186–190. [Google Scholar]
  45. Bakony–Balaton Geopark. Nature Trails and Guided Geotours. 2020. Available online: http://geopark.hu/en/home/bakony-balaton-geopark/geotourism (accessed on 18 November 2020).
Figure 1. The geographical setting of Bakony–Balaton UGGp (red border). The green rectangles mark the sample areas of Figures 3–6. DEM source: [28].
Figure 1. The geographical setting of Bakony–Balaton UGGp (red border). The green rectangles mark the sample areas of Figures 3–6. DEM source: [28].
Ijgi 10 00566 g001
Figure 2. The area of Bakony–Balaton UGGp (red border) is divided into 949 cells, each with an area of 2 × 2 km2. Sub-indices and the geodiversity index are calculated for all cells. The green rectangles mark the sample areas of Figure 3, Figure 4, Figure 5 and Figure 6.
Figure 2. The area of Bakony–Balaton UGGp (red border) is divided into 949 cells, each with an area of 2 × 2 km2. Sub-indices and the geodiversity index are calculated for all cells. The green rectangles mark the sample areas of Figure 3, Figure 4, Figure 5 and Figure 6.
Ijgi 10 00566 g002
Figure 3. Geological sub-indices (first line: SNDXcell, second line: SNDXnorm,cell) in Tapolca Basin, around the volcanic witness mountains.
Figure 3. Geological sub-indices (first line: SNDXcell, second line: SNDXnorm,cell) in Tapolca Basin, around the volcanic witness mountains.
Ijgi 10 00566 g003
Figure 4. Pedological sub-indices (first line: SNDXcell, second line: SNDXnorm,cell) in Tapolca Basin, around the volcanic witness mountains.
Figure 4. Pedological sub-indices (first line: SNDXcell, second line: SNDXnorm,cell) in Tapolca Basin, around the volcanic witness mountains.
Ijgi 10 00566 g004
Figure 5. Geomorphons around Zirc (Northern Bakony). The base data is the MERIT DEM, the analysis was performed with GRASS GIS.
Figure 5. Geomorphons around Zirc (Northern Bakony). The base data is the MERIT DEM, the analysis was performed with GRASS GIS.
Ijgi 10 00566 g005
Figure 6. Hydrological index determined from the generated watercourse network, categorised using Strahler’s order (around Herend). The small digits are the stream hierarchical levels while the big ones are the hydrological indices of the corresponding cells.
Figure 6. Hydrological index determined from the generated watercourse network, categorised using Strahler’s order (around Herend). The small digits are the stream hierarchical levels while the big ones are the hydrological indices of the corresponding cells.
Ijgi 10 00566 g006
Figure 7. Index values visualised in QGIS: (a) geological, (b) pedological, (c) geomorphological (relief + hydrography), (d) palaeontological, (e) mineral occurrences (with mineral waters and springs), and (f) normalised geodiversity index (by summing up the 5 sub-index).
Figure 7. Index values visualised in QGIS: (a) geological, (b) pedological, (c) geomorphological (relief + hydrography), (d) palaeontological, (e) mineral occurrences (with mineral waters and springs), and (f) normalised geodiversity index (by summing up the 5 sub-index).
Ijgi 10 00566 g007
Figure 8. Some examples of geosites from high-geodiversity areas: (a) Koloska Cliffs in the Balaton Uplands near Balatonfüred; (b) a key section in the Tihany Peninsula near Sajkod; (c) panorama of Lake Balaton from the lookout tower of the Mencshely Volcano Educational Trail and (d) panorama with lavender and vineyards from Aszófő.
Figure 8. Some examples of geosites from high-geodiversity areas: (a) Koloska Cliffs in the Balaton Uplands near Balatonfüred; (b) a key section in the Tihany Peninsula near Sajkod; (c) panorama of Lake Balaton from the lookout tower of the Mencshely Volcano Educational Trail and (d) panorama with lavender and vineyards from Aszófő.
Ijgi 10 00566 g008
Figure 9. Comparison of the two geodiversity index maps. The most significant differences (Keszthely Mountains, Tapolca Basin and Séd valley) are marked with circles/ellipsoids. Categories on both maps were created using the natural breaks method.
Figure 9. Comparison of the two geodiversity index maps. The most significant differences (Keszthely Mountains, Tapolca Basin and Séd valley) are marked with circles/ellipsoids. Categories on both maps were created using the natural breaks method.
Ijgi 10 00566 g009
Figure 10. The spatial distribution of the 77 most important and visited geosites of the Geopark (red dots). Most of them are located in cells with high geodiversity values.
Figure 10. The spatial distribution of the 77 most important and visited geosites of the Geopark (red dots). Most of them are located in cells with high geodiversity values.
Ijgi 10 00566 g010
Table 1. Inequality between the sub-indices: the ‘Minimum’ and ‘Maximum’ columns show the smallest and largest values of sub-indices and the overall geodiversity index occurring in the cells without normalisation. ‘Max. % of geodiversity’ and ‘Median % of geodiversity’ show the possible proportion of sub-index values in the final value of geodiversity index without normalisation.
Table 1. Inequality between the sub-indices: the ‘Minimum’ and ‘Maximum’ columns show the smallest and largest values of sub-indices and the overall geodiversity index occurring in the cells without normalisation. ‘Max. % of geodiversity’ and ‘Median % of geodiversity’ show the possible proportion of sub-index values in the final value of geodiversity index without normalisation.
Indices (without Normalisation)MinimumMaximumMedianMax. % of GeodiversityMedian % of Geodiversity
Geological021748.8%33.3%
Pedological05211.6%9.5%
Geomorphological113930.2%42.8%
Palaeontological012127.9%4.8%
Mineral occurrences05111.6%4.8%
Geodiversity14321100%100%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pál, M.; Albert, G. Refinement Proposals for Geodiversity Assessment—A Case Study in the Bakony–Balaton UNESCO Global Geopark, Hungary. ISPRS Int. J. Geo-Inf. 2021, 10, 566. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10080566

AMA Style

Pál M, Albert G. Refinement Proposals for Geodiversity Assessment—A Case Study in the Bakony–Balaton UNESCO Global Geopark, Hungary. ISPRS International Journal of Geo-Information. 2021; 10(8):566. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10080566

Chicago/Turabian Style

Pál, Márton, and Gáspár Albert. 2021. "Refinement Proposals for Geodiversity Assessment—A Case Study in the Bakony–Balaton UNESCO Global Geopark, Hungary" ISPRS International Journal of Geo-Information 10, no. 8: 566. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10080566

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop