Urban Gully Development in Southwestern Nigeria

This study examined the development of urban gullies in a part of Southwestern Nigeria. This was done with a view to determine the volume and amount of sediment loss from the 1 and 2 order gullies in a cascade system. Twenty (20) towns were randomly selected from where thirty (30) gully systems comprising forty (40) 1 order and five (5) 2 order gullies were measured with Tape and Total Station (TS). The gully dimensions measured include gully lengths, depths and cross sectional areas. The core samples were also taken from the study gullies for the determination of the soil bulk density. The results showed that the mean volume of soil loss from the 1 order and (2) order gullies in the study area are 1612.633 m and 1629.922 m, respectively. The weight of the estimated sediment loss from the 1 order and 2 order gullies are 2661.621 and 2657.318 tonnes, respectively. The result further showed that the total soil loss of 119661.442 tonnes was obtained from the study gullies and that the 1 order gullies accounted for 87% (106,375 tonnes) of the total soil loss. This was an indication of soil degradation in the area.


Introduction
Gullies have been recognized as an important environmental threat in many parts of the world [1].This phenomenon has been identified as a major factor in soil degradation, water quality deterioration and changes in channel morphology in the humid tropics [2].Gullying could become even more severe in areas of high population growth especially in an urban centre and high rainfall intensity [3].A number of gullies developed because of a decrease in the erosional resistance of the land surface or an increase in the erosional forces acting on land [4].
Soil erosion as related to anthropogenic factors in Nigeria is well documented in agriculture [5].But, however, the fact that several Nigerian towns and cities are still experiencing erosion is not much known/documented.Though, urban gulling affects a relatively small part of the total land area of Nigeria, the land value of the portion of the land and infrastructure destroyed make it by far more harmful than in agricultural land [6].Soil erosion in urban setting in Nigeria, has not been adequately studied or documented.For instance, [7][8][9] and [10], among others have focused mainly on erosion processes and the nature of the earth materials removed within some urban centres in Southern Nigeria, while [11,12] examined off site effects of erosion on urban stream channels.Also, past studies, for example, [13,14] in AdoEkiti; [15] and [16,17] in Ilorin, focused on the contribution of rainfall parameters on sediment generation and distribution from various land use surfaces in the urban areas.These studies showed that increase in volume of runoff generated from large areas of impermeable urban surfaces coupled with loose and incoherent nature of the soil on steep surfaces and from unpaired drains contributed largely to gully formation/development in urban areas [see 10].Much as these studies contributed to the existing knowledge on urban gullying, however, they are based on very few gullies in very few urban centres.Urban gullying is becoming a topical issue worldwide and the outcry for prompt and proper management of this phenomenon call for in-depth knowledge and under-standing of gully initiation and development.Hence, this study attempts to assess the development of gullies from twenty towns/ cities in Ondo and Ekiti States of southwestern Nigeria.The main objective of this study is to document sediment loss from 1 st and 2 nd order gullies in urban areas of southwestern Nigeria.constitute the study area.The Ondo State is partially on sedimentary formations while Ekiti is mainly on Basement complex rocks.The choice of the two states was as a result of their true representation of southwest Nigeria (excluding Lagos State) in terms of geology and physiographical characteristics.The study area which is part of the eastern section of Southwest Nigeria lies between latitude 5˚57'N and 9˚12'N and longitude 2˚40'E and 6˚03'E (see Figure 1).
The population of the study area is 5,825,236 (Ondo-3,441,024 and Ekiti-2,384,212 (NPC, 2006).The people of the region are predominantly Yorubas comprising sub-groups such as Akures, Ekiti, Ondos, Ikales and Ilajes.Culturally, before the 19 th century, the Yorubas settled in agglomerations or clusters either in towns or villages.Prominent among these towns/cities are Akure, Owo, Ondo, Ikare, Okitipupa, Ado-Ekiti, Ikere, Ijero and Efon Alaaye among others.As observed by [7], the development of most of these towns and villages occurred without any systematic planning.Even today, the poor The study area is underlain especially at the southern part by Post-Jurassic sedimentary rocks and alluvium known as Coastal plain sands while the northern and central parts are mainly underlain by quartzite, quartzschist, magnetite, gneisses, granite gneisses, charnockites and coarse porphyritic-biotite granite rocks of the Precambrian Basement Complex [18].Specifically, Oke-Imesi, Ijero, Efon-Alaaye, Ilawe axis of the study area was underlined by major rocks of the Ife-Ilesa Schist belt which comprises of Amphibolites complex schist's and Quartzite sequence [19].These rocks constitute the parent materials of the Oke-Messi soil Association which comprised of different soil series.The most important are the gravelly, sandy and hill drift series [9].On the whole, 70% of the studied gully systems fall within the various aspect of Basement Complex while 30% are within Neocene in the Sedimentary formations.
The area is characterized by humid tropical climate (Am of Koppen's climatic classification) with mean annual temperature of 27˚C and a mean annual rainfall ranging between 1250 mm and 1400 mm and distributed between the months of March and October with peaks in July and September, and short dry spell in August, though, this varying from year to year in its occurrence [20].The rainfall effectiveness is between 6 -9 months in the year [21].The onset and withdrawal of the rains are marked by thunderstorms accompanied by high rainfall intensity [22].
Also, the forest and derived savanna communities of the study area are derivatives of the tropical rain and tropical deciduous forests which developed in a climatic region characterized by the dominance of humid over arid tropical conditions [23].Gallery forests are also common along stream banks.However, this has also been severely altered as a result of encroachment from urban land use [24].

Methods and Materials
All settlements that fall within the demographic standard of the United Nation definition of urban towns of 20,000 and above inhabitants were compiled on state basis using 2006 population census figures on town basis.However, some towns that failed the above definition were included if they served as administrative headquarters or as university towns.The towns were further classified on the basis of those with gullies and without gullies.In all, 25 towns (13 in Ekiti and 12 in Ondo States) were identified to have serious gully occurrences.From the list, 20 towns were randomly selected using table of random numbers.This was done in order to have reasonable representation of the gullies in the study area.From each selected town, at least one gully site was selected for detail measurement of gully dimensions using line tape and Global, GPS Positioning system.In all, thirty (30) gully systems involving forty (40) 1 st order and five (5) 2 nd order gullies were selected and studied.
The processed data from the GPS and Total Station (TS) were stored in the Note pad of Microsoft office.The values of X and Y coordinates were later exported into Auto CAD environment from which gully widths, lengths and surface areas were determined using the geostatistical tool of the software.The average of the gully widths in the gully channel represents the gully width and was used in the determination of the cross sectional areas of the study gullies.
The cross sectional area for each point along the gully length was calculated using the equation described by [25,26] and later expressed as where, A = Cross sectional area (m 2 ); w = width of the gully (m); d = depth (mean) depth of the gully (m).
The values of the cross sectional areas obtained were used to determine the volume of soil loss from the study gullies.Two methods available for the determination of volume of soil loss are End Area and Prismodal formula [see 27,28].The End Area formula was used for the calculation of volume of soil loss because of its ability to calculate volume irrespective of whether the number of cross sectional areas is odd or even.The equation for the End Area is expressed as: where, d = the distance between cross sectional areas y 1 = first sectional area; y 2 = second sectional area.Soil samples were taken from the gully shoulders (tops) and floors.The number of samples from each gully depends on the length and slope gradient [29].Based on this principle, 60 samples were collected for textural properties analysis.In addition, core samples were also taken using McCauley sampler for the determination of soil bulk density which was used to estimate the weight of sediment loss (tonnes) from the study gullies.The soil particle size analysis was determined using hydrometer method [30].The details on the methods adopted in this study are documented elsewhere [31].

Soil Physical Characteristics
The soil physical characteristics of the study gullies are shown in Table 1.As evident from the table, the values of bulk density obtained at the top soil, ranged between 1.29 g/cm 3 at Ijero and 1.82 g/cm 3 at Oke Ibedo, Ilawe with a mean of 1.53 ± 0.112 g/cm 3 while the gully floor has values which ranged from 1.49 g/cm 3 at Odo Owa to 1.96 g/cm 3 at Usiri, Ikole and a mean of 1.72 ± 0.135 g/cm 3 .As shown in the table, all the soils sampled irespective of depth exhibited high bulk density values which are far above 1.0 -1.3 g/cm 3 considered for well aggregated forest or grassland vegetation [32,33].Also, the values obtained in this study was in the same range for the values obtained for urbanized city of Miami-Dade and Gainesville in Florida which were averaged at 1.63 g/cm 3 and 1.52 g/cm 3 , respectively [34,35].The implication of high values of bulk density in the study area is that there would be reduction in the infiltration capacities of the soils, thus making the soil prone to greater runoff and soil loss.When this happens, more water would be generated as surface runoff in the case of a steep slope terrain and flooding as in the case of lowland area.When the run off flows downslope, there is tendency for it to create rills which may later converge down slope, especially if it is concentrated, to form gullies.The above process explains gully formation at Idogun and Orunbemiku in Ode Irele, Adu in Akure and Odo Ese 1 st order gullies II and III in Oke Imesi where the greater proportion of the surface runoff flows over heavily compacted urban roads from upslope of the study catchments.
Further, the soils in the study area are predominantly sandy clay with mean % sand (53.43%), while clay and silt are 32.45% and 14.12%, respectively at the gully top.At the gully floor, the clay content exhibits highest mean percent (44.88%), followed by sand (42.64%) and silt (12.48%)Table 1 also showed that sand content decreases with soil depth while clay content increases with soil depth.Ordinarily, the high proportion of sand at the upper top soil could not have supported gully initiation but other factors such as high volume of run off as a result of increase in impermeable surfaces and other anthropogenic factors, most of the urban bare soils are easily susceptible to action of running water (runoff).The soils in this study exhibited higher clay content and higher bulk density at the gully floor.This further confirmed the finding [35] in their study of urban gullies in Sao Luis City, Brazil, where it was observed that soil with higher proportion of sand content associated with lower bulk density while soil with higher clay content have higher bulk density.At the gully top, the values of clay ratio ranged between 1.00 and 3.60 at Oke-Ibedo, Ilawe and Ajeka, Okitipupa, respectively with a mean value of 2.30.The values at the gully floor was between 0.49 at Ikalako, Aramoko and 4.06 at Odo Osun, Igbara Oke with a mean value of 1.48.The mean value of clay ratio (3.30) as an index of soil erobidibility further revealed the susceptibility of the top soil to gully initiation by running water in the study area.

Gully Morphometry and Development
Tables 2 and 3 showed the morphometric properties, volume of soil loss and total sediment loss from the study gullies.As evident from Table 2, the gully lengths ranged between 15m at Ijemoji, Igede and 485.38 m at Ode Ose in Owo with the mean length of 131.44 ± 90.70 m.The gully depths ranged from about 0.6 m to 7.0 m with the shallowest depth at Odo pupa in Igede and the deepest at Ita Osun in Efon Alaaye.It has a mean value of 1.66 ± 1.20 m.Also, gully head depth varied from 0.13 m (1 st order) at Odo Ese to 6.70 m at Ilutitun with the mean value of 1.76 m ± 1.33 m.Fifteen of the 1 st order gullies had head depths ≥ 2 m, fourteen of which were formed as a result of abrupt termination half way of concrete channels.Gully widths ranged from 1.25 m at Idogun in Ode Irele to 17.09 m at Ode Ose in Owo with the mean value of 5.12 ± 3.76 m. based on the order, the mean depth values of 5.075 ± 3.93 m and 5.506 ± 2.284 m are obtained for the 1 st and 2 nd order, respectively (see Table 3).All these values of gullies dimensions suggest that gullies in the study area are moderate [36] when compared with observations by [25] in a part of southeastern Nigeria.Also, the deep gullies in the area exhibited unstable side wall as a result of sloughing and basal undermining.Sloughing is usually proceeded by side failure.Examples of gullies formed in this way are those at Okitipupa, Ilutitun and Ita Osun in Efon Alaaye.
The slope gradients of the study gully catchments ranged from 1˚10' to 8˚25', ordinarily under vegetal cover, these slope gradients would not enhance erosion processes but due to exposure to direct raindrop impact and human activities coupled with poor civil engineering work to control the runoff, gully has become pronounced virtually in all the catchments.For instance, building of drains along urban streets and concrete channels is a deliberate attempt to control surface runoff but the improper handling of such projects often accelerated or triggered gully formation and growth at the receiving location of the runoff.Specifically, many of the channels and drains constructed are at times abandoned or terminated half way before getting to the natural drainage channels (water bodies) (see Plates 1(a) and (b), Plates 2(a) and (b)).A total number of 24 (60%) 1 st order gullies were identified to be triggered or accelerated in the study area due to abrupt termination of concrete channels half way.Gullies in this category include Ibaka in Akungba, Ikese in Oka-Akoko, Osinle in Akure, ItaOsun in Efon Alaaye etc. Also, as evident from   systems in the entire study area.This value was in the same range with the sediment loss of 116,402 tonnes obtained from a single urban gully system at Queen Ede College Road, Benin City in 1992 (see [9]).Out of this total sediment loss, 106,375 tonnes (89%) was estimated from the 1 st order study gullies.However, it is of interest to note that sediment loss from Efon Alaaye gullies (10697.80tonnes) was drastically reduced when compared with 31,150 tonnes obtained for the same gullies in 1995 [9].This decline was as a result of aggressive engineering erosion control measures adopted by the government.

Conclusions
This study investigated the development of urban gullies in southwestern Nigeria.Gullies from twenty towns randomly selected from both Ondo and Ekiti States were measured for determination of their dimensions (lengths, widths, depths and surface areas) as well as volume of soil loss and weight of sediment loss.It was established that gullies in the area were of 1 st and 2 nd order, with 1 st order dominating (85%) the gully systems.The gully dimensions measured revealed a mean width ranged from 1.25 m at Idogun to 17.09 m at Ose in Owo with overall mean depth of 5.12 ± 3.76 m.The  mean depth ranged approximately from 0.6 m to 7.0 m with shallowest gullies at Odo pupa, Igede and deepest at Ita Osun, Efon Alaaye with 1.66 ± 1.20 m.The gully head depth varies from 0.13 m at Odo Ese, Oke Imesi to 6.7 m at Eddy Hotel, Ilu titun with a mean value of 1.756 ± 1.33.It must be noted that the deep gullies exhibited unstable sidewall as a result of sloughing involving basal undermining and slumping.Also, the mean cross sectional area (CSA) ranged between 1.05 m 2 at a 1 st order gully, Idogun in Irele and 59 m 2 at Eddy Hotel in Ilutitun.The highest soil loss of 24,232 tonnes was obtained at Ose in owo.The second and third highest values of 17,545 and 11,546.142tonnes were obtained from Eddy Hotel, Ilutitun and Oke-Ibedo, Ilawe, respectively.Most of the 1 st order gullies development or formation may not be unconnected to uncontrolled street surface runoff.It was also revealed that 14 out of 15 gullies with the gully head depth ≥ 2 m were formed as a result of uncontrolled surface runoff.
On the whole, gully initiation and development in the study area (Ondo and Ekiti States) have contributed immensely to environmental degradation in many of the towns studied.Also, the values of gully dimensions (widths, depths and lengths) indicated that the development has not reached a disaster level, but the restoration of damage done require assistance from the outside especially from the government and Non-Government Organizations (N-GOs).

Figure 1 .
Figure 1.Map of Southwestern Nigeria showing the Study Area (Ondo and Ekiti States).

Plate 1 .
(a) & (b) Typical drains constructed to conduct surface runoff at a street in Ore.
Plate 2. (a) A Gully Caused by the Half Way Termination of a Concrete Channel at Ose Gully in Owo; (b) A Gully Caused by the Half Way Termination of a Concrete Channel at the Back of the Gospel Faith Mission Church, Oke Ila, Ado Ekiti.

Table 3 ,
a total of 119661.442onnes of sediment loss was obtained from all the gully t

Table 1 . Soil Physical Properties of Gully Sites in the Study Area.
Source: Fieldwork, 2011; C = Clay; SC = Sandy clay; SL = Sandy loam, SCL = Sandy clay loam; TOP = Soil properties at gully top; FLO = Soil properties at gully floor.

Table 2 . Morphometric Attributes of the Studied Gullies.
The gullies are 2 nd order and they have no gully heads.