Rainfall Trend and Its Relationship with Normalized Di ﬀ erence Vegetation Index in a Restored Semi-Arid Wetland of South Africa

: Clearance of terrestrial wetland vegetation and rainfall variations a ﬀ ect biodiversity. The rainfall trend–NDVI (Normalized Di ﬀ erence Vegetation Index) relationship was examined to assess the extent to which rainfall a ﬀ ects vegetation productivity within Nylsvley, Ramsar site in Limpopo Province, South Africa. Daily rainfall data measured from eight rainfall stations between 1950 and 2016 were used to generate seasonal and annual rainfall data. Mann-Kendall and quantile regression were applied to assess trends in rainfall data. NDVI was derived from satellite images from between 1984 and 2003 using Zonal statistics and correlated with rainfall of the same period to assess vegetation dynamics. Mann-Kendall and Sen’s slope estimator showed only one station had a signiﬁcant increasing rainfall trend annually and seasonally at p < 0.05, whereas all the other stations showed insigniﬁcant trends in both rainfall seasons. Quantile regression showed 50% and 62.5% of the stations had increasing annual and seasonal rainfall, respectively. Of the stations, 37.5% were statistically signiﬁcant at p < 0.05, indicating increasing and decreasing rainfall trends. These rainfall trends show that the rainfall of Nylsvley decreased between 1995 and 2003. The R 2 between rainfall and NDVI of Nylsvley is 55% indicating the inﬂuence of rainfall variability on vegetation productivity. The results underscore the impact of decadal rainfall patterns on wetland ecosystem change.


Introduction
Vegetation is critical for assessing the conditions of terrestrial wetland ecosystems impacted by global climate changes [1]. Wetland Vegetation plays crucial roles in habitat provision for wildlife, preventing soil loss, slowing water movement, and in the process of removing pollutants and purifying water [1][2][3]. Climate change warms the atmosphere, changes rainfall patterns and timing, leaving some parts of the world dryer and others wetter than before [4][5][6][7], and this impacts on wetland ecological processes [7]. The potential of climate change together with its effects on rainfall distribution is a threat to vegetation development. In the continent of Africa, rainfall amounts are predicted to decrease, especially in Sub-Saharan Africa, forcing already high temperatures to further increase [8][9][10]. The more rainfall received at a time implies more vegetation cover [11], especially in semi-arid regions where rainfall is erratic and is exacerbated by global climate change. Rainfall affects both the spatial and temporal patterns of water availability on the earth's surface and is used in climate science to trace the magnitude and extent of variability of climate change [12][13][14]. The emphasis is on investigating the We hypothesised that there is a relationship between the NDVI and rainfall patterns in Nylsvley. The findings provide additional insights for understanding climate variability within the Nyl River system (NRS), including the restored Nylsvley wetland. More light will also be shed on our understanding of the climate change impacts on vegetation growth in wetlands of semi-arid regions and assists with proper management and monitoring of this and similar wetlands for early warning against negative effects. Thus, an attempt was made to investigate the diurnal seasonal and annual rainfall in the NRS to determine the significant trends in rainfall from the year 1950 to 2016. It also explored the seasonal rainfall-NDVI relationship in Nylsvley between 1983 and 2003 as a contribution towards understanding the local-level impacts of climate variations in the restored wetland vegetation ecosystem. In addition, there is little evidence of studies that show long-term trend changes in rainfall and their ecological impacts on vegetation growth within the Nylsvley wetland. Restoration activities are costly hence the need to conduct such a study to fully understand climate change impacts within this wetland.

Study Area
Nylsvley floodplain (3975 hectares in size) is located to the east of Waterberg between Modimolle and Mookgophong in Limpopo Province of South Africa. Figure 1 shows the location of the Nylsvley Nature Reserve within the NRS [39,45]. The flood plain is located at 28 • 42.00 East and 24 • 39.00 South, within the Nyl River system which is managed by the Limpopo Conservation Authority. It was declared a Ramsar site of international conservation importance in 1998, and the wetland is largely recognized as it represents one of the most important avian conservation areas in South African wetlands [29,39,44]. The Nylsvley freshwater wetland is part of the Nyl River system situated in the semi-arid savannah region in Limpopo province of South Africa. This freshwater ecosystem is the largest inland floodplain covering an area of about 16,000 ha when fully inundated. The wetland's ability to provide groundwater recharge throughout the year makes it valuable and prompts the need for its conservation [45]. Rowberry et al. [52] note that the ecological status of the freshwater wetland has The floodplain houses over 100 aquatic bird species, 57 of which breed in this wetland. This makes it an important nesting and feeding site for migratory birds and several endangered species. There are 23 waterfowl bird species found in this wetland that are listed in the International Union for Conservation of Nature Red data book of threatened species [46,47]. The wetland, apart from its conservation value, is also agriculturally important as it provides animals with extensive grasslands for year-round grazing.
Nyslvley is in a semi-arid region characterized by hot-wet and warm-dry seasons with an average annual rainfall of 630 mm [35,45,46]. The NRS catchment is underlain by sedimentary rocks of the Waterberg, mainly composed of sandstone, grit, mudstone, siltstone, shale, and conglomerate [48]. The soils are predominantly acidic sands, which are highly variable, gravels and loams occurring at 1200 mm depth with altitude ranging from 1180 to 1508 m [48]. The Nyslvley flood plain, located in the middle of NRS is characterised by both well and poorly drained soils, black clay soils, and alluvial soils [1,2]. Nylsvley has a history of being grazed by domestic and non-domesticated animals.
There is highly seasonal rainfall, 95% of it occurs between November and March [35] with a mid-season dry spell during significant periods of growth. Winters (May to July) are dry and springs (September) have little to no rainfall (10 mm on average). The mean annual temperature is 18 • C with temperatures ranging from −6 • C (winter) to 39 • C (summer) [49]. The catchment area lies in~90% natural savanna bushveld and 10% croplands [49]. Landuse activities include agriculture (growing of crops and rearing of animals), game farming, and settlements. The Nyl River systems (NRS) consist of small tributaries that lie in the Waterberg plateau [50].
Different wetland vegetation species cover the park, including those in salt marshes, floodplains, and grasslands, being mostly Leersia hexandra, Eleocharis dregeana, and the grass Oryza longistaminata, which is restricted to Nylsvley only [44]. The most dominant species found in Nylsvley freshwater swamps are reeds (Phragmites australis). Rich underlying geological formations are responsible for supporting the unique biodiversity of plant and animal life [51].
Nylsvley Nature Reserve (NNR), which is a Ramsar site of international importance, has been under threat for the past two decades. Studies show that the water quality is poor and the population of aquatic bird species visiting NNR wetland has reduced drastically by 80-90% [28,30,36] The Nylsvley freshwater wetland is part of the Nyl River system situated in the semi-arid savannah region in Limpopo province of South Africa. This freshwater ecosystem is the largest inland floodplain covering an area of about 16,000 ha when fully inundated. The wetland's ability to provide groundwater recharge throughout the year makes it valuable and prompts the need for its conservation [45]. Rowberry et al. [52] note that the ecological status of the freshwater wetland has been threatened by the increasing developments within upstream catchments affecting the availability of water and vegetation in the floodplain. Sustainable conservation and management of wetlands require regular monitoring and assessments to check on the status and changes of vegetation in wetland zones [53]. Wetland assessment involves the inventory and mapping of vegetation species within the wetlands [54,55]. Identification, mapping, and assessing the condition of freshwater wetlands is of paramount importance for various management and planning activities at global, regional, and local scales [53,56,57]. This offers information on wetland vegetation types, their distribution and phenology, and acts as a good indicator in cases of assessing wetland deterioration [53,58,59].
Originally, the NNR was part of a privately-owned property that was active in agriculture. Due to the reserve's unique landscape and biodiversity, it was commissioned to the government in 1974, and the other portion, Vogelfontein was added in 1997. The nature reserve harbors South Africa's largest freshwater flood plain, a variety of birds, and many other faunal species. Since the year 2000, the improvement of ecological integrity has become the key objective with a clear shift from projects at a habitat scale to integrated projects of the entire ecosystem [60]. Simenstad et al. [61] view the 21st century as a period for fixing global environmental problems that were caused by human ignorance. In this century, the desire is to restore the disturbed areas and assess the progress of such restoration activities, as ecosystem functioning is impacted by alterations or changes in the landscapes [61]. It is vital to note that restoration plays a critical role in environmental management and policy development for monitoring purposes [62]. This helps in understating the benefits of most restoration efforts done today especially the development of vegetation in semi-arid regions under climate change projected period.

Data Availability
To assess and understand the trends in rainfall within the NRS, daily rainfall data measured from eight rainfall stations ranging from 27 to 66 years between the years 1950 and 2016 were studied. The selected long-term period provided a common data set for rainfall and allowed for reliable climatic variations in this study. The rainfall data set was obtained from the South African Weather Service (SAWS) and Agriculture Research Council (ARC), two and six stations, respectively. This data has the advantage of quality control with information on missing or unaudited data explained. However, it should be noted that climate time series data are affected by several non-climatic factors which include changes in instruments, human errors, or different observation practices [63]. Hence more rainfall stations were studied before relating vegetation of the wetland to Nylsvley wetland rainfall station to check on the reliability of this station data. Total daily rainfall values were used to come up with mean monthly, seasonal, and annual rainfall. Mean monthly rainfall was correlated with mean NDVI in Nylsvley from 1984 to 2003. Nylsvley rainfall station data was limited up to the 2003 period forcing the analysis to also end in the same period. Seasonal rainfall time series were selected to show changes that may not be depicted by total annual rainfall. Mazvimavi [64] highlighted that overtime changes in rainfall at smaller spatial scales have the potential to show climate change impacts on a localized scale.
NDVI data were derived from Landsat 5 images downloaded from USGS through the earth explorer website (https://earthexplorer.usgs.gov/) and have been geometrically and radiometrically corrected. The imagery at 30 m resolution with a temporal resolution of 16-days (row 170/077, UTM projection), acquired in the summer season for selected years between 1984 to 2003, were considered in this study, see Appendix A for the selected images. The selected sensor is the best available cloud-free image in the study area for the period. Images from the beginning of the season were selected as most were free from cloud cover (October). The winter vegetation coverage of semi-arid regions was mostly poor, hence we focused on summer seasons (October to March), selecting cloud-free images. NDVI values were thus averaged and utilized for analysis in this study. The detailed properties of NDVI are provided in many other studies including [8,26,65,66]. The research was done in semi-arid regions in the 1980s and the 1990s, demonstrating the seasonal relationship between NDVI and rainfall variations [67]. NDVI enabled us to examine ecosystem dynamics and long term trends in vegetation [68].

Data Analysis
To analyze trends in rainfall data, statistical methods have been developed and applied in studies over many years [42,63,69,70]. Most of these studies have used linear regression models, based on ordinary least-squares regression, which were used mostly in linear trend studies. These models have limitations surrounding providing information on the mean condition of indices, ignoring aspects of the lower or upper tails of index distributions that are more important than the mean trends in climate trend studies [71]. Odiyo et al. [42] highlighted that using more than one method in trend analysis has the advantage of improving results. Thus, to boost the accuracy of and confidence in the study results, the Mann-Kendall trend test (MK), Sen's slope estimator, and quantile regression (QR) were used to analyze annual and seasonal rainfall data of NRS and final trend station results were compared. Tareghian and Rasmussen [71] and Lee et al. [72] noted that the QR method has recently been introduced in climate-related studies, and it is new in rainfall application studies in Southern Africa [64,73].
The non-parametric MK method has been preferred over the parametric test method in analyzing data trends as it also estimates trend direction [69]. The Sen's slope method is used for estimating the Sustainability 2020, 12, 8919 6 of 22 slope of a linear trend. It has been widely used for determining the trend magnitude in hydrological and meteorological time series data [14,74,75]. Asfaw et al. [10] note that non-parametric statistical tests are suitable for data that is not normally distributed and with missing data which is a common phenomenon in hydrological time series data. Shapiro-Wilk Normality Test (R-Studio software) was used for testing the normality distribution of the data of all the rainfall stations. MK was used in detecting the significance of the trend in rainfall time series and it works well with data that has outliers. The MK test is based on the statistic S defined in the following formula: where N stands for the number of data points, x j and x i represent the values at time j and i(j > i), respectively, sgn stands for sign [14]. The Sen's slope estimator determined the magnitude of the trend in rainfall data. Pohlert [76] states that Sen's slope test computes the linear rate of change and intercept as follows: where d = slope, k = median X denotes the variable, n = number of data and i, j = indices. Sen's slope is therefore calculated as the median from all slopes: b = Median d k . These intercepts are computed for each time step t, and a is intercept as shown in the formula below: This formula computed the upper and lower confidence limits for Sen's slope.
The results of Sen's slope that had a negative sign represent the decreasing trend, while a positive sign showed the increasing trend [77]. Long-term trend analysis had an advantage, as it can reveal the starting point of the trend year and changes occurring over time, as well as abrupt detections of changes in time series. The Sen's slope estimator embedded within MK was used to determine the magnitude of trends in time series. MK, a non-parametric test, was ideal for analyzing hydrological and climatological data which were characterized by missing or non-normally distributed data as this was a common feature [78].
The QR is an extension of the least-squares regression method [63] and examines functional relations between variables for all portions of a probability distribution and estimates multiple rates of change (slopes) from the minimum to maximum responses [79,80]. It provides a complete picture of the relationships among variables when both the lower and upper quantiles are of interest in a study [81]. Changes in rainfall regimes may not equally affect all rainfall percentile values [82] and hence the need to examine changes in several quantiles in detecting trends among different rainfall distributions. In this study, it was used in the estimation of slopes of trends in quantiles 0.1, 0.2, 0.5, 0.7, and 0.9 of the total recorded daily rainfall values which were used to compute seasonal and annual rainfall for NRS. Quantiles 0.1 and 0.2 represented low rainfall, 0.5 and 0.7 median, and 0.9 indicated high rainfall events [79,80,83]. The first order QR model had X which was defined as covariant of Y and τ and X = x, as follows: where Y is a random variable of the total rainfall time series and θ has values ranging from 0 ≤ θ ≤ 1 with Y probability being equal or less than τ which signifies the θ th quantile. The β(θ) 1 is the coefficient, the β(θ) 0 is the intercept depending on the θ th quantile value and ξ signifies an error with exception of zero [63,64]. The QR estimates also the p-values, t-statistics, standard errors, and confidence intervals. The negative or positive trend detection that is different from zero indicates that the θ th quantile of Sustainability 2020, 12, 8919 7 of 22 rainfall is either decreasing or increasing signifying low or high rainfall, respectively. Each rainfall station was considered independently in this study. The QR estimates and MK results were compared through values of slope coefficients and p-values across quantiles [84]. The Statistical Analysis System (SAS) 'quantreg' package software (SAS v94), was run for QR analysis and MK for trend analysis. The SAS v94 software version can show upper and lower tail distributions that provide additional crucial information on changes in all parts of the distribution of a series [63]. A regression model with heterogeneous variance implies that there is not a single rate of change that characterizes variation in probability distributions. To focus only on analyzing changes on means may under or overestimate or even fail to distinguish real non-zero changes in heterogeneous distributions [85]. Trends in all quantiles were compared using p-values and slope coefficients to investigate their similarities and differences in explaining different quantile regression slopes.
The acquired Landsat 5 images were imported in the ArcGIS 10.7 environment, for analysis. The study area was extracted by clipping the image with the Nylsvley boundary in ArcGIS 10.7. All downloaded images were refined using Erdas 2018 and ArcGIS 10.7 tools and the images were clipped to the study area shapefile boundary before extracting the mean NDVI values. The study area was an inland nature reserve with limited air pollution, with no specific intra-class vegetation analysis needing to be done, and atmospheric corrections for haze was not done as cloud-free day images were used. The NDVI as described in [26,86] was used for developing NDVI maps using the spatial analysis tool within ArcGIS 10.7 software. NDVI calculation uses the Near InfraRed (NIR) and the RED bands as vegetation strongly reflect in the NIR. The NDVI formula is as follows: NDVI values range between −1 to 1, with positive values indicating high amounts of vegetation and negative values responding to poorly or sparsely vegetated areas. Using the image analysis toolbar in ArcGIS 10.7, NDVI maps showing areas of healthy and not healthy vegetation were produced. Mean NDVI values were created using Zonal statistics as a table tool in ArcGIS 10.7. The correlation coefficient which measures the strength of the linear relationship of two variables was computed as postulated by Onema and Taigbenu [18]. The magnitude of the correlation between NDVI and rainfall in Nylsvley was measured in the period 1984 and 2003, with statistical significance at a confidence level of 95% (p = 0.05). The R-squared (R 2 ) revealed the goodness of fit in a linear regression model was used in the analysis of Rainfall-NDVI data. The correlation formula used is as follows: where r = correlation coefficient of variables rainfall and NDVI, x i = i th of the first variable (rainfall), y i = i th of the second variable (NDVI), x = mean for rainfall, y = mean for NDVI. The extracted NDVI values were exported to Excel for the creation of graphs with output NDVI and rainfall trend patterns. The study hypothesized that there is a trend in both NDVI and rainfall time series.

Rainfall Variability Across Stations
The eight stations had rainfall data observations ranging between 27 to 66 years with Groenfontein and Modimolle stations having least and most observations, respectively (Table 1). Positive and negative Sen's slopes indicated an increase and decrease in annual rainfall, respectively. The analysis of the annual rainfall showed four (50%) of the eight stations (Mookgophong, Modimolle, Baffelspoort, and Elandspoort) with increasing rainfall trends and the other four (50%) (Groenfontein, Moorddrift, Nylsvley, and Polokwane) with decreasing rainfall trends over the years under study. Seasonally, an additional station, Groenfontein, showed an increasing rainfall trend over the years, making 62.5% of the stations to be with increasing rainfall and 3 stations (37.5%) with decreasing rainfall (Table 1). The QR model trend analysis highlighted annual and seasonal rainfall results ( Table 2). The negative slope coefficients showed decreasing trends in all quantiles of either annual or seasonal rainfall. Similarly, positive slope coefficients showed increasing rainfall trends in both seasons studied. Polokwane was the only rainfall station showing decreasing seasonal rainfall trends in all quantiles. Table 2 clearly shows the stations with increasing or decreasing annual and seasonal rainfall trends in respective quantiles in NRS.
The results thus show that 37.5% of the rainfall stations were statistically significant in specific quantiles with either increase or decrease in annual or seasonal rainfall. Nylsvley station showed a decrease in all and most quantiles for annual and seasonal rainfall, respectively, whilst Elandspoort showed an increase in most and all quantiles for annual and seasonal rainfall, respectively. Groenfontein showed a decrease in most quantiles of annual rainfall and middle and upper quantiles of seasonal rainfall whilst lower quantiles (0.1, 0.2) showed increasing seasonal rainfall. The mean and median trends as described by Fan and Chen [63], reflect long-term trends in either lower or upper tails of the conditional distributions of the annual and seasonal rainfall. The trends in the median are like those in the mean, unlike the lower and upper tails that show different trends from the mean slopes (Table 1 for slope distribution of the plotted stations). The annual and seasonal rainfall of the eight rainfall stations analysed between 1950 and 2016 are presented in Figures 2 and 3.
The spatial distribution of stations with a statistically significant change in rainfall are shown in Figure 4. Elandspoort rainfall station, upstream of the Nylsvley Nature Reserve, had a statistically significant increase in rainfall (Figure 4). Nylsvley station, located inside the Nature Reserve, had reported a decrease in rainfall as revealed by rainfall data recorded over 53 years in this study. Whilst Groenfontein station, which had statistically significant decreasing annual rainfall, is located upstream of the Nature Reserve.  The spatial distribution of stations with a statistically significant change in rainfall are shown in Figure 4. Elandspoort rainfall station, upstream of the Nylsvley Nature Reserve, had a statistically significant increase in rainfall (Figure 4). Nylsvley station, located inside the Nature Reserve, had reported a decrease in rainfall as revealed by rainfall data recorded over 53 years in this study. Whilst Groenfontein station, which had statistically significant decreasing annual rainfall, is located upstream of the Nature Reserve.

Validation of Rainfall Results
Mann-Kendall trend test showed only one station (Elandspoort), showing a significant increase in both annual and seasonal rainfall at p < 0.05, Table 1. QR results, with p-values that were greater than α (0.05) showed non-statistically significant increasing or decreasing rainfall trends in the stations studied. Elandspoort station was statistically significant with increasing annual and seasonal rainfall at median quantiles, 0.5 and 0.7. The statistically significant stations at the 5% significance level varied per quantile. For instance, Groenfontein station was statistically significant in lower quantile (0.2) while Nylsvley station was also statistically significant in lower quantiles (0.1 and 0.2) whilst Elandspoort was statistically significant in median quantiles (0.5 and 0.7) for annual rainfall trends. Statistically significant stations are shown in Table 2 in bold and Figure 5 shows QR slopes of these statistically significant stations.

Validation of Rainfall Results
Mann-Kendall trend test showed only one station (Elandspoort), showing a significant increase in both annual and seasonal rainfall at p < 0.05, Table 1. QR results, with p-values that were greater than α (0.05) showed non-statistically significant increasing or decreasing rainfall trends in the stations studied. Elandspoort station was statistically significant with increasing annual and seasonal rainfall at median quantiles, 0.5 and 0.7. The statistically significant stations at the 5% significance level varied per quantile. For instance, Groenfontein station was statistically significant in lower quantile (0.2) while Nylsvley station was also statistically significant in lower quantiles (0.1 and 0.2) whilst Elandspoort was statistically significant in median quantiles (0.5 and 0.7) for annual rainfall trends. Statistically significant stations are shown in Table 2 in bold and Figure 5 shows QR slopes of these statistically significant stations.

Mean Monthly Rainfall and NDVI Results
The results of the average monthly summer rainfall and mean NDVI are as shown in Figures 6  and 8. The average summer season rainfall from 9 different months of separate years between 1984 and 2003 explains that there was a positive significant relationship between rainfall and NDVI in

Mean Monthly Rainfall and NDVI Results
The results of the average monthly summer rainfall and mean NDVI are as shown in Figure 6 and Figure 8. The average summer season rainfall from 9 different months of separate years between 1984 and 2003 explains that there was a positive significant relationship between rainfall and NDVI in Nylsvley shown by, correlation of 0.74 at p < 0.02 which was less than alpha of 0.05. The regression, R 2 explains that 55% of the data fits well on the model, shown in Figure 7. Onema and Taigbenu [18] highlighted that correlation measures the strength of two variables in a linear relationship and in this study, the strength of NDVI and rainfall was measured using rainfall and NDVI data for 9 years.

Mean Monthly Rainfall and NDVI Results
The results of the average monthly summer rainfall and mean NDVI are as shown in Figures 6  and 8. The average summer season rainfall from 9 different months of separate years between 1984 and 2003 explains that there was a positive significant relationship between rainfall and NDVI in Nylsvley shown by, correlation of 0.74 at p < 0.02 which was less than alpha of 0.05. The regression, R 2 explains that 55% of the data fits well on the model, shown in Figure 7. Onema and Taigbenu [18] highlighted that correlation measures the strength of two variables in a linear relationship and in this study, the strength of NDVI and rainfall was measured using rainfall and NDVI data for 9 years.   Figure 8 shows the trend of mean NDVI and average monthly rainfall in the summer season in Nylsvley wetland. Examining the relationship between NDVI and summer rainfall as shown in Figure 8, indicated that there was a general agreement between rainfall amount and availability of vegetation. Where rainfall increased, NDVI tended to increase too, see Figure 8, between the 1986 and 1991 period. The years 1993 and 1995 saw a decrease in both rainfall and NDVI. From the period 1997 as rainfall decreased, NDVI continued to drop despite rainfall increases as of 1998. This shows  Figure 8 shows the trend of mean NDVI and average monthly rainfall in the summer season in Nylsvley wetland. Examining the relationship between NDVI and summer rainfall as shown in Figure 8, indicated that there was a general agreement between rainfall amount and availability of vegetation. Where rainfall increased, NDVI tended to increase too, see Figure 8, between the 1986 and 1991 period. The years 1993 and 1995 saw a decrease in both rainfall and NDVI. From the period 1997 as rainfall decreased, NDVI continued to drop despite rainfall increases as of 1998. This shows that vegetation takes time to recover and regrow. By 2003, both rainfall and NDVI showed a downward trend, Figure 8. The mean monthly rainfall and mean NDVI values were high mostly in the months of November-January, indicating the strong relationship between the two. The beginning of season months, such as October, received lower rainfall and had low vegetation indicated by low NDVI.  Figure 8 shows the trend of mean NDVI and average monthly rainfall in the summer season in Nylsvley wetland. Examining the relationship between NDVI and summer rainfall as shown in Figure 8, indicated that there was a general agreement between rainfall amount and availability of vegetation. Where rainfall increased, NDVI tended to increase too, see Figure 8, between the 1986 and 1991 period. The years 1993 and 1995 saw a decrease in both rainfall and NDVI. From the period 1997 as rainfall decreased, NDVI continued to drop despite rainfall increases as of 1998. This shows that vegetation takes time to recover and regrow. By 2003, both rainfall and NDVI showed a downward trend, Figure 8. The mean monthly rainfall and mean NDVI values were high mostly in the months of November-January, indicating the strong relationship between the two. The beginning of season months, such as October, received lower rainfall and had low vegetation indicated by low NDVI.

Discussion
This study used the MK and Sen's slope, non-parametric trend statistic methods to determine the trend significance and magnitude, ideal for the comparison of total annual and seasonal rainfall variations within the Limpopo Province [87]. The study observed that both MK and Sen's slope methods gave similar trends for the rainfall series of all the eight stations studied (Table 1). From annual and seasonal analyses, four and five of the eight stations showed increases in annual and seasonal rainfall, respectively, in the NRS. Using QR, the results show that about half of rainfall stations showed decreasing annual rainfall trends in most quantiles while 37.5% of rainfall stations showed decreasing seasonal rainfall trends. Thus, most of the rainfall stations (62.5%) showed increasing seasonal rainfall trends in most of the quantiles. There are two rainfall stations (Moorddrift and Nylsvley), with constantly decreasing annual and seasonal rainfall trends on all and most of the quantiles, respectively (Table 2).
New et al. [27] in their study of precipitation from 1961-2000 highlighted that few stations in 14 African nations studied showed decreases in rainfall; however, the depicted trends by these authors were not statistically significant. This agrees with the studies by Hulme et al. [88] and Kruger [87], who reported that most stations studied in South Africa between 1901-1995 showed no significance in trends. This study, even though it was done at a localised scale and with data length varying (27-66 years), agrees with previous studies that an equal number of stations at the annual level, and fewer stations at the seasonal, in the region showed non-statistically significant decreases in annual and/or seasonal rainfall, respectively [87]. The results also show the significant increasing trends in both annual and seasonal rainfall at the Elandspoort station, whilst the rest of the stations with increasing trends have statistically non-significant increasing trends.
The decreases and increases in trends shown in the data could be because of increases and decreases in the extreme and inter-annual inconsistency of rainfall within some years in some provinces of South Africa [87]. This can explain why Elandsport and Groenfontein stations, both upstream of Nylsvley (Figure 2), have increasing or decreasing rainfall, respectively. Odiyo et al. [42] investigated variability in rainfall using long-term data (1931/1932-2005/2006) of the Luvuvhu River Catchment (LRC) of South Africa. Two of the six studied rainfall stations showed decreasing trends in 5-and 10-years mean rainfall with decadal fluctuation in a 10-year mean of daily rainfall. Their study concluded that the variability of rainfall has increased in LRC for over 86 years. Various other studies on long-term rainfall changes in South Africa include studies on analyses of rainfall data from 157 stations between 1880 and 1972 [89]. A study by Nicholson [90] reveals no trends in mean annual rainfall over Southern Africa between 1900-1970 and 1961-1990 periods. Kruger [87] reports significant decreases in annual precipitation in some parts of South Africa that include Northern Limpopo, Western KwaZulu Natal, North Eastern Free State, and the South-Eastern regions of Eastern Cape, and increases in rainfall were projected in Northern, Western, and Eastern Cape Provinces. Lynch [91] notes gradual rainfall increases in the Potchefstroom area from 1925 to 1998 whilst Warburton and Schulze [92] report a decrease in median annual rainfall in the Limpopo Province over the latter half of the 20th century.
All rainfall stations located in the upper section of the NRS, namely, Mookgophong, Modimolle, Elandspoort, and Buffelspoort have generally received increased rainfall in both seasons that were studied. In the middle section of NRS, two stations, Nylsvley and Polokwane, received decreased rainfall both annually and seasonally. The decrease in rainfall in Nylsvley over the years contributes to less water available in the Nylsvley floodplain. Moorddrift station situated in the lower section of the NRS has also faced decreased rainfall both annually and seasonally under the study (Figures 2 and 3). These differences can be explained by either natural climatic factors or differences in microclimates [93,94]. The main causes of long-term rainfall variability across the Southern Africa region and the continent at large are complex as they are influenced by the changes taking place in the oceans together with atmospheric circulation. These include the ENSO (El Nino southern oscillation cycle) [23] and the fluctuations in sea surface temperature (SST) which has no relationship with the ENSO [95] and human anthropogenic activities on the environment leading to land cover changes [88]. The results of this study and previous studies on rainfall trends show a clear indication of how variable rainfall has been over the centuries.
The study showed no widespread evidence that rainfall in NRS had decreased or increased due to climate variations. The study based its evidence on results from the eight stations studied with annual and seasonal rainfall and showed no strong evidence that seasonal or annual rainfall has declined or increased since only three of the eight stations had statistically significant results. The study, therefore, rejects the null hypothesis as most of the stations (62.5%) showed a statistically insignificant change over time. Using two methods revealed that rainfall trends from three stations are statistically significant. MK had depicted the Elandspoort station as having a statistically significant increasing rainfall trend. Additional stations, Nylsvley and Groenfontein were depicted by the QR method as having statistically significant decreasing trends in annual rainfall.
The Nylsvley rainfall station, located within the flood plain with over 53 years of recorded annual rainfall data, has received a maximum and minimum annual rainfall of 733 mm in 1988 and 231 mm in 2001, respectively, shown in Figure 2. The Nylsvley rainfall station is crucial as it provides rainfall statistics of the flood plain and the Nature Reserve. The trend of low annual rainfall in Nylsvley has been recorded from the years 1995 to 2003 which received between 231 mm and 483 mm. This decrease could impact on vegetation growth and the availability of flooding water within the flood plain and its biological diversity. This decrease agrees with the results of the study that the flood plain has consistently received low rainfall as shown by the decreasing trends in both rainfall and NDVI, as shown in Figure 6. The decrease in rainfall in the Nylsvley agrees with the ideology that the flood plain is experiencing dry years as exacerbated by low rainfall being received in some of the years.
Though many stations showed decreasing rainfall trends, they were mostly not statistically significant or proved and this could imply that there are interannual and seasonal variations of rainfall in NRS. The insignificant trends can also be linked to a relatively high year-to-year variability of rainfall in South Africa (including Limpopo Province) which is associated with some cloud band and ENSO [96]. Mazvimavi [64] highlighted no changes over-time for rainfall in Zimbabwe and stated that global warming may cause changes in rainfall patterns, but these effects were not statistically significant. The author later stipulated that the Southern Africa region is semi-arid to arid with inter-annual rainfall variability. In the region, including Limpopo Province of South Africa, there is need to be vigilant to such decadal or inter-annual changes to be ready to respond to uncertainty raised by climate change and/or anthropogenic factors. New et al. [27] in their study highlighted the decrease in precipitation of Southern Africa even though the study was not statistically significant either.
Onema and Taigbenu [18] also emphasised in their study that the NDVI-rainfall relationship is linked, hence they expected to see a trend in most of their data. This trend has also been shown in the Nylsvley wetland indicating that, when restoration activities take place in the wetland, management should be informed of the water challenges needed for vegetation growth and other life forms including flooding in the wetland.
The major decrease in rainfall in the years between 1991 and 1993 agrees with studies by Nicholson et al. and Richard et al. [97,98], as there was a serious drought in the Southern Africa region that caused significant summer rainfall variations. This could have led to decrease in NDVI in Nylsvley as the area could have received less rainfall to promote plant vigour. The annual and seasonal rainfall trends for the Nylsvley station from 1953 to 2003 have shown a decreasing trend, as seen in Table 1. A study on vegetation greenness (NDVI) and rainfall variability conducted in the Sahel region in the period 1982-2003 concluded that rainfall remains the main causative factor for vegetation greenness increase [99]. Other studies by Martiny et al. and Barbosa et al. [8,66], also support the NDVI-rainfall relationship, although anthropogenic factors impact negatively too [26]. The decrease in NDVI within the study area could also be attributed to other factors such as anthropogenic activities, for example, agriculture in the area [45]. Agriculture leaves most surfaces as bare ground with infertile soils. As such, vegetation in these old fields takes a long time to naturally recover [100], which could be the case in some parts of Nylsvley.
Nylsvley has a history of grazing both from domesticated and wild animals [45]. The growing of crops ceased in part of the wetland in 1974 and the remaining part in 1997, as it became a protected area. The degradation of the wetland through the removal of vegetation has forced the government, through its project implementer (Working for Wetlands), to embark on restoring the wetland so that it can return closer to its natural state. The decrease in NDVI due to less rainfall has implications in ecosystem development and services. For instance, in the Nylsvley wetland, game animals obtain their food through grazing and less vegetation has implications when planning for these grazers. Vegetation also provides a habitat for migratory birds that visit Nylsvley, and this attracts tourists to visit for bird watching and is a source of income to this region and employment creation in the tourism sector. Thus, the rainfall-NDVI relationship is critical as it helps managers to prepare for early negative climatic effects.

Conclusions and Recommendations
This study addressed the historical rainfall trends in NRS and rainfall-vegetation relationships to assess climate effects on wetland vegetation growth. NDVI was thus used as a proxy of vegetation and showed that as rainfall decreased or increased over the years under study, NDVI seemingly followed suit (Figures 6 and 8). The study utilised two methods in assessing historical rainfall trends which revealed that three stations (Groenfontein, Elandspoort, and Nylsvley) have fluctuating rainfall patterns with Nylsvley consistently reporting decreasing rainfall trends. Three out of the eight rainfall stations studied demonstrated evidence of the changes in annual and seasonal rainfall patterns. The decrease in rainfall in Nylsvley has negatively impacted the floodplain wetland with implications on tourists intending to visit it for bird watching and other activities.
The trend of mean NDVI and monthly rainfall shows a positive relationship in some years under study with other years showing negative NDVI increase. This may suggest a combination of years of drought and the effects of negative human activities such as grazing, which is still active in Nylsvley.
The changes seen in this study area require proper management and monitoring for sustainable wetland protection of this Ramsar site. This study also suggests that the knowledge of changes in rainfall and NDVI patterns could be used in planning for the effective utilisation of water resources throughout the Nyl River system. Since the study showed varying stations indicating decrease and/or increase in rainfall trends, there is a need to investigate the effects of anthropogenic changes such as land cover/use changes which affect the availability of water resources and the success of natural vegetation restoration within Nylsvley Acknowledgments: The authors thank South Africa Weather Services and Agricultural Research Council for providing long term rainfall data. GreenMatter is also acknowledged for the financial support in the form of a bursary to the student that made it possible to accomplish this research task.

Conflicts of Interest:
The authors declare that there are no conflicts of interests.