Assessing Freshwater Changes over Southern and Central Africa (2002–2017)

: In large freshwater river basins across the globe, the composite inﬂuences of large-scale climatic processes and human activities (e.g., deforestation) on hydrological processes have been studied. However, the knowledge of these processes in this era of the Anthropocene in the under-studied hydrologically pristine South Central African (SCA) region is limited. This study employs satellite observations of evapotranspiration (ET), precipitation and freshwater between 2002 and 2017 to explore the hydrological patterns of this region, which play a crucial role in global climatology. Multivariate methods, including the rotated principal component analysis (rPCA) were used to assess the relationship of terrestrial water storage (TWS) in response to climatic units (precipitation and ET). The use of the rPCA technique in assessing changes in TWS is warranted to provide more information on hydrological changes that are usually obscured by other dominant naturally-driven ﬂuxes. Results show a low trend in vegetation transpiration due to deforestation around the Congo basin. Overall, the Congo (r2 = 76%) and Orange (r2 = 72%) River basins maintained an above-average consistency between precipitation and TWS throughout the study region and period. Consistent loss in freshwater is observed in the Zambezi ( − 9.9 ± 2.6 mm/year) and Okavango ( − 9.1 ± 2.5 mm/year) basins from 2002 to 2008. The Limpopo River basin is observed to have a 6% below average reduction in rainfall rates which contributed to its consistent loss in freshwater ( − 4.6 ± 3.2 mm/year) from 2006 to 2012.Using multi-linear regression and correlation analysis we show that ET contributes to the variability and distribution of TWS in the region. The relationship of ET with TWS (r = 0.5) and rainfall (r = 0.8) over SCA provides insight into the role of ET in regulating ﬂuxes and the mechanisms that drive precipitation in the region. The moderate ET–TWS relationship also shows the effect of climate and anthropogenic inﬂuence in their interactions.


Introduction
Terrestrial water storage (TWS) is a function of changes in the earth's gravity fields. Its relation to hydrological applications is characterized by the sum of water stored underground, in the soil, as snow, over the land surface, and canopy. The changes in any of the above components remain one of the most critical elements in the balance of the hydrological cycle. Monitoring the changes that occur within these TWS components is an important hydrological exercise, as it helps in determining the cause of climate change, thus making a cogent recommendation as to the required environmental response necessary to mitigate potential climatic hazards at the basin or regional scale. With the advent of satellite remote sensing and altimetry, extensive data collection and observation that have been recorded in hydrology, especially in the monitoring of water storage and fluxes in a changing world, thus enabling a universal evaluation that spans political boundaries [1][2][3].
An important intention of this study is to understand the hydrological response of the South Central African (SCA) region to climatic conditions. Amidst other factors that contribute to climate change, the study focuses on two notable water cycle variants; evapotranspiration (ET) and precipitation. The net-zero water cycle is achieved majorly as a result of these two factors. Spatial and temporal variability in rainfall amount and evapotranspiration indexes as a result of global climate change is posing a serious danger to water resources management, food security, rain-fed agriculture, and poverty reduction, especially in developing countries such asthe SCA region. Therefore, understanding the variability, magnitude, patterns, and drivers of precipitation and evapotranspiration is crucial for studying the climate systems of the SCA region. Following substantial environmental shifts such as deforestation, increased temperature patterns due to anthropogenic climate change, fluctuations in precipitation patterns amongst others being reported around the SCA region [4][5][6][7][8][9][10][11][12][13], a potential climatic backlash is inevitable. We therefore attempt to comprehend the region's response to these conditions through a short-term trend analysis between 2002 and 2017 to make informed decisions about managing the expected regional climatic shifts.
Southern and Central Africa is a region of strong climatic hotspots that play key roles in global climate. For instance, the Zambezi, Orange, Okavango and Limpopo basins in Southern Africa discharge approximately 400 cubic km of water every year, while the Congo basin's rainfall climatology in Central Africa overshadows global tropical rainfall during transition seasons, thereby contributing to the global climatic impacts, as well as water resources and socio-economic changes [14][15][16][17][18]. Additionally, an estimated 25% to 40% of ET from the Congo River basin is recycled as local rainfall, and model simulations show that variations in ET rates within the basin area affect the moisture cycling phase across the African continent [4]. The temporal evaluation by Mohammed and David (2019) [19] between 2002 and 2017 showed a significant variation in TWS over the Congo basin in Central Africa, which extends to the southwestern part of Sudan. Furthermore, intense water deficit anomalies in Central Africa between 2005 and 2007 [14,20], and the recent progressive negative trend in TWS of SCA (Mohammed and David, 2019; Hua et al., 2016;Sorí et al., 2017) [10,20,21], are undoubtedly clear impacts of climate variability, as well as pointers to climate change in the Central and Southern African region.
Assessing and monitoring variations in hydrological patterns at a global scale is practically difficult with the use of insitu observations alone, especially due to the high cost of installing and maintaining station controls, as well as the existence of gaps in gauged data (Rodell et al., 2018) [22]. With the advent of satellite remote sensing in monitoring water storage and fluxes, several studies have shown how satellite observations can be used to improve understanding of globalhydrology (Burnett et al., 2020;Ahmed et al., 2019) [23,24]. The launch of the Gravity Recovery and Climate Experiment (GRACE) satellite mission made this possible, including tracking of ice-sheet and glacier ablation, monitoring of groundwater depletion, and impacts of climate change on freshwater variations (e.g., Tapley et al., 2004;Wahr et al., 1998;Rodell and Famiglietti, 1999;Swenson et al., 2006;Cazenave et al., 2010) [25][26][27][28][29].
The trend study carried out by Ahmed et al., (2019) [19] and Rodell et al., (2018) [22] used the breakpoint algorithm and apparent trend analysis, respectively. This breakpoint algorithm, which identifies the structural changes in the regression relationship and splits the time series into different segments, suffers from a high snapping tendency causing instability in the run-time system in the event of inexact medians [67,68]. In a bid to improve our understanding of the land water storage over SCA and to monitor the changes and hydro-geological fluctuations at the major river basins in the region, this study attempts to emphasize the recent annual and seasonal inconsistencies of the region's TWS within the period between 2002 and 2017 by estimating the variability trends in TWS using the rotated principal component analysis. The rPCA technique is an improvement on the normal PCA technique which assesses extra elements in the reduction of the dimensionality of a dataset, as well as discovering hidden structures in the dataset. Since the rPCA technique is a non-parametric method, it can be applied todata with any distribution. Analogous to the study carried out by Ujeneva and Abiodun, (2015) [52], the rotated varimax option is implemented to distinguish spatial coherence (modes) in the temporal variability of the TWS. Its prevalence over other methods of trend analysis (breakpoint and apparent trend analysis) is due to its dedication to the varied retained variability modes, as well as the ease in correlating against them   [69,70]. Besides academic curiosity, this work equally attempts to assess the merit of the study undertaken by Rodell et al., (2018) [22] who used the apparent trend analysis to observe a significant increase and substantial wetting trend within the SCA region between 2002 and 2016. The study also reported a negative trend along the coast of Southeastern Africa (housing the Limpopo River basin) of −12.9 ± 2.3 Gt year −1 , reflecting a recent severe drought in the region (Figure 1g,h). The intended approach is based on the analysis of the relationship and existing variability between the GRACE-derived TWS changes and other relevant remote sensing datasets applied in this study. This region-based investigation attempts to quantify and assess the short-term trends in GRACE-derived TWS over SCA as well as explore the factors (natural or anthropogenic) contributing to the trend. Therefore, with the implementation of a rotated principal component analysis (rPCA), and multiple linear regression analysis (MLRA), this study looks into inter-annual and seasonal fluctuations in TWS variations using GRACE and climate drivers over SCA, i.e., precipitation and evapotranspiration patterns in the region. This study specifically seeks to (i) distinguish trends and TWS variability hotspots in relation to rainfall and ETa over SCA (ii) improve understanding ofthe seasonal cycle of TWS and attempt to make future predictions of hydrological state and patterns (iii) study the rate of evapotranspiration as a result of temperature fluctuations as well as its relationship to TWS variability.
The remainder of this study is segmented as follows; Section 2 provides a brief introduction of the study area and analysis coverage zone; Section 3 gives a thorough discussion on the data and methodology employed. This is followed by the analysis and discussion of the results in Section 4, and finally the conclusion of the study is presented in Section 5.
correlating changes in GRACE-derived TWS with changes in climate derived output (e.g., TWS vs. ETa).
The remainder of this study is segmented as follows; Section 2 provides a brief introduction of the study area and analysis coverage zone; Section 3 gives a thorough discussion on the data and methodology employed. This is followed by the analysis and discussion of the results in Section 4, and finally the conclusion of the study is presented in Section 5.

Location
Our study area covers the southern and central African region between 8.5° E to 40.8° E and 11.5° N to 35° S. The study domain is characterized by a latitudinal gradient that includes some major basins (For example, Congo, Zambezi, Limpopo, Okavango and Orange basins) (see Figure 1), whose combined effect has contributed significantly to the climatic fluctuations of the region, and by extension the entire continent.

River Basins
The Congo basin being the third largest in the world plays a significant role in global climate change. It is located in central-equatorial Africa, consisting of several African countries including, the people's Republic of Congo, the Democratic Republic of Congo (DRC), the Central African Republic, and parts of Zambia, Cameroon, Angola, and Tanzania (Sorí et al., 2017) [21]. Despite the role of the Congo basin at global and regional

Location
Our study area covers the southern and central African region between 8.5 • E to 40.8 • E and 11.5 • N to 35 • S. The study domain is characterized by a latitudinal gradient that includes some major basins (For example, Congo, Zambezi, Limpopo, Okavango and Orange basins) (see Figure 1), whose combined effect has contributed significantly to the climatic fluctuations of the region, and by extension the entire continent.

River Basins
The Congo basin being the third largest in the world plays a significant role in global climate change. It is located in central-equatorial Africa, consisting of several African countries including, the people's Republic of Congo, the Democratic Republic of Congo (DRC), the Central African Republic, and parts of Zambia, Cameroon, Angola, and Tanzania (Sorí et al., 2017) [21]. Despite the role of the Congo basin at global and regional scales in terms of both water discharge and drainage basin size, it is reported to have a historical evolution of hydro-meteorological events, climatological, flood and drought episodes (Zhou et al., 2014;Hua et al., 2016;Ndehedehe et al., 2018a) [10,12,16]. The Zambezi river basin is situated at the southern end of the Congo river basin, with an experience of immense precipitation between June and September every year. The drainage is the major supplier of electricity, freshwater and fish to the populace of Angola, Botswana, Namibia, Tanzania, Zimbabwe, Zambia, Mozambique, and Malawi. Amidst housing an extensive wet plain, the Zambezi River basin is also accountable for the climate regulation Remote Sens. 2021, 13, 2543 5 of 30 of a prolific ecosystem of savannas and humid forests surrounding it. The Zambezi River basin maintains an eastern outlet, into the Indian ocean. Towards the southwestern region of the Zambezi River basin lies the Okavango delta, draining the wet highlands of Angola. The Okavango basin is an endorheic river basin covering most of Namibia and Botswana, including parts of South Africa, Zambia, and Angola. It maintains an outstanding physical feature which houses a large semi-arid sandy savannah. During the fall and winter months, water evaporates from the Okavango delta (a vast inland delta formed from the basin's downstream end) situated in northern Botswana. Towards the southeastern part of the Okavango and the southern part of the Zambezi lies the Limpopo River basin. This basin is endowed with groundwater resources capable of supplementing surface water resources in parts of Botswana, South Africa, Mozambique, and Zimbabwe. The Orange River basin also classified as a passive volcanic margin basin, makes its way from the highlands in the east through the Kalahari depression in the west, to empty into the South Atlantic Ocean. It has maintained a steady precipitation decline around the southern and eastern parts compared to the northern part. The rainfall season takes place around December to February within the summer months, while the dry season peaks around November, within the late spring period. Due to the hydrological complexity and lack of monitoring stations in the northern part of the Orange River basin, a large seasonal wetland in the region is left unmonitored.  [81,82]. A large overestimation in the tropical rainforest region of Africa ( Figure 1) over the TRMM precipitation radar is determined in December-January-February and in March-April-May and a minimized bias in June-July-August and September-October-November using zonal mean analysis. However, the bias which is extensively high for all models in the dry seasons when rainfall is minimal is less prominent in the dry seasons of the Southern and Central African climatic regions. Additionally, in order to preserve and maintain a consistent spatial resolution with the other datasets used in this study (such as the GRACE-derived TWS solutions), it was necessary to resample the TRMM/IMERG precipitation dataset.

Evapotranspiration (ET)
Remote sensing products have long been identified as the most feasible and practicable means of providing spatially distributed ET information over the terrestrial surface. The evapotranspiration dataset used in this study was provided by the Numerical Terradynamic Simulation Group (NTSG), University of Montana. MODIS16A2 global ET monthly datasets from 2002 to 2017 were obtained over an 8-day temporal period at a 0.5 km pixel resolution. The calculation of this ET data is based on the Penman-Monteith equation (Monteith, 1965) [83] and identifies the evaporation from the soil, wet canopy, as well as plant transpiration at dry canopy surfaces (Mu, Zhao, and Running, 2011) [84]. However, these data were gridded again so as to correspond with other datasets used in this study.

Rotated Principal Component Analysis (rPCA)
This study applied the rPCA technique on the TWS dataset to be able to assess extra variability dimensions, as well as uncover hidden structures in the dataset ( Figure 2). The idea of the rPCA technique is hinged on the PCA technique. The principal component analysis-based cumulant scheme is basically used to statistically decompose a given dataset into diverse dimensions while retaining interpretability and diminishing data loss (e.g., Ndehedehe et al., 2016b;Ziehe, 2005;Jolliffe, 2003) [15,69,70].The simplicity and proficiency in the explanation of optimized variance parameters, with a minimum number of principal components, contributed to the popularity and relevance of this statistical data analysis tool in climate science (Ndehedehe et al. 2016a [14,36].The rPCA gives the best possible illustration of a p-dimensional dataset in q-dimensions (q < p) by attempting to maximize the variance in q dimensions. For regularization, outliers and gross errors were eradicated from the observation matrixŷ k in order to achieve a robust rotated principal component analysis (e.g., Jolliffe, 2003) [69,70]. The need for the identification of periodic signals embedded in multi-resolution datasets is increasing, and has resulted in the application of the robust PCA to several vigorous lines of research such as; image processing, machine learning, bioinformatics, or web data analysis including the regionalization and analysis of global geophysical time series data (e.g., Ndehedehe [85][86][87][88]. Its unique and simple interactions allow for its broad range of applications in spatial science (Westra et al., 2009) [89], the choice of this technique in this study is due to its ability to linearly combine centered (original) variables provided their units are the same. It is common practice to use a predefined percentage of the total variance in deciding the number of PCs to be retained (70% of total unpredictability is normal, if subjective, cut-off point), although the need to graphically represent the information causes the utilization of just the first six rPCs. Regardless, the percentage of total variance accounted for is basic in the determination of the quality of low dimensional visual representation of the dataset (Ndehedehe, 2019) [85]. The quality of any q-dimensional estimation can be measured by the inconsistency of the set of retained PCs. The sum of variances of the p original variables is the trace of the covariance matrix S. The standard measure of the quality of a given PC is represented by a simple matrix theory which attempts to show the proportion of total variance accounted for, where tr(S) represents the trace of S. The inherent incremental nature of the PC shows theproportion of total variance explained by a set U, and is often denoted as a % of total variance accounted for; ∑ j∈U π j × 100%. PCA reduces the dimensions of multivariate data by linearizing functions of the original variables and parsing them into new variables.
where ( ) represents the trace of S. The inherent incremental nature of the PC shows theproportion of total variance explained by a set U, and is often denoted as a % of total variance accounted for; ∑ × 100% ∈ . PCA reduces the dimensions of multivariate data by linearizing functions of the original variables and parsing them into new variables. The linear combinations for K principal components at a given discrete-time T aregiven by: The linear combinations for K principal components at a given discrete-time T are- The state matrixx T,K contains rows corresponding to the time T in months, and the variables K represents the observations, e.g., TWS. In the linear arrangement above (Equation (2)), the y values are orthogonal, while y T,1 to y T,K represents the rPC signifying the variance of the system matrix with the initial rPC explaining the highest variability. The weight of the original variables in the rPCs is provided by the coefficients of the linear observations (eigenvectors) also known as factor loadings. The factor loadings also referred to as the component coefficient is determined by dividing the factor's eigenvalue by the total number of variables. While the normalized eigenvectors (e.g., P 11 , P 12 , P 13 , . . . , P KK ) are the component coefficients, each characteristic root measures the amount of variation in the total sample accounted for each factor (Ndehedehe et al., 2016a) [14]. A factor's characteristic roots may be computed as the sum of its squared component coefficient for all the variables.

Multiple Linear Regression Analysis
The multi-linear regression analysis is a robust statistical technique that applies a least-squares approach in modeling hydrological and climate datasets. In its application, the relationship between a dependent variable (e.g., TWS changes) and an independent variable (e.g., precipitation patterns) is evaluated to model the trends of precipitation time series and GRACE-derived TWS.To assess the seasonal and inter-annual variabilities existent in the TWS, ETa and precipitation datasets at every grid point, the regression model, is matched to the temporal patterns of the data (see Rieser et al., 2010) [90]. Where, D T/R is a time function (t * ) of either the TWS or precipitation dataset, γ 0 and γ 1 refers to the constant offset and linear trends respectively, γ 2 and γ 3 describes the annual signals, while γ 4 and γ 5 accounts for the semi-annual signals, the eccentricity between the derived observations and model outputs are represented by z. The Root Mean Square Error (RMSE), annual and semi-annual amplitudes are given by, where obs D T R and sim D T R represent derived observations and model outputs from Equation (3), respectively, for the monthly interval mts. As a measure to validate the model strength in the simulation of the observed parameters, the phase lag, as well as the correlation of determination for each grid point, was computed. With the use of other assessment parameters such as the p-value, coefficient of determination and so on, the multi-linear regression method maintains a unique view of accuracy assessment for climate and hydrological studies. The use of this technique in modeling the trends, correlation as well as the harmonic constituents of TWS changes (such as ET, and precipitation patterns) over the SCA region provide a distinctive approach in the accuracy assessment of this study.

Variability Index, Trends, and Annual Amplitudes
The need to classify the time series of the changes in TWS and rainfall residuals over the SCA region resulted in the analysis of the monthly variability of trends and seasonal signals for the two products. In order to analyze the trend variations as well as examine the seasonal variations in the SCA region, the variability index from which the semi-annual and annual components were removed was computed for the two products (TWS and rainfall). The time series of changes in the two products comprise trends and seasonal components which can be interpreted into different climatic schemes such as normal, extreme, wet and dry.
where, ∆ TWS and ∆ RF represents the variability index for the TWS and rainfall, respectively, D 1 and D 2 represent the TWS and rainfall data, respectively, averaged over the SCA region. µ and σ are the corresponding mean and standard deviation for the data over the period of study. The variability index for both products was computed as standardized deviations. In addition to this, an average of the TWS and rainfall plots was made, to understand the evolving annual amplitude of rainfall and the extent of TWS changes over the study region through the temporal patterns. Furthermore, based on the least-squares estimation of the regression coefficient, trends in rainfall and TWS were estimated for each grid cell.
where β is the least-squares vector of unknowns, .
A is the design or coefficient matrix, . P is the weight matrix of the observations, L b is the vector of observations for both products. The estimation can be positive denoting a positive increasing trend or negative denoting a decreasing trend. A popular non-parametric method used in testing the significance of trends was employed in this evaluation. Applying the seasonal Mann-Kendall test, the significance of the observed trends was checked (Mann, 1945;Kendall, 1970;Hamed, 2009) [91][92][93]. To test for parametric trends, tests such as the inversion, regression, turning point, and student T-tests are employable options in the examination of the statistical significance of trends. However, factors such as data weight and type tend to affect them, and this could result in a violation of the normality of these models ( Figure 3).
where, ∆ TWS and ∆ RF represents the variability index for the TWS and rainfall, respectively, D1 and D2 represent the TWS and rainfall data, respectively, averaged over the SCA region. and are the corresponding mean and standard deviation for the data over the period of study. The variability index for both products was computed as standardized deviations. In addition to this, an average of the TWS and rainfall plots was made, to understand the evolving annual amplitude of rainfall and the extent of TWS changes over the study region through the temporal patterns. Furthermore, based on the least-squares estimation of the regression coefficient, trends in rainfall and TWS were estimated for each grid cell.
where is the least-squares vector of unknowns, ̇ is the design or coefficient matrix, ̇ is the weight matrix of the observations, is the vector of observations for both products. The estimation can be positive denoting a positive increasing trend or negative denoting a decreasing trend. A popular non-parametric method used in testing the significance of trends was employed in this evaluation. Applying the seasonal Mann-Kendall test, the significance of the observed trends was checked (Mann, 1945;Kendall, 1970;Hamed, 2009) [91,92,93]. To test for parametric trends, tests such as the inversion, regression, turning point, and student T-tests are employable options in the examination of the statistical significance of trends. However, factors such as data weight and type tend to affect them, and this could result in a violation of the normality of these models ( Figure 3). This however resulted in the predilection of non-parametric methods (such as the Mann-Kendall, Hotelling-Pabst, and Sen tests amongst others) over their parametric counterparts especially given their significance in trend detection (Gocic and Trajkovic, 2013;Gentilucci et al., 2021) [94,95].
where, MK is the non-parametric Mann-Kendall statistic, n is the number of data locations, D j and D i represent data values at time j and i (j > i).
The Mann-Kendall test statistic represents the negative and positive transformation for all significant data points. Under the null hypothesis, the statistics mean (E[M]) = 0, and the variance ( Depending on the Z-transformation given below, The MK test statistic is estimated to be Gaussian. The null hypothesis (H 0 ) which signifies no trend was tested at a 95% confidence level. Since the critical value of the standard normal distribution was less than the computed absolute value of the test statistic, the hypothesis of the positive or negative trend was accepted at α = 0.05.

Spatio-Temporal Variability of Rainfall and Terrestrialwater Storage
This section discusses the spatio-temporal evolutions of prominent rainfall and terrestrial water storage modes from 2002 to 2017 over the SCA region. Its relevance is due to the complex hydrological characteristics of the SCA region which in turn plays a key role in the climate variability in Africa and the world at large. Using the multi-linear regression algorithm, strong relationships were observed between TWS and rainfall.

Terrestrial Water Storage
Regarding the analysis of GRACE-derived TWS over the SCA using the rPCA technique, the dominant spatio-temporal patterns of the changes in TWS over the SCA region between the 2002 and 2017 period are discussed here. For ease of interpretation, the TWS grid is decomposed into sets of rotated principal component analysis (rPCA), with their equivalent empirical orthogonal functions (rEOFs, i.e., variations in spatial patterns) with the use of the rPCA covariance technique. The standard deviations of the corresponding rPCA modes were used in the standardization of the EOFs as well as normalizing the rPCs to make them unitless. The statistically significant rPCAs, which are represented by the first six rPCA modes (Figure 3a-f), gave a cumulative variance of 94.51% and were assumed as significant signals denoting over 95% of the total TWS fluctuations over the SCA region. From the rPCA analysis, it is deduced that the highest EOF loadings are observed in a large portion of Central Africa and the northern portion of Southern Africa in the first rotated orthogonal mode (Figure 4).
In the first rPCA mode, the annual variability of TWS changes in the SCA region is represented by a variance of 37.76%. It can therefore be deduced that the strongest annual variability occurred over the Congo River basin which consists of the democratic republic of Congo, Angola, Zambia, Malawi, and the northern part of Mozambique. Especially over Zambia, Malawi and the Democratic Republic of Congo, the significance of rainfall makes the region remarkably humid. Moreover, rEOF1 ( Figure 4a) and its corresponding rPC (Figure 5a) shows an oscillating trend in TWS over the SCA region which peaked in2002, 2008, 2011 and 2017. The analysis also shows a minimum of TWS changes during 2007. The dramatic leap in trend experienced between 2007 and 2008 in TWS fluctuation over the study region from −2 to 1.5 can be attributed to an unusual increase in multi-annual and annual rainfall rate experienced in the region. However, in the seasonal analysis represented in Table 1, we observe that the region had its greatest gain in TWS between March-May and had the greatest loss in TWS between September-November. While these seasons are regarded as wetting seasons throughout the basin region, we observe a significantly higher precipitation rate from March to May. The high rate of TWS loss from September to November is because of the high transfer of latent heat during this period, as well as anthropogenic factors. The second EOF and its corresponding rPC, which explains 19.97% of the total TWS variability represents another annual TWS variation with a significant negative trend over Zambia, Malawi, and a part of Angola. This is well detailed in the downward trend shown in Figure 4b from 2005 to 2011. Thereafter, an average signal was experienced from 2011 to 2016, before the surge in 2017 and 2018. Since the second rPCA mode describes a total TWS variability of 19.97%, the trends in annual variations relative to the corresponding EOF loadings over Zambia, Zimbabwe, Malawi and both positive and negative trends in North and North Western Angola respectively. Our attention is particular drawn to the Zambezi River basin which observed a significant decline in TWS averaged over the study period. With an estimated mean annual rainfall of 100 to 140 mm in the north of the Zambezi River basin, a steady decline in the above figure is observed towards the south and south-eastern portion of the basin. A multi-annual precipitation rate of around four to six months summer rainy season is observed when the Inter-Tropical Convergence Zone (ITCZ) moves over the river basin from the north between October and March. This, therefore, results in a high ET and ETa rate (200−230 mm) as well as significant water loss in swamps and floodplains especially in the south-western portion of the river basin. From Figure 4b, a steady decline around the Zambezi River basin (crossing or either forming boundaries with Angola, Zambia, Namibia, Botswana, Zimbabwe and Mozambique) was observed from 2006 to 2012. This steady decline in TWS can be attributed to the decline in precipitation pattern within this period, high ET rate (Table 1) and the resultant anthropogenic effect in this regard. An interesting fact as pointed out by Rodell et al., 2018 [22] is that within this period, the Okavango delta and areas west of the coast experienced a notable change in theirhydro-climate. Given that the area-averaged annual rainfall from 1979-2015 was 970 mm, the threshold was increased over five times between 2006 and 2011. Initially, a permanent climate shift was speculated due to the significant decrease in annual rainfall rate between 1950-1975 and 1980-2005. Based on these observations and re-analysis, we therefore attribute natural variability to this GRACE period trend.   July-August-September rainfall (Figures 3g and 5c,d) over the study period. Due to the small sizes of these rivers and lakes and their proximity to the Congo River basin, its contribution to global climate change is usually attributed to the Congo River basin during global re-analysis. In the corresponding rPC plot, a steady decline in TWS is observed over these smaller river basins and the Congo River basin between 2002 and 2006. This is followed by an oscillating trend which is averaged as an increase in TWS over the above-mentioned basins (Figure 4c). The third rPCA mode representing another annual variation shows 7.20% of the total TWS variability. This variation is centered over the Congo River basin with a strong focus on smaller river basins west of the Congo River basin. They include the Cross Sanaga Bioko coastal forests, Ogooue River basin, San Benito River, Utamboni River, Mbe and Nyanga River, and Lake Cayo. The trend in TWS over this region is due to the July-August-September rainfall (Figures 3g and 5c,d) over the study period. Due to the small sizes of these rivers and lakes and their proximity to the Congo River basin, its contribution to global climate change is usually attributed to the Congo River basin during global re-analysis. In the corresponding rPC plot, a steady decline in TWS is observed over these smaller river basins and the Congo River basin between 2002 and 2006. This is followed by an oscillating trend which is averaged as an increase in TWS over the above-mentioned basins (Figure 4c).
In the fourth rPCA mode representing another annual variation, a 7.93% total TWS variability is recorded. An increase in TWS is observed over the northeast of the Democratic Republic of Congo, north western Zambia, and south eastern Angola. This maps a portion of the Congo River basin. Additionally, an increase in TWS is observed in Malawi, Mozambique and Zimbabwe which houses other smaller rivers such as Lake Chilwa, Pugwe, Sabi and Buzi. Additionally, two basins of interest (i.e., Zambezi and a portion of the Limpopo River basin) benefit from this increase in TWS changes.
The fifth rPCA presents a bi-modal multi-annual variation, with a record of 6.21% of total TWS variability. The corresponding EOFs (Figure 4e) show a better representation of a positive and negative trend around Angola/Namibia and northern Zambia/Malawi respectively. This indicates the presence of existing water lakes in these regions. Upon close investigation, it was observed that the Angola/Namibia positive trend is attributed to the existence of the Otjikoto lake, Ntwetwe, Sua Pan and smaller basins existent west of the Zambezi River basin. To confirm this, the time series of Figure 4e is correlated against the time series of the Otjikoto Lake fluctuation, to obtain a 0.8-degree correlation (Not shown).
The sixth rPCA represents an annual TWS variation with a 15.44% record of total TWS variability (Figures 3f and 4f). The positive variations occur around the Central African Repulic and the Otjikoto Lake, while strands of negative variations occur around Gabon and Northern Zambia. The reduction in TWS changes observed around Gabon and Northern Zambia stems from the dry January to June periods (Figure 6a,b) predominant in the study region. However, Figure 1d agrees with the study undertaken by Mohammed and David (2019) [19] which reported a negative trend in TWS within the Central African Republic, part of DR Congo, Republic of the Congo and Gabon between 2002 and mid-2007.  Figure A3), precipitation, and evapotranspiration within the five major river basins in SCA from 2002 to 2017. The statistical report of each basin shows the relationship intercept (inter.) between the three products (unit:mm). African Repulic and the Otjikoto Lake, while strands of negative variations occur around Gabon and Northern Zambia. The reduction in TWS changes observed around Gabon and Northern Zambia stems from the dry January to June periods (Figure 6a,b) predominant in the study region. However, Figure 1d agrees with the study undertaken by Mohammed and David (2019) [19] which reported a negative trend in TWS within the Central African Republic, part of DR Congo, Republic of the Congo and Gabon between 2002 and mid-2007.

Rainfall Variability
The seasonal precipitation patterns that exist in the Central African region are relatively higher than in the Southern African region. This could be argued to be because of the presence of significant river basins in the central region such as the Congo River basin. Apart from the OND period (Figure 3h), which suggests the presence of an average amount of rainfall in the Southern African region (i.e., Namibia, Botswana, Lesotho,

Rainfall Variability
The seasonal precipitation patterns that exist in the Central African region are relatively higher than in the Southern African region. This could be argued to be because of the presence of significant river basins in the central region such as the Congo River basin. Apart from the OND period (Figure 3h), which suggests the presence of an average amount of rainfall in the Southern African region (i.e., Namibia, Botswana, Lesotho, Swaziland, South Africa, Zimbabwe, Zambia, Malawi, and Mozambique), other seasons show very little rainfall in the region. However, the occurrence of significant rainfall can be seen from all seasonal periods in the central African region (Figure 3e-g). The mobile hydrological activity of the central African region is seen because of the surfeit of physical phenomena moderating its climate system. The Future Climate for Africa (FCFA) program which is implemented by IMPALA, FRACTAL and the Integrating Hydro-climate science into policy decisions for climate-resilient infrastructure and livelihoods in East Africa (HyCRISTAL) models, for instance, improves continent-wide prediction of climate on a 5 to 40-year scale, and provides insights and more reliable information about climate processes and extremes in the central and southern African region. Given that a vast majority of the populace in the SCA relies on subsistence agriculture for income and survival (e.g., USAID, 2020) [96], it becomes necessary to have accurate predictions of the climate and weather pattern of the region, as the livelihood of numerous people depend on it.
The SCA region is plagued with highly unstable weather patterns which causes a significant effect on industry (mining and tourism), energy (electricity generation and hydropower/ renewable energy), transport (aviation, shipping ports and other transport corridors), agriculture (drought) and hydrology (flooding). Regardless of the vast water resources that exist in the region, which includes rivers, lakes, groundwater systems, and wetlands, the region has a history of fluctuations in climate (e.g., Rodell et al., 2018) [22], specifically changes in rainfall amount, which were linked to the El-Nino southern oscillation phenomenon. These frequent climatic susceptibilities present a significant threat to the available water resources, as well as increase the frequency and intensity of climatic events (e.g., Kristen et al., 2007) [97]. On the spatial and temporal dynamics of rainfall over the study region (Figure 1a), the obtained rPCA cumulative variability accounts for 75.12% for the first six most dominant modes. Since the accuracy of the rPCs are directly proportional to the area of the spatial domain, an improvement of the variances obtained from the six most dominant modes was determined by masking the river basins of interest (Figure 1,2) in order to enhance scaling. In Figure 6, the utmost loadings from the first rE-OFs, characterizing 19.51% of the total variability are concentrated around Angola, Zambia, and Congo, which also houses the Congo River basin (Figure 6a). This rotated orthogonal mode depicts sub-regions with significantly strong annual rainfall variability. The second rEOF mode, which explains 12.16% of the total variability, represents multi-annual signals caused by bimodal rainfall patterns and is dominant around Gabon, parts of Congo and Rwanda and the Central African Republic. The third rEOF mode characterizes 10.51% of the total variability describes the multi-annual rainfall patterns which converge over the Central African Republic, Congo and Gabon. The fourth rEOF mode, explains 11.34% of the total variability depicting annual rainfall patterns dominant around Malawi, Zambia, Botswana, South Africa, Lesotho, Swaziland, Zimbabwe, Mozambique, and parts of Angola. The fifth and sixth rotated orthogonal modes and their corresponding rPCs explain 6.23% and 15.37% of the total variability, respectively, both denoting annual patterns of rainfall over the SCA region. The mobility of these precipitation patterns is majorly due to climate teleconnections, circulation features, annual rainfall cycles, as well as the impacts of the large river basins existent in the region (Ndehedehe et al., 2016a) [14].These rainfall structures and patterns over Central Africa suggest the dominance of annual variability in the FCFA models. This annual variability can be connected to the existent changes in the rainfall pattern of the third quarter of the year (Nicholson, 2014) [98].
A considerable relationship is observed between the dominant patterns of the spatially varying rainfall ( Figure 6) and the seasonal distribution of rainfall (Figure 3e-g). Additionally, the intense rEOF loadings (Figure 7a-c) observed between −20 • N and 10 • N (the Congo River basin)strongly suggest the significant changes in annual and multi-annual precipitation patterns that occurs in the Congo River basin area which spans across Angola, Zambia, Rwanda, Central Africa Republic, Gabon and the Democratic Republic of Congo during the April to December (Figure 6a-d) period. The temporal evolvement of the corresponding leading rEOF loadings illustrates that between the period of August to November, the maximum peaks of precipitation are observed. Additionally, this is in tandem with Figure 3e-h, which shows that the strong spatial patterns of the observed rainfall are usually between March-May and September-November. Some studies have shown that this rainfall mode, found in Figure 7a  Precipitation is extremely variable with a strong presence over all seasonal time scales in the humid central African region (Figure 6a-c). This therefore warranted the inclusion of all monthly precipitation information due to their diversity in local climates.
Basically, the leading rainfall mode observed in the central African region (Figure.  7a-c), which has strong rEOF loadings over SCA sufficiently describes the inter-annual variability of rainfall in the region. This is coherent with previous studies and is generally accepted by decision-makers for planning purposes (Rodell et al., 2018; Mohammed and David, 2019) [19,22]. Between 1999 and 2002 rainfall averaged 4% above the usual Precipitation is extremely variable with a strong presence over all seasonal time scales in the humid central African region (Figure 6a-c). This therefore warranted the inclusion of all monthly precipitation information due to their diversity in local climates.
Basically, the leading rainfall mode observed in the central African region (Figure 7a-c), which has strong rEOF loadings over SCA sufficiently describes the inter-annual variability of rainfall in the region. This is coherent with previous studies and is generally accepted by decision-makers for planning purposes (Rodell et al., 2018; Mohammed and David, 2019) [19,22]. Between 1999 and 2002 rainfall averaged 4% above the usual reading, while it averaged 1% below the usual reading during the remaining period of the GRACE obtained readings in the central African region (Rodell et al., 2018) [22]. Between 2014 and 2015, a dry precipitation pattern was observed over this region (Figure 8a-c) thereby implicating large-scale climatic oscillation as the ultimate driver. This, however, corroborates the report made by Rodell et al., 2018 [22]. The dominant rainfall modes (Figures 6a and 7a) represent the principal climatological properties and changes in the SCA region. For example, some regions of SCA have experienced episodes of drought events (Figure 7e,f, Ujeneva and Abiodun, 2015) [52], and are marked out in the amplitudes of their temporal variations (Figure 7). Based on the observed peak amplitudes (Figure 7), extreme wet and dry years which are results of strong impacts of El Nino and La Nina cycles (Ujeneva and Abiodun, 2015) [52] can also be recognized in the temporal evolutions of the analyzed rainfall modes (Figure 7a-d). Therefore, the annual and multi-annual variability based on the bi-modal rainfall patterns can be said to be major drivers of land water storage in the SCA region.

Trends in TWS and Rainfall
To spatially average monthly values over the SCA region, it was pertinent to estimate the seasonal cycle of rainfall and changes in TWS (Figure 3). While the averaged changes in TWS show maximum seasonal variations in the January-February-March periods (Figure 2a), the averaged values for the precipitation seasonal windows display maximum rainfall between the mid-August (Figure 2g) to mid-November periods also (Figure 2h). This evidently shows that in as much as the time lag between TWS and precipitation in most regions may be less than one month, there still exists an approximately four to five-month lag of precipitation behind TWS (Figure 8a,b) in most regions. The time response of TWS to rainfall is further verified through the temporal analysis of the first two dominant TWS and precipitation rotated orthogonal modes which clearly indicate little or no lag exists in Figures 4a, 5a, 6a and 7a, but an approximately four-to-five-month lag exists in Figures 4b, 5b, 6b and 7b. Most importantly, the rainfall seasonal periods are homogenous with the classification of major FCFA models outlined in recent studies (Mamane and Malam, 2015) [99]. A significant observation made from the TWS/rainfall variability index is the correspondence of dry and wet seasons. This therefore suggests precipitation as a significant correspondent to the hydrological flux in the SCA region (Rodell et al., 2018) [22]. Additionally, it purports observed change in precipitation patterns in the SCA region as likely to be the most substantial incentive for surface mass variations, hence the

Trends in TWS and Rainfall
To spatially average monthly values over the SCA region, it was pertinent to estimate the seasonal cycle of rainfall and changes in TWS (Figure 3). While the averaged changes in TWS show maximum seasonal variations in the January-February-March periods (Figure 2a), the averaged values for the precipitation seasonal windows display maximum rainfall between the mid-August (Figure 2g) to mid-November periods also (Figure 2h). This evidently shows that in as much as the time lag between TWS and precipitation in most regions may be less than one month, there still exists an approximately four to five-month lag of precipitation behind TWS (Figure 8a,b) in most regions. The time response of TWS to rainfall is further verified through the temporal analysis of the first two dominant TWS and precipitation rotated orthogonal modes which clearly indicate little or no lag exists in Figures 4a, 5a, 6a and 7a, but an approximately four-to-five-month lag exists in Figures 4b, 5b, 6b and 7b. Most importantly, the rainfall seasonal periods are homogenous with the classification of major FCFA models outlined in recent studies (Mamane and Malam, 2015) [99]. A significant observation made from the TWS/rainfall variability index is the correspondence of dry and wet seasons. This therefore suggests precipitation as a significant correspondent to the hydrological flux in the SCA region (Rodell et al., 2018) [22]. Additionally, it purports observed change in precipitation patterns in the SCA region as likely to be the most substantial incentive for surface mass variations, hence the dominant patterns perceived in the variations of TWS.

Hydrological Response of TWS to Climate Data
The lack of in-situ hydrological data to assist in adequate monitoring of variations in TWS have resulted in the use of several drought indexes to facilitate the observation of water availability and hydrological specifications (e.g., Laux, 2009;Heim, 2002;McKee et al., 1993McKee et al., , 1995 [100][101][102][103]. Regardless of inherent uncertainties, the use of these indexes combined with hydrological models has shown relatively good performance in water storage predictions (Pedinotti et al., 2012;Xie et al., 2012) [38,39]. Nevertheless, the model output is still plagued with data insufficiency needed for calibration and evaluation purposes, thus hampering data integrity for accurate predictions. This has therefore, resulted inseveral studies attempting to merge remotely sensed observations and model outputs in order to enhance the estimation of TWS variation, as a means to bypass the issue of sparse model data (e.g., Boone et al., 2009;Schuol et al., 2006;Wagner et al., 2009;Leblanc et al., 2006) [104][105][106][107]. This study uses satellite observations of precipitation and ET to determine the climatic conditions of the SCA region, and its impact on the TWS fluctuations of its prevalent five major river basins.

TWS and Links to Evapotranspiration Anomalies (ETa) and Precipitation
TWS as a fundamental component of the global water cycle is imperative for climate studies, water resources, agriculture, and ecosystem experiments. This therefore encourages an improved understanding of the potential relationship between the Terrestrial Water Storage Change (TWSC) and varied atmospheric phenomena, which could result in an enhanced parameterization in climate models. On the other hand, evapotranspiration is an important but less understood component of the terrestrial watercycle. However, its relevance in linking the land and atmospheric branches of the hydrological cycle by dictating moisture and heat exchange between them makes it a relevant factor in climate studies (Ndehedehe et al., 2018b) [17]. ET attempts to utilize a strong influence on the patterns of atmospheric water vapor transfer.
From this study, we are able to understand the direct and/or indirect impact of ETa on TWS fluctuations in the SCA region. Inasmuch as an average correlation exist between the two components (Figure 9a), we find the relationships between Eta ( Figure 10) and TWS from the rPCA plots. The improved algorithm developed by Mu et al. (2011) [84] provides insight into the necessity of segmenting ET components given their temporal varied response to climatic episodes. This improved algorithm helps to isolate and refine episodes such as water, snow, and ice evaporation (WSI ET), soil evaporation (ET) and even vegetation transpiration (Veg. Trans) to include both day and nighttime readings. The highest ET rate is experienced in the Congo River basin, recording an average of 222.9 mm between 2002 and 2017 (Table 2, Figure 11). The same region also experienced the highest average precipitation of 157.305 mm from 2002 to 2017 (Table 3). The strong relationship between the ET and precipitation in the region is evident in Figure 12a. However, the extremely high precipitation prevalent in this region (Figures 3e-h and 12a, Table 2) has resulted in a steady decline of WSI ET based on the PE-based trend in Figure 10a (Table 4) and stabilizes between 2007 and 2010 (Figure 1c,d). In total, the Congo River basin gained TWS of 30.00 mm between the 2002 and 2017 period. This same region however experienced the strongest variability in TWS when the rPCA technique was applied (Figure 4a). The positive TWS trend observed around the Congo River basin in Figures 4a and 4c can be said to be a result of the high precipitation patterns in the corresponding years ( Figure 12, Table 3).     Figure 11. Mean evapotranspiration of the five major river basins in SCA from July 2002 to December 2017. In order to improve data spread in the basins, the gamma stretch technique is applied. Figure 11. Mean evapotranspiration of the five major river basins in SCA from July 2002 to December 2017. In order to improve data spread in the basins, the gamma stretch technique is applied.  As the amount of water leaving the region (ET) was higher than the amount of water coming in (precipitation), the positive trend in TWS can therefore be attributed to the extreme humid nature of the Congo River basin (Figure 1b,c). The relatively high ETa rate can be attached to the elevated rate of solar radiation absorption and disposition, low vegetation transpiration due to the ongoing deforestation, and the optimum distribution and movement of moisture within the soil of the region (Figure 10a, White and Toumi, 2014;Samba and Nganga, 2012) [108,109]. Evapotranspiration affects energy balance which means that sunlight also plays a key role in its ultimate effects. This same result  As the amount of water leaving the region (ET) was higher than the amount of water coming in (precipitation), the positive trend in TWS can therefore be attributed to the extreme humid nature of the Congo River basin (Figure 1b,c). The relatively high ETa rate can be attached to the elevated rate of solar radiation absorption and disposition, low vegetation transpiration due to the ongoing deforestation, and the optimum distribution and movement of moisture within the soil of the region (Figure 10a, White and Toumi, 2014;Samba and Nganga, 2012) [108,109]. Evapotranspiration affects energy balance which means that sunlight also plays a key role in its ultimate effects. This same result was seen in the rest of the river basins, each with a corresponding TWS variability. For example, the Zambezi River basin maintained a relatively high ET rate with an average of 148.188 mm/year between 2002 and 2017 (Table 2, Figure 11), and an average precipitation rate of 81.161 mm/year between the same period ( Table 3) (Table 3). This indicates a strong rate of solar radiation absorption and disposition throughout the region. A lot of the water gained in this basin is lost through ET (Figure 12c). In total, the Okavango basin gained a TWS of 0.254 mm between the 2002 and 2017 period (Figure 4e). Among all the basins used in this study, the Limpopo basin recorded the highest amount of temporal TWS loss (Table 4). It maintained an average ET and precipitation rate of about 86.58 mm and 44.06 mm, respectively. The ET rate is almost twice the precipitation rate of this region. It therefore follows that without runoff, precipitation and groundwater, this basin will cease to exist in a few years. Figure 12d shows the direct relationship between these two climatic drivers. The TWS of the Limpopo basin is represented as a bi-modal multi-annual variation with a negative trend from 2002 to 2012 (Figure 1h). Nevertheless, a TWS of 0.084 mm was gained during this period. This same case is similar in the Okavango basin which asides incoming precipitation rate, can be argued to have sufficient groundwater to augment the increasing trend in WSI ET (Figure 12c (Table 4).
Due to the complexities that surround ET relationships to TWS fluctuations, few studies have attempted to link them. This study directly shows the impact of ET rate on TWS and suggests that a high ET rate complemented by increased WSI ET, Soil ET and vegetation transpiration can evidently affect TWS anomalies as seen in the Congo River basin (Tables 2-4, Figure 1c,d, Figures 4a, 10a, 12a and A1).

Influence of Evapotranspiration Anomalies on Rainfall Frequency during the 2002-2017 Period
The presence of surface waters (e.g., rivers, lakes, and reservoirs), soil moisture, and groundwater in the south and central African region (especially over the major river basins used in this study, Figure 1) are major constituents of the GRACE derived TWS, which are driven by precipitation and evapotranspiration rate. Precipitation and evapotranspiration represent the volume of inbound and outgoing water in the hydrological budget which is considered to have a net-zero loss. A low rate of evaporation in the cloud can lead to a sudden drop in evapotranspiration. This in turn results in low humidity which directly translates to a reduction in the prevalence of rainfall in that region (Table 1). Without ET, circulation and rainfall patterns would be quite different. From Figure 8b, a high correlation coefficient between TWS and rainfall is observed (Table 1). Rainfall in the Congo River basin is distributed bi-modally in the annual cycle, with two wet seasons March-May and September-October, and two dry seasons June-August and December-February (Table 1, Figures 3e-h and 12a). During the first wet season (MAM) in the Congo River basin (see Figure A2a), an estimated ETa rate of 271.33 mm is recorded in the heart of the basin (Table 1), while the Northern and Eastern parts recorded a relatively lowerETa rate of 239.75 mm. The effect of the variation in ETa rate can be seen in Figure 3e,h, where the domination in rainfall pattern is shown in the Northern and South eastern part of the river basin. Between the end of summer and the beginning of the winter period, a significant reduction in ETa is recorded over the basin (Figure A2b), due to the drop in ETa drivers during this period. This therefore resulted in a significant shift in rainfall pattern (Table 1) represented in Figure 3f.
This implies that the break between the summer and winter periods generates significant fluctuations that affect ETa and rainfall patterns. Between mid-June to mid-November, a consistent average of 231.06 mm was recorded in the upper part of the basin, while an average of 204.86 mm was recorded in the lower part of the basin. A reflection of this ETa fluctuation is seen in Figure 3f,g where a greater wetting period is observed in the upper basin part. The Zambezi River basin, on the other hand, recorded a more stable ETa throughout the summer (Figures A1d and A2e-h), including the break between summer and winter. The highest ETa reading was recorded between December-February, with an average of 236.73 mm/year in the northern and eastern parts of the basin, and 198.4 mm throughout the remaining basin ( Figure 11b, Table 1). In the Limpopo River basin, the strongest ETa fluctuation was recorded between December and February which is the hottest period in the region (Table 1, Figure A2l). It maintained an average ETa rate of 157.21 mm/year in the upper region, and 133.58 mm/year in the lower basin region (Table 1, Figure 11d). However, between June and August ( Figure A2j), the lowest ETa rate of about 43.49 mm was recorded throughout the basin ( Table 1). The Orange River basin generally maintained a low ETa rate due to its relatively small size. However, it recorded the lowest ETa rate between mid-February and mid-June ( Figure A2n) of about 26.76 mm/year (Table 1) during the wet season. Contrary to this observation, an ETa rate of 37.90 mm/year was recorded between mid-august and mid-December ( Figure A2p). The other seasons around the basin maintained an average ETa rate of 31.20 mm/year (Table 1).
In an observation made in this study, the ETa fluctuations around the seasons were not consistent as we can see that the lowest ETa rate for the Congo, Okavango and Limpopo River basins was recorded between June and August ( Figure 12, Table 1), whereas the lowest ETa rate for the Zambezi and Orange River basin varied between September and November and between July and August ( Figure 11, Table 1).
The variations in the seasonal ETa rates of these river basins ( Figure A2a-p) signify the presence of contrasting factors that influence these changes (e.g., Figure 10). The Congo River basin maintained the highest ETa fluctuations, while the Orange River basin has the lowest ETa fluctuations over the seasons. As much as size plays a role in the region's climatic balance, other existing factors including the presence of a relatively strong vegetation index, the region's average surface temperature and the rate of transpiration and latent energy flow each contributing significantly to the varied fluctuations determined here.

Conclusions
Evapotranspiration is an important component of the hydrological cycle over land. It returns water to the atmosphere and supplies run off. Global Climate Models (GCM) now attempt to integrate climatic details through model-coupling of the land and ocean surface atmospheric models. For the atmosphere, these models provide solutions to the 3D fluid conservation equation for heat and momentum, as driven by atmospheric radiation and latent heat release. The water returned by ET to the atmosphere helps drive circulation patterns both in models and reality. Therefore, understanding the usefulness of ET in undertaking regional and global climate investigations is crucial. For this to be achieved, we utilized several algorithms to easily identify the cause analysis for the SCA region. Given the increasing interest in the concept that deforestation directly impedes access to freshwater, this study corroborated this claim as observed in the Congo River basin. However, the rate of deforestation in the Central African region suggests an imminent negative effect on freshwater and the regional ecosystem of the residents of this region. This study employed the use of the rotated principal component analysis and multi-linear regression analysis in a bid to understand the TWS changes of the Southern and Central African region using five major river basins (Congo, Zambezi, Limpopo, Orange, and Okavango) within the study region. Results from our analysis show that: (i) High annual variability of GRACE-derived surface mass variations is observed in the Democratic Republic of Congo, Central African Republic, Zambia, Malawi, Republic of the Congo, and Gabon. This suggests the huge impact of the Congo basin in regulating TWS fluctuations in this region. Additionally, the Congo River basin, which is largely dominated by multi-annual signals was identified from the rPCA result. (ii) Precipitation over the region is dominated by annual and bi-modal multi-annual patterns influenced by circulation features, vegetation height, and climate teleconnections. Even though an overall TWS loss in the studied river basins was not observed, episodes of estimated losses in TWS in the northern part of the Central African region which is consistent with past studies are noted. This loss in TWS seems to be caused by natural inter-annual variability, although past studies insinuated that the surface runoff rate is enhanced by deforestation. This means that an increase in WSI ET is expected. The region also experienced an increase in precipitation rate between mid-2013 and mid-2016. Due to the consistency in the postulated negative correlation between TWS in the Congo and Amazon River basins, we therefore claim large-scale climatic oscillation as the ultimate driver of the TWS fluctuation in this region amidst other factors. (iii) The Limpopo River basin experienced a negative trend of about −4.6 ± 3.2 mm in TWS between 2006 and 2012. The 6% below average reduction in rainfall rate is also a contributor to this situation. The center of the region is dominated by Lake Malawi which maintained a correlation of 0.88 with the regional TWS. The lake declined at an estimated rate of 77 mm/year during this period, thus contributing to the observed trend in TWS. Therefore, we attribute the cause of this apparent trend to natural and climatic variability, although a 7% decrease in precipitation is predicted during this century. Additionally, the Orange River basins had an average TWS gain of 0.34 ± 9.1 mm/year regardless of the high fluctuation rate observed in precipitation patterns. (v) Overall, rainfall leads TWS in approximately five months majorly in the central part of the study region with maximum correlation coefficients (r) ranging from 0.7 to 0.9. TWS in some parts of SCA however, show low and modest correlations with precipitation. In these hydrological regions, precipitation leads TWS with a phase lag ranging from two to five months and shows that besides rainfall, other primary drivers of variations in TWS exist, and are yet another pointer to the role of ETa in TWS fluctuations carried out in this study. Additionally, considering the average r 2 values and uncertainties in modeling precipitation over the humid part of the study area (major central Africa), advancing the multi-linear regression model for better climatology and freshwater prediction could mean inflation of the independent variables to include other germane physical processes. (vi) The ETa fluctuations observed in this study correlated more with rainfall than TWS because of the respective direct and indirect interactions existent between these components. The highest ET rate was observed during the dry period (for example, DJF) because of atmospheric radiation and latent heat release. However, we observed the anomaly spiking in intervals between the summer and winter period as shown in the appendix plots. This evidently shows that the average relationship between ETa and TWS can be sustainable as well as important in analyzing TWS fluctuation patterns during the period of high climatic variations. of the study region with maximum correlation coefficients (r) ranging from 0.7 to 0.9. TWS in some parts of SCA however, show low and modest correlations with precipitation. In these hydrological regions, precipitation leads TWS with a phase lag ranging from two to five months and shows that besides rainfall, other primary drivers of variations in TWS exist, and are yet another pointer to the role of ETa in TWS fluctuations carried out in this study. Additionally, considering the average r 2 values and uncertainties in modeling precipitation over the humid part of the study area (major central Africa), advancing the multi-linear regression model for better climatology and freshwater prediction could mean inflation of the independent variables to include other germane physical processes. (vi) The ETa fluctuations observed in this study correlated more with rainfall than TWS because of the respective direct and indirect interactions existent between these components. The highest ET rate was observed during the dry period (for example, DJF) because of atmospheric radiation and latent heat release. However, we observed the anomaly spiking in intervals between the summer and winter period as shown in the appendix plots. This evidently shows that the average relationship between ETa and TWS can be sustainable as well as important in analyzing TWS fluctuation patterns during the period of high climatic variations.