# Carbonate rocks and karst water resources in the Mediterranean region

Julian Xanke, Nico Goldscheider, Michel Bakalowicz, Juan Antonio Barberá, Stefan Broda, Zhao Chen, Mohamed Ghanmi, Andreas Günther, Andreas Hartmann, Hervé Jourde, et al.

► To cite this version:

Julian Xanke, Nico Goldscheider, Michel Bakalowicz, Juan Antonio Barberá, Stefan Broda, et al. Carbonate rocks and karst water resources in the Mediterranean region. Hydrogeology Journal, 2024, 32 (5), pp.1397-1418. 10.1007/s10040-024-02810-1. hal-04873735

HAL Id: hal-04873735

https://hal.science/hal-04873735v1

Submitted on 8 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Carbonate rocks and karst water resources in the Mediterranean region

Julian Xanke<fnref n="1" />,<fnref n="2" /> ID · Nico Goldscheider<fnref n="2" /> · Michel Bakalowicz<fnref n="3" /> · Juan Antonio Barberá<fnref n="4" /> · Stefan Broda<fnref n="5" /> · Zhao Chen<fnref n="6" /> · Mohamed Ghanmi<fnref n="7" /> · Andreas Günther<fnref n="5" /> · Andreas Hartmann<fnref n="6" /> · Hervé Jourde<fnref n="3" /> · Tanja Liesch<fnref n="2" /> · Matías Mudarra<fnref n="4" /> · Marco Petitta<fnref n="8" /> · Nataša Ravbar<fnref n="9" /> · Zoran Stevanović<fnref n="10" />

Received: 10 October 2023 / Accepted: 3 July 2024 / Published online: 3 August 2024
© The Author(s) 2024

## Abstract

Carbonate rocks in the Mediterranean region form karst landscapes with a variety of morphological and hydrological features, and are of particular interest from a water management perspective as they represent major karst aquifers. The Mediterranean Karst Aquifer Map and Database (MEDKAM) provides a 1:5,000,000 scale map showing the distribution of carbonate and evaporite rocks that can host karst groundwater resources, with additional information on other hydrogeological settings, selected terrestrial and submarine karst springs, caves and karst groundwater-dependent ecosystems. A statistical evaluation shows that carbonate rocks cover ~39.5% of the Mediterranean region within a 250-km focus area from the coastline. North Africa has the largest continuous area of carbonate rocks, while smaller countries in the Middle East and the Dinarides have the largest proportion of carbonate rocks in relation to their total area. Carbonate rocks are also widespread in coastal areas, occurring along ~33.6% (14,000 km) of the total Mediterranean coastline, including large islands such as Crete and Mallorca, and ~25.9% (6,400 km) of the continental coastline. Two additional maps display (1) groundwater recharge, showing a climatic gradient from north to south, and (2) groundwater storage trends, indicating a mean annual karst groundwater loss from 2003 to 2020 of 436 million m<fnref n="3" /> in the 250-km area. This study quantifies the carbonate rocks in the Mediterranean region and shows their importance for groundwater resources. MEDKAM will serve as a basis for further research and improved international cooperation in karst groundwater management.

**Keywords** Hydrogeological mapping · Carbonate rocks · Karst groundwater resources · Groundwater recharge · Mediterranean region

<footnote>☒ Julian Xanke
julian.xanke@tzw.de</footnote>

<footnote>☒ Nico Goldscheider
nico.goldscheider@kit.edu</footnote>

<footnote n="1">¹ TZW: DVGW-Technologiezentrum Wasser (German Water Centre), Section Risk Management, Karlsruher Straße 84, 76139 Karlsruhe, Germany</footnote>

<footnote n="2">² Division of Hydrogeology, Karlsruhe Institute of Technology (KIT), Institute of Applied Geosciences (AGW), Kaiserstr. 12, 76131 Karlsruhe, Germany</footnote>

<footnote n="3">³ HydroSciences Montpellier (HSM), Univ. Montpellier, CNRS, IRD, Montpellier, France</footnote>

<footnote n="4">⁴ Department of Geology and Centre of Hydrogeology, University of Málaga (CEHIUMA), 29071 Málaga, Spain</footnote>

<footnote n="5">⁵ Federal Institute for Geosciences and Natural Resources (BGR), Wilhelmstr. 25–30, 13593 Berlin, Germany</footnote>

<footnote n="6">⁶ Faculty of Environmental Sciences, Institute of Groundwater Management, Technische Universität Dresden, Bergstraße 66, 01069 Dresden, Germany</footnote>

<footnote n="7">⁷ Laboratory of Mineral Resources and Environment (LRME), Department of Geology, Faculty of Sciences of Tunis, University of Tunis El Manar, 2092 Tunis, Tunisia</footnote>

<footnote n="8">⁸ Department of Earth Sciences, Sapienza University of Rome, P.Le A. Moro 5, 00185 Rome, Italy</footnote>

<footnote n="9">⁹ Karst Research Institute ZRC SASU, Titov Trg 2, 6230 Postojna, Slovenia</footnote>

<footnote n="10">¹⁰ Department of Hydrogeology, Centre for Karst Hydrogeology, Faculty of Mining & Geology, University of Belgrade, Djusina 7, 11000 Belgrade, Serbia</footnote>
**Introduction**

Carbonate rocks form karst aquifers that serve as impor-
tant freshwater resources for humans and ecosystems, but
pose an ongoing challenge for research and management
due to their unique hydraulic properties (Goldscheider
and Drew 2007; Stevanović 2019). Recent evaluations
by Chen et al. (2017) indicate that exposed karstifiable
carbonate rocks cover 21.6% of the European land sur-
face; however, the total area of carbonate rock occurrence
is much larger, as it also includes nonexposed carbonate
rocks confined by younger formations. The exposed karst
surfaces correspond to the main recharge zones of karst
aquifers, which are often hydraulically connected over
large areas and across national boundaries (Hartmann
et al. 2014). Karst aquifers are often the source of large
springs (Bonacci et al. 2023; Fan et al. 2023; Stevanović
and Milanović 2023) and have unique underground cave
systems, both of which host diverse ecosystems (Siegel
et al. 2023). Karst landscapes serve as recreational areas,
and many karst features are classified as World Heritage
Sites (Goldscheider 2019; Gunn 2021). Furthermore, car-
bonate rocks serve as important industrial raw material,
host oil and gas reservoirs, contain commercial ore bod-
ies, and provide information about the Earth’s evolution
through biological, chemical, mineralogical, and isotopic
signatures (Mackenzie 1978), and are part of the global
carbon cycle (Hartmann et al. 2009).

Carbonate rocks are a group of sedimentary rocks composed of carbonate minerals. Typical examples are limestone, composed of calcite (CaCO₃), and dolomite (or dolostone), composed of dolomite [CaMg(CO₃)₂] (note that the term “dolomite” is used both for the mineral and the rock). The metamorphism of limestone and dolomite leads to the formation of marble, which has a nearly identical mineralogical and geochemical composition as the parent rock but with larger crystals and less intergranular porosity (Goldscheider and Drew 2007). The peculiarity of carbonate rock is related to its strong solubility, which occurs mainly along fractures and bedding planes due to the flow of (ground)water with CO₂ from the atmosphere and soil. This karstification occurs mainly in the shallow phreatic zone and is a self-reinforcing process leading to the formation of a network of channels (Ford and Williams 2007; Bakalowicz 2005) that are embedded in and interact with a rock matrix that is less karstified. Karst aquifers are characterised by a high degree of hydraulic anisotropy (higher conductivity in the direction of the channels than perpendicular to them), heterogeneity (higher porosity and permeability in the channels than in the matrix) and temporal variability in water availability (variable water level, storage volume

and spring discharge) and water quality (Ravbar 2013;
Stevanović 2015). Managing karst aquifers is, therefore,
challenging.

The distribution of karstifiable rocks is well documented
in many countries and regions; however, on a continental
scale, there are only a few approaches that show the dis-
tribution of carbonate rocks. The global distribution of
karst formations was first shown by Ford and Williams
(1989), revised by Williams and Ford (2006) and further
used in Ford and Williams (2007), albeit highly general-
ized. More recently, the World Karst Aquifer Mapping Pro-
ject (WOKAM) has prepared a much more detailed map
and database that is available in digital and printed form at
scales of 1:25,000,000 and 1:40,000,000 (Chen et al. 2017;
Goldscheider et al. 2020); however, the inconsistent data-
base at the supra-regional level is the difficulty of such a
large-scale map development. In recent decades, large-scale
continuous maps of geologic, lithologic, and water resource
information have also been developed, although the degree
of accuracy generally decreases with the size of the area
under consideration. The Global Lithologic Map (GLiM),
for example, provides lithological information translated
from existing regional geologic maps and literature; how-
ever, the merging of different datasets often involves uncer-
tainties and regional differences (Hartmann and Moosdorf
2012). In order to better harmonize such differences, the
International Hydrogeological Map of Europe (IHME1500)
has been developed since the 1960s under the auspices of
the International Association of Hydrogeologists (IAH) and
with the support of the Commission on the Geological Map
of the World (CGMW). In 2002, the Worldwide Hydrogeo-
logical Mapping and Assessment Programme (WHYMAP)
was launched by the United Nations Educational, Scientific
and Cultural Organization (UNESCO) and the German Fed-
eral Institute for Geosciences and Natural Resources (BGR)
to map invisible subsurface water resources (Richts et al.
2011). In this context, various maps have been developed
that provide global information on groundwater resources,
groundwater vulnerability to flood and drought events
(Richts and Vrba 2016), and there is also the already men-
tioned WOKAM (Chen et al. 2017).

The objective of this study is to provide a compendium of
the developed Mediterranean Karst Aquifer Map and data-
base (MEDKAM; Xanke et al. 2022), based on the Interna-
tional Hydrogeological Map of Europe (IHME1500). The
mapping approach follows WOKAM. MEDKAM was devel-
oped at a digital scale of 1:1,500,000 and printed at a scale
of 1:5,000,000, providing more detailed, more accurate and
new information on the following:

* Quantification of the spatial distribution of different karstifiable rocks in the Mediterranean countries and

<figure><img src="image_2.png" /></figure>
other basic information on other aquifer types depend-
ing on their formation and aquifer properties.

* Quantification of the distribution of coastal karst in the Mediterranean region and the bordering states.

* Terrestrial and submarine karst springs and information on their range of discharge as well as karst caves and their spatial extent.

* Distribution of the long-term mean karst groundwater recharge as well as the trends in karst groundwater storage.

MEDKAM is a contribution to the World-wide Hydroge-
ological Mapping Programme (WHYMAP) of the UNESCO
Intergovernmental Hydrological Programme (IHP) and is
available for download as a digital version at the product
center of German Federal Institute for Geosciences and
Natural Resources (BGR; Xanke et al. 2022).

Geographical setting

The evolution of the Mediterranean karst landscape is
closely linked to geotectonic processes since the late Meso-
zoic (Duggen et al. 2003; Cavazza and Wezel 2003), with
today’s landscape strongly impacted by the collision of Eura-
sia and Africa since the late Cretaceous that created several
mountain belts stretching from North Africa and the Iberian
Peninsula across Europe, the Balkan region and the Cauca-
sus to the Middle East. These characteristic landscapes are
further shaped by local climatic conditions—for example,
higher precipitation rates at higher elevations combined with
the tectonically highly deformed carbonates. This leads to
stronger karstification (development of, e.g. sinkholes, kar-
ren, dolines) than in the lower elevations with a predomi-
nantly dry climate and only slightly tectonically stressed ter-
rain (hilly terrain with deep dry valleys). A partial drying of
the Mediterranean Sea during the late Miocene (Messinian

salinity crisis), 6–5 million years ago (Jolivet et al. 2006;
Roveri et al. 2014), and subsequent water level changes in
the following millions of years until today, have led to deep
karstification at various levels, especially in the coastal area.
As a result, some of the karst systems emerge at the cur-
rent sea level, while some are also found inland above, or
below, which then emerge as submarine springs near the
coast (Bakalowicz 2018).

A clear delineation of the Mediterranean region is dif-
ficult because all common classifications—such as cli-
matic (Koeppen 1936; Peel et al. 2007) or Mediterranean
vegetation cover use (Olson et al. 2001) or geographically
by surface catchments or bordering states—include some
shortage in coverage. In each of these classifications, certain
areas, often near the coast, are not included (e.g. large parts
of Northern Africa when using the climatic classification
approach), while areas far from the Mediterranean Sea are
included (e.g. the entire basin of the Nile when using the
hydrologic basin approach); also, in these cases, only surface
characteristics are decisive for the delineation. In the case
of the hydrogeological consideration of the Mediterranean
area and in particular of the karstifiable rocks, a focus area
of ~250 km around the Mediterranean Sea is suggested by
Siegel et al. (2023), which on one hand includes the most
important karst landscapes along the coast, considers nar-
row countries like Italy as a whole and includes the Medi-
terranean area of the aforementioned criteria. However, the
decisive factor for this type of delineation is that it includes
the unique deep karstification caused by the Messinian salin-
ity crisis, which has partly reached far into the inland around
the Mediterranean Sea (Audra et al. 2004; Fleury et al. 2007;
Bakalowicz 2014).

The focus of the map development of Xanke et al. (2022) and the extension of the karst feature database for karst springs and caves was on the entire Mediterranean region, including countries within the extent of the printed MEDKAM (Fig. 1). The creation of the database on karst

<figure><img src="image_3.png" /><figcaption>Fig. 1 Location of the Mediterranean region, which includes all colored regions/countries and the Mediterranean focus zone (250 km), extension of the printed MEDKAM and the regions and countries of Europe, North Africa and the Middle East evaluated in this study</figcaption></figure>
groundwater-dependent ecosystems (Siegel et al. 2023) was
on the 250 km focus area around the coastline (Fig. 1).

The statistical evaluation of karst distribution in this
manuscript is done for (1) the whole Mediterranean region
(which includes all colored regions/countries as shown in
Fig. 1), (2) the main geographical areas of Europe, Mid-
dle East and North Africa, (3) each individual country that
lies to a significant extent within the 250 km focus area, as
well as Portugal, which is often considered to be part of
the Mediterranean region, and (4) the Mediterranean focus
area (250 km zone). Countries that were not included in the
statistical analysis due to their small size (<250 km²) but
are located in the study area are Monaco (MC), San Marino
(SM), and the Vatican City State (VA).

**Materials and methods**

**Map of karstifiable rocks and other aquifer types**

**Basic mapping approach and legend**

The basis for the MEDKAM is the International Hydro-
geological Map of Europe (IHME1500), which consists of
30 hydrogeological map sheets (five of which are unpub-
lished) covering almost the entire European continent
and parts of the Middle East (BGR and UNESCO 2019;
Günther and Duscher 2019; Duscher et al. 2015) and has
recently been extended by four additional map sheets to
the North African region (unpublished). The development
of the IHME map series began in the 1960s under the
auspices of the International Association of Hydroge-
ologists (IAH) and with the support of the Commission
on the Geological Map of the World (CGMW). The sci-
entific editing, cartography, printing and publication of
the map sheets and explanatory notes were directed and
financially supported by the German Federal Institute for
Geosciences and Natural Resources (BGR) and the United
Nations Educational, Scientific and Cultural Organization
(UNESCO). The map has been available since 2013 in
printed form and as a vector dataset in Shapefile format
for geographic information system (GIS) visualization and
data processing (Duscher et al. 2015). The hydrogeologi-
cal content is based on contributions from the respective
mapped countries, coordinated on a transnational basis.
The attributes of the IHME contain different levels of
lithological description that allow the assignment of the
formations to the MEDKAM units. To fill in the missing
parts in North Africa and the Arabian Peninsula, complete
the Mediterranean mapping region and enable statistical
calculation of the entire countries, the digital version of
the British Geological Survey’s quantitative maps of Afri-
can groundwater resources (MacDonald et al. 2012) and

the United States Geological Survey's bedrock geology of
the Arabian Peninsula (Pollastro 1998) were used. They
have been carefully analyzed with regard to the MEDKAM
lithology and have undergone the same review process as
the IHME; however, for further use of MEDKAM, only
the IHME-based shapefile is available for download (see
section "Map availability").

The proposed map provides information on the distri-
bution of carbonate rocks and evaporites in the Mediter-
ranean and in other hydrogeological settings; however, the
term “karstifiable rocks” mainly refers to carbonate rocks
and evaporites, which consist of easily soluble minerals
and are major karst aquifers. In the following, carbonate
rocks are, therefore, referred to as karst aquifers. Carbonate
rocks generally consist of more than 75% carbonate miner-
als (Ford and Williams 2007), although this is usually not
indicated in geological and lithological maps and must be
derived from general nomenclature. In the case of IHME,
this nomenclature and categorization is given and is, there-
fore, regarded as a reliable basis—for example, limestone,
dolomite or chalk are considered to be relatively pure car-
bonate rocks, as defined above. However, a mixture with
or an alternation between sedimentary and carbonate rocks
is often found in the description and requires classification
by regional experts or derivation from literature and more
detailed regional maps. In this assessment, it is assumed that
the exposed carbonate rocks are karstified to some extent,
without specifying the extent. Therefore, a lithological clas-
sification was developed following the World Karst Aquifer
Mapping Project (Chen et al. 2017) and the geological units
were grouped into four hydrogeologically meaningful units:

1. Karst aquifers in sedimentary and metamorphic carbonate rocks

2. Karst aquifers in evaporite rocks

3. Various hydrogeological settings in other sedimentary and volcanic formations (karst aquifers are possibly present at depth)

4. Local, poor and shallow aquifers in other metamorphic rocks and igneous rocks (no karst aquifers present at depth)

The improvement over WOKAM is that it is no longer
necessary to distinguish between areas with continuous
and discontinuous carbonate (or evaporite) rocks because
the more detailed scale of IHME (1:1,500,000) enables a
more precise delineation of the carbonate (and evaporite)
rock areas. The unit "various hydrogeological settings in
other sedimentary and volcanic rocks" includes a wide
variety of geological formations with different hydrau-
lic properties, but also indicates the possible presence of
confined karst aquifers at greater depth. "Local, poor and
shallow aquifers in other metamorphic and igneous rocks"

<figure><img src="image_4.png" /></figure>
are grouped into one unit. As crystalline rocks generally form the geologic basement, karst aquifers are usually not present at depth in these regions. This improvement allows a more detailed statistical analysis of the distribution of karst aquifers in the Mediterranean and an indication of areas where karst aquifers may be present at depth, covered by other formations.

Similar to WOKAM, the map unit “Karst aquifers in sedimentary and metamorphic carbonate rocks” includes sedimentary limestones, dolomites and chalks, which are mostly biochemical sediments, as well as limestones and dolomites that have undergone some degree of metamorphism, e.g. into marble or limestone shale. Since it is very difficult to distinguish between sedimentary diagenesis and metamorphism, these rock types are combined into a single unit as they are usually referenced by their content of carbonate minerals (>75%; Ford and Williams 2007), although this information is not usually detailed on maps or in publications (nor in the IHME), it is a common approach to rock mapping and characterization. The “Karst aquifers in evaporite rocks” mapping unit, by definition, includes chemical sedimentary rocks that consist mainly of highly soluble minerals such as gypsum, anhydrite or halite. These types of rocks typically form in basins where arid climates cause the water to become supersaturated, resulting in the precipitation of these minerals. The unit of “various hydrogeological settings in other sedimentary and volcanic formations” is probably the largest group and includes unconsolidated sediments of various compositions, including alluvial fans, sandy and gravelly sequences as well as silty and clayey low-conductivity layers. Also included here is the group of consolidated sedimentary rocks such as shale or sandstone. All of these sedimentary and volcanic rocks can be deposited over karst formations. “Local, poor, and shallow aquifers in other metamorphic and igneous rocks” are the most clearly defined. These formations generally do not include carbonate rocks and are usually dominated by low hydraulic conductivity and, hence, low aquifer productivity.

The approach to identifying the nonexposed carbonate rocks is the same as in Chen et al. (2017). Data on covered rock formation are typically poor and must, therefore, be inferred from the extent of individual geosystems (USGS 2021; Pawlewicz et al. 2002; Persits et al. 1997; Pollastro et al. 1999) or other data sources. If karst aquifers are exposed within such a delineated geosystem and are adjacent to sedimentary rocks, it is reasonable to assume that the karst is dipping beneath them. However, the focus here was on large continuous karst areas and not on small-scale sequences of subdivided rock formations as occurs, for example, in mountain ranges due to the strong thrusting and folding. Information on

nonexposed carbonate rocks is important because they may contain significant karst groundwater bodies, which are often present in confined conditions.

## Reclassification, evaluation and iterative improvement of the map

For the reclassification of the lithological units from the IHME to the MEDKAM classification, first, an assignment to the rock groups in the attribute table was made. If the assignment was not conclusive, the more detailed lithological descriptions in the IHME were used as a reference. In a second step, official regional karst and lithological maps were used for comparison with the IHME and, if necessary, further adjustments were made by reclassification. However, the reference maps sometimes have a different scale and, therefore, differ in the topology of the individual geological polygons, and in this case, do not always match the polygons of the IHME. In some regions, a compromise had to be made between accurate selection of individual polygons and consistent classification of karst areas over a larger area, often resulting in generalization from continuous small karst outcrops. As noted previously, the greatest difficulty is in assigning lithological descriptions of sequences of carbonate rocks and sedimentary formations to either unit—for example, marlstone and sandstone intercalations have been assigned to sedimentary formations in many areas. However, when additional carbonate intercalations are described and regional maps also indicate carbonate rocks, an assignment to the karst aquifer unit has been made.

In some regions, the various rock descriptions of the IHME are not very precise. In particular, the alternation of limestone, sand, marl or claystone is difficult to classify as pure karst because of the lack of information on the proportions of each type of rock. For this reason, an attempt has been made to make a classification based on existing geological and lithological maps and descriptions, as well as on regional experts. Nevertheless, it cannot be excluded that, especially in the case of the mixed rock sequences, the assignment could have been made to another hydrogeological unit—for example, shale is by definition a sedimentary rock, but often undergoes some degree of metamorphism (e.g. slate) and could also be classified as a metamorphic rock. In this study, however, slates are classified as sedimentary rocks. The combination “slates, quarzites, sandstones, shales” is assigned to sedimentary rocks, while the combination “slates, phyllites, mica schists and graywackes of various metamorphic degrees, partly gneissic” is assigned to metamorphic rocks. Since these descriptions are not always precisely distinguishable, by adapting them to regional
lithological maps, it is possible that areas with the same IHME classification are assigned to a different MEDKAM unit.

## Database on karst springs, caves and groundwater-dependent ecosystems

The database collects and standardizes existing relevant information related to karst aquifers. In particular, the data on karst springs and caves established in the WOKAM project (Chen et al. 2017) are extended and insufficient data from countries and are being improved and completed. Another new development is the inclusion of selected karst groundwater-dependent ecosystems (KGDE) compiled by Siegel et al. (2023). Due to the large amount of data and the larger scale of MEDKAM, it is possible to show different types of karst features such as the distinction between continental and submarine springs and different types of KGDEs. The supplementary data are based on literature research and data provided by regional experts. In this study, only a selection of the karst features collected for MEDKAM and stored in the database is presented. The criteria used to select the karst features shown in this paper are location within the Mediterranean area (250-km zone), regional importance such as spring flow, or length of cave. Consideration has also been given to consistent documentation of large regions; therefore, a few smaller springs or caves are also shown. It should be noted that the number of karst features, springs, caves, KGDEs stored in the MEDKAM database is only a selection of the actual number of existing features. In many countries, sufficient data sets are not publicly available, or if they are, they are not of sufficient accuracy and quality.

Karst springs are documented with their low-flow and high-flow discharge values, which represent the range of reported discharge data, but are not necessarily the minimum and maximum values. Where only one value was available, it was assigned to the appropriate category (low or high); in some cases, only the mean discharge was available. Special attention is given to coastal and submarine springs. Karst caves are documented with their depth and length where this information was available. A new and important aspect of MEDKAM is the karst groundwater-dependent ecosystems, which provide an indication of the state of water quality. Although these ecosystems are only briefly mentioned, they play an important role locally, often providing habitat for rare species and plants. A multidisciplinary characterization of selected KDGEs is documented by Siegel et al. (2023).

## Assessment of karst groundwater recharge

An existing karst-specific simulation model (VarKarst-R; Hartmann et al. 2020) was used to determine the 30-year average (1990–2019) of groundwater recharge in the

MEDKAM region, at a spatial resolution of 0.25° and using global climate data products as input (precipitation, temperature: GLDAS, Rodell et al. 2004; potential evapotranspiration: GLEAM, Martens et al. 2017). For the MEDKAM region, the recharge map from Hartman et al. (2020) was resolved at a higher spatial resolution and interpolated to include previously unmapped areas. This was done for simplification and presentation reasons and is, therefore, subject to some uncertainty in ~18% of the area, especially in the Middle East and North Africa. From the distribution of mean annual karst groundwater recharge, a country-specific average was calculated.

## Assessment of trends in karst groundwater storage

The study of trends in karst groundwater storage was based on the work of Xanke and Liesch (2022). They investigated the extent of changes in groundwater storage (GWS) over the period 2003–2020 for the Euro-Mediterranean region using the latest data from the Gravity Recovery and Climate Experiment (GRACE/GRACE-FO; Save et al. 2016; Save 2020) satellite mission and recently reanalyzed ERA5 land climate data from the European Centre for Medium-Range Weather Forecasts (Hersbach et al. 2020). For this study, these trends in groundwater storage were masked by the MEDKAM carbonate outcrops. It should be emphasized that the GRACE data are based on an accuracy of 0.5°. The values should, therefore, always be considered as a product of their environment and only apply to the entire vertical groundwater column. Overlapping aquifers and catchments are, therefore, included in the signal and, therefore, do not apply exclusively to the karst, but are projected onto it. From the distribution of mean annual trends in karst groundwater storage, a country-specific average was calculated.

## Results

### Distribution of carbonate and evaporite rocks

The analysis of the spatial distribution of the different hydrogeological units (Fig. 2) shows that 26.1% of the whole Mediterranean region is covered by carbonate rocks corresponding to more than 3.2 million km², and only 0.9% consists of exposed evaporite rocks (Table 1; Fig. 3). Much more prevalent are the various hydrogeological settings in other sedimentary and volcanic formations, covering 64.2% of the total Mediterranean region. Local poor and shallow aquifers in other metamorphic and igneous rocks are much less prevalent and account for 8.8%, much of which is found in European countries.

Within the Mediterranean focus area, 39.5% is covered by carbonate rocks and 1.7% by evaporite rocks. This shows

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 2** The Mediterranean Karst Aquifer Map (MEDKAM) with references to other hydrogeological settings, as well as information on the borders between exposed and non-exposed karst, and a selection of karst springs (see table in section “Karst springs”), submarine springs (see table in section “Submarine karst springs”) and caves (see table in section “Karst caves”). An indication of the level of karstification around the Mediterranean basin is given by the Mediterranean Sea level during the Messinian salinity crisis and the almost drying up of the sea (Xanke et al. 2022). The map was created for the entirety of the countries for which this section of the figure provides only partial coverage (e.g. countries of North Africa and the Middle East)</figcaption></figure>

<figure><img src="image_7.png" /></figure>
<caption>**Table 1** Size of the hydrogeological settings and their percentages in the main geographical regions and in the corresponding whole Mediterranean region as well as in the Mediterranean focus area (250 km zone).</caption>

Information on the size of the countries in this study is taken from the International Hydrogeological Map of Europe (IHME1500)

<table><thead><tr><th rowspan="2">Geographical area</th><th rowspan="2">Total area<sup>a</sup></th><th colspan="2">Karst aquifers in sedimentary and metamorphic carbonate rocks</th><th colspan="2">Karst aquifers in evaporite rocks</th><th colspan="2">Various hydrogeological settings in other sedimentary and volcanic formations</th><th colspan="2">Local, poor and shallow aquifers in other metamorphic and igneous rocks</th></tr><tr><th>1,000 km<sup>2</sup></th><th>1,000 km<sup>2</sup></th><th>%</th><th>1,000 km<sup>2</sup></th><th>%</th><th>1,000 km<sup>2</sup></th><th>%</th><th>1,000 km<sup>2</sup></th><th>%</th></tr></thead><tbody><tr><td>Europe</td><td>2,019</td><td>521</td><td>25.8</td><td>23</td><td>1.1</td><td>1,075</td><td>53.2</td><td>400</td><td>19.9</td></tr><tr><td>Middle East</td><td>1,108</td><td>338</td><td>30.5</td><td>28</td><td>2.5</td><td>634</td><td>57.2</td><td>108</td><td>9.8</td></tr><tr><td>North Africa</td><td>5,537</td><td>1,400</td><td>25.3</td><td>25</td><td>0.5</td><td>3,860</td><td>69.7</td><td>252</td><td>4.5</td></tr><tr><td>Mediterranean region</td><td>8,664</td><td>2,259</td><td>26.1</td><td>76</td><td>0.9</td><td>5,569</td><td>64.2</td><td>761</td><td>8.8</td></tr><tr><td>Mediterranean focus area (250-km zone)</td><td>3,123</td><td>1,233</td><td>39.5</td><td>52</td><td>1.7</td><td>1,541</td><td>49.3</td><td>297</td><td>9.5</td></tr></tbody></table>

that groundwater resources in karst areas are quite significant in the Mediterranean region. The various hydrogeological settings in other sedimentary and volcanic formations cover 49.3%, while local poor and shallow aquifers in other metamorphic and igneous rocks only account for 9.5% (Table 1; Fig. 3).

Looking at the main geographical areas (as defined in Fig. 1, including the entire territory of each country, even if not fully covered by Fig. 2), the proportion of carbonate rocks is the highest in the countries of the Middle East, at 30.5% (Table 1; Fig. 3), although the proportions vary greatly between the individual countries, ranging from 22.0% in Turkey to 72.7% in Jordan. The proportion of karst in European countries reaches 25.8%, with a wide range of proportions when looking at individual countries. Among the four largest European countries, the proportion ranges from 16.7% in Spain, 22.0% in Italy and 26.9% in Greece to 36.9% in France (Table 2). Much higher proportions are found in smaller countries such as Bosnia and Herzegovina (42.0%) or Montenegro (66.1%).

The North African countries, with 25.3%, have the lowest proportion of carbonate rocks of the three major regions, but they also have the largest and, in some cases, the most extensive continuous carbonate areas. Libya tops the list with 31.0% karst, corresponding to an area of 504.9 thousand km². In Egypt, the proportion is 38.9%, with 402.5 thousand km² of karst, and in Algeria, 301.5 thousand km² of karst make up ~13.0% of the total area, although these are not large continuous carbonate areas as in Libya and Egypt, but numerous small areas.

Evaporite rocks are much rarer, occurring in only 2.5% of Middle Eastern countries, 1.1% of European countries and 0.5% of North African countries. Turkey and Spain have by far the largest areas of evaporite rocks, with 20.7 and 20.4 thousand km² respectively, while Tunisia has the highest proportion at 5.5%.

In general, the smaller Mediterranean countries sometimes have relatively high proportions of karst in relation to their total area. Malta stands out as an island state with 100% karst. The Palestinian Territories (West Bank and Gaza Strip) also have a high proportion of karst (78.8%), as do Lebanon

<figure><img src="image_8.png" /><figcaption>**Fig. 3** Hydrogeological settings and their percentages in the main geographical regions (see colored areas in Fig. 1) and in the corresponding whole Mediterranean region as well as in the Mediterranean focus area (250 km zone)</figcaption></figure>

<footnote><sup>a</sup>Water surfaces of lakes not included</footnote>
<table><caption>**Table 2** Size of countries and percentage of different hydrogeological settings. Information on the size of countries in this study is taken from the International Hydrogeological Map of Europe (IHME1500)</caption><thead><tr><th rowspan="2">Country</th><th rowspan="2">Total area<sup>a</sup><br>1,000 km<sup>2</sup></th><th colspan="2">Karst aquifers<br>in sedimentary<br>and metamor-<br>phic carbonate<br>rocks</th><th colspan="2">Karst aquifers<br>in<br>evaporite<br>rocks</th><th colspan="2">Various hydro-<br>geological set-<br>tings in other<br>sedimentary<br>and volcanic<br>formations</th><th colspan="2">Local poor<br>and shallow<br>aquifers in<br>other meta-<br>morphic and<br>igneous rocks</th></tr><tr><th>1,000<br>km<sup>2</sup></th><th>%</th><th>1,000<br>km<sup>2</sup></th><th>%</th><th>1,000<br>km<sup>2</sup></th><th>%</th><th>1,000<br>km<sup>2</sup></th><th>%</th></tr></thead><tbody><tr><td>Albania</td><td>28</td><td>8</td><td>26.9</td><td>0</td><td>0.0</td><td>17</td><td>58.8</td><td>4</td><td>14.3</td></tr><tr><td>Algeria</td><td>2,316</td><td>301</td><td>13.0</td><td>6</td><td>0.3</td><td>1,902</td><td>82.1</td><td>107</td><td>4.6</td></tr><tr><td>Andorra</td><td>0</td><td>0</td><td>2.3</td><td>0</td><td>0.0</td><td>0</td><td>45.5</td><td>0</td><td>52.2</td></tr><tr><td>Austria</td><td>84</td><td>21</td><td>25.4</td><td>0</td><td>0.0</td><td>38</td><td>45.5</td><td>24</td><td>29.1</td></tr><tr><td>Bosnia and Herzegovina</td><td>51</td><td>21</td><td>42.0</td><td>0</td><td>0.0</td><td>24</td><td>47.7</td><td>5</td><td>10.3</td></tr><tr><td>Bulgaria</td><td>113</td><td>32</td><td>28.3</td><td>0</td><td>0.0</td><td>56</td><td>49.3</td><td>25</td><td>22.4</td></tr><tr><td>Croatia</td><td>56</td><td>22</td><td>39.2</td><td>0</td><td>0.0</td><td>34</td><td>59.5</td><td>1</td><td>1.3</td></tr><tr><td>Cyprus</td><td>9</td><td>2</td><td>24.9</td><td>0</td><td>0.0</td><td>5</td><td>55.0</td><td>2</td><td>20.1</td></tr><tr><td>Egypt</td><td>1,034</td><td>402</td><td>38.9</td><td>9</td><td>0.8</td><td>534</td><td>51.7</td><td>89</td><td>8.6</td></tr><tr><td>France</td><td>561</td><td>207</td><td>36.9</td><td>0</td><td>0.0</td><td>248</td><td>44.2</td><td>106</td><td>18.9</td></tr><tr><td>Greece</td><td>28</td><td>8</td><td>26.9</td><td>0</td><td>0.0</td><td>17</td><td>58.8</td><td>4</td><td>14.3</td></tr><tr><td>Israel</td><td>22</td><td>10</td><td>44.9</td><td>0</td><td>0.0</td><td>12</td><td>54.7</td><td>0</td><td>0.3</td></tr><tr><td>Italy</td><td>302</td><td>66</td><td>22.0</td><td>2</td><td>0.8</td><td>192</td><td>63.7</td><td>41</td><td>13.6</td></tr><tr><td>Jordan</td><td>89</td><td>64</td><td>72.7</td><td>0</td><td>0.0</td><td>23</td><td>25.9</td><td>1</td><td>1.4</td></tr><tr><td>Lebanon</td><td>10</td><td>7</td><td>68.0</td><td>0</td><td>0.0</td><td>3</td><td>32.0</td><td>0</td><td>0.0</td></tr><tr><td>Libya</td><td>1,626</td><td>505</td><td>31.0</td><td>0</td><td>0.0</td><td>1,113</td><td>68.4</td><td>8</td><td>0.5</td></tr><tr><td>Malta</td><td>0</td><td>0</td><td>100.0</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td></tr><tr><td>Montenegro</td><td>14</td><td>9</td><td>66.1</td><td>0</td><td>0.0</td><td>5</td><td>32.8</td><td>0</td><td>1.1</td></tr><tr><td>Morocco</td><td>406</td><td>152</td><td>37.4</td><td>2</td><td>0.5</td><td>204</td><td>50.3</td><td>48</td><td>11.8</td></tr><tr><td>North Macedonia</td><td>25</td><td>4</td><td>14.0</td><td>0</td><td>0.0</td><td>10</td><td>41.3</td><td>11</td><td>44.7</td></tr><tr><td>Palestinian Territories</td><td>6</td><td>5</td><td>78.8</td><td>0</td><td>0.0</td><td>1</td><td>21.2</td><td>0</td><td>0.0</td></tr><tr><td>Portugal</td><td>90</td><td>3</td><td>3.5</td><td>0</td><td>0.0</td><td>42</td><td>47.3</td><td>44</td><td>49.2</td></tr><tr><td>Serbia</td><td>88</td><td>11</td><td>12.1</td><td>0</td><td>0.0</td><td>62</td><td>69.8</td><td>16</td><td>18.1</td></tr><tr><td>Slovenia</td><td>20</td><td>10</td><td>47.3</td><td>0</td><td>0.0</td><td>10</td><td>50.4</td><td>1</td><td>4.3</td></tr><tr><td>Spain</td><td>506</td><td>84</td><td>16.7</td><td>20</td><td>4.0</td><td>296</td><td>58.5</td><td>105</td><td>20.8</td></tr><tr><td>Switzerland</td><td>41</td><td>13</td><td>30.3</td><td>0</td><td>0.0</td><td>19</td><td>45.4</td><td>10</td><td>24.2</td></tr><tr><td>Syria</td><td>208</td><td>82</td><td>39.4</td><td>7</td><td>3.3</td><td>119</td><td>57.3</td><td>0</td><td>0.0</td></tr><tr><td>Tunisia</td><td>155</td><td>39</td><td>25.2</td><td>8</td><td>5.5</td><td>108</td><td>69.3</td><td>0</td><td>0.0</td></tr><tr><td>Turkey</td><td>774</td><td>171</td><td>22.0</td><td>21</td><td>2.7</td><td>476</td><td>61.5</td><td>107</td><td>13.8</td></tr></tbody></table>

(68.0%) and Montenegro (66.1%). The proportion of outcropping karst areas below 10% can be found only in Portugal with 3.5%; however, this country is mostly outside the 250-km area of focus. In many countries, the exposed karst does not always represent the actual regional karst groundwater resources, as karst often extends underground (Fig. 2).

The results show that the Mediterranean region, with 2.26 million km² of karst or 39.5% of its area, is one of the most karstified regions compared to other large countries and regions of the world. Goldscheider et al. (2020) show in the WOKAM map that China has a slightly larger karst area of 2.54 million km², corresponding to 26.5% of its total land area, while Europe has ‘only’ 21.8%, North Africa

19.6% and Asia 18.6%. However, the proportions of karst in MEDKAM and WOKAM differ because of the different observation and classification scale (WOKAM distinguishes between continuous and discontinuous karst) and partly also because of the different lithological basis.

The presence of karst aquifers and their percentage of the total area can be an indication of the presence of karst groundwater resources and the dependency of humans and the environment on these water resources. However, climatic conditions, i.e. groundwater recharge (see section “Karst groundwater recharge”), as well as the thickness, functioning and storage capacity of karst aquifers, determine the actual availability of these karst groundwater resources.

<footnote><sup>a</sup>Water surfaces of lakes not included</footnote>
<table><caption>**Table 3** Length of the coastline of individual states, separated into continental coastline and islands, with the corresponding proportion of coastal karst. Information on the length of coastlines in this study is taken from the International Hydrogeological Map of Europe (IHME1500)</caption><thead><tr><th rowspan="2">Country</th><th rowspan="2">Length of coastline including islands</th><th colspan="2">Length of coastal karst, including islands</th><th rowspan="2">Length of continental coastline</th><th colspan="2">Length of continental coastal karst</th></tr><tr><th>km</th><th>%</th><th>km</th><th>%</th></tr></thead><tbody><tr><td>Albania</td><td>512</td><td>201</td><td>39.2</td><td>498</td><td>190</td><td>38.2</td></tr><tr><td>Algeria</td><td>1,296</td><td>399</td><td>30.8</td><td>1,288</td><td>399</td><td>31.0</td></tr><tr><td>Bosnia and Herzegovina</td><td>20</td><td>20</td><td>100.0</td><td>20</td><td>20</td><td>100.0</td></tr><tr><td>Croatia</td><td>3,982</td><td>3,491</td><td>87.7</td><td>1,492</td><td>1,156</td><td>77.5</td></tr><tr><td>Cyprus</td><td>646</td><td>63</td><td>9.8</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Egypt</td><td>2,216</td><td>205</td><td>9.2</td><td>2,401</td><td>205</td><td>8.5</td></tr><tr><td>France</td><td>1,765</td><td>182</td><td>10.3</td><td>919</td><td>161</td><td>17.5</td></tr><tr><td>Greece</td><td>11,712</td><td>5,016</td><td>42.8</td><td>4,129</td><td>1,246</td><td>30.2</td></tr><tr><td>Israel</td><td>189</td><td>9</td><td>5.0</td><td>189</td><td>9</td><td>5.0</td></tr><tr><td>Italy</td><td>6,971</td><td>1,367</td><td>19.6</td><td>3,765</td><td>792</td><td>21.0</td></tr><tr><td>Lebanon</td><td>225</td><td>118</td><td>52.4</td><td>224</td><td>118</td><td>52.7</td></tr><tr><td>Libya</td><td>1,835</td><td>509</td><td>27.7</td><td>1,831</td><td>504</td><td>27.5</td></tr><tr><td>Malta</td><td>134</td><td>134</td><td>100.0</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Montenegro</td><td>218</td><td>148</td><td>67.9</td><td>214</td><td>144</td><td>67.3</td></tr><tr><td>Morocco</td><td>461</td><td>122</td><td>26.5</td><td>461</td><td>122</td><td>26.5</td></tr><tr><td>Palestinian Territories</td><td>40</td><td>0</td><td>0.0</td><td>40</td><td>0</td><td>0.0</td></tr><tr><td>Slovenia</td><td>35</td><td>0</td><td>0.0</td><td>35</td><td>0</td><td>0.0</td></tr><tr><td>Spain</td><td>2,556</td><td>675</td><td>26.4</td><td>1,711</td><td>143</td><td>8.4</td></tr><tr><td>Syria</td><td>192</td><td>46</td><td>24.0</td><td>192</td><td>46</td><td>24.0</td></tr><tr><td>Tunisia</td><td>1,619</td><td>86</td><td>5.3</td><td>1,325</td><td>86</td><td>6.5</td></tr><tr><td>Turkey</td><td>5,107</td><td>1,219</td><td>23.9</td><td>3,990</td><td>1,069</td><td>26.8</td></tr><tr><td>Total</td><td>41,731</td><td>14,010</td><td>33.6</td><td>24,725</td><td>6,409</td><td>25.9</td></tr></tbody></table>

## Mediterranean coastal karst

Data on the length of coastlines vary considerably depending on the scale and source of the data. Therefore, information on country size and coastline length from other data sources may differ from the data used in this study (IHME1500), which is based on the NATO VMAP data (NIMA 2001). The total length of the continental coastline and mapped islands is ~41,600 km, of which ~24,700 km is the continental coastline and ~16,900 km is the cumulative coastline of the numerous Mediterranean islands. To assess the occurrence of both carbonate and evaporite rocks on the coast, the term coastal karst is used here to encompass both hydrogeological units.

The presence of these rocks on the coast usually implies the presence of submarine karst springs or caves, which makes these coastal sections very vulnerable to saltwater intrusion (Fleury et al. 2023). The coastal karst is, therefore, particularly worthy of protection and requires correspondingly adapted management, especially in the case of islands. The statistical analysis shows that 33.6% of the total Mediterranean coastline, including the islands, is karst, corresponding to a length of ~14,000 km. Of this, the continental coastline alone accounts for ~6,400 km of karst, or 25.9%

of the total coastline (Table 3). The country with the longest coastline, including the numerous islands, is Greece with a total length of ~11,700 km, of which 42.8% is karst. If only the continental coastline of ~4,100 km is considered, 30.2% is karst. Italy follows with a total coastline of ~6,970 km, of which 19.6% is karst, and Turkey with ~5,100 km total coastline and 23.9% of karst (Table 3). Countries with a very high proportion of coastal karst are Croatia with 87.7%, Montenegro with 67.9% and Lebanon with 52.4%. Some smaller countries or countries with a short coastline, such as Malta and Bosnia and Herzegovina have a karst coverage of 100% in the studied mapping scale and resolution.

Many of the larger Mediterranean islands have a high proportion of coastal karst, such as the Balearic Islands: 43.5% of the coast of Ibiza consists of karst, Menorca has 60.7% and Mallorca has 78.0% coastal karst. The islands of the Dinarides, in particular, have a high proportion of coastal karst, with some islands consisting almost entirely of carbonate rock, such as Cres, Hvar and Dugi Otok (Fig. 4b) or Korčula or Brač (Table 4). The many Greek islands are also strongly characterised by karst, and some have a high proportion of coastal karst such as Kefalonia (87.8%), Lefkada (81.9%) or Zakynthos (60.7%; Fig. 4c). The largest Greek island, Crete, also has a high proportion of its coastline covered by karst, at 65.1% (Fig. 4d).
<figure><img src="image_9.png" /><figcaption>**Fig. 4** Selected examples of Mediterranean islands with important absolute or relative occurrence of coastal karst: a the Balearic Islands Ibiza, Mallorca and Menorca, b islands along the northern part of the</figcaption></figure>

Dinaric coast of Croatia, c the southern Italian island Sicily, d the Greek islands Crete, e Ionian islands Lefkada, Kafelonia and Zakynthos

## Karst springs

Some of the world's largest springs in the Mediterranean region occur in the Dinarides, such as the Buna spring (Q12 in Table 5) and Vrela Trebišnjice spring (Q49 in Table 5) in Bosnia and

Herzegovina, with maximum discharge of 380 and 219 m³/s, respectively. Springs in the eastern Mediterranean also show discharge maxima of several tens of cubic meters per second, such as the Afqa spring (Q1) in Lebanon with 65.5 m³/s, which is partly used for local domestic water supply (Schuler and
<table><caption>**Table 4** Selection of the larger Mediterranean islands that are not independent states, with indication of the length of the coastline and the proportion of karst. Cyprus and Malta are included in Table 3. Information on the size of countries and the length of coastlines in this study is taken from the International Hydrogeological Map of Europe (IHME1500)</caption><thead><tr><th>Name</th><th>Country</th><th>Length of coastline (km)</th><th>Length of coastal karst (km)</th><th>Percentage of karst on coastline (%)</th></tr></thead><tbody><tr><td>Brač</td><td>Croatia</td><td>113.5</td><td>113.5</td><td>100.0</td></tr><tr><td>Chios</td><td>Greece</td><td>173.4</td><td>101.6</td><td>58.6</td></tr><tr><td>Corse</td><td>France</td><td>763.9</td><td>11.8</td><td>1.5</td></tr><tr><td>Cres</td><td>Croatia</td><td>188.5</td><td>188.5</td><td>100.0</td></tr><tr><td>Crete</td><td>Greece</td><td>845.7</td><td>550.3</td><td>65.1</td></tr><tr><td>Dugi Otok</td><td>Croatia</td><td>134.1</td><td>134.1</td><td>100.0</td></tr><tr><td>Euboea</td><td>Greece</td><td>582.1</td><td>274.6</td><td>47.2</td></tr><tr><td>Hvar</td><td>Croatia</td><td>177.1</td><td>177.1</td><td>100.0</td></tr><tr><td>Ibiza</td><td>Spain</td><td>138.7</td><td>60.3</td><td>43.5</td></tr><tr><td>Kefalonia</td><td>Greece</td><td>216.5</td><td>190.0</td><td>87.8</td></tr><tr><td>Korčula</td><td>Croatia</td><td>122.6</td><td>122.6</td><td>100.0</td></tr><tr><td>Korfu</td><td>Greece</td><td>198.6</td><td>159.5</td><td>80.3</td></tr><tr><td>Krk</td><td>Croatia</td><td>137.2</td><td>103.2</td><td>75.2</td></tr><tr><td>Lefkada</td><td>Greece</td><td>103.2</td><td>84.5</td><td>81.9</td></tr><tr><td>Lesbos</td><td>Greece</td><td>306.9</td><td>40.5</td><td>13.2</td></tr><tr><td>Mallorca</td><td>Spain</td><td>389.9</td><td>304.2</td><td>78.0</td></tr><tr><td>Menorca</td><td>Spain</td><td>171.7</td><td>104.2</td><td>60.7</td></tr><tr><td>Pag</td><td>Croatia</td><td>209.2</td><td>137.8</td><td>65.9</td></tr><tr><td>Rab</td><td>Croatia</td><td>73.3</td><td>45.0</td><td>61.4</td></tr><tr><td>Rhodos</td><td>Greece</td><td>205.3</td><td>42.0</td><td>20.4</td></tr><tr><td>Sardegna</td><td>Italy</td><td>1,180.3</td><td>166.5</td><td>14.1</td></tr><tr><td>Sicilia</td><td>Italy</td><td>990.5</td><td>265.2</td><td>26.8</td></tr><tr><td>Zakynthos</td><td>Greece</td><td>118.7</td><td>71.4</td><td>60.2</td></tr></tbody></table>

Margane 2013), or the Kırkgöz spring (Q31) in Turkey with 62.8 m³/s. Many of the springs play an important role in regional water supply, but are also life-giving for unique ecosystems—for example, Montpellier is supplied by the Lez spring (Q37; Bakalowicz 2011), which hosts an important ecosystem and feeds the Lez River with a maximum discharge of 11.9 m³/s (Siegel et al. 2023). An example of a large karst spring on an island is the Almyros of Heraklion (Q5) in Crete with a maximum discharge of 30.0 m³/s. However, the countless small karst springs around the Mediterranean, which are only partially documented in this study, also serve the local population for domestic, agricultural, livestock and energy purposes in large parts of the Mediterranean region. Examples are the Maqar spring (Q39), one of many in the Al Jabal Al Akhdar region of northern Libya (Hamad and El Hasia 2022), or the Ain El Gudeirat spring (Q3) in the north-eastern Sinai region of Egypt, which also provides water for irrigation and local water supply (LaMoreaux and Tanner 2001).

## Submarine karst springs

A selection of 39 submarine springs in the Mediterranean are presented (Table 6), most of which do not have information on their discharge due to the nature of their submarine location. For this reason, information is only available for a few springs such as for the Chekka spring (S15) in Lebanon, where a maximum discharge of 60 m³/s was measured, or for the Source de Port-Miou (S29) in France, where a maximum discharge of 50 m³/s was recorded. Additionally, the Moraig and Toix spring (S23) in Spain has a maximum discharge of 9.0 m³/s, while the Uji Ftohte (S38) in Albania measured 4.3 m³/s. A well-known and studied spring is Source de la Vise (S36) in France, which still has a maximum discharge of 0.45m³/s. These examples show that submarine springs can indeed release significant amounts of freshwater and, therefore, represent a potentially usable water resource despite the technical challenges. A special feature is the Argostoli system on the island of Kafelonia (Fig. 4e), where a submarine swallow hole infiltrates the seawater on the west coast and discharges ~2 weeks later ~15 km away on the east coast, partly through coastal springs but also through a submarine spring (proven by tracer test).

## Karst caves

Table 7 gives an overview of 30 selected caves from the MEDKAM database that represent the diversity and size of the known caves in the Mediterranean focus area. The longest is the Schönberg cave system (C26) in the Alps, which is 140 km long and has a remarkable depth of 1,061 m. It is mentioned here because it is still within the 250-km zone, but is not strictly speaking part of the classic Mediterranean. Other large cave systems include the Hammam Trozza (C13) in Tunisia and the Coume Ouarnède cave system (C4) in France, which are both over 100 km long. Other deep caves include the Čehi 2 cave in Slovenia (C3) with 1,502 m depth, the Clot d'Aspres cave system (C4) in France with a depth of 1,066 m and the Željezna Jama cave (C30) in Montenegro, which is 1,027 m deep. An example of a cave on an island is the Cuevas del Drach (C7, Fig. 4a) on Mallorca, which is located on the south-eastern side and is ~1.2 km long.

## Karst groundwater recharge

Karst groundwater recharge shows a clear trend from north to south (Fig. 5). Countries in the Alpine region and parts of the Dinarides have the highest average karst groundwater recharge of more than 800 mm/year in Switzerland, Austria and Slovenia. In Croatia, Bosnia and Herzegovina and Albania, the recharge is still between 500 and 800 mm/year. Lower average recharge rates of 300 to 500 mm/year are found in Italy, Serbia, France, Portugal and Northern
<table><caption>**Table 5** Selected springs from the MEDKAM database with estimations for low-flow and high-flow discharge (m³/s)</caption><thead><tr><th>ID</th><th>Name</th><th>Country</th><th>Longitude</th><th>Latitude</th><th>Low (m³/s)</th><th>High (m³/s)</th></tr></thead><tbody><tr><td>Q1</td><td>Afqa</td><td>Lebanon</td><td>35.893600</td><td>34.067550</td><td>0.1</td><td>64.5</td></tr><tr><td>Q2</td><td>Ain Al Laf'a</td><td>Morocco</td><td>-6.462580</td><td>33.847490</td><td>0.2<sup>b</sup></td><td>ND</td></tr><tr><td>Q3</td><td>Ain El Gudeirat</td><td>Egypt</td><td>34.404390</td><td>30.640440</td><td>0.02</td><td>0.2</td></tr><tr><td>Q4</td><td>Ain Ribaa</td><td>Morocco</td><td>-5.232118</td><td>33.757743</td><td>1.0<sup>b</sup></td><td>ND</td></tr><tr><td>Q5</td><td>Almyros of Heraklion</td><td>Greece</td><td>24.872480</td><td>35.410210</td><td>3.3</td><td>30</td></tr><tr><td>Q6</td><td>Anjar</td><td>Lebanon</td><td>35.939060</td><td>33.775486</td><td>0.48</td><td>9.2</td></tr><tr><td>Q7<sup>a</sup></td><td>Ayn Zayanah</td><td>Libya</td><td>20.159586</td><td>32.211600</td><td>5.0<sup>b</sup></td><td>ND</td></tr><tr><td>Q8</td><td>Basso Bussento group</td><td>Italy</td><td>15.540000</td><td>40.137200</td><td>0.2</td><td>60</td></tr><tr><td>Q9</td><td>Bistrica (Syri I Kalter)</td><td>Albania</td><td>20.195200</td><td>39.926440</td><td>13.5</td><td>36</td></tr><tr><td>Q10</td><td>Boka</td><td>Slovenia</td><td>13.485874</td><td>46.321800</td><td>0</td><td>100</td></tr><tr><td>Q11</td><td>Bournillon (Source du)</td><td>France</td><td>5.434349</td><td>45.054841</td><td>0</td><td>80</td></tr><tr><td>Q12</td><td>Buna</td><td>Bosnia and Herzegovina</td><td>17.908740</td><td>43.257010</td><td>3</td><td>380</td></tr><tr><td>Q13<sup>a</sup></td><td>Capo Pescara springs</td><td>Italy</td><td>13.821573</td><td>42.163795</td><td>5.5</td><td>7.5</td></tr><tr><td>Q14</td><td>Divje jezero</td><td>Slovenia</td><td>14.032789</td><td>45.982689</td><td>ND</td><td>60</td></tr><tr><td>Q15</td><td>Düdenbaşı</td><td>Turkey</td><td>33.726680</td><td>36.952010</td><td>0.03</td><td>16.7</td></tr><tr><td>Q16</td><td>Dumanlı<sup>c</sup></td><td>Turkey</td><td>31.553452</td><td>36.921409</td><td>25</td><td>50</td></tr><tr><td>Q17</td><td>Figeh</td><td>Syria</td><td>36.161543</td><td>33.611977</td><td>2.1</td><td>23.4</td></tr><tr><td>Q18</td><td>Font Estramar</td><td>France</td><td>2.958325</td><td>42.859012</td><td>0.8</td><td>20</td></tr><tr><td>Q19</td><td>Fontaine-de-Vaucluse</td><td>France</td><td>5.132778</td><td>43.917778</td><td>3.1</td><td>81.4</td></tr><tr><td>Q20<sup>*</sup></td><td>Fontaine de Nîmes</td><td>France</td><td>4.348898</td><td>43.840438</td><td>0.09</td><td>18.38</td></tr><tr><td>Q21<sup>a</sup></td><td>Foux de la Vis</td><td>France</td><td>3.483763</td><td>43.900002</td><td>1.2</td><td>> 30</td></tr><tr><td>Q22</td><td>Gillardes (Source des)</td><td>France</td><td>5.888000</td><td>44.759500</td><td>0.5</td><td>60</td></tr><tr><td>Q23</td><td>Glava Zete</td><td>Montenegro</td><td>18.996960</td><td>42.675130</td><td>2</td><td>100</td></tr><tr><td>Q24</td><td>Gökova springs</td><td>Turkey</td><td>28.426054</td><td>37.037400</td><td>10.2</td><td>13.5</td></tr><tr><td>Q25</td><td>Gökpınar. Mugla</td><td>Turkey</td><td>38.036874</td><td>37.024109</td><td>6.2</td><td>15</td></tr><tr><td>Q26</td><td>Gradac</td><td>Serbia</td><td>19.856170</td><td>44.210440</td><td>0.4</td><td>10</td></tr><tr><td>Q27</td><td>Hubelj</td><td>Slovenia</td><td>13.918356</td><td>45.899864</td><td>0.04</td><td>60</td></tr><tr><td>Q28</td><td>Izvor Krke</td><td>Croatia</td><td>16.240090</td><td>44.042570</td><td>1.2</td><td>100</td></tr><tr><td>Q29</td><td>Izvor Kupe</td><td>Croatia</td><td>14.694140</td><td>45.491510</td><td>1.1</td><td>100</td></tr><tr><td>Q30</td><td>Izvori Cetine</td><td>Croatia</td><td>16.437270</td><td>43.965270</td><td>0.6</td><td>60</td></tr><tr><td>Q31</td><td>Kamniška Bistrica</td><td>Slovenia</td><td>14.594674</td><td>46.327850</td><td>2.1</td><td>62.6</td></tr><tr><td>Q32</td><td>Kırkgöz</td><td>Turkey</td><td>33.582590</td><td>37.058190</td><td>7.1</td><td>62.8</td></tr><tr><td>Q33</td><td>Klokot</td><td>Bosnia and Herzegovina</td><td>15.806750</td><td>44.824030</td><td>2.3</td><td>70</td></tr><tr><td>Q34</td><td>Krka</td><td>Slovenia</td><td>14.776623</td><td>45.889899</td><td>0.8</td><td>80</td></tr><tr><td>Q35</td><td>Kumanitsa</td><td>Bulgaria</td><td>24.708350</td><td>42.753750</td><td>0.5</td><td>3</td></tr><tr><td>Q36<sup>a</sup></td><td>Lez (Source du)</td><td>France</td><td>3.844122</td><td>43.718235</td><td>0</td><td>11.9</td></tr><tr><td>Q37</td><td>Ljubljania</td><td>Slovenia</td><td>14.300385</td><td>45.950265</td><td>1.5</td><td>120</td></tr><tr><td>Q38</td><td>Mareza</td><td>Montenegro</td><td>19.182530</td><td>42.479920</td><td>2.0</td><td>20</td></tr><tr><td>Q39</td><td>Mqar</td><td>Libya</td><td>22.520500</td><td>32.623000</td><td>ND</td><td>ND</td></tr><tr><td>Q40<sup>a</sup></td><td>Nacimiento del rio Mundo</td><td>Spain</td><td>-2.438531</td><td>38.455362</td><td>0</td><td>100</td></tr><tr><td>Q41</td><td>Olukköprü spring</td><td>Turkey</td><td>34.182060</td><td>37.179670</td><td>7</td><td>30</td></tr><tr><td>Q42</td><td>Ombla</td><td>Croatia</td><td>18.141550</td><td>42.675570</td><td>4</td><td>138</td></tr><tr><td>Q43</td><td>Ouyisse (Source de l')</td><td>France</td><td>1.584203</td><td>44.790271</td><td>0.5</td><td>150</td></tr><tr><td>Q44</td><td>Riječina</td><td>Croatia</td><td>14.427135</td><td>45.425682</td><td>0</td><td>150</td></tr><tr><td>Q45</td><td>Rižana</td><td>Slovenia</td><td>13.890424</td><td>45.528058</td><td>0.3</td><td>80</td></tr><tr><td>Q46</td><td>Ses Font Ufanes</td><td>Spain</td><td>2.961845</td><td>39.802997</td><td>0</td><td>30</td></tr><tr><td>Q47</td><td>Soča</td><td>Slovenia</td><td>13.728532</td><td>46.412272</td><td>1.1</td><td>100</td></tr><tr><td>Q48</td><td>Unica</td><td>Slovenia</td><td>14.250353</td><td>45.820185</td><td>0.2</td><td>90</td></tr><tr><td>Q49</td><td>Vrela Trebišnjice<sup>d</sup></td><td>Bosnia and Herzegovina</td><td>18.428360</td><td>42.853590</td><td>2.5</td><td>219</td></tr><tr><td>Q50</td><td>Zarka spring (El Assi)</td><td>Lebanon</td><td>35.944326</td><td>34.136104</td><td>6.9</td><td>24.3</td></tr></tbody></table>

ND not defined

<footnote><sup>a</sup>Connected to an important karst groundwater-dependent ecosystem (Siegel et al. 2023)</footnote>

<footnote><sup>b</sup>Only mean discharge was available and classified as “low” discharge</footnote>

<footnote><sup>c</sup>Formerly one of the largest karst springs in the world, but submerged in a reservoir since 1982</footnote>

<footnote><sup>d</sup>These springs were submerged at the end of the 1960s during the construction of the Bileća reservoir and are its main water source</footnote>
<table><caption>**Table 6** Submarine karst springs with estimations for low-flow and high-flow discharge (m³/s)</caption><thead><tr><th>ID</th><th>Name</th><th>Country</th><th>Longitude</th><th>Latitude</th><th>Low</th><th>High</th></tr></thead><tbody><tr><td>S1</td><td>Aguadulce</td><td>Spain</td><td>-2.249925</td><td>36.792473</td><td>ND</td><td>ND</td></tr><tr><td>S2</td><td>Aigua Dolç and la Falconera</td><td>Spain</td><td>1.849875</td><td>41.215732</td><td>ND</td><td>ND</td></tr><tr><td>S3</td><td>Anavalos Kiveri</td><td>Greece</td><td>22.737265</td><td>37.518283</td><td>ND</td><td>ND</td></tr><tr><td>S4</td><td>Antalya</td><td>Turkey</td><td>30.6788347</td><td>36.7905208</td><td>ND</td><td>ND</td></tr><tr><td>S5</td><td>Antheron</td><td>Greece</td><td>23.51455</td><td>38.482027</td><td>ND</td><td>ND</td></tr><tr><td>S6</td><td>Arbanija and Slanita</td><td>Croatia</td><td>16.325141</td><td>43.519898</td><td>ND</td><td>ND</td></tr><tr><td>S7<sup>a,b</sup></td><td>Argostoli (a=swallow hole, b=spring)</td><td>Greece</td><td>20.470881</td><td>38.175685</td><td>ND</td><td>ND</td></tr><tr><td>S8</td><td>Bali</td><td>Greece</td><td>24.78715</td><td>35.41252</td><td>ND</td><td>ND</td></tr><tr><td>S9</td><td>Benissa and Calpe</td><td>Spain</td><td>0.096236</td><td>38.666637</td><td>ND</td><td>ND</td></tr><tr><td>S10</td><td>Brač Island</td><td>Croatia</td><td>16.458537</td><td>43.294363</td><td>ND</td><td>ND</td></tr><tr><td>S11</td><td>Bue Marino</td><td>Italy</td><td>9.623155</td><td>40.246866</td><td>ND</td><td>ND</td></tr><tr><td>S12</td><td>Cannes</td><td>France</td><td>7.084438</td><td>43.542896</td><td>ND</td><td>ND</td></tr><tr><td>S13</td><td>Castellamare</td><td>Italy</td><td>12.877892</td><td>38.049807</td><td>ND</td><td>ND</td></tr><tr><td>S14</td><td>Cefalù</td><td>Italy</td><td>14.044757</td><td>38.037825</td><td>ND</td><td>ND</td></tr><tr><td>S15</td><td>Chekka</td><td>Lebanon</td><td>35.72068</td><td>34.329819</td><td>0.2</td><td>ND</td></tr><tr><td>S16</td><td>Chidro</td><td>Italy</td><td>17.683414</td><td>40.303202</td><td>2.5<sup>b</sup></td><td>ND</td></tr><tr><td>S17</td><td>Denia</td><td>Spain</td><td>0.129757</td><td>38.841706</td><td>ND</td><td>ND</td></tr><tr><td>S18</td><td>Gaeta</td><td>Italy</td><td>13.569334</td><td>41.231956</td><td>ND</td><td>ND</td></tr><tr><td>S19</td><td>Gurdić and Škurda</td><td>Montenegro</td><td>18.745531</td><td>42.477181</td><td>ND</td><td>ND</td></tr><tr><td>S20</td><td>Latakia</td><td>Syria</td><td>35.787184</td><td>35.48411</td><td>ND</td><td>ND</td></tr><tr><td>S21</td><td>Lera Pass</td><td>Albania</td><td>19.611853</td><td>40.148594</td><td>0.4<sup>b</sup></td><td>ND</td></tr><tr><td>S22</td><td>Makarska</td><td>Croatia</td><td>17.009245</td><td>43.297702</td><td>ND</td><td>ND</td></tr><tr><td>S23</td><td>Moraig and Toix</td><td>Spain</td><td>0.168832</td><td>38.709284</td><td>0.3</td><td>9</td></tr><tr><td>S24</td><td>Mortola</td><td>Italy</td><td>7.554824</td><td>43.777485</td><td>0.02</td><td>0.2</td></tr><tr><td>S25</td><td>Nabrežina</td><td>Italy</td><td>13.6672236</td><td>45.7404762</td><td>ND</td><td>ND</td></tr><tr><td>S26</td><td>Novi Vinodolski</td><td>Croatia</td><td>14.794591</td><td>45.122805</td><td>ND</td><td>ND</td></tr><tr><td>S27</td><td>Ovacik</td><td>Turkey</td><td>33.198631</td><td>36.126952</td><td>1.0<sup>b</sup></td><td>ND</td></tr><tr><td>S28</td><td>Polla di Cadimare</td><td>Italy</td><td>9.871071</td><td>44.062526</td><td>ND</td><td>ND</td></tr><tr><td>S29</td><td>Port-Miou (Source sous-marine de)</td><td>France</td><td>5.513319</td><td>43.205184</td><td>3.0</td><td>50</td></tr><tr><td>S30</td><td>Pozzalo</td><td>Italy</td><td>14.84062</td><td>36.71742</td><td>ND</td><td>ND</td></tr><tr><td>S31</td><td>Roquebrune-Cap Martin</td><td>France</td><td>7.485171</td><td>43.746165</td><td>0.1</td><td>0.2</td></tr><tr><td>S32</td><td>Rosas springs</td><td>Spain</td><td>3.242482</td><td>42.222494</td><td>ND</td><td>ND</td></tr><tr><td>S33</td><td>San Remo</td><td>Italy</td><td>7.787234</td><td>43.815296</td><td>ND</td><td>ND</td></tr><tr><td>S34</td><td>Senj</td><td>Croatia</td><td>14.897833</td><td>44.991176</td><td>ND</td><td>ND</td></tr><tr><td>S35</td><td>Šibenik</td><td>Croatia</td><td>15.839243</td><td>43.717598</td><td>ND</td><td>ND</td></tr><tr><td>S36</td><td>Source de la Vise</td><td>France</td><td>3.670188</td><td>43.440172</td><td>&lt;0</td><td>0.45</td></tr><tr><td>S37</td><td>Taranto</td><td>Italy</td><td>17.223062</td><td>40.448685</td><td>ND</td><td>ND</td></tr><tr><td>S38</td><td>Uji Ftohte</td><td>Albania</td><td>19.439414</td><td>40.385471</td><td>1.8</td><td>4.3</td></tr><tr><td>S39</td><td>Volovsko</td><td>Croatia</td><td>14.318947</td><td>45.342365</td><td>ND</td><td>ND</td></tr></tbody></table>

ND = not defined; <0 = source of estavelle type that functions as a sink for marine water in low-flow conditions

Macedonia, while recharge rates of less than 100 mm/year are found in the MENA region and the Middle East, with Libya and Egypt having the lowest average recharge rates of 13.5 and 1.6 mm/year respectively, apparently due to the

extensive desert regions in their central and southern parts (Table 8).

The average annual karst recharge volume in the Mediterranean region is 740,339 million m³, while within the

<footnote>a a = swallow hole, b = spring</footnote>

<footnote>b Only one discharge value was available and classified as “low” discharge</footnote>
<table><caption>**Table 7** Selected karst caves from the MEDKAM database of the Mediterranean region</caption><thead><tr><th>ID</th><th>Name</th><th>Country</th><th>Longitude</th><th>Latitude</th><th>Length (km)</th><th>Depth (m)</th></tr></thead><tbody><tr><td>C1<sup>a</sup></td><td>Ayalon</td><td>Israel</td><td>34.92798</td><td>31.91032</td><td>2.7</td><td>100</td></tr><tr><td>C2</td><td>Ayvaini Cave</td><td>Turkey</td><td>28.42250</td><td>40.07350</td><td>5.5</td><td>ND</td></tr><tr><td>C3</td><td>Čehi 2</td><td>Slovenia</td><td>13.51470</td><td>46.36730</td><td>5.5</td><td>1,502</td></tr><tr><td>C4</td><td>Clot d'Aspres C.S</td><td>France</td><td>5.57460</td><td>45.02210</td><td>40</td><td>1,066</td></tr><tr><td>C5</td><td>Coume Ouarnède C.S</td><td>France</td><td>0.86640</td><td>42.97250</td><td>105.8</td><td>975</td></tr><tr><td>C6<sup>a</sup></td><td>Cueva de Nerja</td><td>Spain</td><td>-3.84483</td><td>36.76199</td><td>4.8</td><td>158</td></tr><tr><td>C7<sup>a</sup></td><td>Cueva del Gato</td><td>Spain</td><td>-5.23820</td><td>36.72721</td><td>4.5</td><td>ND</td></tr><tr><td>C8<sup>a</sup></td><td>Cuevas del Drach</td><td>Spain</td><td>3.33028</td><td>39.53583</td><td>1.2</td><td>25</td></tr><tr><td>C9</td><td>Dupnisa Cave</td><td>Turkey</td><td>27.33220</td><td>41.50260</td><td>3.2</td><td>ND</td></tr><tr><td>C10<sup>a</sup></td><td>Frasassi caves</td><td>Italy</td><td>12.96195</td><td>43.40033</td><td>20</td><td>ND</td></tr><tr><td>C11<sup>a</sup></td><td>Friouatto</td><td>Morocco</td><td>-4.07234</td><td>34.10468</td><td>0.27</td><td>ND</td></tr><tr><td>C12<sup>a</sup></td><td>Grotta Zinzulusa</td><td>Italy</td><td>18.43074</td><td>40.01204</td><td>0.26</td><td>ND</td></tr><tr><td>C13<sup>a</sup></td><td>Grotte ain Dhab</td><td>Tunisia</td><td>9.48143</td><td>35.90235</td><td>3</td><td>ND</td></tr><tr><td>C14</td><td>Hammam Trozza</td><td>Tunisia</td><td>9.58000</td><td>35.55000</td><td>120.4</td><td>23</td></tr><tr><td>C15</td><td>İnsuyu Cave</td><td>Turkey</td><td>30.22270</td><td>37.39340</td><td>8.1</td><td>ND</td></tr><tr><td>C16</td><td>Jeita Grotto</td><td>Lebanon</td><td>35.38289</td><td>33.56362</td><td>9</td><td>ND</td></tr><tr><td>C17<sup>a</sup></td><td>Karst en yesos de Sorbas</td><td>Spain</td><td>-2.10807</td><td>37.09194</td><td>50</td><td>ND</td></tr><tr><td>C18</td><td>Kef el Kaous</td><td>Algeria</td><td>1.38000</td><td>35.12000</td><td>4.1</td><td>ND</td></tr><tr><td>C19</td><td>Kef Toghobeit</td><td>Algeria</td><td>5.08270</td><td>35.05200</td><td>3.9</td><td>722</td></tr><tr><td>C20</td><td>Kita Gaćešina—Draženova puhaljka</td><td>Croatia</td><td>15.85710</td><td>44.26920</td><td>26.4</td><td>737</td></tr><tr><td>C21<sup>a</sup></td><td>Oued Chaara</td><td>Morocco</td><td>-4.30768</td><td>33.86006</td><td>7.7</td><td>ND</td></tr><tr><td>C22</td><td>Pınargözü Cave</td><td>Turkey</td><td>31.18270</td><td>37.41480</td><td>16</td><td>720</td></tr><tr><td>C23<sup>a</sup></td><td>Postojnska Jama</td><td>Slovenia</td><td>14.20821</td><td>45.78293</td><td>24.12</td><td>115</td></tr><tr><td>C24<sup>a</sup></td><td>Rhar Bou Ma'za</td><td>Algeria</td><td>-1.25400</td><td>34.70900</td><td>18.4</td><td>35</td></tr><tr><td>C25</td><td>Saint-Marcel d'Ardéche Cave</td><td>France</td><td>4.54100</td><td>44.33210</td><td>51.2</td><td>233</td></tr><tr><td>C26</td><td>Schönberg-Höhlensystem</td><td>Austria</td><td>13.77230</td><td>47.70150</td><td>140.1</td><td>1,061</td></tr><tr><td>C27</td><td>Tilkiler Cave</td><td>Turkey</td><td>31.51067</td><td>36.95344</td><td>6.8</td><td>159</td></tr><tr><td>C28</td><td>Migovec System</td><td>Slovenia</td><td>46.25000</td><td>13.76000</td><td>43</td><td>972</td></tr><tr><td>C29</td><td>Yenesu Cave</td><td>Turkey</td><td>27.56320</td><td>41.37090</td><td>1.6</td><td>ND</td></tr><tr><td>C30</td><td>Željezna Jama</td><td>Montenegro</td><td>19,22890</td><td>42.74340</td><td>17,5</td><td>1,027</td></tr></tbody></table>

Mediterranean focus area (250 km zone), it is 244,569 million m³. At the level of individual countries (the whole area of countries is considered), France has the highest mean annual karst recharge with 83,435 million m³, followed by Turkey with 39,081 million m³ and Italy with 32,418 million m³. Due to the simplified method of calculating mean annual recharge, this is only a solid estimate that should be validated and improved for some regions.

Data on recharge rates or total renewable groundwater from other sources, such as the WaterGAP model (Müller Schmied et al. 2021) or the Food and Agriculture Organisation's AQUASTAT database (FAO 2023), show different values (even when downscaled to karst areas). This is due to the different calculation methodology and database, illustrating the challenge of obtaining accurate and reliable data on groundwater recharge.

## Trends in karst groundwater storage from GRACE

It can be seen from Fig. 6 that karst aquifers in all parts of the Mediterranean are affected by a decrease in groundwater storage, although in some areas there is also a slight increase. In the whole Mediterranean karst region, a mean trend of -1.1 mm/year is detected corresponding to an annual lost volume of 2,544.3 million m³. Within the Mediterranean focus area, the mean trend is -0.3 mm/year, summing up to a total annual loss of 436.2 million m³. Of the 28 countries studied, 23 show an average negative trend in karst groundwater storage change, while only 5 show a positive trend. The strongest negative trends are observed in Malta with -19.3 mm/year, followed by Cyprus with -9.3 mm/year and Austria with -6.0 mm/year. Slightly negative trends are observed in

<footnote><sup>a</sup>Connected to an important karst groundwater-dependent ecosystem (Siegel et al. 2023)</footnote>
<figure><img src="image_10.png" /><figcaption>**Fig. 5** Groundwater recharge, expressed as a 30-year average (1990–2019), derived from the karst-specific simulation model VarKarst-R (Hartmann et al. 2020)</figcaption></figure>

large parts of France, particularly in the Paris region and the southern part of the country; northern Italy, the entire Middle East and parts of the Algerian highlands are particularly affected by strong negative trends. The Dinarides, some areas of eastern Turkey and the coastal fringes of the major karst platforms in North Africa are also affected; however, the latter region also has extensive areas of slightly increasing groundwater storage, as do small areas in southern Spain, southern Italy, Greece and throughout Turkey (Fig. 6). Portugal with 3.0 mm/year, Egypt with 2.1 mm/year, Libya with 1.4 mm/year, Greece with 0.3 mm/year and Lebanon with 0.03 mm/year show on average positive trends. This could either be due to a small increase in natural karst groundwater recharge, or return flows from agricultural irrigation could also be responsible for increasing groundwater storage, which would be plausible in countries with generally very low karst groundwater recharge rates. However, this would not explain small positive trends in nonarable desert regions such as Libya and Egypt. As a result of the different proportions of karst areas in the countries, the annual water loss rates are also very different—for example, Turkey has an average annual loss of 693.8 million m³, followed by France with 645.0 million m³ and Algeria with 558.8 million m³. The positive GRACE trends in Egypt and Libya, with their extensive karst aquifers, add up to a mean annual increase of 856 million and 728.8 million m³ respectively, while in

countries with smaller karst areas, the increase is only a few million m³ (Table 9). The results show the high pressure on the karst groundwater resources in the Mediterranean region and, thus, also on the people and the environment that depend on them. Sustainable use is, therefore, particularly important in arid countries, especially in the context of climate change, which is predicted to be more pronounced in the Mediterranean (Zittis et al. 2019).

## Conclusion and outlook

This study provides an improved detailed assessment of the spatial distribution of carbonate and evaporite rocks as karst aquifers and other hydrogeological formations in the Mediterranean region and the Mediterranean focus area defined by a 250-km zone. It includes karst formations both on the continent and on the main islands and provides the first statistical evaluation of coastal karst. It also presents a selected dataset of a unique combination of terrestrial and submarine karst springs and their range of discharge, of karst caves and their spatial extent. In addition, the results of existing studies on karst groundwater recharge and on the trends in karst groundwater storage were applied to MEDKAM, which provides more accurate and quantified information on potential karst groundwater resources.

The main findings of this study are:

<figure><img src="image_11.png" /></figure>
<table><caption>**Table 8** Country-specific mean value of long-term mean groundwater recharge in carbonate and evaporite rocks (1990–2019)</caption><thead><tr><th rowspan="2">Country</th><th rowspan="2">Total area of carbonate and evaporite rocks<br>1,000 km²</th><th colspan="2">Country-specific mean value of long-term groundwater recharge in carbonate and evaporite rocks (1990–2019)</th></tr><tr><th>mm/year</th><th>Million m³/year</th></tr></thead><tbody><tr><td>Albania</td><td>8</td><td>503.7</td><td>3,861</td></tr><tr><td>Algeria</td><td>308</td><td>50.8</td><td>15,641</td></tr><tr><td>Andorra</td><td>11</td><td>292.0</td><td>3.2</td></tr><tr><td>Austria</td><td>21</td><td>841.6</td><td>17,922</td></tr><tr><td>Bosnia and Herzegovina</td><td>21</td><td>658.0</td><td>14,091</td></tr><tr><td>Bulgaria</td><td>32</td><td>225.2</td><td>7,213</td></tr><tr><td>Croatia</td><td>22</td><td>598.5</td><td>13,219</td></tr><tr><td>Cyprus</td><td>2</td><td>180.5</td><td>408</td></tr><tr><td>Egypt</td><td>411</td><td>1.6</td><td>648</td></tr><tr><td>France</td><td>207</td><td>403.1</td><td>83,435</td></tr><tr><td>Greece</td><td>8</td><td>225.8</td><td>1,731</td></tr><tr><td>Israel</td><td>10</td><td>64.1</td><td>622</td></tr><tr><td>Italy</td><td>69</td><td>472</td><td>32,418</td></tr><tr><td>Jordan</td><td>64</td><td>19.8</td><td>1,272</td></tr><tr><td>Lebanon</td><td>7</td><td>215.8</td><td>1,485</td></tr><tr><td>Libya</td><td>505</td><td>13.5</td><td>6,825</td></tr><tr><td>Malta</td><td>0.03</td><td>65.0</td><td>0.02</td></tr><tr><td>Montenegro</td><td>9</td><td>672.6</td><td>6,130</td></tr><tr><td>Morocco</td><td>154</td><td>102.4</td><td>15,726</td></tr><tr><td>North Macedonia</td><td>4</td><td>359.1</td><td>1,278</td></tr><tr><td>Palestinian Territories</td><td>5</td><td>105.8</td><td>501</td></tr><tr><td>Portugal</td><td>3</td><td>398.6</td><td>1,257</td></tr><tr><td>Serbia</td><td>11</td><td>422.4</td><td>4,540</td></tr><tr><td>Slovenia</td><td>10</td><td>841.8</td><td>8,068</td></tr><tr><td>Spain</td><td>105</td><td>179.4</td><td>18,793</td></tr><tr><td>Switzerland</td><td>13</td><td>934.4</td><td>11,699</td></tr><tr><td>Syria</td><td>89</td><td>78.7</td><td>6,992</td></tr><tr><td>Tunisia</td><td>48</td><td>67</td><td>3,185</td></tr><tr><td>Turkey</td><td>191</td><td>204.4</td><td>39,081</td></tr><tr><td>Mediterranean region</td><td>2,345</td><td>317.2</td><td>740,339</td></tr><tr><td>Mediterranean focus area (250-km zone)</td><td>1,284,500</td><td>190.4</td><td>244,569</td></tr></tbody></table>

* 39.5% of the Mediterranean focus area (250-km zone) is covered by carbonate rocks, corresponding to more than 1.2 million km², while only 1.7% consists of exposed evaporite rocks.

* Looking at the total area of each geographical region, Middle Eastern countries have the highest percentage of carbonate rocks (30.5%), followed by European countries (25.8%) and North African countries (25.3%). For evaporite rocks, the same order is found for the Middle East (2.5%), followed by the European countries (1.1%) and the North African countries (0.5%).

* The largest and most extensive continuous carbonate rocks can be found in North African countries with

301,000 km² in Algeria, 402,000 km² in Egypt and 505,000 km² in Libya.

* The highest proportions of carbonate rocks in relation to the total area of the countries (apart from countries <20,000 km² such as Montenegro or Lebanon) can be found in Jordan with 72.7%, Slovenia with 47.3% and Israel with 44.9%.

* Coastal karst (carbonate and evaporite rocks) is found along ~14,000 km of coastline around the Mediterranean Sea, including islands, which is 33.6% of the total coastline and ~25.9% of the continental coastline.

* At the country level, the longest discontinuous coastal karst (carbonate and evaporite rocks) is in Greece with
<figure><img src="image_12.png" /><figcaption>**Fig. 6** Extent of trends in groundwater storage (GWS) in karst areas over the period 2003–2020 for the Euro-Mediterranean region using the latest data from the Gravity Recovery and Climate Experiment (NASA 2021) satellite mission and recently reanalyzed ERA5 land climate data (Muñoz Sabater 2019) from the European Centre</figcaption></figure>

for Medium-Range Weather Forecasts (modified after Xanke and Liesch 2022). The apparently coarser resolution of the GRACE data in North Africa results from the stereographic projection used, which leads to a stronger distortion in the direction of the equator

~5,016 km, followed by Croatia with ~3,491 km and Italy with ~1,367 km, while Malta and Bosnia and Herzegovina have a 100% proportion of karst on their coast, followed by Croatia with 87.7% and Montenegro with 67.9%.

* The karst groundwater recharge is higher in the northern part of the Mediterranean region than in the southern part and amounts to a total of ~244,600 million m³/year in the Mediterranean focus area (250-km zone). In terms of volume, France is the country with the largest average annual karst groundwater recharge (83,435 million m³).

* Negative trends in karst groundwater storage can be observed in many Mediterranean regions, with only a few areas showing slightly increasing trends such as parts of the Atlas Mountains in Morocco, Algeria, and Tunisia. Also, extensive areas in Libya and Egypt show slightly positive trends. For the Mediterranean focus area (250-

km zone), an average annual loss of karst groundwater resources of 436.2 million m³ is calculated.

Carbonate and evaporite rocks and, thus, karst groundwater resources, are widespread in the Mediterranean region and play an important role in supplying freshwater to people and the environment. In the vast majority of Mediterranean countries, karst aquifers cover more than a quarter of the land area, and in some countries well over 50%, which is remarkable compared to other large regions of the world, as statistical evaluations by Goldscheider et al. (2020) show. This wide distribution also means that karst is found along the coast, especially on some of the many islands, most of which are karstified, emphasizing the vulnerability of these areas to saltwater intrusion and making balanced groundwater use particularly important there (Fleury et al. 2023).

<figure><img src="image_13.png" /></figure>
<table><caption>**Table 9** Country-specific trends in karst groundwater storage in carbonate and evaporite rocks (2003–2020)</caption><thead><tr><th>Country</th><th>Total area of carbonate and evaporite rocks<br>1,000 km²</th><th>Average of mean annual trend in groundwater storage in carbonate and evaporite rocks<br>(2003–2020)<br>mm</th><th>Average of mean annual trend in groundwater storage in carbonate and evaporite rocks<br>(2003–2020)<br>million m³</th></tr></thead><tbody><tr><td>Albania</td><td>8</td><td>-1.62</td><td>-12.4</td></tr><tr><td>Algeria</td><td>308</td><td>-1.82</td><td>-558.8</td></tr><tr><td>Andorra</td><td>11</td><td>-1.76</td><td>-0.019</td></tr><tr><td>Austria</td><td>21</td><td>-6.01</td><td>-128.0</td></tr><tr><td>Bosnia and Herzegovina</td><td>21</td><td>-2.86</td><td>-61.3</td></tr><tr><td>Bulgaria</td><td>32</td><td>-1.82</td><td>-58.2</td></tr><tr><td>Croatia</td><td>22</td><td>-3.06</td><td>-67.6</td></tr><tr><td>Cyprus</td><td>2</td><td>-9.26</td><td>-20.9</td></tr><tr><td>Egypt</td><td>411</td><td>2.08</td><td>856.3</td></tr><tr><td>France</td><td>207</td><td>-3.12</td><td>-645.3</td></tr><tr><td>Greece</td><td>8</td><td>0.33</td><td>2.5</td></tr><tr><td>Israel</td><td>10</td><td>-2.54</td><td>-24.6</td></tr><tr><td>Italy</td><td>69</td><td>-2.18</td><td>-149.4</td></tr><tr><td>Jordan</td><td>64</td><td>-4.85</td><td>-312.6</td></tr><tr><td>Lebanon</td><td>7</td><td>0.03</td><td>0.2</td></tr><tr><td>Libya</td><td>505</td><td>1.44</td><td>728.8</td></tr><tr><td>Malta</td><td>0.03</td><td>-19.26</td><td>-0.01</td></tr><tr><td>Montenegro</td><td>9</td><td>-2.57</td><td>-23.4</td></tr><tr><td>Morocco</td><td>154</td><td>-1.57</td><td>-241.8</td></tr><tr><td>North Macedonia</td><td>4</td><td>-0.84</td><td>-3.0</td></tr><tr><td>Palestinian Territories</td><td>5</td><td>-1.64</td><td>-7.8</td></tr><tr><td>Portugal</td><td>3</td><td>2.96</td><td>9.3</td></tr><tr><td>Serbia</td><td>11</td><td>-3.10</td><td>-33.3</td></tr><tr><td>Slovenia</td><td>10</td><td>-3.99</td><td>-38.3</td></tr><tr><td>Spain</td><td>105</td><td>-1.91</td><td>-200.0</td></tr><tr><td>Switzerland</td><td>13</td><td>-5.51</td><td>-69.0</td></tr><tr><td>Syria</td><td>89</td><td>-4.87</td><td>-432.8</td></tr><tr><td>Tunisia</td><td>48</td><td>-3.01</td><td>-142.9</td></tr><tr><td>Turkey</td><td>191</td><td>-3.63</td><td>-693.8</td></tr><tr><td>Mediterranean region</td><td>2,345</td><td>-3.01</td><td>-2,544.3</td></tr><tr><td>Mediterranean focus area<br>(250 km zone)</td><td>1,285</td><td>-0.34</td><td>-436.2</td></tr></tbody></table>

Karst aquifers in the Mediterranean region are under increasing pressure, as indicated by the negative trends in groundwater storage, which has implications for spring discharge and thus for rivers and groundwater-dependent ecosystems, as well as for regional drinking water supplies. Recent climate projections for the Mediterranean are not optimistic, expecting rising temperatures and reduced water

availability, although with regional differences (Pal et al. 2004; Mariotti et al. 2015; Lionello and Scarascia 2018; Zit-tis et al. 2019). According to Xanke and Liesch (2022), the causes of negative trends are diverse and can be attributed regionally to high, unsustainable groundwater abstraction as well as declining natural groundwater recharge. Positive trends, on the other hand, might regionally be associated
with return flows from irrigation and not necessarily with increased groundwater recharge.

MEDKAM serves as a basis for various research and management approaches at different scales to keep the karst groundwater resources in the Mediterranean region usable for future generations. The spatial and statistical distribution of karst aquifers provides valuable information for scientific and regulatory purposes and can be used in the context of similarly scaled applications and data sets related to climatic, hydrological or geographical issues such as studies by Zhang et al. (2023), who conducted a global analysis of land use changes in karst areas and assessed their impacts on water resources. Additionally, it can be used to upscale or transfer local and regional management options to a supra-regional scale, which may involve, for example, karst-specific protection measures such as vulnerability mapping approaches (Zwahlen 2003), flood prevention and flood regulation approaches (Stevanović 2010; Jourde et al. 2014), or active karst groundwater management through managed aquifer recharge and storage (Xanke 2017). Further development of MEDKAM could include a more detailed scale, or a further subdivision into different types of carbonate rocks and their genesis or the degree of karstification, both in the subsurface and in the epikarst.

## Map availability

The PDF version of the printed MEDKAM and more information on it, as well as a link to download the shapefile, which is free for noncommercial use, can be found here: https://doi.org/10.25928/MEDKAM.1.

**Acknowledgements** This study is a contribution to the project "Karst Aquifer Resources availability and quality in the Mediterranean Area" (KARMA)

**Funding** Open Access funding enabled and organized by Projekt DEAL. Financial support of KIT researchers through the German Federal Ministry of Education and Research (BMBF) and the European Commission through the Partnership for Research and Innovation in the Mediterranean Area (PRIMA) program under Horizon 2020 (grant agreement no. 01DH19022A) is greatly appreciated, as well as funding of our KARMA project partners by their respective national funding agencies (MITECO for Spain, ANR for France, MIUR for Italy, CNRS-L for Lebanon, ANPR for Tunisia). Open Access funding enabled and organized by Project DEAL. The authors thank Salah Hamad for contributing data on karst springs for Libya.

## Declarations

**Conflict of Interest** The authors declare that they have no conflict of interest.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

Audra P, Mocochain L, Camus H, Gilli É, Clauzon G, Bigot JY (2004) The effect of the Messinian Deep Stage on karst development around the Mediterranean Sea: examples from southern France. Geodinam Acta 17(6):389–400. https://doi.org/10.3166/ga.17.389-400

Bakalowicz M (2005) Karst groundwater: a challenge for new resources. Hydrogeol J 13(1):148–160. https://doi.org/10.1007/s10040-004-0402-9

Bakalowicz M (2011) Management of karst groundwater resources. In: van Beynen P (eds) Karst Management. Springer, Dordrecht. https://doi.org/10.1007/978-94-007-1207-2_12

Bakalowicz M (2014) Karst at depth below the sea level around the Mediterranean due to the Messinian crisis of salinity. Hydrogeological consequences and issues. Geol Belg 17(2014):96–101. https://popups.uliege.be/1374-8505/index.php?id=4435

Bakalowicz M (2018) Coastal karst groundwater in the Mediterranean: a resource to be preferably exploited onshore, not from karst submarine springs. Geosciences 8(7):258. https://doi.org/10.3390/geosciences8070258

BGR and UNESCO (2019) International hydrogeological map of Europe 1:1,500,000 (IHME1500). Digital map data v1.2. BRG, Hannover, Germany; /UNESCO, Paris

Bonacci O, Roje-Bonacci T, Vrsalović A (2023) Kuk K (2023) What happened to the karst spring Ombla water temperature? Carbon Evapor 38:68. https://doi.org/10.1007/s13146-023-00892-8

Cavazza W, Wezel FC (2003) The Mediterranean region: a geological primer. Episodes 26(3):160–168

Chen Z, Auler AS, Bakalowicz M, Drew D, Griger F, Hartmann J, Jiang G, Moosdorf N, Richts A, Stevanovic Z, Veni G, Goldscheider N (2017) The World Karst Aquifer Mapping Project: concept, mapping procedure and map of Europe. Hydrogeol J 25:771–785. https://doi.org/10.1007/s10040-016-1519-3

Duggen S, Hoernle K, Van den Bogaard P, Rüpke L, Phipps Morgan J (2003) Deep roots of the Messinian salinity crisis. Nature 422(6932):602–606. https://doi.org/10.1038/nature01553

Duscher K, Günther A, Richts A, Clos P, Philipp U, Struckmeier W (2015) The GIS layers of the “International Hydrogeological Map of Europe 1: 1,500,000” in a vector format. Hydrogeol J 23(8):1867–1875. https://doi.org/10.1007/s10040-015-1296-4

Fan X, Goeppert N, Goldscheider N (2023) Quantifying the historic and future response of karst spring discharge to climate variability and change at a snow-influenced temperate catchment in central Europe. Hydrogeol J 31:2213–2229. https://doi.org/10.1007/s10040-023-02703-9

FAO (2023) AQUASTAT database. AQUASTAT website. https://www.fao.org/aquastat/statistics/query/index.html. Accessed Mar 2023

<figure><img src="image_14.png" /></figure>
Fleury P, Bakalowicz M, de Marsily G (2007) Submarine springs and coastal karst aquifers: a review. J Hydrol 339(1–2):79–92. https://doi.org/10.1016/j.jhydrol.2007.03.009

Fleury P, Pistre S, Bakalowicz M (2023) Coastal karst aquifers and submarine springs: what future for their water resources? C R Géosci 355(S1):1–14. https://doi.org/10.5802/crgeos.168

Ford D, Williams PD (2007) Karst hydrogeology and geomorphology. Wiley, Chichester, UK

Ford DC, Williams PW (1989) Karst geomorphology and hydrology, vol 601. Hyman, London

Goldscheider N (2019) A holistic approach to groundwater protection and ecosystem services in karst terrains. Carbon Evapor 34(4):1241–1249. https://doi.org/10.1007/s13146-019-00492-5

Goldscheider N, Chen Z, Broda S, Auler AS, Bakalowicz M, Drew D, Hartmann J, Jiang G, Moosdorf N, Stevanovic Z, Veni G (2020) Global distribution of carbonate rocks and karst water resources. Hydrogeol J 28(5):1661–1677. https://doi.org/10.1007/s10040-020-02139-5

Goldscheider N, Drew D (eds) (2007) Methods in karst hydrogeology. International Contributions to Hydrogeology, 26. CRC, Boca Raton

Gunn J (2021) Karst groundwater in UNESCO protected areas: a global overview. Hydrogeol J 29:297–314. https://doi.org/10.1007/s10040-020-02206-x

Günther A, Duscher K (2019) Extended vector data of the International Hydrogeological Map of Europe 1:1,500,000 (Version IHME1500 v1.2). Technical note, BGR, Hanover, Germany

Hamad SM, El Hasia A (2022) An overview of Karst groundwater springs in Al Jabal Al Akhdar region (North East Libya). In: Groundwater for Sustainable Livelihoods and Equitable Growth. Taylor and Francis, London, pp 179–194. https://doi.org/10.1209/9781003024101

Hartmann J, Moosdorf N (2012) The new global lithological map database GLiM: A representation of rock properties at the Earth surface. Geochem Geophys Geosyst 13(12). https://doi.org/10.1029/2012GC004370

Hartmann J, Jansen N, Duerr HH, Kempe S, Koehler P (2009) Global CO2-consumption by chemical weathering: what is the contribution of highly active weathering regions? Glob Planet Chang 69(4):185–194. https://doi.org/10.1016/j.gloplacha.2009.07.007

Hartmann A, Goldscheider N, Wagener T, Lange J, Weiler M (2014) Karst water resources in a changing world: review of hydrological modeling approaches. Rev Geophys 52:218–242. https://doi.org/10.1002/2012RG000443

Hartmann A, Liu Y, Olarinoye T, Berthelin R, Marx V (2020) Integrating field work, large sample hydrology and modeling to inform (inter) national governance of karst water resources. https://doi.org/10.31223/osf.io/8cde3

Hersbach H, Bell B, Berrisford P, Hirahara S, Horányi A, Muñoz-Sabater J, et al (2020) The ERA5 global reanalysis. Q J Royal Meteorol Soc 146(730):1999–2049. https://doi.org/10.1002/qj.3803

Jolivet L, Augier R, Robin C, Suc JP, Rouchy JM (2006) Lithospheric-scale geodynamic context of the Messinian salinity crisis. Sed Geol 188–189:9–33. https://doi.org/10.1016/j.sedgeo.2006.02.004

Jourde H, Lafare A, Mazzilli N et al (2014) Flash flood mitigation as a positive consequence of anthropogenic forcing on the groundwater resource in a karst catchment. Environ Earth Sci 71:573–583. https://doi.org/10.1007/s12665-013-2678-3

Koeppen W (1936) Das geographische System der Klimate [The geographical system of climates]. In: Koeppen W, Geiger G (eds) Handbuch der Klimatologie [Manual of climatology]. Borntraeger, Berlin, 44 pp

LaMoreaux PE, Tanner JT (2001) Springs and bottled waters of the world: ancient history, source, occurrence, quality and use. Springer, Heidelberg, Germany

Lionello P, Scarascia L (2018) The relation between climate change in the Mediterranean region and global warming. Reg Environ Change 18(5):1481–1493. https://doi.org/10.1007/s10113-018-1290-1

MacDonald AM, Bonsor HC, Dochartaigh BÉÓ, Taylor RG (2012) Quantitative maps of groundwater resources in Africa. Environ Res Lett 7(2):024009. https://doi.org/10.1088/1748-9326/7/2/024009

Mackenzie FT (1978) Carbonate mineralogy and geochemistry. In: Middleton GV, Church MJ, Coniglio M, Hardie LA, Longstaffe FJ (eds) Encyclopedia of sediments and sedimentary rocks. Encyclopedia of earth sciences series. Springer, Dordrecht, The Netherlands. https://doi.org/10.1007/978-1-4020-3609-5_35

Mariotti A, Pan Y, Zeng N, Alessandri A (2015) Long-term climate change in the Mediterranean region in the midst of decadal variability. Clim Dyn 44(5–6):1437–1456. https://doi.org/10.1007/s00382-015-2487-3

Martens B, Miralles DG, Lievens H, Van Der Schalie R, De Jeu RA, Fernández-Prieto D, Beck HE, Dorigo WA, Verhoest NE (2017) GLEAM v3: Satellite-based land evaporation and root-zone soil moisture. Geosci Model Dev 10(5):1903–1925. https://doi.org/10.5194/gmd-10-1903-2017

Müller Schmied H, Cáceres D, Eisner S, Flörke M, Herbert C, Niemann C, Peiris TA, Popat E, Portmann FT, Reinecke R, Schumacher M, Shadkam S, Telteu CE, Trautmann T, Döll P (2021) The global water resources and use model WaterGAP v2. 2D: model description and evaluation. Geosci Model Dev 14(2):1037–1079. https://doi.org/10.5194/gmd-14-1037-2021

Muñoz Sabater J (2019) ERA5-Land hourly data from 1981 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). https://doi.org/10.24381/cds.e2161bac. Accessed Mar 2021

NASA (2021) Jet Propulsion Laboratory. GRACE & GRACE-FO - data months/days. https://grace.jpl.nasa.gov/data/grace-months/. Accessed 12 Nov 2021

NIMA (National Imagery and Mapping Agency) (2001) Performance Specification Vector Smart Map (VMap) Level 0: MILPRF-89039, Amendment 2

Pal JS, Giorgi F, Bi X (2004) Consistency of recent European summer precipitation trends and extremes with future regional climate projections. Geophys Res Lett 31(13). https://doi.org/10.1029/2004GL019836

Pawlewicz MJ, Steinshouer DW, Gautier DL (2002) Map showing geology, oil and gas fields, and geologic provinces of Europe including Turkey. US Geol Surv Open-File Rep 97-470-I, 14 pp. https://doi.org/10.3133/ofr97470I

Peel MC, Finlayson BL, McMahon TA (2007) Updated world map of the Köppen-Geiger climate classification. Hydrol Earth Syst Sci 11(5):1633–1644. https://doi.org/10.5194/hess-11-1633-2007

Persims FM, Ahlbrandt TS, Tuttle ML, Charpentier RR, Brownfield ME, Takahashi KI (1997) Maps showing geology, oil and gas fields and geological provinces of Africa. US Geol Surv Open-File Rep 97-470-A. https://doi.org/10.3133/ofr97470A

Pollastro RM, Karshbaum AS, Viger RJ (1999) Maps showing geology, oil and gas fields and geologic provinces of the Arabian Peninsula. US Geol Surv Open-File Rep 97-470-B, 14 pp. https://doi.org/10.3133/ofr97470B.

Pollastro RM (1998) Bedrock geology of the Arabian Peninsula and selected adjacent areas (geo2bg). US Geol Surv Data Release

Ravbar N (2013) Variability of groundwater flow and transport processes in karst under different hydrologic conditions. Acta Carsol 42(2–3). https://doi.org/10.3986/ac.v42i2.644
Richts A, Vrba J (2016) Groundwater resources and hydroclimatic extremes: mapping global groundwater vulnerability to floods and droughts. *Environ Earth Sci* 75(10):1–15

Richts A, Struckmeier WF, Zaepke M (2011) WHYMAP and the groundwater resources map of the world 1: 25,000,000. In: Sus-
taining groundwater resources (pp 159–173). Springer, Dordrecht,
The Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10

Rodell M, Houser PR, Jambor U, Gottschalck J, Mitchell K, Meng CJ, Arsenault K, Cosgrove A, Radakovich J, Bosilovich M, Entin JK, Walker JP, Lohmann D, Toll D (2004) The global land data assimila-
tion system. Bull Am Meteor Soc 85(3):381–394

Roveri M, Flecker R, Krijgsman W, Lofi J, Lugli S, Manzi V, Sierro JF, Bertini A, Camerlenghi A, De Lange G, Govers R, Hilgen JF, Hübscher C, Meijer TP, Stoica M (2014) The Messinian salinity crisis: past and future of a great challenge for marine sciences. Mar Geol 352:25–58. https://doi.org/10.1016/j.margeo.2014.02.002

Save H (2020) CSR GRACE and GRACE-FO RL06 Mascon Solutions v02. https://doi.org/10.15781/cgq9-nh24

Save H, Bettadpur S, Tapley BD (2016) High resolution CSR GRACE RL05 Mascons. J Geophys Res Solid Earth 121:7547–7569. https://doi.org/10.1002/2016JB013007

Schuler P, Margane A (2013) Water balance for the groundwater contribution zone of Jeita Spring using WEAP: including water resources management options & scenarios. Technical Coopera-
tion Project Protection of Jeita Spring, German-Lebanese Techni-
cal Cooperation Project Protection of Jeita Spring, BGR, Orleans,
France

Siegel L, Goldscheider N, Petitta M, Xanke J, Andreo B, Bakalow-
icz M, Barberá JA, Bouhlila R, Burg A, Doummar J, Ezzine I,
Fernández-Ortega J, Ghanmi M, Jourde H, Marín AI, Mhimdí
A, Pipan T, Ravbar N, Stevanović AM, Stevanović Z (2023) Dis-
tribution, threats and protection of selected karst groundwater-
dependent ecosystems in the Mediterranean region. Hydrogeol J.
https://doi.org/10.1007/s10040-023-02711-9

Stevanović Z (2010) Intake of Bolje Sestre karst spring for the regional water supply of the Montenegro coastal area. In: Kresic N, Stevanović Z (eds) Groundwater hydrology of springs. Elsevier, Amsterdam, pp 457–477. https://doi.org/10.1016/B978-1-85617-502-9.00023-2

Stevanović Z (ed) (2015) Karst aquifers-characterization and engineer-
ing. Springer. https://doi.org/10.1007/978-3-319-12850-4

Stevanović Z (2019) Karst waters in potable water supply: a global scale overview. Environ Earth Sci 78:662. https://doi.org/10.1007/s12665-019-8670-9

Stevanović Z, Milanović P (2023) South-eastern Dinaric karst: con-
trasts in water treasury. Environ Earth Sci 82:215. https://doi.org/
10.1007/s12665-023-10904-4

USGS (2021) Geologic provinces of the world. https://certmapper.cr.usgs.
gov/data/apps/world-maps/. Accessed 25 May 2021

Williams PW, Ford DC (2006) Global distribution of carbonate rocks. Zeitschr Geomorphol 147(Suppl):1–2

Xanke J (2017) Managed aquifer recharge into a karst groundwater system at the Wala reservoir, Jordan. PhD Dissertation, Karlsruher Institut für Technologie (KIT). https://doi.org/10.5445/IR/1000068919

Xanke J, Liesch T (2022) Quantification and possible causes of declin-
ing groundwater resources in the Euro-Mediterranean region from
2003 to 2020. Hydrogeol J 30:379–400. https://doi.org/10.1007/
s10040-021-02448-3

Xanke J, Goldscheider N, Bakalowicz M, Barberá JA, Broda S, Chen Z, Ghanmi M, Günther A, Hartmann A, Jourde H, Liesch T, Mudarra M, Petitta M, Ravbar N, Stevanović Z (2022) Mediterranean Karst Aquifer Map (MEDKAM), 1:5,000,000. https://doi.org/10.25928/MEDKAM.1

Zhang J, Liesch CZ, Goldscheider N (2023) Global analysis of land-use changes in karst areas and the implications for water resources. Hydrogeol J 31(5):1197–1208. https://doi.org/10.1007/s10040-023-02650-5

Zittis G, Hadjinicolaou P, Klangidou M, Proestos Y, Lelieveld J (2019) A multi-model, multi-scenario, and multi-domain analy-
sis of regional climate projections for the Mediterranean. Reg
Environ Change 19(8):2621–2635. https://doi.org/10.1007/
s10113-019-01565-w

Zwahlen F (ed) (2003) Vulnerability and risk mapping for the protec-
tion of carbonate (karst) aquifers. Office for Official Publications
of the European Communities, Luxembourg

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.