Identification of Groundwater Potential Recharge Zones in Flinders Ranges, South Australia Using Remote Sensing, GIS, and MIF Techniques

In Australia, water resource management is a major environmental, biological, and socio-economic issue, and will be an essential component of future development. The Hawker Area of the central Flinders Ranges, South Australia suffers from a lack of reliable data to help with water resource management and decision making. The present study aimed to delineate and assess groundwater recharge potential (GWRP) zones using an integration between the remote sensing (RS), geographic information system (GIS), and multi-influencing factors (MIF) approaches in the Hawker Area of the Flinders Ranges, South Australia. Many thematic layers such as lithology, drainage density, slope, and lineament density were established in a GIS environment for the purpose of identifying groundwater recharge potential zones. A knowledge base ranking from 1 to 5 was assigned to each individual thematic layer and its categories, depending on each layer’s importance to groundwater recharge potential zones. All of the thematic layers were integrated to create a combined groundwater potential map of the study area using weighting analysis in ArcGIS software. The groundwater potential zones were categorized into three classes, good, moderate, and low. The resulting zones were verified using available water data and showed a relative consistency with the interpretations. The findings of this study indicated that the most effective groundwater potential recharge zones are located where the lineament density is high, the drainage density is low, and the slope is gentle. The least effective areas for groundwater recharge are underlain by shale and siltstone. The results indicated that there were interrelationships between the groundwater recharge potential factors and the general hydrology characteristics scores of the catchment. MIF analysis using GIS mapping techniques proved to be a very useful tool in the evaluation of hydrogeological systems and could enable decision makers to evaluate, better manage, and protect a hydrogeological system using a single platform.


Introduction
In Australia, groundwater is vital for domestic, agricultural, mining, and industrial purposes [1]. It makes up approximately 17% of the available water resources and over 30% of the total water consumption and is accessible to 60% of the continent's total area [2]. In last decades, most of the continent has faced water scarcity due to the overuse of groundwater that led to an unacceptable decline in water levels causing serious issues such as increased groundwater salinity, high pumping costs, seawater intrusion, and loss of access to groundwater by users and ecosystems [1]. Traditional methods such as geophysical and

Geology and Hydrogeology
The study area is a part of the Adelaide Fold Belt, which consists of Paleoproterozoi to Mesoproterozoic basement rocks extending beyond a thick sedimentary sequence tha is Neoproterozoic to Cambrian [29][30][31]. The crystalline basement and the overlying sedi mentary rocks were deformed in an orogenic setting in the late Cambrian-early Ordovi cian Delamerian Orogeny to form the typical strike ridges of the Flinders Ranges ( Figur  3). Structurally, the study area is fragmented by an anticline-syncline NW and NNW fold system resulting in a series of parallel ridges separated by relatively wide valleys [30]. Th ridges are dominated by high resistance Cambrian and Neoproterozoic quartzites, and

Geology and Hydrogeology
The study area is a part of the Adelaide Fold Belt, which consists of Paleoproterozoic to Mesoproterozoic basement rocks extending beyond a thick sedimentary sequence that is Neoproterozoic to Cambrian [29][30][31]. The crystalline basement and the overlying sedimentary rocks were deformed in an orogenic setting in the late Cambrian-early Ordovician Delamerian Orogeny to form the typical strike ridges of the Flinders Ranges ( Figure 3). Structurally, the study area is fragmented by an anticline-syncline NW and NNW fold system resulting in a series of parallel ridges separated by relatively wide valleys [30]. The ridges are dominated by high resistance Cambrian and Neoproterozoic quartzites, and limestones and calcareous meta-siltstones (Wilpena, Umberatana, and Hawker groups), while alluvial deposits are outcropped in valleys and low topographic areas.
Most of the information regarding the hydrogeological setting was sourced from mineral or hydrocarbon exploration reports [34]. Groundwater of variable quantity and quality occurs throughout the Quaternary sediments and the underlying fractured metasediments [28]. Throughout the study area, the ridges are considered to be groundwater divides recharging the more distal parts in the low topographic areas. The thickness of the Quaternary sediments is variable ranging from 100 m in the northwest, to around 20 m in the southeast [34].
Most of the information regarding the hydrogeological setting was sourc mineral or hydrocarbon exploration reports [34]. Groundwater of variable quan quality occurs throughout the Quaternary sediments and the underlying f metasediments [28]. Throughout the study area, the ridges are considered to be  [31][32][33].

Processing and Extraction of Thematic Layers
In this study, groundwater recharge potential zones were delineated and evaluated by preparing and integrating RS and GIS data. Five influencing factors including lithology, lineament, drainage, slope, and landcover were considered. Various thematic layers of theses influencing factors were analyzed using the GIS and MIF approach to assess the GWRP zones. The methodology is given in the flowchart of Figure 4.

Multi-Influencing Factors Approach (MIF)
The MIF approach has been extensively and successfully applied by many researchers in groundwater recharge potential mapping studies [3,35]. In the MIF method, all thematic maps and their individual classes are weighted based on the literature review and according to their hydrological significance for the groundwater recharge potentiality [35]. Figure 5 illustrates the conceptual diagram of the interrelationships among the different influencing factors. From the figure, it is evident that the magnitude of the influence of each factor on the groundwater recharge potential can be computed and evaluated from the interrelationships among the different influencing factors (major and minor). A score of 1.0 was assigned for a strong relationship between the parameters, while a score of 0.5 was given if a minor effect was found, reflecting weak relationships between the factors. For instance, major interrelationships were observed for lineaments on landcover and drainage. Therefore, its evaluated weight was 2.0. This high weight value means that the factor significantly influences the recharge potentiality. The methodology starts with identifying and assigning data sources of the various conditioning factors relevant to the groundwater potential (Table 1) followed by preprocessing these sources to ensure a uniform projection and resolution. Data extraction and analysis were carried out using ArcGIS 10.7. The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER GDEM) was used with topographic and geologic maps as base maps for the delineation and evaluation of the watershed (Table 1). GDEM is processed in GIS environment by filling the sinks, flow direction, and flow accumulation grids to extract streams, watersheds, and topographic and slope features.
Landsat Enhanced Thematic Mapper (ETM+) imagery with free cloud cover scenes for the year of 2015 (Path = 98 and row = 81) during summer season were best to distinguish the spectral signatures of the different landcover types. The images were preprocessed in the ENVI 5.4 environment to reduce or eliminate haze effects before mosaicking and sub-setting. Enhancement image techniques including RGB band combinations, band rationing, and decorrelation stretch of the highest contrast, were applied to the Landsat imagery. The Line Module of PCI Geomatica 2020 Software Package was then used for the detailed automatic extraction of the lineaments. The lineament density map was derived from the Landsat images based on grid cell size of 5.0 km 2 . The landuse/landcover (LULC) features were classified using the maximum likelihood classifier in ArcGIS 10.7. The LULC categories were further verified using image interpretation elements such as drainage, tone, texture, and relief for the satellite images. The average annual rainfall data for a period of 30 years was collected from the Bureau of Meteorology, Australia from the grid file. Then, it was exported as a raster format using GIS and processed to extract the rainfall values. To overcome the low resolution of the grid data compared to the other factors' thematic layers, an equal distribution of 400 points was created over the study area, and the rainfall value was extracted and interpolated using the inverse distance weighting (IDW) method to 15 m resolution. In addition, groundwater data were retrieved from the hydrogeological database of Primary Industries and Resources, South Australia (https://map.sarig.sa.gov.au/) (accessed on 12 July 2021). This database contains information recorded by well drillers about well location, water depth, water yield, and rock type. The well database for the Flinders Ranges contains information for approximately 90% of the wells in the region. All of the thematic layers were then stored and classified in a geo-database as GIS layers. Analyses and interpretations were conducted directly from the created geo-database for the influencing factors.

Multi-Influencing Factors Approach (MIF)
The MIF approach has been extensively and successfully applied by many researchers in groundwater recharge potential mapping studies [3,35]. In the MIF method, all thematic maps and their individual classes are weighted based on the literature review and according to their hydrological significance for the groundwater recharge potentiality [35]. Figure 5 illustrates the conceptual diagram of the interrelationships among the different influencing factors. From the figure, it is evident that the magnitude of the influence of each factor on the groundwater recharge potential can be computed and evaluated from the interrelationships among the different influencing factors (major and minor). A score of 1.0 was assigned for a strong relationship between the parameters, while a score of 0.5 was given if a minor effect was found, reflecting weak relationships between the factors. For instance, major interrelationships were observed for lineaments on landcover and drainage. Therefore, its evaluated weight was 2.0. This high weight value means that the factor significantly influences the recharge potentiality. A cumulative score weight of both major and minor effects in terms of groundwater recharge potentiality was calculated to obtain the relative score for the different influenc ing factors (Table 2). This was further used to estimate the proposed score of each influ encing factor, as shown in the following equation where "S" is the proposed score of every selected influencing factor, and "A" and "B" are the major and the minor interrelationships between each two factors, respectively. The proposed score of every influencing factor was then divided and assigned to every reclas sified sub-factor (Table 3). All the parameters, with their potential weights, were reclassi fied and integrated into five classes by applying weighted overlay analysis in the GIS en vironment to determine the GWRP zones. To obtain the GWRP map, rasters with the influencing factors were overlaid using the following formula GWRP = Li × Li + Lin × Lin + As × As + Lu × Lu + Dr × Dr + Rf × Rf (2 where Li is the lithology, Lin is the lineament, As is the average slope, LU is the land use Dr is the drainage density, Rf is the rainfall, r is the rating of each class of each thematic layer, and w is the weight of each thematic layer. The obtained values for the GWRP were classified into three categories from 1 to 3. A cumulative score weight of both major and minor effects in terms of groundwater recharge potentiality was calculated to obtain the relative score for the different influencing factors (Table 2). This was further used to estimate the proposed score of each influencing factor, as shown in the following equation where "S" is the proposed score of every selected influencing factor, and "A" and "B" are the major and the minor interrelationships between each two factors, respectively. The proposed score of every influencing factor was then divided and assigned to every reclassified sub-factor (Table 3). All the parameters, with their potential weights, were reclassified and integrated into five classes by applying weighted overlay analysis in the GIS environment to determine the GWRP zones.

∑10.5 ∑100
To obtain the GWRP map, rasters with the influencing factors were overlaid using the following formula GWRP = Li r × Li w + Lin r × Lin w + As r × As w + Lu r × Lu w + Dr r × Dr w + Rf r × Rf w (2) where Li is the lithology, Lin is the lineament, As is the average slope, L U is the land use, Dr is the drainage density, Rf is the rainfall, r is the rating of each class of each thematic layer, and w is the weight of each thematic layer. The obtained values for the GWRP were classified into three categories from 1 to 3.

Verification of the GWRP
In the present study, the final map obtained using the MIF approach was validated using an independent dataset that was not used to construct the map. The observed wellyields of the water wells were obtained from the Department for Environment and Water, South Australia. The water wells were referred to as high, medium, and low water yield and used to validate the accuracy of the final groundwater recharge potential map.

Influencing Factors
The analysis of seven major factors was performed using GIS to demonstrate and evaluate their effect on the groundwater recharge potentiality.

Drainage and Drainage Density
The nature, shape, and density of drainage are essential factors that can affect the runoff, flow, and recharge of groundwater [26]. The results indicated that the watershed shows a dendritic drainage pattern; the total area (A) of the watershed is 916.4 km 2 and length (LB) is 49.5 km (Figure 6a). Using the method proposed by Strahler (1964), the watershed was found to comprise a six-order pattern with a total of 2646 streams, of which nearly 50.36% are first order, 23.15% are second order, 12.71% are third order, 7.22% are fourth order, 3.63% are fifth order and 0.03% are sixth order streams (Table 2). Lineament density is considered a useful tool to identify "hot spots" of groundwater recharge potentiality [42][43][44]. Lineament density (Ld) is defined as the total lineament length in a unit area. It is expressed as where Σ L denotes the total length of lineaments and A is the unit area. The highest densities were identified in the northeast and southwest, and around the Hawker town (Figure 6b). The lineament densities in this study area were quite high, which can be attributed to the geologic setting of the study area.

Topographic and Slope Features
The relief ratio (Rr) was introduced to express the difference between the highest and the lowest elevation of the basin [45]. The maximum elevation of the basin was 1000 m and the lowest was 250 m; therefore, the relief of the basin was 750 m. The relief ratio for this area was 0.02, which indicated a moderate topographic nature. The relief ratio and the elongated ratio of the watershed implied low resistance and more erosion and weathering processes, which are positive for groundwater potential [36].
Basin slope (Sb) is another influencing factor, which is mainly controlled by climate processes in areas outcropped by different lithologies [46,47]. In the present study, based on the degree of the slope, the watershed was classed into gentle (<10°), moderate (10-20°), and steep (>20°). Gentle slopes were classified as "good" for groundwater recharge as nearly flat terrain is the most favorable for infiltration (Figure 7a). Moderate slopes were categorized as "moderate" for groundwater recharge due to their slightly undulating topography, which permits partial runoff. Steep slopes allow high surface runoff and a negligible amount of infiltration and, therefore, were assigned to the "poor" category for groundwater recharge. The total length of all streams was 2589.6 km and ranged from 1280 km to 71 km for the first order and sixth order, respectively (Table 4). It was evident that the first order streams had the greatest total length of stream segments and, consequently, the length decreased as the stream order increased. Length of streams is a significant hydrological factor as it reflects the runoff characteristics and, hence, is related to recharging aquifers. Streams of smaller length are characterized by steep slopes, good runoff, and low recharge, while streams with longer lengths are commonly indicative of a smoother slope, low runoff, and good recharge. This change in stream length for the different orders indicated the impact of altitude, lithological variations, and slopes. More importantly, it indicated the hydrological characteristics of the underlying rock of the different stream orders in the study area [36]. Drainage density (D) was introduced by Horton (1945) as the total length of the stream per unit area. Drainage density provides a quantitative measure of the average length of streams in the whole basin. Greenbaum (1985) defined the drainage density (D d ) as the total drainage length of the drainage in a unit area. It is expressed as where ∑ i=n i=1 S i denotes the total length of the drainage, and A is the unit of the area. The calculated drainage density of the study area was 2.83 km/km 2 (Figure 6a). Drainage density values could be linked to the permeability of the dominant landcover, topography, and vegetation [37][38][39]. Based on the findings, the drainage density in the study area was classified into five categories: very good, good, moderate, poor, and very poor density. Generally, low drainage density is more likely to dominate in highly permeable, dense vegetation, and low relief areas, while high drainage density occurs in low permeable, less vegetation, and high relief areas [40]. Accordingly, high drainage areas have less water infiltration, whereas low-density areas have high water infiltration to the aquifers [41].

Lineament and Lineament Density
Two main directions, NE-SW and E-W, were identified from the images. The orientations of the lineaments matched the faults-oriented NE-SW throughout the study area. The lineament length distribution obtained from Landsat images showed that the lineament lengths were highly variable (SD = 775 m), with a mean length of about 1338 m (Figure 6b). The lineament length distribution reflects high variability in lineament length and number, which can be attributed to the geological characteristics of the study area (Figure 6b).
Lineament density is considered a useful tool to identify "hot spots" of groundwater recharge potentiality [42][43][44]. Lineament density (L d ) is defined as the total lineament length in a unit area. It is expressed as where ∑ i=n i=1 L i denotes the total length of lineaments and A is the unit area. The highest densities were identified in the northeast and southwest, and around the Hawker town (Figure 6b). The lineament densities in this study area were quite high, which can be attributed to the geologic setting of the study area.

Topographic and Slope Features
The relief ratio (Rr) was introduced to express the difference between the highest and the lowest elevation of the basin [45]. The maximum elevation of the basin was 1000 m and the lowest was 250 m; therefore, the relief of the basin was 750 m. The relief ratio for this area was 0.02, which indicated a moderate topographic nature. The relief ratio and the elongated ratio of the watershed implied low resistance and more erosion and weathering processes, which are positive for groundwater potential [36].
Basin slope (Sb) is another influencing factor, which is mainly controlled by climate processes in areas outcropped by different lithologies [46,47]. In the present study, based on the degree of the slope, the watershed was classed into gentle (<10 • ), moderate (10-20 • ), and steep (>20 • ). Gentle slopes were classified as "good" for groundwater recharge as nearly flat terrain is the most favorable for infiltration (Figure 7a). Moderate slopes were categorized as "moderate" for groundwater recharge due to their slightly undulating topography, which permits partial runoff. Steep slopes allow high surface runoff and a negligible amount of infiltration and, therefore, were assigned to the "poor" category for groundwater recharge.
water (<1%). Due to the high potentiality to recharge and the storage of groundwater, water bodies were assigned the highest weightage. Lower weightage was allocated to industrial lands due to their poor water holding capacity and less infiltration. Agricultural areas were given moderate weightage factors and native covers were allocated the highest weightage due to their ability to reduce the speed of water flow and increase the time duration for infiltration.

Rainfall
In arid areas, rainfall is considered as one of the most influential factors in groundwater recharge capability [49]. Processes such as precipitation, evapotranspiration, and percolation can control the water availability at the land surface [35]. As mentioned earlier, rainfall is the main source for the aquifers' recharge in the study area. The climate of the area is classified as arid to semi-arid with a high variation in rainfall. The average annual rainfall of the study area ranges from 35 to 50 mm. Accordingly, the annual rainfall is classified into five classes. Higher rainfall potentially infiltrates into the soil and is stored in the aquifers (Figure 7d). The rainfall distribution map indicates that the annual rainfall is higher in the high topographic areas than in the low topographic areas. Rainfall

Lithology and Landcover
In the Hawker Area, while the lithological log data for the wells was not available, the lithological grouping was decided based on detailed available geological maps. The geological map of the study area revealed that the high elevation features are covered by metasediments, quartzite, sandstone, limestone, shale, and siltstone, and the lower elevation areas are dominated by recent alluvial sediments (Figure 7b). Different rock formations were assigned rankings based on their recharge potential properties. Gener-ally, massive unfractured lithologic units have little influence on groundwater recharge potentiality, while alluvial and sandstones with secondary porosity form a good potential groundwater zone.
The landuse/cover is a significant factor that can influence the recharge, infiltration, and runoff process of the groundwater regime [48]. In Figure 7c, the primary landuse/cover is composed of grazing native vegetation (68.55%), agriculture and plantations (22.29%), conservation and natural environments (7.26%), industrial (1.90%), and water (<1%). Due to the high potentiality to recharge and the storage of groundwater, water bodies were assigned the highest weightage. Lower weightage was allocated to industrial lands due to their poor water holding capacity and less infiltration. Agricultural areas were given moderate weightage factors and native covers were allocated the highest weightage due to their ability to reduce the speed of water flow and increase the time duration for infiltration.

Rainfall
In arid areas, rainfall is considered as one of the most influential factors in groundwater recharge capability [49]. Processes such as precipitation, evapotranspiration, and percolation can control the water availability at the land surface [35]. As mentioned earlier, rainfall is the main source for the aquifers' recharge in the study area. The climate of the area is classified as arid to semi-arid with a high variation in rainfall. The average annual rainfall of the study area ranges from 35 to 50 mm. Accordingly, the annual rainfall is classified into five classes. Higher rainfall potentially infiltrates into the soil and is stored in the aquifers (Figure 7d). The rainfall distribution map indicates that the annual rainfall is higher in the high topographic areas than in the low topographic areas. Rainfall is typically regular in the hillslope areas in contrast to the foothill areas [50].

Recharge Potential Model
By classifying and implementing the MIF approach, the final map of the GWRP was produced, and grouped into three classes (good, moderate, and low) (Figure 8a). It was observed that about 28.97%, 58%, and 13.03% of the areas fall under good, moderate, and low groundwater potential zones, respectively (Figure 8b). Good groundwater potential zones were observed in the central, western, and small parts of the eastern areas due to the distribution of low slope, high drainage density, alluvial sediments, high lineament density, and natural and agricultural land. Medium and low groundwater potential recharge zones existed mainly in elevated areas in the southern and eastern parts of the study area where higher slopes and siltstone are the main lithology.

Verification of the Recharge Potential Model
For the purpose of verification, validation of the groundwater potential map was undertaken using water yield data from the study area (Figure 8b,c). The water yield data was used to statistically demonstrate the failure and success relative to each obtained class of recharge potential. The results showed that 29.29% of the wells with a high water yield (mean = 1.04 L/S) are sited on good areas, 62.86% with moderate yield (mean = 0.80 L/S) are located on moderate areas, and 7.86% are on poor groundwater potential recharge zones (Figure 8

Discussions
Influencing factors related to groundwater recharge potentiality were selected and weighted in order to identify favorable zones for recharge using the RS, GIS, and MIF approach. The findings indicated that the groundwater regime of arid regions is mainly controlled by the topography, geology, fracture systems, and climatic conditions. As shown in this study, the various groundwater potential recharge zones are distributed along the study area. The northern and western parts near lineaments with low drainage density, and where secondary structures and alluvial lithology are located, are characterized as good groundwater potential recharge areas (Figure 9). The central areas, which are dominated by a low lineament density and high slope features, are classified as low recharge potential areas. This variability in groundwater potentiality can be attributed to

Discussions
Influencing factors related to groundwater recharge potentiality were selected and weighted in order to identify favorable zones for recharge using the RS, GIS, and MIF approach. The findings indicated that the groundwater regime of arid regions is mainly controlled by the topography, geology, fracture systems, and climatic conditions. As shown in this study, the various groundwater potential recharge zones are distributed along the study area. The northern and western parts near lineaments with low drainage density, and where secondary structures and alluvial lithology are located, are characterized as good groundwater potential recharge areas (Figure 9). The central areas, which are dominated by a low lineament density and high slope features, are classified as low recharge potential areas. This variability in groundwater potentiality can be attributed to many geological factors [8,10,51]. Lineaments, as a structural weakness, could play an essential role in groundwater potential recharge zones by increasing the infiltration capacity of the soil, and facilitating the groundwater flow and movement [26,49]. In addition, the occurrence of lithologic units with high permeability such as alluvial sediments could increase groundwater storage and yields [52]. These directly impact water infiltration and percolation [23,53]. This infiltration capacity could also be influenced by the dominant slope gradient. Gentle slopes promote groundwater recharging potentiality by providing sufficient time for infiltration, while areas with steep slopes permit low recharging due to rapid flow and low infiltration capacity [15]. Drainage density is considered another significant factor in identifying recharge sites because of its influence on surface runoff. The occurrence of a dense drainage system indicates high water permeability zones in the region [54]. many geological factors [8,10,51]. Lineaments, as a structural weakness, could play an essential role in groundwater potential recharge zones by increasing the infiltration capacity of the soil, and facilitating the groundwater flow and movement [26,49]. In addition, the occurrence of lithologic units with high permeability such as alluvial sediments could increase groundwater storage and yields [52]. These directly impact water infiltration and percolation [23,53]. This infiltration capacity could also be influenced by the dominant slope gradient. Gentle slopes promote groundwater recharging potentiality by providing sufficient time for infiltration, while areas with steep slopes permit low recharging due to rapid flow and low infiltration capacity [15]. Drainage density is considered another significant factor in identifying recharge sites because of its influence on surface runoff. The occurrence of a dense drainage system indicates high water permeability zones in the region [54].
The effects of the combination of these influencing factors on controlling the groundwater regime contribute to the recharge potentiality of different aquifers. Areas with good to moderate groundwater potential zones can attribute this to the lithology, slope, and drainage density territories. This classification could, thus, provide preferred targets for groundwater exploration since it combines several hydrogeological features that could affect groundwater availability. The findings from this study are comparable to the findings of other recharge potential studies in many areas around the world [15,20,24,27,55]. However, it is important to note that the delineation of the recharge potential map using the MIF method might not be a better reflection of the realities on the ground. This The effects of the combination of these influencing factors on controlling the groundwater regime contribute to the recharge potentiality of different aquifers. Areas with good to moderate groundwater potential zones can attribute this to the lithology, slope, and drainage density territories. This classification could, thus, provide preferred targets for groundwater exploration since it combines several hydrogeological features that could affect groundwater availability. The findings from this study are comparable to the findings of other recharge potential studies in many areas around the world [15,20,24,27,55].
However, it is important to note that the delineation of the recharge potential map using the MIF method might not be a better reflection of the realities on the ground. This study was limited by the absence of recent drilling flows to validate the sensitivity of the results. However, it appeared that the drilling program was not optimized and suffered from a lack of preliminary investigation of groundwater availability. This was evident when comparing the locations of the wells throughout the study area. The siting of boreholes in relation to the local topography and, occasionally, to local fractures is the most common drilling method. In many cases, however, social factors such as the location of the house to be supplied or the boundary of the property determine the siting of the borehole. As a result, the validation of the groundwater potential map should be improved using recent drilling data in order to better evaluate the realities on the ground. In addition, the delineation of a good recharge potential zone does not assure a 100% success rate in drilling, even if the prospect appeared to be promising and verified by available water yield data. In-depth investigation is required within the proposed approach to overcome the limitation encountered from the hypothetical interpolations of surface data.

Conclusions
This study obtained a groundwater recharge potential map of the Hawker Area in the Flinders Ranges, South Australia, using a GIS, RS, and MIF approach. Analysis of the results provides ample information and analysis for decision makers. Integration of thematic maps such as drainage density, lithology, slope, landcover, rainfall, and lineament density could present valuable information for local authorities and planners regarding suitable areas for groundwater recharge potentiality in the Flinders Ranges.
The results showed that the most promising groundwater recharge potential zones are located in areas with high lineament density, low drainage density, and gentle slopes. Furthermore, the data analysis suggested that alluvial sediments and sandstone have high infiltration capacity. The areas that are least effective for groundwater recharge are underlain by shale and siltstone, and have low lineament density. This was validated by the correlation between water yield and the recharge potential categories. Through the integration of GIS, RS, and MIF, along with other hydrological characteristics of the drainage basin, the strategy of identifying and selecting groundwater recharge zones could provide better groundwater management plans.
Based on the present study, a number of future investigations are suggested. To start with, lineament zones need to be investigated in more detail. The focus should be on classifying, categorizing, and exploring these lineaments by means of drilling and field mapping both across their width and along the strike. Additionally, lithological data for most of the wells are not available; therefore, further data are required to develop a feasible interpretation of the recharge potential mapping. Moreover, a detailed investigation of the water chemistry and water level should be undertaken to understand the impacts of lithology and water rock interactions to give more meaningful results. This can be achieved by installing new piezometers in the monitoring boreholes to obtain hydraulic head data for a better understanding of the flow directions. Finally, building on different applied approaches and models such as the analytical hierarchy method (AHP) and frequency ration (FR), and optimizers such as UCODE or hydroPSO, as well as sensitivity analysis, could provide more meaningful results for practical applications [56][57][58].