Hydroclimatic Extremes in the Limpopo River Basin, South Africa, under Changing Climate

: This research study evaluated the projected future climate and anticipated impacts on water-linked sectors on the transboundary Limpopo River Basin (LRB) with a focus on South Africa. Streamﬂow was simulated from two CORDEX-Africa regional climate models (RCMs) forced by the 5th phase of the Coupled Model Inter-Comparison Project (CMIP5) Global Climate Models (GCMs), namely, the CanESM2m and IPSL-CM5A-MR climate models. Three climate projection time intervals were considered spanning from 2006 to 2099 and delineated as follows: current climatology (2006–2035), near future (2036–2065) and end of century future projection (2070–2099). Statistical metrics derived from the projected streamﬂow were used to assess the impacts of the changing climate on water-linked sectors. These metrics included streamﬂow trends, low and high ﬂow quantile probabilities, the Standardized Streamﬂow Index (SSI) trends and the proportion (%) of dry and wet years, as well as drought monitoring indicators. Based on the Mann-Kendall (MK) trend test, the LRB is projected to experience reduced streamﬂow in both the near and the distant future. The basin is projected to experience frequent dry and wet conditions that can translate to drought and ﬂash ﬂoods, respectively. In particular, a high proportion of dry and a few incidences of wet years are expected in the basin in the future. In general, the ﬁndings of this research study will inform and enhance climate change adaptation and mitigation policy decisions and implementation thereof, to sustain the livelihoods of vulnerable communities.


Introduction
Adverse repercussions attributed to climate change and variability on socio-economic sectors, such as water resources, energy, agriculture, health and ecosystem services, are gradually becoming a major concern to various communities, locally, nationally, regionally and even globally [1]. Consequently, there are growing calls for government institutions to revise and enhance the existing mitigation and adaptation planning strategies, to prepare the most vulnerable communities exposed to stressors of climate change [1]. Hydrological extremes such as floods and drought cause more increased havoc yearly-these impacts include loss of life and agricultural produce as well as the destruction of infrastructure at various spatial and temporal scales [2]. Such extremes are also attributed to the vulnerability of groundwater resource systems as well as the sensitivity of such systems to natural disasters through impacts on the changes in climatic variables such as precipitation, evapotranspiration and streamflow [3,4]. Groundwater systems form part of the hydrological cycle and are too affected by climate change and variability, coupled with socio-economic activities [5,6]. For instance, variations in rainfall patterns coupled with land-use change effects can intensify soil erosion affecting renewable groundwater resources. Such variations can have an impact on groundwater-surface systems, including wetlands and other groundwater dependent ecosystems. Consequently, groundwater systems are essential for effective monitoring, planning and management of water resources, under the changing climate [3,4].
Several research studies suggest that the observed historical trends in the frequency of heavy to extreme hydrological events are associated with extreme precipitation events [7], increases in global average air, ocean and near-surface temperatures, widespread melting of snow and ice, as well as rising global average sea level [1, [8][9][10][11][12][13]. Furthermore, future climate projections have indicated a decrease in groundwater recharge by the 2050s, considered as a limiting factor for the sustainable use of groundwater [14], and increase in most climatic variables (e.g., global mean temperature [15][16][17], heat stress and frequent and prolonged heat waves [18,19] and highly variable precipitation [17]) in most parts of Africa. Such changes in projected climate variables are likely to exacerbate the frequent occurrences of hydrological extremes thereby elevating the strain on the already stressed water resources for agriculture, industry and households, thus resulting in a significant impact on local populations, as well as their livelihoods [20].
The unprecedented upsurge of climate change impacts calls for continuous monitoring and appraisal of metrics of climate change across all climate-sensitive sectors. This includes the monitoring of climatological variables such as precipitation, evaporation and air temperature, among others, as well as hydrometric variables such as streamflow (discharge). There are various ways in which climatological and hydrometric variables can be monitored. These include the use of integrated models to model interactions between natural and human processes [3], the assessment of trends in hydrological statistical parameters (e.g., annual mean, median, minimum and maximum) of the climate variables [21][22][23][24], General Extreme Value (GEV) analysis to determine the probability distribution and return periods of extreme events [25,26] and drought monitoring based on drought indices and drought monitoring indicators [27][28][29][30][31]. Research studies on monitoring of hydroclimatic variables are an essential supportive tool for policy makers in disaster preparedness, particularly in the most vulnerable communities, who are anticipated to always respond in earnest and continuously propose climate-related policies that support societal well-being, protection of infrastructure and all sectors of the economy.
In general, research studies in the monitoring of hydroclimatic variables and inherent extremes at catchment, national regional and global levels are well documented in the literature. However, there have been insufficient studies of future projections of hydroclimatic extremes, particularly in South Africa, more so at the catchment level, despite the country's vulnerability to climate change. Given that the socio-economic impacts of weather and climate extremes are likely to increase as the extremes become more frequent [20], simulations of reliable future climate changes will aid our understanding of the projected climate and the inherent impacts on water-linked sectors. In this regard, Water 2020, 12, 3299 3 of 20 this study aims to evaluate the projected future climate and anticipated impacts on water-linked sectors in the Limpopo River Basin (LRB) with a focus on South Africa. To achieve this, streamflow simulated from the Coordinated Regional Downscaling Experiment (CORDEX) model projections are assessed at three different time intervals. Statistical metrics such as trends in projected streamflow, the Standardized Streamflow Index (SSI), the proportion of dry and wet years and drought monitoring indicators are used to assess the impacts of future climate on water-linked sectors in the LRB.

Study Area
The LRB, shown in Figure 1, is the fourth largest transboundary river basin in Southern Africa spanning over 416,296 km 2 . The basin spreads over four countries and flows from the border between South Africa and Botswana, through South Africa and Zimbabwe and into Mozambique where it enters the Indian Ocean at Zongoene [32]. There are about 14 million people who are supported by the basin and slightly over 50% of the population is rural, however, the rate of urbanization in the basin has also increased over time [33]. A significant proportion of the rural population is involved in rainfed agriculture to sustain their livelihoods [34]. The basin has a wide range of wildlife and biodiversity which sustains tourism in the region. Currently, almost all the water in the upper LRB is allocated to different water users in the social and economic sectors such as agriculture, industry, power stations and municipalities [35].
Water 2020, 12, x FOR PEER REVIEW 4 of 21 Figure 1. The Limpopo River Basin and its attributes. The dots represent the grid centers of the Coordinated Regional Downscaling Experiment (CORDEX) models at 50 km spatial resolution.

Observational Datasets
Monthly observed rainfall and minimum and maximum temperature data were obtained from the Hadley Center Climate Research Unit (CRU) TS2.0 monthly 50 km resolution observed dataset [37]. Daily streamflow observational data for the study site was obtained from the Department of Water and Sanitation of South Africa. Four downstream streamflow gauge stations were selected for the study. The selected stations have been operational for over three decades; hence, they provide reliable and continuous datasets. Two of the downstream streamflow stations were selected for the period, spanning from 1976 to 2005. The averaged streamflow data from the selected stations were utilized to bias correct the simulated streamflow across the four analyzed periods. The climate in the LRB is generally subtropical dry savannah with highly variable rainfall across the basin, ranging between 200 mm/year in the western drier parts of the basin and 600 mm/year in the eastern parts close to the Indian Ocean. The south middle parts of the catchment receive above 1500 mm/year, and more than 90% of the rain is received between October and April [32]. The LRB is also prone to extreme weather events such as floods and droughts, highlighting the need for disaster preparedness and efficient early warning systems [33]. The climate is influenced by, among other things, Water 2020, 12, 3299 4 of 20 the continental tropical, equatorial convergence zone and moist, marine, subtropical eastern and marine western Mediterranean air masses [36]. On the one hand, the lower Limpopo River for example from Chokwe to the Guija district to the Mozambican coast at Xai Xai are low lying areas of the floodplain that have rich soils but are more susceptible to flooding and loss of agricultural produce and other assets as a result of high rainfall and tropical cyclones in the area as well as the upper catchments [33]. Consequently, despite the rich soils, poverty, loss of biodiversity, huge sediment load and food and water insecurity are some of the persistent challenges faced by communities in these parts of the basin.

Observational Datasets
Monthly observed rainfall and minimum and maximum temperature data were obtained from the Hadley Center Climate Research Unit (CRU) TS2.0 monthly 50 km resolution observed dataset [37]. Daily streamflow observational data for the study site was obtained from the Department of Water and Sanitation of South Africa. Four downstream streamflow gauge stations were selected for the study. The selected stations have been operational for over three decades; hence, they provide reliable and continuous datasets. Two of the downstream streamflow stations were selected for the period, spanning from 1976 to 2005. The averaged streamflow data from the selected stations were utilized to bias correct the simulated streamflow across the four analyzed periods.

CORDEX Models and Sensitivity Analysis
The present research study used the CORDEX-Africa daily climate simulation data at a spatial resolution of 50 km to force an Agricultural Catchment Research Unit (ACRU) hydrological model. The CORDEX-Africa simulations were obtained from the CORDEX Rossby Centre regional model (RCA4), e.g., the Regional Climate Models (RCMs) downscaled using a dynamic downscaling approach. In this technique, climate data from Global Climate Models (GCMs) are used as the main inputs that drive the RCMs. In the current analysis, the CORDEX climate simulation data from the RCA4 RCMs were driven by nine GCMs of the 5th phase of the Coupled Model Inter-Comparison Project (CMIP5) [38]  The nine GCMs were first assessed through sensitivity analysis. The purpose of sensitivity analysis was to select the best model, that is, the model that was able to simulate the climate conditions in the LRB. Finding a model that reproduces the general features of the catchment's climate correctly is always a challenge due to various factors, which include model parametrization, the spatial resolution of the models and topographic features, including land-surface heterogeneities. Given that rainfall is the main driver of streamflow, RCMs' simulated rainfall was evaluated against the observations. Models that performed well in simulating spatial distribution of annual rainfall were chosen to drive the hydrology model. As shown in Figure 2, the spatial distribution of annual rainfall totals is well captured by CanESM2 and IPSL-CM5A-MR RCMs. In the basin, many of areas in the western part receive a relatively low annual rainfall (about 500-600 mm/year), and rainfall increases towards the eastern part with some areas reaching a maximum of 1000-1500 mm/year. The two models were therefore selected and utilized to assess the historical and projected streamflow in the LRB and the inherent impacts on water-linked sectors. The ensemble model calculated from six GCMs was also included in the analysis. Results from the two models and the ensemble were compared with the aim of assessing the importance of using a correct model, particularly at a regional level, where the resolution of the models is always problematic.
analysis was to select the best model, that is, the model that was able to simulate the climate conditions in the LRB. Finding a model that reproduces the general features of the catchment's climate correctly is always a challenge due to various factors, which include model parametrization, the spatial resolution of the models and topographic features, including land-surface heterogeneities. Given that rainfall is the main driver of streamflow, RCMs' simulated rainfall was evaluated against the observations. Models that performed well in simulating spatial distribution of annual rainfall were chosen to drive the hydrology model. As shown in Figure 2, the spatial distribution of annual rainfall totals is well captured by CanESM2 and IPSL-CM5A-MR RCMs. In the basin, many of areas in the western part receive a relatively low annual rainfall (about 500-600 mm/year), and rainfall increases towards the eastern part with some areas reaching a maximum of 1000-1500 mm/year. The two models were therefore selected and utilized to assess the historical and projected streamflow in the LRB and the inherent impacts on water-linked sectors. The ensemble model calculated from six GCMs was also included in the analysis. Results from the two models and the ensemble were compared with the aim of assessing the importance of using a correct model, particularly at a regional level, where the resolution of the models is always problematic.

Hydrological Model
The ACRU model, developed at the former School of Bioresources Engineering and Environmental Hydrology of the University of KwaZulu-Natal, South Africa [39], was used in this study to simulate streamflow across the selected CORDEX models and time intervals. This

Hydrological Model
The ACRU model, developed at the former School of Bioresources Engineering and Environmental Hydrology of the University of KwaZulu-Natal, South Africa [39], was used in this study to simulate streamflow across the selected CORDEX models and time intervals. This hydrological model is a multi-purpose and integrated multi-soil-layered agro-hydrological modelling system, which runs on a daily time-step. ACRU is configured to be hydrologically sensitive to catchment land uses and the corresponding changes, hence, it has the ability to simulate catchment hydrological responses relating to land management, including the impacts of land-cover/land-use changes on water flows [39,40]. The model requirements include inputs of known and measurable spatio-temporal variable factors that characterize the watershed. The minimum catchment information required to operate ACRU is the daily precipitation as well as daily maximum and minimum temperature datasets. Output information derived from ACRU include daily streamflow, separated into stormflow and baseflow, peak discharge, reservoir status, recharge to groundwater, irrigation water supply and demand, crop yields and sediment yield. The ACRU model has been widely used in various hydrological studies across Southern Africa [41][42][43] as well as in other parts globally, including Canada [44], New Zealand [45] and Ghana [46], among others. For a detailed description and theoretical background of the ACRU model, the reader is referred to [39]. In this study, ACRU was used in a lumped mode to simulate streamflow for four different time intervals, namely, the historical (1976-2005), present (2006-2035), near future (2036-2065) and distant future (2070-2099).

Streamflow Simulations
In this study, the impacts of climate change in the LRB were based on the analysis of streamflow simulated from CORDEX-Africa, forced by two CMIP5 GCMs under the RCP4.5 and RCP8.5 scenarios. In this regard, changes in the streamflow were analyzed based on 30-year future projection periods, Water 2020, 12, 3299 6 of 20 arbitrarily delineated in a specific time scale, from near (present), mid-century and distant future periods. Changes in the streamflow under the classified periods were compared with a 30-year reference (historical) period [47,48], which in this study was set from 1976 to 2005. The assessment of the impacts of climate change on water resources and water-linked sectors in the LRB was based on the evaluation of streamflow trends of the annual mean, minimum and maximum hydrological statistical parameters. The low-, mid-and high-flow (quantiles) characteristics were evaluated given that they relate to dry, normal and wet conditions, respectively. In addition, the trends in the drought duration and severity, calculated from the SSI, were analyzed across different periods.

Mann-Kendall Trend
The Mann-Kendall (MK) test, also considered as the original form of the MK test [49,50], was applied to the streamflow and SSI time series across the study periods to examine and detect if there is a monotonic increasing or decreasing linear or non-linear trend of the variables over time. The MK test is non-parametric in nature and is considered as a suitable method to detect shifts and trends of climatic variables [51]. Consequently, the MK trend has been widely used to detect trends of climatic variables in meteorology, hydrology and climatology research fields [52][53][54]. The test is robust against non-normal distribution datasets and it is suitable for data with missing values. Other modified MK tests reported in the literature, e.g., the modified MK test with trend-free pre-whitening [55], the modified MK test that considers the presence of autocorrelation structures in a time series [56] and the modified MK test that addresses the effects of the presence of trends and persistence, cover both the short-and long-term persistence in the trends [57]. In this study, statistically significant trends were computed for the 95% confidence level, thus, the null hypothesis (stating that there is no significant trend in the assessed time series) was rejected for every p-value less than 0.05.

Extreme Value Analysis of Streamflow
The GEV distribution, characterized by three parameters, namely, the location, scale and the shape, was utilized in this study to assess the features of simulated streamflow across the CORDEX models, RCPs and timescales. Daily simulated streamflow was aggregated into yearly maximum time series for each RCP scenario and timescale. The yearly maximum streamflow time series was fitted to the GEV distribution, resulting in an estimation of the location, scale and shape parameters that are suitable to explain the probability distribution of streamflow extremes. The estimated parameters were assessed based on the three classes of the GEV family of distribution, namely, the Gumbel, Frechet and Weibull (also referred to as type I, type II and type III extreme value distributions, respectively) [26]. The differences between these families of distributions depend on the value of the shape parameter. In particular, the GEV distribution corresponds to type I, type II and type III extreme value distributions, when a shape parameter is equal to zero, greater than zero or less than zero, respectively. Consequently, the shape parameter was used in this study to identify the extreme value distribution that best fits/describes streamflow extremes in the LRB across all the projected timescales.

Standardized Streamflow Index
This study adopted the principle of the Standardized Precipitation Index (SPI) [58] for the calculation of the SSI at two selected accumulation periods, with particular interest to the periods relating to hydrological drought. In this regard, simulated streamflow was aggregated into monthly streamflow data over 6-and 12-month accumulation periods. The gamma probability distribution (mostly recommended for SPI calculation) was utilized to fit the aggregated monthly streamflow time series [59]. The SSI was calculated across the selected models and the study periods, and results were interpreted based on the classification scheme of SPI, according to the World Meteorological Organization standards [60]. The resulting SSI time series was used to calculate drought characteristics, e.g., the frequency of dry/wet years, drought duration and severity (Drought Monitoring Indicators, DMIs) and the corresponding trends based on the MK trend test.

Trends in Historical Streamflow
The impacts of climate change on water resources in the LRB were investigated based on the assessment of trends using three hydrological statistical parameters, namely, the annual mean, annual minima and annual maximum of streamflow, computed from the observed and simulated streamflow from the ACRU hydrological model based on the CORDEX models, viz, CanESM2, IPSL-CM5A and Ensemble. The resulting MK trends in the selected hydrological statistical parameters are shown in Figure 3, which shows a comparable trend pattern between the observations and model simulations. In particular, a significant matching pattern is observed in the annual mean and annual maximum trends, where trends in annual streamflow across all observed (with exception to stream gauge station A6H021) and simulated streamflow depicts a negative trend in the annual mean and maximum streamflow. On the other hand, two stream gauge stations and the ensemble model depict positive trends in the annual minima, while the rest of the stations and CORDEX models exhibit negative trends, although in some cases, such trends are extremely minimal for the observed streamflow time series. Statistically significant trends at a 95% significance level are indicated in yellow. In general, trends across the hydrological statistical parameters and across the flow gauge stations and model simulations suggest that the LRB experienced reduced streamflow between 1976 and 2005.

Trends in Projected Streamflow
Impacts of projected future climate on water resources in the LRB were based on trends analysis of projected streamflow simulated from the CORDEX models under RCP8.5 and RCP4.5. The results of the selected hydrological statistical parameters are shown in Table 1 for three different time intervals, namely, the present period (defined here as the period from 2006 to 2035), near future, spanning 2036-2065 and the distant future covering the 2070-2099 period. In this study, both the historical  and present climate are considered as baselines to assess the impacts on the future climate on water security in the LRB. The RCP8.5 scenario, also referred to as "business as usual", relates to the expected outcome of climate impacts in a case where no efforts are made to reduce greenhouse gas emissions. The results presented in Table 1 depict negative trends across the three selected hydrological statistical parameters and periods under the RCP8.5 scenario, although most of the observed trends are not statistically significant at a 95% significance level. Consequently, there is a decrease in streamflow during 2006-2035, and such a pattern is likely to continue in both the near and distant future under the "business as usual" scenario. The RCP4.5 is a long-term stabilization scenario that considers the global emissions of greenhouse gases and land-use-landcover in a global economic framework, with the assumption of enforcing emissions mitigation policies [61,62]. This scenario covers almost all the sectors of the economy, including agriculture and land-use emissions [61]. The added features of the RCP4.5 scenario result in a slight change in trends of projected streamflow in the LRB. In particular, trends in the annual mean and annual minimum and maximum of the projected streamflow exhibit an increasing pattern across the three time intervals considered in the current research study.

Trends in Projected Streamflow
Impacts of projected future climate on water resources in the LRB were based on trends analysis of projected streamflow simulated from the CORDEX models under RCP8.5 and RCP4.5. The results of the selected hydrological statistical parameters are shown in Table 1 for three different time intervals, namely, the present period (defined here as the period from 2006 to 2035), near future, spanning 2036-2065 and the distant future covering the 2070-2099 period. In this study, both the historical  and present climate are considered as baselines to assess the impacts on the future climate on water security in the LRB. The RCP8.5 scenario, also referred to as "business as usual", relates to the expected outcome of climate impacts in a case where no efforts are made to reduce greenhouse gas emissions. The results presented in Table 1 depict negative trends across the three selected hydrological statistical parameters and periods under the RCP8.5 scenario, although most of the observed trends are not statistically significant at a 95% significance level. Consequently, there is a decrease in streamflow during 2006-2035, and such a pattern is likely to continue in both the near and distant future under the Water 2020, 12, 3299 8 of 20 "business as usual" scenario. The RCP4.5 is a long-term stabilization scenario that considers the global emissions of greenhouse gases and land-use-land-cover in a global economic framework, with the assumption of enforcing emissions mitigation policies [61,62]. This scenario covers almost all the sectors of the economy, including agriculture and land-use emissions [61]. The added features of the RCP4.5 scenario result in a slight change in trends of projected streamflow in the LRB. In particular, trends in the annual mean and annual minimum and maximum of the projected streamflow exhibit an increasing pattern across the three time intervals considered in the current research study. Flow quantiles are essential in hydrological studies as they relate to conditions of water resources and, more generally, provide information on the potential status of hydrological extremes such as drought and floods. In the current study, trends in the 0.1, 0.5 and 0.9 flow quantiles were calculated across the CORDEX models and periods of analysis, with results given in Table 2. The 0.1, 0.5 and 0.9 quantiles relate to dry, normal and wet conditions, respectively. Persistent dry and wet conditions over time are likely to lead to drought and floods (flash floods), respectively. The results in Table 2 show a decreasing trend in the low, medium and high quantiles across the CORDEX GCMs under the RCP8.5 scenario and the projected time intervals (with exceptions to the distant future projections). On the contrary, an increasing pattern in the low, medium and high quantiles is observed for model simulations under the RCP4.5 emission scenario for the 2006-2035 projected period. The near-future projections (2036-2065) depict negative trends in low, medium and high flow quantiles across the CORDEX models under the RCP4.5 emission scenarios. For the distant future, low flow is projected to decrease under the RCP4.5 scenario, whereas high flow is projected to increase under the same emission scenario. In general, the results suggest that while the LRB is likely to experience frequent dry conditions that can translate to drought, the region is also likely to experience wet conditions that can result in (flash) floods. Table 2. Trends (m s −1 ) in flow quantiles (low = Q 0.1 ; medium = Q 0.5 and high = Q 0.9 ) under the representative concentration pathways 8.5 and 4.5; bold indicates statistically significant trends.

Extreme Value Analysis
Results for the GEV distribution analysis are given in Table 3, with the respective estimated parameters, viz, the location, scale and shape, across the models and timescales. The location parameter of streamflow extremes ranges between 4.5 to 13.7 and 6.1 to 11.9 for model simulations under the RCP8.5 and RCP4.5 scenarios, respectively. The scale ranges from the lowest, 3.9, to the highest, 15.2, whereas the shape parameter of streamflow extremes is between 0.8 and 1.9 across the RCP scenarios and the time intervals. Comparison of the values of the shape parameter with the classes of GEV family of distributions suggests that the simulated streamflow extremes across the CORDEX models under both the RCPs and timescales are best described by the Fréchet type II extreme value distribution. Table 3. Parameters (location, scale and shape) of the generalized extreme value distribution, across the timescales and CORDEX models under the representative concentration pathways 8.5 and 4.5.

Assessment of Hydrological Drought in the Limpopo River Basin
The SSI is used in this study to evaluate projected hydrological drought conditions in the LRB. The study considered only 6-and 12-month accumulation periods, since these periods are related to hydrological drought, with inherent impacts manifested in water reservoirs. The results (only for SSI-6 are presented) for the three CORDEX models under the RCP8.5 and RCP4.5 scenarios are presented in Figures 4 and 5, respectively. In each figure, the positive SSI values correspond to wetter conditions, whereas negative values represent drier conditions. The overall interpretation of the SSI values is similar to the SPI [60]. The figures depict the frequency of occurrence of projected dry and wet conditions in months within the LRB. The SSI-6 under both RCP scenarios demonstrates a significant variation in seasonal to annual streamflow, which fluctuates between wet and dry months. Drought events are likely to be more frequent, but in most cases at shorter time scales. Furthermore, the LRB is projected to experience intra-annual moderate to very wet conditions that are expected to last for shorter periods. In general, the SSI-6 computed from the CORDEX model simulations under the "business as usual" scenario indicates that the LRB is likely to experience reduced streamflow in the near future as well as moving into the distant future. Furthermore, the region is likely to experience significant wet conditions, categorized as ranging from moderate to very wet. Such wet conditions are expected to make the LRB susceptible to flash flooding.

Proportion of Projected Dry/Wet Years
The SSI-12 time series computed from the CORDEX models under RCP8.5 and the stabilized RCP4.5 scenarios were used to assess the proportion of projected dry (drought) and wet (excess rainfall) years in the study area, focusing on the near and distant future periods. There is no universally scientifically accepted definition of either a dry year or a wet year. Subsequently, this study adopted definitions set by [63] whereby a drought year was defined as a year exhibiting an SSI

Proportion of Projected Dry/Wet Years
The SSI-12 time series computed from the CORDEX models under RCP8.5 and the stabilized RCP4.5 scenarios were used to assess the proportion of projected dry (drought) and wet (excess rainfall) years in the study area, focusing on the near and distant future periods. There is no universally scientifically accepted definition of either a dry year or a wet year. Subsequently, this study adopted definitions set by [63] whereby a drought year was defined as a year exhibiting an SSI value less than −1.00, whereas a wet year was defined as one having an SSI value of greater than 1.00. Any other years outside these definition terms (e.g., SSI ranging between −0.99 and 0.99) were classified as years exhibiting near-normal conditions. Furthermore, in the present study, the dry/wet episodes in a year are determined as a percentage of counts of the respective dry and wet occurrences (as per the SSI thresholds stated in the foregoing) over a 12-month accumulation period. The results of this analysis are illustrated in Figures 6 and 7 for the near-and distance-future time intervals, respectively. Based on model simulations under the RCP8.5 scenario, given in Figure 6a

Trends in Drought Monitoring Indicators
Drought monitoring indicators-the drought duration and drought severity, were computed using the SSI-6 and -12 time series across the CORDEX models under the RCP8.5 and RCP4.5 scenarios. The DMIs' trends across the models, RCP scenarios and study periods are presented in Table 4. While simulations under the RCP8.5 scenario depict zero trends in drought duration (DD), for SSI-6, across the models and study time intervals, negative and positive trends are observed under the RCP4.5 scenario, for the near and distant future, respectively. Moreover, negative trends in the drought severity are observed across the models, RCPs and periods of study, with exceptions for simulations under the RCP4.5 stabilized scenario, for the distant future. Furthermore, negative trends in DD (SSI-12) are observed for the near-future epoch and no trends for the distant future under both RCP scenarios. A similar trend pattern is observed for drought severity (DS) across the RCPs with the exception of the distant future, which shows positive trends in DS. Only a few of the observed positive (negative) trends are statistically significant at a 95% significance level (see the trends in boldface). In general, drought conditions are projected to persist in the LRB during the entire period, (near to distant), as confirmed by the drought episodes in Figures 4 and 5, although such conditions are likely to be less severe. Table 4. Trends in drought monitoring indicators based on standardized streamflow index 6 and standardized streamflow index 12; bold represents statistically significant trends at 95% confidence level. DD-drought duration, DS-drought severity.

Discussion
The LRB has, in the past, attracted interest from various stakeholders, e.g., from the science community, in-country and transboundary hydro-governance structures and institutions, as well as the general society. Such interests are attributed to, among other things, the significant socio-economic role the basin plays, e.g., to sustain the livelihoods of the 14 million people across the four co-riparian states of Botswana, Mozambique, South Africa and Zimbabwe. Despite its critical importance, the LRB is considered as an arid to semi-arid region and one of the most stressed water basins in Africa [64]. Consequently, water resources are a limiting factor for development in the region, and this is exacerbated by climate change and variability. Water demand for South Africa from the LRB is currently in the order of 3000 Mm 3 /year [64], of which irrigation uses the most water, followed by urban supply, mining and power production. The impacts of climate change and variability, coupled with other factors, such as population growth, urbanization, industrial development and increasing agricultural activities [65], are likely to exacerbate the strain on water resources in the basin, with variable impacts on different water-linked sectors, as well as on socially and economically disadvantaged groups, who depend directly on these resources for their livelihoods. Climate change is expected to bring changing rainfall patterns for the LRB, resulting in more frequent droughts and floods, which will impact agriculture, food security and other water-linked sectors.
The results reported in this study show that water resources in the LRB are already stressed, as supported by the observed decreasing trends in historical streamflow from 1976 and 2005, which are in accordance with studies reported by [66] on annual variations of rainfall, minimum and maximum temperature climate variables and [67][68][69] on the frequent occurrences of drought over Southern Africa. Such trends are likely to continue considerably in the future, as a result of the projected high variability of rainfall patterns and increases in hydro-climatic extremes within the region. These changes are expected to reduce the water available for irrigation (which is already in high demand) and dryland farming in the basin, with impacts on food security, decrease in crop yields, the shift in planting dates and cropping patterns, changes in the suitability of crops and crop varieties, loss of wildlife, livestock and income. To alleviate future climate impacts on agricultural systems, there is a need to consider an integrated approach of managing the groundwater and surface water systems, as an alternative option to managing water resources under the changing climate [4].
Based on climate projection models, the LRB is likely to continue to experience reduced streamflow in the near and distant future projected periods; these findings corroborate the results of [70], although the projections were made up to 2030. In particular, projected hydroclimatic extremes reported in this study, e.g., reduced streamflow, observed low and high flow patterns, sporadic dry/wet episodes in the region, as well as the general persistence of drought (the duration and severity), will have more pronounced impacts on the different water-linked sectors of the economy, in the near future. The sporadic wet episodes will lead to flash floods especially in flood-prone areas such as the lower Limpopo River floodplain, and this is a challenge in itself because the most vulnerable often resist relocation due to their dependency on agricultural activities along the floodplain. Unless efficient early warning systems are put in place to support a timely response to impending floods, the most vulnerable communities will continue to suffer from the harshness of climate change and variability.
The likelihood of occurrences of floods in the LRB will also increase the incidence of vector-borne disease, which affects both humans and animals. Furthermore, the projected dry episodes under climate change in the LRB will exacerbate the frequent occurrences of drought. It is well known that malaria is prevalent in the LRB, and it is characterized by diverse malaria transmission stabilities where the highland regions have no malaria transmission, the downstream to midland regions have unstable transmission and the downstream to lowland regions have stable malaria transmission [71]. Consequently, the projected dry episodes will lead to a geographical shift of malaria within the different regions of the basin, whereas flash floods will increase the risk of malaria [72], cholera, diarrhea and bilharzia diseases, which are a notable concern in the LRB [33]. Moreover, projected hydroclimatic extremes in the LRB will lead to ill-health related diseases, such as nutrition-related, water-related, airborne and dust-related and vector-borne diseases and mental health disorders, among others [73]. Furthermore, the projected dry episodes under climate change in the LRB will exacerbate the water demand for power generation, taking into account that the already existing coal-based power plants in the Limpopo Province require a considerable amount of water. The use of renewable sources such as solar photovoltaic systems that require very minimal amounts of water may be beneficial for the LRB under changing climatic conditions.

Conclusions
This contribution has used CORDEX-Africa model simulations to assess the impacts of future climate change on water-linked sectors in the LRB of South Africa. Based on the analysis of projected streamflow statistical metrics, the findings of the study can be summarized as follows: • The LRB is experiencing frequent dry conditions (decrease in streamflow) under the current climatology. The conditions are projected to continue in both the near future and towards the end of the century time intervals. Prolonged dry conditions can translate to drought in the basin.

•
The region is also likely to experience wet conditions at an average that can result in (flash) floods. • Significant dry and wet years are projected in the LRB in the near future and as we move into the distant future. The CORDEX model analysis under the RCP8.5 and RCP4.5 scenarios project that the LRB will experience, on average, 64% (46%) and 45% (55%) of dry (wet) years, respectively in the near future. As we move towards the end of the century period, the basin is projected to experience an average of 48% (52%) and 55% (45%) of dry (wet) years when considering GCM simulations under the RCP8.5 and RCP4.5 scenarios, respectively.

•
The hydrological extremes are best represented by the Fréchet generalized extreme value distribution over the LRBs. This GEV family could be ideal for calculating the return periods of extreme events in the LRB.
The findings from this study can be used to support water resource management, agricultural planning and effective decision making to attain sustainable practices regarding water resources in the LRB. The results contribute towards enhancing the ability of farmers, catchment management agencies and extension services within the basin to develop sustainable adaptation strategies that can be used to alleviate the impacts of climate change.