ABSTRACT
Soil attributes are easily modified for different agricultural purposes, requiring the adoption of appropriate practices, according to the local particularities in order to maintain its production capacity. Given the above, this study aimed to assess soil physical attributes under different agricultural uses conditions in an Oxisol. The survey was conducted in Brazil with the following treatments: ​​Cassava monoculture area, beans monoculture and native forest. Results were interpreted with multivariate analysis. Native forest differed from the other agricultural uses due to surface consolidation, lack of soil preparation practices and increased organic matter input, having higher water retention. The conclusion was that there were statistical irregularities in soil attribute tendencies. Different agricultural uses interfered in the physical attributes when compared to the native forest, which produced better results, followed by bean and cassava monocultures. The most significant physical attributes to distinguish agricultural uses were: Aggregation; weighted average diameter; soil density; particle density; total porosity; field capacity and available water.
Key words: Soil properties, density, porosity, water retention.
Soils are three-dimensional, natural and dynamic bodies. They are formed on the Earth’s surface through climate and organisms environmental factors action on the source material in function of relief, time action and
spatial variation, according to these factors combination possibilities of these factors (Kämpf and Curi, 2012).
Under natural conditions the soil is in dynamic equilibrium and is resistant to changes (Goedert and Oliveira, 2007). However, human action has promotes attribute changes (Silva et al., 2007), causing soil natural fertility decline (Cordeiro et al., 2004), as its structure is easily changed. Thus, other properties are also degraded, such as soil density, micro, macro and total porosity, organic matter, water retention and infiltration in the soil profile (Reichert et al., 2009; Ferreira et al., 2010; Cunha et al., 2011).
Therefore, inappropriate land use causes changes in physical and hydraulic properties. The magnitude of these changes magnitude varies depending on the soil type, weather and agricultural crops management (Kay, 1990).
In general, arid and semiarid regions soils are naturally fragile because of their low water storage. In addition, plant root system growth is also physically limited in these soils (Cardoso, 2002), increasing its susceptibility to human action improper use practices.
Soil studies in Rio Grande do Norte state western region and in the municipality of Martins, RN state, are scarce. Therefore, physical attributes quantification is necessary to establish appropriate farming practices establishment, in order to seek better conditions in relation to the soil water dynamics, aeration and structure, which are dynamic attributes with spatial and temporal variability.
Given the above, this study aimed to assess soil physical properties under different agricultural uses in Martins, RN. Multivariate analysis was used as a tool to interpret results.
Research was conducted in Martins, RN, Brazil. Martins is located in the Rio Grande do Norte state western meso region and in the Umarizalmicro region, in the following geographical coordinates: 6° 05´ 16´´ South, 37° 54´ 40´´ West. It is located in the Borborema plateau. However, its relief further comprises the Depressão Sertaneja, covering an area of 169.47 km². Climate is classified as tropical rainy (Aw) according to Köppen classification, with average rainfall of 1133.8 mm, whose rainy season is from January to June, the average annual temperature is of 25°C and subperennial forest vegetation type, which is associated with the Caatinga vegetation.
The following treatments were assessed: cassava monoculture (CASS), bean monoculture (BEAN) and native forest (NF).
The cassava monoculture area covered about 1 ha. The site was cleared in 2003. Plant residues were burned in small piles, with plowing and harrowing being conducted afterwards. These practices were always adopted before cassava planting, about once every two years. Sampling was performed after burning.
The bean monoculture area covered 1 ha, and like the cassava area, it was deforested in the same year (2003). Plant residues were burned and soil tillage was held afterwards through plowing and harrowing, the latter being held once a year, before sowing.
Through visual analysis by the farmer, cattle manure was applied to improve soil fertility. However, such practice was not performed during sampling.
The native forest area was used as a reference to compare with the other agricultural uses. The vegetation is of subperennial type, consisting of broad-leaved trees, with relatively slender and dense trunks, which is typical of tropical rainy climate zones. The most commonly found species are: Psiium firmum, Talisia esculenta, Enterolobium contortisiliquum, Pithe colobium polycephalum, Sideroxylon obtusifolium, Helicostylistomenosa, Hymenaea courbaril, Dipteryx odorata, Mimosa sepiaria, Copaifera langsdorffii, among others.
Agricultural uses textural classification was: sandy clay to native forest (NF), cassava monoculture (CASS) and bean monoculture (BEAN). Two soil profiles were opened in a representative area of the study sites. Soil samples with deformed structure in their respective horizons were collected for physical analysis, as well as diagnostic horizons identification and classification, according to the Brazilian Soil Classification System (Donagema et al., 2011; Santos et al., 2013).
Profile 1 showed a textural classification ranging from sandy clay on A and BA horizonsto clay in the Bw horizon. Profile 2 ranged between sandy clay in the A horizon to clay in the Bw horizon. Both were classified as Oxisol.
In order to assess soil physical properties, four representative points in each area were sampled to collect disturbed and undisturbed samples, in depths from 0.00 to 0.30 m. Subsequently, two representative profiles of the areas under study were opened for soil classification of soil, according to methodology described by the manual of methods and soil analysis (Donagema et al., 2011).
In order to collect soil samples with undisturbed structure an Uhland type device was used, with cylindrical metal sampler dimensions of 0.05 m diameter and 0.05 m height was used. Physical attributes, retention curve, soil density, total porosity, macro and microporosity were assessed. Eight samples were collected by point and 32 samples were collected for each area under study, totaling 128 samples. After collection, samplers were coated with aluminum foil and taken to the laboratory with care, in order to maintain the soil structure and moisture.
Deformed structure samples collection were carried out in four representative points of each area, with three repetitions, totaling 48 samples, Collection was conducted with the aid of a straight shovel, with samples being transferred to properly labeled plastic bags. Subsequently, samples were used to determine textural classification and particle density.
Textural classification (granulometric analysis) carried out by pipette method using chemical dispersant (sodium hexametaphosphate). Soil density (SD) was determined by volumetric ring method and particle density (PD) was determined by volumetric flask method, using alcohol. These determinations were performed following Embrapa soil analysis methods manual (Donagema et al., 2011).
In order to determine aggregation and average diameter, wet screening method was used. Five samples of 25 g for each agricultural use were used, with a sample being used for moisture determination (correction factor). Each sample was transferred to a filter paper disk located on the upper sieve (2.00 mm) and spread over the entire filter paper surface. Moistening was conducted through capillarity (4 min) in a four sieves set with the following mesh diameters: 2.00; 1.00; 0.50 and 0.25 mm. After the required time, the filter paper was removed. Sieve sets agitation was performed in avertical oscillation apparatus (42 oscillations/min) during four minutes inside barrels with water. Fractions retained on each sieve were transferred to beckers that were previously treated and identified, and sent to an oven at 105°C for 24 h. After dry matter obtainment and sand content (sodium hexametaphosphate) deduction, water stable aggregates (aggregation) percentage and weighted average diameter (WAD) were determined for each class. WAD was obtained through Equation 1, proposed by Castro Filho et al. (1998).
.png)
Where: WAD–weighted average diameter, mm; xi–classes average diameter, mm; wi - each class proportion in relation to the total.
For TP, macroporosity and microporosity determination, tension table was used. For microporosity, a tension of 6 kPa was used. Macroporosity was calculated through the difference between total porosity and microporosity.
The graphical relation between matric potential and soil water content is calledsoil water characteristic curve, or moisture retention curve, that is, matric potential (jm) x water content (θ). For retention curve obtainment, tensions of 0; 2; 6; 10; 33; 100; 300; 500 and 1500 kPa were used. In the application of 0; 2; 6 and 10 kPa tensions, the tension table was used. Low tension cameras were used for 33 and 100 kPa, and high tension cameras were used for 300; 500 and 1500 kPa.
Soil water retention curves adjustment was made based on van Genuchten’s equation (1980), (Equation 2) using the SWRC software, developed by Dourado Neto et al. (1990). The equation considers the matric potential (jm) as an independent variable, and water content (q) as dependent variable:
Where: qr -residual water content, m3 m-3; qs- saturated water content, m3 m-3; |jm|- matric potential, kPa a, m, n –Equation empirical parameters.
As it is a theoretical parameter that varies between 10 kPa (sandy soils) and 33 kPa (clay soils), field capacity value (FC) was determined according to the textural classification. Permanent wilting point (PWP) was obtained with 1500 kPa tension, and available water (AW) was obtained through the difference between FC (at 10 kPa tension) and PWP.
Soil physical attributes average values under different agricultural uses (NF, CASS and BEAN) were interpreted through ANOVA, which was performed with Tukey’s test (p<5%), descriptive analysis (average, minimum, maximum, variance, standard deviation and coefficient of variation) In addition, multivariate statistical analysis tool was used in Principal Component Analysis (PCA), using the GENES program.
Soil physical attributes were measured in different unit systems, with data standardization being necessary, since variance is influenced by the attributes in the measurement units in question. Correlation matrix was established after data standardization, in order to verify the percentage and the degree of importance of these correlations, with values that were higher or equal to 0.7 being considered.
Through soil physical attributes analysis of variance under different agricultural uses (AU analysis of variance, it was shown that there was no statistical difference for AGREG, SD, PD, MACRO, MICRO, FC and AW there was no statistical difference. WAD statistically differed at1% probability and TP and PWP were statistically different (p<0.05) (Table 1). These differences are probably due to soil structural arrangement being a dynamic, functional, functional and complex property that is easily modified property by land use (Carneiro et al., 2009; Kämpf and Curi, 2012).
Evaluating SD, PD, TP and PWP attributes coefficients of variation (CV) showed low values (CV lower than 10%); as for aggregation, WAD, MACRO and average FC, CV was between 10 and 20%. For MICRO and AW, CV was very high (higher than 30%), according to Pimentel (2009).
The correlation matrix obtained with the physical attributes studied in the different agricultural uses (NF, CASS, BEAN), in the municipality of Martins-RN were highly correlated (values above 0.7), showing interrelation in the different agricultural uses (Table 2). High and positive correlations between aggregation (AGREG), weighted average diameter (WAD) and microporosity (MICRO) were found.
High positive correlations were found between SD and TP for MACRO, FC and AW (except for TP, in which high positive correlation was found for MICRO). These correlations can be justified due to SD and TP not being limiting factors for other attributes, with soil density values ranging from 1.14 to 1.28 g cm-3 and total porosity ranging between 37 and 47% (Table 3). According to Borges et al. (1997), SD values above 1.62 g cm-3 were partial obstacles for root growth while working in an Oxisol with medium texture at the Triângulo Mineiro region.
In the native forest, higher AGREG, WAD, SD, TP, MACRO, FC and AW average values were found (Table 3), what can be attributed to its surface higher consolidation due to tillage practices absence and higher organic matter input provided higher average values. Data tendency is in accordance with Corrêa (2002), who found increased aggregates and weighted average diameter innative forest, with reduction as intensive land use practices were carried out. Different agricultural management practices interfered in aggregates formation because it modifies soil organic matter dynamics (Zanatta et al., 2007) and conditions for microorganism activity (Vargas and Scholles, 2000).
Considering the principal component 1 (PC1), which showed 55.26% of total data variation, physical attributes that had the highest weights (in module) were: soil density (SD), total porosity (TP), field capacity (FC) and available water (AW). In CP2, which explained 41.46% of data explanation, the largest eigenvectors were: Aggregation (AGREG.), weighted average diameter (WAD) and particle density (PD) (Table 4).
In the graphic representation (Figure 1), the two first principal components scores (PC1 and PC2) were considered for groups’ interpretation. Two different groups were formed; the first refers to cassava (CASS) and beans (BEAN) monoculture, indicating that practices adopted in these areas have not contributed to physical attributes differentiation. Analyzing the second group, native forest (NF), it differed from the first group (CASS and BEAN), being higher regarding soil physical properties. Albuquerque et al. (2005), Souza et al. (2005), Carneiro et al. (2009), Corrêa et al. (2010) obtained similar results, as agricultural use has changed soil physical properties, when compared to the native vegetation area.
Assessing soil water retention curve under different agricultural uses (NF, CASS and BEAN), which was adjusted according to the mathematical model proposed by Van Genuchten (1980), (Figure 2) the different agricultural uses showed uniformity, except for cassava monoculture, at 0 to 100 kPa tensions, which were below the others, due to lower total porosity (Table 3).
Water retention process occurs in unsaturated soils as a result of capillarity. Adsorption forces operating in the soil matrix vary depending on the texture, mineralogy, depth, structure, porosity, organic matter and agricultural uses (Libardi, 2005; Reichardt and Timm, 2004).
Viscosity is a property that reflects the ease with which the molecules or particles slide over the other, being directly proportional to the volume of the particles varying with the temperature. It is affected by solutes type and concentration of solutes (Reichardt and Timm, 2004).
Native forest (NF) showed higher water retention, total porosity, field capacity and available water. This can be justified by OM increase due to vegetation cover maintenance, soil non-disturbance and iron aluminumoxides presence, which give higher aggregate stability. Results obtained by (Roth et al., 1991; Bertol et al., 2000; Giarola et al., 2002) corroborate with this study.
Soil and agricultural crops management modify the soil structure, and consequently, its physical properties, such as density, porosity (Faria et al., 1998; Pires et al., 2012) and soil water retention (Ramos et al., 2013). Therefore, management based on monoculture, compared to native forest, negatively influenced on physical and hydraulic attributes (Table 3 and Figure 2). In Brazil Portela et al. (2001), among others, Cintra et al. (2009), found modifications in soil pores quantity and quality. The soils of both studies were classified as Oxisol in water retention, in a tableland ecosystem soil under different land use.
Different agricultural uses influence soil physical attributes mean values. The most significant physical attributes to distinguish agricultural uses were: aggregation; weighted average diameter; soil density; particle density; total porosity; field capacity and available water. Native forest obtained better results, followed by cassava and beans monoculture. Soil water characteristic curves showed uniformity among different agricultural uses with higher water retention in the native forest.
The authors have not declared any conflict of interests.
REFERENCES
|
Albuquerque JA, Argenton J, Bayer C, Wildner L do P, Kuntze MAG (2005). Relação de atributos do solo com a agregação de um Latossolo Vermelho sob sistemas de preparo e plantas de verão para cobertura do solo. Rev. Bras. Ciênc. Solo 29:415-424.
Crossref
|
|
|
|
Bertol I, Cogo NP, Levien R (1997). Erosão hídrica em diferentes preparos do solo após as colheitas de milho e trigo, na presença e ausência dos resíduos culturais. Rev. Bras. Ciênc. do Solo 21:409-418.
Crossref
|
|
|
|
|
Borges EM Lombardi Neto F, Correa GF, Costa LM (1997). Gesso e matéria seca vegetal na floculação de argila e na produção de soja em um latossolo vermelho-escuro com camada subsuperficial compactação. Rev. Bras. Ciênc. Solo 21:125-130.
|
|
|
|
|
Cardoso EL (2002). Solos do Assentamento Urucum - Corumbá, MS: caracterização, limitações e aptidão agrícola. Corumbá: Embrapa Pantanal. Documentos, 30. P. 35.
|
|
|
|
|
Carneiro MAC, Souza ED, Reis EF, Pereira HS, Azevedo WR (2009). Atributos físicos, químicos e biológicos de solo de cerrado sob diferentes sistemas de uso e manejo. Rev. Bras. Ciênc. do Solo 33:147-157.
Crossref
|
|
|
|
|
Castro Filho C, Muzilli O, Podanoschi AL (1998). Estabilidade dos agregados e sua relação com o teor de carbono orgânico em um Latossolo Roxo Distrófico, em função de sistemas de plantio, rotações de culturas e métodos de preparo das amostras. Rev. Bras. Ciênc. do Solo 22:527-538.
Crossref
|
|
|
|
|
Cintra FLD, Resende RS, Silva MLL, Portela JC (2009). Efeito de volumes de água de irrigação no regime hídrico de solo coeso dos tabuleiros e na produção de coqueiro. Rev. Bras. Ciênc. Solo 33:1041-1051.
Crossref
|
|
|
|
|
Cordeiro FC, Dias F de C, Merlim A de O, Correia MEF, Aquino AM, Brown G (2004). Diversidade da macrofauna invertebrada do solo como indicadora da qualidade do solo em sistema de manejo orgânico de produção. Revista Universidade Rural: Série Ciênc. Vida 24:29-34.
|
|
|
|
|
Corrêa JC (2002). Efeito de sistemas de cultivo na estabilidade de agregados de um Latossolo Vermelho-Amarelo em Querência, MT. Pesquisa Agropecuária Bras. 37:203-209.
Crossref
|
|
|
|
|
Corrêa M, Freire MBG, Ferreira RLC, Silva JAA, Pessoa LGM, Miranda MA, Melo DVM (2010). Atributos físicos de solos sob diferentes usos com irrigação no semiárido de Pernambuco. Rev. Bras. Eng. Agríc. Ambiental 14:358-365.
|
|
|
|
|
Cunha EQ, Stone LF, Moreira JAA, Ferreira EPB, Didonet AD, Leandro WM (2011). Sistemas de preparo do solo e culturas de cobertura na produção orgânica de feijão e milho. I - Atributos físicos do solo. Rev. Bras. Ciênc. Solo 35:589-602.
Crossref
|
|
|
|
|
Dourado Neto, D, Van Lier QJ, Botrel TA, Libardi PL (1990). Programa para confecção da curva de retenção de água no solo utilizando o modelo de van Genuchten. Eng. Rural 1:92-102.
|
|
|
|
|
Donagema GK, Campos DVB, Calderano SB, Teixeira WG, Viana JHM (2011). (Org.). Manual de métodos de análise de solos. 2.ed. rev. Rio de Janeiro: Embrapa Solos. P. 230. Documentos 132.
|
|
|
|
|
Faria JC, Schaefer CER, Ruiz HA, Costa LM (1998). Effects of weed control on physical and micropedological properties of Brazilian Ultisol. Rev. Bras. Ciênc. Solo 22:731-741.
Crossref
|
|
|
|
|
Ferreira RRM, Tavares Filho J, Ferreira VM, Ralisch R (2010). Estabilidade física de solo sob diferentes manejos de pastagem extensiva em Cambissolo. Semina: Ciênc. Agrárias 31:531-538.
|
|
|
|
|
Giarola NFB, Silva AP, Imhoff S (2002). Relações entre propriedades físicas e características de solos da região Sul do Brasil. Rev. Bras. Ciênc. Solo 26:885-893.
Crossref
|
|
|
|
|
Goedert WJ, Oliveira AS (2007). Fertilidade do solo e sustentabilidade da atividade agrícola. In: Novais RF, Alvarez VVH, Barros NF, Fontes
|
|
|
|
|
RLF, Cantarutti RB, Neves JCL (2007). Fertilidade do solo. Viçosa: Soc. Bras. Ciênc. Solo Cap. 18:991-1017.
|
|
|
|
|
Kämpf N, Curi N (2012). Conceito de solo e sua evolução histórica. In: Ker JC, Curi N, Schaefer CEGR, Torrado PV (2012). Pedologia: Fundamentos. Viçosa: Soc. Bras. Ciênc. Solo. Cap. 1:1-20.
|
|
|
|
|
Kay BD (1990). Rates of change of soil structure under different cropping systems. Adv. Soil Sci. 12:1-41.
Crossref
|
|
|
|
|
LIBARDI PL (2005). Dinâmica da água no solo. São Paulo: Editora da Universidade de São Paulo (EDUSP). P. 344.
|
|
|
|
|
Pimentel GF (2009). Curso de estatística experimental. 15. Ed. Piracicaba: FEALQ. P 451.
|
|
|
|
|
Pires BS, Dias Junior MS, Rocha WW, Araújo Júnior CF, Carvalho RCR (2012). Modelos de capacidade de suporte de carga de um Latossolo Vermelho-amarelo sob diferentes usos e manejos. Rev. Bras. Ciênc. Solo 36:635-642.
Crossref
|
|
|
|
|
Portela JC, Libardi PL, van Lier Q de J (2001). Retenção da água em solo sob diferentes usos no ecossistema tabuleiros costeiros. Rev. Bras. Eng. Agríc. Ambient. Campina Grande 5(1):49-54.
|
|
|
|
|
Ramos BZ, Pai PSM, Freitas WA, Dias Junior M de S (2013). Avaliação dos atributos físico-hídricos em um Latossolo Vermelho distroférrico sob diferentes sistemas de manejo-Lavras/Minas Gerais/Brasil. Rev. Ciênc. Agrárias 36:340-346.
|
|
|
|
|
Reichardt K, Timm LC (2004). Solo Planta e atmosfera: conceitos, processos e aplicações. São Paulo. P 478.
|
|
|
|
|
Reichert JM, Kaiser D R, Reinert DJ, Riquelme UFB (2009). Variação temporal de propriedades físicas do solo e crescimento radicular de feijoeiro em quatro sistemas de manejo. Pesquisa Agropecuária Bras. 44:310-319.
Crossref
|
|
|
|
|
Roth CH, Castro Filho C, Medeiros GB (1991). Análise de fatores físicos e químicos relacionados com a agregação de um Latossolo roxo distrófico. Rev. Bras. Ciênc. Solo 15:241-248.
|
|
|
|
|
Santos HG, Jacomine PKT, Anjos LHC, Oliveira VA, Oliveira JB, Coelho MR, Lumbreras JF, Cunha TJF (2013). Sistema Brasileiro de Classificação de Solos. 3. Ed. Brasília: Embrapa P. 353.
|
|
|
|
|
Silva MB, Kliemann HJ, Silveira PM, Lanna AC (2007). Atributos biológicos do solo sob influência da cobertura vegetal e do sistema de manejo. Pesquisa Agropec. Bras. 42:1755-1761.
Crossref
|
|
|
|
|
Souza ED, Carneiro MAC, Paulino HB (2005). Atributos físicos de um Neossolo Quartzarênico e um Latossolo Vermelho sob diferentes sistemas de manejo. Pesquisa Agropecuária Bras. 40:1135-1139.
Crossref
|
|
|
|
|
Van Genuchten MT (1980). A Closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44:892-898.
Crossref
|
|
|
|
|
Vargas LK, Scholles D (2000). Biomassa microbiana e produção de C-CO2 e N mineral de um Podzólico Vermelho-Escuro submetido a diferentes sistemas de manejo. Rev. Bras. Ciênc. Solo 24:35-42.
Crossref
|
|
|
|
|
Zanatta JA, Bayer C, Dieckow J, Vieira FCB, Mielniczuk J (2007). Soil organic carbon accumulation and carbon costs related to tillage, cropping systems and nitrogen fertilization in a subtropical Acrisol. Soil Till. Res. 94:510-519.
Crossref
|
|
.png)
.png)
