Total Water Storage Change in Cameroon: Calculation, Variability and Link with Onset and Retreat Dates of the Rainy Season

Total water storage change (TWSC) was calculated using CRU (Climatic Research Unit) monthly gridded data for the period 1962–1993 over Cameroon. Investigations were conducted to link its annual cycle with both the beginning and the end of the rainy season. A method was derived as an alternative to determine onset and retreat dates of the rainy season. Two methods were used for the calculation of TWSC. The first method used potential evapotranspiration (PET) from the Thornthwaite formula (PET TH) and the second, CRU gridded PET data estimated from the Penman–Monteith formula (PET PM). A comparative study of the corresponding TWSC, namely TWSC TH and TWSC PM , respectively, was done. According to the preliminary results, the study area is classified as humid below latitude 8 • N and semiarid above. The results of the spatial and temporal variations showed a close correlation between the two methods, but with a slight gap between their different values, those of TWSC PM being larger and fluctuating less. The annual cycles of TWSC and PR generally showed similar patterns, and their intensities decreased from the southern part of the area (Equatorial forest zone) to the northern part (Sahelian zone). For mean TWSC = 0, two different points were identified: the first and the second corresponding dates matching the onset and retreat months of the rainy season, respectively, except in the arid area (Sahelian zone), where only the retreat month of the rainy season was perfectly determined. The delay observed in the determination of rainfall onset date in that area is assigned to PET formulas that are defined only for humid areas and to the influence of high temperature just before the beginning of the rainy season, promoting the rapid evaporation of soil water immediately after the first rains. Application of the same method (TWSC = 0) for the individual year showed similar performances. Although TWSC is always negative in Zone 3 and positive in Zones 1 and 2, the study of the interannual variabilities revealed an overall declining trend due to a stronger decrease in precipitation compared with PET. Moreover, the decrease during dry months is more remarkable than during wet months.


Introduction
In many African countries, ignorance of total water storage change is the main cause of the depletion of agricultural productivity and mismanagement of hydro-electric energy.Urbanism, civil engineering works and many other activities depending on hydrology also suffer severe damage.According to Indrani and Abir [1], any change in the hydrological cycle directly affects the availability and quality of fresh water, and the hydrological cycle is driven by the variety of complex processes and interactions closely linked to the atmospheric circulation and to temperature.Small changes in precipitation or evapotranspiration also imply changes in streamflows, water quality and ecosystem processes [2].Thus, as declared by Alexandris et al. [3], knowledge about the period of available soil/ground water and its quantity can determine optimum sowing and groundwater exploitation times.Given such interest, many scientists are committed to the study of water storage and related parameters like potential evapotranspiration, actual evaporation, total water storage, etc, and their environmental impacts under a changing climate [4][5][6][7][8][9][10].
Potential evapotranspiration is water likely to be lost in conditions where water is abundant in the soil [11].It is calculated using mathematical formulas [12][13][14][15].Among these formulas, that of Penman-Monteith [16] shows good results in several areas [13] and is recommended by the Food and Agriculture Organisation (FAO) [13,14].The disadvantage of this method is that it requires a significant number of climatic data, which is not always available in several regions of poor countries.In the contrary, the empirical method of Thornthwaite, which is still used, is one that needs only the temperature data.Actual evapotranspiration is the water actually dissipated into the atmosphere as vapor [11].It is difficult to measure accurately, but can be estimated through empirical mathematical equations [17,18] and Satellite estimates [19,20].According to Brutsaert [21], actual evapotranspiration and potential evapotranspiration have the same value under conditions of ample water availability at the surface.The total water storage change, which is defined in the simplest form as precipitation minus evapotranspiration, is difficult to measure due to the absence of observations of evapotranspiration [22].However, it can be estimated through satellite observations [5,[23][24][25], which were demonstrated to be efficient [26,27].However, despite this efficiency, it was found that the spatial concentration and distribution of water masses (reservoirs, lakes, glaciers and small-scale groundwater systems with spatial scales well below GRACE nominal resolution) affect GRACE (Gravity Recovery and Climate Experiment) water storage estimates [23].Because the installation and maintenance of evapotranspiration measurement networks are very expensive for poor countries, empirical formulas are often used to estimate the total water storage change.The most used include the atmospheric water balance [28,29], land water balance [17] and model simulations [19].Some authors [1] underline that the management of water resources is crucially important for seasonally arid countries, especially those whose rainfall is driven by monsoon because they are likely to experience dramatic changes in land use patterns.Therefore, the investigation on total water storage change over recent decades is necessary for understanding the large-scale water storage variation process.Using the approximation that inter-basin groundwater flow is limited, Cao et al. [17] modified the basin-scale water balance equation [28,29] to evaluate the evapotranspiration variations and the feedback of groundwater storage during the past 50 years in the North of China Plain.
In the scientific literature in general, studies have not been carried out on the link between total water storage change and the dates of the beginning and end of rains.The objectives of this work are to study the variabilities of total water storage change and to investigate its link with the onset and retreat dates of the rainy season.We focus our study on the Cameroon territory, one of the Sub-Saharan countries that has experienced problems of water scarcity during recent decades, and its northern part (arid-area), moreover affected by frequent extreme climate events.It is hoped that the results of this work will contribute to the understanding of the large-scale water storage variation process and provide guidance to farmers, hydrologists, agronomists, civil engineers and researchers.The next section (Section 2) describes the study area, and Section 3 presents the data used and defines the methodology.The results are shown in Section 4, and the work ends with the discussion and conclusion in Section 5.

Study Area
The study area chosen is in Central Africa and comprises Cameroon and some neighboring areas delimited by the coordinates 0 • -14 • N and 5 • -20 • E (Figure 1).This area was classified into three climatic zones according to Guenang and Mkankam [30]: the Equatorial forest zone (Zone 1), the midland zone (Zone 2), which predominantly covers highlands, and the Sahelian zone (Zone 3).In each of these zones, all stations have roughly the same onset and retreat dates of the rainy season.Rainfall onset and retreat dates in the area follow the northward (from Zone 1 to Zone 3) and the southward (from Zone 3 to Zone 1) direction, respectively [30].These movements are connected to those of the inter-tropical discontinuity (ITD) [30].The climate in Cameroon varies greatly with terrain, ranging from tropical rain along the coastal region through the hot semiarid northern plains to the arid Sahelian region in the far North.Most of the high intensity rainfall events in Cameroon are from convective clouds.Annually, there are several cases of extreme precipitation events.Major consequences are heavy rains, floods and landslides.Agriculture is the main economic activity (generally done at subsistence levels), and electricity is mostly produced by hydro-electric dams.The productivity of these activities is strongly dependent on the available soil moisture and water level of the rivers, respectively.

Data Used
The data used in this study are monthly precipitation (PR), surface temperature (TAS), potential evapotranspiration (PET), mean temperature (T mean ) and minimum and maximum temperatures (T min and T max ) covering the time period 1962-1993.They are from CRU (Climatic Research Unit) grid data Version 3.0 [14] with a 0.5 • × 0.5 • longitude/latitude spatial resolution.All of these data are downloadable free of charge from the CRU website (http://badc.nerc.ac.uk/data/cru/).CRU precipitation data have already been applied in the area.They were used to compute the multi-scalar Standardized Precipitation Index (SPI), and the results were compared to those from observation data [31].The results showed that CRU precipitation quite well corroborates the observations.

Methodology
Calculations are made for the three different zones of the study area.Beforehand, the aridity index (I a ) adopted by the United Nations [32,33] is calculated to characterize each part of the area.It is defined as the ratio of the climatological annual mean values of precipitation (PR) and potential evapotranspiration (PET) [6,32,34,35]: Depending on the values of this index, the areas are classified as follows [32,33]: hyper-arid for I a < 0.05, arid for 0.05 < I a < 0.2, semi-arid for 0.2 < I a < 0.5, dry sub-humid for 0.5 < I a < 0.65 and humid for I a > 0.65.
Under the hydrological assumption that runoff in the basin has been reduced to an insignificant level [36], Cao [17] proposes the following balanced equation to define the total water storage change ∆S (including soil water storage change and groundwater storage change): where ET is the evapotranspiration.
In the current study, the absence of evapotranspiration data in the study area forces us to use potential evapotranspiration instead of actual evapotranspiration.This substitution is rational given that the largest part of the study area satisfies conditions so that the two parameters match.This assumption is discussed in the manuscript.Therefore, we define total water storage change (TWSC) as: Monthly surface temperatures were used to compute potential evapotranspiration (PET) according to the Thornthwaite equation (PET TH ) [37]: where: PET TH is the estimated potential evapotranspiration (mm•month −1 ); T a is the average daily temperature (degrees Celsius; if this is negative, use 0) of the month being calculated; N is the number of days in the month being calculated; L is the average day length (hours) of the month being calculated.
is a heat index which depends on the 12 monthly mean temperatures T ai .For a comparative study, monthly CRU PET calculated from Penman-Monteith equation (PET PM ) [3,13,16] was used.PET PM or ET 0 is defined by the following formula: ET 0 is also called the reference evapotranspiration (mm•day For more details, see Allen et al. [16]. Removing PET TH and PET PM from PR leads to defining two parameters, namely TWSC PM (PR − PET PM ) and TWSC TH (PR − PET TH ), respectively.In a given area, the computation of each parameter was done for each grid point, per month and per year.Results were spatially averaged giving one data time series.The mean annual cycles of all parameters were computed, and an investigation was undertaken on the link between TWSC annual cycle and the onset and retreat dates of the rainy season computed by Guenang and Mkankam [30] using the cumulative percentage mean rainfall amount method [38].The study of interannual variabilities was also part of this work.
For each year, the onset and retreat dates were determined as those corresponding to the first and the second intercept, respectively, between the curve of TWSC and the reference line for TWSC = 0. Figure 2 illustrates this interpolation where m I is the unknown month.Considering the property of Thales in the approximative triangles II 1 A 1 and II 2 A 2 , the following relation can be written: This leads to the following equation: Finally, the date is deduced as followed:

Classification of Parts of the Study Area According to the Aridity Index
The following table (Table 1) shows the results of the calculation of the aridity index for the three different zones of the study area.From these results, Zones 1 and 2 below 8 • N are classified as humid, while Zone 3 above this latitude is classified as semiarid.

Annual Cycle of PR
The annual cycle of PR for all three zones of the study area is shown in Figure 3; the onset and retreat months of the rainy season are shaded.
Overall, PR shows three different regimes for the three zones.In Zone 1 (Figure 3a), PR shows two modes, the first one during the period MAMJ (March-April-May-June) (a peak observed in May) and the second one during ASON (August-September-October-November) (a peak observed in October).
These two modes are alternated with a dry season DJF (December-January-February) (the shortest of the study area) where the most severe period is from December-January and an intermediate season (June-July-August) recording rainfall ranging between 100 and 140 mm•month −1 .We note that during the year, the minimum PR is 50 mm•month −1 (Figure 3a).Furthermore, it is shown that the soil of this zone is mainly yellow ferralitic, characterized by a green forest vegetation and low water holding capacity, which is however not an important limitation since it rains almost all year long, and many rivers cross the area [39,40].Thus, the soil is always wet and promotes evapotranspiration.In Zone 2 (Figure 3b), PR only shows one mode ranging from April-October.It is the longest mode of the area (seven months long).Furthermore, the highest rainfall is recorded during August and September.The soil is slightly wet during the dry period NDJFM (November-December-January-February-March) compared with Zone 1.
The abundance of PR in Zones 1 and 2 is due to the fact that they are influenced by the West African monsoon activity and are bordered in their western parts by the Atlantic coast.Zone 3, a part of the Sahelian area is the driest.It has the shortest rainy season May-June-July-August-September (with a peak observed during August) followed by a dry period, the longest of the area during which the period of DJF does not record any rain.This lack of PR is an unfavorable condition for the application of PET formulas that are established on the basis of a soil still provided with water.For each zone, there are two gray vertical dashed lines representing the mean onset and retreat month of the rainy season [30].These months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).PET are from the Thornthwaite and Penman-Monteith methods.The labels on the x-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.

Annual Cycle of PET
Annual cycles of PET TH and PET PM are shown in Figure 3. Their patterns are similar in each zone, but with a gap of about 20 mm•month −1 , the values of PET TH being the largest.Overall, PET increases from Zone 1 (southern part of the area) to Zone 3 (northern part of the area).

Annual Cycle of TWSC and Its Link with Onset and Retreat Months of the Rainy Season
Figure 4 shows the mean monthly annual cycle of TWSC for all zones of the study area.TWSC TH and TWSC PM show similar patterns with the best correlation coefficients in Zones 1 and 2 (r = 0.993007).Overall, each pattern shows two points where TWSC is zero.These two points correspond to those of Figure 3 where there is an intersection between the PR curves and those of PET.The corresponding months of these two points match in many cases the months of beginning and end of the rains.The following table (Table 2) shows the lag between the months when TWSC is zero and those of the onset and retreat months of the rainy season.The estimates are for TWSC TH and TWSC PM .and (c) for the Sahelian zone named Zone 3 (north of study area).For each zone, there are two gray vertical dashed lines representing the mean onset and retreat month of the rainy season [30].These months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).The labels on the x-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.r indicates the Spearman rank order correlation coefficient between TWSC TH and TWSC PM .
Table 2. Errors on the estimated onset and retreat month of the rainy season using the TWSC method.

Zone Error on Onset Month Error on Retreat Month
Using TWSC TH Using TWSC PM Using TWSC TH Using TWSC PM It is observed that the two methods give similar results, but with a slight difference as a consequence of the difference between the values of PET TH and PET PM .For TWSC = 0 and comparing the corresponding first and second months to those of the onset and retreat months of the rainy season, respectively, we accurately obtain the same results for the onset month of the rainy season in Zone 2 and the retreat month in Zone 3.For the others, the onset and retreat months of the rainy season are obtained with a difference less than a unit of time in data (one month), except in Zone 3, where the difference between the two onset months is greater.The first cause of these errors is attributed to PET formulas that are limited to areas where the soil should be wet enough.These formulas become less efficient in arid areas and during dry periods, as is the case in the Sahelian zone (Zone 3) and especially during the period NDJFMA (November-December-January-February-March-April) when the deficit in soil water (TWSC < 0) is high, hence the delay observed in the onset determination of the rainy season.These formulas put things right when the rains begin to achieve a perfect determination of the retreat month of the rainy season in the arid area (Figure 4c).Second, the effect of temperature is to be taken into account in these errors.Indeed, Figure 5 showing the annual cycles of average maximum and minimum temperatures in the three zones highlights the highest temperatures in the Sahelian zone (Figure 5c).For each zone, there are two gray vertical dashed lines representing the mean onset and retreat month of the rainy season [30].These months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).The labels on the x-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.
Analysis of Figures 3 and 5 shows a strong correlation between PET and extreme temperatures (especially maximum temperature).They follow the same patterns.Consequently, the rise of air temperature contributes mostly to soil heat and therefore facilitates soil water evaporation.

Favorable Periods for Rainwater-Dependent Activities
The above analysis shows that months with TWSC equal to zero (PR − PET = 0) are very close to the onset and retreat months of the rainy season.When TWSC is greater than zero (PR − PET > 0), it is the rainy season, and the amount of rainwater is important enough to sustain evaporation and provide a water reserve to the area.This represents favorable conditions for plants to grow and for activities depending on the availability of soil/ground water storage.For TWSC less than zero (PR − PET < 0), rainwater is not sufficient to assure potential evapotranspiration, then the ground water reserve is requested to contribute for evaporation and will consequently decrease during the corresponding period (dry season).
Table 3 recapitulates the means of TWSC for each zone and for the two methods used.The differences between results from the two methods (TWSC PM − TWSC TH ) are added.
These results reveal that the two methods agree in determining humid periods except in Zone 2, where TWSC PM indicates eight months long against six for TWSC TH .In fact, the longest period with TWSC > 0 is observed in Zone 1, while Zone 3 shows the shortest one.These periods extend from March-November (nine months) in the Equatorial forest zone (Zone 1), from March-October (eight months) and from May-October (six months) according to TWSC PM and TWSC TH , respectively, in Zone 2 (transition zone), and from July-September (three months) in the Sahelian zone (Zone 3).Standard deviations on TWSC in most cases are higher in arid areas (Zone 3) and during dry months, while those from TWSC TH are most often slightly larger.This indicates stronger inconsistency of the Thornthwaite formula in determining PET in arid areas comparatively to the Penman-Monteith formula.The main cause is the unsatisfied potential conditions in arid areas and during dry periods in such a way that PET formulas become less efficient.For TWSC PM > 0, the difference between TWSC PM and TWSC TH (Table 4) shows a gap ranging from 19.58-32.55mm•month −1 , 10.64-20.61mm•month −1 and 10.73-17.05mm•month −1 in Zones 1, 2 and 3, respectively.In Zone 1, TWSC PM is greater than TWSC TH along the year.In Zones 2 and 3, it is greater during the rainy season or when it approaches.

Onset and Retreat Dates of the Rainy Season for an Individual Year Applying the TWSC Method
It is important to check how the method of TWSC behaves for a given year so that it can be used for projection, yearly planning activities or decision making.Figure 6 shows time series of the onset and retreat dates of the rainy season in each of the three zones applying the TWSC method.Onset and retreat dates of the rainy season for an individual year and using the TWSC method.TWSC TH and TWSC PM were considered for the three defined zones: (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and (c) for the Sahelian zone named Zone 3 (north of study area).The labels on the y-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.For each zone, there are two horizontal dashed lines representing the mean onset and retreat month of the rainy season [30] and labeled by one and two stars, respectively.These reference months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).
Results show that the dates obtained from TWSC TH = 0 and TWSC PM = 0 corroborate in most years, except very few cases where the discordance is height (retreat date in Zone 1 in 1964Zone 1 in , 1981Zone 1 in , 1983Zone 1 in and 1990).In general, the calculated dates are consistent with the reference months and are slightly earlier for onset and later for retreat, except with onset in Zone 3 (Sahelian zone), where the dates are the latest.The reasons for these weakness were mentioned in Section 4.2.3.

Interannual Variability of TWSC
Figure 7 shows the annual cycle of box-plots carrying some statistical information on the values of TWSC (TWSC TH and TWSC PM ) for the different zones of the study area.These boxes provide information about the TWSC mean, standard deviation and extreme values.The box-plots representing TWSC TH and TWSC PM have similar patterns, both following the shape of annual precipitation cycles, except in Zone 3, where TWSC shows an unexpected convex shape during the dry period.This deformation highlights the inconsistency of PET formulas for arid areas and especially during dry periods.From these figures, it is also noted that deficit in precipitation increases (TWSC decreasing) northward.The standard deviation between the interannual values of TWSC is generally greater during the rainy season in view of the importance of the height of box-plots during this season.The intersection with the reference line (TWSC = 0) is observed during the onset months of the rainy season in Zones 1 and 2 and during the retreat month in Zone 3. The other intersection points with the reference line match the onset and retreat months if the extreme values of TWSC are considered.In any case, mean TWSC TH is closer to the reference line than mean TWSC PM , except for onset date in Zones 2 and 3, where TWSC PM outperforms.
Figure 8 shows the time series of TWSC TH and TWSC PM for each zone of the study domain.The goal is to examine the effects of climate change on the interannual values of TWSC.Once again, we note that patterns of both TWSC TH and TWSC PM corroborate and that TWSC PM values remain somewhat larger over the years.The trends are all negative, with amplitudes increasing from Zone 1 to Zone 3, showing a decrease in TWSC with the magnitudes increasing northward, TWSC PM underestimating the decline in Zone 3 (Figure 8c).With regard to the positive values of TWSC in Zones 1 and 2, and negative in Zone 3 for most years, one is tempted to conclude that Zones 1 and 2 have not witnessed deficits in soil water (although declining) and that Zone 3 has almost always been in deficit; worse still, these deficits continue to decline.
In order to understand the reason why TWSC, although decreasing, was always positive in Zones 1 and 2 and negative in Zone 3, interannual PR, PET TH and PET PM were studied (Figure 9).The following statements were deduced from this figure: (a) precipitation amounts (PR) were always higher than PET TH and PET PM in Zones 1 and 2 and lower in Zone 3; (b) In the three zones, PR decreased with slopes in Zones 1 and 2 (−0.43 and −0.52, respectively) stronger than in Zone 3 (−0.37).In the meantime, trends in PET TH and PET PM , although increasing in most cases, remained too weak.Consequently, the difference PR − PET = TWSC decreased in all zones and remains positive in Zones 1 and 2 and negative in Zone 3. A careful examination of interannual TWSC at the monthly scale (for each month of the year) (Figure 10) reveals more details.
Zones 1 and 2 have indeed actually witnessed some deficits (TWSC < 0) in the months of January, February, March (less frequent in Zone 1), April and November (Zone 2 only) and December.These deficits that have been of low frequency and amplitude are masked by the high frequency and amplitude of the interannual accumulation of water in the soil during the other months.The opposite situation is observed in Zone 3 where overall soil water deficit outweighs soil water accumulation.The months when TWSC is positive are July, August and September, corresponding to the rainy season.

Discussion and Conclusions
We have calculated total water storage (TWSC) in three defined climatic zones over Cameroon by subtracting potential evapotranspiration from precipitation.Data used were CRU precipitation (PR), CRU PET from the Penman-Monteith formula (PET PM ) and PET TH that we calculated from the Thornthwaite formula using CRU precipitation and temperatures.Therefore, two parameters defined as TWSC PM (PR − PET PM ) and TWSC TH (PR − PET TH ) were derived and analyzed based on annual cycle and interannual variability for each climatic zone of the study area.
The aridity index helped us to determine that the area is humid below 8 • N of latitude and semiarid above.This classification led to better appreciation of the correlation that exists between PET calculated from mathematical formulas and actual evapotranspiration.Indeed, both values are equal in humid areas, while PET becomes greater than actual evapotranspiration when tending to arid areas.
It was found that TWSC and PR both decreased northward and showed similar patterns.The values of TWSC PM remained slightly higher and more consistent.This consistency of the Penman-Monteith PET formula confirms its efficiency compared to other formulas [13].By setting TWSC equal to zero, a close link was observed between the corresponding months and onset and retreat dates of the rainy season calculated by Guenang and Mkankam [30] using the cumulative percentage mean rainfall amount method [38], except in the Sahelian area (Zone 3), where a lag of about 1.5 months is observed in rainfall onset month determination.This delay is due firstly to PET formulas, which are defined only for humid areas.Therefore, they become less efficient in arid areas (Sahelian zone, e.g., Zone 3) and during the dry period when the soil is in deficit of water (TWSC < 0).These formulas put things right when the rains begin to achieve a perfect determination of the retreat date of the rainy season.All in all, the success of the method in Zones 1 and 2 comes from the fact that the condition of sufficient water availability at the surface is satisfied in such a way that the PET is equal to actual evapotranspiration [21], whereas that condition is not satisfied in semiarid area (Zone 3), especially during the dry season.Secondly, it is due to strong evaporation in arid areas (e.g., Zone 3) as the consequence of high temperatures during dry periods, especially at the approach of the rainy season.In fact, this high difference between the two dates occurred only in the hottest part of the study area (Zone 3), where monthly temperatures were generally higher and oscillated between 15 • C and 38 • C.During the period of the lag, intense soil heat facilitates strong evaporation of soil water after the first rains, letting the soil dry.Furthermore, it was observed that the interannual variabilities of rainfall onset in this zone are higher [30].
The study of the interannual variabilities revealed an overall declining trend of TWSC change as the result of the sharp interannual decrease (with a strong contribution of the decline during dry months) of yearly precipitation compared to PET.In many areas of the world, negative trends were also found in TWSC change estimated from different methods.For instance, Moiwo et al. [25] used satellite-based (GRACE and GLDAS) and ground-measured (precipitation) hydrological data to show water storage depletion in the North China Region [25].A study using a simple empirical climate model (ECM) and the Gravity Recovery and Climate Experiment (GRACE) satellite data found an overall negative storage change in the Himalayas and Tibetan Plateau [19]; they attributed the possible causes of this storage loss to human activities (e.g., farm irrigation) and climate change (e.g., temperature rise and global warming).Agriculture is expected to suffer from this change given that the availability of water resources is closely linked to precipitation and therefore to the rainy season.
This study brings more useful information for water resources management in the sense that it specifies the periods when the water can be available in the soil/ground.Taking these periods into account can avoid seed losses and may be of benefit to farmers.This can also help to plan activities that are rain-fed or hydrologically dependent, like agriculture, civil engineering works (bridge and road construction, etc.) or optimization of the hydro-electric energy production.The study also found that setting TWSC to zero and determining the corresponding dates can be used as an alternative method for the onset and retreat dates' determination of the rainy season.However, the use of potential evapotranspiration instead of actual evapotranspiration to calculate TWSC leads to a biased onset date in arid areas.In future investigations, it would be interesting to focus on the correction of this bias and to go further into the establishment of the link between TWSC and the onset and retreat dates of the rainy season; for instance, the use of data at finer scale (weekly means, pentad means or daily means) in order to increase the accuracy of the connection between these parameters.It would also be of great interest to use actual evapotranspiration if available.

Figure 1 .
Figure 1.Study area with CRU grid points.The topography is shaded, and the area is divided into three climatic zones as classified by Guenang and Mkankam [30].

Figure 2 .
Figure 2. Illustration of the interpolation method used to deduce onset and retreat dates of the rainy season applying the TWSC method.

Figure 3 .
Figure 3. Annual cycles of mean monthly PR and PET (expressed in mm•month −1 ) for the period 1962-1993: (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and (c) for the Sahelian zone named Zone 3 (north of study area).For each zone, there are two gray vertical dashed lines representing the mean onset and retreat month of the rainy season[30].These months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).PET are from the Thornthwaite and Penman-Monteith methods.The labels on the x-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.

Figure 4 .
Figure 4. Annual cycles of mean monthly TWSC TH and TWSC PM (expressed in mm•month −1 ) for the period 1962-1993: (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2;and (c) for the Sahelian zone named Zone 3 (north of study area).For each zone, there are two gray vertical dashed lines representing the mean onset and retreat month of the rainy season[30].These months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).The labels on the x-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.r indicates the Spearman rank order correlation coefficient between TWSC TH and TWSC PM .

Figure 5 .
Figure 5. Annual cycles of average, minimum and maximum temperatures (T mean , T min and T max respectively, expressed in • C) for the period 1962-1993: (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and (c) for the Sahelian zone named Zone 3 (north of study area).For each zone, there are two gray vertical dashed lines representing the mean onset and retreat month of the rainy season[30].These months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).The labels on the x-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.

Figure 6 .
Figure 6.Onset and retreat dates of the rainy season for an individual year and using the TWSC method.TWSC TH and TWSC PM were considered for the three defined zones: (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and (c) for the Sahelian zone named Zone 3 (north of study area).The labels on the y-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.For each zone, there are two horizontal dashed lines representing the mean onset and retreat month of the rainy season[30] and labeled by one and two stars, respectively.These reference months are absolutes and correspond to those containing onset and retreat dates in terms of pentads and their variabilities (not more than six pentads for the highest range).

Figure 7 .
Figure 7. Box-plots showing annual cycles of TWSC TH and TWSC PM (expressed in mm•month −1 ) for the period 1962-1993: (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and (c) for the Sahelian zone named Zone 3 (north of study area).The labels on the x-axis (J, F, ..., D) indicate the months of the year (January, February, ..., December) respectively.The months marked by one and two stars indicate the onset and retreat months of the rainy season, respectively.The parameters used to build a box are: mean TWSC (mid-value of box), mean + standard deviation (top value of box), mean − standard deviation (bottom value of box) and minimum and maximum TWSC (bottom and top values, respectively).

Figure 8 .
Figure 8.Time series of mean monthly TWSC (expressed in mm•month −1 ), calculated using PET obtained from the Thornthwaite and Penman-Monteith methods: (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and (c) for the Sahelian zone named Zone 3 (north of study area).The regression lines are shaded and trend coefficients indicated by b1 and b2 for TWSC TH and TWSC PM , respectively.

Figure 9 .
Figure 9.Time series of mean monthly PR, PET TH and PET PM per zone (expressed in mm•month −1 ): (a) for the Equatorial forest zone named Zone 1 (south of study area); (b) for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and (c) for the Sahelian zone named Zone 3 (north of study area).The regression lines are shaded and trend coefficients indicated by b0, b1 and b2 for PR, PET TH and PET PM , respectively.

Figure 10 .
Figure10.Time series of mean TWSC (expressed in mm•month −1 ) for each month (a-l) of the year and for each zone, calculated using PET obtained from the Thornthwaite and Penman-Monteith methods: for the Equatorial forest zone named Zone 1 (south of study area); for the transition zone between the Equatorial forest zone and the Sahel zone named Zone 2; and for the Sahelian zone named Zone 3 (north of study area).The trend coefficients are indicated by b1 and b2 for TWSC TH and TWSC PM , respectively.

Table 1 .
Aridity index calculated for each zone of the area and from PET using the Thornthwaite and Penman-Monteith mathematical formulas.

Table 3 .
Total water storage change (TWSC) from PET calculated using the Thornthwaite and Penman-Monteith formulas.Positive and negative values indicate the availability and deficit of rainwater in the area, respectively.Standard deviations are added to each value.

Table 4 .
Difference between TWSC PM and TWSC TH .