Introduction

Groundwater is an essential lifeline for various purposes, especially in arid and semi-arid regions1. Kuwait faces persistent water resource shortages and water quality issues in an arid region. Kuwait faces one of the highest per capita consumption rates in the world. The country grapples with acute water scarcity, leading to significant annual expenses on direct water supply and distribution2. Kuwait’s renewable water resources per person have experienced a gradual decline, decreasing from 17.3 m3/year in 1977 to 4.8 m3/year in 20182. As a result of this scarcity, the country heavily depends on desalination plants to manage its freshwater requirements.

The availability of usable groundwater in Kuwait is limited, with the water mostly being brackish to salty. Unfortunately, even these available scarce fresh groundwater lenses are at risk due to mismanagement (excessive pumping), leading to a decline in water quality (salinity deterioration) through excessive pumping3. An example can be found in Al-Rawdatain and Umm Al-Aish areas, where substantial fresh groundwater reserves provided approximately 20 mm3/year during the early 1960s. However, recent mismanagement has significantly reduced this potential to merely 3 mm3/year. Moreover, improper placement of well screens between transition zones of different water bodies or poor design in brackish-saline zones3 can exacerbate salinity issues, negatively affecting Kuwait’s water resources. Even brackish groundwater with total dissolved solids (TDS) ranging between 1000 and 10,000 mg/l, which is still considered strategically crucial for Kuwait, has not been spared from the detrimental effects of over-pumping. Salinity levels have surged from 4000 to 20,000 mg/l following excessive pumping activities.

Addressing the depletion of groundwater resources in Kuwait requires effective management of pumping practices, comprehensive plans for groundwater exploration, and the strategic siting of groundwater boreholes based on well-planned studies. Recently, remote sensing and geoinformatics have demonstrated their effectiveness in identifying potential groundwater zones4,5,6,7,8,9,10,11. Additionally, diverse multi-criteria decision-making methods such as regression models12, frequency ratio13, influencing factor14, and Analytic Hierarchy Process (AHP)6,15 have introduced good results regarding groundwater potentiality mapping.

In Table 1; Fig. 1, a literature review is presented. Literature is organized to show which thematic layers different authors used to identify groundwater potential zones in other parts of the world. Based on Fig. 1, it is clear that among 15 studies, all of them used slope in their research. Most of them used land use and land cover (14), lineament density (14), drainage density (13), geology/lithology and rainfall (12), and geomorphology (9). Other thematic layers, such as the digital elevation model, normalized differential vegetation index, and topographic wetness index, are used in a few research studies. Indices such as NDWI, LST study, and DTW are not frequently used, which suggests that our research has a holistic approach toward our aim.

Fig. 1
figure 1

Overview of the number of thematic layers used in literature. GL geology and lithology; G geomorphology; LULC land use & land cover; R rainfall; LD lineament density; DD drainage density; S slope; NDVI normalized difference vegetation index; NDWI normalized difference water index; LST land surface temperature; TWI topographic wetness index; DTW depth to water table.

To tackle this pressing issue for Kuwait, the current research endeavors to develop a solution by modeling groundwater potential zones in Kuwait, integrating multiple parameters using Fuzzy-AHP and GIS techniques. The key novelty of this research is that it marks the first comprehensive study to characterize Kuwait’s groundwater resources by integrating a wide range of datasets, including geological, hydrological, and climatic information. By employing advanced multicriteria analysis alongside cutting-edge web-based tools such as Google Earth Engine (GEE), this work sets a new precedent in groundwater assessment for the region. The main contribution of the current research is the holistic application of all relevant factors after conducting a comprehensive literature review to build a groundwater potential map for Kuwait as a whole. Additionally, the current study integrates GEE analysis of climatic conditions and various indices with geological, morphological, and hydrogeological factors to provide suitable recommendations for the entire groundwater system in Kuwait. These approaches have shown promising results in offering viable recommendations for water-related challenges in arid and semi-arid surroundings.

Table 1 Literature review of the theme used to identify groundwater potential zones.

Description of the study area

Site description

Kuwait is located in the northwestern part of the Arabian Gulf, situated between latitudes 28⁰ 30′ and 30⁰ 05′ to the north and longitudes 33⁰ 46′ and 30⁰ 48′ to the east Fig. 2. It shares its southern and southwestern borders with Saudi Arabia, while its northwestern border is with Iraq. Kuwait’s eastern coastline borders the Arabian Gulf. Figure 1 illustrates that the country includes several islands, where the largest is Bubiyan Island, followed by Warba Island, and the nearest to the coast, Failaka Island. The total land area covers 17,818 km2.

Fig. 2
figure 2

The study area’s location map (Kuwait) shows its main features and the drainage network over the state.

Generally, Kuwait’s topography is characterized by a mostly flat terrain with a somewhat uneven desert consisting of gently rolling plains covered in gravel. The land gradually inclines from sea level along the eastern coast of the Arabian Gulf towards the west and southwest. Based on Fig. 3A,B. The highest point in the southwestern corner reaches 284 m above sea level. Small hills are scattered across Kuwait, including those along the Jal Al-Zour ridge, which rises to 145 m and overlooks the northern coast of Kuwait Bay. This ridge stretches from the northeast to the southwest, near Jahra. Additionally, there are hills like Al-Laiyah and Keraa Al-Marw, which stand at 138 m in height, and Al-Ahmadi Hill in the southern region of Kuwait, reaching an elevation of 125 m. These hills are more extensive and elongated than the earlier ones, and they are covered with pebbles of various sizes.

Kuwait also features numerous valleys, lowlands locally known as Al-Khubarat, and dunes. Among the notable valleys in Kuwait are the Al-Batin Valley, which runs along the country’s western borders, and the Al-Sheqaq Valleys, situated in the northwestern part of the country. Al-Khubarat are distributed in various areas, with the most significant being Al-Rawdatain and Umm Al-Aish, both located in the northern region26. The central geological units covering Kuwait are introduced in Fig. 3C. Figure 3D illustrates the extensive array of geomorphological features covering Kuwait’s landscape. These features consist of desert plains and dunes distributed across Kuwait’s south and southeastern areas and shaped by the effects of wind, resulting in impressive Aeolian desert formations. The primary geomorphic characteristics include sandy areas (mobile sand sheets), stony surfaces (containing calcrete and gypcrete deposits), minor ridges, and Karst depressions. Additionally, the region encompasses coastal plains and tidal flat features, with the most significant geomorphic elements being recurved spits, sandy beach deposits, silt and clay-based evaporite deposits, and sabkhas. Various igneous and metamorphic rocks also characterize gravel ridges (comprising pebbles and gravel). Inter-ridge surfaces and playa deposits are found across the Fan Delta surface, primarily influenced by the ancient river system along Wadi Al-Batin26,27. Figure 3E shows the general level of the water table within Kuwait, and the distribution of the total dissolved solids is presented in Fig. 3F.

Fig. 3
figure 3

Research area: (A) Digital elevation model29; (B) Aspect; (C) General geology; (D) Geomorphology30; (E) Depth to water31; (F) Total Dissolved solids31.

Climatic conditions

The climate of the State of Kuwait is characterized by long, hot, dry summers and short, warm winters with occasional rain. During the summer months, dust-laden winds are common, and humidity in the air increases32. The highest temperature recorded throughout history was 51 °C in July 1978, and the lowest reached − 4 °C in January 1964. According to information from the General Directorate of Civil Aviation, average temperatures usually range between 47 °C in summer and 6 °C in winter. These significant temperature differences vary with changes in annual precipitation rates, ranging from 22 mm to 352 mm per year33. Figure 4A presents average precipitation. The “mean daily maximum” (solid red line) shows the maximum temperature of an average day every month in Kuwait. Likewise, the “mean daily minimum” (solid blue line) shows the average minimum temperature. Hot days and cold nights (dashed red and blue lines) show the average hottest day and coldest night of each month of the last 30 years. Temperatures, and be prepared for hotter and colder days. Wind speeds are not showed per default but can be enabled at the bottom of the graph. Figure 4B presents wind rose in Kuwait.

Fig. 4
figure 4

(A) Average monthly precipitation and temperature for Kuwait; (B) Rose wind34.

Geology and hydrogeology

Figure 3C illustrates the exposed rock formations in the study area, originating from sedimentary deposits from the Palaeocene to the recent period. These formations can be categorized into Al-Hasa and Kuwait35. The Al-Hasa group comprises limestone, dolomite, anhydrite, and clay. It is divided into three formations: Umm er Radhuma from the Palaeocene to the Middle Eocene, Rus in the Lower Eocene, and Damman in the Middle Eocene. The Damman Formation unconformably overlays the Middle Eocene dolomite-limestone of the Dammam Formation. It varies in thickness from 390 m in the northeast to just 8 m in the southeast of Kuwait36.

On the other hand, the Kuwait Group is composed of fluvial sediments such as sand, gravel, calcareous sand, sandstone, clays, gypsum, limestone, and marl. It consists of three formation units, which were formed in Ghar in the Miocene, Fars in the Pliocene, and Dibdibba in the Pleistocene. These formations originated from torrential floods that transported rock formations from the Arabian Peninsula and deposited them in northern Kuwait. The lower member of the Dibdibba Formation is extensively exposed in Kuwait’s western and north-central parts, consisting of weathered medium to very coarse-grained calcretic sandstones with grit and pebbles. The upper member of the Dibdibba Formation can be found beneath the undulating plain in the north and northeast Kuwait, extending south-westward southwestward along the characteristic curved ridges of the northwestern area. This upper member consists of a sequence of ungraded, often cross-bedded sands and gravels with occasional intercalations of sandy mud37.

Additionally, the study area includes sedimentary formations from the Quaternary era, specifically the Pleistocene and Holocene. These formations encompass land deposits of gravel and sand and coastal deposits of calcareous sandstone, sand, and mud. The area is characterized by multiple lineaments with various orientations, including N-NE/S-SW, NE/SW, and N-NW/S-SE38. Experience rapid reductions in floodwater flow due to quick infiltration into the dry and porous substrate. In the past, river networks supplied water to the desert, eventually forming groundwater. However, Kuwait’s groundwater is considered fossil water due to challenges in replenishing and maintaining it under current climatic conditions. Kuwait has two main groundwater layers. The upper layer, known as the Kuwait Group aquifer, consists of greenish gravelly sands, estimated to be about 7,000–10,000 years old, coinciding with the northward incursion of the summer monsoon in northern Arabia. Intense summer monsoons likely contributed to forming active stream networks in northeastern Arabia and Kuwait, carving the landscape, and forming the Kuwait Group aquifer. The lower layer, known as the Dammam formation, comprises chalky limestone and dolomite. The Kuwait Group Aquifer is unconfined, allowing water to infiltrate from the surface, while the Dammam Formation is confined, with impermeable layers above and below it39,40. As illustrated in Fig. 3E, groundwater flows from higher elevations in the southwest, approximately 90 m above sea level, toward the Gulf Coast in the northeast, consistent with Kuwait’s topography. This understanding is supported by hydrogeological studies that confirm the vulnerability of the Kuwait Group to contamination and the salinity issues in the Dammam Formation resulting from over-extraction. Groundwater quality varies from moderately saline in the southwest to highly saline in the northeastern region, as shown in Fig. 3F.

Freshwater lenses are found in the north, such as in Al-Roudtain and Umm Al-Aish. Several geological and rock conditions contributed to their formation, including heavy rainfall exceeding the infiltration capacity into surface depressions, the creation of local surface pools, and the rapid infiltration of accumulated rainwater into the groundwater. These freshwater resources contain less than 1000 mg/L of total dissolved solids and have been used for drinking and bottled water production, with their extraction for drinking purposes ceasing due to the development of desalination plants31.

Data used

Elevation

The region’s altitude is a basic aspect of its topography and plays a crucial role in evaluating groundwater possibilities. Utilizing elevation data from ALOS PALSAR, we examined the topography of the area, visually depicted in Fig. 3A, B.

Slope

The slope is pivotal in influencing various factors, including surface runoff, groundwater recharge, etc. Steeper slopes hinder groundwater infiltration with regards to the rapid increase in surface runoff during the rainfall and heightened soil erosion rates. Conversely, areas with gentler slopes promote more significant groundwater recharge, allowing rainwater to remain longer and gradually seep into the ground6. Consequently, the hill demonstrates an inverse relationship between infiltration rates and surface water infiltration while positively correlating with surface runoff41. This aspect can prolong the time water infiltrates the soil in areas with gentle flooding, ultimately raising the groundwater Tables19,21,42,43,44,45. To determine the surface slope, spatial analysis tools within the ArcGIS environment were employed, utilizing high-precision ALOS PALSAR digital elevation models (DEM) with a 12.5 meters spatial resolution as shown in the Fig. 5A46,47.

Drainage density

Drainage density is of a importance in assessing groundwater movement and recharge as it directly influences the speed of water flow and its ability to penetrate the aquifer. The relationship between drainage density and surface runoff, which subsequently affects water infiltration into the ground, is inversely related to permeability5,6,48. It is vital to identify potential groundwater areas since the drainage density influences the groundwater recharge rate in a specific region. High drainage density areas generally have a lower potential for groundwater recharge due to increased surface runoff rates. Conversely, regions with lower drainage density often exhibit higher groundwater potential as they offer more significant opportunities for groundwater replenishment. ALOS PALSAR DEM data was used to create a drainage network map shown in Fig. 5C, followed by the Line Density spatial analysis tool in ArcGIS to produce a drainage density map18,20,21,44,49,50,51.

Lineament density

Lineaments play a crucial role in replenishing underground groundwater. They are responsible for enhancing the porosity and permeability of the essential groundwater indicators when rainwater infiltrates the soil. Lineaments also assist in identifying potential groundwater locations. The presence and flow of groundwater depend on the infiltration facilitated by these contours. A higher density of lineaments indicates a more significant potential for groundwater recharge, while a lower density suggests a lower likelihood of finding a potential groundwater source. According to Magowe and Carr21,52,53 there is a positive correlation between increased lineament density and potential groundwater areas. To implement this theory, a numerical score was assigned to lineament density. The identification of areas with the greatest groundwater potential was based on their proximity to regions characterized by high lineament density, which were subsequently assigned the highest grades, as demonstrated in Table A (1). Utilizing the “kernel density” method within a Geographic Information System (GIS) framework, a lineament density map was generated to visualize these findings as shown in the Fig. 5D.

Rainfall

Rainfall is indispensable for the hydrological cycle as it serves as the main source of freshwater and directly impacts the replenishment and augmentation of groundwater that percolates into the earth’s subsurface layers through fissures and soil or via water flowing through the drainage network, depending on topographical conditions. The region relies on seasonal rainfall received from the northeastern and southwestern monsoons. Researchers have highlighted the crucial role of rainfall density, intensity, and duration in influencing the percolation process. High-intensity, short-duration rainfall leads to reduced percolation and increased surface runoff, whereas low-intensity, long-duration rainfall contributes to a higher percolation rate than surface runoff. Consequently, the quantity of rainfall in the area matters, and the duration and density of rainfall play a crucial role in groundwater recharge. Various researchers have emphasized a strong positive relationship between rainfall and groundwater likelihood, indicating that heavy rainfall suggests a more significant potential for groundwater recharge. Conversely, reduced rainfall has a lower chance of reaching the groundwater18,21,43,50. Rainfall data from the IMERG dataset for the period 2000–2022, obtained from the Precipitation Measurement Missions https://www.chrsdata.eng.uci.edu, (PMM) at the University of California, Irvine, with a spatial accuracy of 4 km was utilized in this study. Rainfall data was imported into ArcGIS as raster data and employed to create a rainfall map as shown in the Fig. 5E.

Topographic wetness index

The Topographic Wetness Index (TWI) plays major role within the hydrological system, providing as a reflection of potential groundwater seepage affected by the area’s topographical features. It has proven crucial in pinpointing the extent and location of saturated sources adding to surface runoff generation. Therefore, areas with the highest TWI values are labeled as holding the most substantial groundwater potential. Many researchers have utilized TWI as a key factor to delineate groundwater potential areas45,46,51. As highlighted in a study by Beven and Kirkby52, TWI connects ascending areas, acting as a measure of water flow towards a specific point (i.e., the local slope), thereby indicating subsurface lateral transmissivity. The TWI is determined using the following equation.

$$\:TWI=\text{ln}\left(\frac{As\:}{tan\beta\:}\right)$$

In this context, “As” represents the specific catchment area, and “tan β” corresponds to the slope angle at the given location.

Groundwater depth

Groundwater depth presents the upper boundary of the saturated zone; therefore, the depth of the water table is a key element of the hydrological cycle in the region, alongside topography and rainfall patterns. Changes in groundwater pressure result in fluctuations in the water table level.

Land use and land cover

Land use and land cover (LULC) are essential factors affecting the groundwater recharge rate, surface runoff, and evaporation12,36. This usage provides essential information on various topics such as leakage and recharge, soil moisture, and the distribution of surface and groundwater and their occurrence (da Costa et al., 2019). The LULC map data was obtained from the Global Land Cover database37 with a 10-meter resolution as shown in the Fig. 6A. Subsequently, adjustments were made by cross-referencing field data and aerial photographs through ArcGIS 10.8 to ensure accuracy.

Geomorphology

The region’s geomorphology is a fundamental element for dividing groundwater recharge areas, as it provides evidence of various shapes and processes, such as temperature fluctuations, groundwater movement, and geochemical interactions within the Earth’s subsurface, and more24,56,57,58. Seepage capabilities will affect the potential for groundwater recharge in different geomorphological units59. Barchan Sand Dunes, Active Sand Sheets, Desert Floor Deposits, Siliciclastic Granule Lag, Fall Dunes, Deflated Rugges Sand Sheets, Smooth Sand Sheets, Gravel Lag, Barchanoid Ridges, Playa Deposits, Urban areas, Coastal Plain Deposits, Desert Floor Deposits Calcretic Debris are among the vital geomorphological features in the study area as shown in Fig. 3D.

Geology

The potential existence of groundwater major relies to geology features of the area. Geological features heavily influence groundwater’s presence and distribution across an area. Rainwater struggles to seep through without suitable rock formations, and groundwater finds storing underground tricky. Water movement and storage capacity within groundwater layers vary based on rock characteristics like porosity and permeability6,44. Different geological indicators results in a significant role in the movement and accumulation of groundwater, necessitating careful consideration in groundwater research. For instance, rocks like sandstone, which possess high porosity, efficiently retain water, whereas low-porosity rocks such as clay and oil-bearing formations have limited water storage capacity. Moreover, fractures and voids within rocks serve as vital conduits for guiding groundwater movement and facilitating its flow between rock layers60.

Normalized differential vegetation index (NDVI)

NDVI is a reliable gauge of groundwater levels, particularly in these specific regions and ecosystems. Robust, lush vegetation with a high greenness index boosts groundwater recharge. This is achieved through the plant’s canopy and root systems, intercepting rainfall, and slowing surface water runoff. Consequently, more water can seep into the soil, replenishing groundwater. Thus, a high greenness index can significantly increase groundwater recharge in vegetated areas. To assess Groundwater Potential (GWP) and identify land cover in the study area, a pixel-based NDVI within the range of -1.0 to + 1.0 was computed in this research as shown in the Fig. 6B, utilizing the equation from Taloor et al.6,16,61,62 applied to Sentinel 2 based Google Earth Engine data.

NDVI = (NIR - RED) / (NIR + RED), where RED is B4, 664.5 nm, NIR is B8, 835.1 nm.

Normalized differential water index (NDWI)

The NDWI is considered an effective tool for exploring potential groundwater areas, as the presence of water in the soil or vegetation serves as an indirect indicator of groundwater availability. When plants contain a high percentage of water, it increases the vegetation cover density, allowing more water to penetrate the soil, thus replenishing groundwater. Higher NDWI values are associated with increased groundwater nourishment, and areas with consistently high NDWI values are considered reliable groundwater sources. The role of the NDWI is more prominent in areas with high soil moisture content, indicating the presence of water in the soil. This indirectly serves as evidence of groundwater potential in the region. The NDWI in this study was calculated using the equation proposed by Liu et al.63 and applied to data from the Google Earth engine based on Landsat 8. The formula used is (NIR - SWIR) / (NIR + SWIR), where NIR represents the near-infrared range and SWIR represents the shortwave infrared range. Generally, higher NDWI values indicate higher water content in the observed area as shown in the Fig. 6C.

Land surface temperature

Land surface temperature (LST) is a valuable indicator for assessing groundwater potential in a specific area, considering various environmental factors. Lower temperatures, often linked to vegetation and specific ground covers, may suggest conducive conditions for groundwater recharge. In urban settings, impervious surfaces like asphalt and concrete contribute to elevated ground surface temperatures, reducing groundwater recharge potential. Conversely, vegetation aids in cooling the land surface through shade and evaporation, enhancing groundwater possibilities in densely vegetated regions. In our study, Landsat 8 data was utilized to derive LST values, which were applied to Google Earth Engine data based on thermal bands from satellite images (bands 10 and 11). The process involved converting Digital Numbers (DN) to radiance using sensor calibration parameters, subsequently calculating brightness temperature per Planck’s law (Eq. 1), followed by conversion to Land surface temperature using the Stefan-Boltzmann law (Eq. 2).

$$\:{T}_{b}=\frac{{k}_{2}}{{ln}\left(\frac{{K}_{1}}{L}+1\right)}$$

where: Tb is the brightness temperature, k2and K1 are sensor-specific calibration constants, L is the radiance.

$$\:LST={\left(\frac{{T}_{b}}{\epsilon\:}\right)}^{\frac{1}{4}}$$

Where: LST is the Land Surface Temperature, ε is the emissivity.

Table 2 The input data, spatial accuracy, and sources used in this study.
Fig. 5
figure 5

Research area: (A) Slope; (B) Flow accumulation; (C) Drainage density; (D) Lineament density; (E) Annual Rainfall65; (F) Topographic wetness index.

Fig. 6
figure 6

Research area: (A) Land use and land cover; (B) Normalized differential vegetation index; (C) Normalized differential water index; (D) Land surface temperature.

Methodology

Mapping the groundwater potential zone with conventional techniques is quite difficult. The lack of thematic layers in the traditional approach—such as data on soil, groundwater level, and geomorphology—causes problems. The absence of historical maps or data from pertinent departments and organizations makes the process of gathering data even more difficult. This makes it more difficult to validate the groundwater potential zone mapping’s correctness and dependability. These difficulties bring to light the shortcomings and inadequacies of the conventional approaches used to carry out these mapping exercises. The creation of thematic layers using satellite data has become more commonplace recently for a variety of uses, such as mapping groundwater potential zones. Satellite imagery is a useful tool for obtaining information about vegetation indices, land cover, and other pertinent factors. Thematic layers obtained from satellite data can be used as a basis for groundwater potential analysis and assessment. This method enhances the final map’s accuracy and dependability. Satellite data were used to produce several thematic layers, such as: drainage density, slope, elevation, groundwater level, land use and land cover (LULC), geomorphology, soil, and drainage density.

Figure 7 illustrates the conceptual flowchart outlining the research design and methodology employed in this study. The flowchart delineates the automation process facilitated by Remote Sensing (RS) and Geographic Information System (GIS) software, accommodating spatial data variation between maps and remotely sensed data as shown in the Table 2 such as Alos PALSAR-DEM and satellite imagery. OpenStreetMap (OSM) serves as a valuable resource, offering user-acquired data freely accessible for download. Integration of OSM can enhance base maps while amalgamating diverse maps generated through remote sensing techniques. Furthermore, the proposed map delineates the preparation of data in a suitable format. According to the characteristics and multiplicity of spatial datasets used, the research methodology handle it while preprocessing the data using the data management resample tool in the Arc-GIS program; that tools allow controlling the cell size of the data that comes out of it, where all outputs were standardized to become in the cell size of 10 m, another preprocessing of data are conversion to raster tools (as polygon to raster or interpolation tools). The detailed steps of this methodology are elaborated in the subsequent sub-sections.

Fig. 7
figure 7

Conceptual flowchart of this research is connected.

Working with topographic maps and satellite image analysis

The topographic map underwent geometric correction/georeferencing utilizing ESRI ArcGIS software, employing the Orthorectification process to mitigate spatial distortion and align the X and Y dimensions to a specific geographic reference system, the Universal Transverse Mercator (UTM) coordinate system—Datum: World Geodetic System (WGS 1984), Zone 39 North. Subsequently, the map was digitized onscreen, and its spatial data were saved in a vector formats as a feature classes geodatabase (Geo-Db). The Geo-Db was constructed using the ArcGIS program (Arc Catalog module) to generate various vector data types of the study area features: points, lines, and polygons. Identification and attribute data were then inputted. Following this, the spatial sub-setting process was applied to select and extract the entire raster data from the topographic map sheet using the spatial analyst toolbox/extraction toolset/extract by mask tool. Additionally, a more recent mapping of roads and other features within the studied area was conducted using ESRI Satellite Images.

Working with digital elevation model

Studying terrain characteristics is crucial for the park’s development process. A Digital Elevation Model (DEM) with a spatial resolution of  12.5 m was prepared for the study using data from the ALOS PALSAR Terrain. These data underwent preprocessing in ArcGIS software and were re-projected to align with the geographic coordinate system (UTM – Datum: WGS 1984 zone 39 N) for consistency with maps. The fill enhancement function was utilized to rectify errors in the digital elevation model (DEM). Additionally, a hydrology model was employed to generate the drainage network and basin watershed and boundary for the studied area.

Google earth engine

Scalable, cloud-based geospatial retrieval and processing are made possible by Google Earth Engine (GEE). Additionally, GEE offers free cloud-based computing capability for geospatial data analysis as well as access to the overwhelming majority of publicly accessible, multi-temporal Remotely Sensed data of satellite imagery and geospatial datasets (such as Landsat, Sentinel, and MODIS, …etc.)68. The geospatial community has embraced GEE recently and facilitated various local, regional, and global environmental applications. Maps depicting the world’s changing forests are among the significant environmental applications for global urban change58,59; GEE offers JavaScript and Python APIs for communicating with the Earth Engine servers to do computations. indices.

One sensor’s worth of images is provided by Earth Engine and is compiled into a “collection” called image collection. Using a straightforward graphical user interface (GUI) and the Function.filter (), users can also define a point or polygon anywhere on Earth and filter image collection data from each satellite for user-defined time frames, area extent or location, and the percentage of an image pixel that is contaminated by cloud interference. They can also choose from a variety of predefined or custom image band combinations that align with the objective and research question. A script was created to filter date, check the cloud-free availability of images, preprocessing and calculation of used indices: NDVI, NDWI, and LST.

Weighting methods of criteria

This section offers an introduction to the Fuzzy Analytical Hierarchy Process (FAHP) and entropy-based aggregation methods employed to prioritize influential criteria from various experts’ viewpoints.

This method can assist policymakers in setting priorities and producing appropriate outcomes. It is a helpful tool for resolving difficult decision-making challenges. By substituting a series of pairwise comparisons and findings for composite judgments, the AHP supports both subjective and objective aspects of a choice. Moreover, decision-making bias can be eliminated by using the AHP as a useful tool to evaluate the correctness of decisionmakers’ assessments. It is predicated on the pairwise comparison between the selected criteria that is produced by the experts’ assessment. To determine the relative importance of the most significant themes and their characteristics, views were solicited. There may be some disagreement because the comparisons are based on subjective or individual viewpoint. One of the most important benefits of the AHP is that it can be used to quantify the degree of consistency among paired comparisons of various criteria. This consistency ratio was found to be one way to verify that the perceptions are consistent. The AHP covers and integrates spatial data as input and produces decisions as output. By creating a pairwise comparison matrix with Saaty’s scale, the qualitative data of many themes and attributes are converted into quantitative data69,70,71,72. Weights were adopted that were consistent with the natural and environmental conditions of the study area, and criteria that were not included in similar areas and were applied in other areas of a different nature were modified to match the conditions of the study area according to consultation with experts and specialists.

Trapezoidal interval type II fuzzy Analytical Hierarchy process

Gupta and Lee59 provided the trapezoidal Interval Type-2 Fuzzy Analytical Hierarchy Process (IT2 FAHP) approach. This method is used to provide the weights of the criteria in the analysis. Experts’ responses, based on AHP are documented in decision matrices. These decision matrices are transformed into trapezoidal fuzzy numbers73. Following this, the aggregated comparison matrix is constructed by calculating the arithmetic mean of the respondents’ responses according to Eq. (1):

$$\:downright{Z}_{q}={\left({{\stackrel{\sim}{\stackrel{\sim}{x}}}_{ij}}^{q}\right)}_{m\times\:n}=\left[\begin{array}{ccc}{{\stackrel{\sim}{\stackrel{\sim}{x}}}_{11}}^{q}&\:\cdots\:&\:{{\stackrel{\sim}{\stackrel{\sim}{x}}}_{1n}}^{q}\\\:&\:\ddots\:&\:\\\:{{\stackrel{\sim}{\stackrel{\sim}{x}}}_{m1}}^{q}&\:\cdots\:&\:{{\stackrel{\sim}{\stackrel{\sim}{x}}}_{mn}}^{q}\end{array}\right],\:1\le\:i\le\:m,\:1\le\:j\le\:n,\:1\le\:q\le\:k$$
(1)

Where; \(\:{\stackrel{\sim}{\stackrel{\sim}{x}}}_{ij}\) represents the aggregated matrix of the upper and lower IT2 fuzzy sets, which should be calculated: \(\:{\stackrel{\sim}{\stackrel{\sim}{\:x}}}_{ij}=\left(\frac{{{\stackrel{\sim}{\stackrel{\sim}{x}}}_{ij}}^{1}+{{\stackrel{\sim}{\stackrel{\sim}{x}}}_{ij}}^{2}+\dots\:+{{\stackrel{\sim}{\stackrel{\sim}{x}}}_{ij}}^{k}}{k}\right)\). Besides, \(\:k\:\)refers to the number of respondents. Equations (23) can be selected to obtain the upper (\(\:{P}^{U}\)) and lower (\(\:{P}^{L}\)) fuzzy matrices:

$$\:{P}^{U}=\left[\begin{array}{ccc}p({{\stackrel{\sim}{A}}_{1}}^{U}\ge\:{{\stackrel{\sim}{A}}_{1}}^{U})&\:\cdots\:&\:p({{\stackrel{\sim}{A}}_{1}}^{U}\ge\:{{\stackrel{\sim}{A}}_{n}}^{U})\\\:&\:\ddots\:&\:\\\:p({{\stackrel{\sim}{A}}_{n}}^{U}\ge\:{{\stackrel{\sim}{A}}_{1}}^{U})&\:\cdots\:&\:p({{\stackrel{\sim}{A}}_{n}}^{U}\ge\:{{\stackrel{\sim}{A}}_{n}}^{U})\end{array}\right],\:1\le\:i\le\:n$$
(2)
$$\:{P}^{L}=\left[\begin{array}{ccc}p({{\stackrel{\sim}{A}}_{1}}^{L}\ge\:{{\stackrel{\sim}{A}}_{1}}^{L})&\:\cdots\:&\:p({{\stackrel{\sim}{A}}_{1}}^{L}\ge\:{{\stackrel{\sim}{A}}_{n}}^{L})\\\:&\:\ddots\:&\:\\\:p({{\stackrel{\sim}{A}}_{n}}^{L}\ge\:{{\stackrel{\sim}{A}}_{1}}^{L})&\:\cdots\:&\:p({{\stackrel{\sim}{A}}_{n}}^{L}\ge\:{{\stackrel{\sim}{A}}_{n}}^{L})\end{array}\right],\:1\le\:i\le\:n\:$$
(3)

The ranking value of the upper membership function\(\:Rank\left({{\stackrel{\sim}{A}}_{i}}^{U}\right)\), lower membership function\(\:Rank\left({{\stackrel{\sim}{A}}_{i}}^{L}\right)\), and fuzzy set\(\:Rank\:\left({\stackrel{\sim}{\stackrel{\sim}{A}}}_{i}\right)\)are computed using Eqs. (46):

$$\:Rank\left({{\stackrel{\sim}{A}}_{i}}^{U}\right)=\frac{1}{n(n-1)}\left(\sum\:_{k=1}^{n}p\left({{\stackrel{\sim}{A}}_{i}}^{U}\ge\:{{\stackrel{\sim}{A}}_{k}}^{U}\right)+\frac{n}{2}-1\right),\:1\le\:i\le\:n$$
(4)
$$\:Rank\left({{\stackrel{\sim}{A}}_{i}}^{L}\right)=\frac{1}{n(n-1)}\left(\sum\:_{k=1}^{n}p\left({{\stackrel{\sim}{A}}_{i}}^{L}\ge\:{{\stackrel{\sim}{A}}_{k}}^{L}\right)+\frac{n}{2}-1\right),\:1\le\:i\le\:n$$
(5)
$$\:Rank\:\left({\stackrel{\sim}{\stackrel{\sim}{A}}}_{i}\right)=\frac{Rank\left({{\stackrel{\sim}{A}}_{i}}^{U}\right)+Rank\left({{\stackrel{\sim}{A}}_{i}}^{L}\right)}{2},\:1\le\:i\le\:n$$
(6)

Where; \(\:\sum\:_{i=1}^{n}Rank\:\left({\stackrel{\sim}{\stackrel{\sim}{A}}}_{i}\right)=1\)

Finally, Eqs. (7, 8) are used to normalize ranking values of the upper membership function (\(\:{{d}_{ij}}^{U}\)), lower membership function (\(\:{{d}_{ij}}^{L}\)). To calculate the relative weights, the fuzzy set’s average normalized ranking value (\(\:{d}_{ij}\)) is computed using Eq. (9):

$$\:{{d}_{ij}}^{U}=\frac{Rank\left({{\stackrel{\sim}{A}}_{i}}^{U}\right)}{\sum\:_{j=1}^{J}Rank\left({{\stackrel{\sim}{A}}_{i}}^{U}\right)},\:1\le\:i\le\:n,\:1\le\:j\le\:J$$
(7)
$$\:{{d}_{ij}}^{L}=\frac{Rank\left({{\stackrel{\sim}{A}}_{i}}^{L}\right)}{\sum\:_{j=1}^{J}Rank\left({{\stackrel{\sim}{A}}_{i}}^{L}\right)},\:1\le\:i\le\:n,\:1\le\:j\le\:J$$
(8)
$$\:{d}_{ij}=\frac{{{d}_{ij}}^{U}+{{d}_{ij}}^{L}}{2}$$
(9)

Where; \(\:{d}_{ij}\) indicates the preference rating of the criteria and it shall be noted that\(\:\sum\:{d}_{ij}=1\).

MCDA suitability model according to entropy-based aggregation approach

Gupta and Lee59 provided the trapezoidal Interval Type-2 Fuzzy Analytical Hierarchy Process (IT2 FAHP) approach. We used this approach to calculate the weights of the criteria in the analysis. The weighted average of each criterion (\(\:{{w}_{i}}^{avg}\)) can be defined using Eq. (10):

$$\:{{w}_{i}}^{avg}=\frac{\sum\:_{j=1}^{m}{{w}_{i}}^{j}\times\:{S}_{j}}{\sum\:_{j=1}^{m}{S}_{j}}$$
(10)

Where; \(\:{S}_{j}\) refers to the number of respondents in the jth group.

Measure the entropy value (\(\:{E}_{j}\)) for each group of respondents using Eq. (11):

$$\:{E}_{j}=\sum\:_{i=1}^{n}\left(\frac{-1}{\text{l}\text{n}\left(2\right)}\left({{w}_{i}}^{j}\times\:\text{ln}\left({{w}_{i}}^{j}\right)+{{w}_{i}}^{avg}\times\:\text{l}\text{n}\left({{w}_{i}}^{avg}\right)\right)\right)$$
(11)

Where n refers to the number of criteria.

The diversification value (\(\:{div}_{j}\)) for each corespondent and normalize these values (\(\:{\stackrel{-}{w}}_{j}\)) to determine the significance level of each group using Eqs. (12, 13):

$$\:{div}_{j}=n-{E}_{j}$$
(12)
$$\:{\stackrel{-}{w}}_{j}=\frac{{div}_{j}}{\sum\:_{j=1}^{m}{div}_{j}}$$
(13)

Calculate the compromise weight (\(\:{cw}_{i}\)) of criteria and rank the estimated weights in descending order to determine the most influential criteria using Eq. (14):

$$\:{cw}_{i}=\sum\:_{j=1}^{m}{\stackrel{-}{w}}_{j}\times\:{{w}_{i}}^{j}$$
(14)

Finally, Table A (1) presents assigned weights and ratings based on each factor, while Table A(2) presents mutual independence and ratings among criteria with a consistency check of 10%.

Results

A comprehensive set of thirteen thematic maps was constructed to evaluate the groundwater potential in Kuwait. These maps encompassed data related to various factors such as geology, slope, rainfall, elevation, Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Land Surface Temperature (LST), drainage density, Land Use and Land Cover (LULC), lineament density, topographic wetness index (TWI), geomorphology, and the water table. Creating these datasets involved combining remote sensing techniques and traditional data sources, utilizing software tools like ERDAS and ArcGIS.

Among these thematic maps, topographic elevation, slope, drainage density, and TWI maps were derived from ALOS PALSAR RTC data. The remaining maps were developed using data from various sources, including Landsat-8, SENTINEL-2 satellites, and conventional resources like geology, rainfall, water table, NDVI, NDWI, LST, and LULC maps. All converted layers underwent reclassification to ensure data consistency, aligning them with the study site’s specific requirements. They were then resampled to establish a uniform cell size of 10 m. Reclassification simplified the interpretation of raster data by merging ranges of cell values into a unified value based on a consistent ranking system, with object (using value) in the ramp reclassify tool (spatial analysis tool-box in Arc-GIS software) is used to specify how to reclassify values of the input raster. This streamlined the process of comparing and ranking sites based on their suitability. Additionally, this study analyzed the relationship between the variables employed for predicting groundwater potential during the weight calculation and the actual groundwater potential.

Our DEM elevation analysis revealed that the highest point reached 290.94 m above sea level (ASL), and the lowest point descended 22 m below sea level. According to a study by Melese et al.52 and Hasanuzzaman et al.25, coastal regions generally exhibit lower elevations, significantly impacting groundwater availability. Additionally, it is observed that lower elevations have a more significant influence on groundwater presence compared to higher elevations, leading to variations in climatic conditions. These variations can affect rainfall patterns, soil quality, vegetation types, land use, and coverage. To delve deeper into this analysis, as depicted in SM 1 (A), the terrain is categorized into nine classes based on elevation, assigning the lowest rank to nearly flat areas and the highest rank to those with relatively high runoff elevations. In summary, elevation emerges as a critical determinant of groundwater potential, holding a secondary position with a normalized weight of 0.066 about its impact on groundwater potential areas, as illustrated in Table A (1) based on the chosen factors in this study.

Slope results (Fig. 5A)) revealed that the study area predominantly exhibited slopes between 0° and 63°. Consequently, the surface slope was categorized into five classes on the slope map. The class with the lowest value was assigned a higher ranking due to its nearly flat terrain, while the class with the highest value was assigned a lower ranking owing to its relatively high runoff levels. The slope is significant, with a normalized weight of 0.046, based on its relative influence on identifying potential groundwater areas among the selected parameters in this study as presented in Table A (1). Utilizing standard and customized weight values, the study area’s slope was categorized into nine groups, with slopes less than 5 degrees indicating good potential. In contrast, those exceeding 40 degrees signify significantly low probability, as SM1 (B) illustrates.

Besides the slope, a comprehensive drainage analysis highlighted the runoff tributaries and the drainage density. Our results highlighted that drainage density values range from 0 to 17.6 m2/km2, as depicted in Fig. 5C. Among the selected parameters in this study, drainage density holds the third position with a normalized weight of 0.136 due to its relative importance in determining potential groundwater areas. Consequently, areas with low drainage density received higher scores, while regions with dense drainage received lower scores, as shown in Table A (1). The study area underwent reclassification into nine distinct categories for analytical purposes, unveiling a notable groundwater recharge potential despite its low drainage density. Regions characterized by low density exhibit a moderate capacity for recharge, whereas those with medium density demonstrate a lower recharging capability. Conversely, areas marked by high drainage density exhibit a notably low potential for groundwater recharge, as depicted in SM1 (C), Table A (1).

A structural analysis was performed through lineaments analysis, shown in Fig. 5D. The latter revealed that lineament density varied from 0 to 4.27 m/km2. Lineament density emerged as the fourth most significant parameter in delineating groundwater potential zones in this study as indicated in Table A (1), with a normalized weight value of 0.037. The line density map revealed that the central part of the study area is highly suitable for groundwater due to the density of lineaments. In contrast, the western and the eastern regions have a lower potential for groundwater recharge as the lineament density decreases from the center towards the east and west as depicted in Fig. 5D. The study area underwent reclassification into nine clearly defined categories, revealing that regions characterized by low line density display diminished recharge capacity. Conversely, areas featuring high line density demonstrate a notable potential for groundwater recharge, as illustrated in SM 1(D) and Table A (1).

The rainfall results are one of the most critical factors in the current research. Our findings revealed that the average annual precipitation varies between 39.8 mm and 107.1 mm. A map with satisfactory outcomes was created using the inverse distance weighting (IDW) method. Consequently, regions in the north, receiving the most rainfall, were rated the highest, while ratings diminished towards the southern areas. Figure 5E illustrates the rainfall map categorization, adjusted to a 10 m2 spatial resolution, and segmented into nine subclasses through a natural fraction classification technique. Regions with scant precipitation were depicted as having minimal recharge potential, whereas those experiencing abundant rainfall were deemed to possess high recharge potential, as shown in SM 1(E). Given the substantial influence of rainfall on potential areas, it was attributed to a high weight of 0.042. as determined among the selected parameters in this study (Table A1).

Besides the rainfall results, TWI was calculated to enhance our results. Figure 5F visually depicts the range of Topographic Wetness Index (TWI) values in our investigation, which spanned from − 13 to 27 within the study area. These values were subsequently categorized into nine groups for further analysis. Higher weights were assigned to high topographic wetness index values and, conversely, lower significance to lower values, as illustrated in SM1 (F). Given the significant impact of the topographic wetness index on potential groundwater areas, a substantial weight of 0.206 was allocated to it. as determined among the selected parameters in this study Table A (1).

Rainfall is also a significant contributor to fluctuations in the groundwater level, as it serves as a source of recharge to the groundwater layer, replenishing it from rainfall. Thus, a depth-to-water table analysis was performed based on previous studies. It is well-known that changes in the water table within the region’s topography primarily result from variations in rainfall patterns, groundwater recharge, and discharge34,35. Groundwater in Kuwait is predominantly characterized by two distinct aquifers: the Kuwait Group aquifer and the Dammam aquifer. As illustrated in Fig. 3E, the water table depth in Kuwait ranges from 2 to 30 m, with the map displaying a gradual decline in the water table depth towards the outlet in the study area. Higher weights were assigned to lower water table depths and vice versa. Values were reclassified into nine categories, as shown in SM2(A). Given the significant impact of the water table depth on potential groundwater areas, a substantial weight of 0.125 was assigned to it, as determined among the selected parameters in this study Table A(1).

LULC analysis is considered one of the most common investigations in groundwater potentiality mapping in recent years. Our results revealed that The LULC map for the studied area consists of seven main categories, which are Water, Trees, Flooded vegetation, Crops, Built Area, Bare ground, and Rangeland, as shown in Fig. 6A land use patterns in the study area. Water bodies and agricultural land are excellent sources for groundwater recharge because they favor rainwater filtration and an increased rate of groundwater replenishment. In contrast, arid and built-up areas are considered less important in this context. The data indicates their contribution to increased surface runoff, reduced leakage rates, and decreased groundwater replenishment. For this reason, classifications have been assigned according to the available water quantity for replenishment based on land surface characteristics. Agricultural lands were given the highest rating of 0.00%, paved built-up areas were given the lowest rating of 0.25%, and water bodies were assigned a classification coefficient of 0.11%. Higher weights were allocated to these categories in the following order: Vegetation, Trees, Flooded vegetation, waterbody, Rangeland, Bare land, Built Up Area, and Clouds, respectively, as shown in SM2 (B) in terms of their impact on groundwater replenishment Table A (1). LULC is ranked with a natural weight of 0.021 based on its relative impact on potential groundwater areas, as determined among the selected parameters in this study Table A (1).

Besides the previous factors, geomorphological and geological features were analyzed to highlight recharge capabilities within the study area, as shown in SM2(C). The latter denotes grades assigned based on the recharge capabilities of different geomorphological areas. Subsequently, subjective weights were assigned based on the characteristics of each layer. Coastal Plain Deposits and Sand and Gravel deposits were rated highly due to their excellent permeability, receiving a grade of nine. Conversely, Playa Deposits and Urban areas, known for their high water-holding capacity, were given a lower grade of one. Highlands and plateaus influenced this decision and typically facilitate surface runoff. Geomorphology is positioned with a natural weighting of 0.052, reflecting its relative impact on potential groundwater zones within the specified parameters of this study Table A (1). Geologically, our results and findings depend mainly on the rock units’ distribution in Fig. 3C. Utilizing the Analytic Hierarchy Process (AHP) framework, the resultant geological map underwent reclassification into nine primary categories based on their hydrogeological characteristics, as depicted in SM 2 (D). Geology is positioned with a natural weighting of 0.026, reflecting its relative impact on potential groundwater zones within the specified parameters of this study Table A(1).

One of the most critical factors affecting GWPZ is NDVI. Our results of NDVI through GEE highlighted that the resulting NDVI map values (see Fig. 6B) span from − 0.22 to 0.43, with values exceeding 0.4, signifying the presence of dense vegetation. Values lower than 0.15 indicate an absence of vegetation, such as barren land and sand, while values below zero indicate urban areas and water bodies, as noted by Rulinda et al.39. The NDVI images were subsequently classified into nine groups, revealing differences in the distribution of plants between parts of the study area, and they are spread in the Abdali area in northern Kuwait, Al-Wafra in southern Kuwait, and the Sulaibiya agricultural area.as shown in SM3 (A). NDVI is positioned with a natural weighting of 0.121, reflecting its relative impact on potential groundwater zones within the specified parameters of this study Table A(1).

Additionally, and using GEE, NDWI was calculated to strengthen our results. The resulting NDWI map values range from − 0.58 to 0.50, as shown in Fig. 6C. Subsequently, the NDWI images were classified into nine groups, revealing differences in surface water distribution among different parts of the study area. These differences are apparent in the eastern part of the country (See SM3 (B)). NDWI is positioned with a natural weighting of 0.046, reflecting its relative impact on potential groundwater zones within the specified parameters of this study Table A(1).

Besides all the previous factors, LST was studied as an indicator for GWPZ, where the lower the LST, the higher the probability of closer groundwater zones to the surface. The resulting LST map values range from − 0.58 to 0.50. as shown in Fig. 6D. Subsequently, the LST images were classified into six groups, revealing differences in surface water distribution among different parts of the study area. These differences are apparent in the eastern part of the country (See SM3 (C)). LST is positioned with a natural weighting of 0.076, reflecting its relative impact on potential groundwater zones within the specified parameters of this study Table A(1).

Groundwater potential zones (GWPZs)

The results of the previous factors were spatially analyzed to check the GWPZs within the study area. Our spatial overlay analysis and weighted analysis findings are introduced in Fig. 8; Table 3. Based on these two supported results, groundwater potential zones in the north and east parts of Kuwait are very highly suitable. However, very high suitability is less than 1.5%, and suitable zones with high and moderate degrees cover 12% and 37%, respectively. Based on input data, these are logical results since rainfall, NDVI, NDWI, and LST impact results in this area. Lower parts of the country have higher chances of groundwater potential. Shallow The low potential covers nearly 50% of the country’s area, covering more than 7500 km2.

Fig. 8
figure 8

Mapping of Groundwater potential zone levels in Kuwait.

Table 3 Suitability degree.

Verification for GWPZs

Field verification and statistical validation (kappa coefficient) were integrated to check the reliability of our methodology and confirm our research findings (Locations of GWPZs). The Kappa coefficient, commonly used to estimate interrater reliability, can be employed in the context of test–retest. In the test–retest, the Kappa coefficient indicates the extent of agreement between frequencies of two sets of data collected on two different occasions. Using this formula:

$$\:K=\frac{{P}_{agree}-{P}_{chance}}{1-{P}_{chance}}$$
(15)

P agree shows the proportion of trials in which judges agree, and P chance shows the proportion of trials in which agreement would be expected due to chance74. If the results are from 0.81 to 1.00, it shows almost perfect agreement, while 0.61 to 0.80 shows substantial validation quality. Based on the mentioned Table 4, the Kappa coefficient of the suitability scale is presented. Results indicate very high agreement results.

Table 4 Kappa coefficient.

Field verification is presented in Fig. 9. According to results shown in Figs. 8 and 9, and Table 4 (kappa coefficient), it is clear how the utilized methodology delivered successful findings through the excellent matching between the locations of water wells and highly potential GWZs. Based on Table 4, a high suitability scale of 0.77 indicates 49 wells in that area. If we compare the results with Fig. 9A,B, most of the wells are found in very high groundwater potential zones (beige color), and significant parts are in high groundwater potential zones (yellow). Table 4, a very high suitability scale, shows the density ratio of 0.076 wells per km2, suggesting that the highest value of existing wells is on the high suitability scale. This could mean that the kappa coefficient of 0.76 was not correctly calculated. However, the suitable area size for a suitable area size for a high suitability scale, suppose we check the appropriate area size for a high suitability scale. In that case, the results show that the area size suits a high suitability scale. We get logical results, and it can be concluded that our research is defined and conducted correctly.

Fig. 9
figure 9

Field validation points (water wells) zones match the high potential zones.

The minimum highest areas (central west area) are the highest. The minimum distribution of wells is located in the highest areas (central west area), consistent with the research area’s geology, geomorphology, and digital elevation model. Also, it indicated low shallow, shallow groundwater potential zones (Figs. 8 and 9). These results are in harmony with Tables 3 and 4. Based on the utilized approach and validation methods, our results show a considered show a considered show a considered shows how a considered show a considerable degree of consistency. Furthermore, our methodology could be applied not only in Kuwait but also in similar research areas in the world.

Discussion

This section will focus on addressing the comparison of our findings with literature that studied groundwater potential sin kuwait focusing mainly on their applied methodology and findings as shown in Table 5. The table highlights the utilized methods and findings for several hydrological and hydrogeological studies dealing with the characterization of water resources in Kuwait over the last 3 decades. This table highlights the novelty of our research on Kuwait, as we did not find any recent studies characterizing groundwater resources in Kuwait by integrating several datasets (geological, hydrological, climatic, etc.) using advanced multicriteria analysis and implementing web-based tools, e.g., GEE.

Table 5 Overview of Groundwater studies and findings in Kuwait.

Extensive efforts were made to compile detailed information on groundwater yield and the primary agricultural fields reliant on groundwater resources. Our findings confirm significant groundwater usage in the regions identified as high-potential zones. Figure 10a–d highlighted the distribution of agricultural fields and their associated irrigation wells. These wells are predominantly situated within the high-potential groundwater zones identified in our analysis, providing a robust cross-validation of our results. This integration further substantiates the accuracy and relevance of the identified groundwater potential zones for agricultural applications.

Fig. 10
figure 10

Distribution of agricultural fields and associated irrigation wells.

Our research offers valuable insights into managing groundwater resources in Kuwait’s semi-arid regions, crucial for sustainable agriculture. By identifying high-potential groundwater zones using Google Earth Engine and MCDM, we provide a framework that can guide water resource allocation and agricultural planning. The findings suggest targeted groundwater extraction in areas such as the north and east, which would optimize water use for crops. Policy recommendations include promoting efficient irrigation practices, integrating groundwater with other water sources like treated wastewater, and enforcing groundwater monitoring to prevent overuse. These strategies align with Kuwait’s Vision 2035 for enhancing food security and agricultural sustainability.

Future research should consider the impact of climate change on groundwater availability and explore additional datasets to refine predictions. Technological advancements, such as AI-driven models, could further improve groundwater management. Limitations include data accuracy and the regional specificity of the study, which may require adjustments for broader applications. Despite these, our study provides a practical, data-driven approach for optimizing water resources, with immediate implications for policymakers in Kuwait.

Conclusions and recommendations

In a Semi-Arid Region of KUWAIT, groundwater potential zone was detected using Multi-Criteria Decision Making. Analyzing several previous types of research and all sorts of studies evaluating groundwater potential zones, thirteen criteria were studies evaluating groundwater potential zones; thirteen criteria were selected to incorporate multicriteria decisions. The highest assigned weight is on rainfall, 0.162; topographical wetness index, 0.13; and NDVI, 0.122. Based on different input weights, different results emerge. However, the main idea is to find the perfect combination of parameters and weights. The main idea is to see an ideal combination of parameters and weights to evaluate groundwater potential. Comparing our results to other studies is almost impossible due to the different parameters used, the various parameters used, and the area’s specific characteristics. However, some similarities can be drawn to other authors on similar (larger scale) areas.

The topographic wetness index with rainfall has the most effect on the research and detects the most impact on the research study of groundwater potential in semi-arid regions. Our research shows that combining more factors can increase precision while ensuring the stability of the detection process of groundwater potential zones. Another indication of quality research is. It shows the kappa values of suitability scale ranks and the kappa values of suitability scale ranks, which show substantial kappa coefficient values of 0.76 and 0.77, respectively. Our research concluded that a high degree of suitability for groundwater potential in Kuwait’s semi-arid region is 242 km2.

Locations, areas, and water wells of very high and high potential zones are located in the northeast part of the research area and also in the south part of the research area in Kuwait. This is presented in our research and ground-validated and ground-validated. The latest research findings offer valuable insights for both public and commercial entities involved in groundwater management, land use planning, and environmental conservation within the researched area. Furthermore, the methodology employed in this study can be adapted to analyze similar watersheds sharing comparable hydrogeologic, topographic, and climatic attributes. Various strategies, such as enhancing groundwater replenishment, augmenting water supply, and implementing rainwater harvesting initiatives in urban settings, can improve water resource management. Additionally, using advanced geospatial techniques is imperative to manage water resources better, curb overexploitation, and explore untapped areas effectively. A more comprehensive multidisciplinary approach is advocated to refine groundwater potential zone mapping techniques, thereby fostering sustainable water resource management. Future research endeavors need to scrutinize the effects of climate change and urbanization on groundwater recharge and potential. Incorporating advanced remote sensing data like LIDAR, machine learning algorithms, real-time monitoring, and an integrated approach is recommended for future groundwater potential zone mapping, ensuring the system’s sustainability. Emphasis is placed on conducting further research encompassing socioeconomic dimensions, long-term monitoring, field validation, and extending the scope to different geographical regions. Such future undertakings are pivotal in promoting sustainable water resource management, meeting the escalating water demands across various industries, and advancing global understanding of groundwater potential.