Effect of Desertification on Some Selected Soil Properties in a Semi Arid Part of Northwestern Nigeria

Gudu, being a fringe local government area of Sokoto State of Nigeria, is in the heart of desertification zone of Nigeria where agriculture as the sole means of livelihood for the ever increasing population is endangered. The extent of the effect of desertification on the soil quality is unknown, and for sustainable livelihood, timely information on vegetation and soil resources of the area is of essence. It was therefore hypothesized that desertification has effect on the Soil Organic Carbon (SOC) and particle size. Objectives include assessment of vegetation condition of the area, vegetation change detection and identification of hotspots and assessment of soil properties in the hotspots. Biennial time series (1998 and 2014) of 1-km 10-day SPOT Normalized Difference Vegetation Index (NDVI) was subjected to Time Series Analysis (TSA). Vegetation change detection was carried out to isolate hotspots (improved and degraded areas) from which soil samples were taken. ANOVA was used to test for difference in the measured soil properties at P < 0.05. Results show a general reduction vegetation cover, typical of a semi arid environment. SOC was generally low but higher in improved vegetation areas. While the degraded areas were purely sandy, improved vegetation was either Sand clay loam or Loamy sand, indicating a relationship between SOC and soil texture. It was concluded that desertification had effect on the selected soil properties in the study area. Embarkation on massive afforestation programme and protection of existing forest reserves coupled with environmental education for the farmers were recommended.


Introduction
The United Nations Convention to Combat Desertification [1] defined desertification as land degradation in arid, semiarid and dry sub humid areas resulting from climatic variations and anthropogenic activities.It is synonymous with land degradation in dry lands, manifesting in vegetation loss, water and wind erosion, salinization of irrigated and soil compaction [2]- [4].It is also seen as a change in soil properties, vegetation or climate, which results in a persistent loss of ecosystem services that are fundamental to sustaining life [5].
Dry lands of the world are of critical concern as they cover over 41% of the surface of the earth with over 2 billion inhabitants [3].The Africanarid and semiarid ecosystems are sensitive and dynamic, responding to fluctuations in climate and anthropogenic perturbations [6].In the early 1970s, countries in the West African Sahel suffered prolonged drought years, caused by long-term rainfall fluctuations [6] [7].The drought led to deaths of livestock and many inhabitants of the region [2].Vast acreage of land which previously was productive in the northern Nigeria during the past three decades has been ravaged by desertification, resulting in farmers abandoning their farmlands due to loss of sustainable livelihood [5] [8].It was affirmed that Nigeria loses 35,000 ha of land to desertification a year [9].
Soil quality and vegetation degradation in arid and semi-arid areas are the precursors and corollaries of desertification, and both are known to be highly correlated in most landscapes as their dynamics are being influenced by climate and anthropogenic activities [3] [10] [11].Although recent studies have discovered some level of greening in the Sahel proving possible resuscitation of the ecosystem [12], a hallmark of desertification is its irreversibility or difficulty of rehabilitation of the affected area as the degradation process continues even when external environmental stresses, such as heavy overgrazing, are removed [13].
Desertification consequences are far-reaching, and combating the menace demands proactive approaches [14].One of such approaches is to monitor the phenomenon via measuring land degradation and desertification processes [15] [16].Satellite images have been widely used to monitor vegetation production for the purposes of mapping land degradation [11].Studies have NOAA's operational environmental satellites and conventional data have been successfully used for drought detection and impact assessment in 20 countries [17].The study discovered that drought can be detected weeks earlier before arrival.Similarly, the nature of vegetation trend across the Sahel was investigated [12].The study successfully linked vegetation change (greening) of the study area to climate, land use management and migration.The footprint of vegetation dynamics in Kenya was traced using Advanced Very High Resolution Radiometer (AVHRR)/Moderate Imaging Spectro-radiometer (MODIS) NDVI data [18].Results proved that vegetation trend monitoring with remotely sensed data can give accurate indication of climate change.A time series of 1-km SPOT Normalized Difference Vegetation Index (NDVI) was used to simulate the spatial trend of desertification in Sokoto and its environs, Nigeria [19].Results showed the direction of desertification movement and that the inter-annual vegetation vigour exhibited a diminishing trend over the time series.
It has been affirmed that between 50% to 75% of the Border States of Nigeria are affected by desertification [9].Gudu, being a fringe local government area of Sokoto State, is in the heart of desertification zone of Nigeria.Incidentally, it is one of the areas that benefitted from the afforestation scheme of the Colonial government.Today, the area, including the forest reserves, has largely been affected by desertification.Consequently, agriculture as the sole means of livelihood for the ever increasing population is endangered.
Desertification is expected to get worse in Nigeria with the growing population and climate change [9] [20].Thus, for sustainable livelihood, timely information on vegetation and soil resources of the area is of essence.Remote sensing studies on desertification effect on the soil of the area are scanty in the literature.It is therefore hypothesized in this paper that desertification has effect on the Soil Organic Carbon (SOM) and particle size.Objectives include assessment of vegetation condition of the area, vegetation change detection and identification of hotspots and assessment of soil properties in the hotspots.

Study Area
The study area is located in Gudu Local Government Area of Sokoto State, which is situated in the extreme northwest of Nigeria and bordered by Niger Republic in North-West (Figure 1).The annual rainfall is around 500 mm experienced between June and September.The hot season peaks between March and April with a temperature of about 43˚C.The cold season which is accompanied by dust and sand laden winds is experienced between November and February, with the lowest temperature around 12˚C [20] [21].These climatic records have naturally conditioned the study area to Sudan savannah vegetation type composed primarily of shrubs and grasses.However, a large part of it is forested, gazetted since the colonial period.
Two major geological formations characterize the area: basement complex and sedimentary rocks.The earlier comprises old volcanic and metamorphic rocks which are mainly granite and metasediments, while the latter is the Gwandu formation overlying basement rocks in the lullmeden basin, extending from Sokoto, Niger to Mali.The sedimentary rocks comprise sandstones and clays with a high potential for ground water.The soil is characterized by sandy topsoil of Aeolian origin while the subsoil is clayey, derived from marine deposit.The area is located in the popular Sokoto Plains characterised by monotonous planation of 300 m (above sea level) with isolated mesas.
Based on the estimate of 2006, the population of the area was put at 42,172 people.Rain-fed agriculture is the main source of livelihood, and crops grown are largely grains.This also goes with light grazing by the inhabitants but the area is largely exploited by nomadic livestock rearing across the borders.Agriculture, bush burning, fuel-wood extraction, among other activities are a threat to the already fragile vegetation and soil resources.

Data Collection and Preparation
Two types of data were used for this study: 1 km SPOT Normalized Difference Vegetation Index (NDVI) archived as 10-day composites (Figure 2) and data on soil organic carbon (SOC) and particle size of the area.The NDVI composites covering a period between April 1998 and December 2014 were freely downloaded via www.spot-vegetation.com in zipped format.The zipped SPOT NDVI composites were extracted after which the study area was sub-mapped from each scene.Maximum Value Composite (MVC) was used to recompose the ten-daily composites first into monthly composites and later into annual composites as the process has the capability of further correcting the inherent atmospheric effect [22].Some composites could not be retrieved during data download due to technical failure; thus, annual composites were recomposed into biennial composites beginning from 1998 to 2014 (Figure 2).
Time Series Analysis was performed on the NDVI data series to determine the Principal Component (PC) and .By this method, the condition of the vegetation over time can be inferred.The OLS regression is a generalized linear modelling technique that may be used to model a single response variable (vegetation) which has been recorded on at least an interval time scale [25].It is represented mathematically as: Y = a + bX, where Y is the continuous response variable (vegetation), X is the continuous explanatory variable (time), a is the value of Y when X is zero, and b is the regression coefficient that describes the change in Y that is associated with a unit change in X [25].The extreme values of OLS image (both positive and negative) represent the hotspots, and NDVI profiles were calculated at selected points of such areas to show vegetation trends.Two types of hotspots were marked in the area: areas of increasing greening and areas of vegetation degradation.This classification was carried out using three standard deviation to the mean: =μ ± 3σ, where μ is the mean and σ the stan-dard deviation.The area coverage of these categories were estimated and compared.Soil samples were taken from these hotspots, and soil properties were compared using independent t-test.Data on soil was obtained by sampling using paired-plot approach [26] [27].Four well-spaced pixels were purposively selected from areas of improved vegetation and four other pixels selected in the degraded areas such that the latter pixels were closely juxtaposed to those in the improved vegetation areas as indicated by the analysis performed on the NDVI images.Thus, eight pixels were sampled in all.In each pixel (1 km 2 ), samples of soil were collected at four points and homogenized, out of which two handfuls were taken to the laboratory for analysis.Soil samples were taken to 30 cm depth using soil augers, and sampling was done at the time when biological activity was low [28], corresponding to the dry season.
Collected samples were bagged, labelled and transported to the laboratory for analysis.The samples were airdried for 48 hours, crushed and sieved with a 2-mm sieve.The samples were later used for testing for Soil Organic Carbon (SOC) and particle size (Sand, Silt and Clay).Particle size was analysed using the Bouyoucos hydrometer method [29] and classification was carried out using the textural triangle.The SOC was determined using the Walkley-Black method [30], while rating of SOC level was based on the scale in Table 1 [31].

Vegetation Condition of the Study Area
The condition of vegetation describes the overall health of the vegetation as enhanced or constrained by environmental and anthropogenic factors.This was calculated using Time Series Analysis (TSA) to generate the Principal Component (PC) (Figure 3) which presents the average vegetation condition of the area [23].In other words, it gives an illustration of the general spatio-temporal trend of the vegetation.Figure 4 shows the Loadings which describes the relationship between the PC and the time series images in Figure 2. The Loadings shows that the PC carries 92.32% of the information contents embedded in the time series, indicating a very high representativeness of the area.It also shows that there was a rising improvement in vegetation condition between the year 1998 and 2002, but the trend began to head down afterwards.This presents a general form of decline in vegetation which is likely a response to anthropogenic or climatic influence or both.[23].Both degraded and green areas (Figure 5) represent hotspots from which soil samples were collected.Two profiles, one for each of the degraded and green areas, were computed (Figure 6 and Figure 7).Figure 6 shows increase in vegetation, while Figure 7 indicates a downward trend of vegetation change.Area estimates of Figure 5 show that the Green Areas cover 10,600 hectares while areas of No-change cover 148,800 hectares; about 80,600 hectares of vegetation are areas of vegetation degradation within the time period.

Soil Property Status within the Hotspots
Table 2 presents the ratings of soil properties at the selected hotspots, while Table 3 represents the descriptive summaries and ANOVA test result of the soil properties at the sampled hotspots.Although values are rated low in both hotspots classes (Table 2), it can be seen that SOC is higher in A (Improved vegetation) than in C, and the difference is statistically significant (p < 0.05) (Table 3).Classification shows that the general texture in the degraded class is sand, while texture varies between Sandy clay loam and loamy sand in the green areas.The ANOVA test also indicates that that the textural differences are significant (p < 0.05).These findings show that improved vegetation can positively affect SOC and texture.

Vegetation Condition of the Study Area
Findings show that 92.32% information was contained in the Principal Component (PC) image which shows a very high representativeness of the vegetation of the study area.PCs account for majority of the variability in a dataset, giving room for the description of information with few variables [32] [33].PCA has been used to maximize change information in high variance components [34].A PC of 72% information content was observed while assessing vegetation change in Nigeria with advanced very high resolution radiometer NDVI [35].Further, a general form of decline in vegetation was discovered.It was argued that the north-western Sokoto which includes the study area is increasingly becoming desertified, manifesting a diminishing trend of interannual vegetation vigour [19].

Vegetation Change Detection
Results also show that more areas exhibited vegetative changes than vegetation improvement.Between 50% to 75% of the Border States of Nigeria are reportedly affected by desertification [9].Although some recent studies have discovered some level of greening in the Sahel proving attributed to management, migration, among others [12], desertification is expected to get worse in Nigeria with the growing population and climate change [9] [20].

Soil Properties at Hotspots
Results show that SOC is generally low in the study area and that SOC is higher in the improved vegetation than in the degraded areas.SOC is directly related to the amount of Soil Organic Matter (SOM) supplied to the soil [36] and this is related to the quantity and quality of the SOM as well as the environmental condition [37] [38].
The area is generally plagued by overgrazing, bush burning, fuel-wood extraction, near total removal of crop residue, naturally low vegetation cover, rapid mineralization due to high temperature and moisture conditions, intensive cultivation, among others [29] [39]- [41] such that even the forest reserves have become scrubby.All these potentially reduced the vegetation biomass and the addition of SOM to the soil.SOM improves the physical properties of soil, increases the Cation Exchange Capacity (CEC) and water-holding capacity of sandy soil, and it contributes to the structural stability of clay soils by helping to bind particles into aggregates [42]- [46].Conversely, unfriendly agricultural practices impact negatively on the terrestrial capability to sequester CO 2 [47].Further, results show that the general texture in the degraded class is sand, while texture varies between Sandy clay loam and loamy sand in the green areas, suggesting a mutual relationship between SOC and texture.Although the study area is characteristically sandy in nature consequent upon aeolian deposition, the loamy sand and sandy clay loam in the green areas is due to the enrichment from the soil organic matter which in turn is a function of vegetation improvement.
Studies have demonstrated soil textural changes in similar semi-arid environments [48] [49].Different LULC types reportedly have different effects on textural change through erosion prevention [50].This is plausible in a vegetated environment such as the green areas in this study.Soil texture is an important physical characteristic of soil which has great significance in land use management [51].Soils dominated by sand particles drain water very quickly but have lower fertility [38] [52], and this is where the soil of the study area belongs.

Conclusion and Recommendations
Based on the findings of this study, it is concluded here that the study area is getting dersertified and this has effect on the SOC and texture quality of the study area.Critical factors for soil property improvement in the desert area include increases in aboveground and belowground biomass through vegetation restoration with attendant improvements in soil chemical and physical properties and soil organic matter contents [53].Thus, it is recommended that government should embark on an afforestation drive to take care of the non-forested areas.Grazing reserves should be established to curtail wanton destruction of vegetation.Alternative fuel such as kerosene and

Figure 1 .
Figure 1.Map of the study area.

Figure 2 .
Figure 2. Biennial 1 km SPOT NDVI composites of the study area from 1998 to 2014.slope coefficient of the Ordinary Least Squares (OLS) regression between the values of each pixel overtime.The PC presents the overall vegetation condition of the area across the time series indicating the trend of vegetation overtime[23] [24].By this method, the condition of the vegetation over time can be inferred.The OLS regression is a generalized linear modelling technique that may be used to model a single response variable (vegetation) which has been recorded on at least an interval time scale[25].It is represented mathematically as: Y = a + bX, where Y is the continuous response variable (vegetation), X is the continuous explanatory variable (time), a is the value of Y when X is zero, and b is the regression coefficient that describes the change in Y that is associated with a unit change in X[25].The extreme values of OLS image (both positive and negative) represent the hotspots, and NDVI profiles were calculated at selected points of such areas to show vegetation trends.Two types of hotspots were marked in the area: areas of increasing greening and areas of vegetation degradation.This classification was carried out using three standard deviation to the mean: =μ ± 3σ, where μ is the mean and σ the stan- The slope coefficient of the Ordinary Least Squares (OLS) image shows pixels with different values, representing different degrees of vegetation change.Negative values represent areas where degradation took place in the time

Figure 3 .
Figure 3. Principal component image showing the average vegetation condition of the area.

Figure 4 .
Figure 4. Component loadings showing a downward relationship between the PC and each time series image.

Figure 5 .
Figure 5. Reclassified OLS image showing clearly patches of change classes.series while positive values show areas of vegetation improvement.Reclassification of these values using threestandard-deviation-to-the-mean algorithm (as the values are normally distributed) produced Figure 5.This map provides information on the spatial distribution of changes in vegetation.Green Areas are those places where vegetation appreciated overtime while No change areas are those places of minor vegetation variation.The degraded areas are places of negative vegetation change[23].Both degraded and green areas (Figure5) represent hotspots from which soil samples were collected.Two profiles, one for each of the degraded and green areas, were computed (Figure6and Figure7).Figure6shows increase in vegetation, while Figure7indicates a downward trend of vegetation change.Area estimates of Figure5show that the Green Areas cover 10,600 hectares while areas of No-change cover 148,800 hectares; about 80,600 hectares of vegetation are areas of vegetation degradation within the time period.

Figure 6 .
Figure 6.Profile of OSL at A.

Figure 7 .
Figure 7. Profile of OSL at C.

Table 2 .
Ratings of SOC and particle size at the hotspots.

Table 3 .
Descriptive summary and ANOVA result.