Research Article - (2018) Volume 6, Issue 1
This study was about Identifying climate change, by adopting the downscaling techniques in Umuahia metropolis. The study adopted ex-post-facto research design and data for maximum and minimum temperatures were collected from the archive of NIMET and for the period 1986-2015. On the other hand large scale predictor’s data were collected from the archive of HadCM3 for these periods 1960-2001 (NCEP) and 1960-2099(HadGCM3). Analyses were done using SDSM, ANOVA, PPMC and MLR. The study unraveled that, Mean sea level pressure (MSLP), Relative humidity at 500 hpa (r500), Relative humidity at 850 hpa (r850), and Temperature at 2 meters above sea level (temp) are the predictors of minimum and maximum temperature in the area. This also showed statistical significance at P<0.05. During validation the monthly sub-model performed better by using these indices for minimum and maximum temperatures respectively R20.85 & 0.70; RMSE 2.14 & 2.72; Rs 0.88 & 0.56; Value 0.00 & 0.00. Conversely, minimum and maximum temperature showed temporal variation for the period 1960-2080 at P<0.05 (F, 284.1) & p<0.05 (F, 227.1) respectively, therefore indicating significant change in temperature characteristics. The study strongly advocate assembling a working group that will work on a regional downscaling project, forging a synergy between Nigerian meteorological agency and the working group, producing a localized GCM, and the need to carry out similar study across the other regions of Nigeria.
Keywords: Climate change; Signals; Down scaling; Models; Normal
Climate is expressed as weather elements’ characteristics observed and documented for a period of 30-35 years [1]. Therefore any change observed in a normal in comparison with the previous normal is termed climate change [2-9]. However, Ayoade [10,11], has argued that, to adequately affirm that, there is change in climate, the best practice is to compare results from three or more normal. This is because the change between two normal may just be ‘noise’. Thus to be able to explain climate change, the intergovernmental panel on climate change has always encouraged the building of General Circulation Models (GCMs) [12]. This is so because it has been verified [12-17] to be an adequate means of explaining climate change issues and forecasting [18] the likelihood of change in climate, and what effects climate change is going to have at regional and global scales [16]. However, the scales at which the GCMs are built are so large, thus the need to perform downscaling [19,20].
Downscaling has emerged as a means (in climatological studies) of resolving what is known at a larger scale to what is obtainable at the local or even point scale [5-8,9,20-22]. It is usual to carry out downscaling using either or both of the following techniques: dynamic and statistical downscaling [18] although this study utilized statistical downscaling. Nevertheless, whether or not the resolution of the various GCMs is improved on, there will still be need to downscale [21]. This is because of the following advantages that the downscaling technique possess: it links the large scale predictors with the predict and (s); identifies the major controls of weather parameters; reproduces the predictand(s) datum or data with high precision; identifies scenario based weather information of the area for which downscaling investigation is carried out [12,20].
All the same, in the developed world historical data of weather parameters are largely available [1] for climate explanation. While in the developing world the case is not the same [10] as station density is relatively sparse [10] incidences of wars, poor training of weather observer, period of weather station establishment etc. all of which affects weather data stock available. As such, identifying historical characteristics of weather for several places in the developing world is almost impossible. This is the case in Umuahia. Weather, data recording, started in the weather station in 1977, thus, trying to identify climate change from such data situation is difficult or spurious. Downscaling techniques therefore offers a scientific solution to this lacuna. This study therefore, sets out to unravel signals of climate change from downscaled data of Umuahia in Abia state.
Study area
Umuahia is capital of Abia State, South-eastern Nigeria, and has been growing in leaps and bounds since becoming state capital on August 27, 1991. The influx of civil and public servants, traders and contractors with the new status brought about rapid population growth and urbanization of the town. The area is a typical humid tropical area with its location at the intercept of latitude 050 29’N and longitude 7033’E, and altitude, at the highest point of about 205m (Figure 1). Annual rainfall totals vary from about 1900 mm to 2400 mm. This rainfall is of the double maxima type with peaks in the months of July and September. In recent times, September has being the month with highest amount of rainfall as against July whereas December has lowest amount (NRCRI, 2016). Air temperature ranges between 21°C and 34°C, and relative humidity of 60% to 83% (NRCRI, 2016).
In terms of methods the study adopted the ex-post-facto research design. Data for temperature parameters (maximum and minimum temperatures) were collected from the archive of the Nigerian meteorological agency (Abuja) for the period 1986-2015. On the other hand large scale predictor’s data were collected from the archive of HadCM39 (Table 1).
Predictor | SDSM code | Predictor | SDSM code |
---|---|---|---|
Mean sea level pressure | MSLP | Surface airflow | P_f |
Surface zonal airflow | P_U | Surface meridional velocity | p_v |
Surface vortices | P_z | Surface wind direction | p-th |
Surface divergence | P_zh | 500-hpa air flow strength | P5_f |
500-hpa zonal velocity | P5-U | 500-hpa meridional | P5_v |
500-hpa vortices | P5_Z | 500-hpa geopotential height | P500 |
500-hpa wind direction | P5th | 500-hpa divergence | P5zh |
850-hpa airflow strength | P8_f | 850-hpa zonal velocity | P8_U |
850 meridional velocity | P8_v | 850-hpa vortices | P8-z |
850 geopotential height | P850 | 850-hpa wind direction | P8th |
850-hpa divergence | P8zh | Relative humidity at 500-hpa | R500 |
Relative humidity at 850-hpa | r850 | Near surface relative humidity | rhum |
Surface specific humidity | Shum | Mean temperature at 2m | temp |
Table 1: Large scale predictor’s data were collected from the archive of HadCM39.
To identify the appropriate grid cell of the hadCM3 that corresponds to the area under investigation, the hadCM3 grid cell was super-imposed on the map of Niger-delta region and the area was picked out within the map. Niger delta region because, the HadCM3 resolution is too large (412 km2) for the area under review, as such trying to identify the grid cell corresponding to the area needs to incorporate adjoining areas to the study area (Abia). Therefore the grid 3,31 of the HadCM3 GCM was selected. Having selected the grid cell, data for large scale predictors were downloaded from http://www.cics.uvic.ca/scenarios/sdsm/select.cgi. However, analyses were done in two categories. This first was achieved via the downscaling techniques and in the environment of the SDSM and the second form of analysis was achieved via the analysis of variance (ANOVA) and Pearson’s product moment correlations (PPMC) in the environment of the SPSS.
Conversely, the SDSM software operates in these itemized steps: quality control (this ensures that data are properly arranged and that they are complete as calibrated), predictors screening (this identifies the various predictors that would likely predict the predict and based on the r value, partial r value, and the p value realized), model graduation or calibration (this step computes a MLR between the predictors and predict and also tells the proportion of influence in the predicted that can be explained by the predictors), weather generation (this regenerates the weather of the area by creating up to 20 ensembles or likelihood of occurrence in the predicted, graphing model out puts, statistical analysis, & scenario generation (here the predictors of the Hadcm3 is used to reproduce the weather element(s) of interest in area where the downscaling study is being carried out and in this case Umuahia metropolis). Generally, temperature is deemed continuous in their distribution, its uncertainty (minimum and maximum temperatures) were computed by comparing the monthly means of the observed data with those of the downscaled using the Wilcoxon signed rank test at 95% level of significance. This uncertainty analysis is necessitated because there was need to find out the ability of the calibrated model to type the real world scenario which also occurs in a dynamic manner. Having downscaled, the minimum and maximum temperatures of the area, ANOVA was then applied on the downscaled temperature to find out changes in the temperature characteristics using the epochs under consideration (1960-2080).
In Table 2, the predictors used for this study are presented and were found to be valid at p<0.05. The result also showed that predictors valid for minimum and maximum temperatures downscaling in the area of study are; Mean sea level pressure (MSLP), Relative humidity at 500 hpa (r500), Relative humidity at 850 hpa (r850), and Temperature at 2 meters above sea level (temp). The chosen predictors in Table 2 also exhibit some physical association with the predict ands (minimum and maximum temperatures) and have been identified by Wilby et al, as basic predictors of temperature characteristics in the tropics and by extension have also been used by Mahmoud and Babel [13], in similar studies with varying levels of success (Table 2).
Predictors selected | SDSM Code | Predict and | Correlation coefficient | Partial correlation |
---|---|---|---|---|
Mean sea level pressure | MSLP | T min & T max | 0.87 & 0.73 | 0.63 & 0.61 |
Relative humidity at 500-hpa | r500 | T max | 0.92 | 0.86 |
Relative humidity at 850-hpa | r850 | T max | 0.81 | 0.71 |
Temperature at 2 meters above sea level | Temp | T min & T max | 0.88 & 0.94 | 0.74 & 0.89 |
Table 2: Selected predictors for model calibration.
During the model calibration three performance indices were used and include P value, Pr value, R value; and two indices were used for checking of model performance, that is, explained variance (EV) and the standard error (SE). Therefore the mean explained variance analyzed ranged between 63-76% and 65 and 83% for minimum and maximum temperatures respectively, while standard error fell between 0.8 and 1.25°C for both minimum and maximum temperatures respectively. This result showed consistency with that of previous related works (Souvignet & Heinrich, 2011; Rashimd & Mukand, (2014). In terms of validation, daily temperature data for the period 1992-2001 was generated by the SDSM and matched with observed minimum and maximum temperature for same period. Performance indicators for the validation were R2 values; root means square error (RMSE) and ratio of the standard deviation (Rs) (Table 3).
Ncep/Hadcm3 | R2 T max & T min | RMSE T max & Tmin | Rs T max & T min | P- Value at 0.05 T max, T min |
---|---|---|---|---|
Daily | 0.27 & 0.21 | 4.35 & 4.52 | 1.00 & 1.02 | 0.00 0.00 |
Monthly | 0.85 & 0.70 | 2.14 & 2.72 | 0.88 & 0.56 | 0.00 0.00 |
Seasonal | 0.52 & 0.61 | 3.72 & 3.58 | 0.94 & 0.87 | 0.00 0.00 |
A2 experiment | ||||
Daily | 0.26 & 0.23 | 4.26 & 4.56 | 1.00 & 1.02 | 0.00 0.00 |
Monthly | 0.72 & 0.81 | 2.12 & 1.94 | 0.92 & 0.86 | 0.00 0.00 |
Seasonal | 0.51 & 0.59 | 3.74 & 3.59 | 0.96 & 0.98 | 0.00 0.00 |
Table 3: Validation of the SDSM data (1991-2001).
The indicators in (Table 3) were computed by synthesizing NCEP & A2 scenario predictor’s data with predict and (T min & T max). However, NCEP performance indicators derived, performed better than that of the A2 scenario experiment predictors, reason being that the model was calibrated from the NCEP predictors. Conversely the model performed better at the monthly sub model scale than at other periods of both seasonal and daily scales (Table 3). This finding is inconsistent with the finding of Babel and Mahmood (2014), whose models performed better at the daily sub model scale for temperature parameters investigated in their study. Adduced reason for this variability in model performance Levels include variation in the environments in which the studies were carried out (i.e. semi-arid & tropical environments) (Table 4).
Umuahia met station max temp | Umuahia met station min temp | |||
---|---|---|---|---|
Months of the year | NCEP | A2 | NCEP | A2 |
January | 0.543 | 0.143 | 0.111 | 0.343 |
February | 0.723 | 0.723 | 0.123 | 0.823 |
March | 0.543 | 0.243 | 0.058 | 0.443 |
April | 0.123 | 0.051 | 0.123 | 0.351 |
May | 0.343 | 0.803 | 0.343 | 0.903 |
June | 0.423 | 0.723 | 0.113 | 0.921 |
July | 0.543 | 0.543 | 0.543 | 0.513 |
August | 0.12 | 0.723 | 0.22 | 0.829 |
September | 0.443 | 0.543 | 0.143 | 0.663 |
October | 0.723 | 0.7 | 0.313 | 0.553 |
November | 0.523 | 0.053 | 0.123 | 0.663 |
December | 0.013 | 0.023 | 0.153 | 0.153 |
Table 4: Uncertainty analysis for temperature downscaling in the study area.
However, there exist some identified uncertainties in Table 4. In the table, NCEP predictors revealed some level of uncertainty for the month of December as the P value for that month is less than 0.05. The same uncertainties are traceable to the A2 scenario for T max for the months of November and December respectively. However for T min, NCEP showed that there are little or no uncertainties between the modeled data mean and the observed data means for all months except for the month of March, which experiences the beginning of rainfall in the area and probably affected the simulation. However, in the A2 experiment for T min, December showed some uncertainty between modeled and observed data at P<0.05 (Table 5).
Sum of Squares | Df | Mean Square | F | Sig. | |
---|---|---|---|---|---|
Between Groups | 299.86 | 3 | 199.953 | 227.138 | .000 |
Within Groups | 159100.2 | 43796 | 13.683 | ||
Total | 159400.1 | 43799 |
Table 5: Downscaled maximum temperature.
Furthermore, Table 5 explains the temporal variation in maximum temperature from downscale data for the area under investigation. It is evident from the table that there is a significant difference in temperature for the period 1960 and 2080. This finding is also valid at p<0.05 (F, 227.1), however, the Duncan statistics in Table 4 revealed that the difference exists in three different categories and the first category is 1961-1990 & 1991-2020 with mean maximum temperature of 31.6Cand 32.5C respectively, while the other two categories are 2021-2050 and 2051-2080 with mean maximum temperature of 33.7C and 34.2C respectively (Table 6).
max_Asnr | ||||
Duncana | N | Subset for alpha = 0.05 | ||
Identifier | ||||
1 | 2 | 3 | ||
1961-1990 | 10950 | 31.6414 | ||
1991-2020 | 10950 | 32.5075 | ||
2021-2050 | 10950 | 33.6456 | ||
2051-2080 | 10950 | 34.2057 | ||
Sig. | 0.317 | 1.000 |
Table 6: Duncan statistics for maximum temperature.
Indications are clear that temperature in the area has continued to rise since 1960 and will continue to rise even till 2080 (Figure 2), and by extension indicates climate change (Tables 7 and 8).
Sum of Squares | Df | Mean Square | F | Sig. | |
---|---|---|---|---|---|
Between Groups | 499.661 | 3 | 166.554 | 284.141 | 0 |
Within Groups | 28308.9 | 43796 | 11.655 | ||
Total | 28808.56 | 43799 |
Table 7: Downscaled minimum temperature for Abia State.
max_Asnr | |||||
Duncana | N | Subset for alpha = 0.05 | |||
Identifier | |||||
1 | 2 | 3 | 4 | ||
1961-1990 | 10950 | 23.0329 | |||
1991-2020 | 10950 | 23.8058 | |||
2021-2050 | 10950 | 24.6860 | |||
2051-2080 | 10950 | 24.9150 | |||
Sig. | 1.000 | 1.000 | 1.000 | 1.000 |
Table 8: Duncan statistics for minimum temperature.
Table 8 reveals that minimum temperature in the area is significantly different over time in the area. This statement is also valid at P<0.05 (F, 284.1). Nevertheless, the Duncan post-hoc test in Table 8 shows that the mean temperature of the four epochs under investigation using downscaled data are significantly different from each other at P<0.05. The indications are also clear that the areas minimum temperature has witnessed significant change of an upward trend (Figure 3). This also indicated that the temperature of the area has been and is still going to change if nothing is done to reduce GHGs emission in the area. Therefore adequate and well-structured action plans are required for reparations of the environment and by extension the meteorological characteristics of the study area in particular and Nigeria in general.
In conclusion, evidence from downscaled data of Umuahia in Abia state indicates that temperature (minimum and maximum) have changed in the area. Nevertheless, the identified change is also attributed to assumptions and factors of climate change as advocated by HadGCM3 (which include but not restricted to, increased and uncontrolled urbanization, GHGS emissions proliferation, and increase population). However, the temporal patterns of temperature indicated that there are changes from normal to normal for current and future temperature patterns although with some uncertainties in the forecasts. This is probably already causing some severe environmental impacts and the fear is that if necessary palliative steps are not taken the environment may witness irreparable change. Nonetheless, the predicted changes in temperature pattern are seen to be continuously on the increase (although with some uncertainties), if anthropogenic conditions are sustained in the area.
However, the researchers having noted that the study adopted only one GCM data set, advices that the assertions of this study be utilized with caution, this is one limitation of the current study. A suggestion to this is that there is need to utilize more GCMs in the area for a more robust study and also to expand the study beyond the boundaries of the current study to account for the boundary conditions which may have not been captured in the current study, since the intension was to capture only the temporal information of the temperature characteristics and to find out how well the HadCM3 captured the climate information of the area. Also worthy of note is the expansion of predict ands, so that a more all-encompassing statement can be made about the climate situation in the study area and beyond in Nigeria. Therefore there is need to do the following:
• Assemble a working group that will work on the suggestions made above
• A synergy be established between Nigerian meteorological agency and the working group
• There is need to produce a localized GCM, which for now is not available
• Also there is need to carry out similar study across the other regions of Nigeria
• The state legislature is by this study, implored to review critically the existing laws on GHGs emission and other anthropogenic activities that exacerbate GHGs emission, for to improve on them and create stiffer sanctions for erring members of the polity.