MULTIVARIATE ANALYSIS OF CHEMICAL AND PHYSICAL ATTRIBUTES OF QUARTZIPSAMMENTS UNDER DIFFERENT AGRICULTURAL USES

The inadequate agricultural management in the semiarid region of northeastern Brazil, associated with the fragility of most soils, has caused degradation of the quality of this resource. This study aimed to evaluate the quality of Quartzipsamments under different uses through the multivariate analysis of chemical and physical attributes. Five areas were selected: four of them with agricultural use (onion, banana, cassava, and corn) and the other with Caatinga. Soil samples were collected at depths of 0.00–0.10, 0.10–0.20, and 0.20–0.40 m for chemical and physical analyses. The data were submitted to descriptive analysis and then to the multivariate analysis using the principal component analysis (PCA) technique and Ward’s clustering method. PCA and clustering techniques allowed the identification of different production systems in subsurface layers, suggesting the effect of fertility on the differentiation of areas with onion and banana, as well as the clustering of cassava and corn cultivations. Physical attributes and sodium content of the surface layer of the area with onion showed a tendency of compaction and salinization process of this environment. In general, fertility was the main responsible for the improvement of soil quality with agricultural use when compared to the Caatinga.


INTRODUCTION
The implementation of agricultural exploitation systems may cause negative impacts on the soil due to the adoption of a model of agriculture that does not prioritize the rational use of natural resources.The occupation of soils considered of low aptitude without the previous assessment of their use capacity can lead to the adoption of inadequate management systems and their degradation, especially when fragile soils, such as sandy soils, are inserted into the productive process (Sales et al., 2010;Mota & Valladares, 2011).However, the agricultural use of Quartzipsamments is common in the semiarid region of northeastern Brazil.According to Jacomine (1996), this soil class covers around 9.3% of the area of the region.
Although there is information on soils from the semiarid region, such as exploratory surveys and other studies of this theme (Jacomine, 1996;Santos et al., 2012), it is always valid to carry out studies on the characterization of different representative soils on a more detailed scale (Mota & Valladares, 2011).Soil characterization studies in regions little explored allows providing more accurate information on the different soil classes throughout the Brazilian territory, as well as systematizing this information on soil properties, which may be a subsidy for the development of management practices, sustainable use, and recovery of degraded areas (Santos et al., 2012).
The detailing of potentialities and fragilities is of extreme importance to the planning to minimize impacts observed continuously in the region (Corrêa et al., 2010).These impacts are considered the result of a combination of agriculture management and environmental characteristics (Corrêa et al., 2010;Mota & Valladares, 2011;Ursulino & Moreno, 2014), which together lead to and accelerate the development of physical and chemical degradation.These processes are mainly due to the use of irrigation systems (furrows and flooding) of low efficiency, overuse of chemical fertilizers, and, in many cases, unfavorable natural drainage conditions.These factors, coupled with the high evapotranspiration demand and low rainfall index (Corrêa et al., 2010), Engenharia Agrícola, Jaboticabal, v.39, n.4, p.457-465, jul./aug. 2019 have accelerated the salinization process of these areas, negatively affecting changes in physical and chemical attributes.
Thus, it is necessary to carry out studies to better understand changes caused by soil use, aiming at maximizing production and avoiding the degradation of agricultural soils since they generate essential information for the management and conservation of tropical sandy soils incorporated into intensive productive processes (Sales et al., 2010).In this sense, according to Corrêa et al. (2010), changes in soil attributes can be observed after agricultural use and agricultural uses can be compared.Also, sustainability can be evaluated from use evaluations with soil attributes as indicators, aiming at the evaluation of productive systems to adapt systems or propose more sustainable uses.
In this context, the multivariate analysis technique can be used as a tool to explain the maximum correlation between soil characteristics and indicate those that most contribute to soil characterization and/or soil alteration (Freitas et al., 2014;Oliveira et al., 2015).Several studies have applied the multivariate technique for soil quality analysis (Pragana et al., 2012;Freitas et al., 2014;Arcoverde et al., 2015;Oliveira et al., 2015).
Therefore, this study aimed to evaluate the quality of Quartzipsamments under different uses through the multivariate analysis of chemical and physical attributes.

MATERIAL AND METHODS
The studied Quartzipsamments are located in the municipalities of Casa Nova, Remanso, and Sento Sé, on the banks of Sobradinho Lake in the state of Bahia.According to Köppen classification, the climate BSwh′ (warm and semiarid) predominates in 75% of the area, with annual precipitation varying from 500 to 900 mm; the arid climate can be found in the other 25% of the area, with annual precipitation lower than 500 mm.The predominant vegetation is the hyperxerophilic Caatinga.The surroundings of the Sobradinho Lake is characterized by intense agricultural activity, standing out the irrigated agriculture and vegetable cultivation, mainly onion.
Five areas with soils classified as Quartzipsamment, located on the banks of the lake and under different uses, were used in this study: a) onion cultivation (O), b) banana cultivation (B), c) cassava (CA), d) corn cultivation (CO), and e) Caatinga (CAA) (Table 1).Soil tillage for planting in the areas with agricultural use was performed with a disc plow and heavy disc harrow; liming and fertilization was also carried out in the areas.The five areas presented Quartzipsamments with a textural class predominantly loamy sandy at depths of 0.00-10, 0.10-0.20,and 0.20-0.40m, whose mean particle size is shown in Table 2. Disturbed and undisturbed soil samples were collected at depths of 0.00-0.10,0.10-0.20,and 0.20-0.40m.Each area was divided into three sub-areas of sampling to collect ten simple disturbed soil samples, which composed a composite sample, totaling three composite samples per depth and sub-area.Trenches were opened in the middle of each sub-area to collect undisturbed soil samples.
The mean and standard deviation were calculated by the classical statistics for each studied attribute.The multivariate technique was used through the principal component analysis (PCA) and hierarchical cluster analysis.
PCA was performed aiming at reducing a large number of variables into a more significant set (represented by principal components), identifying the variables belonging to the components, and determining how much each variable explains each component to evidence attributes associated with the changes in the areas (Freitas et al., 2014).The minimum of principal Engenharia Agrícola, Jaboticabal, v.39, n.4, p.457-465, jul./aug.2019 components with at least 70% of the total data variance was adopted (Burak et al., 2012).Also, the principal components were those with a value of a correlation coefficient between the principal component and the variable above 0.65 in absolute value, as described in Arcoverde et al. (2015).
The hierarchical cluster analysis was performed by calculating the Euclidean distance between the accessions for the set of 12 attributes using Ward's algorithm to obtain clusters of similar accessions.This analysis was used to verify the similarities between analyzed attributes and studied areas from homogeneous clusters represented in a similarity dendrogram.

RESULTS AND DISCUSSION
The agricultural uses promoted modifications in the mean values of the physical and chemical attributes of Quartzipsamments at the three studied depths (Table 3).Engenharia Agrícola, Jaboticabal, v.39, n.4, p.457-465, jul./aug.2019 Soil physical attributes showed higher Ds values at the three studied depths in the soil with Caatinga; the area with onion showed the highest values among the evaluated production systems (Table 3).Lopes et al. (2012) verified higher values of soil density in the surface soil layer of the Caatinga area, which was attributed to the presence of crop residues in this area.
Soil compaction may occur in the area under onion cultivation due to more continuous use and because it is an irrigated crop with a short cycle.Besides, mechanical management practices such as plowing alter the original soil structure (Cortez et al., 2011;Dantas et al., 2012).Thus, systems that favor soil structuring, such as those that increase organic matter contents, should be adopted in areas with degradation signs (Stefanoski et al., 2013;Sales et al., 2016).
All soil uses had ratios of micropores higher to those of macropores at the three depths (Table 3).Despite this, ratios lower than 2:1 were verified in banana and corn cultivations and Caatinga area at depths of 0.00-0.10 and 0.10-0.20 m, and in banana, cassava, and corn cultivations at a depth of 0.20-0.40m (Table 3).This result is possibly related to the predominance of fine and very fine sand fractions over the coarser fractions, a striking feature in some sandy soils of the semiarid region (Santos et al., 2012).
The depth of 0.0-0.10m showed acidification and loss of exchangeable bases (K, Ca, and Mg), mainly in corn, cassava, and banana cultivations.The OM contents were lower when compared to the Caatinga regardless of the soil use, while P contents remained close, except for cassava (Table 3).
The depth of 0.10-0.20 m presented exchangeable base contents close to those observed in the soil under Caatinga, except for Ca (Table 3).In addition, OM contents were similar to that found in the Caatinga only for onion and banana cultivations.A significant increase in assimilable P contents was observed in soils cultivated with these crops (Table 3).
The depth of 0.20-0.40m had an increase in K, Ca, and Mg contents in relation to the Caatinga, with the maintenance of pH values (Table 3).It is probably due to the management used in each crop, with low or without the addition of inputs in corn and cassava cultivations and high addition of manure and nutrients, mainly phosphate, in the areas with onion and banana.These factors, together with the irrigation (surface) and inadequate management may have contributed to nutrient losses due to export and leaching, taking into account the sandy soil texture.Morais et al. (2014) also observed pH values above neutrality in soils under irrigated banana cultivation and native forest (Caatinga), with low availability of phosphorus and micronutrients under these conditions.These authors also emphasized that the increased pH in cultivated areas is closely related to the fertilizer hydrolysis process, such as urea, which is widely used in onion cultivation.
Increased P in cultivated areas may be explained by the residual effect of phosphate fertilizers over the years although its high concentrations in these areas are associated with overestimation of P availability because this form of P is not readily available to plants, except for those that acidify their rhizosphere, as reported by Morais et al. (2015).According to these authors, the increased P in depth, especially in surface soil layers (Table 3), can be attributed to the high moisture and descending flow of water, which promote an increase in the diffusion and percolation of P.
In general, low CECeff values were observed in the soils due to their sandy constitution and low OM contents (Table 3).According to Santos et al. ( 2012), it favors strong leaching of elements during rainy seasons, leading to their lower exchangeable contents in the soils.Besides these factors, the use of low-efficiency irrigation systems (furrows and flooding) in the region promotes higher losses of soil elements by their removal.
Soils presented, in general, low to medium contents for exchangeable Ca, Mg, and K, mainly under cassava and corn cultivations (Table 3), as they are areas of family production with a low or no contribution of external inputs.Similarly, assimilable P contents in areas with cassava and corn were considered low, unlike that observed in the area with onion at a depth of 0.10-0.20 m and with banana at depths of 0.10-0.20 and 0.20-0.40m (Table 3).
The principal component analysis (PCA) was carried out aiming at evidencing the attributes associated with changes in the areas due to the used agricultural management.Two components were extracted, which accounted for 68.67, 77.47, and 71.58% of the total variability of the data at depths of 0.00-0.10,0.10-0.20,and 0.20-0.40m, respectively ( The principal component 1 (PC1) is more associated with groups of chemical attributes, while the principal component 2 (PC2) is more related to soil physical attributes when considering the three depths (Table 4).It means that soil chemical characteristics were more changed by the management, being more sensitive and having the possible use as quality indicators.The PC1, generated by soil attributes at a depth of 0.00-0.10m (Table 4), explained 50.36% of the total variance of the studied attributes, being those with the highest load (≥ 0.65 in absolute value) an indication soil compaction (Ds), soil reaction (pH), nutrient availability (OM, K, Ca, Mg, CECeff), and salinization (Na) (Table 4).In other words, these attributes discriminated better the differences between areas.
These results corroborate a discussion of Table 3, which shows that, in general, the management of the Quartzipsamment provided a reduction in its natural density and fertility, as well as a slight increase in the exchangeable sodium (Na) content of onion cultivation.According to Dantas et al. (2012), the decreased soil density (Ds) in uses other than onion cultivation can be attributed to a lower soil disturbance and higher return of residues to the soil.
On the other hand, the high Ds values observed in the Caatinga area indicate an advance of degradation problems.Ursulino & Moreno (2014) evaluated the physical quality of an Entisol in the semiarid region of northwestern Brazil and found mean Ds values close to 1.50 Mg m −3 , which is considered high.This result was attributed to deforestation in the area, which causes pluvial erosion effects, causing surface sealing, which increased soil density.It is similar to that observed in the Caatinga area of the present study (Table 3).
Also, a higher Na content was observed in the surface soil layer of the area with onion when compared to the Caatinga area, with a mean value of 125% (Table 3).Na accumulation is mainly related to the management of irrigation water (Morais et al., 2015).Thus, a continuous accumulation should be prevented to avoid affecting the yield of cultivated areas, being required the effective management of the irrigation water depth (Dantas et al., 2012).
Soil porosity (Ma and Mi) and chemical quality were more sensitive in discriminating banana, cassava, and corn uses possibly due to the management systems used in these areas, which contributed to an intermediate change.
In this case, the values of fertility-related attributes were usually similar and thus, different from those observed in onion cultivation and Caatinga (Tables 3 and 4).
Therefore, according to the results found in this layer (Tables 3 and 4), soils cultivated with onion and native vegetation (Caatinga) were strongly discriminated by Ds, OM, CECeff, K, Ca, Mg, and Na when compared to soils with banana, cassava, and corn.In addition, macroporosity was more sensitive in discriminating soils with cassava and Caatinga, while soil microporosity (Mi) influenced soils cultivated with banana, onion, and corn (Tables 3 and 4).
At a depth of 0.10-0.20 m, PC1 explained 54% of the total variance of the data and attributes related to soil macroporosity (Ma), nutrient availability (P, K, Ca, Mg, and CECeff), and accumulation of salts (Na) were more sensitive in distinguishing the areas (Table 4).PC2 explained 23.47% of the total variance in the data, being composed by indicator variables of compaction (Ds and Mi) and soil reaction (pH and H+Al) (Table 4).
The attributes Ds, Ma, and pH strongly discriminated soils with cassava and Caatinga, with a higher proportion of macropores in relation to micropores and lower values of pH and Ds in cassava cultivation, but with opposite results in the Caatinga (Tables 3 and 4).On the other hand, the higher proportion of micropores in relation to macropores discriminated soils with onion, banana, and corn, while P, K, Ca, Mg, Na, and CEC contents strongly discriminated onion cultivation and moderately discriminated the soil with banana (Tables 3 and 4).
Engenharia Agrícola, Jaboticabal, v.39, n.4, p.457-465, jul./aug.2019 At a depth of 0.20-0.40m, PC1 explained 44.19% of the total variance of the data and was composed by the variables OM, P, Ca, Mg, and CECeff, while PC2 explained 27.39% of the total variance and was composed of physical (Ds, Ma, and Mi) and chemical variables (K and H+Al) (Table 4).High P contents were mainly observed in this layer, while Ca and Mg discriminated the area with banana and base accumulation, related to an increase in Ca and Mg contents, discriminated cassava and corn areas (Table 3).
The variables Ds, Mi, and H+Al were strongly sensitive in discriminating soils cultivated with onion and Caatinga, also moderately influencing banana cultivation (Tables 3 and 4).Ma significantly discriminated the uses of cassava and corn and moderately banana cultivation, which was grouped alone in cluster analysis (Tables 3 and 4).
Caatinga and onion areas presented similarities regarding their physical and chemical attributes (Table 3).These areas had the highest Ds values and higher proportions of micropores in relation to macropores at the three depths, as well as the highest values of OM, CECeff, and nutrient contents in the surface soil layer, but reducing in depth (Table 3).Corn and cassava areas were usually more homogeneous, mainly due to the lower mean Ds values and low OM values in the first two layers, which may have promoted the lower fertility (Table 3).Finally, in isolation, soil use with banana, mainly in subsurface layers, presented the highest P and OM values, which are possibly associated with the higher fertility (Table 3).
These results confirmed the importance of the physical attributes Ds, Ma, and Mi as soil quality indicators of the northeastern semiarid region, as reported by Arcoverde et al. (2015).These authors highlighted Ma as a sensitive attribute for discriminating the physical quality at a depth of 0.00-0.10 and 0.10-0.20 m, while Ds and Mi were more important at a depth of 0.20-0.40m, thus being considered good soil quality indicators.Lopes et al. (2012) observed higher OM values in the surface soil layer of a Caatinga area when compared to an area with melon, which was attributed to the presence of dead coverage in this area as depth decreased.Also, a significant increase in OM was observed at a depth of 0.20-0.40m when compared to the soil under Caatinga, which is similar to the result obtained in the area with banana of this study.
The results in the corn and cassava areas, especially in the surface soil layers, are possibly due to the inadequate irrigation management associated with the sandy texture and fertilization, which was insufficient to meet crop demands.Thus, the low values of exchangeable bases are due to the leaching, which promotes the accumulation of these elements at a depth of 0.20-0.40m, acidifying the surface layers, revealed by very low pH values.Xavier et al. (2012) observed that the low pH values in the surface soil layer are due to the low recharge of soluble salts from fertilizer application or even the entry of these ions from irrigation water.
From the confirmation of the soil attributes that promoted differences between areas through the principal component analysis, the hierarchical cluster analysis was performed, according to Ward's method, aiming at classifying the soil attributes and principal components into homogeneous clusters (Freitas et al., 2014).In this sense, the cluster analysis technique for the set of studied variables allowed complementing efficiently the results obtained with the use of the principal component analysis (Figure 1).
With the use of significant attributes in PC1, Ward's method was used with the similarity measure given by the Euclidean distance aiming to group them, defining the cut-off distance equal to 50% of the maximum Euclidean distance.In this analysis, onion, banana, cassava, corn, and Caatinga areas were grouped according to their degree of similarity, classifying them into homogeneous clusters (Oliveira et al., 2015).
The combined use of physical and chemical attributes at a depth of 0.00-0.10m (Figure 1A) allowed the data to be sorted into two clusters: the first one encompassing cassava, corn, and banana areas and the second cluster encompassing the data of Caatinga and onion areas.Production systems grouped into the first cluster are those with the lowest contribution of inputs and those with the highest chemical changes in relation to the Caatinga soil, with reductions in pH values and exchangeable base contents.
The clustering using the data from a depth of 0.10-0.20 m (Figure 1B) provided the formation of two clusters: the first one formed by Caatinga, corn, cassava, and banana areas, and the second one with the onion area.The separation of the onion area from the other production systems is mainly due to differences in P content observed at this depth.
Engenharia Agrícola, Jaboticabal, v.39, n.4, p.457-465, jul./aug.2019 The cluster analysis from the data of soil chemical and physical attributes at a depth of 0.20-0.40m allowed a clear division into two clusters.The first cluster consisted of the data of the area with banana, while the second cluster was composed of the data of the other areas (Figure 1C).In this case, soil P content was probably the attribute that differentiated the banana production system, in addition to Ca and K contents.
A higher similarity was observed between areas with onion and Caatinga at depths of 0.00-0.10 and 0.20-0.40m, with the smallest changes in soil fertility due to onion management.Areas with corn and cassava cultivation showed similarity at the three depths because the management systems are similar, with no replacement of nutrients removed from the soil by crops, causing acidification in the most superficial layer and base accumulation in the deepest soil layer (Figure 1).
The results of the cluster analysis of the three layers for the analyzed physical and chemical attributes (Figure 1) showed that their changes allowed a marked differentiation of clusters, with clear distinction between the onion cultivation of the other areas at a depth of 0.10-0.20 m and the area with banana at a depth of 0.20-0.40m, suggesting a higher effect of management on the subsurface soil layers.
Thus, PCA and cluster analysis showed clear changes in soil physical and chemical attributes due to the specific management adopted in the studied areas, thus being efficient techniques to study the quality of Quartzipsamments from the identification of indicators sensitive to the management.
In general, the differentiation of the areas according to the management was evident despite being very similar due to the predominance of loamy sandy texture.Thus, areas with cassava and corn cultivation had a higher deficiency and limitation from the point of view of soil fertility in the surface layers when compared to the Caatinga area, although with base accumulation due to leaching in depth.This behavior was related to the management of fertilization carried out in these areas.On the other hand, banana cultivation had an improvement in soil quality in depth as fertility increased probably due to higher P, Ca, Mg, and OM contents, which allowed differentiating this use from the others.
Changes in soil physical attributes and Na accumulation in the surface layer show the beginning of degradation processes related to soil compaction and salinization, despite the similarity regarding soil quality of Caatinga and onion areas.

CONCLUSIONS
Principal component analysis and hierarchical cluster techniques allowed the identification of different production systems in subsurface soil layers, suggesting the effect of soil fertility on the differentiation of areas with onion at a depth of 0.10-0.20 m and banana at a depth of 0.20-0.40m, as well as the clustering of cassava and corn cultivations at the three studied layers.Therefore, soil fertility is the main responsible for soil quality.
The physical attributes Ds, Ma, and Mi and Na content of the surface layer of the area with onion showed a tendency of compaction and salinization.The former variables were also sensitive, together with soil fertility, in differentiating the studied areas.
In general, an improvement in soil quality with agricultural use was observed in depth when compared to the Caatinga.

TABLE 2 .
Particle size composition of Quartzipsamments by depth.