Modeled Implications of Global Warming Conditions on Precipitation Totals in Valley and Highland Areas of Nigeria ()
1. Introduction
The study on rainfall response to global warming in low and highland areas of Nigeria is significant because these marginal areas are important agricultural zones of the country; hence an understanding of the variation in rainfall behaviour is vital for agricultural development of the physiographic regions. Agricultural practices in Nigeria are predominantly rain-fed and rainfall characteristics determine the agricultural calendar ( Olayide et al., 2016 ). Changing rainfall totals as a response to global warming will affect agricultural activities and production. The valley areas are flood plains which are used for the cultivation of rice, sugar cane and vegetables while on the highlands crops like Irish potatoes and many varieties of fruits that thrive in cooler environments are cultivated. This is because the low lying valleys used for flood plain cultivation can be inundated reducing the usage of these marginal lands for agriculture. Similarly, erosion on slopes of the highland areas may be aggravated thus destroying farmlands or reducing productivity through fertility loss by erosion.
Many settlements are situated in the low lying valleys of rivers Niger and Benue. Response by increasing rainfall total will worsen the occurrences of flooding which is already prone in these settlements ( Ikusemoran et al., 2013 ). Also, other environmental disasters like erosion and landslides may be intensified especially in the highlands ( Aliba et al., 2015 ). Decreasing response of rains will lead to reduction in fishing activities and hydro-electric power generation which will adversely affect standard of living in the study area. Changing precipitation will worsen issues of migration and conflict in the region. Clashes between cattle herders and farmers have become pronounced due to competition for climate based resources and the situation will get worse if conditions become drier ( Chukwuma, 2020 ). Modeled projection of the implications of the global warming conditions on precipitation totals in the valley and highland areas will avail information on the anticipated future rainfall occurrence in the study area. This will help to plan ahead the expected vulnerability and climate risk, thus embarking on timely adaptation measures.
Topography can affect precipitation in at least three important ways. Firstly, it can intensify the solar heating over land by providing an elevated heat source. This magnifies the regional-scale land-sea temperature contrast and facilitates the onset and maintenance of monsoon regimes that produce convective rainfall and latent heating. Variations of this heating can affect weather and climate in other tropical areas, as well as in the extra tropics via atmospheric teleconnections. Secondly, topography can contribute to orographic rainfall and additional latent heating. That is, as moist air is forced up a slope, the air cools, moisture condenses and rainfall and latent heating occur with consequences similar to those associated with convective heating. Thirdly, topographic barriers in the tropics obstruct the large-scale air flow and force some air to go around rather than over topography. This mechanical effect disturbs the flow downstream and sets up large-scale circulation patterns that extend into the extra tropics via planetary long waves ( Meehl, 1992 ).
Topography also has specially marked effects on weather elements such as air humidity, temperature and rainfall. These effects result in concomitant changes in soil and vegetation distribution over hilly or mountainous areas ( Korner, 2004 : Viviroli & Weingartner, 2004 ). The mountain ecosystems help in stabilizing atmospheric circulation by creating barriers to free movement of winds.
The temperature trend in Nigeria since 1901 (Figure 1) shows increasing pattern, the mean air temperature in Nigeria between 1901 and 2005 was 26.6˚C while the temperature increase for the 105 years was 1.1˚C ( Odjugo, 2010 ). This is obviously higher than the global mean temperature increase of 0.74˚C recorded since 1860 when actual scientific temperature measurement started ( Spore, 2008 ; IPCC, 2007 ). Similarly, according to the Federal Ministry of Environment (2014) , temperature has been on the increase in Nigeria in the last five decades and has been very significant since 1980s. An increase in global precipitation is anticipated but this increase will be uneven across different regions, these calls for regional analysis of precipitation characteristics. Rainfall has however shown a general decreasing trend in Nigeria since 1901 (Figure 2) but with regional differences. This study therefore aims at modeling the response of precipitation to global warming conditions on different physiographic areas of lowland and highland in the Middle Belt Region of Nigeria. The specific objectives are to determine grand average of monthly rainfall and harmonics curves for the
selected stations in lowland and the highland from 1961 to 2016 and model the implications of the global warming conditions on rainfall occurrences in the lowland and highland areas.
This study is limited to the valley area of Makurdi and the highland of the Jos Plateau which are located in the Middle Belt Region of the country within latitudes 6˚24'N to 11˚42'N and longitude 2˚42'E to 13˚6'E covering a total area of 320,000 km2, about 35% of the total land area of Nigeria (Figure 3). The choice
of Makurdi and Jos is because the two stations are the ones with the required recorded data with highest difference in elevation; Makurdi (92 m) and Jos (1286 m) and both stations experience rainfall abnormality in the country by recording rainfall amount that does not coincide with their latitudinal values. The period of study (1961-2016) 56 years is chosen to cover the recent five decades with reported cases of intense global warming and extreme rainfall events worldwide. The base year of 1961 is decided because of data availability in both stations of Makurdi and Jos.
2. Research Methodology
To model the implications of the global warming conditions on rainfall totals in the valley and highland areas, the Harmonic analysis technique was used. Harmonic analysis of a time series uses a Fourier series to analyze periodic fluctuations. It is a particularly useful tool in studying annual precipitation patterns as it reveals the spatial variation of various precipitation characteristics ( Isikwue et al., 2013 ). Harmonic Analysis, which is commonly applied to study periodic variations, decomposes a time series into its constituent parts if the time series represents a periodic phenomenon. It transforms a complex time series to a sum of many sinusoidal functions or harmonics ( Park, 2008 ). The Harmonic analysis is based on a mathematical principle that a curve, viewed as a function, may be represented by a series of trigonometric functions. This is as given by Isikwue et al. (2013) :
(1)
where;
;
(2)
is the observed value at time, t,
is the arithmetic mean, number of observations, i is the number of harmonics and P is the period of observation. In other words, the time series equals the mean plus the sum of all N/2 harmonics. The Equation (1) above can be rewritten as;
(3)
The type of variation dominating the curve is revealed by a comparison of the sizes of the amplitudes Ci, where;
is the amplitude of the ith harmonic and
, arc sin (Ai/Ci) is the time at which the i harmonic has a maximum. It can also be expressed in percentage. A harmonic with overwhelming contribution would definitely account for most of the periodic variation in the data, while the contributions of the other harmonics would be considered negligible. Large first harmonic amplitude suggests strong annual variation, while comparatively large second harmonic amplitude points to strong semiannual variation.
Autocorrelation plots of the data were made in order to determine whether significant relationship in the data. Paired sample T-test will be used to determine whether there is a significant relationship between the estimated rainfall values of the years 1961 and 2016 and their corresponding observed values. This serves as a test for the efficiency of the harmonic analysis model in estimating the rainfall amount for 20 and 50 years ahead in Makurdi and Jos representing two different altitudes.
Harmonic analysis on the average monthly rainfall observations was performed for the 56 years under study using the Turbo Pascal for windows programming language in order to implement the computation. This was carried out by fitting a periodic function of sinusoidal character to enhance the determination of the contribution of each harmonic (expressed as a percentage of total variation in the rainfall measurements it accounts for), the amplitude of each harmonic and the time at which each harmonic is maximum. The rainfall data was partitioned into monthly totals and each year was considered on the basis of the twelve months that make up the year. The 12 months coding is shown in Table 1.
3. Results and Discussion
The harmonic plots of both the valley and highland stations of Makurdi and Jos respectively show that the observed multiple peaks of the total rainfall in the time series plot are now reduced to bimodal patterns, by the amplitudes of the oscillations of the troughs and valleys of the harmonics. Hence, the harmonic analysis could have filtered or smoothened out some noise or spikes in the data. The harmonics plot of average monthly rainfall (mm) for the interval of 12 months for Makurdi and Jos is given in Figure 4 and Figure 5 respectively.
The result of harmonic plots of the grand or net average monthly rainfall (mm) for the interval of 12 months for the valley area of Makurdi and highland area of Jos show that there exists a similarity in the trend of variations of harmonic series plots of the data. Both plots exhibit single rainfall maxima or peak as shown in Figure 6 and Figure 7 respectively. That is, rainfall patterns are mono-modal in nature. This indicates that in a particular year, the rainfall peaks are observed in one single month of August for the valley of Makurdi and either July or August for the highland area of Jos. This implies that rainfall peak is attained earlier in the highland areas of Jos than the valley areas of Makurdi. This result partially agrees with the findings of Ogungbenro & Morakinyo (2014) in their work on rainfall distribution and change detection across climatic zones in Nigeria which revealed that rainfall peaks in July for areas within the Guinea Savanna Zone; but the result refutes their conclusion on the occurrence of August break and a double rainfall maximum within the belt. This difference could be attributed to the unique physiographic attributes of the valley and highland and their influence on rainfall behaviour as opposed to the other regions within the Savanna vegetation belt.
Results extract from the run of harmonic analysis on average monthly rainfall (mm) for the interval of 12 months for Makurdi and Jos are summarized in Table 2 and Table 3 respectively. Table 2 shows that the first harmonic dominates
Table 1. Code for the month axis along with the year
Source: Authors’ fieldwork, 2018.
Source: Authors’ fieldwork, 2018.
Figure 4. Harmonic plots of average monthly rainfall for Makurdi.
Source: Authors’ fieldwork, 2018.
Figure 5. Harmonic plots for average monthly rainfall for Jos.
Source: Authors’ fieldwork, 2018.
Figure 6. Grand average of monthly rainfall and harmonics curves for Makurdi.
the periodic components in the monthly average rainfall of Makurdi as it has the highest percentage contribution of 92.27%, indicating that the monthly rainfall in the station actually fluctuates. The large amplitude of the first harmonic indicates strong annual variation of rainfall in the station as in the works of Wong et al. (2009) and Kirkyla & Hammed (1989) . Results from Table 3 shows that the
Source: Authors’ fieldwork, 2018.
Figure 7. Grand average of monthly rainfall and harmonics curves for Jos.
first harmonics also dominates the periodic components in the monthly average rainfall of Jos as it has the highest percentage contribution of 91.67, indicating fluctuation or periodicity of the monthly rainfall in the station. The periodic fluctuation is dominated by halve a cycle as is evident in the grand average monthly plots in Figure 6.
Using the grand average monthly rainfall obtained as 105.08 and 105.41 mm for Makurdi and Jos respectively as model estimates from the harmonic analysis, the period of 12 months, and the sine and cosine coefficients in Table 2 and Table 3, the periodic function Xt for the monthly average rainfall of Makurdi and Jos were obtained using extracts from equations
and
. Where Ai and Bi are coefficients
of sine and cosine respectively and i’s are integers ranging from 1 to 3 as given in Table 2 and Table 3. The equations were then implemented in the harmonic analysis program to make a short consecutive five year and long five years interval forecasts of monthly average rainfall measurements for Makurdi and Jos. The result of short consecutive five year and long five years interval forecasts of monthly average rainfall is presented in Table 4 and Table 5 for Makurdi and Table 6 and Table 7 for Jos respectively.
Modeled result of the short consecutive five years and long five years interval forecast of average monthly rainfall for the low land area of Makurdi is presented in Table 4 and Table 5 respectively. The result shows an increase in rainfall occurrences from the month of January to August. The forecast however, shows a decrease in rainfall from the month of September to December. The implication of the result is that there will be more rainfall earlier in the season in the lowland areas and lesser of the rainfall towards the end of the season. In other words, more rains should be expected immediately after onsets and lesser rains should be anticipated before cessations. The effect of this rainfall pattern on agriculture is that the second phase of crops that mature towards the end of
Table 2. Result extract from the run of harmonic analysis program on average monthly rainfall for Makurdi, Nigeria.
Source: Authors’ computation, 2018.
Table 3. Result extract from the run of harmonic analysis program on average monthly rainfallfor Jos, Nigeria.
Source: Authors’ computation, 2018.
Table 4. A display of five (5) years forecasts of average monthly rainfall for Makurdi.
Source: Authors’ computation, 2018.
the rainy season will experience moisture deficiency. Also, decrease in rainfall towards the end of the rainy season will reduce the recharge of underground aquifers and worsen the perennial problem dry season water scarcity in the area.
The occurrence of more rainfall earlier in the season in the low land area may
Table 5. A display of five (5) years interval forecasts of average monthly rainfall for Makurdi.
Source: Authors’ computation, 2018.
be explained by the fact that since air descent is more pronounced in the low land (a condition that inhibits rainfall formation), most of the rainfall in the low land region are formed outside of the region but drifted into the area by regional wind movement most especially the squall line movement earlier in the season during the northward movement of the Inter Tropical Discontinuity Zone (ITDZ).
The modeled result of the short consecutive five years and long five years interval forecast of average monthly rainfall for the high land area of Jos presented in Table 6 and Table 7 respectively shows a decrease in the occurrence of rainfall in the months of January and February (Out of season rains). Meanwhile, the forecast shows an increase in rainfall occurrences from the months of March to December. Since the rainy season normally begins in Jos area in April, the implication of the result is that the high land area increasingly becomes drier during the dry season and wetter during the wet season. Also, the increasing rainfall during the wet season will aggravate the ability of running water to detach, erode and transport soil down slope on the steep slopes of the plateau; this in the long run will impoverish the soils, destroy buildings and other structures.
A validation of the forecast results was done. Observations show that the equations exhibit a good fit to the average monthly rainfall as they produced very close estimates of the actual monthly average rainfall. They yielded the same means (average monthly rainfall) as that of the actual data; and very close standard deviations of the actual and that of the model estimates as shown in Table 8, which also display the results of the forecasts with the corresponding actual and estimates of monthly average rainfall measurements. These estimates are
Table 6. A display of five (5) years forecasts of average monthly rainfall for Jos.
Source: Authors’ computation, 2018.
Table 7. A display of five (5) years interval forecasts of average monthly rainfall for Jos.
Source: Authors’ computation, 2018.
known as model estimates. Table 8 displays both the measured and the estimated rainfall values in Makurdi and Jos respectively. The level of deviation measured from the estimated rainfall values were determined using Paired sample T-test analysis. The results show that there are no significant differences at
Table 8. Validation results for Makurdi and Jos stations.
Source: Authors’ computation, 2018.
95% confidence level between the measured and the estimated rainfall in both years as can be seen from the T and P values. Hence, the predictions for the other years ahead are valid as far as rainfall variations in these stations are concerned.
4. Conclusion and Recommendation
Climate change and global warming are topical issues that have come to occupy environmental studies worldwide. The effect of global warming impinges on all facets of life including precipitation. Precipitation is one of the most important elements of weather and climate and in the tropics where agriculture is predominantly rain fed, water supply is majorly from natural sources, an understanding of the direction of response of precipitation to warming global conditions becomes paramount. This study shows the direction of response of precipitation to global warming on the unique marginal physiographic regions of the Benue valley and Jos plateau. The study recommends that refined modeling techniques that bring high quality term forecasts should be developed and used in immediate regions as well as diffusing of information to help interpret forecasts in terms of their agronomic and economic implications to users in local environments.