Pastures Degradation and the Relation with Pedo-Geomorphological Attributes in Watershed

The objective of this study was to map the levels of degradation of pastures and its association with chemical properties of the soil and relief. We selected a sub-basin located in the municipality of Viçosa, Minas Gerais, and we mapped the degradations through Ikonos images that were submitted for classification. Soil samples were collected and analyzed chemically; the results were transformed in spatial data using cokriging, based on variables of the digital elevation model. We associated the maps (cokriging) with data from exposure of the hillside, levels of solar radiation and levels of the pastures’ degradation. The pastures were mainly classified in moderate (44%) and very strong (30%) degradation levels. The areas with mild degradation tend to increase when there is a higher sum of bases, cationic exchange capacity and organic matter. The portions facing North and the increase of solar radiation showed more advanced stages of degradation.


INTRODUCTION AND OBJECTIVES
Although studies and field experiences are fundamental for the knowledge and mapping of soil attributes and classification (Buol et al., 2011), the line of research using pedrometry is gaining focus in soil studies due to the operational agility, low cost and reliable results (Cardoso et al., 2016;Ceddia et al., 2014;Hengl et al., 2015;McBratney et al., 2003). The quantitative techniques of spatial prediction are generally derived from geostatistics or from pedo-transference functions , McBratney et al., 2003Webster, 1994) and these were used in the mapping of attributes and soil classification in studies executed by Carvalho et al. (2002), Minasny & Hartemink (2011), Moore et al. (1993), Souza et al. (2010) and Souza et al. (2015). However, only a few studies correlate the results of the mapping with data of environmental issues.
From the most concerning environmental issues in Brazil, pastures degradation deserves special attention, and although there is not a study that quantifies such degradations in a national context, it is estimated that this problem affects approximately 70% of the pastures nationwide (Dias-Filho, 2011, 2014Peron & Evangelista, 2004). However, the recovery/regeneration process of the pastures depends on the recognition of several factors, for example: the degree of degradation, topography and the physical and chemical properties of the soil (Dias-Filho, 2014). Besides, the recognition of the spatial distribution of the pastures degradation is a basic step for its recovery (Nascimento et al., 2006) and remote sensing helps in the mapping process, with its effectiveness already proven in studies from Chagas et al. (2009), Eckert et al. (2015), Lamb & Brown (2001) and Nascimento et al. (2006).
The predisposition of the pastures to degradation may be related to natural relief conditions or attributes of the soils (Braz et al., 2013;Muller et al., 2001). Regarding the topography, the position of the surface of the hillside can have a negative influence on the characteristics of the vegetation, since in these particular areas there might be higher levels of light incidence, temperature and less water availability (Lima, 1986;Miranda et al., 2011). Studies conducted by Oliveira et al. (1995) and Matos et al. (2008), when analyzing hillsides, observed that the surfaces turned to the north in the southern hemisphere had a higher incidence of light, which led to lower humidity and smaller sized vegetation.
Hence, this study identifies levels of pastures degradation and their relation to spatial data of soil attributes and relief characteristics: (i) surface of exposure of the relief and (ii) levels of solar radiation.

Study area
The study area comprehends a sub-basin of Equídeo, 211 ha, located in the city of Viçosa, Minas Gerais (MG). The climate of the region is of the type Cwa (mesothermic, humid with rainy summers and dry winters), according to the Köppen classification, featuring water scarcity between the months of May and September. The annual averages of precipitation, relative humidity and air temperature are, respectively, 1,340 mm, 80 % and 19 °C (Paula et al., 2002).
There is a dominance of Gneiss rocks, that geomorphologically form elongated tips and a wavy relief. In general, the relief is very dissected, with ridges and narrow valleys. According to the American soil taxonomy system (Usda, 2014), the watershed has a predominance of Inceptisols + Ultisols combination in the hill slopes; Oxisols with variation of red and yellow oxisols based on the Brazilian soil classification (Santos et al., 2018) both dystrophic, especially in the hilltops; and Inceptisols + Aquents in the lower areas ( Figure 1). used was e-Cognition 3, based on the criteria established by Bhaskaran et al. (2010), Burnett & Blaschke (2003) and Monteiro et al. (2009). The recognition of the pattern of pasture degradation was based on proof gathered from field work, in which the levels of degradation were observed and associated to its spectral behavior on the image, the recognized spots were demarcated using a GPS. The levels of degradation were based on criteria determined by Spain & Gualdrón (1991) (Table 1).
The evaluation of the performance of the classifier was based on the Kappa (Ka) index (Congalton, 1991), which is based on the hypothesis test that performs a comparison between the categories of the reference map (validation samples) and the results of the classification. The validation samples were 150 spots randomly distributed on the Ikonos image, along with an observation of the category use for each spot. The results are summarized in the form of a significance matrix, in which the elements of the main diagonal indicate the significance of the classification and allows the calculation of the Ka for a given category, levels of quality were provided for each Ka value (Congalton, 1991;Landis & Koch, 1977).

Methodology procedures
The mapping of the levels of degradation was based on an object-oriented process of classification, performed using image bands (objects), which are grouped in categories with similar characteristics, such as: level of grey, spectral characteristics, texture, shape, compactness, information of the surrounding objects context, and that are organized according to hierarchy and semantics (Baatz & Shäpe, 2000, Benz et al., 2004. Ikonos images were used for the classification, with a spatial resolution of 1 m from 2015. The resolution of the image was changed to 10 m, based on the bilinear interpolation method (Schneider et al., 2003), to eliminate the influence of small objects with spatial signature similar to the pastures degradation levels, and also for ease of comparison with geomorphometric data of the Digital Elevation Model (DEM) of the SRTM (Shuttle Radar Topography Mission) images.
Initially, a visual interpretation of the Ikonos image was carried out to determine the thematic categories and further ahead these categories were segmented using the multiresolution algorithm, creating objects in different sizes, according to the shape, color and uniformity. The software process. The prediction was based on covariates from SRTM: altitude; slope inclination; wetness index; sediment transportation capacity; relief curvature; roughness index; surface exposure of the hillside; and global solar radiation levels (direct and dispersed). The geomorphometrical data were obtained using the Saga (System for Automated Geoscientific Analyses).
The information obtained from the soil analysis was used in the cokriging process, where geomorphometrical data from SRTM DEM was used as predictive covariables. The results were evaluated using a semivariogram, that indicates the greatest distance of which the data still has spatial correlation and a cross-validation process with five samples. After the cokriging process, the semivariogram model was checked to make sure it was adequate according to the normality, steady state, tendency and anisotropy of the samples. Afterwards, the model was evaluated using a rotation estimation technique (Goovaerts, 1997;Souza et al., 2015). The processes were performed with the R software using an R studio interface (Team R Core, 2015).
The only attributes used in the cokriging process and comparison with the levels of pasture degradation were the content of organic matter, the sum of bases and the effective cation exchange capacity -Ecec, which represent all the soil sorption complex. A Pearson correlation analysis was applied between Ecec, SB (sum of bases) and OM (organic matter) with the different levels of pasture degradation, whereby the different degradations were classified with values ranging from 1 (lowest level of degradation) to 4 (highest level of degradation).
The data of the surface exposure of the hillside along with the map of solar radiation were used in the analysis of the distribution of pasture degradation levels, whereby the levels were reclassified in three distinct radiation levels, in MWh/m 2 : low (0.371-0.580), medium (0.580-0.709) and high (0.709-0.817). The flat areas were eliminated due to their high radiation index and they were also not liable for comparison with the hillsides. The procedures used are summarized in Figure 2.
The SRTM DEM with a spatial resolution of 30 m was used to obtain geomorphometrical variables, which were processed to eliminate null and inaccurate values and then interpolated to 10 m, with the use of the tool Topo to Raster from the software ArcGIS®. Topo to Raster interpolation process was based on the algorithm defined by Hutchinson (1989), tested and evaluated in studies of Frey & Paul (2012) and Jing et al. (2014).
The interpolation technique (Topo to Raster) was created with the specific goal to covert vector data into hydrological models of consistent elevation. The method uses the computer efficiency for local interpolation, considering the reverse of the distance, still maintaining the superficial continuity of the global interpolation methods, such as cokriging and spline. The technique is efficient because it allows the use of archives of points with altitude values, contour lines and drainage, which assist in the reduction of errors and provides a hydrologically consistent model (Carmo et al., 2015;Hutchinson, 1989). Thus, data from the grid of points resulting from the DEM itself, contour lines (20 m) and drainage with a corrected direction of flow were considered (IBGE, 1979).
The soil sampling was collected in the field, and spots for sample collecting were determined through a hexagonal dot grill using the ArcGIS software. The software distributed the spots among the area considering an equal distance of 200 m, adding up to a total of 39 spots. Two soil samples were collected from each point (composite sample), in the ranges of 0 cm to 10 cm and 10 cm to 20 cm. The samples were chemically characterized according to the methods of Embrapa (1997), in which the following parameters were analyzed: pH in H 2 O, Ca 2+ , Mg 2+ , K + , Al 3+ , (H + Al) and organic matter. The data gathered was submitted to descriptive statistical analysis to represent the general characteristics of the soils.
Predictive covariables were obtained from the DEM and used to assist spatial prediction of soil attributes (Caten et

RESULTS AND DISCUSSION
The pastures from the studied area occupy a total of 133 ha and occur specially in plain and downslope areas, representing 60 % of the soil use of the basin. Typically, molasses grass (Melinis minutiflora Beauv.) prevails in the study area, considered a plant tolerant to low soil fertility and to moderate water scarcity (Nascimento et al., 2006). In a close-up field analysis, it was observed that a big part of the basin pastures is degraded at some level, especially due to neglect and, in fact, very few areas are actually used as pastures. The remains of the forest (secondary's successional state of the Atlantic rainforest) are associated to higher altitudes and hillsides ( Figure 3).
Ikonos satellite image was able to identify the degraded areas and the classification results were compatible with those recognized in field. Kappa index showed high quality standard when used to evaluate the levels of degradation (Table 2). Among them, the lowest index was designated to very high level of degradation (P4) and 6/13 Floresta e Ambiente 2019; 26(4): e20170616 Souza CMP, Francelino MR, Costa LM, Fernandes Filho EI other categories of urban use, whereby the reflectance of these areas may influence the classifier performance (Nascimento et al., 2006).
Considering specifically the pasture category, it was observed that 44 % was in a moderate degradation stage (P2) and 30 % in a very strong degradation stage (P4), which represents 27 % and 18 % of the basin drainage area, respectively (Table 3). The substantial area with a high degradation stage is possibly because this stage is closely related to the reflectance of the soil exposed areas or areas with similar reflectance (urban areas), consisting of a classifying error (Nascimento et al., 2006).
In general, the basin soils are acid, with values of pH ranging from 5.5 to 4.5 and they have low quantities of cations that are exchangeable such as K + , Ca 2+ and Mg 2+ , which in turn causes a low base saturation in the soil, averaging about 25 %. The levels of organic matter were also low, with an average value in the range of 4.0 dag kg -1 (Table 4). Table 2. Classification level of exactness, according to Kappa index value (Landis & Koch, 1977).  The maps resulted from the cokriging process show that the areas with a higher sum of bases are located in areas with a Inceptisols + Aquents association, presumably related to a bigger accumulation of organic matter in these areas, which is a characteristic of alluvial regions (Assine & Soares, 2004). Also, the soils of the region presented a dominance of kaolinite mineralogy (Campos et al., 2003), a factor that contributes for a bigger charge density and may justify the higher values for the SB. The effective cation exchange capacity Ecec revealed that the lowest values are those of the areas located in highest altitudes of the watershed, which is also the areas of Oxisols occurrence. The geomorphometrical aspects obtained from the SRTM assisted in obtaining the results of spatial attributes of the soil (Figure 4). The validation procedure indicated that the results were satisfactory for interpolation and that the soil sampling associated to the predictive covariables are adjusted to the land variability. The values of the Mean Squared Error (MSE) are close to zero, the values for the Root Mean Square Standardized Effect (RMSSE) are near 1 (Goovaerts, 2001), the coefficient of determination (R 2 ) for the samples of the cokriging process are above 0.54 and the spatial degree of dependency was considered strong (< 0.25) for the sum of bases and Ecec and moderate for organic matter (Table 5).   (Table 6), it is possible to observe a positive correlation between Ecec and the sum of bases, indicating that the increase of one attribute promotes the increase of the other, in a mutual manner. However, the organic matter (OM) attribute showed a negative correlation with the values of the sum of bases, but it is known that the OM may undertake a major role in the increase of negative charge in the soil (Oades, 1988).

Degradation and classes of land use Kappa index Variation (Kappa) Quality
The soil attributes and the pasture levels presented a negative correlation, indicating that an increase of degradation is followed by a decrease in the SB, Ecec and OM values. According to Heringer & Jacques (2002), the reduction of soil fertility is one of the factors of pasture degradation, as a result of inappropriate soil tillage, either due to an excess of pastures or a lack of liming and soil fertility and, at times, associated with the natural degradation tendency of the region (Dias-Filho, 2011). As an addition, the soils of the region naturally have low fertility, a fact that also affects pasture degradation.
Considering the relation between the soil attributes (SB, Ecec and OM) and the level of pastures degradation (Table 7), it is observed that there was a prevalence of milder degradation levels in the areas with a bigger sum of bases and also Ecec. In regard to the organic matter, it is known that it has a close relation with the pastures (D' Andréa et al., 2004;Feigl et al., 1995;Freitas et al., 2000), since these supply a large amount of carbon due to the root system. In the studied area, as the levels of OM in the soil increase, the pasture degradation levels decrease. This result is more evident on the stage P4 of degradation, which had a 72 % reduction in the degradation among the organic matter levels OM1 and OM2 (Table 7).
In general terms, the mild levels of degradation tend to increase, whereas the areas with higher levels of degradation tend to decrease. These variations in area occur as positive attribute variations are noticed.  below the tropics the hillsides of the relief facing north receive more energy compared to those facing south, resulting in a greater evapotranspiration potential, lower humidity content and, while studying pastures, Passos et al. (2017) noticed a higher degradation in these areas associated to an increase of the soil temperature.
Regarding the level of solar radiation on the study area, it was observed that as the radiation increases, the pastures degradation intensifies (Table 9). However, the biomass production of the pastures is dependent on the level of radiation, air temperature and water availability (Teixeira et al., 2011), in adequate amounts, and the type of fodder and soil must also be considered, since studies show that pastures completely exposed to the sunlight may present reduced productivity compared to partially shaded areas or seasons with lower radiation intensity (Paciullo et al., 2011;Santos et al., 2011). The relationship between the levels of degradation of the pastures with the exposed surfaces of the relief (Table 8) showed that 21.20 % of the degradation levels (mild and very strong) were found on the north of the surface. The pasture with a very strong degradation level (P4) had the second largest area in the surface of the relief. The hillsides on the southern hemisphere facing north and west tend to receive more irradiations, presenting lower levels of soil humidity, microbial activity and organic matter content, even if the soil's categories are not altered (Chagas et al., 2013;Passos et al., 2017).
The reduction of soil attributes and negative influence on plants characteristics in areas of the relief facing north and west have been noticed in studies by Lima (1986), Oliveira et al. (1995), Miranda et al. (2011). The authors explain that this may be because Table 8. Degradation levels of the exposed surfaces of the relief.

Levels
North  Table 9. Degradation levels of pastures in areas exposed to solar radiation. The result found indicates that soil conditions and the slope inclination of the studied area are influenced by the increase of temperature and has a direct impact on the sustainability of the pastures, since the increase of solar radiation may have a negative influence on photosynthesis and cause a rise of the plants respiration rate, reducing the net photosynthesis and biomass production.

Degradation levels
The areas with smooth and plain surfaces were not considered in the analysis, because they are exposed to solar radiation during most of the day. The study only contemplated the pastures located in slopes and hillsides, and these areas may be less fertile and have low humidity, which in turn prompts the development of fodder (Santos et al., 2014) and consequently pasture degradation.
On a local scale, solar radiation is a factor that has a very close relation to the topography, the radiation is modified according to inclination, surface orientation and the shape of the land (Kumar et al., 1997, Passos et al., 2017. Thus, the surfaces facing north, probably exposed to a higher level of solar radiation, presented higher levels of pasture degradation, confirmed by the fact that the increase in radiation also leads to an increase of pasture degradation with the soil characteristics of the studied area.

CONCLUSIONS
The classification of the pastures' degradation levels using Ikonos images was satisfactory, validated by field recognition and the Kappa index, which ranged from Very good to Excellent. The geostatistical procedures of cokriging used to map the soil attributes were effective because they correlated with predictive variables coming from the DEM, allowing to estimate values for nonsampled areas that had similar geomorphometrical characteristics to those of sampled soils.
The mapped attributes showed high relation to the soil categories. The areas with Oxisols featured lower values of SB, Ecec and OM. The lower altitude areas presented better chemical characteristics. It was noted that higher levels of pastures degradation had a direct association with lower values of Ecec, SB and OM.
North-facing slope surfaces presented pastures with a higher level of degradation when compared to those in other positions of the relief. The higher levels of solar radiation were also followed up with higher levels of pasture degradation. Considering the soil characteristics of the studied area (slopes, low soil fertility, precipitation, fodder types), such results indicate that the long solar incidence due to relief position or the radiation incidence itself are aspects that have a negative influence on the pastures.