The sustainable management of groundwater resources in developing countries is often challenging due to limited measurement and monitoring infrastructure to collect data necessary for decision support. To make a contribution towards addressing these challenges, this study investigated the use of Internet of Things (IoT) technology and low-cost sensors to collect the required groundwater-level data and develop a model to map the recharge potential with stormwater. The study focused on two stormwater ponds located in a highly urbanised area in Cape Town, South Africa. A combination of Geographic Information System and analytic hierarchy process was integrated to generate a groundwater recharge potential zone map of the study area. The IoT-based data were used to develop and calibrate a numerical groundwater model in MODFLOW. The study determined that retrofitted stormwater ponds are potential groundwater augmentation zones and can provide opportunity for stormwater recharge in urban areas. Overall, this study highlights the potential of IoT to collect hydrogeological data with low-cost sensors. Data can be collected at high temporal resolution, and the spatial scale can be increased due to availability of low-cost sensors.

  • Internet of Things (IoT)-based data were leveraged to address challenges in groundwater management.

  • Rapid battery drain powering the IoT system was mitigated by reducing the data collection frequency and dead sleep mechanisms.

  • Geographic Information System–based analysis and analytic hierarchy process were employed to map the recharge potential in a highly urbanised area in Cape Town.

  • This study demonstrated that stormwater ponds have the potential to recharge groundwater aquifers through infiltration.

Most countries in Southern Africa largely rely on surface water. With prolonged droughts in many countries that are expected to reoccur frequently in the future, there is a major shift to alternative water sources such as groundwater. The main challenge with the utilisation of these alternative water sources is the high monetary cost of quantity and quality monitoring coupled with limited human and financial resources. The novelty in this study is in the application of Internet of Things (IoT) for measurement and analytic hierarchy process (AHP) as a mapping approach for sustainable management of groundwater resources in developing countries. First, the application of IoT for measurement and monitoring is relatively new, while the novelty is the application in the Southern Africa context in the face of major challenges such as availability of energy and Internet; access to IoT components, i.e., sensors and gateways; safety of individuals installing the devices in the field; and vandalism.

Groundwater is the most extensive freshwater resource that is also resilient to climate change and less susceptible to pollution compared to sources such as surface water. To develop a sustainable and resilient water supply system, countries such as South Africa that are located in semi-arid areas and largely depend on surface water need to diversify water sources to include groundwater to meet future demands (DWA 2010). The need to diversify water sources is more urgent in urban areas to meet the rapidly increasing demand attributed to population growth and improving living standards due to rising incomes. In the city of Cape Town (CCT) in South Africa, identifying opportunities for managed aquifer recharge is crucial to sustain groundwater systems (Mauck 2017). Although groundwater resources are extensive and their resilience to climate change is well known as climate change projections for Southern Africa indicate an expectation of hotter and drier weather, particularly in the western part of the region, they need to be managed sustainably and preserved to meet the water demands of the future (DWA 2010). Unfortunately, efforts to sustainably manage groundwater resources are often hindered by limited data and inadequate monitoring capacity to support decision-making (Gaffoor et al. 2020). High-quality data are essential to determine the variability and viability of the groundwater resource (Cobbing & Hiller 2019). IoT technology has provided the opportunity to remotely collect data in real time with high spatial and temporal resolution that can support decision-making and enhance the management of the groundwater resources (Gaffoor et al. 2020). In this study, an investigation was undertaken to assess the application of IoT-based data and the use of AHPs to map groundwater recharge potential around two stormwater ponds in a small, urbanised area in Cape Town. The study involved the collection of groundwater-level data using an IoT monitoring system and then using the Geographic Information System (GIS) and AHP to generate groundwater recharge potential zones in the study area. Using a groundwater model in MODFLOW, the study also estimated the recharge rates and processes in the aquifer.

Study area and data collection

The study focused on two stormwater ponds located south of the CCT, i.e., the ‘School Pond’ at 34°2′3.0″S; 18°35′4.7″E and the ‘Green Dolphins Pond’ at 34°2′14.08″S, 18°35′11.57″E (Figure 1).
Figure 1

Highly urbanised land use in the study area.

Figure 1

Highly urbanised land use in the study area.

Close modal

The topography is relatively flat (less than 3% average slope) with an unconfined aquifer that offers opportunity for groundwater recharge (Adelana 2010; Mauck 2017). The study area is 0.315 km2 comprising the catchment area draining into the stormwater ponds and surrounding areas.

Selected stormwater ponds

The stormwater ponds were mainly selected as suitable sites for the study due to retrofitting works undertaken to enhance groundwater recharge with stormwater in a project linked to this research (Okedi 2019; Tanyanyiwa et al. 2023). Both studies investigated the dynamic aspects of surface–groundwater interaction associated with managed aquifer recharge (MAR) in the study area. This investigation builds on these studies with focus on the role of IoT for measurement and AHP for mapping MAR. The geometry and land use around the two stormwater ponds is also shown in Figure 1. The catchment area of the stormwater ponds was delineated using a 0.5 m high-resolution digital elevation model (DEM) created from a Lidar dataset. It was determined that the ‘School Pond’ and ‘Green Dolphins Pond’ had catchment areas of 0.155 and 0.048 km2, respectively. However, the dataset also captured the building tops with the highest elevation of 47.81 m. A corrected DEM was developed using the Golden Surfer Software and the Kriging method in the ArcGIS software using the contours of the study area in order to obtain the ground-level values used in the modelling (Figure 2).
Figure 2

DEM showing the corrected ground-level values.

Figure 2

DEM showing the corrected ground-level values.

Close modal

Climatic conditions of the study area

Daily rainfall and temperature data were collected from two weather stations near the study area, i.e., Cape Town Automatic Weather Station (CT-AWS) and Kirstenbosch, to model the hydrology of the study area. The CT-AWS station was chosen for further analysis due to its proximity to the study area. The mean monthly rainfall and temperatures and the annual precipitation recorded are shown in Figures 3 and 4.
Figure 3

Mean monthly temperature and rainfall (2002–2021) at the CT_AWT station.

Figure 3

Mean monthly temperature and rainfall (2002–2021) at the CT_AWT station.

Close modal
Figure 4

Annual precipitation (2002–2021) at the CT_AWT Station.

Figure 4

Annual precipitation (2002–2021) at the CT_AWT Station.

Close modal
For estimation of evaporation, the evaporation pan and the Food Agriculture Organisation (FAO) Penman–Monteith equation methods were considered. A weather station near the study area was appropriate as it measured data using the evaporation pan. The station collected data remotely through a cellular network, downloaded daily at midnight, and stored in a server database. The data were utilised to convert evaporation to evapotranspiration (ETo) using Equation (1).
(1)
where is the evapotranspiration (mm/day), is the Class A pan coefficient, and is the evaporation data (mm/day).
The pan evapotranspiration method connects evaporation to ETo and is suitable for regions with limited meteorological data. Doorenbos & Pruitt (1977) presented the Kp values for various weather conditions based on different pan locations as shown in Table 1. The mean daily wind speed and the relative humidity data recorded at the Cape Town Airport located about 7.2 km northeast of the study area were obtained from a study by Okedi (2019) as shown in Figures 5 and 6. It is observed from the mean wind speed during May, June, July, and August that it ranged between 175 and 425 km/day, and the mean relative humidity was >70%. A Kp value of 0.8 was used to estimate the ETo, given a windward distance of 1,000 m. The estimated mean ETo was determined as shown in Table 2.
Table 1

Pan coefficient (Kp) for Class A pan (Doorenbos & Pruitt 1977)

Class A panPan placed in short green cropped area
Wind speed at a height of 2 m (km/day)Windward side distance of green crop (m)Mean daily relative humidity (%)
< 40 (40*)40–70 (55*)> 70 (70*)
<175 (175*) 0.55 0.65 0.75 
10 0.65 0.75 0.85 
100 0.70 0.80 0.85 
1,000 0.75 0.85 0.85 
175–425 (300*) 0.50 0.60 0.65 
10 0.60 0.70 0.75 
100 0.65 0.75 0.80 
1,000 0.70 0.80 0.80 
425–700 (562*) 0.45 0.50 0.60 
10 0.55 0.60 0.65 
100 0.60 0.65 0.70 
1,000 0.65 0.70 0.75 
> 700 (700*) 0.40 0.45 0.50 
10 0.45 0.55 0.60 
100 0.50 0.60 0.65 
1,000 0.55 0.60 0.65 
Class A panPan placed in short green cropped area
Wind speed at a height of 2 m (km/day)Windward side distance of green crop (m)Mean daily relative humidity (%)
< 40 (40*)40–70 (55*)> 70 (70*)
<175 (175*) 0.55 0.65 0.75 
10 0.65 0.75 0.85 
100 0.70 0.80 0.85 
1,000 0.75 0.85 0.85 
175–425 (300*) 0.50 0.60 0.65 
10 0.60 0.70 0.75 
100 0.65 0.75 0.80 
1,000 0.70 0.80 0.80 
425–700 (562*) 0.45 0.50 0.60 
10 0.55 0.60 0.65 
100 0.60 0.65 0.70 
1,000 0.65 0.70 0.75 
> 700 (700*) 0.40 0.45 0.50 
10 0.45 0.55 0.60 
100 0.50 0.60 0.65 
1,000 0.55 0.60 0.65 

*The mean value.

Table 2

Estimation of the reference evapotranspiration using A-pn evaporation

MonthMonthly A-pan evaporation (mm)Mean monthly A-pan evaporation (Epan) (mm/day)Pan coefficient (Kp)Mean monthly evapotranspiration (ETo) (mm/day)
May 44.5 1.44 0.80 1.15 
June 52.8 1.76 0.80 1.41 
July 44.8 1.45 0.80 1.16 
August 64.0 2.06 0.80 1.65 
Average 1.34 
MonthMonthly A-pan evaporation (mm)Mean monthly A-pan evaporation (Epan) (mm/day)Pan coefficient (Kp)Mean monthly evapotranspiration (ETo) (mm/day)
May 44.5 1.44 0.80 1.15 
June 52.8 1.76 0.80 1.41 
July 44.8 1.45 0.80 1.16 
August 64.0 2.06 0.80 1.65 
Average 1.34 
Figure 5

Mean daily wind speeds at Cape Town Airport (Okedi 2019).

Figure 5

Mean daily wind speeds at Cape Town Airport (Okedi 2019).

Close modal
Figure 6

Mean daily relative humidity at Cape Town Airport (Okedi 2019).

Figure 6

Mean daily relative humidity at Cape Town Airport (Okedi 2019).

Close modal
Figure 7

Reference evapotranspiration for Cape Town Airport (Adelana 2010).

Figure 7

Reference evapotranspiration for Cape Town Airport (Adelana 2010).

Close modal
Adelana (2010) calculated the daily ETo in Cape Town from 2000 to 2009 using the FAO Penman–Monteith equation as follows:
(2)
where ETo is the reference evapotranspiration (mm day−1), Rn is the net radiation at the crop surface (MJ m−2 day−1), G is the soil heat flux density (MJ m−2 day−1), T is the mean daily air temperature at 2 m height (°C), u2 is the wind speed at 2 m height (m s−1), es is the saturation vapour pressure (kPa), ea is the actual vapour pressure (kPa), esea is the saturation vapour pressure deficit (kPa), Δ is the slope vapour pressure curve (kPa °C−1), and Γ is the psychrometric constant (kPa °C−1). The estimated average monthly values are as presented in Figure 7.

The maximum daily ETo value was obtained as 3.28 mm/day (3.8 × 10−8m/s). This ETo value was preferred to the pan evapotranspiration value due to limited data caused by vandalism of the weather station. An average ETo value of 2.59 mm/day (3.0 × 10−8 m/s) was used in the model.

Mapping groundwater recharge potential

A combination of GIS and AHP was used to delineate and map the groundwater recharge zone. Seven thematic maps were considered the control factors of groundwater flow and storage in the area. These influencing factors were weighted according to their reaction to groundwater occurrence, expert opinion, and consideration from the reviews from previous studies. The factors and their importance are compiled from previous literature, and appropriate scores were assigned to each factor to produce a final groundwater recharge potential map of the area (Yeh et al. 2016). Various studies have been carried out globally using the integration of geospatial data, multi-criteria decision analysis techniques, remote sensing, and GIS technology to identify suitable sites for groundwater recharge potential through MAR suitability maps (Sallwey et al. 2019; Achu et al. 2020; Ahmed et al. 2021). The assigned values were chosen in comparison to studies like Sener & Davraz (2013) and Arulbalaji et al. (2019). Expert opinion is valuable but introduces subjectivity; therefore, the methodology was subjected to peer review for transparency and credibility.

The process required slope information that was generated from the Shuttle Radar Topography Mission (SRTM) data, drainage density, and lineament density maps that were generated from the DEM created using Lidar dataset. Land use/land cover (LULC) was determined through supervised image classification and the geology, lithology, and soil shapefiles of the Cape Flats Aquifer used to extract the shapefile of the study area using the extraction by mask tool in ArcGIS software. From these data, seven thematic layers/factors, including geology, lineament density, lithology, slope, soil, LULC, and drainage density, were analysed, weighted according to their reaction to groundwater occurrence, and assigned weights from Saaty's scale of relative importance (Table 3).

Table 3

Saaty's scale of relative importance (Saaty 1990)

Intensity of importanceDefinitionExplanation
Equal importance Two elements contribute equally to the objective 
Moderate importance Experience and judgement slightly favour one element over another 
Strong importance Experience and judgement strongly favour one element over another 
Very strong importance One element is favoured very strongly over another, its dominance is demonstrated in practice 
Extreme importance The evidence favouring one element over another is of the highest possible order of affirmation 
Intensity of importanceDefinitionExplanation
Equal importance Two elements contribute equally to the objective 
Moderate importance Experience and judgement slightly favour one element over another 
Strong importance Experience and judgement strongly favour one element over another 
Very strong importance One element is favoured very strongly over another, its dominance is demonstrated in practice 
Extreme importance The evidence favouring one element over another is of the highest possible order of affirmation 

Note: 2, 4, 6, and 8 can be used to express intermediate values.

The thematic layers were compared with each other using a pair-wise comparison matrix as shown in Table 4.

Table 4

Pair-wise comparison matrix between applied parameters for the AHP model

GeologyLithologySlopeDrainage densityLand useLineament densitySoil
Geology 
Lithology 
Slope 1/3 1/3 
Drainage density 1/3 1/3 
LULC 1/5 1/5 1/3 
Lineament density 1/5 1/5 1/3 1/2 
Soil 1/5 1/3 1/5 1/3 1/3 
Total 10 16 18 21 
GeologyLithologySlopeDrainage densityLand useLineament densitySoil
Geology 
Lithology 
Slope 1/3 1/3 
Drainage density 1/3 1/3 
LULC 1/5 1/5 1/3 
Lineament density 1/5 1/5 1/3 1/2 
Soil 1/5 1/3 1/5 1/3 1/3 
Total 10 16 18 21 

A normalised pair-wise comparison matrix was created and the normalised weights for each component were calculated using the average of each row as shown in Table 5.

Table 5

Normalised pair-wise comparison matrix and weights of each factor

GeologyLithologySlopeDrainage densityLand useLineament densitySoilWeight
Geology 0.31 0.29 0.34 0.31 0.31 0.28 0.24 0.30 
Lithology 0.31 0.29 0.34 0.31 0.31 0.28 0.14 0.28 
Slope 0.10 0.10 0.11 0.10 0.18 0.17 0.24 0.14 
Drainage density 0.10 0.10 0.11 0.10 0.06 0.11 0.14 0.10 
LULC 0.06 0.06 0.04 0.10 0.06 0.06 0.14 0.07 
Lineament density 0.06 0.06 0.04 0.05 0.06 0.06 0.05 0.05 
Soil 0.06 0.10 0.02 0.03 0.02 0.06 0.05 0.05 
Total 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 
GeologyLithologySlopeDrainage densityLand useLineament densitySoilWeight
Geology 0.31 0.29 0.34 0.31 0.31 0.28 0.24 0.30 
Lithology 0.31 0.29 0.34 0.31 0.31 0.28 0.14 0.28 
Slope 0.10 0.10 0.11 0.10 0.18 0.17 0.24 0.14 
Drainage density 0.10 0.10 0.11 0.10 0.06 0.11 0.14 0.10 
LULC 0.06 0.06 0.04 0.10 0.06 0.06 0.14 0.07 
Lineament density 0.06 0.06 0.04 0.05 0.06 0.06 0.05 0.05 
Soil 0.06 0.10 0.02 0.03 0.02 0.06 0.05 0.05 
Total 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 

A parameter with a high weight illustrates a layer with a high impact on groundwater recharge and vice versa (Arulbalaji et al. 2019). The higher the normalised principal eigenvector, the greater the influence the parameter has on the groundwater recharge. The matrix was checked for consistency by determining the consistency ratio (CR).
where CI is the consistency index given by Equation (3) and RCI is the random consistency index.
(3)
where is the principal eigenvalue and n is the number of parameters.

A pair-wise matrix is considered consistent when is equal to or more than the number of the layers being examined. Given that of 7 and n of 7 were obtained, a CI of 0 was calculated. The RCI was obtained from Table 6 and the CI was calculated as 0. Therefore, this analysis was considered since Saaty (1990) stated that a CR of equal to or less than 0.1 for larger matrices is acceptable to continue the analysis.

Table 6

Saaty's ratio index for the different values of N

The consistency indices of randomly generated reciprocal matrices
Order of the matrix
N123456789101112
RCI value 0.00 0.00 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51 1.48 
The consistency indices of randomly generated reciprocal matrices
Order of the matrix
N123456789101112
RCI value 0.00 0.00 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51 1.48 

The sub-classes of the layers were re-classified using the natural breaks classification method for assigning weights and ranked according to their relative influence on groundwater development. The thematic layers were integrated with the weighted overlay analysis method using Equation (4).
(4)
where GWPZ is the groundwater potential zone, X is the weight of the thematic layers, Y is the rank of the sub-classes of the thematic layer, A is (1, 2, 3 … X), and B is (1,2,3 … Y).

The final weight and ranking for the different sub-classes of the thematic layers are given in Table 7.

Table 7

Weight and ranking for the different sub-classes of the thematic layers

ParameterSub-classesWeightInfluence (%)RankArea (km2)
Geology Quaternary 0.30 30 0.315 
Lithology Sand 0.28 28 0.315 
Slope (degrees) 0–0.15 0.14 14 0.161 
0.16–0.47 0.007 
0.48–0.82 0.078 
0.83–1.2 0.048 
1.3–1.5 0.020 
Drainage density (km/km20–27.52 (low) 0.10 10 0.234 
27.53–55.04 (very low) 0.057 
55.05–82.56 (moderate) 0.018 
82.57–110.1 (high) 0.005 
110.2–137.6 (very high) 0.001 
LULC Bare ground 0.07 0.008 
Storm water pond 0.019 
Vegetation 0.002 
Built-up area 0.223 
Road 0.064 
Lineament density (km/km20–15.92 (low) 0.06 0.261 
15.93–31.85 (very low) 0.033 
31.86–47.77 (moderate) 0.016 
47.78–63.7 (high) 0.005 
63.71–79.62 (very high) 0.001 
Soil Sandveld 0.05 0.315 
ParameterSub-classesWeightInfluence (%)RankArea (km2)
Geology Quaternary 0.30 30 0.315 
Lithology Sand 0.28 28 0.315 
Slope (degrees) 0–0.15 0.14 14 0.161 
0.16–0.47 0.007 
0.48–0.82 0.078 
0.83–1.2 0.048 
1.3–1.5 0.020 
Drainage density (km/km20–27.52 (low) 0.10 10 0.234 
27.53–55.04 (very low) 0.057 
55.05–82.56 (moderate) 0.018 
82.57–110.1 (high) 0.005 
110.2–137.6 (very high) 0.001 
LULC Bare ground 0.07 0.008 
Storm water pond 0.019 
Vegetation 0.002 
Built-up area 0.223 
Road 0.064 
Lineament density (km/km20–15.92 (low) 0.06 0.261 
15.93–31.85 (very low) 0.033 
31.86–47.77 (moderate) 0.016 
47.78–63.7 (high) 0.005 
63.71–79.62 (very high) 0.001 
Soil Sandveld 0.05 0.315 

IoT-based groundwater-level data

The study area had an existing long-range wide area network (LoRaWAN) coverage; therefore, a SenseCAP LoRaWAN wireless sensor network, with a water-level sensor, node, gateway, and smartphone were used for monitoring and analysing groundwater levels. LoRa technology is distinguished by its extensive coverage, capable of transmitting data over long distances, provided there is a clear line of sight between the gateway and the sensor. LoRaWAN, as a network protocol, leverages LoRa technology for low-power wide area network communications and device management. In rural settings, research indicates that LoRa typically provides coverage of approximately 20 km, whereas in urban areas, LoRa can achieve ranges of up to 5 km (Sanchez-Iborra et al. 2018; Benites et al. 2019; Sendra et al. 2022).

A Seeed SenseCAP LoRaWAN gateway and a Hefei WNK8010 submersible water-level sensor were used. The gateway was installed at The Leadership College-Hyde Park, next to the School pond. The water-level sensor was installed in the secure CCT MON71 borehole at the Green Dolphins Pond at coordinates 34°2′15.42″S latitude and 18°35′11.42″E longitude, around 0.35 km away from the gateway as shown in Figure 8. The borehole consisted of a 165 mm PVC pipe drilled to the basement of the aquifer at a depth of 45 m and a 305 mm steel casing around the PVC pipe to a depth of top 5 m (Figure 9). The gateway was connected to the Internet using a 4G LTE micro-SIM card with a monthly subscription. An uninterruptible power supply was connected to the gateway as a power backup to ensure continuous communication between the sensor and the gateway.
Figure 8

Location of gateway and water-level sensor.

Figure 8

Location of gateway and water-level sensor.

Close modal
Figure 9

MON71 borehole.

The gateway was registered on the ThingSpeak website that stores and displays data using a LoRa application, allowing tracking, analysing, and visualising water-level data in real time on the cloud. The sensor node comprised four components, i.e., microcontroller unit (MCU), sensor, printed circuit board (PCB), and battery as the power supply. All the components were placed in an enclosed compartment for safety and waterproofing purposes. The submersible water-level sensor was used to send raw data to the gateway. The sensor comprises a 6 m long cable and a corrosion-resistant stainless steel housing material. It was connected to the node using the terminal blocks on the Open Water Network PCB to create a circuit for the measurement and transfer of data from the sensor to the gateway where it was stored and made available for processing. The sensor has a measurement range of 0–5 m and converts analogue current to digital current that the MCU can process. A code was developed in MicroPython and used to process the data. The sensor has an onboard voltage divider that converts electrical resistance into water-level depth. It transfers the pressure of the water to the current signal with its reliable amplified circuit and precision temperature compensation using Equations (5) and (6).
(5)
(6)
where P is the liquid pressure, is the liquid density, g is the gravity, h is the depth (distance from surface of the liquid to the sensor), and Po is the atmospheric pressure at liquid surface. The voltage values were converted into depth values using Equation (7).
(7)
where h is the liquid level (m), V is the voltage output signal (volts), is the voltage output corresponding to a depth of 0 m, is the voltage output corresponding to a depth of 5 m, is the density of the liquid (kg/m3), and R is the maximum range of height (m).

The water-level data on The Things Network (TTN) were analysed and visualised using MATLAB and the ThingSpeak application. A channel was created on the ThingSpeak website to generate a right application programme interface (API) to connect the TTN and ThingSpeak and perform an in-depth analysis of the data, view/graph the necessary variables, send alerts, and export the data to Microsoft Excel. The IoT-based data were validated using the manually recorded water-level data at the MON71.

Groundwater model

A numerical model was developed using MODFLOW NWT, with ModelMuse 4 as the graphical and visual user interface to study groundwater flow and water levels. MODFLOW NWT was used with the Upstream-Weighting (UPW) Package. The conceptual model of the unconfined aquifer was divided into two layers based on the subsurface geology of the study area. The first layer (Layer 1) was extended to 20 m below the surface and the second layer (Layer 2) from 20 to 45 m. The top surface of the double-layered model represented the land surface, while the elevation of the bottom cells was determined using the thickness of the layers. The layers were assigned the same aquifer parameters except for the hydraulic conductivity that varied between layers. Different Kx values were assigned as 1.2 × 10−4 and 6.0 × 10−4 m/s for Layers 1 and 2, respectively, given that the coarse sand in Layer 2 has a higher conductivity than the fine to medium sand in Layer 1. As initial conditions, the conductivity in the z-direction (Kz) was assumed to be one-tenth of the x direction (Kx), and that in the y-direction (Ky) was made equal to Kx. The Sy was estimated based on the Theis solution by an excel spreadsheet and Aqua test analysis by Gxokwe & Xu (2017) that determined average storativity of the Cape Flat Aquifer (CFA) to range from 0.004 to 0.2 with the Phillipi borehole having storativity values of 0.1. The IoT-based water-level data collected were used to specify the initial water levels of the model. The groundwater-level data collected for almost 5 months was displayed as graphs on the ThingSpeak channel in real time. The Hydrognomon software was used to plot a time series of the IoT-based groundwater levels from July to December 2022 as shown in Figure 10. Hydrognomon is a feely available software tool designed for the analysis of hydrological time series data. This software operates under the GNU General Public License version 3 (GNU GPLv3) and is integrated into the openmeteo.org framework, with accessible source code (Kozanis et al. 2010).
Figure 10

IoT-based groundwater-level time series (m.a.s.l).

Figure 10

IoT-based groundwater-level time series (m.a.s.l).

Close modal

The dipped groundwater-level data recorded monthly at the MON71 borehole in August, September, and October 2022 were also obtained and compared with the IoT-based groundwater levels. The data showed an acceptable difference of less than 0.1 m. It was determined that IoT is reliable in collecting large amount of data while maintaining the desired accuracy levels, which solves the problem of insufficient hydrogeological data. It is possible to employ IoT-based systems to collect high-resolution data for use in developing and calibrating models for groundwater studies.

In the first 28 days, data were collected at 15-min intervals, which caused the battery voltage to drop below the required threshold of 3.5 mV, leading to a communication failure on 16 August 2022. A new battery was installed on 22 August 2022, and data collection continued at a 24-h resolution as the groundwater levels did not exhibit significant changes within a sub-daily time step. The battery ran out on 26 October 2022 and a replacement battery was installed on 15 November 2022, resulting in a noticeable increase in battery voltage and a corresponding jump in the groundwater data. It was determined that reducing the temporal resolution and frequency of data collection prolongs the battery life.

The water-level depths were converted to groundwater levels and imported into Hydrognomon software to generate time series plots. The average water levels were considered in the steady state flow regime model simulation and calibration. The results from the steady state were used as the initial conditions of the transient state simulation. The model employed the constant head boundary (CHB) and the general head boundary (GHB). The western boundary was set as the CHB with a boundary head of 31 m, while the eastern and northern boundaries were set as the GHB with a head of 35 m based on the DEM elevations. Xu & Beekman (2019) presented the typical results of recharge studies, as shown in Figure 11. In this study, the estimation was based on a study by Van Wyk et al. (2011) that provided results of recharge rates in Southern Africa. The recharge rate from precipitation, estimated as 11.3 mm/year (3.6 × 10−10m/s) using the average annual precipitation of 434 mm, calculated from the selected time series, was assigned as the initial value using the recharge package (RCH).
Figure 11

Typical results of recharge rates in Southern Africa (Van Wyk et al. 2011).

Figure 11

Typical results of recharge rates in Southern Africa (Van Wyk et al. 2011).

Close modal

The drain package (DRN) was used to simulate the amount of stormwater runoff draining into the ponds to facilitate aquifer recharge. The drainage lines shapefile was assigned to the Model_Top with a conductance per unit length of 0.001. The steady state was run with a single time step, while the transient state was run for 20–140 days.

Model calibration and sensitivity analysis

The recharge rates were adjusted according to the study by Xu & Beekman (2019), which classified the recharge rates in Cape Town in the range of 11–54 mm/year. The Kx values of the layers were adjusted within the allowable range of 10–50 m/day as specified by the CFA hydraulic conductivity profile. The evapotranspiration rates were calculated using the FAO Penman–Monteith equation, and the aquifer thickness was adjusted with reference to the CFA thickness. The steady and transient state models were calibrated by visually comparing simulated heads at MON71 and the IoT-based observed data. The root mean square error (RMSE) method was used to determine the accuracy of the model calibration.

The sensitivity analysis was carried out to determine the effect of varying different parameters and stresses on the model calibration process. The model showed that varying recharge rates, hydraulic conductivity, and aquifer thickness influences the fluctuations in the water table and the outflows presented in the water budget tables. An increase in the recharge rate registered an increase in the water table and outflow values; thus, groundwater levels are directly proportional to recharge. The results from the abstraction rate showed that groundwater levels are inversely proportional to the abstraction rate; as an increase in the abstraction rate decreased the water table, more water was taken from the aquifer. A reduction in the hydraulic conductivity in Layer 2 caused a decrease in the water table, while an increase in the hydraulic conductivity in Layer 1 led to an increase in the water table.

According to Dlodlo (2012), linking IoT to sustainable water resources management is very attractive and utilised in the region of water resources management. In this study, IoT-based groundwater-level data were collected using an IoT monitoring system that was easy to use and had the potential to accurately analyse and process large amounts of data used to calibrate the MODFLOW model. Therefore, IoT turned out to be a standout as a groundwater data collection system (Oke et al. 2019).

Zhang et al. (2019) combined AHP and GIS-based analysis with groundwater modelling to identify sites suitable for MAR implementation along the West Coast of South Africa. In this study, the groundwater recharge zone map was prepared by overlaying the cumulative weights of the thematic layers developed with AHP, a multi-criteria decision-making technique widely used for groundwater management (Saaty 1980). The thematic layers are assigned weights to determine the recharge potential score (Pr) and subjected to the weighted overlay analysis in GIS (Yeh et al. 2016; Abijith et al. 2020). It was noted that groundwater levels were highly influenced by the location of the retrofitted ponds and seasonality. High groundwater levels were observed around the stormwater ponds in the winter months of June, July, and August due to the heavy rainfall. Lower water levels were recorded during the summer months from December due to high temperature and low rainfall. As presented in Figure 12, the areas around the two stormwater ponds showed excellent groundwater recharge potential.
Figure 12

Groundwater recharge zone map of the study area.

Figure 12

Groundwater recharge zone map of the study area.

Close modal

Effect of stormwater recharge on groundwater levels

Model calibration is carried out using the trial-and-error method of adjusting the hydrogeological input parameters to represent the field conditions within the acceptable criterion (Khadri & Pande 2016). The steady state model was calibrated to achieve simulated heads close to the field-measured values through the trial-and-error method. The calibration considered the IoT-based groundwater-level data recorded in September as the observed head. The model was run to obtain the results and giving a good relationship between the simulated and observed head at the MON71, with an RMSE value of 0.65. The higher hydraulic head values in Layer 1 were attributed to the lower hydraulic conductivity than those in Layer 2. The higher heads were located in high elevation areas and flow towards the low elevation areas that comprise the stormwater ponds. For Layer 2, the presence of coarser sands contributed to the allocation of a higher hydraulic conductivity than in Layer 1, leading to a lower range of hydraulic heads (Figure 13). The water table is higher in high elevation areas than low elevation areas, especially around the school pond catchment area as shown in Figure 14.
Figure 13

Steady state hydraulic head values for Layer 1 (left) and Layer 2 (right).

Figure 13

Steady state hydraulic head values for Layer 1 (left) and Layer 2 (right).

Close modal
Figure 14

Water table contours around the study area.

Figure 14

Water table contours around the study area.

Close modal
The transient state model considered varying the recharge rates, hydraulic conductivity, aquifer thickness, and stormwater pond influence on the water table and outflows. The model was calibrated using IoT-based water level measured over the wet months of July, August, and September 2022. The steady state model, as the initial condition, was used to determine the reliability of the parameters described. An abstraction rate of 32 L/s, specific yield of 0.01, and a recharge rate of 11.3 mm/year were applied to the 40-m-thick transient model to obtain a percent discrepancy of 0.0%. The difference between the inflow and outflow storage values was attributed to the recharge at the stormwater ponds. The values were checked, and the model was observed to be more stable with lower hydraulic conductivity values in Layer 1 than in Layer 2. The final calibrated model confirmed that the residual value between the observed and simulated heads had an acceptable RMSE value of 0.86. In this study, the acceptable RMSE value was between 0 and 1; otherwise, the input parameters were changed and the model ran until a suitable value was obtained (Pisinaras et al. 2007). The heads for each layer were categorised into five classes. Figure 15 shows the simulated hydraulic heads for Layer 1 (left) and Layer 2 (right). The hydraulic heads reduce towards the location of the stormwater ponds due to elevation differences. The hydraulic head values of Layer 2 indicate that the presence of coarser sands contributed to the allocation of a higher hydraulic conductivity than in Layer 1, leading to the low values of hydraulic heads.
Figure 15

Transient state hydraulic head values for Layer 1 (left) and Layer 2 (right).

Figure 15

Transient state hydraulic head values for Layer 1 (left) and Layer 2 (right).

Close modal

This study showed that IoT-based data can be leveraged to address challenges and limitations prevalent in groundwater resource management in Southern Africa. IoT data can be used to effectively address the lack of high-quality groundwater data, enabling accurate collection, analysis, and processing of groundwater levels. Battery drain was mitigated by reducing data collection frequency, ensuring prolonged sensor battery life. The study also utilised GIS-based analysis, AHP techniques, and MODFLOW software to identify and map the recharge potential of two stormwater ponds. The difference between the observed and simulated heads at the MON71 borehole was 0.86 m. The results from the modelling exercise confirmed the occurrence of recharge at the stormwater pond, validating the effectiveness of stormwater ponds for groundwater augmentation. It was also determined that comprehensive hydrogeological assessments including GIS and AHP techniques were essential for mapping recharge zones and creating groundwater augmentation suitability maps. Stormwater ponds were classified in the excellent groundwater recharge zone, demonstrating their potential for groundwater replenishment through infiltration. Fluctuations in the water table and outflows were influenced by recharge rates, hydraulic conductivity, and aquifer thickness. Higher recharge rates resulted in increased water table and outflow values, establishing a direct relationship between groundwater levels, recharge rates, and abstraction for effective groundwater resource management.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Abijith
D.
, Saravanan, S., Singh, L., Jennifer, J. J., Saranya, T. & Parthasarathy, K. S. S.
2020
GIS-based multi-criteria analysis for identification of potential groundwater recharge zones – a case study from Ponnaniyaru watershed, Tamil Nadu, India
.
HydroResearch
(
3
),
1
14
.
Achu
A. L.
,
Reghunath
R.
&
Thomas
J.
2020
Mapping of groundwater recharge potential zones and identification of suitable site-specific recharge mechanisms in a tropical river basin
.
Earth Systems and Environment
4
(
1
),
131
145
.
doi: 10.1007/s41748-019-00138-5
.
Adelana
S.
2010
Groundwater Resource Evaluation and Protection in the Cape Flats, South Africa. LAMBERT Academic Publishing GmbH & Co., Germany.
Ahmed
A.
,
Alrajhi
A.
&
Alquwaizany
A. S.
2021
Identification of groundwater potential recharge zones in Flinders Ranges, South Australia using remote sensing, GIS, and MIF techniques
.
13
. https://doi.org/10.3390/ w13182571
.
Arulbalaji
P.
,
Padmalal
D.
&
Sreelash
K.
2019
GIS and AHP techniques based delineation of groundwater potential zones: A case study from Southern Western Ghats, India
.
Scientific Reports
9
(
1
),
1
17
.
doi: 10.1038/s41598-019-38567-x
.
Benites
B.
,
Chávez
E.
,
Medina
J.
&
Vidal R
C. M.
2019
LoRaWAN applied in swarm drones: A focus on the use of fog for the management of water resources in Lima-Peru
. In
Proceedings of the 5th International Conference on Mechatronics and Robotics Engineering
,
Rome, Italy
, pp.
171
176
.
Cobbing
J.
&
Hiller
B.
2019
Waking a sleeping giant: Realizing the potential of groundwater in sub-Saharan Africa
.
World Development
122
,
597
613
.
doi: 10.1016/j.worlddev.2019.06.024
.
Dlodlo
N.
2012
Adopting the Internet of Things technologies in environmental management in South Africa
. In
International Conference on Environment Science and Engineering (IPCBEE)
,
Vol. 32
, pp.
45
55
.
Doorenbos
J.
&
Pruitt
W. O.
1977
Crop water requirements. Irrigation and Drainage Paper
. 24, FAO, Rome, Italy.
DWA
2010
Groundwater Strategy 2010
. p.
64
.
Gaffoor
Z.
, Pietersen, K., Jovanovic, N., Bagula, A. & Kanyerere, T.
2020
Big data analytics and its role to support groundwater management in the Southern African development community
.
Water (Switzerland)
12
(
10
).
doi: 10.3390/w12102796
.
Gxokwe
S.
&
Xu
Y.
2017
Conceptualization of urban hydrogeology within the context of water sensitive urban design: Case study of Cape Flats Aquifer (May). Available from: http://etd.uwc.ac.za/.
Khadri
S. F. R.
&
Pande
C.
2016
Ground water flow modeling for calibrating steady state using MODFLOW software: A case study of Mahesh River basin, India
.
Modeling Earth Systems and Environment
2
(
1
),
1
17
.
doi: 10.1007/s40808-015-0049-7
.
Kozanis
S.
, Christofides, A., Mamassis, N. & Efstratiadis, A.
2010
Hydrognomon – Open source software for the analysis of hydrological data
.
European Geosciences Union General Assembly
12
(
May 2010
),
12419
.
doi:10.13140/RG.2.2.21350.83527
.
Mauck
B. A.
2017
The Capacity of the Cape Flats Aquifer and Its Role in Water Sensitive Urban Design in Cape Town
. PhD thesis available at https://open.uct.ac.za/items/c4fdebdd-4ce4-4d16-87cd-d02323fee53e.
Oke
S. A.
,
Alowo
R.
&
Masinde
M.
2019
Using Internet of Things for sustainable groundwater drought management in the Modder River catchment of South Africa, Open Innovations Conference, OI 2019, October, pp. 63–70. doi: 10.1109/OI.2019.8908246
.
Okedi
J.
2019
The Prospects for Stormwater Harvesting in Cape Town, South Africa Using the Zeekoe Catchment as A Case Study
, p.
227
.
Pisinaras
V.
, Petalas, C., Tsihrintzis, V. A. & Zagana, E.
2007
A groundwater flow model for water resources management in the Ismarida plain, North Greece
.
Environmental Modeling and Assessment
12
(
2
),
75
89
.
doi: 10.1007/s10666-006-9040-z
.
Saaty
T. L.
1980
The Analytic Hierarchy Process: Planning, Priority Setting, Resource Allocation
.
McGraw-Hill International Book Co
,
New York, NY, USA.
Saaty
T. L.
1990
How to make a decision: The analytic hierarchy process. European Journal of Operational Research 48 (1), 9–26
.
Sallwey
J.
, Schlick, R., Valverde, J. P. B., Junghanns, R., López, F. V. & Stefan, C.
2019
Suitability mapping for managed aquifer recharge: Development of web-tools
.
Water (Switzerland)
11
(
11
).
doi: 10.3390/w11112254
.
Sanchez-Iborra
R.
, Sanchez-Gomez, J., Ballesta-Viñas, J., Cano, M.-D. & Skarmeta, A. F.
2018
Performance evaluation of LoRa considering scenario conditions
.
Sensors (Switzerland)
18
(
3
).
doi: 10.3390/s18030772
.
Sendra
S.
, Parra, L., Jimenez, J. M., Garcia, L. & Lloret, J.
2022
LoRa-based network for water quality monitoring in coastal areas
.
Mobile Networks and Applications
28
(
1
),
65
81
.
doi: 10.1007/s11036-022-01994-8
.
Tanyanyiwa
C. T.
,
Armitage
N. P.
,
Okedi
J.
&
Schneuwly
R.
2023
Evaluating the hydrogeological and treatment performance of retrofitted infiltration basins for urban MAR in Cape Town, South Africa
. In
Novatech 2023: 11e Conférence internationale sur l'eau dans la ville
.
Van Wyk
E.
,
Van Tonder
G.
&
Vermeulen
D.
2011
Characteristics of local groundwater recharge cycles in South African semi-arid hard rock terrains: Rainfall-groundwater interaction
.
Water SA
38
(
5
),
747
754
.
doi: 10.4314/wsa.v38i5.14
.
Xu
Y.
&
Beekman
H. E.
2019
Review: Groundwater recharge estimation in arid and semi-arid Southern Africa
.
Hydrogeology Journal
27
(
3
),
929
943
.
doi: 10.1007/s10040-018-1898-8
.
Yeh
H. F.
, Cheng, Y.-S., Lin, H.-I. & Lee, C.-H.
2016
Mapping groundwater recharge potential zone using a GIS approach in Hualian River, Taiwan
.
Sustainable Environment Research
26
(
1
),
33
43
.
doi: 10.1016/j.serj.2015.09.005
.
Zhang
H.
,
Xu
Y.
&
Kanyerere
T.
2019
Site assessment for MAR through GIS and modeling in West Coast, South Africa
.
Water (Switzerland)
11
(
8
).
doi: 10.3390/w11081646
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).