Statistical Analysis of Recent and Future Rainfall and Temperature Variability in the Mono River Watershed (Benin, Togo)

This paper assessed the current and mid-century trends in rainfall and temperature over the Mono River watershed. It considered observation data for the period 1981–2010 and projection data from the regional climate model (RCM), REMO, for the period 2018–2050 under emission scenarios RCP4.5 and RCP8.5. Rainfall data were interpolated using ordinary kriging. Mann-Kendall, Pettitt and Standardized Normal Homogeneity (SNH) tests were used for trends and break-points detection. Rainfall interannual variability analysis was based on standardized precipitation index (SPI), whereas anomalies indices were considered for temperature. Results revealed that on an annual scale and all over the watershed, temperature and rainfall showed an increasing trend during the observation period. By 2050, both scenarios projected an increase in temperature compared to the baseline period 1981–2010, whereas annual rainfall will be characterized by high variabilities. Rainfall seasonal cycle is expected to change in the watershed: In the south, the second rainfall peak, which usually occurs in September, will be extended to October with a higher value. In the central and northern parts, rainfall regime is projected to be characterized by late onsets, a peak in September and lower precipitation until June and higher thereafter. The highest increase and decrease in monthly precipitation are expected in the northern part of the watershed. Therefore, identifying relevant adaptation strategies is recommended.


Introduction
Modifications in the climate as a result of both natural and anthropogenic processes have raised considerable concerns (such as more frequent and intense rainfall, droughts, dry spells, violent winds, etc.), as they induce adverse impacts on several development sectors [1].
In recent decades, weather and climate extremes such as droughts, heat waves, wild fires, floods and storms have increased in frequency and intensity in several regions of the world. In fact, Vincent et al. [2] noticed that the percentage of warm nights is increasing while that of cold nights is decreasing in South-America. In addition, the US Climate Change Science Program underscored the fact that heavy precipitations have become more frequent and intense in Northern America [3]. In central Asia, Savitskaya [4] reported that, during the last 50 years, there was high variability in the pattern of precipitation, whereas winter has become warmer in the entire region.

Data Used
Both observation and projection data were used. Daily observed rainfall and temperature were collected from meteorological institutes of Benin and Togo (DMN, Direction de la Météorologie Nationale) for the period 1981-2010. That period is the current normal used for climatological analysis, and this study aims at taking it into account, as previous studies have already accounted for other normals [15,17,18]. As presented by Figure 1, rainfall data were collected at 24 rain gauges within and around Mono watershed (not farther than 25 km). As for temperature data, they were collected from three synoptic stations located in the watershed: Tabligbo, Atakpamé and Sokodé.
Projected data were provided by the regional climate model (RCM), REMO (Table 1). Rainfall and temperature data were extracted for the period 2018-2050 and under two representative concentration pathways: RCP4.5 (intermediate pathway) and RCP8.5 (most extreme pathway). Akinsanola et al. [19] have already noticed that REMO fairly simulates rainfall in West Africa and concluded that it can be used for future climate projections in the region. However, the use of a multimodel ensemble approach (more than one model) is recommended in order to better estimate the actual climate and improve the robustness of climate change projections [20][21][22].

Data Used
Both observation and projection data were used. Daily observed rainfall and temperature were collected from meteorological institutes of Benin and Togo (DMN, Direction de la Météorologie Nationale) for the period 1981-2010. That period is the current normal used for climatological analysis, and this study aims at taking it into account, as previous studies have already accounted for other normals [15,17,18]. As presented by Figure 1, rainfall data were collected at 24 rain gauges within and around Mono watershed (not farther than 25 km). As for temperature data, they were collected from three synoptic stations located in the watershed: Tabligbo, Atakpamé and Sokodé.
Projected data were provided by the regional climate model (RCM), REMO (Table 1). Rainfall and temperature data were extracted for the period 2018-2050 and under two representative concentration pathways: RCP4.5 (intermediate pathway) and RCP8.5 (most extreme pathway). Akinsanola et al. [19] have already noticed that REMO fairly simulates rainfall in West Africa and concluded that it can be used for future climate projections in the region. However, the use of a multi-model ensemble approach (more than one model) is recommended in order to better estimate the actual climate and improve the robustness of climate change projections [20][21][22]. These data were accessed online (https://www.cordex.org) in the context of the Coordinated Regional Climate Downscaling Experiment (CORDEX) over Africa at 0.44 • resolution.

Trend Analysis
Trends and breakpoints were assessed using the non-parametric test of Mann-Kendall first, followed by the Pettitt and SNH tests The Mann-Kendall test is used in order to establish whether there is a trend (increasing or decreasing) in the time series. It is done with a confidence level of 95%, and the hypotheses are: H 0 : there is no trend in rainfall time series; H 1 : there is a trend in rainfall time series.
For both Pettitt and SNH tests, significance level α = 0.05, and the hypotheses are: H 0 : there is no change in annual rainfall data; H 1 : there is a date at which there is a change in the data.
Moreover, in order to assess the trend of annual rainfall at watershed scale, and because rainfall data are not measured in every single grid of the watershed, spatial interpolation was required.
In the scope of this study, ordinary kriging (OK) was chosen over other methods-such as arithmetic mean, Thiessen polygon, inverse distance weighting-because (i) it takes into account not only the distance between observation stations and estimation point but also the distance between stations taken two by two; (ii) it is a stochastic method which provides the best linear unbiased predictions; (iii) the interpolation error can be estimated [23]. Nonetheless, one of the limitations of kriging is that it is not suitable when there are few observation points. Kriging was basically developed for geostatistics purposes [24] but is widely used in climatology. It is worth noting that the 'backbone' of kriging is the variogram which explains the variance of the studied variable with respect to distance between observation points. Equation (1) presents the formula of variogram where G(h) is the variogram, N(h) the number of coupled points separated by the distance h, z(p i ) the observed rainfall at location p i , and z(p i + h) the observed rainfall at location p i + h. Furthermore, considering the fact that rainfall regime in Mono watershed is not homogenous, rainfall trend analysis is carried out with respect to three latitude-based regions, as done in previous studies [15,18]. The regions are defined as follows: latitude < 7, 7 ≤ latitude ≤ 8 and latitude > 8. Hereinafter, these regions are respectively referred to as southern part, central part and northern part of the Mono watershed.
Analysis of temperature trends over the watershed was conducted on the arithmetic mean from the stations of Tabligbo, Sokodé and Atakpamé.

Interannual Rainfall Variability Analysis
The standardized precipitation index (SPI) [25] is a tool recommended by the World Meteorological Organization (WMO) and widely used for quantifying the precipitation deficit over different timescales (3 to 48 months). For the selected timescale, rainfall records are fitted with a probability distribution which is then transformed into a normal distribution so that the mean SPI for the location and desired period is zero. Hence, this method improves the common anomaly method (Equation (2)), which does not take into account the fact that rainfall is typically not normally distributed for a cumulative period of 12 months or less.
where I(i), x i , x m and σ are respectively the standardized index of year i, the value for the year i, the average and the standard deviation of the time series.
In the present study, the SPI 12 (for 12 month's timescale) is used to assess rainfall deficit or excess on a yearly basis. Moreover, SPI 12 is the one recommended for watershed analysis [26]. Table 2 presents the guidelines for analyzing SPI values [25,26]. Raw outputs from RCMs must be corrected prior to local impact studies because of the bias they encompass. There are several bias correction methods but in this study, the methods of delta, linear scaling and empirical quantile mapping (EQM) are used because they have produced satisfactory results in previous studies carried out in similar climatic regions [27][28][29]. The results of Ntcha M'po et al. [27], Essou and Brissette [30] and Speth et al. [31], who bias-corrected REMO data in the Ouémé watershed (Benin), guided the choice of correction methods in this study. Rainfall was corrected with delta method in the south, multiplicative scaling in the central part and EQM in the north. As for temperature data, they were corrected using only EQM.

Percentage of Relative Change in Rainfall Seasonal Cycle
It is computed using Equation (3).
where, P c,m , R proj,m and R obs,m are respectively the percentage of change for mth month, average projected rainfall for month m, and the average rainfall of month m during observation period.

Present Rainfall Pattern
The pattern of monthly rainfall varies from one part of the watershed to the other. The southern and northern parts are respectively characterized by bimodal and unimodal rainfall regime, whereas, a transitory (neither bimodal nor unimodal) regime was found in the central region ( Figure 2). In the south, two peaks are respectively recorded in June and September, while the unique peak in the northern part occurs in August. On the other side, the rainy season in the central part lasts from March to September. These results are in line with previous research findings [15,17,18]. Thus, the bimodal or unimodal characteristics of the precipitation regime seem not to change over the recent periods compared to the historical period.
The results of the Mann-Kendall test (Table 3) underscored the fact that rainfall in the three regions of Mono watershed had an increasing trend during the period 1981 to 2010.
As presented in Table 3 and Figure 3, the SNH test detected break points in the time series while the Pettitt test did not. The change noted in the center and the north may be related to the well-known 1970s and 1980s droughts which affected many West-African countries and was documented by several authors such as Le Barbé and Lebel [32], Le Barbé et al. [33] and Le lay and Galle [34]. In the south, two peaks are respectively recorded in June and September, while the unique peak in the northern part occurs in August. On the other side, the rainy season in the central part lasts from March to September. These results are in line with previous research findings [15,17,18]. Thus, the bimodal or unimodal characteristics of the precipitation regime seem not to change over the recent periods compared to the historical period.
The results of the Mann-Kendall test (Table 3) underscored the fact that rainfall in the three regions of Mono watershed had an increasing trend during the period 1981 to 2010. As presented in Table 3 and Figure 3, the SNH test detected break points in the time series while the Pettitt test did not. The change noted in the center and the north may be related to the well-known 1970s and 1980s droughts which affected many West-African countries and was documented by several authors such as Le Barbé and Lebel [32], Le Barbé et al. [33] and Le lay and Galle [34].  According to the baseline period, the longest dry period is 1981-1986 in the south, 1981-1983 in the central part and 1981-1984 in the north. In addition, the driest year is 1992 in the south, and 1983 for both the center and north. As for years of highest excess, it is 2010 in the south and 1995 in central and northern parts. It is worth noting that, in the three regions, the longest dry period falls in the 1980s drought events. After this specific period no extreme drought occurred (except in 1992 in the south), and there were more years above normal than below. Overall, the period after 1990 is characterized by more wet years, and it explains the trends highlighted by statistical tests performed above. Similar results have been reported in other watersheds in West Africa by several authors, such as Adeyeri et al. [35] in Komadugu-Yobe basin, Nicholson et al. [36] over West Africa, Ozer et al. [37] over the Sahelian region, and by Lawin [38] and Attogouinon et al. [39] in the upper Ouémé river valley. However, this shift to wetter condition is region dependent, because other studies reported a decreasing trend in rainfall patterns over West-Africa [40,41].  According to the baseline period, the longest dry period is 1981-1986 in the south, 1981-1983 in the central part and 1981-1984 in the north. In addition, the driest year is 1992 in the south, and 1983 for both the center and north. As for years of highest excess, it is 2010 in the south and 1995 in central and northern parts. It is worth noting that, in the three regions, the longest dry period falls in the 1980s drought events. After this specific period no extreme drought occurred (except in 1992 in the south), and there were more years above normal than below. Overall, the period after 1990 is characterized by more wet years, and it explains the trends highlighted by statistical tests performed above. Similar results have been reported in other watersheds in West Africa by several authors, such as Adeyeri et al. [35] in Komadugu-Yobe basin, Nicholson et al. [36] over West Africa, Ozer et al. [37] over the Sahelian region, and by Lawin [38] and Attogouinon et al. [39] in the upper Ouémé river valley. However, this shift to wetter condition is region dependent, because other studies reported a decreasing trend in rainfall patterns over West-Africa [40,41].

Future Rainfall Pattern
The results of break-point detection and the Mann-Kendall test performed on annual rainfall under RCP4.5 and RCP8.5 are summarized in Table 4.

Future Rainfall Pattern
The results of break-point detection and the Mann-Kendall test performed on annual rainfall under RCP4.5 and RCP8.5 are summarized in Table 4.
The results of homogeneity tests and the Mann-Kendall test suggest that, all over the watershed, there is neither break-point nor a linear trend in annual rainfall time series, under emission scenarios RCP4.5 and RCP8.5. Statistically, rainfall time series are homogenous and present no trend. Nonetheless, some variabilities are observed. Specifically, in the northern part and for RCP8.5, rainfall might decrease by 2035 and increase thereafter. The increase in annual precipitation over recent decades in Mono river watershed seems not to be maintained in the future. Future pattern of rainfall may be marked by high variabilities. Such an absence of significant trend in rainfall is also reported by N'Tcha M'Po et al. [40] in Ouémé river basin by 2050, using a REMO model. Similarly, Lawin et al. [42] reported no trend in rainfall pattern in the Imbo north plain region in Burundi under central Africa climatology using an ensemble of eight regional climate models. The SPI computed for each region emphasizes this variability ( Figure 5). Under RCP4.5, the number of projected deficit years in the watershed increases slightly from south to north. The years 2020 and 2024 are projected to be extremely wet in the south, whereas 2031 is expected to be extremely dry in the north. However, under RCP8.5, the projected number of deficit years decreases from south to north. The year 2033 is projected to be extremely dry in the north, and again, years 2020 and 2024 for this scenario are expected to be extremely wet in the south. The agreement of both scenarios for years 2020 and 2024 shows that those years will potentially be characterized by extreme precipitations. Furthermore, from south to north, RCP8.5 projects more extremely wet years than RCP4.5.
Furthermore, the pattern of the seasonal cycle of rainfall is expected to undergo some modifications. In the southern part (Figure 6a), rainfall seasonal cycle is projected to keep its bimodal pattern under RCP4.5 and RCP8.5. In addition, both scenarios project almost the same pattern. As in the normal period, the first peak is recorded in June but with a slightly lower amount. The second peak, which normally occurs in September, is expected to extend to October with a higher value.
In the central and northern parts (Figure 6b,c), and under RCP4.5 and RCP8.5, the rainfall regime is projected to be characterized by late onsets and lower precipitation until June, compared to observations, and higher thereafter. Both scenarios converge on the fact that rainfall peak will probably occur in September. The northern part is expected to keep its unimodal pattern under both scenarios, whereas a shift from a transitional regime to a unimodal one is expected in the central region.
Projected future rainfall regimes under the two scenarios are quite similar, apart from in August, where RCP4.5 predicted a slightly larger amount than RCP8.5, and an inverse proportion was predicted in June Moreover, on a monthly scale, Figure 7 depicts how the rainfall seasonal cycle is expected to change under RCP4.5 and RCP8.5 compared to the baseline period.
Under RCP4.5, the relative change in monthly rainfall varies from −5.5% to 8.4% in the southern part, −29.9% to 22.2% in the central part and −39% to 91.4% in the northern part. For RCP8.5, the expected change ranges from −3.5% to 5.8% in the southern part, −55% to 20% in the central part and −64% to 85.9% in the northern part. Therefore, the biggest changes (both increase and decrease) in monthly rainfall are expected in the northern part of the watershed, regardless of the scenario considered. According to RCP4.5 and with respect to the observation period, the highest rainfall decrease during the period 2018-2050 is expected to occur in February, whereas January and November will record the highest increases. rainfall may be marked by high variabilities. Such an absence of significant trend in rainfall is also reported by N'Tcha M'Po et al. [40] in Ouémé river basin by 2050, using a REMO model. Similarly, Lawin et al. [42] reported no trend in rainfall pattern in the Imbo north plain region in Burundi under central Africa climatology using an ensemble of eight regional climate models.
The SPI computed for each region emphasizes this variability ( Figure 5). Under RCP4.5, the number of projected deficit years in the watershed increases slightly from south to north. The years 2020 and 2024 are projected to be extremely wet in the south, whereas 2031 is expected to be extremely dry in the north. However, under RCP8.5, the projected number of deficit years decreases from south to north. The year 2033 is projected to be extremely dry in the north, and again, years 2020 and 2024 for this scenario are expected to be extremely wet in the south. The agreement of both scenarios for years 2020 and 2024 shows that those years will potentially be characterized by extreme precipitations. Furthermore, from south to north, RCP8.5 projects more  Under RCP8.5, the highest decrease is projected to affect rainfall in February, whereas the highest increase is expected in November. Thus, the two scenarios project a more than 70% increase of rainfall in the month of November by 2050, compared to the observation period. In addition, the highest increase is projected by RCP4.5 and the highest decrease by RCP8.5. Globally substantial changes are expected prior to and at the end of rainy seasons.  Under RCP8.5, the highest decrease is projected to affect rainfall in February, whereas the highest increase is expected in November. Thus, the two scenarios project a more than 70% increase of rainfall in the month of November by 2050, compared to the observation period. In addition, the highest increase is projected by RCP4.5 and the highest decrease by RCP8.5. Globally substantial changes are expected prior to and at the end of rainy seasons. Under RCP8.5, the highest decrease is projected to affect rainfall in February, whereas the highest increase is expected in November. Thus, the two scenarios project a more than 70% increase of rainfall in the month of November by 2050, compared to the observation period. In addition, the highest increase is projected by RCP4.5 and the highest decrease by RCP8.5. Globally substantial changes are expected prior to and at the end of rainy seasons.

Present Temperature Change
The two homogeneity tests performed on annual temperature revealed the presence of break-point in the time series (Figure 8).

Present Temperature Change
The two homogeneity tests performed on annual temperature revealed the presence of breakpoint in the time series (Figure 8). The Pettitt test indicates that since 1996 mean annual temperature has increased by 0.51 °C in Mono watershed compared to the period 1981-1995, whereas the results of the SNH test implies an increase of 0.52 °C from the period 1981-1996 to 1997-2010. In addition, the Mann-Kendall test suggested a significant increasing trend (p-value = 3.457 × 10 −6 and τ = 0.6).
Anomalies computation revealed that, from 1981 to 1997, temperature was globally below and near normal, but since 1998 it has stayed above the normal (Figure 9).  The Pettitt test indicates that since 1996 mean annual temperature has increased by 0.51 • C in Mono watershed compared to the period 1981-1995, whereas the results of the SNH test implies an increase of 0.52 • C from the period 1981-1996 to 1997-2010. In addition, the Mann-Kendall test suggested a significant increasing trend (p-value = 3.457 × 10 −6 and τ = 0.6).
Anomalies computation revealed that, from 1981 to 1997, temperature was globally below and near normal, but since 1998 it has stayed above the normal (Figure 9).

Present Temperature Change
The two homogeneity tests performed on annual temperature revealed the presence of breakpoint in the time series (Figure 8). The Pettitt test indicates that since 1996 mean annual temperature has increased by 0.51 °C in Mono watershed compared to the period 1981-1995, whereas the results of the SNH test implies an increase of 0.52 °C from the period 1981-1996 to 1997-2010. In addition, the Mann-Kendall test suggested a significant increasing trend (p-value = 3.457 × 10 −6 and τ = 0.6).
Anomalies computation revealed that, from 1981 to 1997, temperature was globally below and near normal, but since 1998 it has stayed above the normal (Figure 9).  It therefore corroborates the outputs from homogeneity tests. In addition, these results are in line with previous studies which noted similar increasing trend of temperature in West Africa (Badjana [43] in the Kara river basin of Togo, Kabo-Bah et al. [44] in the Ghana part of Volta river basin, Oguntunde et al. [45] at Ibadan and Collins [5] over the West-African region).

Future Temperature Change
Under RCP4.5, the homogeneity tests detected break-points at different dates, but they got an agreement under RCP8.5 ( Figure 10 It therefore corroborates the outputs from homogeneity tests. In addition, these results are in line with previous studies which noted similar increasing trend of temperature in West Africa (Badjana [43] in the Kara river basin of Togo, Kabo-Bah et al. [44] in the Ghana part of Volta river basin, Oguntunde et al. [45] at Ibadan and Collins [5] over the West-African region).

Future Temperature Change
Under RCP4.5, the homogeneity tests detected break-points at different dates, but they got an agreement under RCP8.5 ( Figure 10  Regardless of the scenario used, an overall significant increasing trend in temperature is expected by 2050 (Table 5). Regardless of the scenario used, an overall significant increasing trend in temperature is expected by 2050 (Table 5). Figure 11 depicts anomalies of temperature under RCP4.5 and RCP8.5. For both emission scenarios, the number of years above normal is higher compared to the number of years below normal. Therefore, future climate in Mono river watershed is projected to be warmer by 2050. Such an increasing trend by 2050 has been reported by Nelson et al. [46] in Togo. Similarly, Oyerinde [47] reported a consistently increasing trend in temperature over the Niger River Basin, 5% to 10% under RCP4.5 and 5% to 20% under RCP8.5, using an ensemble model from eight regional climate models. In the Massili basin of Burkina Faso, Bontogho [29] reported an increase in temperature by 1.8 • C (RCP4.5) and 3.0 • C (RCP8.5) from 1971 to 2050 using the regional model HIRHAM5. Badou [48] reported a temperature increase of up to 0.48 • C under RCP4.5 and up to 0.45 • C under RCP8.5 using the REMO model in the Benin part of Niger River basin. Overall, all models and scenarios considered by several authors converge to a moderate to high increase of temperature all over the world.   For both emission scenarios, the number of years above normal is higher compared to the number of years below normal. Therefore, future climate in Mono river watershed is projected to be warmer by 2050. Such an increasing trend by 2050 has been reported by Nelson et al. [46] in Togo. Similarly, Oyerinde [47] reported a consistently increasing trend in temperature over the Niger River Basin, 5% to 10% under RCP4.5 and 5% to 20% under RCP8.5, using an ensemble model from eight regional climate models. In the Massili basin of Burkina Faso, Bontogho [29] reported an increase in temperature by 1.8 °C (RCP4.5) and 3.0 °C (RCP8.5) from 1971 to 2050 using the regional model HIRHAM5. Badou [48] reported a temperature increase of up to 0.48 °C under RCP4.5 and up to 0.45 °C under RCP8.5 using the REMO model in the Benin part of Niger River basin. Overall, all models and scenarios considered by several authors converge to a moderate to high increase of temperature all over the world.

Conclusions
This paper examined the trend in annual rainfall and annual temperature of Mono river watershed over the observation period 1981-2010 and by 2050 using the regional model REMO under RCP4.5 and RCP8.5. It also assessed the monthly pattern of rainfall over the same periods. During the last three decades, rainfall and temperature have been increasing all over the Mono River watershed. By 2050 and under emission scenarios RCP4.5 and RCP8.5, annual rainfall is projected to be characterized by high variability, whereas a significant increasing trend is projected for annual temperature (warmer future climate). For each of the three defined regions (south, center and north), and under both emission scenarios, the seasonal cycle of rainfall is expected to change: In the southern part, the first peak is projected to reduce slightly, whereas the second peak is expected to increase and shift to October; however, in the central and northern parts, it is expected that there will be late onset of rainfall and higher peaks. In addition, the seasonal cycle of rainfall in the central part is

Conclusions
This paper examined the trend in annual rainfall and annual temperature of Mono river watershed over the observation period 1981-2010 and by 2050 using the regional model REMO under RCP4.5 and RCP8.5. It also assessed the monthly pattern of rainfall over the same periods. During the last three decades, rainfall and temperature have been increasing all over the Mono River watershed. By 2050 and under emission scenarios RCP4.5 and RCP8.5, annual rainfall is projected to be characterized by high variability, whereas a significant increasing trend is projected for annual temperature (warmer future climate). For each of the three defined regions (south, center and north), and under both emission scenarios, the seasonal cycle of rainfall is expected to change: In the southern part, the first peak is projected to reduce slightly, whereas the second peak is expected to increase and shift to October; however, in the central and northern parts, it is expected that there will be late onset of rainfall and higher peaks. In addition, the seasonal cycle of rainfall in the central part is expected to shift from a transitional regime to a unimodal one. Moreover, on a monthly scale, the northern part of the watershed is expected to record the highest increase and decrease in rainfall regardless of the emission scenario considered. These considerable changes in the monthly rainfall of the northern part are expected to occur globally in the dry season, thus indicating potential extreme events. Considering the projected trends and patterns for rainfall and temperature over Mono watershed by 2050, it is recommended that experts identify and implement relevant adaptation strategies.