5adc42b3-104c-49cf-8fde-f24999ad217e

Geoconservation Research

Volcanism & Geoconservation

2023, Volume 6 / Issue 1 / pages(1-17)

Cretaceous Volcanic Rock Geosites of the Papuk UNESCO Global Geopark (Croatia): Scientific Aspect of Geoheritage in Geoeducation,

Geotourism and Geoconservation

Dražen Balen1, Petra Schneider1,*, Zorica Petrinec1, Goran Radonić2, Goran Pavić2

1Department of Geology, Faculty of Science, University of Zagreb, Horvatovac 95, 10000 Zagreb, Croatia

2Papuk Geopark / Papuk Nature Park, Trg Gospe Voćinske bb, 33522 Voćin, Croatia

Corresponding Author:

Petra Schneider Department of Geology Faculty of Science University of Zagreb, Croatia

Research over decades confirms the geological values of the Papuk UNESCO Global Geopark (Croatia) as a unique place in the regional frame where several orogenic events left their traces through the formation of diverse lithologies. The important part of the geological mosaic, at least in the western part of the Geopark, is the variety of igneous (sub)volcanic rocks. Albite rhyolite at Rupnica and Trešnjevica geosites formed in the Late Cretaceous (~81 Ma), recording the geological event(s) associated with the closure of the Neotethys Ocean. At that time, acidic silicate melt rose fast from the deep crustal levels to the near surface, where cooling caused regular cracking and the development of columnar jointing. Today, these geosites attract the attention of visitors and therefore they are important landmarks that contribute to local (geo)tourism. They are also used as educational sites for both higher education and schoolchildren with Rupko’s Geological School, in which the development of columnar jointing is explained popularly, further enhancing public awareness of the geodiversity and geoheritage of the Mt. Papuk area. The recently opened Geo-info Center in Voćin significantly enhances the geoheritage presentation at the Geopark.

Keywords: Papuk UNESCO Global Geopark, Rupnica, Trešnjevica, Albite rhyolite, Columnar jointing

Introduction

Papuk Geopark (Croatia) joined the Europe- an Geoparks Network and UNESCO Global Geoparks Network in 2007. It encompasses the most prominent mountain in the Slavonian Moun- tains, Mt. Papuk. The Slavonian Mountains con- sist of several mountains or hills up to 1000 m

In the northern part of the Slavonian Mountains (including Mt. Papuk), metamorphic and igneous rocks are well exposed. They represent one of the surface exposures of the crystalline basement of

Figure 1. A) Location of the Papuk UNESCO Global Geopark. Upper left inset with simplified map of the Dinaride-Al- pine-Pannonian region shows position of the mountain belts, Pannonian Basin and Tisia Mega-Unit (blue shaded area). The position of Slavonian Mts. is marked by blue rectangle. B) Simplified geological map of Croatia M 1:300.000 (Croatian Geological Survey, 2009) with positions of the key volcanic localities (Rupnica and Trešnjevica). Green line marks the bor- der of Papuk Nature Park, while gray line marks the border of Papuk UNESCO Global Geopark. Legend: 1 – metamorphic rock complex (Precambrian); 2 – progressive metamorphic series (Ordovician-Devonian); 3 – metamorphic rock complex (Ordovician-Devonian); 4 – granitic rocks (Ordovician-Devonian); 8 – Variscan/Hercynian semi-metamorphic complex (De- vonian-Permian); 13b – clastic sedimentary deposits (Late Permian); 14 – clastic sedimentary deposits and carbonates (Early Triassic); 15 – carbonate deposits (Middle Triassic); 44 – clastic sedimentary deposits and carbonates (Miocene, M2,3); 46 – “litavac” (bioclastic beech carbonate) and clastic sedimentary deposits with volcanic rocks (Miocene, M4); 47 – carbonate clas- tic deposits (Miocene, M5,6); 48 – clastic sedimentary rocks and coal (Miocene, M7); 50 – “Paludina beds” (freshwater deposits with Viviparus sp.) (Pliocene); 52 – clastic deposites (Plio-Quaternary); 54a – loess (Pleistocene); 58a – delluvial-proluvial deposits (Holocene); 58b – alluvial deposits (Holocene); α – andesites and rhyolites (Miocene, M3,4).

the Pannonian Basin, which belongs to the Tisia Mega-Unit (e.g., Pamić et al. 2002a, b; in the sci- entific literature also known as Tisza or Tisa). The Tisia Mega-Unit dominates the geology of Mt. Papuk and contributes to the complexity and great diversity of rocks, including granitoids and meta- morphic complexes (Jamičić 1983; Pamić & Lan- phere 1991; Balen et al. 2006, 2013; Horváth et al. 2010) developed through several orogenic cycles of pre-Variscan (Early Paleozoic), Variscan (Late Paleozoic) and Alpine (Cretaceous) ages (Balen et al. 2015, 2018, Biševac et al. 2010, Horvat et al. 2018, Slovenec et al. 2020) and is covered by the Late Permian–Jurassic sedimentary rocks (Jamičić 1983, 1988; Pamić et al. 1988; Pamić & Lanphere 1991). The Slavonian Mountains rise above the lowlands between the Sava and Drava rivers and are thus often referred to as inselbergs (island hills), namely mountains that rise from relatively flat surroundings through tectonics (set of normal faults) and erosion. Therefore, the crystalline basement is exposed in the Slavonian Mountains, surrounded by the Neogene and Qua- ternary sedimentary cover of the Pannonian Basin (Jamičić 1983, 1988; Pamić et al. 1988; Pamić &

Among the numerous geological features of Mt. Papuk, which developed over a vast period of

~500 Ma from the Early Paleozoic (metamorphic rocks) to recent times (on-going formation of tufa rock), the (most) important part of the geological mosaic, at least in western parts, are the Late Cre- taceous volcanic rocks (Schneider et al. 2022). Although volcanic rocks cover only a relatively small area in the western part of Mt. Papuk (~2 % of the surface area), their occurrence and petroge- netic characteristics are important for interpreta- tion of the complex geodynamic evolution of this marginal part of the Pannonian Basin. The best ex- posures of the volcanic rocks are in the vicinity of the town of Voćin (northern slopes of Mt. Papuk; inset in Fig. 1) and are part of the larger volcanic body, Voćin volcanic mass, that covers a total of

of Late Cretaceous age (~81 Ma for the Rupnica geosite, Schneider et al. 2022) that cross cut the Variscan basement rocks (Jamičić & Brkić 1986; Jamičić et al. 1987, 1989; Jamičić 1989; Pamić 1991, 1997). Although the rocks of the Voćin vol- canic mass have limited surface exposure, their representatives (the best exposure of volcanic rocks in the Croatian part of the Pannonian Basin), rarity (columnar jointing phenomenon occurred in the acidic rock), integrity (well-preserved sites), scientific knowledge (scientific publications that emphasize Cretaceous age, high-T conditions, geochemistry corresponding to anorogenic mag- ma, rapid uplift of magma from deeper crustal lev- els, transition from compressional to extensional tectonics) give an opportunity for reconstructions on a regional scale. That is to say, similar intense magmatic activity during the Late Cretaceous also occurred in the adjacent geotectonic unit called the Sava-Vardar Zone (Pamić 1993, 2002) or Sava Zone (Schmid et al. 2008, 2020), making the Cre- taceous volcanic rocks of the western part of Mt. Papuk an important part of the still fragmented geological record of the geodynamic evolution in this region. Decades of geological research in the western part of Mt. Papuk have provided a solid scientific base for the reconstruction of its com- plex geologic evolution, highlighting the current scientific importance of the Papuk Geopark geo- logical inventory.

Papuk UNESCO Global Geopark and its Key Geosites

The area of Mt. Papuk is characterized by rich bio- diversity and historical heritage, but also by ex-

ceptional geodiversity, providing the grounds for the declaration of the Papuk Nature Park (Papuk NP) in 1999, and defining Javna ustanova Park prirode Papuk (Public Institution Papuk Nature Park) as the institution that governs the park, covering an area of 336 km2 (Official Gazette of the Republic of Croatia, NN 45/99 and 96/99). In 2007, the geodiversity and geoheritage of Mt. Papuk were internationally recognized, making it the first protected area in Croatia to receive the status of a UGGp (the area extended to 524 km2) and it became a member of both the European and Global Geopark Networks under the auspices of UNESCO (Division of Earth Sciences).

From a geological point of view, Mt. Papuk is the most diverse mountain in Croatia. Geological di- versity is represented by both geochronological diversity (rocks formed from the early Palaeozoic to the present) and lithological diversity (various sedimentary, metamorphic and igneous rocks). Igneous rocks of the Voćin volcanic mass are presented to the public in two out of 69 geosites of the Papuk UGGp – Rupnica and Trešnjevica. The main features of other geosites are related to tectonic, geomorphological, sedimentological, paleontological, stratigraphic and hydrogeologi- cal features (Public Institution Papuk Nature Park 2011).

The unique geological features of the Rupnica geosite, the most famous geologic symbol of the Papuk UGGp, were well exposed through mining, mainly related to the construction of railway and road infrastructure in the 19th and 20th centuries. The area of the former quarry on the northwestern flank of Mt. Papuk is easily accessible to visitors, and it is the aesthetic component that primarily attracts the visitors’ attention – the morpholog- ical feature of the columnar jointing of volcanic rocks, popularly called organ pipes. However, the presentation of the geosite as a locality where one can observe the effects of the “frozen” fiery forces of nature has promoted its geological significance to the public over the years, making it a landmark of the geological heritage of Papuk UGGp. The nearby abandoned Trešnjevica quarry is another geosite related to the same volcanic episode. Here

a swarm of effusive rocks of basic to acidic com- position is exposed, cross cutting the older igne- ous and metamorphic rocks.

Rupnica Geosite

The more or less regular prismatic fracturing of volcanic rocks into columns has been observed since ancient times and reported from many world-renowned localities. Such descriptions in the literature date back to the 17th century, re- porting the Giant’s Causeway in Northern Ireland (Bulkeley 1693, cited in Tomkeieff 1940). Basal- tic columnar jointing is common in Central and Eastern Europe and superb exposures are known from Austria, Bohemia, Hungary, Poland, Roma- nia, and Slovakia. Most of them form the core of geosites and protected areas, and some of them are even part of the UNESCO Global Geoparks Net- work, such as the Bakony-Balaton Geopark (Mar- tin & Németh 2004) and the Novohrad-Nógrád Geopark (Hetenyi et al. 2012).

Although the exposure of columnar jointing at Rupnica geosite (Fig. 2) is not as impressive in size as the Giant’s Causeway or some other noted sites with the same feature (e.g., Devil’s Postpile; Columbia River Basalts, USA), it is a scientifical- ly significant geosite because of the type of rock in which the columnar jointing developed. Gen- erally, columnar jointing develops readily in less viscous basaltic magmas, as seen for example at the type locality of basalt, the Stolpen Castle Hill (recently corrected to be basanite to nephelinite), which shows basaltic columnar jointing (Tietz et al. 2018). However, columnar jointing is today known to develop in almost all rock types includ- ing welded ignimbrites (Ahn 2018, Lebti et al. 2006), and in various geological environments, including even glaciovolcanic (Tuffen et al. 2001) and subaqueous conditions (Mueller & Donaldson 1992).

The rocks exposed at the Rupnica are acidic and classify as albite and aegirine-albite rhyolites (Ta- jder 1956, 1960). The geodynamics of the origin and emplacement of acidic rocks in the western area of Mt. Papuk (e.g., Balen & Petrinec 2014,

Schneider et al. 2022) place the Rupnica geosite alongside the best columnar rhyolitic rock forma- tions in the world, such as the Early Cretaceous High Island Formation of the Hong Kong UNES- CO Global Geopark (Zhao et al. 2017).

The rocks of Rupnica have a porphyritic tex- ture dominated by albite phenocrysts hosted in a matrix of albite microcrysts, quartz and devit- rified volcanic glass. Accessory minerals include magnetite, apatite, zircon and occasionally alka- li clinopyroxene (aegirine, aegirine-augite). The development of columnar jointing in these rocks, as well as the morphometric characteristics of the columns, were studied by Balen & Petrinec (2014)

using a statistical evaluation of crack counting in which numerous geology and environmental sci- ence students participated. The research revealed that relatively narrow, predominantly four- and five-sided columns, tightly packed striae, and small amplitude wavy sets of columns dominate the exposed part (Figs. 2 and 3). Features quan- tified through the statistical approach, together with petrography, geochemistry, and geology of the geosite surroundings, provided information about the speed of cooling and depth of the (sub) surface acidic volcanic body. Columnar joint- ing developed under the combined influence of large-scale supercooling and thermal contraction. This analysis at Rupnica and at the nearby local-

Figure 2. Panoramic view of the Rupnica geosite. The polygonal columnar pattern developed in the albite rhyolite at Rupnica geosite. The person in the bottom right as a scale.

ities of Omerovac, Jovanovica and Turski grad (our unpublished data), reveals the operation of the “T-junction” formation mechanism that pro- duced predominantly four-sided columns (Fig. 3), indicating rapid cooling at low pressure, with a

near-surface crustal level of emplacement of the subvolcanic body (Balen et al. 2010; Balen 2011; Balen & Petrinec 2014).

ies revealed that the melt originating from mid- to-lower-crustal depths reached the shallow em- placement level as a high-temperature magma (780–900 °C). The source of the parental A-type magma is related to the lower crust in an anorogen- ic tectonic setting with peraluminous, alkali-calcic to alkalic and ferroan character (Balen et al. 2020; Schneider et al. 2022). Recent studies of inclu- sions captured in accessory zircon (i.e., anatase, fluorapatite and hematite) indicate crustal source contribution, rapid uplift and cooling of an oxi- dized magma (Schneider & Balen 2022; Schnei- der et al. 2022). Igneous rocks, now exposed at the Rupnica geosite, formed by cooling and (near) surface crystallization of acidic magma during the Late Cretaceous (~81 Ma, Schneider et al. 2022),

Trešnjevica Geosite

Although Mt. Papuk crystalline is mostly com- posed of pre-Alpine metamorphic and igne- ous rocks, at the Trešnjevica geosite, where the pre-Alpine gneisses and granitoid rocks are ex- posed, swarms of younger igneous rocks can be found (Fig. 4). They crosscut the igneous rocks of the older complex and each other. Many vari- eties of volcanic rocks are found here, including various subtypes of rhyolite, andesite, basalt, me- tabasalt, dolerite (“diabase”), as well as pyroclas- tic rocks and volcanic breccia. The rocks of older complexes, such as granite, pegmatite, gneiss and migmatite, are present at this geosite as well. Geo-

Figure 3. The columnar jointing developed at the Rupnica geosite. A) Columns - side view with propagating cracks and striae. A pen as a scale. B) Polygonal columnar pattern visible at Rupnica geosite, dominated by four-sided columns.

logical relationships are clearly visible and the site thus provides an opportunity to study the different types of igneous rocks, their mineralogy, textures, modes of occurrence, and their original geologi- cal setting based on geochemical data (Balen et al. 2009).

The dominant igneous rock in the dykes on this geosite is rhyolite, which does not differ signifi- cantly from that at the Rupnica geosite. The rhyo- lite varieties include albite, anorthoclase, and ae- girine-albite rhyolite (Tajder 1956, 1960). Other varieties of dyke rocks include andesite which is a dark-gray compact rock (Kišpatić 1887, 1916; Ta- jder 1956). Dark igneous rock is also basalt with

phenocrysts of feldspars (labradorite) up to 4–7 mm in size that can be found in the holocrystal- line matrix together with clinopyroxene (augite), while minerals in the matrix include submillime- ter grains of plagioclase, pyroxene, magnetite, il- menite, apatite and secondary chlorite and calcite. The basalts are porphyritic to ophitic in texture. Besides fresh basalts, there are also altered vari- eties (metabasalts, Pamić 1991), the latter being more common. The texture of the metabasalts is homogeneous and contains up to 20 vol. % of amygdales filled with calcite, chlorite, pumpel- lyite and with minor prehnite, quartz and zeolite. Pyroclastic rocks, including coarse-grained variet- ies and tuffs, cover roughly the same area as effu-

Figure 4. Trešnjevica geosite. A) Panoramic view of the dyke swarm. B) Details of the dyke penetrating the light-coloured rocks of the Variscan igneous and metamorphic complex. Thickness of the dyke ~1 m.

sive rocks (Jamičić et al. 1987; Pamić 1991). The pyroclastic rocks contain basaltic and rhyolitic clasts. Tuffs have a massive texture and comprise basalt lithoclasts, with thicknesses varying from the cm- to deca-meter thick intercalations inside pyroclastic rocks. Lithoclasts of basalt, crystalo- clasts of feldspars and clinopyroxene are found inside the tuff, while devitrified volcanic glass is usually present in the matrix. Tuffs are classified as lithoclastic to vitroclastic (Pamić 1991).

The igneous rocks of the Voćin volcanic mass show transitional affinity, ranging from sub-alka- line to alkaline (Fig. 5A). Incompatible trace ele-

ments, trace element ratios, and apparent enrich- ment of light rare earth elements (REE) together with a moderate to low abundance of high field strength elements, a variable abundance of large ion lithophile elements and a negative europium anomaly for acidic varieties, suggest melting of mid- to low-crust material. The distribution of REE in the studied samples (Fig. 5B) supports this conclusion. The high magmatic temperatures are in accordance with the inferred within-plate and volcanic arc tectonic setting. The preliminary age dating of zircon from albite rhyolite using the LA- ICP-MS method revealed a Late Cretaceous age (our unpublished data) and corresponds well with

the age determined for albite rhyolite from Rupni- ca (~81 Ma) and the regional thermal event (Balen et al. 2020, Schneider et al. 2022). The age of the albite rhyolite of Mt. Papuk is coeval with geo- chemically similar igneous rocks from Mt. Požeš- ka Gora (Balen et al. 2020), Mt. Prosara and Mt. Kozara (Ustaszewski et al. 2009, 2010; Cvetković et al. 2014), which are situated within an import- ant Cretaceous suture zone (Sava Zone) after the closure of Neotethys Ocean. These indicate sim- ilar and contemporaneous geological processes that occurred in a long collisional zone between the ancient margins of Europe and Africa.

Host Rocks of the Cretaceous Volcanic Rocks:

Voćin – Zvečevo – Kamenska Section

Although Cretaceous volcanic rocks are the exclu- sive feature of the western part of Mt. Papuk, sig- nificantly older rocks of the Variscan orogenic cy- cle were recognized early as the most widespread and represent a well-documented geodynamic ep- isode in the area (Pamić & Lanphere 1991). The

Variscan rocks are related to geological processes that occurred in the deeper crustal levels, tens of kilometers deep. The origin of these rocks thus differs in terms of age and geological processes from the Cretaceous igneous rocks from the Rup- nica and Trešnjevica geosites.

Outcrops of Variscan granite, migmatite and gneiss can be found along the road Voćin–Zveče- vo–Kamenska. At the Čarugin Kamen geosite, partially melted crustal rocks (migmatites) with typical zebra-like textures are exposed togeth- er with gneisses. Near the geosite, a wide vari- ety of rock textures and structures are observ- able relating to variation in the volume content of key mineral constituents (feldspar, quartz and mica). Granite with its igneous structures forms large irregular bodies that are also present with- in the gneiss-migmatite zone, although they are extensively tectonically disturbed and fragment- ed. They are characterized by a massive texture (visible at the Zvečevo Potočan geosite) and the

Figure 5. A) TAS diagram (Le Bas et al. 1986) illustrating geochemical variation of the Voćin volcanic rocks. Circles represent analyses from the literature (Tajder 1956, 1960; Lugović 1983; Jamičić et al. 1987; Pamić 1991, 1997), while our data are represented by squares for Trešnjevica geosite and triangles for Rupnica geosite.

B) Spider-diagram of rare-earth element abundances normalized on chondrite after Boynton (1984). Gray shaded

mineral assemblage of quartz, feldspar (albite and microcline) with micas (biotite, muscovite). Complex geological relations among various rock types in the western Papuk area are presented on the informative panel at the geo-information point at the Zvečevo site (Fig. 6).

Linking Geoheritahge to Geoeducation, Geotour- ism and Geoconservation

The basic document of the Papuk NP was the Man- agement Plan of the Papuk Nature Park (Public Institution Papuk Nature Park 2011) that, among others, set aims and activities for protection and

preservation of geological and geomorphological diversity and visiting, recreation, education and the marketing of the destination. However, the Study of the management of the number of visitors with an action plan for the management of visits in the Papuk Nature Park (2020–2029) (Public In- stitution Papuk Nature Park 2020) was published as a part of the project Geo Stories of the UNES- CO Geopark. The document contains descriptive

and analytical parts related to visits to the Papuk NP and a visitor management plan with detailed analysis of the structure and the number of visi- tors, visitor monitoring, visitor and interpretation infrastructure, educational and recreation pro- grams, communication and cooperation, together with capacities for visitor management. Accord- ing to this document, geotouristic development of the western area of Mt. Papuk is planned through

Figure 6. Geo-information point Zvečevo by the Voćin – Kamenska road (Western geological profile) with charac- teristic geological section through igneous, volcanic and metamorphic rocks of the area. Left-hand side of profile shows Voćin volcanic mass.

2020–2029, following the establishment of the Geo-info Center in Voćin. Initial visitor numbers to the Geo-info Center in 2022 (18,500 visitors) speak in favor of such predictions.

Besides the scientific value of the previously de- scribed key volcanic sites of the Papuk UGGp, the Rupnica geosite attracts the attention of visitors primarily by its aesthetic value and is therefore an important landmark that contributes to local tourism. The Rupnica geosite is also an import- ant example of the successful revitalization and presentation of mining heritage in the Republic of Croatia as it relates industrial (mining) and natural (geological) heritage. The research conducted by

Briševac et al. (2021) at four Croatian geoheritage sites (including Rupnica) indicated strong and weak points of the revitalization of the geoheritage sites, as well as future opportunities and possible threats to which the natural heritage like Rupni- ca’s columnar jointing might be exposed. Accord- ing to the authors, local authorities in cooperation with local tourist offices have crucial roles in the preservation and management of geological heri- tage, largely assisted by European funding.

Rupnica is a prime example of how the application of geologically informative principles of geotour- ism (Dowling 2011) can be a driving force in its development. Education about basic earth science

principles and their interpretation at the geosite has recently been upgraded by the opening of the Geo-info Center at Voćin, providing the opportu- nity for visitors to develop their knowledge about natural processes and awareness about the need for protection and management of natural and historic heritage of the region and the role they may take in it. There is a short education trail on the site with four explanatory educational panels showing dif- ferent types of volcanos with a general description of volcanic rocks and processes, and three view- ing platforms for observing the jointing phenome- non from different perspectives. One of the panels is dedicated to the columnar jointing developed at the site and includes a short note about the history of protection of the site. A guided tour around the Rupnica geosite is an optional feature of some of the educational programs. The geosite is located on the road leading along the Western geological profile and is easily accessible by car or bus. A few kilometers further on is the Trešnjevica geo- site with three explanatory education panels about the geological significance of this site and a simple experimental setup simulating a volcanic eruption through the classical chemical reaction between an acid and a carbonate.

The area around the Zvečevo locality was one of the most visited excursion zones on Mt. Papuk, with the most significant accommodation facilities in the area. However, it lost its role as a base from which to visit Papuk after the devastation during the Croatian War of Independence in the 1990s. The site is known for its important biodiversity, with rare wet grassland areas, home to an endan- gered marsh butterfly. In the immediate vicinity of this grassland, the geo-information point Zvečevo displays a panel showing a simplified cross-sec- tion through the Western geological profile (Fig. 6).

Although the Cretaceous volcanic occurrences near Voćin and the section through the oldest rock formations of the western part of Mt. Papuk attract more visitors each year and the presentation of the geo-educational content has improved, significant touristic infrastructure is still lacking in the area. The opening of the modern Geo-info Center in

Voćin in 2022 may thus represent a turning point towards the development of sustainable tourism in the western part of the Papuk UGGp.

Geo-info Center Voćin

In early 2022, a visitor center (Fig. 7) was opened in the town of Voćin, about 1 km from the Rup- nica geosite. It is the main presentation center of the Papuk UGGp. In addition to basic information about Papuk NP, the new Geo-info Center offers educational resources in the fields of geology, seismology, biology and archaeology, all based on scientific research. The educational content presented at the center is highly interactive and takes advantage of modern multimedia technolo- gy. From the main entrance, visitors are directed to the interactive panels (touch screens) that cover a variety of geodiversity topics relating to the UN- ESCO Global Geopark Network. The next theme is presented in a cinema that simulates a journey through the geological evolution of the present Papuk area. Visitors can then explore the large hall with Origin of Life, Wall of Evolution and Solar system exhibitions, along with a round table with a relief and an interactive geological map of Mt. Papuk with rock samples, additional interactive panels, the earthquake table and the Earth globe. On the upper floor, biological, archeological and geological collections are presented, together with a small workshop/laboratory and interactive rooms covering bio- and geo-diversity of the area. There is also specific content dedicated to the ge- ology of the nearby Rupnica and Trešnjevica geo- sites. From a small café-restaurant in the Geo-info Center, one can enjoy the view of the medieval fortress above Voćin (Turski grad fortress), built on a volcanic hill above the confluence of the Jo- vanovica and Voćinska Rijeka rivers, which cut through the Voćin volcanic mass revealing nu- merous smaller outcrops of volcanic rocks near the Geo-info Center, such as Jovanovica geosite. The building of the Geo-info Center is designed so that its external towers resemble the four-sided columns of the Rupnica (Fig. 7).

The Geo-info Center Voćin was built as a part of the project Geo Stories of the UNESCO Geopark (https://geoinfocentar.com/geo-price-unes-

co-geoparka/) that was co-financed from the Op- erational Program Competitiveness and Cohesion 2014–2020, submitted to the call for Promotion of sustainable use of nature heritage in national parks and nature parks of the Ministry of Regional De- velopment and European Union Funds. The proj- ect goals were aimed at the development of tour- istic potential in the continental part of Croatia, including improvement in tourist infrastructure and visitor management system of the Papuk NP

as a part of promoting sustainable socio-econom- ic development. The project holder was the Pub- lic Institution Papuk Nature Park, and the project partner was the County Road Administration of the Virovitica-Podravina County.

Educational Programme Rupko’s Geological School

The Rupnica and Trešnjevica geosites are valuable resources for introducing the public to geology.

Figure 7. Geo-info Centre Voćin. The towers on the right-hand side of the photo symbolically represent Rupnica geosite columns. (Photo: K. Toplak)

They illustrate the diverse and interesting geology of the Papuk UGGp with special attention devot- ed to the volcanic phenomena. In addition, these sites are also used as educational sites for both higher education students and schoolchildren. The remarkable example of the variability of volcanic rocks and processes is explained with educational panels on site and with Rupko’s Geological School educational content at both geosites. Rupko’s Geo- logical School is a Papuk UGGp programme that teaches primary school children the basic concepts of geology using examples from the Papuk area. Papuk UGGp tours are available on request and

are conducted by trained guides. As part of the same programme, a geological guide for children has been published in Croatian and English to ac- company the geological education content.

Simulation of Columnar Jointing Process (Starch-Water Experiment)

The more or less regular prismatic fracturing of volcanic rocks into columns at the Rupnica geo- site and many more famous localities is a univer- sal natural feature. In fact, columnar jointing is not restricted to volcanic rocks, but also develops in other materials. The science that describes the

behavior of different materials is reflected in the various branches of physics, materials science, cryoscience and even planetary science (e.g., Bud- kewitsch & Robin 1994, Goehring 2013).

The Geo-info Center Voćin hosts a small workshop as part of Rupko’s Geological School programme, in which a simple experiment is performed, sim- ulating the development of columnar jointing using non-expensive and everyday ingredients (i.e., water and starch). The regularity behind the processes can be demonstrated and the scientific facts behind the phenomenon of columnar jointing (Toramaru & Matsumoto 2004) can be easily com- municated to students through an experiment with a starch-water mixture (i.e., drying a starch-water slurry or cake in a glass jar). The experiment (Fig. 8A) is easily conducted using an ordinary light- bulb as a heating device. Natural processes in the cooling of the silicate melt can be represented by variations in the speed of drying of the starch-wa- ter slurry, leading to the development of different types of junctions (Fig. 8B). A faster starch-wa- ter mixture drying rate (equivalent to the cooling rate in the melt) causes cracking of the mass with prevailing T-junctions (Fig. 8C). Such cracks run almost perpendicular to each other and form a more or less regular rectangular network. In 3D, the cracks extend from the surface into the depth of the material and form four-sided columns at the base of the experimental volume. On the other hand, a slower drying rate will induce Y-junctions where the cracks form a six-sided pattern and six-sided columns are formed in the starch exper- iment (Fig. 8C). This regularity can be extended to the silicate melt (Budkewitsch & Robin 1994, Goehring & Morris 2005, 2008; Goehring et al. 2006, 2009; Jagla & Rojo 2002; Weinberger 2001; Calderón & Sánchez 2023), which cools faster (favoring T-junctions) or slower (producing more Y-junctions).

Statistical evaluation of columnar jointing at the Rupnica geosite (Balen & Petrinec 2014) pro- vides the basis for interpreting the observed and measured cracking phenomena, according to the theory of rapid melt uplift and cooling. Most com- monly observed are T-junctions and four-sided

columns (Figs. 3B and 8B). This simple experi- mental setup thus provides an excellent example of a large-scale natural phenomenon developed at the nearby Rupnica geosite but also at Trešnjevica geosite.

Summary

The main tasks of a geopark are geoconservation with geoheritage protection, geoeducation and geotourism. Such tasks and related activities in the Papuk UGGp are performed in close collaboration with the local community, government offices (Ministry of Economy and Sustainable Develop- ment), European and Global Geopark Networks but also with scientific community. The develop- ment of infrastructure and increased availability of information about 69 documented geosites, as well as educational trails, numerous workshops, educational programs and facilities at geo-info centers are all in progress in the Papuk UGGp area. The Geo-info Center Voćin opened in 2022 has improved the educational offer. A key part of the geological mosaic, at least in the western parts of Papuk UGGp, is the Cretaceous volcanic area of the Voćin volcanic mass, representing a fun- damental part of the geopark and its geodiversity with Rupnica and Trešnjevica as key geosites.

The acidic silicate melts of the Voćin volcanic mass in the Papuk UGGp originated from mid- to-lower-crustal depths and reached a shallow emplacement level as a high temperature (up to 900 °C) magma. Geochemistry points to a paren- tal A-type magma formed in an anorogenic tec- tonic setting. Studies of accessory minerals such as zircon and apatite, as well as studies of mineral inclusions captured in accessory zircon (i.e., ana- tase, fluorapatite and hematite) indicate possible contributions from deeper crustal levels, fast uplift and cooling of magma. The volcanic rocks were formed in the Late Cretaceous (~81 Ma) through cooling at low pressure, i.e., a shallow position within the igneous body. The relatively fast cool- ing rate caused the development of T-junctions in the columnar jointing process. This resulted in predominantly four-sided columns that are well exposed at Rupnica geosite.

Figure 8. A) Columnar jointing experiment at the Rupko’s Geological School using starch-water slurry. B) Compari- son with a real volcanic rock outcrop at Rupnica. C) Sketch of four-, five- and six-sided columns with real examples at Rupnica geosite and sketch of T- and Y-junctions.

The Papuk UGGp systematically encourages re- search activities by providing logistical support in the field, as well as modest financial support under its own funding program. Systematic collection of scientific data over decades and constant valori- zation of the geoheritage, the efforts of the Papuk UGGp staff and the accompanying scientists are rewarded through the development and enhance- ment of geoeducation, geotourism and geoconser- vation efforts of the area as well as with the inter- national recognition.

Conflict of Interest

References

Ahn KS (2018). Formation mechanism of columnar joints at the Sanin Kaigan Geopark in Japan. Journal of the Korean Earth Science Society. 39(6): 575–592. https://doi.org/10.5467/JKESS.2018.39.6.575

Balen D (2011). Evolution of columnar joint polygonal patterns at Rupnica Geosite (Papuk Geopark, Croatia). 10th Alpine workshop “CorseAlp2011”, Abstract vol-

Balen D & Petrinec Z (2014). Development of co- lumnar jointing in albite rhyolite in a rapidly cooling volcanic environment (Rupnica, Papuk Geopark, Cro- atia). Terra Nova. 26: 102–110. https://doi.org/10.1111/ ter.12075

Balen D, Horváth P, Tomljenović B, Finger F, Humer B, Pamić J & Árkai P (2006). A record of pre-Variscan Barrovian regional metamorphism in the eastern part of the Slavonian Mountains (NE Croatia). Mineralogy and Petrology. 87: 143–162. https://doi.org/10.1007/ s00710-006-0120-1

Balen D, Radonić G & Pavić G (2009). Trešnjevica geosite – “volcano traces from the past”. 8th European Geoparks Conference “New Challenges with Geotour- ism” (pp. 175–175). Geopark Naturtejo, Idanha-a-No- va, Portugal.

Balen D, Radonić G & Pavić G (2010). Order and dis- order in columnar joints of Rupnica Geosite (Papuk Geopark). 4th Croatian Geological Congress, Abstracts Book (pp. 137–138). Croatian Geological Survey.

Balen D, Horvath P, Finger F & Starijaš B (2013). Phase equilibrium, geothermobarometric and xenotime age dating constraints on the Alpine metamorphism re- corded in chloritoid schists from the southern part of the Tisia Mega-Unit (Slavonian Mts., NE Croatia). In- ternational Journal of Earth Sciences. 102: 1091–1109. https://doi.org/10.1007/s00531-012-0850-8

Balen D, Massonne H-J & Petrinec Z (2015). Collision related Early Palaeozoic evolution of a crustal frag- ment from the northern Gondwana margin (Slavoni- an Mts., Tisia Mega-Unit, Croatia): Reconstruction of the P-T path, timing and paleotectonic implications. Lithos. 232: 211–228. https://doi.org/10.1016/j.lith- os.2015.07.003

Balen D, Massonne HJ & Lihter I (2018). Alpine metamorphism of low-grade schists from the Slavoni- an Mountains (Croatia): new P-T and geochronological constraints. International Geology Review. 60(3): 288-

Balen D, Schneider P, Massonne H-J, Opitz J, Luptako- va J, Putiš M & Petrinec Z (2020). The Late Cretaceous A-type alkali-feldspar granite from Mt. Požeška Gora (N Croatia): Potential marker of fast magma ascent in the Europe–Adria suture zone. Geologica Charpath- ica. 71(4): 361–381. https://doi.org/10.31577/Geol- Carp.71.4.5

Biševac V, Balogh K, Balen D & Tibljaš D (2010). Alpine (Cretaceous) very low- to low-grade metamor- phism recorded on the illite-muscovite-rich fraction of metasediments from South Tisia (eastern Mt Papuk, Croatia). Geologica Carpathica. 61: 469–481. https:// doi.org/10.2478/v10096-010-0029-9

Boynton WV (1984). Geochemistry of the rare earth elements: meteorite studies. In: Henderson P (ed), Rare earth element geochemistry (pp. 63–114). New York: Elsevier.

Briševac Z, Maričić A & Brkić V (2021). Croatian Geoheritage Sites with the Best-Case Study Analyses Regarding Former Mining and Petroleum Activities. Geoheritage. 13: 95. https://doi.org/10.1007/s12371-

Budkewitsch P & Robin PY (1994). Modelling the evolution of columnar joints. Journal of Volcanology and Geothermal Research. 59: 219–239. https://doi. org/10.1016/0377-0273(94)90092-2

Bulkeley RB (1693). Part of a letter concerning the Gi- ants Causeway in the County of Atrim in Ireland. Philo- sophical Transactions of the Royal Society London. 17: 708–710. (cited from Tomkeieff, 1940).

Calderón CE & Sánchez JJ (2023). Geometría y di- mensiones de estructuras de disyunción columnar en rocas volcánicas de Colombia. Boletín de Geología. 45(1): 37–51. https://doi.org/10.18273/revbol.v45n1- 2023002 (in Spanish with English abstract)

Croatian Geological Survey (2009). Geological Map of the Republic of Croatia. 1:300 000. Croatian Geologi- cal Survey, Department of Geology (in Croatian).

Cvetković V, Šarić K, Grubić A, Cvijić R & Milošević A (2014). The Upper Cretaceous ophiolite of North

Kozara – remnants of an anomalous mid-ocean ridge segment of the Neotethys? Geologica Carpathica. 65(2): 117–130. https://doi.org/10.2478/geoca-2014-0008

Dowling RK (2011). Geotourism’s global growth. Geo- heritage 3(1): 1–13. https://doi.org/10.1007/s12371-

Goehring L & Morris SW (2005). Order and disorder in columnar joints. Europhysics Letters. 69: 739–745. https://doi.org/10.1209/epl/i2004-10408-x

Goehring L & Morris SW (2008). The scaling of co- lumnar joints in basalt. Journal of Geophysical Re- search B. 10203: 113–129.

Goehring L, Lin Z & Morris SW (2006). An experimen- tal investigation of the scaling of columnar joints. Phys- ical Review E. 74: 036115. https://doi.org/10.48550/ arXiv.cond-mat/0606221

Goehring L, Mahadevan L & Morris SW (2009). Nonequilibrium scale selection mechanism for co- lumnar jointing. PNAS. 106(2): 387–392. https://doi. org/10.1073/pnas.0805132106

Hetenyi G, Taisne B, Garel F, Medard E, Bosshard S & Mattsson HB (2012). Scales of columnar jointing in igneous rocks: field measurements and controlling fac- tors. Bulletin of Volcanology. 74(2): 457–482. https:// doi.org/10.1007/s00445-011-0534-4

Horvat M, Klötzli U, Jamičić D, Buda G, Klötzli E & Hauzenberger C (2018). Geochronology of granit- oids from Psunj and Papuk Mts., Croatia. Geochrono- metria. 45(1): 198–210. https://doi.org/10.1515/geo- chr-2015-0099

Horváth P, Balen D, Finger F, Tomljenović B & Krenn E (2010). Contrasting P–T–t paths from the basement of the Tisia Unit (Slavonian Mts., NE Croatia): appli- cation of quantitative phase diagrams and monazite age dating. Lithos. 117: 269–282. https://doi.org/10.1016/j.

Jagla EA & Rojo AG (2002): Sequential fragmentation: The origin of columnar quasihexagonal patterns. Phys- ical Review E. 65: 026203. https://doi.org/10.1103/ PhysRevE.65.026203

Jamičić D (1983). Structural fabric of the metamor- phosed rocks of Mt. Krndija and the eastern part of Mt. Papuk. Geološki vjesnik Zagreb. 36: 51–72. (in Croa- tian)

Jamičić D (1988). Tectonics of the Slavonian Mts. Dis- sertation, University of Zagreb. (in Croatian)

Jamičić D (1989). Basic Geological Map of SFRY. 1:100 000, Daruvar sheet L 33-95: Geološki zavod, Za- greb, Savezni geološki zavod, Beograd. (in Croatian)

Jamičić D & Brkić M (1986). Basic Geological Map of SFRY. 1:100 000, Orahovica sheet L 33-96: Geološki zavod, Zagreb, Savezni geološki zavod, Beograd. (in Croatian)

Jamičić D, Brkić M, Crnko J & Vragović M (1987). Sheet description of the Basic Geological Map of SFRY 1:100 000, Orahovica sheet L 33-96: Geološki zavod, Zagreb, Savezni geološki zavod, Beograd. (in Croatian)

Jamičić D, Vragović M & Matičec D (1989). Sheet de- scription of the Basic Geological Map of SFRY 1:100 000, Daruvar sheet L 33-95: Geološki zavod, Zagreb, Savezni geološki zavod, Beograd. (in Croatian)

Kišpatić M (1887). Voćinski augitni andezit (trahit). Rad JAZU. 83(8): 179–185. (in Croatian)

Kišpatić M (1916). Augitandesit (Trachyt) von Voćin. Izvješća JAZU. 8: 311–312. (in German)

Le Bas MJ, Le Maitre RW, Streckeisen A, Zanettin B (1986). A chemical classification of volcanic rocks based on total alkali – silica diagram. Journal of Pe- trology. 27: 745–750. https://doi.org/10.1093/petrolo- gy/27.3.745

Neogene and Quaternary ignimbrites in the area of Arequipa, Southern Peru: Stratigraphical and petro- logical correlations. Journal of Volcanology and Geo- thermal Research. 154(3–4): 251–275. https://doi. org/10.1016/j.jvolgeores.2006.02.014

Lugović B (1983). Extrusive rocks from the NW part of Mt. Papuk (Croatia, Yugoslavia). Geološki vjesnik. 36: 131–156. (in Croatian, with an English summary)

Martin U & Németh K (2004). Mio/Pliocene phreato- magmatic volcanism in the western Pannonian Basin. Geologica Hungarica, Serie Geologica. 26: 192–192.

Mueller W & Donaldson JA (1992). A Felsic Dyke Swarm Formed under the Sea - the Archean-Hunter- Mine-Group, South-Central Abitibi Belt, Quebec, Can- ada. Bulletin of Volcanology. 54(7): 602–610. https:// doi.org/10.1007/BF00569944

Official Gazette of the Republic of Croatia (1999). Act of Proclamation of Papuk Nature Park NN 45/1999. (in Croatian).https://narodne-novine.nn.hr/clanci/sluzbe- ni/1999_05_45_871.html

Official Gazette of the Republic of Croatia (1999). Regulation on the establishment of the Public In- stitution Nature Park Papuk NN 96/1999. (in Cro- atian).https://narodne-novine.nn.hr/clanci/sluzbe- ni/1999_09_96_1622.html

Pamić J (1991). Upper Cretaceous basaltoid and pyro- clastic rocks from the Voćin volcanic mass on the Pap- uk Mt. (Slavonija, Northern Croatia). Geološki vjesnik. 44: 161–172. (in Croatian, with an English summary)

Pamić J (1993). Eoalpine to Neoalpine magmatic and metamorphic processes in the northwestern Vardar Zone, the easternmost Periadriatic Zone and the south- western Pannonian Basin. Tectonophysics. 109: 273–

Pamić J (1997). Volcanic rocks of the Sava-Drava in- terfluve and Baranja (Croatia). Zagreb: Nafta. (in Cro- atian, with English summary)

Pamić J & Lanphere M (1991). Hercynian granites and metamorphic rocks from the Papuk, Psunj, Krndi- ja and the surrounding basement of the Pannonian Basin (Northern Croatia, Yugoslavia). Geologija Lju- bljana. 34: 81–253. https://doi.org/10.5474/geologi- ja.1991.004

Pamić J, Lanphere M & Mckee E (1988). Radiomet- ric ages of metamorphic and associated igneous rocks of the Slavonian Mountains in the southern part of the Pannonian Basin. Acta Geologica. 18: 13–39.

Pamić J, Balen D & Herak M (2002a). Origin and geo- dynamic evolution of Late Paleogene magmatic asso- ciations along the Periadriatic-Sava-Vardar magmatic belt. Geodinamica Acta. 15(4): 209–231. https://doi.or g/10.1080/09853111.2002.10510755

Pamić J, Balen D & Tibljaš D (2002b). Petrology and geochemistry of orthoamphibolites from the Variscan metamorphic sequences of the South Tisia in Croatia – an overview with geodynamic implications. Interna- tional Journal of Earth Sciences (Geologische Rund- schau). 91: 787–798. https://doi.org/10.1007/s00531- 002-0258-y

Public institution Papuk Nature Park (2011). Manage- ment Plan of the Papuk Nature Park. (in Croatian)

Public institution Papuk Nature Park (2020). Study of the management of the number of visitors with an action plan for the management of visits in the Papuk Nature Park. (in Croatian)

Schmid SM, Bernoulli D, Fügenschuh B, Matenco L, Schefer S, Schuster R, Tischler M & Ustaszewski K (2008). The Alps–Carpathians–Dinarides orogenic sys- tem: correlation and evolution of tectonic units. Swiss Journal of Geosciences. 101: 139–183. https://doi. org/10.1007/s00015-008-1247-3

Schmid SM, Fügenschuh B, Kounov A, Maţenco L, Nievergelt P, Oberhänsli R, Pleuger J, Schefer S, Schus- ter R, Tomljenović B, Ustaszewski K & Van Hinsber- gen DJJ (2020). Tectonic units of the Alpine collision

zone between Eastern Alps and western Turkey. Gond- wana Research. 78: 308–374. https://doi.org/10.1016/j. gr.2019.07.005

Schneider P & Balen D (2022). Rapid uplift of Late Cretaceous acidic magma from northern Croatia de- ciphered by studying inclusions in zircon using Ra- man spectroscopy. EGU General Assembly 2022 Conference Abstracts (pp. EGU22-3895). https://doi. org/10.5194/egusphere-egu22-3895

Schneider P, Balen D, Opitz J & Massonne H-J (2022). Dating and geochemistry of zircon and apatite from rhyolite at the UNESCO geosite Rupnica (Mt. Pap- uk, northern Croatia) and the relationship to the Sava Zone. Geologia Croatica. 75(2): 249–267. https://doi. org/10.4154/gc.2022.19

Slovenec D, Belak M, Mišur I, Šegvić B & Schuster R (2020). The early Paleozoic cumulate gabbroic rocks from the southwest part of the Tisza Mega-Unit (Mt. Papuk, NE Croatia): evidence of a Gondwana suture zone. International Journal of Earth Sciences (Ge- ologische Rundschau). 109: 2209–2233. https://doi. org/10.1007/s00531-020-01896-8

Tajder M (1956). The albite volcanic rocks of the Voćin area and their origin. Acta geologica 1 (Prirodoslovna istraživanja JAZU). 27: 35–48. (in Croatian with En- glish summary)

Tajder M (1960). Anorthoclase-aegirine-rhyolite from Rupnica brook near Voćin. Acta geologica 2 (Priro- doslovna istraživanja JAZU). 29: 95–101. (in Croatian with English summary)

Tomkeieff SI (1940). The basalt lavas of the Giant’s Causeway district of Northern Ireland. Bulletin of Volcanology. 6: 89–143. https://doi.org/10.1007/ BF02994875

Toramaru A & Matsumoto T (2004). Columnar joint morphology and cooling rate: A starch-water mixture experiment. Journal of Geophysical Research. 109: B02205. https://doi.org/10.1029/2003JB002686

Tuffen H, Gilbert J & McGarvie D (2001). Products of an effusive subglacial rhyolite eruption: Blahnukur, Torfajokull, Iceland. Bulletin of Volcanology. 63(2–3): 179–190. https://doi.org/10.1007/s004450100134

Ustaszewski K, Schmid SM, Lugović B, Schuster R, Schaltegger U, Bernoulli D, Hottinger L, Kounov A, Fügenschuh B & Schefer S (2009). Late Cretaceous in- tra-oceanic magmatism in the internal Dinarides (north- ern Bosnia and Herzegovina): implications for the col- lision of the Adriatic and European plates. Lithos. 108: 106–125. https://doi.org/10.1016/j.lithos.2008.09.010

Ustaszewski K, Kounov A, Schmid SM, Schaltegger U, Krenn E, Frank W & Fügenschuh B (2010). Evolu- tion of the Adria–Europe plate boundary in the north- ern Dinarides: from continent–continent collision to back-arc extension. Tectonics. 29: TC6017. https://doi. org/10.1029/2010TC002668

Weinberger R (2001). Evolution of polygonal pat- terns in stratified mud during desiccation: The role of flaw distribution and layer boundaries. GSA Bulletin. 113(1): 20–31. https://doi.org/10.1130/0016-7606(200

Zhao L, Wang L, Tian M & Wu F (2017). Geochem- istry and zircon U-Pb geochronology of the rhyolitic tuff on Port Island, Hong Kong: Implications for early Cretaceous tectonic setting. Geoscience Frontiers. 8(3): 565–581. https://doi.org/10.1016/j.gsf.2016.05.009