Karst System Research Articles

Assessment of diffusive recharge at regional scale in karst systems, Amazon basin, Ecuador

Karst systems play a crucial role in drinking water supply and biodiversity. Therefore, their comprehensive characterization of their recharge mechanisms is essential for safeguard water resource quality and sustainable management. This study utilizes the open-access thematic cartography of Ecuador to analyze the lithological, geomorphological, and hydrographic boundaries of groundwater karst systems. Then, employing the APLIS method, the diffusive recharge for each system is calculated. As a result, five hydrogeological karst systems and twelve small isolated karst bodies were delineated within the Amazon Region of Ecuador, collectively occupying a total area of 11112 km2. The effective infiltration of precipitation within these systems was found to vary between 11.11% and 77.11%, with a predominant low infiltration range (20%-40%) covering 48% of the karst area. One of the multiple parameters used in water balance calculations is defined by this research, serving as the initial step in creating conceptual models for evaluating karst aquifer systems in Ecuador from a hydrological perspective. The findings not only offer insights into the vulnerability and characteristics of these systems but also establish the foundation for informed decision-making in sustainable groundwater resource management and protection efforts.

Defining the functioning of developed karst systems through extreme climatic event analysis, application to the Krupaja Spring, Serbia

Defining the functioning of developed karst systems through extreme climatic event analysis, application to the Krupaja Spring, Serbia

Groundwater Model for Karst and Pelitic Aquifer Systems from a Semi-Arid Region Under Climate Change Scenarios: A Case Study in the Vieira River Watershed, Brazil

Water scarcity is a global issue, especially in semi-arid and arid regions where precipitation is irregularly distributed over time and space. Predicting groundwater flow in heterogeneous karst terrains, which are essential water sources, presents a significant challenge. This article integrates geology, hydrology, and water monitoring to develop a pioneering conceptual and numerical model of groundwater flow in the Montes Claros Region (Vieira River Watershed, Brazil). This model was evaluated under various climate change scenarios, considering changes in rainfall, groundwater consumption, and population growth over the current century. The results indicate that a decline in water table levels is inevitable, primarily driven by population growth and high pumping rates rather than rainfall fluctuations. This underscores the urgent need for improved monitoring, model upgrading, and more importantly, targeted water resource management for Montes Claros.

A cave with remarkably high subterranean diversity in Africa and its significance for biodiversity conservation

Aziza cave, which is also known as kef Aziza or Tazouguert cave, represents an important and large karstic system that consists of more than 3.5 km of surveyed galleries, standing as the fifth most extensive cave system in Morocco and one of the ten largest in North Africa. This study unveils Aziza cave as an important spot of subterranean diversity in Africa. Here, we provide the first checklist of subterranean fauna in this cave, with 26 taxa, comprising 22 troglobiotic and 4 stygobiotic species. Of this total, eight species still require further confirmation of their status. The richest taxa include Coleoptera (5 species), Araneae (4 species), Entomobryomorpha (3 species), and Isopoda (2 species). However, it is noteworthy that only around 34.6% of the cave-restricted species found in the cave have been formally described to date. Additionally, the biodiversity of large system areas remains to be discovered as these areas need to be further explored. Furthermore, this paper highlights the broader conservation challenges faced by subterranean habitats in Morocco, particularly considering human-induced impacts on these remarkable ecosystems. We aim to draw attention to the crucial ecological role of subterranean environments and their extraordinary biological diversity. By doing so, we aim to inspire increased research and conservation initiatives, not just in this area but across Africa.

Modelling hazards impacting the flow regime in the Hranice Karst due to the proposed Skalička Dam

Abstract. This study examines the hydrogeological hazard associated with the construction of the proposed Skalička Dam in the vicinity of the Hranice Karst. Prompted by the catastrophic regional floods in 1997 and 2010, the design of the dam aims to mitigate floods along the Bečva River downstream of the reservoir. However, concerns have been raised regarding the potential disturbance of the natural groundwater regime in the Hranice Karst and the source of mineral waters for the Teplice spa. This is due in particular to the dam's location in an area with limestone outcrops potentially susceptible to surface-water infiltration. Previous studies have also highlighted the strong correlation between the water level in the Bečva River and the water level in karst formations such as the Hranice Abyss, Zbrašov Aragonite Caves, and other caves in the locality. To address these concerns, a nonlinear reservoir-pipe groundwater flow model was employed to simulate the behaviour of the Hranice Karst aquifer and specifically the effects of the dam reservoir's impoundment. The study concluded that the lateral variant of the dam would have a practically negligible impact on the karst water system, with the rise in water level being only a few centimetres. The through-flow variant was found to have a more significant potential impact on water levels and the outflow of mineral water to the spa, with a piezometric rise of about 1 m and an increase in the karst water discharge to the Bečva River of more than 50 %. Based on these results, recommendations for further investigations concerning the design of the dam and its eventual construction have been formulated to reduce geological uncertainties and to minimize the potential impact of the hydraulic scheme on the hydrogeology of the karstic system.

Solute Transport in a Multi-Channel Karst System with Immobile Zones: An Example of Downtown Salado Spring Complex, Salado, Texas

To investigate the influence of flow rate increment on the solute transport parameter of immobile zones in a karst system, a dye tracer test was conducted in the Downtown Salado Spring Complex (DSSC) comprising three springs: Big Boiling, Anderson, and Doc Benedict springs. The Multiflow two-region nonequilibrium model (2RNE) was used to simulate the breakthrough curve (BTC) of the springs, and changes in the solute transport parameters in response to flow rate increment were observed. The simulation result showed that the 2RNE model was capable of reproducing the BTC of all the DSSC springs, with an R-squared value greater than 0.9 in all flow rate increment scenarios. The research demonstrates that a positive correlation will exist between the flow rate and solute transport parameter of the immobile zones if the tracer transport to the spring is truly influenced by immobile zones. In contrast, a negative correlation will exist between the flow rate and mass transfer coefficient if the immobile zone has less influence. Overall, the research provides insights into contaminant movement in karst by documenting how tracers are retained in the immobile fluid zone.

Co-tracing of groundwater emerging and conventional contaminants in karst area of Southwest China: Distribution, source, and source-specific health risks

Groundwater pollution, especially the co-contamination of emerging contaminants (ECs) and conventional contaminants, is increasingly recognized as a critical global concern. Karst aquifers are found throughout the world, providing critical water sources and sustaining rivers and ecosystems, and are particularly vulnerable to pollution. In this study, the characteristics of the contamination, potential sources, and source-specific health risks related to ECs and trace elements (TEs) in groundwater were investigated in Guiyang (a typical karst city in southwest China). The distribution of TEs followed a descending pattern from urban to suburban to rural areas. Interestingly, the median concentrations (217 ng/L) of Caffeine (CAF) in rural areas surpassed those in urban areas (145 ng/L) and even exceeded these in densely populated areas. Additionally, ECs and TEs in groundwater are significantly influenced by population density and land use types within a 1000-meter buffer zone around the sampling sites. Hierarchical cluster analysis and positive matrix factorization suggested that five principal sources of contamination were identified: domestic wastewater (29.7 %), industrial production (25.2 %), agricultural activities (16.3 %), natural and industrial sources (14.9 %), and pharmaceuticals (13.9 %). Monte Carlo simulations revealed that the non-carcinogenic risks to children in urban areas are significant, with 30.0 % of these risks exceeding the safety threshold. Source-specific health risk assessments indicated that industrial production is the primary source of risk. Furthermore, this study introduces a co-tracing method to trace the sources of groundwater contamination in karst systems. The co-tracing method based on the ECs and TEs allows for a deeper understanding of contaminant transfer and its fate in karst groundwater.

The importance of fault control on the geomorphological development of a deep karst system: the case of Sima GESM, a Z-supercave in southern Spain

The mountainous karst massif in the Sierra de las Nieves (southern Spain) and the adjacent large outcrops of subcontinental mantle peridotites provide this area with a distinctive geodiversity that has contributed to the recent declaration of the zone as the new Sierra de las Nieves National Park. In addition, the Alpine tectonic interaction between the peridotites and the carbonate formations originated an important fault, the Turquillas fault, that has conditioned the existence of two tectonic blocks in the karst massif, one uplifted with respect to the other. Also, in the uplifted block, the fault has contributed to the generation of Z-supercaves, which are defined in this paper as caves with more than 1000 m of development in the vertical or Z dimension. Sima GESM with 1100 m of depth is the only Z-supercave known in the giant geologic structure of the Gibraltar Arc, which speaks on the infrequent and exceptionality of the generation of these underground geomorphological systems. The task of this paper is to present a detailed study of the Turquillas fault and its influence in the surface and underground geomorphology and hydrogeology and in particular its influence and implications in the development of the deep karst system defined as the Z-supercave. It is hoped that this study will contribute to identify the necessary conditions for the development of these extraordinary geomorphic systems which are of great international relevance in the major mountain ranges of the world with high relief karst and deep karst systems.

Capacity and current use of the Formoso river watershed, bonito municipality, Mato Grosso do Sul, Brazil: conflicts and land use recommendations

The lack of planning and feasibility of controlling actions over land use and land cover lead to relationship problems between the support capacity that environmental systems present and their current use. The goal of this study is to figure out how much land can be used based on its natural and physical features and to find problems that come up because of the heavy use in the Formoso River watershed. It is one of the most recognized Brazilian watersheds, due to its turistic, karst, and environmental characteristics. The methodology consists in handling and generating data in a GIS environment by performing a synthesis map, as well as adopting geoprocessing criteria and techniques in the field to validate the collected data. The results indicate that the watershed presents large areas without conflicts due to its numerous conservation units. However, major conflicts are identified in other areas, affecting the balance and dynamics of the landscape and its water resources. Such conflicts arise because soybean crops advance in fragile and karst areas. The highest concentration of confit is only 1.61%, particularly near the Formoso River marsh. Therefore, the present work considers a methodology that assesses the capacity of use through a method that has produced a summary document that is applicable and compatible with the karst system.

Utilizing a multi-tracer method to investigate sulphate contamination: Novel insights on hydrogeochemical characteristics of groundwater in intricate karst systems

Karst environments, especially in Mediterranean area, are highly vulnerable to natural and anthropogenic contamination. This study presents a comprehensive hydrogeochemical assesment of surface water and groundwater across a 2300 km2 catchment area spanning Southern Dalmatia (Croatia) and Western Herzegovina (Bosnia and Herzegovina).For the first time in the study area, data were collected over six years integrating ion analysis, sulphur isotope (δ34S) composition, and physical-chemical analysis of water from 30 locations. The research identified four hydrogeochemical facies (carbonate, sulphate, mixed carbonate/sulphate and chloride), influenced by seawater intrusion, carbonate dissolution, evaporite presence, and human activities.Elevated sulphate levels, often exceeding 250 mg/L, were a main focus of the study due to their potential risks to drinking water quality. The study developed a conceptual model to explain the distribution of sulphates, underscoring the importance of evaporite diapirism and δ34S analysis in tracing sulphate origins. These findings contribute to an improved understanding of karst systems and offer essential data for groundwater protection and legislative measures in the Mediterranean region.

Modelling the effect of free convection on permafrost melting rates in frozen rock clefts

Abstract. This research develops a conceptual model of a karst system subject to mountain permafrost. The transient thermal response of a frozen rock cleft after the rise in the atmospheric temperature above the melting temperature of water is investigated using numerical simulations. Free convection in liquid water (i.e. buoyancy-driven flow) is considered. The density increase in water from 0 to 4 °C causes warmer meltwater to flow downwards and colder upwards, resulting in significant enhancement of the heat transferred from the ground surface to the melting front. Free convection increases the melting rate by approximately an order of magnitude compared to a model based on thermal conduction in stagnant water. The model outcomes are compared qualitatively with field data from the Monlesi ice cave (Switzerland) and confirm the agreement between real-world observations and the proposed model when free convection is considered.

High resolution paleo-hydrological record from carbonate deposits in the Roman aqueduct of Traconnade (Aix-en-Provence, SE France)

Carbonate deposits in aqueducts are potential valuable archives of the paleo-hydrological record for the Roman period. Here we present a new high-resolution record from carbonate deposits in the aqueduct of Traconnade built for Aix-en-Provence, the Roman city of Aquae Sextiae. As its name indicates, Aquae Sextiae was built around an abundant source of water but its development, like other ancient cities, required good management of this resource. Through a combination of micromilling and in situ analysis (LA–MC–ICPMS) on thin sections, we measured geochemical variations (stable isotopes and trace elements) and performed sampling for absolute UTh dating. Geochronological results based on classical correction and isochron along with archaeological observations confirm that the carbonate concretion was deposited in Roman times on the walls of the Traconnade aqueduct, built no later than 140 CE. Carbonate δ18O records document the presence of clear geochemical cycles that were mainly controlled by seasonal variations in temperature and δ18O from water circulating in the aqueduct. Minima and maxima in δ18O values were attributed to summer and winter conditions, respectively in agreement with the modern spring’s stable isotope values. The other tracers (Mg/Ca, Sr/Ca and δ13C) also show marked cyclicities with maximum values occurring in summer. The enrichment in trace elements and δ13C were linked to processes occurring during dry periods (1) prior precipitation of calcite in the epikarst, and/or (2) a greater degree of CO2 degassing. The geochemical variations in Ba/Ca, U/Ca and Rb/Ca behave differently and their co-variations may reflect the water-rock interactions in the karst system. We highlight that these tracers can be used as a proxy for intense precipitation events, occurring frequently during autumn. The Traconnade aqueduct carbonate provides sufficient temporal resolution to explore environmental variations on a seasonal scale. These results suggest the occurrence of a double rainy season (spring and autumn) in the Aix-en-Provence region, similar to present-day climatic conditions. Our study confirms the high potential of aqueduct concretions for paleohydrological studies.

Genesis of honeycomb buildups in the Permian Zechstein Group, Southern North Sea

Seismic interpretation has revealed a hitherto unreported honeycomb pattern of carbonate buildups within the Orchard Platform (Southern North Sea). The Z2 Stassfurt Halite Fm. onlaps the southern margin of the Orchard Platform and is also found infilling Z2 intra-platform lagoons to form salt lakes. Post Z2 evaporation, the deeper Z3 water column drowned the Orchard Platform inhibiting the platform recovery attempted by the Z3 Plattendolomit Fm. The palaeobathymetric variability of the drowned Orchard Platform was sufficient to bring parts of the seafloor into the photic zone allowing for the sporadic growth of the Z3 Plattendolomit Fm. However, the palaeobathymetric lows remained beneath the photic zone ensuring an incomplete regeneration of the Orchard Platform with the creation of a high-frequency network of intra-platform lagoons which mimic the polygonal texture of a honeycomb. Whilst previously accepted as collapse structures or karst systems, this study correlates the development of the honeycomb buildups to variations in seafloor palaeobathymetry which in turn mimic the structural lineaments of the Zechstein subcrop. Syn-depositional instability in the Zechstein subcrop caused the topsets of the Z2 salt lakes to become warped. The warped halite provided seed points for Z3 Plattendolomit Fm. growth which allowed for linear ridges of carbonate to traverse the Z2 salt lakes and eventually connect with the honeycomb buildups. Deposition in the Mesozoic lead to loading of the Zechstein. Halite-filled Z3 lagoons accommodated this loading, which caused a pinching effect on the Z3 honeycomb buildups. The sedimentological understanding provided by this study not only de-risks frontier exploration but also provides insight into carbonate growth in restricted platform recovery scenarios.

Fluid flow zones along fracture corridors inferred from collapse dolines in carbonates of the Irecê Basin, Brazil

AbstractDolines are a widespread karst landform, and understanding their genesis can provide critical information about the evolution of subsurface karstification and groundwater flow zones. Swarms of collapse dolines characterise the karst in the southern portion of the Irecê Basin, northeastern Brazil. The basin comprises Neoproterozoic carbonate rocks approximately 100 m thick, bounded by less soluble older quartzites. The distribution of dolines is not uniform, as they cluster along NW‐SE axes and are aligned with cave systems and permeability (fracture) corridors mostly associated with the Almas and Água de Rega ephemeral rivers. Detailed analysis of sectors with a higher density of dolines was performed using LiDAR and unmanned aerial vehicle (UAV) imagery at two different scales to obtain morphometric and orientation patterns. As expected, scale dependency is observed, with an increase in density and circularity index at a higher resolution scale, but a decrease in mean area and perimeter. However, Voronoi polygon indices are similar at both scales (6.86 ± 1); the departure from the more stable number of six due to the clustering of dolines associated with NW‐SE oriented structural trends. Doline density trends closely follow the orientation of caves, but they are not always spatially connected, since doline evolution may either occlude the caves or be related to a non‐identified void. The Irecê Basin karst system displays a continuum of landforms evolving initially from high‐density fracture zones (fracture corridors), which work as favourable dissolution zones, grading into more integrated flow routes, cave systems and eventually clusters of collapse dolines. Collapse dolines are widespread karst landforms worldwide, and their genesis requires the presence of dissolution voids at depth, which may pass undetected from the surface. In some sectors of Irecê Basin, the existence of doline alignments but the absence of known caves indicates the potential of collapse dolines to be used as proxies for subsurface karstification. Such an approach may allow the prediction of the location and orientation of subsurface zones of enhanced dissolution that can potentially represent productive aquifers or favourable zones for geofluids.

Flooding and flood water storage in karst systems of the Mediterranean region

Flooding is a recurring natural phenomenon that can have both life-giving and destructive aspects. In natural environments, floods are often an important element of the seasonal hydrologic cycle that provides water and nutrients to soil, supporting unique, rich and diverse ecosystems. However, flood events can also represent a destructive force that can endanger lives and cause significant damage in urban areas. Karst areas, in particular, are unique because of their special hydraulic characteristics in terms of flood occurrence, the dependence of ecosystems on such events, and attempts to actively store and manage floods. In this article, the hydraulic response of karst aquifers to heavy precipitation events, flood generation, and engineering interventions for flood control are discussed using several examples from karst areas in the Mediterranean region. Flooding mechanisms and regulatory structures in karst poljes are considered using several typical examples from the Dinaric mountain range. In addition, different variants of groundwater abstraction for increasing storage capacity and flood control are presented using examples from France and Montenegro. Managed aquifer recharge in karst areas and adjacent aquifers is demonstrated with examples from Jordan and Algeria. Finally, failed attempts at flood storage in karst reservoirs are presented with examples from Spain and Montenegro. These examples of flood retention in karst areas show the wide range of planning and technical measures and remind us of possible risks and failures in implementation as well as some positive and negative impacts on the environment and especially on ecosystems.

Virtual paleoanthropology in karstic environments. The challenging case of the Neanderthal skeleton from Altamura (southern Italy)

The so-called “Altamura Man” is a Neanderthal skeleton found in 1993 in the Lamalunga karstic system (southern Italy). The skeleton, dated between 172 and 130 ka, still lies deep within a spectacular cave setting, partly embedded within speleothems and extensively covered by coralloid concretions. These peculiar location and conditions would entail major difficulties in excavating and preserving ex-situ this skeleton. Moreover, a decision-making impasse among the institutions involved complicates the problem. Yet, this specimen is one of the best preserved human fossils known. From 2016 to 2020 we conducted extensive in situ investigations of the specimen. Digital techniques and equipment specifically designed for or adapted to the difficult underground context were used in carrying out the most detailed assessment possible of the skeleton. The data thus collected provide an insight into human evolution in Europe, and a first glance at the amount of information this amazing specimen could provide. In this contribution, we want to update the information on the problematic context of Lamalunga, and address the innovative solutions put in place in 2016–2020 to study and preserve the Neanderthal skeleton in this fragile environment.

Influence of structural properties and connectivity of initial fracture network on incipient karst genesis

A karst system is generally difficult to characterize especially regarding its conduits pattern and morphology. Even though speleology allows to map the geometry of cave systems, a large part of karst conduit networks remains unreachable underground. In this study, we use a reactive transport model that simulates dissolution processes in fracture networks to investigate the effect of discrete fracture network (DFN) properties (i.e., anisotropy, connectivity, fracture intersection types) on incipient karst conduit geometry under different boundary conditions (constant hydraulic head versus constant flow rate conditions, directional versus concentrated recharge conditions). We focused on single fracture (with/without cementation) and ladder like fracture networks which are common on limestone. Results show that fracture network properties, play a major role in the final geometry and topology of the incipient karst conduit patterns. Also, the presence of cement on a single fracture favors tortuous and curvilinear wormholes, while the anisotropy of the fracture network favors the formation of anisotropic incipient karst geometries. Finally, the type of hydraulic BC remains a major factor that governs dissolution processes and thus incipient karst conduit geometry.

Water quality management in a tropical karstic system influenced by land use in Chiapas, Mexico

Urbanization and agricultural activities are increasingly threatening karstic systems and the water resources they provide, which are crucial for the livelihood of rural communities. Over the past decade, the Río Grande de Comitán-Lagos de Montebello (RGC-LM) watershed in Chiapas, Mexico, has experienced significant deterioration in water quality. The objective of this study was to analyze the water quality dynamics within the RGC-LM watershed and their potential interaction with land use change. We conducted assessments of the physicochemical and microbiological parameters across the lotic and lentic systems within the watershed over eight monitoring campaigns (from 2013 to 2020). Cluster Analysis (CA) and Principal Component Analysis (PCA) statistical approaches were performed to identify the main drivers of surface water quality variability. This approach aims to enhance water quality assessments within environmental management strategies, with the goal of promoting sustainable land-use practices aimed at protecting and improving water resources within the watershed. Our findings indicate that there is a substantial difference between the upper and lower watershed water quality. The good water quality conditions observed in the higher altitude lakes (mountain lakes) can be attributed to favorable topographical features, specific land use patterns, and effective conservation measures. In contrast, the lotic system in the middle watershed and the lentic water bodies located at lower altitudes (plain lakes) exhibit signs of water quality degradation due to inefficient wastewater treatment and agricultural runoff. This study highlights the importance of implementing adequate environmental management strategies to address these water quality challenges effectively. Given the pressing issue of water quality degradation in karstic environments, we emphasize on the need for increased monitoring frequency of both physicochemical and microbiological parameters to capture seasonal variations and to further understand the vulnerability of water resources in karstic environments. Future assessment of nutrient concentrations, pesticide levels, and blue-green algae populations will also be crucial for evaluating the trophic state at a watershed scale.

Morphogenesis of a new landform responsive to superimposed sulphate and carbonate buried karsts in western Canada

AbstractThis case study documents the morphogenesis of a new karstic landform consisting of a calcareous tufa mound with a gypsum caprock that developed at La Saline Lake, an oxbow of the Athabasca River in northeast Alberta, Canada. The development of this novel tiered carbonate‐sulphate mound architecture resulted from emplacement above the largest known buried salt dissolution collapse‐subsidence karst system, which was superimposed by an overlying buried carbonate karst. Dissolution of halite‐anhydrite beds of the Middle Devonian Prairie Evaporite, 200 m below, resulted in extensive subsidence and brecciation of the overlying Upper Devonian carbonate karst. The development of the calcareous tufa mound resulted from meteoric‐driven groundwater mixed with glacial meltwater directed along the shallowly buried karstic carbonate terrain of Middle‐Upper Devonian limestone, which floors the McMurray Formation, the main reservoir unit of the Lower Cretaceous Athabasca Oil Sands. The flow of the carbonate‐saturated groundwater responsible for the emplacement of the tufa mound at the surface was abruptly redirected along a deeper pathway in response to dissolution collapse‐subsidence adjustments within the buried salt karst of the Prairie Evaporite. The redirected groundwater flow deeper encountered an anhydrite‐dominated segment of the salt dissolution trend within the Prairie Evaporite Formation, 200 m below. The migration up‐section of sulphate‐saturated brine onto the surface of the previously formed tufa mound resulted in the emplacement of a gypsum caprock as a new landform.

The problem of anthropogenic microfibres in karst systems: Assessment of water and submerged sediments

A new worrying micropollutant threathens natural environments: the microfibres (MFs). Natural, regenerated and synthetic MFs have been detected in different environments, as well as in organisms. While synthetic MFs are generally detected in microplastic analyses, natural and regenerated MFs are not taken into account, or are wrongly considered plastics. They are generally considered biodegradable even if their degradation processes in ecosystems are poorly known. Their potential faster degradation could release toxic compounds, and their characteristics could led to a long-term accumulation in the environment. Understanding their dangerousness and the possible impact they could have on ecosystems is fundamental for environment conservation. We collected and investigated water and submerged sediment samples in different caves and springs of the Classical Karst Region (NE Italy), rich in protected habitats and species. MFs were analysed via microscopy and spectroscopy. MFs were found in all samples, highlighting pollution in surface and subterranean habitats of the karst system. MF concentration was higher in submerged sediments respect to waters, highlighting an accumulation of MFs over time. Big microfibres were less abundant, and MF amount increased with the decrease in the considered size. More than 80% of fibres were fluorescent under UV light. Fluorescent MFs were especially transparent, while non-fluorescent ones were mainly black and blue. Most MFs were cellulosic, and synthetic MFs represent only 15–22%, highlighting a significant gap between the MF composition detected in natural environments and the global production of synthetic textiles in recent times. Synthetic MFs were more abundant in waters. Our results improve the knowledge on micropollutants in karst environments, laying the foundations for future research. MF pollution monitoring in karst areas must become a priority for species protection, habitat conservation, and waters management, improving analyses on a larger number of aquatic environments, taking into account the ecological connections between surface and subterranean habitats.