Next Article in Journal
Surface Temperature Trend Estimation over 12 Sites in Guinea Using 57 Years of Ground-Based Data
Next Article in Special Issue
Impact of Climate and Land Use/Land Cover Change on the Water Resources of a Tropical Inland Valley Catchment in Uganda, East Africa
Previous Article in Journal
Stakeholders’ Perception of Climate Actions in Some Developing Economies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Investigation into the Future Changes in Onset and Cessation of Rain and Their Variability over the Aswa Catchment, Uganda

1
Faculty of Engineering and Technology, University of Botswana, Gaborone P/Bag 0061, Botswana
2
Department of Agricultural and Bio Systems Engineering, Makerere University, Kampala P.O. box 7062, Uganda
3
Department of Environmental Science, Faculty of Science, University of Botswana, Gaborone P/Bag 0022, Botswana
*
Author to whom correspondence should be addressed.
Climate 2020, 8(6), 67; https://doi.org/10.3390/cli8060067
Submission received: 28 April 2020 / Revised: 15 May 2020 / Accepted: 19 May 2020 / Published: 29 May 2020
(This article belongs to the Special Issue Climate Change and Water-Related Agricultural Risks)

Abstract

:
Future global warming may result in extreme precipitation events leading to crop, environment and infrastructure damage. Rainfall is a major input for the livelihood of peasant farmers in the Aswa catchment where the future rainfall variability, onset and cessation are also likely to be affected. The Aswa catchment has limited rainfall data; therefore, use of secondary datasets from Tropical Rainfall Measuring Mission (TRMM) is considered in this study, based on the close correlation of the recorded and TRMM rainfall. The latter was used to calibrate the statistical downscaling model for downscaling of two general circulation models to simulate future changes in rainfall. These data were analyzed for trends, wet and dry conditions/variability; onset and cessations of rain using the Mann–Kendall test, Standardized Precipitation Index (SPI) and the cumulative percentage mean rainfall method, respectively. Results show future rainfall is likely to increase, accompanied by increasing variability reaching as high as 118.5%. The frequency of SPI values above 2 (extreme wetness) is to increase above current level during mid and end of the century. The highest rainfall variability is expected especially during the onset and cessation months, which are generally expected to come earlier and later, by up to four and five weeks, respectively. The reliability worsens from the midterm (2036–2065) to long term (2066–2099). These likely changes in rainfall quantities, variability, onset and cessation months are some of the key rainfall dynamics that have implications for future arable agriculture, environment and water resource availability and planning over the Aswa catchment, as is increasingly the case elsewhere.

1. Introduction

Projected future global warming may result in extreme precipitation events causing drought and flooding, which in turn may result in crop yield reduction, increased soil erosion and damage to infrastructure [1,2,3,4,5]. The recent report by the Intergovernmental Panel on Climate Change (IPCC) [6] indicates that global warming above pre-industrial levels (1900) of 1.5 °C and beyond would have adverse impacts on communities that depend on agriculture in the least developed countries. Rainfall is key to livelihoods of peasant farmers who depend on rain fed agriculture [7,8,9,10,11]. Declines in rainfall, increased variability and/or uncertainty would increase the risk of hunger and poverty especially for the peasantry. Fishman [12] argues that an increase in rainfall variability amidst increased annual rainfall overrules the benefit expected from the increase on crop yields. Increased future rainfall variability is, therefore, a great concern for the socio-economic wellbeing of the communities living in Aswa catchment and other areas with similar climatic conditions.
Over the country of Uganda as a whole, more than 72 % of the economically active population derive their livelihood from agriculture, which contributes about 24.6 % of the Gross Domestic Product (GDP) [13]. The study area of Aswa catchment offers a typical Ugandan and African setting where the future socio-economic wellbeing of the community who mostly depend on agriculture is uncertain due to impacts of climate change. More than two million people, of which more than 70% are peasants, live in the Aswa catchment [14,15]. The catchment is characterized by periodic and extended drought spells, leading to chronic food insecurity and famine. With an annual population growth rate of over 3%, and within the context of climate change, the future is increasingly becoming more uncertain. Detailed research on variability in the future rainfall changes and their onset and cessation at the catchment and local level is, therefore, required for effective crop, environment and water resource planning and management.
The Aswa catchment also lacks adequate recorded climate data, as is generally the case across most of the African continent and other developing countries in the world. This challenge, among others, could be addressed by using secondary data, such as the Tropical Rainfall Measuring Mission (TRMM) [15,16]. Such data can be used to downscale General Circulation Models (GCMs) to appropriate local catchment resolutions. GCMs are used to project changes in climatic parameters, such as rainfall under various climate change Representative Concentration Pathway (RCP) scenarios, defined by the Intergovernmental Panel for Climate Change (IPCC). However, these projections are defined at coarse grid resolutions at which the GCMs are configured and, thus, being unable to resolve important sub-grid scale variations. Therefore, the GCMs cannot be used directly for climate change impact studies of small catchments such as Aswa; thus, downscaling is required to generate corresponding data. In addition to GCMs not providing catchment scale information, despite being our major source of knowledge about the future, they can also provide biased information [17,18].
For downscaling GCMs, two downscaling methods, dynamical and statistical, are available. Other nonlinear methods using artificial neural network (ANN) are also being used lately [19,20] for generating long-term rainfall patterns. Statistical Downscaling (SD) is usually preferred over dynamical downscaling as the former produces local weather and climatic time series, which is cheap, readily transferable and computationally less demanding by use of appropriate statistical or empirical relationships with predictor variables and, therefore, has been widely used in climate change studies [21]. Products of statistical downscaling can be limited when there are some systematic errors in input data used to calibrate the downscaling model, however. Systematic biases for the current climate are also unavoidable in GCMs [22]. The biases can be propagated through the downscaling process with far reaching implications on the simulations and any subsequent applications [23,24]. The first step towards reducing the bias and improving the simulations is usually to avoid using unrealistic or inaccurate datasets for calibrating the downscaling. In this regard input data used for downscaling can either be bias corrected before downscaling or the simulations are bias corrected themselves after downscaling [17,24]. Caution is also needed during post processing of the simulations to keep minimal chances that the original climate change signal could be affected as this alteration can substantially affect impact model results [25]. Depending on the level of biasness of downscaled products from observations and the intended use of the products, on a case-by-case basis, bias correction may not be critical.
In Statistical Down Scaling Model (SDSM), a combination of Multiple Linear Regression (MLR) and the Stochastic Weather Generator (SWG) is used for generating climate scenarios for assessment of various parameters of the data [21,26]. This study therefore, uses secondary rainfall data (TRMM) to calibrate SDSM in a downscaling experiment of two General Circulation Models (GCMs)—Hadley Centre Coupled Model version 3 (HadCM3) and Canadian Earth System Model (CanESM2))—for simulation of future rainfalls for the data-constrained, but socio-economically important Aswa catchment located in Uganda, within the context of changing climate, as reported in previous studies [15]. The projected future rainfall is analyzed, using various methods, for changes in variability, onset and cessation. These changes are expected to have a profound impact on socio-economic wellbeing of the farmers and the environment particularly in Aswa catchment, as well as areas and regions that have similar climatic and socio-economic environment.

2. Materials and Methods

2.1. Study Site

The Aswa River catchment, at its confluence with River Nile, lies mainly in the northeastern part of Uganda with a small portion extending into South Sudan. The total catchment area spreads over 31,000 km2 and lies between 31.9° E–34.25° E latitude and 2.05° N–4.09° N longitude. Figure 1 shows the location of the catchment area and the points from which TRMM data were downloaded. Coordinates of the TRMM data locations are shown in Table 1. The catchment has a wide range of annual rainfall amounts and distribution pattern. The northeastern parts have a semi-arid climate with annual rainfall varying between 500 and 800 mm, while the western and southern parts have semi-humid and humid climate with annual rainfall ranging from 800 to more than 1400 mm in some places.

2.2. Rainfall Data

Secondary data from TRMM rainfall were used for this study in place of the scares and inconsistent recorded rainfall data over the catchment. Gridded satellite based TRMM data are available from 2000 to 2015 (16 years). The TRMM rainfall data grid locations used in this study are shown in Table 1. The TRMM grid resolution is 32 by 32 Km, while for the GCM used is within 290 by 435 km. The TRMM data can therefore be used for downscaling GCM data to a finer resolution.

2.3. Data Processing

Limited observed daily rainfall data are available within the vicinity of the Aswa catchment from weather stations such as Lira, Gulu and Kitgum. The relation between the station and secondary data (TRMM) was investigated through correlation to find out if they were closely related enough to give confidence of using TRMM data in calibrating the statistical downscaling model over the catchment. Sometimes secondary data may have interventions due to malfunction of satellite equipment or change of equipment [27] and, therefore, should be tested before use so that the homogenous section is used for the study. Before this validation, TRMM data were put to intervention test analysis to find out if the data set was homogeneous using the standard normal homogeneity test (SNHT) [15].

2.4. Downscaling and Projection of Future Rainfall

TRMM rainfall data were used to calibrate the SDSM for subsequent downscaling of two GCMs (HadCM3 and CanESM2). Many GCMs are in use today, globally, with varying degrees of accuracy. CanESM2 and HadCM3 have been successfully used in the region for climate change studies [28,29] and have, therefore, been adopted for use in this study. The high emission scenario RPC 8.5 (CanESM2), A2 (HadCM3) predictor data variables were used for this study to represent the worst-case scenario. The predictor data for the two GCMs were obtained from two different websites: HadCM3 (http://www.cccsn.ec.gc.ca/?page=pred-hadcm3) and CanESM2 (http://www.cccsn.ec.gc.ca/?page=pred-canesm2).
The predictor variables were screened in order to choose appropriate predictors. This stage is critical because the choice of predictors largely determines the accuracy of the downscaled climate scenario data. The procedure for predictor selection has been well laid by Hasan, et al. [30]. A correlation matrix between the predictor and predictand (TRMM rainfall) was made. The predictor with the highest correlation coefficient was selected and called the super predictor (SP). After choosing the super predictor, the remaining predictors were regressed individually, together with SP, and the absolute correlation coefficient between the predictor and predictand and between individual predictors were obtained together, with their partial correlations and P values. Any predictor with p-value > α (0.05) was removed. The percentage by which partial correlation is lowered in respect to the correlation coefficient PRP (Equation (1)) was calculated for each predictor. The other predictors where selected by calculating the PRP with respect to the absolute correlation using the following equation.
P R P = ( P . r R 1 R 1 )
where:
  • PRP is the percentage by which partial correlation is lowered in respect to the correlation coefficient,
  • P.r is defined as partial correlation coefficient and
  • R1 is the value of correlation coefficient between predictor and predictand.
The predictor that had the minimum PRP was selected as the second most suitable predictor. This same process was done for the third and other following predictors. As a result, the super predictor and the other remaining predictors will have almost zero, or insignificant multi-co-linearity. Five predictors were selected, as can be seen, for example, for the Kitgum location in Table 2, with mean temperature at 2 m as the super predictor. The correlations are statistically significant (p < 0.05).

Calibration and Validation of SDSM Model and Generation of Future Rainfall

SDSM version 4.2.9 was calibrated and validated using the selected appropriate predictors and TRMM rainfall data as predictand for various locations within the Aswa catchment. Such results of predictor selection are shown for Kitgum station in Table 2.
Using the developed model, precipitation data were simulated for 1961–2001(HadCM3) and 1961–2005 (CanESM2). The SDSM was used to generate 20 ensembles. The mean of the ensembles and the TRMM and SDSM weather generator monthly rainfall data were compared. The coefficient of determination, R2, between the (TRMM) and generated monthly data series, Bias (ratio of sum of totals of simulated and TRMM rainfall), F values and the standard error (SE), were used for evaluation of the statistical calibration and validation. The calibration and validation results for various locations within the catchment can be seen in Table 3. The results show that the statistical calibration and validation are good enough to be relied on. The correlation coefficients range from 0.83 to 0.97. The analysis of variance (ANOVA) F values are much greater than the critical values and the standard error for rainfall range from 9.62 to 29.25 mm/month. These values compare well with other studies using SDSM and validation studies in the region [28,29,31]. The bias values are low, with an average value of 1.018.
The calibrated SDSM was used to downscale HadCM3 and CanESM2 for a generation of future rainfall data for the near future (2006–2035), medium term future (2036–2065) and long term future (2066–2099). These data were analyzed for trends and variability following the procedure outlined in Section 2.5.

2.5. Rainfall Trend Analysis

The historical (secondary) TRMM (2000–2015) data were subjected to a homogeneity test in order to identify any interventions in the data. The homogenous portion of the data was then used for further analysis. The simulated near future (2006–2035), medium term (2036–2065) and the long term (2066–2099) future daily and monthly rainfall data were analyzed for trend and variability. The Mann–Kendall test was used for trend analysis and is given by Equations (2) and (3).
τ = S 0.5 n ( n 1 )
S = i = 1 n 1 j = i + 1 n     s i g n    ( x j x i )
where;
  • xj = the observation at time j
  • xi = the observation at time i
  • n = total number of observations
  • Sign (xj − xi) = 1 if computed value > 0
  •       = 0 if computed value is 0
  •       = −1 if computed value < 0
Further analysis was carried out using descriptive statistics such as mean, standard error, standard deviation and coefficient of variability (CV) for the historical (TRMM) and simulated future rainfall data.

2.6. Onset and Cessation of Rainfall Dates

The cumulative percentage mean rainfall method was used to determine the onset and cessation of the historical, near future (2006–35), medium term (2036–65) and long term (2066–99) future simulated rainfall. This method has been used successfully by many researchers [7,32,33,34] for determination of the rainfall onset and cessation dates. The procedure for the method is as follows:
  • the percentage of mean annual rainfall that happens at every 5-day interval is derived;
  • the cumulative percentage of the computed percentage at 5-day intervals is derived;
  • the cumulative percentage at 5-day intervals is plotted for the whole year; and
  • identification of the time of rainfall onset and cessation during the year or season.
The point of first maximum positive curvature and last maximum negative curvature on the graph of the cumulative percentage indicate the mean periods of rainfall onset and cessation, respectively. Another alternative criterion for the onset of the rain is the timing of an accumulated 7% to 8% of the annual rainfall, and the cessation begins after the accumulation of 90% of the annual rainfall.
The cumulative percentage mean rainfall method for determining onset and cessation of rainfall can use either the rainfall amounts or rainy days [7,32,33,35]. Both methods have been applied in this study. There are other methods for determination and definitions of onset and cessation of rain [36,37]. The criterion used for defining onset and cessation of rain can be based on soil water balance/water stress coefficient, the first wet day with accumulated X mm of rainfall during a period of a certain number of days, rainfall–evapotranspiration relation model, and others. The threshold for a rainy day also varies: 0.85 mm, 1 mm and other values [36,38]. The method of cumulative percentage mean rainfall is commonly used by many researchers. It is preferred to other methods for determination of rainfall onset and cessation dates due to its efficiency and freedom from assumption of rainfall threshold values and its reliance on rainfall data alone [39]. The minimum rainfall value for a rainy day was set as 1 mm. This value was adopted from Ngetich, Mucheru-Muna, Mugwe, Shisanya, Diels and Mugendi [36]. The days with rainfall below 1 mm were considered as non-rainy or dry days.

2.7. Standardized Precipitation Index (SPI)

There are several drought indices in use but the most commonly used are the Palmer Drought Severity Index (PDSI) and the Standardized Precipitation Index (SPI) [40,41]. The PDSI is related to a water balance. The input parameters include moisture supply, evaporation, precipitation and runoff. However, the SPI is an index with only precipitation as an input variable [42]. The standardized Precipitation Index (SPI) was used to determine occurrence of extreme rainfall events for the current period and future projections of rainfall in this study. The advantage of SPI over other indices is its robustness to detect different types of droughts such as hydrological, meteorological and agricultural at different time scales (1, 3, 6, 12, 24 up to 60 months) using only rainfall data set.
The SPI computation method was developed by McKee, et al. [43] to evaluate the relative deviations of precipitation from the average normal values and the Gamma distribution is generally used in SPI calculation. It has been used for many studies all over the world [41,44,45,46,47,48] to evaluate occurrences or drought and wet conditions. SPI is based on normalized data and, therefore, does not vary spatially, and as such, droughts and wet occurrences can be assessed in different regions. The SPI is calculated by fitting a certain probability density function to the frequency distribution of rainfall series over the period for the analysis. The resulting fitted function is used to generate cumulative distribution of the data and consequently transformed into the standardized normal distribution used to define SPI.
The Gamma distribution has two main parameters; the probability density function (PDF) and cumulative distribution function (CDF).
The PDF can be expressed as:
g ( X ) = 1 β α . Γ α x α 1 e x / β
where,
  • α = shape parameter (α > 0),
  • β = scale parameter (β > 0),
  • Γα = a gamma function of α and can expressed as;
Γ ( α ) = 0 y α 1 e y d y
The PDF of an observed amount of precipitation for a given month and period of time is expressed as:
G ( x ) = 1 β ^ α ^ . Γ ( α ^ ) 0 x x α 1 e x / β ^ d x
The PDF is transformed to the standard normal variate Z, which is the value of SPI.
S P I = x i x ¯ σ
where, xi is the precipitation of the selected period during the year i, x ¯ is the long term mean precipitation and σ is standard deviation for the selected period. Positive and negative SPI values correspond to wet and drought periods respectively. Table 4 gives the classification of SPI values.
In this study, hydrological drought was evaluated for Aswa catchment. A soft wear SPI generator developed by the National Drought Mitigation Center of University of Nebraska-Lincoln was used to calculate the monthly SPI values for the period 2006–2099.

3. Results

3.1. Correlation of Station and TRMM Rainfall

The TRMM and Observed monthly rainfall data for Lira station were correlated as seen in (Figure 2). The correlation gives coefficient of determination, R2 value of 0.65 and standard error of 47.86 mm. This indicates that the observed and TRMM rainfall data are closely correlated and the TRMM data can be used to represent observed data in case of unavailability of the latter.

3.2. Simulated Future Rainfall

Based on computed test statistics against threshold, homogeneity test results (Table 5) show that the TRMM and weather generated rainfall data series 1975–2005 are homogeneous and, therefore, the whole data series was used to calibrate the SDSM for the study. Future rainfall scenarios were subsequently simulated using the calibrated and validated SDSM for two GCMs (HadCM3 and CanESM2) for various locations within the vicinity of the Aswa catchment. All of the results show future rise in rainfall amounts. Figure 3 shows the graphical simulation results for five locations based on CanESM2 and HadCM3 GCMs.
The annual rainfall at current (2015) level for Agoro is about 800 mm and is expected to rise to about 1000–1200 mm by 2035, 1200–1300 mm by 2065 and up to 1500 mm by 2099 (Figure 3). The annual rainfall at current level for Amuria is about 1500 mm and is expected to rise to about 1700–2000 mm by 2035, 2200–2300 mm by 2065 and up to 3500 mm by 2099 (Figure 3).
The Mann–Kendall test statistics also show that there is going to be positive trend for the simulated future annual rainfall over Aswa catchment as can be seen from Table 5. The Kendall tau values range from 0.172 to 0.683. The increase in positive trend is stepper for the medium term (2036–2065) and long term (2066–2099) periods. This increase in rainfall is also associated with increased variability as we show in the next sections.
The increase in annual rainfall is expected to have impact on agricultural production, water resources and environment management. Some of the impacts may be positive (increased crop yields, water availability) while others may be negative (erosion, flooding, infrastructure damage, increased waterborne diseases).

3.3. Future Rainfall Onset and Cessation Dates

The onset and cessation of rainfall is to vary in the future and depart from the present pattern. In the mid and long term future, onset of rain is expected earlier by 3–4 weeks and rainfall cessation is expected to also commence later than it occurs now, by 1–5 weeks. The rainfall cessation dates are more variable than the onset dates. Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show examples of the mean rainfall onset and cessation dates for present, near, medium and long term future simulations for Agago and Kaabong locations using daily rainfall and rainy days data. These two locations can be seen from Figure 1. Kaabong represents the semi-arid part of the catchment and Agago the humid part. Results for both CanESM2 and HadCM3 predictions are comparably close to each other. However, the future rainfall onset and cessation dates are to become more variable. Some studies from Ghana [32] found that the use of daily rainfall or rainy days for determination of onset and cessation of rainfall give similar results. Odekunle [33] showed that the use of rainy days for determination of onset and cessation of rainfall in Nigeria gave more accurate results than the use of rainfall amounts. However, in our work for Aswa catchment, there is no evidence that this is the case as the two methods gave similar results.
In view of the observations above, it is expected that the future shift of the onset and cessation of rainfall will have impact on the planning and management of farming activities in Aswa catchment. It will also impact crops yields in the catchment. Introduction of new crop varieties, which have short growing periods, irrigation, soil and water conservation may be some of the mitigation measures to be adopted by farmers in the catchment.

3.4. Rainfall Variability Analysis

The simulated future rainfall data series are highly variable. The monthly coefficients of Variation (CV) reach up to 118.5%. Coefficient of variation is a measure of how reliable the rainfall pattern is. A higher variability implies a highly unreliable rainfall pattern. Generally, rainfall with a CV value above 30% is considered to be highly variable [49]. Table 6, Table 7 and Table 8 summarize the results for Agoro, Kaabong and Gulu locations, respectively. Kaabong represents the semi-arid region, Agoro represents the semi-humid region and Gulu represents the humid region of the catchment, as can be seen in Figure 1.
As suggested by the high CV values, the rainfall is expected to be more variable during the onset and cessation months of February–April and September–November, respectively. There is a short dry spell that occurs around the months of June-July. The rainfall projection based on CanESM2 show that the short dry spell of July is to be more severe in the mid and long term future than the current period as shown in Figure 12 and Figure 13 for Agoro and Gulu locations. Extreme heavy rainfall events are expected to increase during February-May, September and November months.
Based on CanESM2 rainfall projection, the three months’ time scale SPI values show that the long term future rainfall pattern is going to be dominated by wet events (severely wet to extremely wet) as compared to the current and near future situation which is dominated by droughts (moderately dry to extremely dry). Figure 14, Figure 15 and Figure 16 show the SPI_3 values of 2006–2099 for Kaabong, Agoro and Gulu locations respectively. Table 9, Table 10 and Table 11 show percentage of variation of SPI values for the period 2006-2099 for Kaabong, Agoro and Gulu locations respectively. The moderately, severely and extremely wet SPI values increase from mid to end of the century compared to the current situation. In Agoro, for example, the percentage of extremely wet SPI values (>2) rise from 0 (2006–2035) to 0.8 (2036–2065) and then to 6.4 (2066–2099). This is in agreement with the high CV values of Kaabong and Agoro, as it is observed in Table 6 and Table 7. These results suggest that the semi-arid and semi-humid section of Aswa catchment will have high rainfall variability in the future. Nicholson, et al. [50] also found that semi-arid climate in Africa are prone to extreme rainfall variability. Batisani and Yarnal [51] and Byakatonda et al. [52] reported increasing rainfall variability in semi-arid Botswana. The expected increased future rainfall variability will have profound impact on various sectors such as agricultural production, water resources, infrastructure and the general environment in Aswa catchment. Droughts will cause crop failure and extreme wet conditions will lead to water logging, flooding and soil erosion, which also affect crop yields negatively. Studies from different parts of the world suggest that high rainfall events can lead to crop yield reduction. Studies from Ghana revealed negative impact on crop yields due to heavy rainfall events [53,54].

4. Discussion

Due to scarce and inconsistent station rainfall data in Aswa catchment, secondary data from TRMM was used to downscale two GCMs (CanEMS2 and HadEM3) for investigating likely future rainfall quantities, variability onset and cessation of rainfall months for this study. However, secondary data from satellite sometimes present challenges of non-homogeneity [27] and week correlation with the actual recorded data. The intervention tests show that the TRMM data is homogeneous. The correction between TRMM and the limited recorded rainfall data in Aswa catchment show a coefficient of determination, R2 value of 0.65 (R = 0.806), which gives confidence of using the TRMM rainfall data as proxy for station data. The biases are generally low with an average of 1.018. TRMM rainfall validation studies carried out in Ghana [31] revealed similar bias values ranging from 1.047 to 1.479. The good correlation between station and TRMM rainfall and minimal bias between downscaled and TRMM rainfall are deemed sufficient to allow for use of TRMM rainfall in the absence of station rainfall data for simulation of future scenarios. The downscaled simulated future rainfall using SDSM, CanEMS2 and HadEM3 CGM predictors show increased rainfall in the future. In reference to the SDSM calibration results, we expect CanESM2 (with higher R2 values) based simulated future rainfall data for Aswa catchment to be more accurate than HadCM3 based data. The correlation coefficients for CanESM2 are marginally higher and the standard errors are lower than for HadCM3 (Table 3).
The results of Man Kendal tau values are positive, as can be seen for example, for Agoro and Amuria locations ranging from 0.172 to 0.683 (Table 5). This suggests a strong positive trend in the future. Agoro represents the semi humid section and Amuria represents the humid section of the catchment. This trend is expected to become steeper as we move from mid (2036–2065) to end of the century (2066–2099). These results are also supported by the SPI values, which can be seen from Figure 14, Figure 15, Figure 16, Table 9, Table 10 and Table 11. The percentage of SPI values greater than 1.5 (severely wet and extremely wet) are to increase as we move to mid and end of the century. This suggests that the future rainfall is not only expected to increase but also have extremely wet patterns, which are likely to affect crop production, infrastructure and the environment negatively.
The expected increase in rainfall trend could be good for improving crop yields in an area such as the Aswa catchment, which has been suffering from droughts and water stress. However, analysis of the rainfall patterns show increased and high monthly CV values in the future. CV values above 30% suggest high variability and poor rainfall distribution. The CV values, as can be seen from Table 6, Table 7 and Table 8 for Agoro, Kaabong and Gulu locations respectively, are high reaching up to 118.5%. The onset and cessation months of February–April and October–November respectively show the highest rainfall variability. The rainfall patterns during these months are expected to become more unreliable as we move to mid centaury to end of the century. The CV values for future simulated rainfall data based on CanESM2 are generally higher than for the data sets based on HadCM3. However, looking at the higher correlation coefficients for CanESM2 data sets during calibration and validation, we consider the corresponding CV values for the data sets to be more reliable.
The analysis of future rainfall onset and cessation dates results also support the indications of the CV. Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show the variation of future rainfall onset and cessation dates for Agago (humid) and Kaabong (semi-arid) locations for CanESM2and HadCM3 based downscaling using both rainfall and rainy days data. The future rainfall onset and cessation dates are expected to come earlier and later than the current dates by up to four and five weeks respectively. For Agago the future rainfall onset is expected to come in February instead of March and cessation is expected in November instead of late October as is currently being experienced. For Kaabong, future rainfall onset is expected to come late February/early March instead of late March as is currently happening, while the future cessation is expected in November instead of the current late October. The two GCM gave similar results for each location. The change in rainfall onset and cessation dates is likely to affect the farming activities of the peasant farmers and crop yields.

5. Conclusions

In this study, we examined the future changes in rainfall patterns, trends, extreme occurrences, onset and cessation of rainfall over Aswa catchment using SPI, Mann–Kendall statistics, descriptive statistics (coefficient of variability), and cumulative percentage mean rainfall methods to provide insight into near future (2006–2035), midterm (2036–2065) and long term (2066–2100) future periods.
The results suggest that annual rainfall over Aswa catchment is going to increase especially in the mid and long term future by up to 100%, and even more. This could be a good development for the farmers if the rainfall turns out to be reliable and well distributed. However, the future presents a highly variable and unreliable rainfall pattern. This can be seen from the high CV (up by 118.5%) and the rise in percentage of SPI values for extremely wet periods. The study reveals SPI values for extremely wet periods increasing from zero in (2006–2035), to 0.8 in (2036–2065) and to 6.4 in (2066–2099) periods for Agoro which lies in the semi-humid region of the catchment. For Kaabong located in the semi-arid region of the catchment, the rise in percentage of SPI values for extremely wet periods is 0.3, 0.3 and 5.4 respectively. For Gulu located in the humid region of the catchment, is from 0 to 0.1 to 1.9 respectively. This suggests that the semi-arid and semi-humid regions of the catchment will experienced more extremely wet conditions in the future compared to the humid region. Extremely dry periods are, however, expected to reduce although they will continue to be experienced as can be seen from Table 9, Table 10 and Table 11.
The onset and cessation of rain are also expected to shift backwards and forward by up to four and five weeks respectively. This in turn is expected to affect farming activities such as field preparation, planting and harvesting. The future scenarios are increased, highly variable and unreliable rainfall in Aswa catchment, with the drier semiarid section being the most affected. This may result in floods, water logging, severe soil erosion and even droughts with consequent reduction in crop yields and negative impact on the livelihood of the farmers. These changes can affect the future agricultural production, water availability, environment and socio-economic wellbeing over the Aswa catchment. The results of the study can be of use for future catchment water resources planning and management, agricultural and socio-economic planning, and introduction of appropriate mitigation measures such as provision of supplementary irrigation, introduction of new crop varieties/farming system, soil and water conservation measures, drainage and reservoir development.

Author Contributions

Conceptualization, M.I., P.T.O., B.P.P. and D.B.M.; Formal analysis, M.I.; Investigation, M.I.; Methodology, M.I.; Resources, P.T.O., B.P.P. and D.B.M.; Supervision P.T.O., B.P.P. and D.B.M.; Validation M.I., P.T.O., B.P.P. and D.B.M.; Visualization, M.I.; Writing—original draft, M.I., P.T.O., B.P.P. and D.B.M.; Writing—review & editing M.I., P.T.O., B.P.P. and D.B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This study was made possible by the sponsorship of EU’s Intra Africa Mobility to Enhance Training of Engineering Graduates in Africa (METEGA) Program at University of Botswana. The authors wish to acknowledge NASA for providing free TRMM satellite data listed in Table 1 of this paper. We would like to thank the developers of SDSM for providing this model for free. The authors also would like to acknowledge the limited available rainfall data for Aswa catchment provided by the Directorate of Water Resource Development (DWRD), Uganda. We would like to thank the University of Botswana for provision of funds through the Office of Research and Development (ORD) for data acquisition.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mukherjee, S.; Aadhar, S.; Stone, D.; Mishra, V. Increase in extreme precipitation events under anthropogenic warming in India. Weather. Clim. Extremes 2018, 20, 45–53. [Google Scholar] [CrossRef]
  2. Lehmanni, J.; Coumou, D.; Frieler, K. Increased record-breaking precipitation events under global warming. Clim. Chang. 2015, 132, 501–515. [Google Scholar] [CrossRef]
  3. O’Gorman, P.A. Precipitation Extremes Under Climate Change. Curr. Clim. Chang. Rep. 2015, 1, 49–59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Kharin, V.V.; Zwiers, F.W.; Zhang, X.; Wehner, M. Changes in temperature and precipitation extremes in the CMIP5 ensemble. Clim. Chang. 2013, 119, 345–357. [Google Scholar] [CrossRef]
  5. Schewe, J.; Heinke, J.; Gerten, D.; Haddeland, I.; Arnell, N.W.; Clark, D.B.; Dankers, R.; Eisner, S.; Fekete, B.M.; Colon-Gonzalez, F.J.; et al. Multimodel assessment of water scarcity under climate change. Proc. Natl. Acad. Sci. USA 2014, 111, 3245–3250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. IPCC. Summary for Policymakers. In Global Warming of 1.5 °C; An IPCC Special Report on the Impacts of Global Warming of 1.5 °C above Pre-Industrial Levels and Related Global Greenhouse Gas Emission Pathways, in the Context of Strengthening the Global Response to the Threat of Climate Change; Masson-Delmotte, V., Zhai, P., Pörtner, H.O., Roberts, D., Skea, J., Shukla, P.R., Pirani, A., Moufouma-Okia, W., Péan, C., Pidcock, R., et al., Eds.; World Meteorological Organization: Geneva, Switzerland, 2018; p. 32. [Google Scholar]
  7. Recha, C.W.; Makokha, G.L.; Traore, P.C.S.; Shisanya, C.; Lodoun, T.; Sako, A. Determination of seasonal rainfall variability, onset and cessation in semi-arid Tharaka district, Kenya. Theor. Appl. Clim. 2011, 108, 479–494. [Google Scholar] [CrossRef] [Green Version]
  8. Traore, B.; Corbeels, M.; Van Wijk, M.; Rufino, M.C.; Giller, K.E. Effects of climate variability and climate change on crop production in southern Mali. Eur. J. Agron. 2013, 49, 115–125. [Google Scholar] [CrossRef]
  9. Sultan, B.; Roudier, P.; Quirion, P.; Alhassane, A.; Müller, B.; Dingkuhn, M.; Ciais, P.; Guimberteau, M.; Traore, S.; Baron, C. Assessing climate change impacts on sorghum and millet yields in the Sudanian and Sahelian savannas of West Africa. Environ. Res. Lett. 2013, 8, 014040. [Google Scholar] [CrossRef]
  10. Omoyo, N.N.; Wakhungu, J.; Oteng’I, S. Effects of climate variability on maize yield in the arid and semi arid lands of lower eastern Kenya. Agric. Food Secur. 2015, 4, 881. [Google Scholar] [CrossRef] [Green Version]
  11. Adhikari, U.; Nejadhashemi, A.; Woznicki, S.A. Climate change and eastern Africa: A review of impact on major crops. Food Energy Secur. 2015, 4, 110–132. [Google Scholar] [CrossRef]
  12. Fishman, R. More uneven distributions overturn benefits of higher precipitation for crop yields. Environ. Res. Lett. 2016, 11, 24004. [Google Scholar] [CrossRef]
  13. Wanyama, J.; Ssegane, H.; Kisekka, I.; Komakech, A.J.; Banadda, N.; Zziwa, A.; Ebong, T.O.; Mutumba, C.; Kiggundu, N.; Kayizi, R.K.; et al. Irrigation Development in Uganda: Constraints, Lessons Learned, and Future Perspectives. J. Irrig. Drain. Eng. 2017, 143, 04017003. [Google Scholar] [CrossRef] [Green Version]
  14. UBOS. National Ppopulation and Housing Census 2014; Uganda Bureau of Statistics (UBOS): Kampala, Uganda, 2014.
  15. Iwadra, M.; Odirile, P.T.; Parida, B.P.; Moalafhi, D.B. Evaluation of future climate using SDSM and secondary data (TRMM and NCEP) for poorly gauged catchments of Uganda: The case of Aswa catchment. Theor. Appl. Clim. 2018, 137, 2029–2048. [Google Scholar] [CrossRef]
  16. Kenabatho, P.K.; Parida, B.; Moalafhi, D. Evaluation of satellite and simulated rainfall products for hydrological applications in the Notwane catchment, Botswana. Phys. Chem. Earth 2017, 100, 19–30. [Google Scholar] [CrossRef]
  17. Maraun, D. Bias Correcting Climate Change Simulations—A Critical Review. Curr. Clim. Chang. Rep. 2016, 2, 211–220. [Google Scholar] [CrossRef] [Green Version]
  18. Zickfeld, K.; Eby, M.; Weaver, A.J.; Alexander, K.; Crespin, E.; Edwards, N.R.; Eliseev, A.V.; Feulner, G.; Fichefet, T.; Forest, C.E.; et al. Long-Term Climate Change Commitment and Reversibility: An EMIC Intercomparison. J. Clim. 2013, 26, 5782–5809. [Google Scholar] [CrossRef]
  19. Hossain, I.; Rasel, H.M.; Alam Imteaz, M.; Mekanik, F. Long-term seasonal rainfall forecasting using linear and non-linear modelling approaches: A case study for Western Australia. Theor. Appl. Clim. 2019, 132, 131–141. [Google Scholar] [CrossRef]
  20. Abbot, J.; Marohasy, J. Skilful rainfall forecasts from artificial neural networks with long duration series and single-month optimization. Atmos. Res. 2017, 197, 289–299. [Google Scholar] [CrossRef]
  21. Huang, J.; Zhang, J.; Zhang, Z.; Xu, C.-Y.; Wang, B.; Yao, J. Estimation of future precipitation change in the Yangtze River basin by using statistical downscaling method. Stoch. Environ. Res. Risk Assess. 2010, 25, 781–792. [Google Scholar] [CrossRef]
  22. Fotso-Nguemo, T.C.; Vondou, D.A.; Pokam, W.M.; Djomou, Z.Y.; Diallo, I.; Haensler, A.; Tchotchou, L.A.D.; Kamsu-Tamo, P.H.; Gaye, A.T.; Tchawoua, C. On the added value of the regional climate model REMO in the assessment of climate change signal over Central Africa. Clim. Dyn. 2017, 49, 3813–3838. [Google Scholar] [CrossRef]
  23. Mehrotra, R.; Sharma, A. A Multivariate Quantile-Matching Bias Correction Approach with Auto- and Cross-Dependence across Multiple Time Scales: Implications for Downscaling. J. Clim. 2016, 29, 3519–3539. [Google Scholar] [CrossRef]
  24. Moalafhi, D.B.; Sharma, A.; Evans, J.; Mehrotra, R.; Rocheta, E. Impact of bias-corrected reanalysis-derived lateral boundary conditions on WRF simulations. J. Adv. Model. Earth Syst. 2017, 9, 1828–1846. [Google Scholar] [CrossRef] [Green Version]
  25. Sangelantoni, L.; Russo, A.; Gennaretti, F. Impact of bias correction and downscaling through quantile mapping on simulated climate change signal: A case study over Central Italy. Theor. Appl. Clim. 2018, 135, 725–740. [Google Scholar] [CrossRef]
  26. Mahmood, R.; Babel, M.S. Future changes in extreme temperature events using the statistical downscaling model (SDSM) in the trans-boundary region of the Jhelum river basin. Weather. Clim. Extremes 2014, 5, 56–66. [Google Scholar] [CrossRef] [Green Version]
  27. Saha, S.; Moorthi, S.; Wu, X.; Wang, J.; Nadiga, S.; Tripp, P.; Behringer, D.; Hou, Y.-T.; Chuang, H.-Y.; Iredell, M.; et al. The NCEP Climate Forecast System Version 2. J. Clim. 2014, 27, 2185–2208. [Google Scholar] [CrossRef]
  28. Babel, M.S.; Turyatunga, E. Evaluation of climate change impacts and adaptation measures for maize cultivation in the western Uganda agro-ecological zone. Theor. Appl. Clim. 2014, 119, 239–254. [Google Scholar] [CrossRef]
  29. Mekonnen, D.F.; Disse, M. Analyzing the future climate change of Upper Blue Nile River basin using statistical downscaling techniques. Hydrol. Earth Syst. Sci. 2018, 22, 2391–2408. [Google Scholar] [CrossRef] [Green Version]
  30. Hasan, D.S.N.A.B.P.A.; Ratnayake, U.; Shams, S.; Nayan, Z.B.H.; Rahman, E.K.A. Prediction of climate change in Brunei Darussalam using statistical downscaling model. Theor. Appl. Clim. 2017, 133, 343–360. [Google Scholar] [CrossRef]
  31. Amekudzi, L.K.; Osei, M.A.; Atiah, W.A.; Aryee, J.; Ahiataku, M.A.; Quansah, E.; Preko, K.; Danuor, S.K.; Fink, A.H. Validation of TRMM and FEWS Satellite Rainfall Estimates with Rain Gauge Measurement over Ashanti Region, Ghana. Atmos. Clim. Sci. 2016, 6, 500–518. [Google Scholar] [CrossRef] [Green Version]
  32. Amekudzi, L.K.; Yamba, E.I.; Preko, K.; Asare, E.O.; Aryee, J.; Baidu, M.; Codjoe, S.N.A. Variabilities in Rainfall Onset, Cessation and Length of Rainy Season for the Various Agro-Ecological Zones of Ghana. Climate 2015, 3, 416–434. [Google Scholar] [CrossRef]
  33. Odekunle, T.O. Determining rainy season onset and retreat over Nigeria from precipitation amount and number of rainy days. Theor. Appl. Clim. 2005, 83, 193–201. [Google Scholar] [CrossRef]
  34. Odekunle, T.O.; Balogun, E.E.; Ogunkoya, O.O. On the prediction of rainfall onset and retreat dates in Nigeria. Theor. Appl. Clim. 2005, 81, 101–112. [Google Scholar] [CrossRef]
  35. Guenang, G.M.; Kamga, F.M. Onset, retreat and length of the rainy season over Cameroon. Atmos. Sci. Lett. 2012, 13, 120–127. [Google Scholar] [CrossRef]
  36. Ngetich, F.K.; Mucheru-Muna, M.; Mugwe, J.; Shisanya, C.; Diels, J.; Mugendi, D. Length of growing season, rainfall temporal distribution, onset and cessation dates in the Kenyan highlands. Agric. For. Meteorol. 2014, 188, 24–32. [Google Scholar] [CrossRef]
  37. Mugalavai, E.M.; Kipkorir, E.C.; Raes, D.; Rao, M.S. Analysis of rainfall onset, cessation and length of growing season for western Kenya. Agric. For. Meteorol. 2008, 148, 1123–1135. [Google Scholar] [CrossRef]
  38. Akinseye, F.M.; Agele, S.O.; Traore, P.C.S.; Adam, M.; Whitbread, A.M. Evaluation of the onset and length of growing season to define planting date—‘a case study for Mali (West Africa)’. Theor. Appl. Clim. 2015, 124, 973–983. [Google Scholar] [CrossRef] [Green Version]
  39. Odekunle, T.O. Rainfall and the length of the growing season in Nigeria. Int. J. Clim. 2004, 24, 467–479. [Google Scholar] [CrossRef]
  40. Kallis, G. Droughts. Annual Review of Environment and Resources. Annu. Rev. 2008, 33, 85–118. [Google Scholar] [CrossRef]
  41. Tirivarombo, S.; Osupile, D.; Eliasson, P. Drought monitoring and analysis: Standardised Precipitation Evapotranspiration Index (SPEI) and Standardised Precipitation Index (SPI). Phys. Chem. Earth 2018, 106, 1–10. [Google Scholar] [CrossRef]
  42. Kumar, A. Characterizing meteorological drought using standardized precipitation index for Dehradun, Uttarakhand. Indian J. Soil Conserv. 2016, 44, 157–162. [Google Scholar] [CrossRef]
  43. McKee, T.B.; Doesken, N.J.; Kleist, J. The Relationship of Drought Frequency and Duration to Time Scales. In Proceedings of the Eighth Conference on Applied Climatology, Anaheim, CA, USA, 17–22 January 1993. [Google Scholar]
  44. Guenang, G.M.; Kamga, F.M. Computation of the Standardized Precipitation Index (SPI) and Its Use to Assess Drought Occurrences in Cameroon over Recent Decades. J. Appl. Meteorol. Clim. 2014, 53, 2310–2324. [Google Scholar] [CrossRef]
  45. Yusof, F.; Hui-Mean, F.; Suhaila, J.; Yusop, Z.; Ching-Yee, K.; Jamaludin, S.S.S. Rainfall characterisation by application of standardised precipitation index (SPI) in Peninsular Malaysia. Theor. Appl. Clim. 2013, 115, 503–516. [Google Scholar] [CrossRef]
  46. Bari, S.H.; Rahman, T.U.; Hoque, M.A.; Hussain, M. Analysis of seasonal and annual rainfall trends in the northern region of Bangladesh. Atmos. Res. 2016, 148–158. [Google Scholar] [CrossRef]
  47. Shefine, B.G. Analysis of Meteorological Drought Using SPI and Large Scale Climate Variability (ENSO)-A case study in North Shewa Zone, Amhara Regional State, Ethiopia. Hydrol. Curr. Res. 2018, 9. [Google Scholar] [CrossRef]
  48. Karanja, A.; Ondimu, K.; Recha, C. Analysis of Temporal Drought Characteristic Using SPI Drought Index Based on Rainfall Data in Laikipia West Sub-County, Kenya. OALib 2017, 4, 1–11. [Google Scholar] [CrossRef]
  49. Kisaka, M.O.; Mucheru-Muna, M.; Ngetich, F.K.; Mugwe, J.; Mugendi, D.; Mairura, F. Rainfall Variability, Drought Characterization, and Efficacy of Rainfall Data Reconstruction: Case of Eastern Kenya. Adv. Meteorol. 2015, 2015, 1–16. [Google Scholar] [CrossRef] [Green Version]
  50. Nicholson, S.E.; Funk, C.; Fink, A.H. Rainfall over the African continent from the 19th through the 21st century. Glob. Planet. Chang. 2018, 165, 114–127. [Google Scholar] [CrossRef]
  51. Batisani, N.; Yarnal, B. Rainfall variability and trends in semi-arid Botswana: Implications for climate change adaptation policy. Appl. Geogr. 2010, 30, 483–489. [Google Scholar] [CrossRef]
  52. Byakatonda, J.; Parida, B.; Kenabatho, P.; Moalafhi, D. Modeling dryness severity using artificial neural network at the Okavango Delta, Botswana. Glob. Nest. J. 2016, 18, 463–481. [Google Scholar]
  53. Ibn Musah, A.-A.; Du, J.; Udimal, T.B.; Sadick, M.A. The Nexus of Weather Extremes to Agriculture Production Indexes and the Future Risk in Ghana. Climate 2018, 6, 86. [Google Scholar] [CrossRef] [Green Version]
  54. Derbile, E.K.; Kasei, R.A. Vulnerability of crop production to heavy precipitation in north-eastern Ghana. Int. J. Clim. Chang. Strat. Manag. 2012, 4, 36–53. [Google Scholar] [CrossRef]
Figure 1. Location of Aswa River catchment and GIS generated average Tropical Rainfall Measuring Mission (TRMM) annual rainfall (2000–2015) distribution.
Figure 1. Location of Aswa River catchment and GIS generated average Tropical Rainfall Measuring Mission (TRMM) annual rainfall (2000–2015) distribution.
Climate 08 00067 g001
Figure 2. Correlation of TRMM and observed rainfall data for Lira.
Figure 2. Correlation of TRMM and observed rainfall data for Lira.
Climate 08 00067 g002
Figure 3. Simulated future rainfall for Agago, Amuria, Gulu, Agoro and Kaabong locations (a) CanESM2-SDSM (2006–2100) (b) HadCM3-SDSM (2006–2097).
Figure 3. Simulated future rainfall for Agago, Amuria, Gulu, Agoro and Kaabong locations (a) CanESM2-SDSM (2006–2100) (b) HadCM3-SDSM (2006–2097).
Climate 08 00067 g003aClimate 08 00067 g003b
Figure 4. Agago Future Mean Rainfall Onset and Cessation (TRMM, weather generated Rainfall (1975–2001) and projected (HadCM3) 2006–2035, 2006–2065 and 2006–2099).
Figure 4. Agago Future Mean Rainfall Onset and Cessation (TRMM, weather generated Rainfall (1975–2001) and projected (HadCM3) 2006–2035, 2006–2065 and 2006–2099).
Climate 08 00067 g004
Figure 5. Agago Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM Rainfall (1975-2001) and projected (HadCM3) 2006–2035, 2006–2065 and 2006–2099).
Figure 5. Agago Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM Rainfall (1975-2001) and projected (HadCM3) 2006–2035, 2006–2065 and 2006–2099).
Climate 08 00067 g005
Figure 6. Agago Future Mean Rainfall Onset and Cessation (TRMM Rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099).
Figure 6. Agago Future Mean Rainfall Onset and Cessation (TRMM Rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099).
Climate 08 00067 g006
Figure 7. Agago Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM Rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099).
Figure 7. Agago Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM Rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099).
Climate 08 00067 g007
Figure 8. Kaabong Future Mean Rainfall Onset and Cessation Dates (TRMM rainfall (1975–2001) and projected (HadCM3) 2006–2035, 2036–2065 and 2066–2099).
Figure 8. Kaabong Future Mean Rainfall Onset and Cessation Dates (TRMM rainfall (1975–2001) and projected (HadCM3) 2006–2035, 2036–2065 and 2066–2099).
Climate 08 00067 g008
Figure 9. Kaabong Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM rainfall (1975–2001) and projected (HadCM3) 2006–2035, 2036–2065 and 2066–2099).
Figure 9. Kaabong Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM rainfall (1975–2001) and projected (HadCM3) 2006–2035, 2036–2065 and 2066–2099).
Climate 08 00067 g009
Figure 10. Kaabong Future Mean Rainfall Onset and Cessation Dates (TRMM rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099.
Figure 10. Kaabong Future Mean Rainfall Onset and Cessation Dates (TRMM rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099.
Climate 08 00067 g010
Figure 11. Kaabong Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099.
Figure 11. Kaabong Future Mean Rainfall Onset and Cessation for Rainy Days (TRMM rainfall (1975–2005) and projected (CanESM2) 2006–2035, 2036–2065 and 2066–2099.
Climate 08 00067 g011
Figure 12. Agoro CanESM2 future average monthly rainfall variations for (a) maximum monthly rainfall and (b) minimum monthly rainfall.
Figure 12. Agoro CanESM2 future average monthly rainfall variations for (a) maximum monthly rainfall and (b) minimum monthly rainfall.
Climate 08 00067 g012
Figure 13. Gulu CanESM2 future average monthly rainfall variations for (a) maximum monthly rainfall and (b) minimum monthly rainfall.
Figure 13. Gulu CanESM2 future average monthly rainfall variations for (a) maximum monthly rainfall and (b) minimum monthly rainfall.
Climate 08 00067 g013
Figure 14. CanESM2 3-month period SPI values for Kaabong for the duration 2006–2099.
Figure 14. CanESM2 3-month period SPI values for Kaabong for the duration 2006–2099.
Climate 08 00067 g014
Figure 15. CanESM2 3-month period SPI values for Agoro for the duration 2006–2099.
Figure 15. CanESM2 3-month period SPI values for Agoro for the duration 2006–2099.
Climate 08 00067 g015
Figure 16. CanESM2 3-month period SPI values for Gulu for the duration 2006–2099.
Figure 16. CanESM2 3-month period SPI values for Gulu for the duration 2006–2099.
Climate 08 00067 g016
Table 1. TRMM rainfall data location points.
Table 1. TRMM rainfall data location points.
LocationLatitude, °NLongitude, °E
Agoro3.8833.13
Agago2.8833.38
Amuria2.1333.63
Kaabong3.6333.88
Kitgum3.3832.88
Gulu2.6332.63
Lira2.1333.38
Nimule3.6332.13
Source: http://iridl.ldeo.columbia.edu/SOURCES/.NASA/.GES-DAAC/.TRMM_L3/.TRMM_3B42/.v7/%20.daily/ (with a resolution of 0.25° by 0.25° or 32 by 32 Km).
Table 2. Predictor selection for Kitgum Rainfall.
Table 2. Predictor selection for Kitgum Rainfall.
Predictor DescriptionR1P.rPPRP = ((P.r − R1)/R1)P.rPPRP = ((P.r − R1)/R1)P.rPPRP = ((P.r − R1)/R1)P.rPPRP = ((P.r − R1)/R1)Rank
Mean temperature at 2 m0.59 Super Predictor
Surface specific humidity0.540.290.000.460.190.000.640.070.000.870.030.070.94
850 hPa zonal velocity0.520.370.000.29 2nd
Specific humidity at 850 hPa0.500.310.000.390.220.000.55 3rd
Surface zonal velocity0.490.320.000.360.100.000.790.120.000.760.100.000.79
850 hPa airflow strength0.450.310.000.300.070.000.840.090.000.790.080.000.83
Surface airflow strength0.440.260.000.400.090.000.800.110.000.760.090.000.79
Specific humidity at 500 hPa0.430.240.000.440.190.000.550.100.000.760.120.000.745th
precipitation0.370.170.000.530.040.010.890.100.000.72 4th
Key: PRP is the percentage by which partial correlation is lowered in respect to the correlation coefficient, P.r is partial correlation coefficient, R1 is the value of correlation coefficient between predictor and predictand, P is probability value. Values in bold indicate the highest correlation coefficient for RI column or the lowest PRP value in the PRP column.
Table 3. Calibration and validation of Statistical down Scaling Model (SDSM) results for various locations within the Aswa catchment.
Table 3. Calibration and validation of Statistical down Scaling Model (SDSM) results for various locations within the Aswa catchment.
LocationGCMParameterProcessBiasRegression StatisticsANOVA
Rainfall Period RR2Standard ErrorObservationsFSignificance of F
AgagoHadCM3 Calibration1975–1990(16yrs)1.0340.870.7629.251926122.58 × 10−61
Validation 1991–2001(11yrs)0.9880.880.7726.931324304.61 × 10−43
CanESM2 Calibration1975–1990(16yrs)1.0210.970.9513.8419235351.00 × 10−124
Validation1991–2005(15yrs)1.0100.970.9512.6418031672.54 × 10−115
AgoroHadCM3 Calibration1975–1990(16yrs)1.0480.850.7221.031924941.01 × 10−54
Validation 1991–2001(11yrs)0.9960.930.8614.111327837.11 × 10−57
CanESM2 Calibration1975–1990(16yrs)0.9640.950.9012.3319217978.59 × 10−99
Validation1991–2005(15yrs)1.0110.910.8315.431808422.10 × 10−69
AmuriaHadCM3 Calibration1975–1990(16yrs)1.0300.870.7626.821926083.94 × 10−61
Validation 1991–2001(11yrs)1.0550.860.7426.431323772.92 × 10−40
CanESM2 Calibration1975–1990(16yrs)1.0440.940.8819.3619213424.79 × 10−88
Validation1991–2005(15yrs)1.0570.950.9016.1518016447.80 × 10−92
KaabongHadCM3 Calibration1975–1990(16yrs)0.9950.900.8115.441928121.67 × 10−70
Validation 1991–2001(11yrs)1.1360.830.6915.751322889.52 × 10−35
CanESM2 Calibration1975–1990(16yrs)0.9800.950.8911.5919215873.49 × 10−94
Validation1991–2005(15yrs)0.9660.940.899.6218013721.38 × 10−85
GuluHadCM3 Calibration1975–1990(16yrs)1.0190.960.9219.1919222483.12 × 10−107
Validation 1991–2001(11yrs)0.9820.930.8723.211328857.40 × 10−60
CanESM2 Calibration1975–1990(16yrs)1.0190.970.9515.4519235743.73 × 10−125
Validation1991–2005(16yrs)0.9980.960.9218.1818020096.77 × 10−99
Table 4. Classification of SPI values.
Table 4. Classification of SPI values.
SPI ValueClass
2.00 or moreExtremely wet
1.50–1.99Severely wet
1.00–1.49Moderately wet
0–0.99Mildly wet
0 to −0.99Mildly dry
−1.00 to −1.49Moderately dry
−1.50 to −1.99Severely dry
−2 or less Extremely dry
Source: [44].
Table 5. Results of Homogeneity Test and Mann–Kendall Statistics for TRMM and Simulated Rainfall Data for Agoro and Amuria Locations in Aswa Catchment. SNHT = standard normal homogeneity test.
Table 5. Results of Homogeneity Test and Mann–Kendall Statistics for TRMM and Simulated Rainfall Data for Agoro and Amuria Locations in Aswa Catchment. SNHT = standard normal homogeneity test.
Homogeneity Test-SNHTMann–Kendall Statistics
GCM,
Location
Data TypePeriod p-Value (Two-Tailed)AlphaCommentNumber of ObservationsKendall’s TauSen’s SlopeComment
CanESM2
AgoroTRMM and generated monthly Rainfall1975–20050.0990.05p > α (0.05), homogeneous372
TRMM and simulated annual Rainfall2006–2099 940.68318.52Positive trend
TRMM and simulated annual Rainfall2006–2065 600.44910.42Positive trend
TRMM and simulated annual Rainfall2006–2035 300.2006.628Positive trend
HadCM3
AgoroTRMM and simulated annual Rainfall2006–2097 920.6693.219Positive trend
TRMM and simulated annual Rainfall2006–2065 600.4632.339Positive trend
TRMM and simulated annual Rainfall2006–2035 300.1951.426Positive trend
CanESM2
AmuriaTRMM and generated monthly Rainfall1975–20050.2200.05p > α (0.05), homogeneous372
TRMM and simulated annual Rainfall2006–2100 940.64710.52Positive trend
TRMM and simulated annual Rainfall2006–2065 600.4729.072Positive trend
TRMM and simulated annual Rainfall2006–2035 300.1726.286Positive trend
HadCM3
AmuriaTRMM and simulated annual Rainfall2001–2097 920.64718.54Positive trend
TRMM and simulated annual Rainfall2001–2065 600.48213.70Positive trend
TRMM and simulated annual Rainfall2001–2035 300.2189.212Positive trend
Table 6. Agoro TRMM and Simulated (CanESM2and Had3) Monthly Rainfall Coefficient of Variation.
Table 6. Agoro TRMM and Simulated (CanESM2and Had3) Monthly Rainfall Coefficient of Variation.
Jan.Feb.Mar.Apr.MayJun.Jul.Aug.Sep.Oct.Nov.Dec.
CanESM2
1975–2005 TRMM and Weather Generator
Mean24.228.448.2116.974.183.2108.9117.355.983.865.925.0
Standard Deviation18.423.010.624.118.813.711.123.613.415.528.110.0
CV, %76.180.822.020.625.316.510.220.124.018.542.740.1
2006–2035 TRMM and Simulation
Mean21.558.084.0164.2123.679.956.2191.4159.456.6166.352.9
Standard Deviation7.255.838.652.937.423.816.430.051.827.892.720.6
CV33.496.346.032.230.329.729.115.732.549.055.738.9
2006–2065 TRMM and Simulation
Mean24.296.6107.2209.9145.897.247.7195.1197.958.1160.849.8
Standard Deviation8.3110.372.580.640.236.115.733.374.068.891.218.5
CV34.4114.167.738.427.637.132.817.137.4118.556.737.1
2006–2099 TRMM and Simulation
Mean28.7224.0150.0246.5176.3128.840.4206.1271.754.4171.047.9
Standard Deviation14.2255.494.9103.559.960.916.634.5126.762.1101.516.6
CV, %49.4114.063.342.034.047.341.116.746.6114.059.434.7
HadCM3
1975–2001 TRMM and Weather Generator
Mean46.638.649.162.774.173.6101.8101.092.986.374.755.8
Standard Deviation26.521.532.044.842.430.334.726.628.223.835.333.0
CV, %56.955.765.171.357.241.134.026.330.327.647.259.1
2006–2035 TRMM and Simulation
Mean46.638.649.162.774.173.6101.8101.092.986.374.755.8
Standard Deviation26.521.532.044.842.430.334.726.628.223.835.333.0
CV, %56.955.765.171.357.241.134.026.330.327.647.259.1
2006–2065 TRMM ans Simulation
Mean66.973.783.783.889.276.978.969.661.967.372.067.7
Standard Deviation39.753.354.449.045.641.844.342.940.942.445.340.0
CV, %59.472.365.058.551.154.456.161.666.063.062.959.2
2006–2097 TRMM and Simulation
Mean77.767.267.465.676.374.680.585.489.589.890.483.6
Standard Deviation57.952.753.650.752.047.246.057.166.361.564.755.3
CV, %74.578.479.577.368.163.357.166.974.168.571.666.1
The red and black items show the high CV values. The red color gives the highest CV values for rainfall onset and cessation periods.
Table 7. Kaabong TRMM and Simulated (CanESM2and HadCM3) Monthly Rainfall Coefficient of Variation.
Table 7. Kaabong TRMM and Simulated (CanESM2and HadCM3) Monthly Rainfall Coefficient of Variation.
Jan.Feb.Mar.Apr.MayJun.Jul.Aug.Sep.Oct.Nov.Dec.
CanESM2
1975–2005 TRMM and Weather Generator
Mean9.912.447.576.290.451.563.967.028.840.543.924.4
Standard Deviation6.44.616.115.535.514.07.519.69.513.717.616.4
CV, %65.136.833.820.339.327.111.829.333.033.740.167.5
2006–2035 TRMM and Simulation
Mean7.321.070.392.8172.841.7111.161.624.812.258.134.2
Standard Deviation1.29.623.518.541.811.723.119.65.56.615.515.9
CV, %16.245.733.420.024.228.220.831.822.254.126.746.5
2006–2065 TRMM ans Simulation
Mean7.337.277.7101.0239.239.2133.555.124.28.376.338.0
Standard Deviation1.226.132.324.087.810.935.018.85.16.732.117.7
CV, %15.770.141.623.736.727.926.234.121.081.242.146.5
2006–2099 TRMM and Simulation
Mean7.264.591.2114.1322.534.5172.851.623.65.8107.743.2
Standard Deviation1.848.338.531.8139.112.166.318.25.16.567.521.2
CV, %25.274.942.327.843.135.138.435.321.6111.862.649.1
HadCM3
1975–2001 TRMM and Weather Generator
Mean11.214.668.785.3102.656.052.166.425.451.162.827.5
Standard Deviation3.74.320.09.130.510.25.35.43.06.712.612.3
CV, %32.929.229.210.629.718.310.38.111.813.220.044.8
2006–2035 TRMM and Simulation
Mean49.151.266.077.187.690.7102.696.480.474.573.962.4
Standard Deviation25.839.048.662.060.440.334.341.532.430.323.532.8
CV, %52.576.273.780.469.044.433.543.140.440.631.852.6
2006–2065 TRMM and Simulation
Mean88.486.693.988.3103.1101.097.892.294.295.2100.997.8
Standard Deviation58.258.453.952.163.138.535.846.075.771.969.172.2
CV, %65.967.557.459.061.338.136.649.980.375.668.573.8
2006–2097 TRMM and Simulation
Mean120.8109.8127.7127.5135.1135.8130.1128.3120.0117.3123.6132.8
Standard Deviation79.178.7100.2109.696.3115.783.781.080.675.387.689.9
CV, %65.571.778.586.071.385.264.363.167.264.270.867.7
The red and black items show the high CV values. The red color gives the highest CV values for rainfall onset and cessation periods.
Table 8. Gulu TRMM and Simulated (CanESM2and HadCM3) Monthly Rainfall Coefficient of Variation.
Table 8. Gulu TRMM and Simulated (CanESM2and HadCM3) Monthly Rainfall Coefficient of Variation.
Jan.Feb.Mar.Apr.MayJun.Jul.Aug.Sep.Oct.Nov.Dec.
CanESM2
2006–2099 TRMM and Simulation
Mean82.4141.7297.9314.1253.4153.8220.3380.4323.2278.9137.442.8
Standard Deviation86.584.6116.363.463.034.1117.7109.261.288.451.626.6
CV, %105.059.739.020.224.922.253.428.718.931.737.662.3
2006–2065 TRMM and Simulation
Mean44.299.9237.5307.3225.8154.5144.9322.1334.5243.2144.636.2
Standard Deviation41.954.774.764.042.534.946.087.465.069.249.214.8
CV, %94.754.831.420.818.822.631.727.119.428.534.041.0
2006–2035 TRMM and Simulation
Mean27.175.3196.6289.3211.9153.7114.2263.1331.2224.9146.434.3
Standard Deviation14.842.748.571.147.835.618.867.358.363.748.613.1
CV, %54.756.624.724.622.623.116.525.617.628.333.238.3
1975–2005 TRMM and Weather Generator
28.432.7127.7198.0200.4108.6150.5172.9149.1160.679.431.8
Standard Deviation21.920.135.350.535.728.135.438.730.526.629.911.8
CV, %77.261.527.725.517.825.923.522.420.416.637.637.1
HadCM3
2006–2097 TRMM and Simulation
Mean140.5160.5227.2293.8180.5114.4345.7179.0332.0171.154.2105.3
Standard Deviation55.165.273.4101.663.929.2203.774.3119.645.027.542.5
CV, %39.240.632.334.635.425.558.941.536.026.350.840.4
2006–2065 TRMM and Simulation
Mean120.0148.6228.0269.5195.9125.0229.1186.0261.4180.256.983.3
Standard Deviation44.263.775.080.349.925.859.386.780.849.632.732.8
CV, %36.942.932.929.825.520.625.946.630.927.557.439.4
2006–2035 TRMM and Simulation
Mean106.4108.1172.7220.2177.7121.9218.1238.8231.8213.383.376.9
Standard Deviation50.750.558.479.761.325.227.639.941.227.223.434.7
CV, %47.646.733.836.234.520.712.616.717.812.728.145.1
1975–2005 TRMM and Weather Generator
Mean26.437.1132.5222.1208.3112.4170.7175.8158.9161.374.533.8
Standard Deviation19.613.727.031.121.619.59.815.919.418.415.47.9
CV, %74.336.920.414.010.317.35.89.012.211.420.723.3
The red and black items show the high CV values. The red color gives the highest CV values for rainfall onset and cessation periods.
Table 9. Percentage variation of Standardized Precipitation Index (SPI) values for Kaabong over the period 2006–2099.
Table 9. Percentage variation of Standardized Precipitation Index (SPI) values for Kaabong over the period 2006–2099.
SPITime Period
ValueClass2006–20352036–20652066–2099
Percentage of period
≥2Extremely Wet0.30.35.4
1.5 to 1.99Severely Wet0.82.511.0
1.0 to 1.49Moderately Wet0.82.822.8
0.99 to −0.99Mildly wet/dry (near normal)63.483.655.4
−1.0 to −1.49Moderately Dry19.05.82.5
−1.5 to −2.0Severely Dry12.83.62.2
≤−2.0Extremely Dry2.81.40.7
Maximum SPI Value 2.22.53.0
Minimum SPI Value −2.9−2.5−2.8
Table 10. Percentage variation of SPI values for Agoro over the period 2006–2099.
Table 10. Percentage variation of SPI values for Agoro over the period 2006–2099.
SPITime Period
Value Class2006–20352036–20652066–2099
Percentage of period
≥2Extremely Wet0.00.86.4
1.5 to 1.99Severely Wet1.10.69.8
1.0 to 1.49Moderately Wet2.04.424.8
0.99 to −0.99Mildly wet/dry (near normal)62.080.856.4
−1.0 to −1.49Moderately Dry22.96.91.5
−1.5 to −2.0Severely Dry10.14.40.7
≤−2.0Extremely Dry2.01.90.5
Maximum SPI Value 1.873.483.0
Minimum SPI Value −2.48−2.9−2.28
Table 11. Percentage variation of SPI values for Gulu over the period 2006–2099.
Table 11. Percentage variation of SPI values for Gulu over the period 2006–2099.
SPITime Period
ValueClass2006–20352036–20652066–2099
Percentage of period
≥2Extremely Wet0.00.11.9
1.5 to 1.99Severely Wet0.30.44.6
1.0 to 1.49Moderately Wet0.83.110.5
0.99 to −0.99Mildly wet/dry (near normal)64.573.367.3
−1.0 to −1.49Moderately Dry20.414.39.8
−1.5 to −2.0Severely Dry10.16.34.3
≤−2.0Extremely Dry3.92.51.7
Maximum SPI Value 1.92.52.8
Minimum SPI Value −3.3−3.3−3.3

Share and Cite

MDPI and ACS Style

Iwadra, M.; Odirile, P.T.; Parida, B.P.; Moalafhi, D.B. An Investigation into the Future Changes in Onset and Cessation of Rain and Their Variability over the Aswa Catchment, Uganda. Climate 2020, 8, 67. https://doi.org/10.3390/cli8060067

AMA Style

Iwadra M, Odirile PT, Parida BP, Moalafhi DB. An Investigation into the Future Changes in Onset and Cessation of Rain and Their Variability over the Aswa Catchment, Uganda. Climate. 2020; 8(6):67. https://doi.org/10.3390/cli8060067

Chicago/Turabian Style

Iwadra, Michael, P. T. Odirile, B. P. Parida, and D. B. Moalafhi. 2020. "An Investigation into the Future Changes in Onset and Cessation of Rain and Their Variability over the Aswa Catchment, Uganda" Climate 8, no. 6: 67. https://doi.org/10.3390/cli8060067

APA Style

Iwadra, M., Odirile, P. T., Parida, B. P., & Moalafhi, D. B. (2020). An Investigation into the Future Changes in Onset and Cessation of Rain and Their Variability over the Aswa Catchment, Uganda. Climate, 8(6), 67. https://doi.org/10.3390/cli8060067

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop