The assessment of observed and simulated flows over wet and dry sub-periods using the GR4J model and the Nash-Sutcliffe criteria is crucial for understanding the model's performance and reliability. Figure 3 shows the variability of experiential and simulated flows, allowing for a detailed comparison during different sub-periods, such as calibration and validation. By analyzing the wet and dry sub-periods separately, researchers can gain insights into the model's ability to capture the hydrological response under different hydroclimatic conditions. This analysis helps identify any biases or limitations in the model's representation of flow dynamics during wet or dry periods. It provides valuable information for water resource management, as the performance of the model can differ under varying climatic conditions. The use of various Nash-Sutcliffe criteria enhances the evaluation of model performance. These criteria, which measure the contract among experiential and simulated flows, provide a quantitative calculation of the model's accuracy. By applying dissimilar Nash-Sutcliffe criteria, such as the daily, monthly, or annual criteria, researchers can assess the model's performance at different temporal scales. This analysis helps identify potential discrepancies between observed and simulated flows and guides the refinement of the model calibration. The comparison of observed and simulated flows over wet and dry sub-periods provides valuable insights into the GR4J model's performance under different hydrological conditions. It allows for a comprehensive evaluation of the model's ability to capture flow variability and provides a basis for further model improvement. This information is vital for effective water management and decision-making, as it helps understand the model's reliability in simulating flows during wet and dry periods.
4.1. Variability of observed and simulated flows over wet and dry sub-periods in calibration and validation with the GR4J model using the various Nash-Sutcliffe criteria
Figures 3, 4, and 5 provide important insights into the analysis of this basin using GR4J model. Figure 3 depicts the variability of precipitation, which is a crucial input for hydrological modeling. Understanding the temporal and spatial patterns of rainfall is vital for accurately simulating flow in the catchment. The figure allows researchers and water managers to assess the distribution, incidence, and strength of precipitation events, which can have a important effect on the hydrological response in the basin. Figure 4 showcases a comparison between observed and simulated flows. This evaluation is essential for assessing the performance of the GR4J model in replicating the hydrological behavior of the Gambia basin. By visually comparing the observed and simulated flows, researchers can identify any discrepancies or biases in the model outputs. This analysis helps to understand how well the model captures the temporal variations and magnitudes of the actual flows in the basin. The figure provides a comprehensive assessment of the model's capability to reproduce the observed hydrological processes and informs decisions related to water management and planning. Figure 5 presents the probabilities of non-exceedance derived from the GR4J model using the calculated Nash-Sutcliffe criteria. This figure offers valuable information about the reliability and confidence of the model's predictions. It provides a probability distribution of flow values, indicating the likelihood of different flow levels occurring in the future. These probabilities are useful for decision-making processes, as they assist in understanding the range of possible flow scenarios under different conditions. Water managers can utilize this information to evaluate potential risks and develop appropriate strategies for sustainable water resource management. In summary, Figs. 3, 4, and 5 play a crucial role in the assessment and analysis of the Gambia basin using the GR4J model. They provide insights into precipitation variability, model performance in simulating flows, and the probabilities of non-exceedance. These visual representations enhance our understanding of the hydrological processes in the basin, guide decision-making, and contribute to effective water resource management.
According to Fig. 6, the GR4J model exhibited better performance during the validation phases compared to the calibration period. The correlation coefficient (r) values for the validation periods of 1991–2000 and 2001–2010 were 0.523 and 0.542, respectively, while the coefficient for the calibration period was 0.516. It is worth exploring whether enhancing the GR4J model's ability to match a proportional fraction of soil moisture could further improve its performance. However, this assumption requires further investigation in future research. In the GR4J model, the soil moisture fraction is calculated as the difference between available soil moisture and field capacity. In nature, soil moisture reaches saturation levels within 2 to 4 days, after which it gradually approaches field capacity through the process of soil water drainage. The GR4J model does not explicitly require an upper limit for saturation soil moisture, which may contribute to its effectiveness in simulating flow.
Table 1 shwos the outcomes of the model performance criteria during the calibration and validation phases. The criteria used include Nash-Sutcliffe efficiency (Nash), bias in whole flow volume (Q), and efficiency for each objective function applied to the gauging station to assess model performance. The research investgaton of this study is presented in Table 1. The scatterplots of observed and simulated daily hydrographs are shown in Figs. 1, 2, and 3, it is marked that the Nash-Sutcliffe criteria, including Nash-Sutcliffe (Q), Nash-Sutcliffe (VQ), and Nash-Sutcliffe (ln(Q)), yielded acceptable results. A acceptable relationship was observed among the observed and simulated flow data, with the Nash-Sutcliffe (Q) criterion placing significant importance on the differences between simulated and observed flood flows. The values obtained were 0.623, 0.711, and 0.578 during the calibration period (1981–1990) and the validation periods (1991–2000 and 2001–2010), respectively. The Nash-Sutcliffe (ln(Q)) criterion, which effectively represents changes in the hydrological system during low-flow periods and provides better performance for low flows, yielded values of 0.694, 0.711, and 0.737 in the calibration period and the validation periods, respectively. To assign similar weight to the simulation of flood and low-water flows (or the simulation of average flows), the Nash-Sutcliffe (VQ) criterion was utilized, resulting in values of 0.778, 0.827, and 0.800 during the calibration period and the validation periods, respectively during the calibration period (1981–1990) and the validation periods (1991–2000 and 2001–2010).
Table 1
Model performance criteria in the calibration and validation phase
Calibration period (1981–1990)
|
Efficiency criteria
|
Nash-Sutcliffe (Q)
|
0,623
|
Nash-Sutcliffe (VQ)
|
0,778
|
Nash-Sutcliffe (ln(Q))
|
0,694
|
Balance sheet
|
1,314
|
Validation period 1 (1991–2000)
|
Efficiency criteria
|
Nash-Sutcliffe (Q)
|
0,711
|
Nash-Sutcliffe (VQ)
|
0,827
|
Nash-Sutcliffe (ln(Q))
|
0,711
|
Balance sheet
|
0,997
|
Validation period 2 (2001–2010)
|
Efficiency criteria
|
Nash-Sutcliffe (Q)
|
0,578
|
Nash-Sutcliffe (VQ)
|
0,800
|
Nash-Sutcliffe (ln(Q))
|
0,737
|
Balance sheet
|
0,926
|
Calculation of the errors using the observed flow/simulated flow ratio gives an uncertainty of around 1.314 for the calibration period (1981–1990), indicating a slight underestimation of flows by the model compared with observed flows. As for the validation periods, these errors are 0.997 and 0.926 respectively in 1991–2000 and 2001–2010, indicating a slight overestimation of flows by the model compared with observed flows. We note here that the values of r and efficiency are close for the three periods. Overall, flows at the Simenti gauging station are both slightly underestimated and slightly overestimated by the model for all objective functions, depending on the weights used. Compared with all the objective functions, Nash-Sutcliffe (Q) shows the poorest performance in terms of volume bias for the different calibration and validation periods
Comparisons of correlation coefficients, volume bias, and efficiency results indicate that the GR4J model provided a better estimation of simulated flow when the Nash-Sutcliffe (VQ) objective function was used. Similar findings were observed in numerous calibration processes conducted in various catchments in Africa, where the Nash-Sutcliffe (VQ) objective function yielded improved estimates of daily flows, timing, and volume ratios (Vaze et al., 2011). It is noteworthy that the model's performance showed a shift from underestimation during calibration to overestimation during validation at the Simenti station. Interestingly, in absolute terms, the model performed better during validation than during calibration when using the Nash-Sutcliffe (VQ) objective function. Figure 3 depicts a curve indicating that, in certain years, the model failed to capture flows greater than 15 mm during both the calibration and validation periods at the Simenti station. Consequently, larger water volume deficits and even greater increases in water volume were observed during the validation period at the Simenti station. These disparities in volume could be attributed to variations in flow between the calibration and validation periods. Specifically, the mean annual flow recorded in the calibration period was 0.385 mm, while it increased to 0.603 mm in validation period 1 and further to 0.756 mm in validation period 2 (Table 2). As for the average flow simulated by the model, it is 0.142 mm in the calibration period (i.e. 0.142 mm behind the observed flow), 0.626 mm in validation period 1 (i.e. 0.023 mm above the observed flow) and 0.536 mm in validation period 2 (i.e. 0.220 mm behind the observed flow).
Table 2
Average flows observed and simulated (mm) by the GR4J model in the Casamance basin at Simenti
Period
|
Calibration (1981–1990)
|
Validation 1 (1991–2000)
|
Validation 2 (2001–2010)
|
Flow
observed
|
Simulated flow
|
Difference
|
Observed flow rate
|
Simulated flow
|
Difference
|
Observed flow rate
|
Simulated flow
|
Difference
|
|
Average flow
|
0,527
|
0,385
|
0,142
|
0,603
|
0,626
|
-0,023
|
0,756
|
0,536
|
0,220
|
|
High water
|
0,856
|
1,064
|
-0,208
|
1,338
|
1,216
|
0,122
|
1,558
|
1,088
|
0,470
|
|
Low water
|
0,046
|
0,139
|
-0,094
|
0,112
|
0,161
|
-0,049
|
0,178
|
0,137
|
0,041
|
|
Similarly, the mean high-water discharge observed in the calibration period is 0.856 mm, whereas it increases to 1.338 mm in validation period 1 and 1.558 mm in validation period 2. As for the model-simulated mean high-water discharge, it is 1.064 mm in the calibration period (i.e. 0.208 mm greater than the observed discharge), 1.216 mm in validation period 1 (i.e. 0.122 mm less than the observed discharge) and 1.088 mm in validation period 2 (i.e. 0.470 mm less than the observed discharge). For the low-water period, the simulated flow is greater than the observed flow in calibration (with a value of 0.094 mm) and in validation 1 (0.049 mm), while in validation 2, a delay is noted with a value of 0.041 mm.
Table 3
Correlation between observed and simulated flows, average performance and uncertainties in calibration and validation in the Casamance basin at Simenti
Station
|
Timing
|
Validation 1
|
Validation 2
|
Correlation coefficient
|
0,516
|
0,523
|
0,542
|
Average performance (%)
|
69,8
|
75,0
|
70,5
|
Average uncertainties in mm
|
0,142
|
-0,023
|
0,220
|
The correlation between observed flows and simulated flows (Table 3) shows better correlation coefficients for validation than for calibration at Simenti. However, the results obtained show that the GR4J model is an efficient model for simulating daily flows, especially in the Gambia basin at the Simenti station. The simulated values are close to those observed and attest to the model's validity. The analysis of the different performances (Table 3) shows first of all that the validation performances are better than the calibration performances. These performances are 69.8% for the calibration period (1981–1990), 75% for validation period 1 (1991–2000) and 70.5% for validation period 2 (2001–2010) at the Simenti station. Uncertainties are low and are of the order of 0.142 mm over the calibration period (i.e. underestimation of flows by the model), -0.023 over validation period 1 (i.e. overestimation of flows by the model) and 0.220 over validation period 2 (i.e. underestimation of flows by the model). Furthermore, it is important to note the differences in rainfall-runoff ratios during the calibration and validation periods in the basin. The average rainfall for calibration period 1, calibration period 2, and validation period is 1158 mm, 1203 mm, and 1289 mm, respectively. These statistics indicate varying relationships between rainfall and runoff in the basin during these periods. Additionally, there is a notable increase in both rainfall and runoff observed during the validation period, particularly during the high-water period at the Simenti station. Considering these characteristics, it is possible that the parameters calibrated to capture the phenomena during the calibration period might not accurately reproduce the flow during the validation period. One factor that influences this is the parameter associated with groundwater exchange, which affects conveyance. During the calibration period, the precipitation exceeds the flow at the Simenti station. The presence of both negative and positive parameter values, as well as their magnitudes, can be attributed to the differences in the time periods between calibration and validation. This emphasizes the importance of having similar climatic conditions between the calibration and validation periods, as parameters derived from a wetter (or drier) period may not be applicable to a drier (or wetter) period.
4.2 Future hydrological trends
The assessment of climate change impacts on hydrological systems plays a vital role in the development of adaptive strategies for water resource management, flood risk control, mitigation measures, and ecological protection. These assessments typically involve the integration of future projections of climatic factors, such as precipitation and temperature, derived from climate models with hydrological models to simulate and understand the potential changes in hydrological processes under future climate conditions. By combining these models, we can gain insights into how climate change may affect water availability, flood patterns, and overall ecosystem health, enabling us to plan and implement appropriate measures to adapt to these changes effectively.
4.2.1. Future annual hydrological trends
To characterise the annual variability of future flow in the Gambia basin at the Simenti station, characteristic annual flow values are shown in Table 4 and Fig. 7. As observed data for the different hydrological components of the basin were available, they were used as well as the model outputs as baseline or reference data, allowing a comparison of the flow components for the future scenarios with the simulated historical values. At the Simenti station, where the interannual modulus over the reference period (1985–2014) is 168 m3/s for observed flows and 187 m3/s for simulated flows, the future period (2021–2100) recorded a modulus of 98.3 m3/s under the SSP 245 scenario, i.e. an average deviation of -41.5% compared with the observed flow over the reference period, and a modulus of 91.4 m3/s under the SSP 585 scenario, i.e. an average deviation of -45.5% compared with the observed flow over the reference period.
Table 4
Future changes in flow (in %) on an annual scale over the four future periods at the Simenti station
SSP245
|
1985–2014
|
2021–2040
|
Trend
|
2041–2060
|
Trend
|
2061–2080
|
Trend
|
2081–2100
|
Trend
|
Flow rates
|
168
|
108
|
-35,8
|
94,0
|
-44,0
|
98,3
|
-41,4
|
93,2
|
-44,4
|
Precipitation
|
1281
|
1134
|
-147
|
1107
|
-174
|
1109
|
-172
|
1113
|
-168
|
Average temperature
|
28,0
|
29,0
|
0,99
|
29,6
|
1,65
|
30,3
|
2,27
|
30,7
|
2,71
|
FTE
|
1467
|
1511
|
44,4
|
1540
|
73,8
|
1569
|
102
|
1588
|
121
|
SSP585
|
1985–2014
|
2021–2040
|
Trend
|
2041–2060
|
Trend
|
2061–2080
|
Trend
|
2081–2100
|
Trend
|
Flow rates
|
168
|
119
|
-29,1
|
108
|
-35,7
|
74
|
-55,6
|
64
|
-61,6
|
Precipitation
|
1281
|
1141
|
-140
|
1128
|
-153
|
1058
|
-223
|
1031
|
-250
|
Average temperature
|
28,0
|
29,1
|
1,08
|
30,2
|
2,26
|
31,6
|
3,66
|
33,2
|
5,22
|
FTE
|
1467
|
1515
|
48,8
|
1568
|
101
|
1631
|
164
|
1700
|
234
|
The variations in rainfall and temperature that controlled to variations in possible flow are presented in Table 4. The analysis of runoff over the future period shows a decrease in runoff over the different sub-periods under both scenarios, in phase with the decrease in rainfall. There is a constant downward trend in the predictable hydrological variations for entirely scenarios. The decrease in discharge over the upcoming period diverse considerably between the four scenarios. Ranked from smallest to largest over the horizons, this decrease in flow, compared with the flow observed over the reference period, is of the order of -35.8%, -44.0%, -41.4% and − 44.4% respectively for the horizons 2040, 2060, 2080 and 2100 under the SSP 245 scenario. For the SSP 585 scenario, this reduction over these four horizons is respectively − 29.1%, -35.7%, -55.6% and − 61.6%. While for the first two horizons (2040 and 2060), the decline is greater under the SSP 245 scenario, for the last two horizons (2080 and 2100), the decline is greater under the SSP 585 scenario.
4.2.1. Future monthly hydrological trends
The flow characteristics on a monthly scale in basin area at the Simenti station over the situation period and the upcoming period are shown in Table 5 and Figs. 8 and 9. The river regime at this station is therefore characterised by a three-month high-water period and a nine-month low-water period. The regime of the Gambia catchment at Simenti is characterised by a maximum in September (683 m3 /s over the reference period, 401 m3 /s under the SSP 245 scenario and 360 m3 /s under the SSP 585 scenario) and a minimum in May (0.66 m3 /s over the reference period, 6.66 m3 /s under the SSP 245 scenario and 6.52 m3 /s under the SSP 585 scenario): this is a unimodal regime. The high-water period lasts three months, making it a pure tropical regime.
Table 5
Average monthly flow values for the reference and future periods at the Simenti station
|
SSP 245
|
SSP 585
|
|
1985–2014
|
2021–2100
|
Difference (%)
|
1985–2014
|
2021–2100
|
Difference (%)
|
January
|
20,7
|
25,4
|
22,6
|
20,7
|
24,4
|
17,6
|
February
|
8,94
|
16,1
|
80,3
|
8,94
|
15,5
|
73,9
|
March
|
3,00
|
11,4
|
278
|
3,00
|
11,0
|
267
|
April
|
1,63
|
8,48
|
421
|
1,63
|
8,27
|
408
|
May
|
0,66
|
6,66
|
911
|
0,66
|
6,52
|
889
|
June
|
3,81
|
5,51
|
44,5
|
3,81
|
5,41
|
42,1
|
July
|
73,1
|
8,56
|
-88,3
|
73,1
|
8,63
|
-88,2
|
August
|
352
|
104
|
-70,3
|
352
|
87,1
|
-75,2
|
September
|
683
|
401
|
-41,3
|
683
|
360
|
-47,3
|
October
|
420
|
298
|
-29,0
|
420
|
295
|
-29,9
|
November
|
120
|
116
|
-3,52
|
120
|
111
|
-7,62
|
December
|
48,8
|
47,9
|
-1,77
|
48,8
|
45,5
|
-6,68
|
Rainy season
|
306
|
163
|
-46,6
|
306
|
151
|
-50,7
|
Dry season
|
29,1
|
33,1
|
13,7
|
29,1
|
31,7
|
9,00
|
Year
|
168
|
98,3
|
-41,4
|
167,7
|
91,4
|
-45,5
|
The average monthly rates of change in flow below the two dissimilar SSP scenarios for the dissimilar times are presented in Figs. 8 and 9. Overall, there is a descending trend in flow for the months July to December ahead, but the decrease varies considerably from month to month. In general, the greatest decrease is noted over the far upcoming under both scenarios and occurs during the high-water period. On a monthly scale, the decrease in flow is expected to be greatest in August and September in both the near and distant future. Moreover, the maximum value of river discharge has changed from September in the historical period and the near future (2021–2040) and medium future 1 (2041–2060) to October in the medium future 2 (2061–2080) and far future (2081–2100) under the SSP 585 scenario (Fig. 8).
The range of deviations for the change in projected future flow is shown in Fig. 8. A disparate evolution is associated with the change in flow, with a large difference between high-water months with negative deviations (decrease in future flow) and low-water months with positive deviations (increase in future flow). The rainy seasons (high-water period) will record a fall of around − 46.6% under SSP 245 and − 50.7% under SSP 585. Conversely, dry seasons (low-water periods) will increase by 13.7% under SSP 245 and 9.0% under SSP 585.
An analysis of annual flow changes in the Gambia basin during the historical period reveals an uneven distribution of monthly flow at the Simenti station due to the seasonal influence of low flow originating from the water table. This uneven distribution has implications for the hydrological and ecological processes within the basin. The anticipated decrease in flow, coupled with the high uncertainty of future floods, will significantly impact the hydrological resource system and weaken its overall functioning. Understanding the changes in the contribution of different hydrological components to climate change within the river basin is of great importance. This study aimed to simulate changes in the main hydrological components relative to a baseline period (1985–2014) using a model to represent future scenario periods. By analyzing and comparing the simulated flow components, insights can be gained into the potential changes in the hydrological system under different climate change scenarios. Figure 9 shows that in the near, medium and distant future, discharge decreases under SSP 245 and SSP 585, respectively, compared with the control period (1985–2014). The changes and the magnitude of the changes in discharge are very different are consistent in the different scenarios with rainfall. The Gambia basin, located in the West African region, exhibits hydroclimatic variability that aligns with the projections of the Intergovernmental Panel on Climate Change (IPCC). Projections from climate models indicate a decrease in average monthly flows at the Simenti station under two scenarios (SSP 245 and SSP 585) for future time horizons. These findings suggest that surface water resources in the catchment areas of the basin are expected to continue declining throughout the 21st century, particularly under the more severe SSP 585 scenario. The vulnerability of the Gambia basin to climate change is evident, given the projected greater reduction in rainfall. The observed trends in water resources align with similar findings in other regions, highlighting the extreme scenarios presented by CMIP6. Runoff in the Gambia basin is primarily influenced by rainfall, and its seasonal changes correspond to those of precipitation patterns projected by climate models. The results obtained were based on the use of the GR4J global hydrological model, which is relatively simple and relies on four parameters. Utilizing a distributed hydrological model could potentially enhance the accuracy and reliability of these results.