Introduction

Water use in Africa is forecast to dramatically increase (Wada and Bierkens 2014), as more than half of global population growth by 2050 is projected to occur within the sub-Saharan region (UN 2019). Demand will be enhanced by increases in drinking water consumption, industry (Nieuwoudt et al. 2004), and irrigation for food production (Altchenko and Villholth 2015) in a continent where only 5% of arable land is currently irrigated (Siebert et al. 2010).

Groundwater is the largest store of freshwater in sub-Saharan Africa (MacDonald et al. 2012) and offers unrealised potential to contribute to meeting future water demands (Cobbing and Hiller 2019). Additionally, the exploitation of groundwater resources offers a myriad of benefits over surface water alternatives, particularly in this region of pronounced climatic variability (Braune and Xu 2010; Gaye and Tindimugaya 2019) of which 40% is classified as drylands (Cobbing and Hiller 2019).

Long-term groundwater hydrographs provide a direct indicator of groundwater abstraction sustainability from a quantitative perspective: allowing assessments of changes in storage, understanding recharge, and linkages to climate and land use change (Cuthbert et al. 2019). The rarity of such hydrographs has meant assessments of sub-Saharan African water security typically rely on datasets derived from large-scale models (e.g. Döll and Fiedler 2008), with limited validation using field observations (Sood and Smakhtin 2015), or large-scale data reviews (MacDonald et al. 2021). There are examples in the literature of multidecadal African hydrographs which have been used to estimate groundwater recharge and its association with rainfall intensity and climate variability (Kotchoni et al. 2019; Owor et al. 2009; Sibanda et al. 2009; Taylor et al. 2013); or changes in storage relating to land use changes (Favreau et al. 2009), abstraction and managed aquifer recharge (Murray et al. 2018). However, these studies have only examined hydrograph variation at a limited number of locations: from a single wellfield (Taylor et al. 2013) to a limited number of sites spread across an entire country (Kotchoni et al. 2019).

Broader studies include Cuthbert et al. (2019) who suggested three types of annual rainfall-recharge relationship relating to climatic zones from analysis of 14 multi-decadal hydrographs in nine countries across sub-Saharan Africa. Type 1 hydrographs were only present in humid and sub-humid climates where recharge was perennially consistent and uncorrelated to annual rainfall. Type 2 hydrographs showed increasing annual recharge with annual precipitation above governing precipitation thresholds, which tend to be greater as aridity increases, and occur across several climatic zones. Type 3 hydrographs occur in semi-arid to hyper-arid zones and have complex precipitation-recharge relationships, with recharge often only occurring in response to extreme precipitation events. In addition, Ascott et al. (2020) analysed 12 multi-decadal hydrographs from the Burkina Faso national monitoring network and distinguished them into sites showing: (1) long-term decline, (2) intra-annual variability, and (3) multi-decadal variability.

This paper presents an analysis of 21 long-term groundwater hydrographs from two adjacent semi-arid catchments in Limpopo Province, South Africa. This dataset is rare in the context of the previously mentioned studies as it allows us to investigate a larger number of long-term hydrographs distributed across a smaller area (35,100 km2) where climate variation is more limited. The area includes crystalline rock aquifers, which are widespread across the continent (Wright 1992). Furthermore, groundwater use is already extensively utilised here, with well-developed irrigated agriculture comprising an estimated 23% of cropped land in the Province (Cai et al. 2017). The paper aims to: (1) rationalise hydrograph behaviour into distinct types; (2) understand why these types exist spatially in the context of climate, hydrogeology, and land use and abstraction; (3) explore how these types relate to the recharge drivers of streamflow and climate, including the El Niño-Southern Oscillation (ENSO).

Methods

Study area

The study area comprises the Mogalakwena (19,200 km2) and Sand River (15,900 km2) catchments south of the east–west trending Soutpansberg mountain range within the Limpopo River Basin, Limpopo Province, northern South Africa (Fig. 1a,b). This semi-arid area is characterised by mean annual precipitation of 458 and 634 mm at Mara and Bela-Bela meteorological stations (1965–2010, Fig. S1 of the electronic supplementary material (ESM)—respectively, which is much lower than mean annual potential evaporation of 1,456 mm (DWS 2015). Rainfall falls predominantly (93%) in the summer wet season between October, the start of the hydrological year, and April (Fig. 1c). The main river channels in the catchments are usually perennial, but flows are limited in the dry season, being absent during periods of drought (see Fig. S2 of the ESM for monthly flows in the Mogalakwena River from 1965 to 2008). Most tributaries are ephemeral, flowing during the wet season or in response to event-driven intense rainfall (Holland and Witthüser 2011). Streamflow is primarily a result of recent rainfall, with baseflow a minor component as evidenced by limited or absent dry season flow. Elevation ranges between approximately 700 and 2,100 m above sea level.

Fig. 1
figure 1

a Location of study area within South Africa; b Detailed study area map illustrating geology, urban centres (black circles), South African Weather Service meteorological stations (black squares), river gauge (black triangle), and surface water flow direction (blue arrows). Bedrock geology simplified from the 1:1000000 map of South Africa (Council for Geoscience 2019); c Boxplot of monthly rainfall at Mara (1965–2010) showing median (dark blue), interquartile range (IQR) (blue), calculated minima and maxima (1.5 × IQR) (light blue), and outliers excluded

Annual rainfall is highly variable (coefficient of variation 25–35%, Fig. S1 of the ESM), with a consequential cycle of droughts and floods (Wetterhall et al. 2015). ENSO is considered the main cause of interannual variability in rainfall in southern Africa (Kolusu et al. 2019). Positive ENSO events (El Niño) and negative ENSO events (La Niña) have been associated with below-average and above-average summer rainfall, respectively, in South Africa (MacKellar et al. 2014; Reason et al. 2005). Nevertheless, any relationship between ENSO and rainfall is nonlinear and complicated by interactions with other modes of atmospheric variability (Fauchereau et al. 2009) such as the Indian Ocean Dipole (Gaughan et al. 2016).

The geology of the catchments is dominated by: crystalline basement rocks (gneiss/granite); the Northern Limb of the Bushveld Complex; the Soutpansberg and Waterberg Group; and the Karoo Supergroup. Basement rocks are of Archaean age and mainly comprise granite and gneiss. The Bushveld Complex intruded into the basement approximately 2 billion years ago and includes the Lebowa Granites and the mainly gabbroic succession of the Rustenburg Layered Suite. These crystalline rocks are all overlain by a regolith of between 10 and 50 m thick (Holland 2011). The Proterozoic Soutpansberg and Waterberg Group include red quartzitic sandstones and conglomerates with basaltic lavas forming the E–W trending Soutpansberg. The Karoo Supergroup consists of a mixture of sandstone, siltstones, and mudstones, as well as the Lebombo Basalt. There are frequent ENE–NE trending dyke swarms evident in aeromagnetic data (Stettler et al. 1989).

Younger alluvial deposits are found along the main river channels. Localised investigations have demonstrated 4–16 m of clay and sand at the Polokwane sewage works, and 10–14 m of red or sandy clay overlying 4–10 m of sand, gravel and pebbles at Mara within the Sand River Catchment (Vegter 2003). These deposits comprise local aquifers that are completely saturated beneath the rivers during streamflow. On the floodplain of the Nyl River, around 35 m of sediments have accumulated in the Nylsvlei wetland upstream of the confluence with the Mogalakwena River (McCarthy et al. 2011).

Mean transmissivity and mean borehole yields are similar between dominant lithology types in Limpopo Province with variations of 20–44 m2/day and 0.7–1.3 L/s, respectively (Holland 2011; Holland and Witthüser 2011). Basement aquifers are slightly higher yielding than elsewhere in southern Africa and there are localised high yield anomalies, notably air lift yields of 40 L/s are common around Mogwadi (Holland 2011). Water storage and movement is via the matrix and fractures in the weathered basement and Bushveld Complex; fractures in the Soutpansberge/Waterberg and Karoo; and only the matrix in the alluvium (du Toit 2003). Runoff is generally low due to the low topographic gradient and sandy soils (DWS 2015).

There has been intensive investment in pivot irrigation (Cai et al. 2017; Ebrahim et al. 2019), and it is estimated that irrigation for the commercial agricultural sector constitutes around 80% of water withdrawals in both catchments, which is almost exclusively supplied by groundwater (DWAF 2004). The town of Mogwadi and many rural communities are completely groundwater dependent; wellfields also partially supply Polokwane and Mokopane, in addition to some mining operations (DWAF 2004; DWS 2015).

Data and study period selection

Within the study catchments there are at least 21 monitoring boreholes with long-term records (>20 years) monitored by the Department of Water and Sanitation (DWS); the longest record extends from 1965 until 2013. Retrieved DWS data were converted to mean monthly time series from a mix of monthly dipped and hourly logged observations. Linear interpolation was then used to infill any gaps equal to or less than 3 months, given the long hydraulic memory typically observed within aquifers in drylands (Opie et al. 2020). Outstanding gaps were infilled using the hydrograph with the highest correlation coefficient, where overlapping data were available. The first derivative of the correlated hydrograph was used to infill because of discrepancies in the absolute values between sites. All infilling was visually inspected to confirm anomalous data were not introduced. Finally, a 20-year study period was selected where all records were complete following infilling: between November 1980 and November 2000 (Fig. 2). The entire dataset during this period was 94% original, 4% infilled by linear interpolation, and 2% infilled using a correlated hydrograph. There are limited supporting data available for the boreholes: depths are unknown, there are no lithological logs or completion details, site specific aquifer properties are unclear with an absence of pumping test data, and the authors are not aware of any attempts to characterise fractures.

Fig. 2
figure 2

Borehole hydrograph data availability and infilling. Dotted line illustrates the selected study period. Respective DWS names and metadata for the sites are provided in Table 1

Long-term daily rainfall data were retrieved from South African Weather Service (SAWS) climate stations at Mara (1949-present, 900 m above mean sea level (m asl)) and the more elevated Bela-Bela (1937-present, 1130 m asl). No infilling of data was conducted over the study period and any missing data were assumed to be zero. The Bela-Bela record contains 34 daily gaps, 33 of which occurred during the dry season; no rainfall on the single missing wet season day occurred at Mara. The Mara record includes 37 gaps, of which 32 days were in the dry season. The corresponding Bela-Bela data over the five missing wet season days were 0, 0, 1, 4.8, and 33 mm. Rainfall at both stations were summed to monthly data.

Monthly river flow data were retrieved from DWS for the Mogalakwena gauge (A6H009) aggregating across an area of 14,700 km2, almost the entire catchment, with data spanning from 1960 to 2008. During the study period, there were four gaps: three of which occurred towards the end of the dry season in either August or September, when flows are generally low or absent, and one during the drought of 1992, in February, when rainfall was only 24 mm at Mara (Fig. S2 of the ESM). These four gaps were filled by linear interpolation. There were no suitable river flow data covering the study period available for the Sand River Catchment.

The study period commences following successive years of consistently above mean rainfall (Fig. S1 of the ESM)—for example, between 1972 and 1980 at Mara, annual rainfall ranged between 513 and 734 mm, with a mean of 610 mm, >30% above the longer term mean. The study period includes major regional droughts including 1982 and 1992 (Trambauer et al. 2014). Exceptional rainfall also triggered major flood events in both 1996 and 2000 (Crimp and Mason 1999; Dyson and Van Heerden 2001). The 2000 event was the most extreme—for example, 1,388 mm of rainfall fell in Louis Trichardt during February 2000 (Dyson and Van Heerden 2001).

Supporting high spatial resolution (30 s of a longitude/latitude degree) distributed datasets of annual mean rainfall (1970–2000; Fick and Hijmans 2017) and aridity index (Trabucco and Zomer 2018) were retrieved. These climatological datasets were mapped across the study area and individual values extracted for each borehole (Table 1) to supplement the data from the two meteorological stations. Land use data were retrieved from the South African Department of Agriculture, Forestry and Fisheries (DAFF 2015). These data were simplified for visualisation into eight categories: natural vegetation, commercial agriculture (including pivot irrigation), commercial agriculture (irrigated without the use of pivots), subsistence agriculture, mining, town/village, open water, and wetland. All agricultural categories are classified into low, medium, or high intensity in the source data. The Multivariate ENSO Index Version 2 (MEI) (NOAA 2020) was used to investigate the teleconnection between groundwater levels and ENSO.

Table 1 Metadata for each site

Normalising hydrological data

To compare hydrographs across the study area, hydrographs were normalised to the standardised groundwater level index (SGI) following Bloomfield and Marchant (2013) for the study period. SGI is a nonparametric normalisation: data from a site are split into observations from each calendar month; data within each month are ordered, assigned a rank, and an inverse normal cumulative distribution function is applied; finally, the normalised monthly indices are merged to form a continuous SGI time series (Bloomfield et al. 2015). An SGI time series has a mean of zero and a standard deviation of 1. Note that this process inherently removes seasonality from the hydrographs.

Precipitation and river flow data were normalised following the same approach as the SGI, for consistency, and form a version of the standardised precipitation index (SPI; McKee et al. 1993) and standardised streamflow index (SSI; Vicente-Serrano et al. 2012), respectively. This process was undertaken for accumulation periods (aggregated rainfall) of k = 1–36 months to create SPI1–36 and SSI1–36. Note that as the accumulation period increases, month-to-month time-series variability is smoothed.

Statistical analysis

Cluster analysis was used to explore similarity between SGI time series and provide rationalisation of hydrographs into specific types. There are a wide range of clustering algorithms available (Haaf and Barthel 2018), which can be broadly split into hierarchical and non-hierarchical, although the decision over the approach to pursue is rather subjective (Webster and Oliver 1990). The hierarchical hclust function in R was selected for this study; utilising the minimum error sum of squares algorithm of Ward (1963). Initially each time series is allotted to its own cluster, before the algorithm proceeds iteratively joining the two most similar clusters, until a single cluster is formed. At each iteration, distances between clusters are recalculated by the Lance-Williams dissimilarity update formula applied according to Ward’s approach.

The autocorrelation of a time series demonstrates its correlation to a delayed copy of itself as a function of the delay, which was determined for each SGI time series using the Acf function in R. The autocorrelation of a hydrograph is representative of the hydrological memory in the system and was used to define the optimal number of clusters. Four clusters captured the main types of autocorrelation structure, with further clusters being single sites. A mean hydrograph (MH) of all SGI hydrographs within each cluster was produced to rationalise hydrological responses for further investigation.

SGI MHs were cross-correlated with both SPIs and SSI across accumulation periods of 1–36 months and lags of 0–6 months to examine differences in hydrological drivers between cluster types (Barker et al. 2016; Bloomfield and Marchant 2013). To investigate the relationship between ENSO anomalies—represented by MEI.v2 (MEI)—and SGIs, the MEI was averaged over October to April for each hydrological year and correlated by Spearman’s rank-order correlation with the change in SGI over the same period. All statistical analyses were performed in R version 3.4.0 (R Core Team 2017).

Results

Hydrograph classification

The SGI hydrographs were grouped into four distinct clusters (C1–C4; Fig. 3a,b), which each display a unique autocorrelation structure (Fig. 3c). The clusters are numbered sequentially according to the maximum length of their significant mean positive autocorrelation: C1 (28 months), C2 (36 months), C3 (56 months), C4 (70 months). Increasing the number of clusters beyond four produces clusters containing only single sites—for example, cluster 5 is BH 9 and cluster 6 is BH 5. These two sites do have a similar autocorrelation structure to the C2 mean, but they are the two sites with autocorrelation furthest from any of the four cluster means (Fig. 3c).

Fig. 3
figure 3

a Hierarchical cluster dendrogram separated into four clusters (C1–C4); b Mean SGI hydrograph for each cluster; c correlograms for each SGI time series split by hierarchical cluster for lags up to 80 months. The mean correlogram for each cluster is shown in black and the horizontal dotted lines denotes significance (p < 0.05)

Different cluster mean hydrographs (MHs) have characteristic responses to precipitation and length of recessions (Figs. 3b and 4). These observations explain the autocorrelation: increasing autocorrelation is associated with increasingly subdued responses to rainfall and longer periods of groundwater recession. MH C1 fluctuates through multiple periods of above and below normal SGI between 1980 and 1991 before retreating to its lowest levels between 1992 and 1995. During 1996, the SGI rises rapidly and remains above normal through to 2000 where SGIs increase to their maximum. MH C2 interannual fluctuations are similar but more limited than MH C1 between 1980 and 1996, with the time series in recession except during 1987 and 1991. During 1996, the SGI rises rapidly before recessing below normal in 1998 and increasing again in 2000. MH C3 is in recession in nearly all years with the notable exceptions of the extreme rainfall events of 1996 and 2000 occurring during negative ENSO periods (Fig. 4b,c). MH C4 is also in recession in nearly all years, with the only notable exception of 2000 when both SPI and SSI were most positive (Fig. 4c). Both C3 and C4 MHs are at their maxima at the beginning of the study period.

Fig. 4
figure 4

a All SGI time series grouped by cluster (C1–C4); b ENSO anomalies represented by the MEI; and c SPI24 at Mara and Bela (Bela-Bela), SSI24 at Moga (Mogalakwena), and SGI of the MH for each cluster

There is broad consistency of SGI time series within each cluster and with its respective MH (Fig. 4a). The most divergent time series within a cluster are BH 5 and BH 9 within C2, as previously distinguished by their autocorrelation (Fig. 3a,c). These sites have periods of below normal SGI pre-1988 when other time series are at least normal; do not have the same coherent negative SGI in 1992–1995; and are more negative in 1998–1999. Interestingly, BH 20 is the only site in C4 that responds to the 1996 rainfall event (Fig. 4a).

Spatial distribution of cluster types

In the Mogalakwena catchment, there appears a transition from predominantly C2 sites and a single C1 site, to C3 sites, and then C4 sites downstream in the catchment (Fig. 5). The transition to C3 sites occurs abruptly at Mokopane, for example site 13 (C2, upstream of town) is within 5 km of site 14 (C3, within the town). In the Sand River catchment, there is a similar transition downstream along the catchment from site 9 (C2), to site 16 (C3) and finally, site 21 (C4). However, around Polokwane, C1, C3 and C4 sites are all located within 7.5 km, with a transition from C1 to C4 within only 4 km to the south-west of the town.

Fig. 5
figure 5

Spatial distribution of clustered borehole sites in the study area. Geology simplified from the 1:1000000 map of South Africa (Council for Geoscience 2019). Blue arrows indicate surface-water flow direction

Relationships between mean hydrographs and rainfall/river flow

Mean hydrographs C1, C2, C3, and C4 are strongly, moderately, weakly and very weakly correlated, respectively, with SPI for a given accumulation period (Fig. 6; Table 2). For MH C1 the correlation increases with SPI accumulation periods for SPI at Mara up to the maximum of 36 months, and increases with accumulation periods for SPI at Bela-Bela before peaking at 25 months. The correlation coefficient (r) is 0.63 and 0.52 for Bela-Bela and Mara, respectively, by k = 12. For MHs C2 and C3, the correlation with SPI increases with accumulation period up to around 5–10 months, before declining, then increasing again and peaking at 21–35 months. For MH C4, the most correlated accumulation periods are much shorter (2–3 months), although given the weakness of the relationship, it is challenging to contrast this value with the other cluster accumulation periods.

Fig. 6
figure 6

Correlation heat maps for each of the four mean hydrographs (MHs) for a range of accumulation periods and lags of SPI Bela-Bela, SPI Mara, and SSI Mogalakwena. Black dots mark the highest Pearson’s correlation coefficient (r)

Table 2 Maximum correlation (r) between SGI and SPI for a range of accumulation periods (k) and lags up to 6 months

Increasing lag in SPI reduces correlations for C1. The optimal lags for C2 and C3 are a month or two months, respectively, with further lags reducing the correlation. The optimal lags are more inconsistent for C4 at 2–6 months. Maximum correlation coefficients are reasonably consistent for both rainfall sites (Table 2).

MHs C2, C3, and C4 are more strongly correlated with SSI than SPI, with C2 and C4 being the most and least correlated, respectively (Fig. 6; Table 2). MH C1 is similarly correlated with both SSI and SPI. The most correlated accumulation periods and lags for each cluster are similar for both SSI and SPI (Table 2). Therefore, there is evidence that groundwater storage variations are more statistically related to river flow dynamics than rainfall dynamics for C2, C3, and C4, but not C1.

Relationships between mean hydrographs and ENSO

There is a negative relationship between mean wet season MEI and the change in SGI for all clusters (Fig. 7). The greatest increases in SGI generally occurred during the negative ENSO (La Niña) years of 1996 and 2000 and there was no increase in SGI during the most positive ENSO (El Niño) years (MEI > 1.5). Negative MEI does not necessarily relate to rises in the water table over the wet season—for example, in 1989, which had the strongest negative anomaly, rainfall was 316 mm in Mara and groundwater levels did not rise over the wet season in any cluster. There is a more consistent tendency for groundwater levels not to rise during strongly positive ENSO years.

Fig. 7
figure 7

Wet season mean Multivariate ENSO Index (MEI. v2) versus change in SGI for all clusters highlighting two wet season periods of anomalously high rainfall. Correlation is assessed by Spearman’s Rank, n = 20, and * denotes significance at p < 0.05

The nonparametric correlation coefficients are most negative and only significant (p < 0.05) for C3 and C4 (Fig. 7). Wet season increases in SGI at C3 and C4 are only evident when MEI < 0.4. C1 consistently rises across a much broader range of MEI values compared to the other MHs. Additionally, C1 does not respond as appreciably as the other MHs to the 2000 event, with the magnitude of the response being comparable to several other years. The 1996 and 2000 wet season increases in SGI are greatest at C2, but there is no clear relationship with ENSO with those years omitted.

Discussion

Factors controlling cluster types

Rainfall and aridity

Rainfall and aridity are a key control on groundwater hydrographs and recharge across the continent (Cuthbert et al. 2019; MacDonald et al. 2021). There is some evidence for a transition from either C1 or C2 to C3 or C4 down both catchments. Any downstream transition could relate to rainfall variability, particularly given the nonlinearity observed between rainfall and recharge in the continent, whereby small increases in annual rainfall in excess of a given threshold can lead to substantial rises in recharge (Taylor et al. 2013). However, there is minimal evidence for lower annual rainfall at C3 and C4 sites compared with C1 and C2 sites, although lowest annual rainfall does occur at site 19 (C4) (Fig. 8a; Fig. S3 of the ESM). Furthermore, the aridity index displays no appreciable variation between clusters despite varying between 0.24 and 0.38 at the sites (Fig. 8b; Fig. S4 of the ESM). It should also be reiterated that changes in hydrograph type also occur over short distances where any climate variation is negligible. Consequently, there is limited evidence for rainfall and aridity differentiating between the hydrograph clusters, and other local drivers are important.

Fig. 8
figure 8

Hydrogeological similarities between hydrograph clusters: a mean annual rainfall (Fick and Hijmans 2017), b aridity index (Trabucco and Zomer 2018), c distance to surface water, d groundwater level in metres below datum (m bd) during study period. Boxes illustrate the 25 and 75th percentiles, dissected by the median; whiskers indicate the 10th and 90th percentiles; and all outliers are shown as dots

Hydrogeology

In semi-arid climates, groundwater recharge is often postulated to occur beneath riverbeds (Cuthbert et al. 2016; Meredith et al. 2015; Scanlon et al. 2006), therefore, one may expect sites further from rivers to be less dynamic, with longer autocorrelation. However, Van Wyk (2010) demonstrated hydrographs did not adhere to this suggested model in crystalline rock settings in South Africa with more dynamic hydrographs observed away from rivers in more upland settings where soils and regoliths were thinner. In terms of the hydrographs from this study, two of the three sites furthest from a river are in C4 (Fig. 8c). Moreover, site 19, the only C4 type near a river (ca. 400 m), lies close to only a small tributary of the Mogalakwena, which is likely to flow less frequently than the main channel with a comparably smaller catchment area. There is also evidence for a transition from C1 types (sites 4, 2, 1), to C3 type (site 15) and then C4 type (site 20) as you move approximately perpendicular away from the river running through Polokwane (Fig. 5). Nevertheless, there are also inconsistencies in these relationships—for example, site 18 (C3) also lies close to the river in Polokwane, the second furthest site from the river is a C1 type (site 3), and there is no clear variation in distance from nearest surface water between C1, C2, and C3 (Fig. 8c).

A deeper unsaturated zone can smooth groundwater level changes and increase the autocorrelation structure; however, groundwater levels were similar across C1, C2, and C3 types, with the deepest mean groundwater levels across all sites actually at a C2 type (site 8) 41.5 m below datum (bd). There is evidence for generally deeper groundwater levels at C4 types (Fig. 8d), partially a result of the declining trend in levels over the study period (Fig. S5 of the ESM).

There is no relationship between cluster type and geology (Fig. 5). All four cluster types occur in the Archaean gneiss/granite (Fig. 5), including C1, C3 and C4 types in close proximity around Polokwane. Types C2, C3, and C4 also all occur in the Rustenburg Layered Suite. It is unclear whether boreholes are screened within the overburden or bedrock, which could behave differently, although it is common practice to screen boreholes across water strikes in both the regolith and fractured rock. Furthermore, there is uncertainty over variations in hydraulic conductivity and storage between sites, as well as the potential of structural controls, principally dykes and faults, on the hydrogeology, which could restrict groundwater flow (Ebrahim et al. 2019). Therefore, there is potential to have compartmentalised aquifer units (Abiye et al. 2020) that are less dynamic within the same broad geological classification.

Abstraction and land use

There are concerns in both catchments regarding the high abstraction of groundwater (Abiye et al. 2020; Busari 2008; DWAF 2004; DWS 2015; Masiyandima et al. 2002). Abstraction has the potential to introduce downward trends in hydrographs (Oiro et al. 2020) and increase autocorrelation (Wendt et al. 2020). Pivot irrigation is well developed in the study area (Cai et al. 2017; Ebrahim et al. 2019), particularly north of Polokwane in the Sand Catchment (Fig. 9), which is mainly reliant on groundwater. Specifically, around Mogwadi, near site 21 (C4), it is reported that pivot irrigation for commercial agriculture has reduced groundwater levels by 50 m between the 1970s and 2000 (Fallon et al. 2019; Masiyandima et al. 2002) and satellite imagery confirms extensive pivot irrigation encircling the site (Fig. S6a of the ESM). More localised pivot irrigation is also proximal to sites 18 and 16 (C3; Fig. 9). There is a mixture of pivot irrigation and medium- to high-intensive commercial agriculture, presumably irrigated by other means, around site 19 (C4) (Fig. 9; Fig. S6b of the ESM). Some C1 and C2 sites (e.g. 3, 5, 12) are located in areas of commercial agriculture, but this is predominately classed as low-intensive. Therefore, sites in areas of pivot irrigation and medium- to high-intensive commercial agriculture appear to be in C3 and C4.

Fig. 9
figure 9

Simplified land use classification of the study area (DAFF 2015) displaying all sites coloured by cluster type. MAR on the inset figure corresponds to the Polokwane Managed Aquifer Recharge scheme

Groundwater in the Polokwane Municipality is some of the most highly utilised in the study area and a previous water balance by the DWS stated it was >300% overexploited in 2010 (DWS 2016). The municipal water supply to the town has been mainly imported from neighbouring dammed surface water catchments since 1958 because groundwater was unable to sustain the growing demand (Vegter 2003). Groundwater still provided a subordinate component that varied between 2.17 and 4.01 million m3/year during the study period (Vegter 2003). Sites 15 (C3) and 20 (C4) are both located in the Sterkloop Wellfield, whilst sites 18 (C3) and 4 (C1) are just outside the Sand River Wellfield. Sites 1 (C1) and 2 (C1) lie within the urban conurbation (Fig. 9).

There is mining and smelting activity around Polokwane, with site 20 being within 2.5 km of this land use (Fig. 9), although the extent of groundwater use by these industries is unknown. The town also operates a managed aquifer recharge (MAR) scheme, but this is located north of Polokwane (Fig. 9) and produces an undulating water table in boreholes only downstream of the scheme (Murray and Tredoux 2002). The association of types C3 and C4 with the Polokwane wellfields, as opposed to the C1 sites located within only 2 km of these sites and outside the wellfields, indicates municipal abstraction may be an important control on groundwater hydrographs locally.

Groundwater resources are also under pressure around Mokopane, where there is a shift from C2 to C3 sites in the Mogalakwena Catchment. Notably, groundwater levels have been declining around the Mokopane wellfield (Busari 2008), although it is unclear if it is the wellfield itself which has led to the declining water table. Mining is a large employer in the town, including the largest open pit platinum mine in the world north of the town. The extensive workings are large water consumers, although the mine is downstream of the C3 sites. Furthermore, it should be noted that given the water pressures in both catchments, several mines have switched to utilising grey water from both the Mokopane and Polokwane wastewater treatment works.

Evidence for recharge processes

Stronger relationships between groundwater levels and river flow than groundwater levels and rainfall for C2, C3 and C4 may suggest that groundwater recharge is dominated by leakage from surface water at those sites. The Polokwane MAR scheme makes use of the interconnected Sand River, alluvial aquifer, and bedrock aquifers where treated municipal wastewater discharges to the river to recharge the bedrock aquifer that is then tapped for public water supply and agriculture (Murray and Tredoux 2002). However, Walker et al. (2018) concluded the 2-m-thick alluvial deposits beneath a reach of the Molotsi sand river, a bordering catchment to the Sand River with similar geology, were hydraulically unconnected to the underlying fractured basement. Furthermore, along the Nyl River, an upstream tributary of the Mogalakwena (Fig. 5), Tooth et al. (2002) reported that flooding had effectively sealed the floodplain with thin clay layers which limit groundwater recharge. They did note, though, that some recharge does occur at the margins of the floodplain where the superficial deposits are coarser.

Therefore, within the literature, there is a lack of consistent local studies that support widespread surface water and groundwater connectivity, let alone the suggested dominance of focussed recharge from surface waters. An equally plausible explanation of the relationship between SGIs and SSI is that both river flow and groundwater recharge are an expression of soil moisture excess in the catchments and therefore related only indirectly. The SGI and SSI relationship could also theoretically indicate a groundwater-fed river network, but this is unlikely in the local conditions where baseflow is a minor component of streamflow.

Importance of teleconnections and extreme events

The relationship to ENSO is strongest for MHs C3 and C4. It is suggested here that the stronger relationship is a combination of two factors. Firstly, frequent recharge signals that are unrelated to a particularly strong/weak ENSO are masked by abstraction. Secondly, the lowering of the water table over several years by abstraction may increase available subsurface storage resulting in a more notable rise in the water table when extreme ENSO-driven rainfall events occur. Evidence for this latter factor can be drawn from the contrasting responses to the extreme rainfall event of 2000. SGIs in 2000 for C3 and C4 MHs rose substantially but did not achieve their maxima following years of predominantly groundwater recession. On the other hand, the magnitude of the rise in the C1 MH was comparable to multiple other years, despite the extreme nature of the rainfall, and the SGI peaked. Examination of the raw groundwater level data from three sites (1, 2, and 4) of the four C1 sites demonstrates that water levels in the 2000 wet season were around their peaks across their entire records (up to 40 years) and relatively shallow (6.9, 2.9, 2.6 m below datum; Fig. S5 of the ESM). This evidence may indicate the response to the 2000 event at some C1 sites may had been storage limited. It is also possible that at such shallow water depths in these crystalline rock settings, there are layers within the soils of very high permeability that could initiate rapid lateral flows and restrict further rises in levels (Bonsor et al. 2014). Vertical variability in specific yield could be significant in terms of recharge response and the mean water levels prior to the event differed for C1 to C4 at 8.2, 16.1, 12.2 and 29.3 m bd, respectively. However, there is no information on weathering thickness or lithology from the individual boreholes to assess vertical contrasts in specific yield.

Long-term groundwater hydrographs elsewhere in the continent, that tend to be located in wellfields, display similar behaviour to C3 and C4, such as the 60-year hydrographs from the Dodoma wellfield, Tanzania (Taylor et al. 2013). Such datasets from areas of high abstraction may overemphasise the importance of recharge from extreme events, as well as teleconnections, compared to more natural settings. Indeed, recent modelling work by Seddon (2019) for the Dodoma wellfield does show that accounting for the influence of abstraction on groundwater levels hydrographs does reduce the number of years reporting zero recharge. Nevertheless, their analysis still demonstrates heavy rainfalls contribute disproportionately to recharge and extreme events are undoubtedly invaluable from a water security perspective in Africa.

The major recharge events for all hydrographs types are associated with extreme rainfall occurring during negative ENSO years. The 1996 event is highly pronounced within C1, C2, and C3 MHs, and follows the early 1990s drought where groundwater levels had fallen widely. It is unclear why two of the C4 types (sites 19 and 21) do not respond to this event, but they are furthest downstream and distant from the main river channels. The more substantial 2000 event triggered strong responses in all types.

There is consistent evidence across all cluster types that groundwater levels do not rise during extreme positive ENSO years. This is supported by an analysis of ENSO and rainfall in Limpopo Province, which demonstrated a robust positive relationship between ENSO anomalies (Niño 3.4 sea surface temperature) and dry spell frequency between 1979 and 2002 (Reason et al. 2005). ENSO anomalies could also be used to predict groundwater recharge, given that model predictions of rainfall in southern Africa can be improved by the inclusion of ENSO (Kolusu et al. 2019; Landman and Beraki 2012).

Conclusions

Twenty-one groundwater hydrographs in two adjacent semi-arid (0.24–0.38 aridity index) catchments in South Africa, from 1980 to 2000 when records overlapped, were classified into four cluster types (C1–C4). Hydrograph types C1 through to C4 show a transition of increasing autocorrelation and increasingly subdued rainfall responses. C1 type is strongly related to rainfall, has the optimal cross-correlation with the standardised precipitation index (SPI), fluctuates on a typically annual basis, and has significant (p < 0.05) positive autocorrelation of 28 months. C4 type is minimally related to rainfall and characterised by multiannual recessions with a significant positive autocorrelation of 70 months with levels rising only in response to an extreme rainfall event in 2000. C2 and C3 types are intermediates between the C1 and C4 extremes.

C1 type is similarly related to both the SPI and the standardised streamflow index (SSI), whereas C2–C4 are more strongly associated with SSI. These correlations may suggest C2–C4 are more dependent on focussed recharge from riverbeds, though further investigation using environmental tracers is required to confirm this assertion as river flow and diffuse recharge are both driven by an excess of soil moisture.

There is a tendency for substantial or minimal recharge to occur during extreme negative and positive ENSO years, respectively, thus ENSO anomalies could be useful to predict groundwater recharge. Large recharge events occur in 1996 and 2000 during negative ENSO years (La Niña), though large recharge events do not always occur during such years. Declines in the water table are associated with extreme positive ENSO years (MEI > 1.5). Only SGI changes in C3 and C4 types are significantly correlated with wet season El Niño–Southern Oscillation (ENSO) anomalies across the study period.

The range of groundwater hydrograph types, sometimes within only a few kilometres, cannot be attributed to spatial variability in either climate or hydrogeological settings such as distance to surface water, depth to water table, or available geological information. C3 and C4 sites appear to be associated with areas of high groundwater abstraction such as municipal wellfields and intensive commercial agriculture. It is considered that high levels of abstraction near C3 and C4 sites mask frequent small-scale recharge events observed at C1 and C2 sites. Lowering of water levels by abstraction may also increase available storage resulting in greater capture of recharge from extreme events and/or produce contrasting hydrograph responses to recharge due to vertical variability in aquifer properties. Consequently, extreme events associated with positive ENSO years are most notable in the C3 and C4 time series.

Abstraction can bias interpretations of groundwater hydrographs concerning: the regularity of recharge, the relative importance of extreme recharge events, the strength of the relationship with the potential recharge drivers of rainfall and streamflow, and the significance of teleconnections. Therefore, care should be taken when analysing groundwater level data from areas of high abstraction such as within municipal wellfields or near intensively irrigated agriculture.