MORPHOLOGICAL FEATURES AND FORMATION CONDITIONS OF THE ALMOPIA SPELEOPARK CAVES (LOUTRA ALMOPIAS, N. GREECE) MORFOLOŠKE ZNAČILNOSTI IN POGOJI NASTANKA JAM V SPELEOPARKU ALMOPIA (LOUTRA ALMOPIAS, SEVERNA GRČIJA) Georgios LAZARIDIS 1,* & Vasilios MELFOS 2 Abstract UDC 551.435.84(495) Georgios Lazaridis & Vasilios Melfos: Morphological features and formation conditions of The Almopia Speleopark caves (Loutra Almopias, N. Greece) The Almopia Speleopark caves are located at the Almopia ba- sin in northern Greece, at the foothill of Voras Mountain, and are formed in the Maestrichtian limestones of the Pelagonian zone. They are studied on the basis of their meso- and micro- scale morphology as well as their horizontal pattern, in order to investigate the character of the forming aquifer. Emphasis is given on the morphological description of the Loutra Almo- pias Cave. Cave morphology is dominated by the presence of cupolas, rock bridges, ridges and “windows”, abrupt termina- tions of fracture guided passages, pendants, rising channels, pseudonotches, false-floors and spongework. Speleogens indi- cate a speleogenesis due to slowly natural convecting hot water bodies. Phreatic calcite from the Varathron Cave is analyzed on the basis of the fluid inclusions in order to investigate the physicochemical conditions of the convecting water bodies. This has shown that the calcite was formed at temperatures ranging between 120 and 189 ºC, with a peak around 150 ºC. The fluids were dominated by NaCl of very low salinities (0.2- 1.0 wt% NaCl equiv.), showing probably the incorporation of meteoric waters. Key words: Phreatic calcite, thermal water, hypogene speleo- genesis, hydrothermal cave, Loutra Almopias Cave, Greece. Izvleček UDK 551.435.84(495) Georgios Lazaridis & Vasilios Melfos: Morfološke značilnosti in pogoji nastanka jam v speleoparku Almopia (Loutra Almo- pias, severna Grčija) Speleopark Almopia je ob vznožju gore Voras v porečju Almo- pije v Severni Grčiji. Jame so nastale v maastrichtskih apnencih Pelagonske cone. Z analizo jamskih skalnih oblik in tlorisov jam smo raziskovali pogoje njihovega nastanka. Največ raziskav smo naredili v jami Loutra Almopias. Med oblikami so značilne kupole, skalni mostovi, skalni grebeni, okna, nezvezni zaključki rovov nastalih ob razpokah, dvigajoči se kanali, lažne korozij- ske zajede in tla ter prepleti drobnih kanalov. Iz skalnih oblik sklepamo, da so jame nastale s korozijo počasnih konvekcijskih tokov termalne vode. Fizikalno-kemične pogoje speleogene- ze smo določali z analizo tekočinskih vključkov v freatičnem kalcitu iz jame Varathron. Temperatura izločanja kalcita je bila med 120°C in 180°C, z vrhom temperaturne porazdelitve pri 150°C. Nizka slanost tekočinskih vključkov, masni delež soli je med 0,1% in 1%, kaže da tekočinske vključke sestavlja pretežno meteorna voda. Ključne besede: freatični kalcit, termalna voda, hipogena spe- leogeneza, hidrotermalna jama, Loutra Almopias, Grčija. ACTA CARSOLOGICA 50/1, 37-48, POSTOJNA 2021 1,2 Faculty of Geology, Aristotle University, Thessaloniki, GR-54124, Greece, geolaz@geo.auth.gr * Corresponding author Received/Prejeto: 02.09.2019 DOI: 10.3986/ac.v50i1.7592 COBISS: 1.01 INTRODUCTION The Almopia Speleopark comprises a group of caves with geological, paleontological and archaeological interest. The so-called “Bear Cave” (Loutra Almopias Cave, former name Loutra Arideas Cave) is important due to a collection of more than 15,000 specimens of large mammals recovered during excavations since 1992 (Tsoukala & Rabeder 2005). Most of the material is attributed to Ursus ingressus (Rabeder et al. 2004). The associated fauna includes spotted cave hyena, cave lion, leopard, wolf, fox, badger, mustelids and artiodac- tyls of Late Pleistocene age (Tsoukala 1994; Tsoukala et al. 1998, 2001; Tsoukala & Rabeder 2005; Pappa et al. 2005; Tsoukala et al. 2006). Additionally, many thou- sands specimens of micromammals have been identi- fied (Chatzopoulou et al. 2001). Since 1990, when the first explorations took place in the Speleopark, more than ten caves and shel- ter caves have been recorded. In 2004, a new effort be- gun, aiming in interpreting the speleogenesis of these caves through detailed morphological and mineralogi- cal investigations on the basis of a resurvey (Lazaridis 2005). The study of these caves is part of a broader on- going research on hypogene caves and speleogenesis in Greece (Lazaridis 2017). In previous works two development phases of the caves were distinguished based on their morphology. The main phase took place in deep-seated conditions by a slowly convecting and upwelling thermal water body, while in the second one surface water, which en- tered the caves, re-sculptured the walls of some cave passages under turbulent flow conditions (Lazaridis 2006). Furthermore, hypogene speleogenesis is also known from the Provalata Cave in North Macedonia (Temovski et al. 2013; Temovski 2017). The aim of the present study is to define the con - ditions under which the main phase of speleogen- esis took place on the basis of the cave morphology at Varathron Cave, and to determine the physicochemi- cal conditions forming the phreatic speleothems (cave clouds) in the context of the fluid inclusion data. Fur- thermore, a detailed description on the dissolution morphology of the Loutra Almopias Cave is presented and discussed. GEOLOGICAL SETTING The Speleopark is located at the Almopia basin in northern Greece, at the foothill of Voras Mountain (Kato Loutraki village). The area under study belongs to the Pelagonian zone of the Hellenides and is situated not far from the geological boundary with the Almopia zone (Fig. 1). The area consists of alpine metamorphic rocks, Mesozoic carbonates and flysch deposits of Up- per Cretaceous-Palaeocene age. Volcanic rocks of in- termediate composition and high potassium character intrude the Almopia zone (Vougioukalakis 2002). The volcanic activity is dated between 5 Ma (million years) at the eastern and 1.8 Ma at the south-eastern part of the Voras Mountain (Kolios et al. 1980). The caves were formed in the Maestrichtian limestones of the Pelagon- ian zone. Geological formations are dominated tecton- ically by the ENE-WSW striking Loutraki Fault and the NNW-SSE striking overthrust of Almopia zone onto the Pelagonian (Fig. 1). Ore-bearing faults (filled with pyrite, Mn- and Fe-oxides) occur north and north-west of the Kato Loutraki village (Mountrakis 1976). The Loutraki Fault bounds the Voras Mt. against the Al- mopia basin, which is due to the Neogene extensional tectonic regime (Chatzidimitriades 1974; Mountrakis 1976; Eleftheriadis 1977). Three uplifted denudation surfaces of Neogene age and two piedmont surfaces of Villafranchian-Vil- lanyian or Pleistocene age had been referred by Psi- lovikos & Kanetsis (1989) at the northern part of the Pelagonian zone. These uplifted movements resulted in the down-cutting erosion of the Thermopotamos River that flows through the Almopia Speleopark area. Numerous Ca-Mg-HCO 3 hot springs with a tem- perature of 37 ºC are found in the Almopia Speleopark area and are probably related to the volcanic activity of the Pliocene (Dimopoulos 1990; Patras 1990). GEORGIOS LAZARIDIS & V ASILIOS MELFOS ACTA CARSOLOGICA 50/1 – 2021 38 MORPHOLOGICAL FEATURES AND FORMATION CONDITIONS OF THE ALMOPIA SPELEOPARK CAVES (LOUTRA ALMOPIAS, N. GREECE) METHODS Investigation of speleogenesis processes that formed the caves is done with detailed survey of the caves interior and underground small- and medium-scale geomorpho- logical features. In particular the identification of hypo- gene speleogenesis (see Palmer 2000, 2007; Klimchouk 2007, 2009, 2016; Dublyansky 2013, 2014; Klimchouk et al. 2017) is based on a number of criteria proposed and discussed in a number of papers (e.g., Klimchouk, 2007; Audra et al. 2009a, b; Dublyansky 2013; Klimchouk et al. 2014; De Waele et al. 2016). Cave morphology is examined on the basis of the horizontal cave pattern (Palmer 2000, 2005) and by recognizing specific meso- and micro-scale morphology in cave passages (Kempe et al. 1975; Bögli 1978; White & Deike 1989; Lauritzen & Lundberg 2000; Lundberg 2005; Klimchouk 2007), in or- der to investigate the character of the host aquifer (Ford, 2000; Klimchouk, 2007). Basic terminology is given in Fig. 2. Morphometrical parameters are calculated ac- cording to Klimchouk (2007) and Lazaridis (2009). The calcite samples collected from the Varathron Cave were studied by optical microscopy at the Depart- ment of Mineralogy, Petrology, Economic Geology, Aris- totle University of Thessaloniki, for the determination of Fig. 1: Geological sketch-map of NW Almopia Basin (based on Mercier et al. 1988); inset maps: Greece with the Almopia Basin depicted (upper left corner); Almopia Basin and Speleopark (upper middle); 3D view of the geology of the Almopia Speleopark broader area (upper right corner, vertical exaggeration 1.5). Fig. 2: Sketch showing basic geomorphological features and termi- nology used in text. Earlier stage of development is dominated by isolated adjacent cavities that become interconnect with dissolu- tion progress. ACTA CARSOLOGICA 50/1 – 2021 39 the constituent minerals, their texture and the grain size. In addition powders of the samples were analysed by X- ray diffraction (XRD) in order to distinguish calcite from aragonite or dolomite. Microthermometric data were obtained using a LINKAM THM-600/TMS 90 heating-freezing stage cou- pled to a Leitz SM-LUX-POL microscope, housed at the Department of Mineralogy, Petrology, Economic Geol- ogy of the Aristotle University of Thessaloniki. Calibra- tion of the stage was achieved using organic standards with known melting points (chloroform −63.5 ºC, naph- thalene 80.35 ºC, Merck 135 135 ºC, saccharine 228 ºC, Merck 247 247 ºC) and ice (H 2 O)). The precision of the heating and freezing measurements were ±1 °C and ±0.2 ºC, respectively. Fluid inclusion shapes and sizes, spatial relationships among inclusions and minerals and phases within inclusions were observed in three doubly-polished thin sections prepared at the University of Hamburg, Germany. Routine heating-freezing runs were performed on a total of 79 fluid inclusions. The FLINCOR program (Brown 1989) was used to process the fluid inclusion data for calculating salinities and densities. MORPHOLOGY OF THE ALMOPIA SPELEOPARK CAVES Almopia Speleopark caves are of spongework horizon- tal pattern and network mazes or a combination of both patterns. Especially the caves of the northern slope of the Thermopotamos valley are of spongework pattern in con- trast to the network caves of the southern slope. Cave pas- sages are fracture guided and of elliptical shape. In various places assemblages of circular, in cross section, passages that developed along fractures or bedding planes, form boneyard morphology (e.g., Antarton Cave). Cave passag- es are commonly abruptly terminated (“blind” passages). Meso- and micro-scale cave morphology is domi- nated by cupolas, feeders, fissure-like feeders and side-wall feeders, pendants, rock bridges inside chambers and pas- sages, rock ridges that separate chambers, pseudonotches (sensu Osborne 2004), “windows” that connect chambers, in some cases flat passage ceilings, and wide wall pockets flattened upwards. Network mazes present the same spe- leogens except those formed in chambers, such as “win- dows” , rock bridges and bedrock ridges. Exceptionally, in some cave passages scallops and false floors are formed, due to the above-mentioned second phase of develop- ment. The alteration of the original cave morphology by dissolution in the vadose zone is very limited. Cave morphology strongly indicates phreatic speleo- genesis in deep-seated setting by slowly convecting water bodies. However, scallops, detritus sediments, eroded de- posits and false-floors represent invasion by surface wa- ter (Thermopotamos Stream). This second phreatic event followed the main cave development after the uplifting movements of the mountain and the cave drainage. Apart from the speleogens, in an inclined passage of the Varathron Cave calcite formations were found that in- dicate phreatic origin and they are strongly related to the first phreatic event that developed the caves. This deposit covers the perimeter of the cave passage forming in some places cave clouds and is analyzed on the basis of the fluid inclusions in the present paper in order to investigate the physicochemical conditions of the convecting water bodies. DETAILED SOLUTIONAL MORPHOLOGY OF THE LOUTRA ALMOPIAS CAVE The Loutra Almopias Cave is located at the northern slope of the Thermopotamos Valley. The two-dimensional morphometrical parameters of the cave are (Fig. 3): cave area 810 m 2 ; passage length: ~250 m; area of cave field: 1470 m 2 (polygon method, Klimchouk 2007); passage density 168 km -1 ; areal coverage: 55%; V1-score (sensu Lazaridis 2009): -20.69. The limestone is dipping toward the NW . The cave entrance is developed due to valley col- lapse at 540 m above sea level, having no relation to the recent topography. The cave floor at the entrance is about 2 m higher than inside the cave due to collapsed boul- ders. The entrance chamber shows spongework dissolu- tion pattern at the ceiling, a rock bridge and a “window” to the northern part and a pillar that separates a parallel gallery to the rock bridge passage (indicative morpho- logical features from all cave chambers are given in Fig. 4a-g). Blind passages are also present as well as to the rest of the chambers. Two NE-SW striking passages connect the en- trance chamber with the first chamber of the cave (LAC I). The northern one is a narrow and high passage that terminates to bedrock ridges at both ends, which prob- ably constituted a fissure-like feeder (Fig. 5a). The cave floor is about 3 m lower than the entrance chamber and GEORGIOS LAZARIDIS & V ASILIOS MELFOS ACTA CARSOLOGICA 50/1 – 2021 40 MORPHOLOGICAL FEATURES AND FORMATION CONDITIONS OF THE ALMOPIA SPELEOPARK CAVES (LOUTRA ALMOPIAS, N. GREECE) Fig. 4: Dissolution morphology of the Loutra Almopias Cave. a. spongework dissolution; b. com- bination of rock bridge, a window and a bedrock ridge to the back- ground; c. pendants; d. cupolas; e. wall pocket with flat ceiling; f. ris- ing channel; g. small solution pits on the cave-ceiling; h. side-wall feeder. Indicative scale bar: 1m: a, b ,c, e; 0.5m: d; 0.25m: f (Photo: G. Lazaridis). Fig. 3: Geomorphologic sketch- map of the “Bear Cave”. Cave out- line modified after Kampouroglou & Chatzitheodorou 1999. Cupolas are distributed in the whole cave area and only some main features are depicted on map for better ap- pearance. ACTA CARSOLOGICA 50/1 – 2021 41 about 2 m lower than the first chamber. The southern passage is a “window” connection of the two chambers. The chamber LAC I (Fig. 5b) is the larger one divided in three areas (LAC Ia, Ib, Ic). Two of these areas are wall- pockets of about 4-5 m length (or blind passages) con- nected to rising channels. The third area is the southern chamber of the cave. LAC Ic presents blind passages to the eastern side, a small passage connecting the cham- ber to an outlet of the most eastern chamber of the cave (LAC III), cupolas on the ceiling and a ‘window’ to the southern part where the ceiling exerts downwards. LAC I presents boneyard morphology to the northwestern side of the chamber, cupolas on the ceiling, rising chan- nels (LAC Ia and Ib). A small cupola-shaped chamber is located in the central part of LAC I connected with “windows” and rock ridges (Fig. 5). LAC II chamber is connected to LAC I with window and a rock ridge, and to LAC Ia with a “window” near the cave ceiling (about 10 m higher the cave-floor). LAC II presents cupolas and large ceiling pendants as well as a massive ‘old’ flowstone deposited to the eastern blind passage towards the cham- ber (indicative photos of speleothems are given in Fig. 6). LAC III is of similar morphology with cupolas on the ceiling (Fig. 5). To the western side-wall two feeders have been observed, one of which is the connection to LAC Ic. The chamber terminates to the eastern at a blind pas- sage where eroded sensu lato flowstone observed, similar to that of LAC II. The last two chambers are connected through a “window” passage. Cave-walls present in all chambers pockets nearly flat that grade downward into sloping sidewalls. Although, the cave walls and ceilings bear many typical dissolution forms, the features on the cave floor are obscured by sediment fill, which is a common feature of the Almopia Speleopark caves. At least in two caves, the Gremos Cave and the Varathron Cave sediment fill- ing of at least 15-20 m thickness has been removed due to erosion, indicating that a significant vertical cave devel- opment can be considered despite the fact that the caves are significantly sediment filled. However the sediment thickness in Loutra Almopias Cave is unknown. A prob- able feeder close to the eastern sidewall of the chamber LAC I (depicted on Fig. 3), is more than 15 m in depth but not accessible due to the small dimensions of the opening. This feeder could be connected to a lower cave or passage due to seasonal air wind occurrence, indicat- ing a considerable depth regarding the lateral extent of the cave and the sediment thickness. The sediments of the cave are detrital in origin, derived from the wider area of the valley (Tsirambides 1998; Kabouroglou & Chatzitheodorou 1999), with flow- stone intercalations. GEORGIOS LAZARIDIS & V ASILIOS MELFOS Fig. 5: Loutra Almopias Cave: a. The fissure-like feeder passage that connects the entrance chamber with LAC I chamber of the cave. b. View of the main part of LAC I chamber during excavation. The rising channel of LAC Ic is seen on the left side and the passage to LAC II at right side. c. The LAC III chamber of the cave. d. Rock ridge, partially covered by sediments, in the central chamber LAC I (Photo: G. Lazaridis). ACTA CARSOLOGICA 50/1 – 2021 42 MORPHOLOGY, MINERALOGY AND FLUID INCLUSIONS OF PHREATIC CALCITE FROM THE V ARATHRON CAVE The V arathron Cave is located at the northern slope of the Thermopotamos valley (Fig. 1). The two lower entrances of the cave are located at 500 m of altitude and the third one at 520 m. The solutional morphology of the cave is typical for the Almopia Speleopark caves. The higher en- trance is connected to the cave with an inclined passage, whose walls are covered by phreatic calcite (subaque- ous coatings) forming cave clouds (mammilaries, Text. Fig. 7). Cave clouds are speleothems that are common in hypogenic caves. This speleothem has been reported in caves formed by rising thermal water such as Devil’s Hole, Lechuguilla Cave, Carlsbad Caves, Giusti Cave etc. (Hill & Forti 1997) and from the hypogene Petralona Cave in Greece (Lazaridis 2009). The carbonate samples from the speleothems of the Varathron Cave consist mainly of calcite with traces of aragonite (Fig. 8). The calcite is coarse-grained with the maximum grain size reaching up to 1 cm. The carbonate texture is homeoablastic to weakly heteroblastic, demonstrat- ing grains of almost same size, except some which are smaller. Sometimes a slight grain elongation is observed (Fig. 9a). The calcite is relatively free of strain and forms MORPHOLOGICAL FEATURES AND FORMATION CONDITIONS OF THE ALMOPIA SPELEOPARK CAVES (LOUTRA ALMOPIAS, N. GREECE) Fig. 6: Speleothems of the Loutra Almopias Cave. a. and a’. “old” flowstone; b. cave corals; c. cave corals; d. stalactites and flowstone; e. rimstone dams; f. rims on verti- cal side-wall (Photo: G. Lazaridis). ACTA CARSOLOGICA 50/1 – 2021 43 grains that have sutured boundaries, sometimes em- bayed or dentate (Fig. 9b). In contrast, a few anhedral, equigranular crystals have straight boundaries, which in- tersect at 120˚ triple junctions (Fig. 9c). Such equilibrium grain configurations are normally present in thermal an- nealed crystals. In the present study, fluid inclusions data were ob- tained from calcite, which comprises the phreatic speleo- thems in the Varathron Cave. Fluid inclusions were eval- uated using fluid inclusion types and fluid inclusion as- semblages (FIAs). This evaluation ensures that the results were not influenced by different sizes and shapes of fluid inclusions and it helps to eliminate variable data caused by changes after entrapment (Goldstein & Reynolds 1994). Fluid inclusions have regular or irregular shapes, and are isolated or are arranged in clusters and planes. Inclusions with negative crystal or rounded to elongated isometric shapes, which are mostly isolated, often occur along crystal faces of calcite (Fig. 9d). They are assumed to be primary in origin with only a few considered as sec- ondary, according to the criteria of Roedder (1984) and Bodnar (2003). Microthermometric measurements were conducted mainly on primary fluid inclusions; inclusions that have been necked down or were secondary in origin were avoided. Microthermometric results from the fluid inclusions are depicted in Fig. 10. At room temperature, only two phase liquid-vapor inclusions were identified. Inclusions that were analyzed, ranged in diameter between 8 and 95 μm and homogenized into the liquid state. The primary GEORGIOS LAZARIDIS & V ASILIOS MELFOS Fig. 7: Entrance of Varathron Cave and the inclined passaged covered by phreatic calcite coatings. Inlet: Close up of sampling site (Photo: G. Lazaridis). Fig. 8: X-ray diffraction (XRD) diagram of the speleothems at the Varathron cave, consisting mainly of calcite with traces of aragonite. ACTA CARSOLOGICA 50/1 – 2021 44 MORPHOLOGICAL FEATURES AND FORMATION CONDITIONS OF THE ALMOPIA SPELEOPARK CAVES (LOUTRA ALMOPIAS, N. GREECE) Fig. 9: Microphotographs of the calcites from the Varathron Cave, Almopia Speleopark. a. Elongated calcite crystals; b. calcite with sutured boundaries, sometimes embayed or dentate; c. anhedral, equigranular calcite crystals with straight boundaries, which inter- sect at 120º triple junctions; d. a two-phase aqueous fluid inclusion in calcite from the Varathron Cave. L Liquid, V vapor (Photo: V. Mel- fos). Fig. 10: Distribution of the homog- enization temperatures of fluid in- clusions in calcite from the speleo- thems of the Varathron Cave. fluid inclusions have relatively consistent liquid to vapor ratios (15–20 vol % vapor). The variability in homogeni- zation temperature data may be due to real variability in the FIAs or to post-entrapment processes such as thermal re-equilibration or undetectable necking-down (Gold- stein 2003). Homogenization temperatures of all but seven of sixty-nine fluid inclusions range from 120 ºC to 189 ºC, with a peak around 150 ºC (Fig. 13). The other eight in- clusions, possibly secondary, homogenized between 68 ºC and 104 ºC. Eutectic temperatures (Te) of fluids range from −20.5 ºC to −21.5 ºC, suggesting that salts in the flu- ids are dominated by NaCl (Crawford 1981; Shepherd et al. 1985). The final ice melting temperatures in the same inclusions range from -0.1 ºC to -0.6 ºC, corresponding to salinities between 0.2 and 1.0 wt% NaCl equiv. Many inclusions in calcite can be grouped into FIAs which have a constant liquid to vapor ratios (~20%) and variable shapes, clustering along linear trends that do not cross cleavage or growth boundaries. FIAs showed restricted temperature ranges such as 131 ºC to 148 ºC, 148 ºC to 163 ºC, 153 ºC to 161 ºC and 154 ºC to 166 ºC, showing that the assemblages are true FIAs and, therefore, that these ranges probably represent primary Th. In contrast, the low Th fluid inclusions (Fig. 10) were affected by sec - ondary processes, such as necking-down and leakage. ACTA CARSOLOGICA 50/1 – 2021 45 GEORGIOS LAZARIDIS & V ASILIOS MELFOS CONCLUSIONS The Almopia Speleopark comprises a group of caves formed in the Maastrictian limestones of the Pelagon- ian zone, at the southern foothills of the Voras Moun- tain. The intrusion of the volcanic rocks and the thermal gradient are related to a number of Ca-Mg-HCO 3 hot springs, found in the area, with a present temperature up to 37 ºC at the Almopia Speleopark area (Dimopoulos 1990; Patras 1990). The Loutra Almopias Cave and the Varathron Cave are the most important karst features showing various morphologies, which indicate phreatic conditions. Ac- cording to Lazaridis (2006) slowly convecting water bodies, possibly related to the hydrothermal field of the area, reacted with the limestones resulting in a compli- cated cave system. The Loutra Almopias Cave revealed the most of the speleogens that have been recorded in the Almopia Speleopark caves, which are characteristic for hypogene speleogenesis. The fluid inclusion study con- firms the hydrothermal-hypogene origin of the Almopia Speleopark caves, proposed in previous work, based on the cave morphology and the geological setting of the area (Lazaridis 2006). The phreatic calcite crust found in the Varath- ron Cave, forming cave clouds, is composed of coarse- grained calcite related to hydrothermal fluids. The fluid inclusion study has shown that the calcite has been formed at temperatures ranging between 120 and 189 ºC, with a peak at 150 ºC. The fluids were dominated by NaCl with very low salinities (0.2-1.0 wt% NaCl equiv.), showing probably the incorporation of meteoric waters. This is supported by the Cl - content as well as the rela- tion of the stable isotopes (O 18 and D) which according to Dotsika et al. (2006) indicate a strong contribution of meteoric water. Some parts of the caves were subjected to minor modifications by a second phase of speleogenesis that is related to the geomorphological evolution of the area, which comprises mountain uplift, down-cutting erosion and the invasion of the Thermopotamos Stream into the caves. The latter, re-sculptured some cave walls under turbulent flow conditions and deposited detrital sedi- ments as well as the fossils found in the Loutra Almopias Cave. ACKNOWLEDGEMENTS W e thank Professor Stephan Kempe, University of Darm- stadt, Germany, for his thorough and helpful discussions on the caves’ morphology. We also thank the reviewers of an early version L. Plan, P . Voudouris, A. Klimchouk and P . Audra for their fruitful comments and suggestions. REFERENCES Audra, P., Mocochain, L., Bigot, J.Y. & J.C. Nobécourt, 2009a: Morphological indicators of speleogenesis: hypogenic speleogenesis.- In: Klimchouk, A. & D. Ford (eds.) Hypogene Speleogenesis and Karst Hy- drogeology of Artesian Basins. Institute of Speleology and Karstology, pp. 17-22, Simferopol. Audra, P., Mocochain, L., Bigot, J.Y. & J.C. Nobecourt, 2009b: Hypogene cave patterns.- In: Klimchouk, A. & D. Ford (eds.) Hypogene Speleogenesis and Karst Hydrogeology of Artesian Basins. Institute of Spele- ology and Karstology, pp. 23-32, Simferopol. Bodnar, R.J., 2003: Introduction to fluid inclusions.- In: Samson, I. et al. (eds.) Fluid inclusions: analysis and interpretation. Short course series, Vol. 32, Mineral- ogical Association of Canada, pp. 1-8, Québec. Bögli, A., 1978: Karsthydrographie und Physische Speläol- ogie.- Springer-Verlag, pp. 292, Heidelberg. Brown, P .E., 1989: FLINCOR: A microcomputer program for the reduction and investigation of fluid inclusion data.- American Mineralogist, 74, 11, 1390-1393. Chatzidimitriadis, E.A., 1974: Die Problematik der Er- wärmung von Th ermalquellen und ihre Zusammen- hänge zu den geologisch-tektonischen Verhältnis- sen ihrer Umgebungen anhand von zwei Beispielen, nämlich der Heilquelle von Ano Loutraki bei Aridea West Mazedonien und von Trajanoupolis West Th- rakiens.- Scientific Annals, Faculty of Physics and ACTA CARSOLOGICA 50/1 – 2021 46 MORPHOLOGICAL FEATURES AND FORMATION CONDITIONS OF THE ALMOPIA SPELEOPARK CAVES (LOUTRA ALMOPIAS, N. GREECE) Mathematics, Aristotle University of Thessaloniki, 14, 113-128. Chatzopoulou, K., Vasileiadou, A., Koliadimou, K., Tsoukala, E., Rabeder, G. & D. Nagel, 2001: Prelimi- nary report on the Late Pleistocene small mammal fauna from Loutraki “Bear Cave” (Pella, Macedonia, Greece).- Cadernos do Laboratorio Xeoloxico de Laxe, 26, 485-495. Crawford, M.L., 1981: Phase equilibria in aqueous flu- id inclusions.- In: Hollister, L.S. & M.L. Crawford (eds.) Short course in fluid inclusions: applications to petrology. Mineralogical Association of Canada, pp. 75–100, Québec. De Waele, J., Audra, P., Madonia, G., Vattano, M., Plan, L., D'Angeli, I.M., Bigot, J.Y . & J.C. Nobécourt, 2016: Sulfuric acid speleogenesis (SAS) close to the wa- ter table: Examples from southern France, Austria, and Sicily.- Geomorphology, 253, 452–467. http:// dx.doi.org/10.1016/j.geomorph.2015.10.019 Dimopoulos, G., 1990: Hydrogeological and hydro- chemical characteristics of thermal mineral waters of Greece.- In: Proceedings of the 2nd Symposium of mineral waters, 41-57. Thessaloniki, (in Greek). Dotsika, E., Poutoukis, D., Michelot, J.L. & W. Klopp- mann, 2006: Stable isotope and chloride, boron study for tracing sources of boron contamination in groundwater: boron contents in fresh and thermal water in different areas in Greece.- Water, Air, and Soil Pollution, 174, 19–32. https://doi.org/10.1007/ s11270-005-9015-8 Dublyansky, Y .V ., 2013: Karstification by geothermal wa- ters.- In: Shroder, J. & A. Frumkin, (eds) Treatise on Geomorphology. Volume 6, Academic Press, pp. 57–71, San Diego. Dublyansky, Y .V . 2014: Hypogene speleogenesis – discus- sion of definitions.- In: Klimchouk, A. et al. (eds.) Hypogene Cave Morphologies. Special Publication 18, Karst Waters Institute, pp.1-3, Leesburg. Eleftheriadis, G., 1977: Study of the volcanic rocks in South Almopia.- PhD Thesis, Aristotle University of Thessaloniki, pp.173. Ford, C.D., 2000: Speleogenesis under unconfined set- tings.- In: Klimchouk, B.A. et al. (eds.) Speleogen- esis. Evolution of karst aquifers. National Speleologi- cal Society, pp. 319-324, Huntsville. Goldstein, R., 2003: Petrographic analysis of fluid inclu- sions.- In: Samson, I.M. et al. (eds.) Fluid inclusions: analysis and interpretation. Short course series, Vol. 32, Mineralogical Association of Canada, pp. 9–53, Québec. Goldstein, R.H. & T.J. Reynolds, 1994: Systematics of fluid inclusions in diagenetic minerals.- SEPM Short Course 31, SEPM Society for Sedimentary Geology, pp.199, Tulsa. Hill, C.A. & P. Forti, 1997: Cave minerals of the world. 2nd ed., National Speleological Society, pp. 463, Huntsville. Kabouroglou, E. & T. Chatzitheodorou, 1999: Geomor- phological changes and sedimentation of the A cave (Agiasma) of Loutraki (Pella, Macedonia, Greece).- In: Proceedings of the 5th Congress of the Greek Geo- graphical Society, 83-93, Athens (in Greek). Kempe, S., Brandt, A., Seeger, M. & F. Vladi, 1975: "Fac- etten" and "Laugdecken", the typical morphological elements of caves developing in standing water.- Annales de Spéléologie, 30, 4, 705-708. Klimchouk, A.B., 2007: Hypogene Speleogenesis: Hydro- geological and Morphogenetic Perspective.- Special Paper no.1, National Cave and Karst Research Insti- tute, pp. 106, Carlsbad. Klimchouk, A.B., 2009: Morphogenesis of hypogenic caves.- Geomorphology, 106, 1-2, 100-117. https:// doi.org/10.1016/j.geomorph.2008.09.013 Klimchouk, A.B., 2018: Advances in understanding of hypogene karst.- In: Stafford, K. & G. Veni (ed.) Hy- pogene Karst in Texas. Texas Speleological Survey, pp. 1-122, Austin, TX, USA. Klimchouk, A.B., Sasowsky, I.D., Mylroie, J., Engel, S.A. & A. Summers-Engel (eds.), 2014: Hypogene Cave Morphologies. Special publication 18, Selected papers and abstracts of the symposium, 2 nd -7 th February 2014, San Salvador Island, Bahamas. Karst Waters Institute, pp.103, Leesburg. Klimchouk, A., Palmer, A.N., De Waele, J., Auler, A.S. & P. Audra (eds.), 2017: Hypogene Karst regions and caves of the world.- Springer, pp. 911, Cham. Kolios, N., Innocenti, F., Manetti, P., Peccerillo, A. & O. Giulani, 1980: The Pliocene volcanism of the Voras Mts (Central Macedonia, Greece).- Bulletin V olcan- ologique, 43, 3, 553–568. https://doi.org/10.1007/ BF02597692 Lauritzen, S.E. & J. Lundberg, 2000: Solutional and ero- sional morphology.- In: Klimchouk, B.A. et al. (eds.) Speleogenesis. Evolution of karst aquifers. National Speleological Society, pp. 408-426, Huntsville. Lazaridis, G., 2005: Speleological research in the Loutra Arideas area (Macedonia, Greece).- Neue Forschun- gen zum Höhlenbären in Europa, 45, 57-64. Lazaridis, G., 2006: Almopia Speleopark (Pella, Mac- edonia, Greece): Morphology-Speleogenesis of the caves.- Scientific annals, School of geology, Aristotle University of Thessaloniki (AUTH), Special V olume 98, 33–40. Lazaridis, G., 2009: Petralona Cave: Morphological anal- ysis and a new perspective on its speleogenesis.- In: ACTA CARSOLOGICA 50/1 – 2021 47 GEORGIOS LAZARIDIS & V ASILIOS MELFOS Klimchouk, A. & D.C. Ford (eds.) Hypogene Speleo- genesis and Karst Hydrogeology of Artesian Basins. Special Paper 1, Ukrainian Institute of Speleology and Karstology, pp. 233–239, Simferopol. Lazaridis, G., 2017: Hypogene Speleogenesis in Greece.- In: Klimchouk, A. et al. (eds.) Hypogene Karst Re- gions and Caves of the World. Springer, pp. 225-239, Cham. Lunberg, J., 2005: Karren.- In: Culver, D.C. & W .B. White (eds.) Encyclopedia of caves. Elsevier Inc, pp. 315– 321, Amsterdam. Mercier, J.J., V ergely, P . & A. Galeos, 1988: Geological map of Greece, Arnissa sheet, scale 1:50.000.- Athens. Mountrakis, D., 1976: Contribution to the knowledge of the geology of the north boundary of the Axios and Pelagonian zones in the K. Loutraki-Orma area (Al- mopia).- PhD Thesis, Aristotle University of Thes - saloniki, pp. 164. Osborne, R.A.L., 2004: The troubles with cupolas.- Acta carsologica, 33, 2, 9-36. https://doi.org/10.3986/ ac.v33i2.289 Palmer, A.N. 2000: Hydrogeologic control of cave pat- terns.- In: Klimchouk, B. A. et al. (eds.) Speleogen- esis. Evolution of karst aquifers. National Speleologi- cal Society, pp. 77-90, Huntsville. Palmer, A.N., 2005: Passage growth and development.- In: Culver, D.C. & W.B. White (eds.) Encyclopedia of caves. Elsevier Inc, pp. 440-444, Amsterdam: El- sevier. Palmer, A.N., 2007: Cave Geology.- Cave Books, pp. 454, Dayton. Pappa, S., Tsoukala, E., Lazaridis, G. & G. Rabeder, 2005: Milk teeth of Quaternary carnivores from Northern Greek Caves.- Neue Forschungen zum Höhlenbären in Europa, 45, 169-182. Patras, D., 1990: Hydrogeological and hydrochemical research of the thermometallic springs in the area of Loutraki, Aridea (Loutra Pozar). In: Proceedings of the 2nd Symposium of mineral waters, 167–185, Thessaloniki, (in Greek). Psilovikos, Α. & E. Κanetsis, 1989: Distribution and pale- ogeographic significance of the planation surfaces on the mountain blocks of the Pelagonian massif.- Bulletin of the Geological Society of Greece, 23, 1, 279-287. Rabeder, G., Hofreiter, M., Nagel, D. & G. Withalm, 2004: New taxa of Alpine Cave Bears (Ursidae, Camivo- ra).- Cahiers Scientifiques, 2, 49-67. Roedder, E., 1984: Fluid Inclusions.- Reviews in Mineral- ogy, 12, 1-646. Shepherd, T., Rankin, A. & D. Alderton, 1985: A practical guide to fluid inclusion studies.- Blackie and Son, pp. 239, Glasgow. Temovski, M., Audra, P., Mihevc, A., Spangenberg, J.E., Polyak, V., McIntosh, W. & J.Y. Bigot, 2013: Hypo- genic origin of Provalata Cave, Republic of Macedo- nia: a distinct case of successive thermal carbonic and sulfuric acid speleogenesis.- International Jour- nal of Speleology, 42, 3, 235-246. https://dx.doi. org/10.5038/1827-806X.42.3.7 Temovski, M., 2017: Hypogene Karst in Macedonia.- In: Klimchouk, A. et al. (eds.) Hypogene Karst Re- gions and Caves of the World. Springer, pp. 241-256, Cham. Tsirambides, A., 1998: Study of floor sediments from the Agiasma Cave of Loutraki, Pella (Macedo- nia, Greece).- Bulletin of the Geological Society of Greece, 34, 339-349. (in Greek). Tsoukala, E., 1994: Bärenreste aus Loutraki (Mazed., Griechen.).- In: 2 nd Höhlenbären Symposium, 15 th - 18 th September 1994, Corvara. Tsoukala, E. & G. Rabeder, 2005: Cave bears and Late Pleistocene associated faunal remains from Loutra Arideas (Pella, Macedonia, Greece). 15 years of re- search. - Neue Forschungen zum Höhlenbären in Europa, 45, 225-236. Tsoukala E., Rabeder, G. & S. Verginis, 1998: Ursus spe- laeus and associated faunal remains from Loutraki (Pella, Macedonia, Greece). Excavations of 1996. Geological and geomorphological approach.- In: Proceedings of the 4th International Höhlenbären Symposium, Velenje. Tsoukala E., Rabeder, G. & S. Verginis, 2001: Ursus spe- laeus and Late Pleistocene associated faunal remains from Loutraki (Pella, Macedonia, Greece). Excava- tions of 1999.- Cadernos do Laboratorio Xeoloxico de Laxe, 26, 441-446. Tsoukala, E., Chatzopoulou, K., Rabeder, G., Pappa, S., Nagel, D. & G. Withalm, 2006: Paleontological and stratigraphical research in Loutra Arideas Bear Cave (Almopia Speleopark, Pella, Macedonia, Greece).- Scientific Annals, School of Geology, Aristotle Uni- versity (AUTH), special volume, 98, 41-67. Vougioukalakis, G., 2002: Petrological, geochemical, vol- canological study of the Almopias Pliocene volcanic formations and their correlation with the geothermal manifestations in the area.- Ph.D Thesis, Aristotle University of Thessaloniki, pp. 303. White, B.W . & H. G. Deike III, 1989: Hydraulic geometry of cave passages.- In: White, W. & E. White (eds.) Karst Hydrology. Concepts from the Mammoth Cave Area. Van Nostrand Reinhold, pp. 223-258, New Yor k . ACTA CARSOLOGICA 50/1 – 2021 48