RESUMO
Compreender a variabilidade espacial dos atributos químicos e físicos do solo é importante para melhorar as práticas de manejo e conservação do solo. Por sua vez, a variabilidade espacial das propriedades do solo resulta da variação das características morfológicas do relevo. Assim, o objetivo deste estudo foi avaliar a influência da curvatura da superfície na variabilidade espacial dos atributos químicos e físicos do solo na região de Mar de Morros (Pinheiral-Rio de Janeiro, Brasil). Foram selecionadas duas pedoformas adjacentes com curvatura convexa e côncava e coletadas amostras de solo em uma grade com espaçamento fixo de 10 metros entre os pontos. Foram coletadas um total de 56 amostras de solo (0-5 cm de profundidade) nas duas pedoformas e analisados os atributos físicos e químicos do solo. Na análise dos dados foi utilizada estatística descritiva e geoestatísitca. Os atributos químicos do solo apresentaram um padrão aleatório quanto a distribuição espacial em ambas as pedoformas. No entanto, a pedoforma côncava apresentou maiores valores de pH e potássio e valores mais baixos para alumínio quando comparada a pedoforma convexa. Por outro lado, os valores de silte apresentaram dependência espacial em ambas às pedoformas. Já a densidade do solo e o teor de argila apresentaram dependência espacial na pedoforma convexa. Os valores da densidade do solo e silte aumentaram no sentido do terço superior para o terço inferior em ambas as pedoformas. Os resultados mostram que para esta área de estudo, a curvatura da superfície tem maior influência na dependência espacial dos atributos físicos do que para os atributos químicos solo.
Palavras chave:
Relevo; Geoestatística; Fertilidade; Floresta Atlântica
ABSTRACT
Understanding the spatial variability of soil chemical and physical attributes is important for improving management practices and soil conservation. In turn, the spatial variability of soil properties results from variation in morphological relief characteristics. Thus, the purpose of this study was to evaluate the influence of landform curvature on the spatial variability of soil chemical and physical attributes in the Mar de Morros region (Pinheiral-Rio de Janeiro State, Brazil). Two adjacent landforms were selected with convex and concave curvature and sampled in a regularly spaced grid of 10 meters. A total of 56 soil samples (0-5 cm depth) from the two landforms were collected and analyzed for physical and chemical attributes. The data were analyzed using descriptive statistics and geostatistics. All chemical attributes showed random patterns of spatial variability in both landforms. The concave landform had higher values of pH and potassium and lower values of aluminum than the convex landform. In contrast, silt content showed spatial dependence in both the concave and convex landforms. Bulk density and clay showed spatial dependence in the convex landform. Bulk density and silt content increased from the shoulder to the footslope of both landforms. The results show that, for this study area, landform curvature has more influence on the spatial dependence of soil physical attributes than of soil chemical properties.
Keywords:
Relief; Geostatistics; Fertility; Atlantic Forest
INTRODUCTION
In the Mar de Morros region (Pinheiral-Rio de Janeiro State, Brazil), the Atlantic Forest is represented by Seasonal Semideciduous Forest, which historically has been an important area for growing coffee and raising livestock. Areas previously occupied by Atlantic Forest are now predominantly pastures with different degrees of degradation. Some are still used as unmanaged pasture; however, others were abandoned and have reverted to native vegetation in different states of natural succession (MENEZES, 2008MENEZES, C. E. G. Integridade de paisagem, manejo e atributos do solo no Médio Vale do Paraíba do Sul, Pinheiral-RJ. 2008. 175p. Tese de Doutorado. Universidade Federal Rural do Rio de Janeiro, Seropédica, 2008.).
The rate of natural vegetation succession in anthropized abandoned areas is determined by the physical and chemical characteristics of the soil and is fastest in the most fertile environments.
The rate of natural vegetation succession in disturbed abandoned areas is determined by factors, such as seed bank, seedling source, degree of human disturbance, and soil physical and chemical properties. Natural succession is fastest in the most fertile environments (PINTO et al.; 2008PINTO, S. I. C.; MARTINS, S. V.; BARROS, N. F.; DIAS, H. C. T.; KUNZ, S. H. Influence of environmental variables on the shrub and tree species distribuition in two Semideciduous Forest sites in Viçosa, Minas Gerais, Brazil. Revista de Biologia Tropical, v. 56, n. 3, p. 1557-1569, 2008.). As observed by Souza et al. (2012SOUZA, P. B.; LELIS, J. J.; SCHAEFER, C. E. G. R.; SOUZA, A. L.; MEIRA NETO, J. A. A. Distribution of tree species in a geomorphological and pedological gradient of Sub montane Semidecidual Seasonal Forest in the vicinty of Rio Doce state park, Minas Gerais. Revista Árvore , v. 30, n. 4, p. 707-718, 2012.) in Rio Doce State Park in Minas Gerais, the distribution of tree species of Semideciduous Sub-montane Forest was affected by the gradient of soil physical and chemical attributes. The distribution of tree species in the forest was also related to soil drainage regimes (FERREIRA-JUNIOR et al., 2007FERREIRA-JÚNIOR, W. G.; SILVA, A. F.; SCHAEFER, C. E. G. R.; MEIRA NETO, J. A. A.; DIAS, A. S.; IGNÁCIO, M.; MEDEIROS, M. C. M. P. Influence of soils and topographic gradientes on tree species distribution in a brazilian Atlantic Tropical Semideciduos Forest. Edinburgh Journal of Botany, v. 64, n. 2, p. 137-157, 2007.).
The spatial variability of soil physical and chemical attributes is determined by the source material, climate and relief characteristics of the landscape (CEDDIA et al., 2009CEDDIA, M. B.; VIEIRA, S. R.; VILLELA, A. L. O.; MOTA, L. S.; ANJOS, L. H. C.; CARVALHO, D. F. Topography and spatial variability of soil physical properties. Scientia Agricola, v. 66, n. 3, p. 338-352, 2009.; CAMARGO et al., 2010CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Spatial variability of physical atributes of an Alfisol under diferent hillslope curvatures. Revista Brasileira de Ciência do Solo, v. 34, n. 3, p. 617-30, 2010.). Relief forms in these environments vary from convex, through linear, to concave (TROEH, 1965TROEH, F. R. Landform equations fitted to contour maps. American Journal of Science, 263:616-27, 1965.).
Each landform has variability in surface curvature, declivity and altitude, which influence the dynamics (accumulation and dispersion) of water and sediment in the landscape (SANCHEZ et al., 2009SANCHEZ, R. B.; MARQUES JÚNIOR, J.; SOUZA, Z. M.; PEREIRA, G. T.; MARTINS FILHO, M. V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, v. 68, n. 4, p. 1095-1103, 2009.; CAMARGO et al., 2010CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Spatial variability of physical atributes of an Alfisol under diferent hillslope curvatures. Revista Brasileira de Ciência do Solo, v. 34, n. 3, p. 617-30, 2010.). Convex, linear and concave surfaces influence the source material of sediment and the intensity and direction of water flow within the soil profile (ARTUR et al., 2014ARTUR, A. G.; OLIVEIRA, D. P.; COSTA, M. C. G.; ROMERO, R. E.; SILVA, M. V. C.; FERREIRA, T. O. Variabilidade espacial dos atributos químicos do solo, associada ao microrrelevo. Revista Brasileira Engenharia Agrícola e Ambiental, v. 18, n. 2, p.141-149, 2014.). According to Gandolfi (2000GANDOLFI, S. História natural de uma Floresta Estacional Semidecidual no município de Campinas. 2000. 520p. Tese de Doutorado. Universidade Estadual de Campinas, Campinas.), the variation in the surface slope conditions also influence soil humidity and fertility gradients.
Some studies have observed that concave environments show greater spatial variability in soil attributes than do linear environments (SOUZA et al., 2003SOUZA, C. K,; MARQUES JÚNIOR, J.; MARTINS FILHO, M. V.; PEREIRA, G. T. Influência do relevo e erosão na variabilidade espacial de um Latossolo em Jaboticabal (SP). Revista Brasileira de Ciência do Solo , v. 27, n. 1, p. 1067-74, 2003.; MONTANARI et al., 2005MONTANARI, R.; MARQUES JÚNIOR.; PEREIRA, G. T.; SOUZA, Z. M. Forma da paisagem como critério para otimização amostral de Latossolos sob cultivo de cana-de-açúcar. Pesquisa Agropecuária Brasileira, v. 40, n. 1, p. 69-77, 2005.; CAMARGO et al., 2010CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Spatial variability of physical atributes of an Alfisol under diferent hillslope curvatures. Revista Brasileira de Ciência do Solo, v. 34, n. 3, p. 617-30, 2010.). According to Sanchez et al. (2009SANCHEZ, R. B.; MARQUES JÚNIOR, J.; SOUZA, Z. M.; PEREIRA, G. T.; MARTINS FILHO, M. V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, v. 68, n. 4, p. 1095-1103, 2009.), greater variability in the concave environments is linked directly to variability in relief characteristics.
The pattern of soil attributes in different landscape conditions can be understood by characterizing the spatial variability of the soil physical and chemical attributes based on the surface geomorphology. According to Sibaldelli et al. (2015SIBALDELLI, R. N. R.; CARVALHO, J. D. F. C.; DE OLIVEIRA, M. C. N. Uso de geoestatística no estudo da variabilidade espacial da capacidade de troca de cátions do solo. Global Science and Technology, v. 8, n. 1, p. 141-156, 2015.), understanding the spatial and temporal variability of soil attributes is required in studies of agricultural soils and their productive capacities. In non-production forests, understanding the spatial variability of soil properties is important for improving management techniques aimed at recovery from anthropogenic disturbances.
Geostatistics have been used by many researchers to evaluate physical and chemical attributes of the soil (AQUINO et al., 2015AQUINO, R. E.; COSTA CAMPOS, M. C.; OLIVEIRA, I. A.; SIQUEIRA, D. S.; RODRIGUES SOARES, M. D.; FREITAS, L. Técnicas geoestatísticas na avaliação de atributos químicos em Cambissolo com agrofloresta e cana-de-açúcar em Humaitá, Amazonas. Revista Brasileira de Ciências Agrárias, v. 10, n. 4, p. 544-552, 2015.; JUNIOR et al., 2014JUNIOR O. , O. S.; SANTOS, A. C.; SANTOS, P. M.; ROCHA, J. M. L.; ALEXANDRINO, E. Distribuição espacial de atributos químicos do solo em áreas sob diferentes usos agrícolas. Revista Engenharia na Agricultura, v. 22, n. 5, p. 442-455, 2014.; OLIVEIRA et al., 2013OLIVEIRA, I. R. D.; TEIXEIRA, D. D. B.; PANOSSO, A. R.; CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Modelagem geoestatística das incertezas da distribuição espacial do fósforo disponível no solo, em área de cana-de-açúcar. Revista Brasileira de Ciência do Solo , v. 37, n. 2, p. 1481-1491, 2013.). Thus, the objective of this study is to evaluate the influence of landform curvature on the spatial variability of physical and chemical soil attributes in the Mar de Morros region, in the municipality of Pinheiral, Rio de Janeiro State, Brazil.
MATERIAL AND METHODS
The study was conducted in the municipality of Pinheiral, Rio de Janeiro State, in the Middle Paraíba Region. The study site is located in the Mar de Morros geomorphologic region, in the sub-basin of the Ribeirão Cachimbal, which is in the hydrographic basin of the Paraíba do Sul river. Detailed information about the study area and its location is presented in Figure I.
Geographical location of the study area. In blue the municipality of Pinheiral, Rio de Janeiro State.
Two adjacent landforms were selected, one with convex curvature (convex-divergent) and concave curvature (concave-convergent). Both landforms had a predominantly southeastern aspect and current Semideciduous Sub-montane plant cover (CONAMA, 1994BRASIL. Ministério do Meio Ambiente, Resolução Conama 006 de 4 de maio de 1994. Brasília, DF: Ministério do Meio Ambiente. Disponível em: <Disponível em: http://www.mma.gov.br/conama/legiano1.cfm/
>. Acessado em 20 de maio de 2012.
http://www.mma.gov.br/conama/legiano1.cf...
; IBGE, 2012INSTITUTO BRASILEIRO DE GEOGRAFIA E ESTATÍSTICA. Manual Técnico da Vegetação Brasileira: Sistema fitogeográfico inventário das formações florestais e campestres, técnicas e manejo de coleções botânicas, procedimentos para mapeamentos. 2ª ed. Rio de Janeiro. 2012, 271p.). This vegetation has established itself in the site through natural regeneration and has been present in this environment for more than sixty years.
The soils of both landforms were classified as Inceptsol (Cambissolo Háplico). In each landform, soil samples were collected from 0-5 cm depth. We focused on the top 5 cm because this layer is the most sensitive to surface runoff processes, erosion and drainage. In the convex landform, 24 samples were collected in a grid pattern with 10-meter spacing. In the concave landform, 32 samples were collected, and grid spacing was reduced to 8 meters in the central section (Figure 2). This procedure aimed to avoid sampling areas that had surface water flow during rainy periods.
Morphological characteristics (curvature and slope), landscape position and the regular grid of soil sampling points from convex and concave landforms.
A total of 56 soil samples were collected from the two landforms. Bulk density (BD) was measured on undisturbed soil samples collected using an Uhland sampler. Soil texture (clay, silt and sand) and chemical attributes (pH, Ca2+, Mg2+, K+, Al3+, P and C) were measured on deformed samples (EMBRAPA, 1997EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA. Manual de métodos de análise de solo. Rio de Janeiro, EMBRAPA/SNLCS. 1997, 212p.).
Statistical analyses were performed using two approaches. The first analysis consisted of descriptive statistics and geostatistics. This analysis was performed to evaluate whether the soil attributes showed spatial dependence, using GEOSTAT software (VIEIRA et al., 1983). For soil attributes that did not show spatial dependence (i.e., had a completely random pattern), the average value of each attribute was compared between landforms using a Tukey test (α=0.05).
The steps of geostatistical analysis were as follows: 1- descriptive statistical analysis; 2 - determination of spatial dependence (variography), 3 - modeling and validation of experimental semivariograms and interpolation. All steps were performed using the GEOSTAT software (VIEIRA et al, 1983). The semivariances of soil attributes were calculated according to the equation below, where: g(h) - semi variance estimated at a distance h; N(h) - number of pairs of values Z (x i), Z (x i + h)] separated by a vector h; X i - spatial position of variable Z, and Z - values of variables (physical and chemical attributes of soil).
The experimental semivariogram was also classified according to the degree of randomness (DR). The DR is calculated through the following equation 2, where: C0 = intercept or nugget effect and C0 + C1 = Sill.
Based on the methodology proposed by Cambardella et al. (1994CAMBARDELLA, C. A.; MOORMAN, T. B.; NOVAK, J. M.; PARKIN, T. B.; KARLEN, D. L.; TURCO, R. F.; KONOPKA, A. E. Field scale variability of soil properties in central Iowa soils. Soil Science Society of America Journal, v. 58, n. 5, p.1501-1511, 1994.), the following classes were used: Low Degree of Randomness (DR ≤ 25%), Moderate DR (DR between 25% and 75%), High DR (DR ≥ 75%), and finally Complete Randomness (DR = 100%). In the last class, the semivariogram is classified as pure nugget effect.
The interpolation of the soil attributes was performed using Ordinary Kriging (OK). The OK only utilizes primary data, such as clay content measured at sampled locations u to estimate the clay content at unsampled locations. The stationarity of the mean is assumed only within a local neighborhood W(u), centered at the location u being estimated. Here, the mean is deemed to be a constant but unknown value, i.e., m(u´) = constant but unknown, “u´Î W(u). The OK estimator is written as a linear combination of the n(u) data Z(uα) with a single unbiasedness constraint:
withThe unknown local mean m(u) is filtered from the linear estimator by forcing the kriging weights to sum to 1.
RESULTS AND DISCUSSION
Soil physical and chemical attributes for the convex and concave landforms are presented in Table I. Mean soil bulk density was not different between landforms. In contrast, clay and sand contents tended to be higher in the convex landform, while the highest value of silt content was observed in the concave landform. These differences in soil texture are likely due to different patterns of water flow in the two landforms, which in the convex landform are dominated by the process of side outflow and the carrying effect of soil particles (OLIVEIRA et al., 2013OLIVEIRA, D. P.; FERREIRA, T. O.; ROMERO, R, E.; FARIAS, P. R. S.; COSTA, M. C. G. Microrrelevo e a distribuição de frações granulométricas em Cambissolos de origem calcária. Revista Ciência Agronômica , v. 44, n. 4, p. 676-684, 2013.).
Soil pH, Ca+2, Mg+2, P and K+ were higher in the concave landform, while Al+2 and C contents were higher in the convex landform. According to Souza et al. (2004SOUZA, Z. M.; MARQUES JÚNIOR, J.; PEREIRA, G.T. Influência da pedoforma na variabilidade espacial de alguns atributos físicos e hídricos de um Latossolo sob o cultivo da cana-de-açúcar. Irriga, Piracicaba, v. 9, n.1, p.1-11, 2004.), small variations in landform relief are responsible for differences in the variability of chemical attributes in the landscape. Differences in soil chemical properties are also associated with variation in organic matter accumulation processes, nutrient transport and sediment and water dynamics (CAMPOS et al., 2010CAMPOS, M. C. C.; RIBEIRO, M. R.; SOUZA JUNIOR, V. S.; RIBEIRA FILHO, M. R.; OLIVEIRA, I. A. Interferência dos pedoambientes nos atributos do solo em uma topossequência da transição Campos/Floresta. Revista Ciência Agronômica, Fortaleza, v. 41, n. 4, p. 527-535, 2010.).
For soil chemical attributes, pH and potassium mean levels were statistically higher in the concave landform, while aluminum and carbon mean levels were statistically higher in the convex landform (Table I). The influence of relief characteristics on nutrient availability in different landforms was also observed by Montanari et al. (2010MONTANARI, R.; MARQUES JÚNIOR, J.; CAMPOS, M. C. C.; SOUZA, Z. M.; CAMARGO, Z. A. Caracterização mineralógica de Latossolos em diferentes feições do relevo na região de Jaboticabal, SP. Revista Ciência Agronômica , v. 41, n. 2, p. 191-199, 2010.), which associated the variability in chemical attributes with different drainage conditions. The drainage influenced the redistribution of nutrients in the soil layers. Studies have highlighted the role and influence of water flow in the distribution of soil nutrients. Regions with water accumulation tend to have higher soil fertility and higher silt and sand contents (FERREIRA-JUNIOR et al., 2007FERREIRA-JÚNIOR, W. G.; SILVA, A. F.; SCHAEFER, C. E. G. R.; MEIRA NETO, J. A. A.; DIAS, A. S.; IGNÁCIO, M.; MEDEIROS, M. C. M. P. Influence of soils and topographic gradientes on tree species distribution in a brazilian Atlantic Tropical Semideciduos Forest. Edinburgh Journal of Botany, v. 64, n. 2, p. 137-157, 2007.; SOUZA et al., 2012SOUZA, P. B.; LELIS, J. J.; SCHAEFER, C. E. G. R.; SOUZA, A. L.; MEIRA NETO, J. A. A. Distribution of tree species in a geomorphological and pedological gradient of Sub montane Semidecidual Seasonal Forest in the vicinty of Rio Doce state park, Minas Gerais. Revista Árvore , v. 30, n. 4, p. 707-718, 2012.).
The variation in carbon content between the landforms is related to differences in litter decomposition rates. According to Santos (2014SANTOS, G. L. Efeito da Pedoforma no Processo de Sucessão Secundária em Fragmentos Florestais na Região do Médio Vale do Paraíba do Sul, Pinheiral, RJ. 2014. 130p. Dissertação de Mestrado. Universidade Federal Rural do Rio de Janeiro, Seropédica, 2014.), higher decomposition rates were observed in convex landforms than in concave ones. This pattern was associated with variations in environmental conditions.
Studies developed by Souza et al. (2003SOUZA, C. K,; MARQUES JÚNIOR, J.; MARTINS FILHO, M. V.; PEREIRA, G. T. Influência do relevo e erosão na variabilidade espacial de um Latossolo em Jaboticabal (SP). Revista Brasileira de Ciência do Solo , v. 27, n. 1, p. 1067-74, 2003.) and Montanari et al. (2005MONTANARI, R.; MARQUES JÚNIOR.; PEREIRA, G. T.; SOUZA, Z. M. Forma da paisagem como critério para otimização amostral de Latossolos sob cultivo de cana-de-açúcar. Pesquisa Agropecuária Brasileira, v. 40, n. 1, p. 69-77, 2005.) showed that, in comparison to linear landforms, concave landforms had greater variation in granulometric attributes, organic matter content, superficial flow and slope, indicating that soil attributes have different patterns as a function of landform curvature. According to Resende et al. (1997), concave landforms have areas with higher erosion rates in the highest landscape positions and higher sediment accumulation rates in the lowest landscape positions. This likely explains why concave landforms generally show more heterogeneous distributions of soil attributes.
All soil physical attributes in the concave landform followed a normal distribution. In the convex landform, however, clay and sand contents did not follow a normal distribution (Table I). According to Sirtoli et al. (2008SIRTOLI, A. E.; SILVEIRA, C. T.; SILVA, C. R.; MONTOVANI, L. E.; RIBEIRO, S. R. A.; OKA-FIORI, C. Atributos topográficos secundários no mapeamento de pedoformas. Geociências, v. 27, n. 1, p. 763-78, 2008.), in concave environments the processes of sediment loss and accumulation occur simultaneously, which leads to soil physical attributes following a normal distribution., In convex environments, in contrast, sediment loss is the dominant process and causes a non-linear distribution.
Soil carbon and phosphorus also followed a normal distribution in both landforms (Table I), which can be related to the distinctive processes of sediment loss and accumulation for phosphorus and to the contribution of decomposing leaf litter for carbon. According to Santos (2014SANTOS, G. L. Efeito da Pedoforma no Processo de Sucessão Secundária em Fragmentos Florestais na Região do Médio Vale do Paraíba do Sul, Pinheiral, RJ. 2014. 130p. Dissertação de Mestrado. Universidade Federal Rural do Rio de Janeiro, Seropédica, 2014.), which evaluated the effect of landform on secondary succession in forest fragments, the distribution of soil chemical attributes in forest environments was influenced by relief characteristics and litter inputs, both of which drive organic-matter dynamics and nutrient cycling in the soil.
Potassium levels in the convex landform and magnesium levels in the concave landform followed normal distributions, while pH, calcium and aluminum did not follow normal distributions in either landform (Table I). According to Artur et al. (2014ARTUR, A. G.; OLIVEIRA, D. P.; COSTA, M. C. G.; ROMERO, R. E.; SILVA, M. V. C.; FERREIRA, T. O. Variabilidade espacial dos atributos químicos do solo, associada ao microrrelevo. Revista Brasileira Engenharia Agrícola e Ambiental, v. 18, n. 2, p.141-149, 2014.), higher spatial variability in soil chemical attributes results from microrelief characteristics, which influence the direction of soil water flow.
Comparing statistical parameters, the highest values for variance, standard deviation and coefficient of variation were observed for the concave landform; this pattern can be associated with variation in altitude, slope and landscape position (Figure 3). According to Campos et al. (2010CAMPOS, M. C. C.; RIBEIRO, M. R.; SOUZA JUNIOR, V. S.; RIBEIRA FILHO, M. R.; OLIVEIRA, I. A. Interferência dos pedoambientes nos atributos do solo em uma topossequência da transição Campos/Floresta. Revista Ciência Agronômica, Fortaleza, v. 41, n. 4, p. 527-535, 2010.), the distribution of soil physical and chemical attributes are dependent on the conditions of the landforms.
Semivariograms of experimental and respective models of soil physical attributes in convex and concave landforms (0-5 cm depth).
The results of the geostatistical analysis of soil physical and chemical attributes are presented in Table 2 and Figure 3. None of the chemical attributes showed spatial dependence; they had DR values of 100% indicating pure nugget effect. The geostatistical analysis demonstrated that the chemical attributes showed a random pattern in both the convex and the concave landforms.
Consequently, the average values of chemical attributes in each landform could be statistically compared using a Tukey test, since the errors follow a random pattern. Another aspect of the pure nugget effect for chemical attributes is that ordinary kriging is not the best interpolation method for generating maps of these attributes. Ordinary kriging only outperforms other interpolation methods when the attribute shows spatial dependence. Due to the observed lack of spatial dependence, maps of chemical attributes were not generated.
The variability in soil physical attributes and erosion factors in different landforms was evaluated by Sanchez et al. (2009SANCHEZ, R. B.; MARQUES JÚNIOR, J.; SOUZA, Z. M.; PEREIRA, G. T.; MARTINS FILHO, M. V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, v. 68, n. 4, p. 1095-1103, 2009.). The authors observed that the morphological characteristics (relief characteristics and slope) of the landform were responsible for different distribution patterns of soil attributes in the landscape. In addition, the landform factor exerted a greater influence on the distribution of soil attributes than did the erosion factor.
The random pattern of spatial variability observed for soil chemical attributes along the topographic gradient can be explained by the influence of nutrient cycling and soil cover (SANTOS, 2014SANTOS, G. L. Efeito da Pedoforma no Processo de Sucessão Secundária em Fragmentos Florestais na Região do Médio Vale do Paraíba do Sul, Pinheiral, RJ. 2014. 130p. Dissertação de Mestrado. Universidade Federal Rural do Rio de Janeiro, Seropédica, 2014.). The soil chemical attributes in the convex and concave environments can be related to sediment loss dynamics and the distribution of plant and animal material along the landforms. Rassol et al. (2014RASSOL, S. N.; GAIKWAD, S. W.; TALAT, M. A. Relationships between soil properties and slope segments of sallar wullarhama watershed in the Liddar Catchment of Jammu and Kashmir. Asian Journal of Engineering Research, v. 2, n. 2, p. 1-10, 2014.) evaluated the relationships between the soil properties and the slope at the boundary between the Jammu and Caxemira watersheds in Liddar, and observed that the slope and hillside positions significantly affected the circulation and accumulation of water, nutrients and sediments in the soil, which caused variability in soil properties along the landscape.
The lack of spatial dependence in soil chemical attributes may be related to the dynamics of variables, such as relief conditions, vegetation and soil in the different landforms. Some studies, such as Souza et al. (2008SOUZA, G. S.; LIMA, J. S. S.; SILVA, A. S.; OLIVEIRA, R. B. Variabilidade espacial de atributos químicos em um Argissolo sob Pastagem. Acta Scientiarum Agronomy, v. 30, n. 4, p. 589-596, 2008.) and Lima et al. (2010LIMA, J. S. S.; SOUZA, G. S.; SILVA, S. S. Amostragem e variabilidade espacial de atributos químicos do solo em área de vegetação natural em regeneração. Revista Árvore, Viçosa, v. 34, n. 1, p. 127-136, 2010.), evaluating the influence of variability in slope and topographic gradients on soil chemical attributes, observed the spatial dependence of soil attributes based on the type of landform. Their results are different from those observed in this study; it is possible to infer that this difference is related to the size of the study areas and the distribution and number of sampling points.
Soil physical attributes showed spatial dependence, though the details differed between landforms. In the convex landform, all soil physical attributes except sand content had experimental semivariograms that were best fitted by a Gaussian model (Table 2 and Figure 3). The range of spatial dependence varied from 77 to 96 meters for different physical attributes. In the concave landform, only silt content showed spatial dependence, and a Gaussian model was also fitted to the experimental semivariograms. The distribution pattern of soil physical attributes can be explained by variability in the topographical conditions that determine the behavior of water flow and sediments in the study region.
Studying the influence of microrelief on the spatial variability of different granulometric soil fractions in the Apodi Plateau, Oliveira et al. (2013OLIVEIRA, I. R. D.; TEIXEIRA, D. D. B.; PANOSSO, A. R.; CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Modelagem geoestatística das incertezas da distribuição espacial do fósforo disponível no solo, em área de cana-de-açúcar. Revista Brasileira de Ciência do Solo , v. 37, n. 2, p. 1481-1491, 2013.) observed a direct influence of microrelief on the distribution of sand, clay and silt contents, which was associated with the action of pedogenesis in soils due to convergent water flows.
The spatial variability maps of soil physical attributes in the convex landform (soil bulk density, clay and silt contents) and the concave landform (silt content) are presented in Figure 4. In all cases, a Gaussian model was fitted to the experimental semivariograms. In the convex landform, the degree of randomness was low for soil bulk density (DR = 0%) and moderate for clay and silt contents (DR = 36% and 37%, respectively). In the concave landform, the degree of randomness was moderate for silt contents (DR = 37%). Soil bulk density and silt content increased from the shoulder to the footslope of both landforms, while clay content showed the opposite pattern in the convex landform. The influence of relief on the spatial distribution of granulometric fractions in Cambisols was also observed by Oliveira et al. (2013OLIVEIRA, I. R. D.; TEIXEIRA, D. D. B.; PANOSSO, A. R.; CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Modelagem geoestatística das incertezas da distribuição espacial do fósforo disponível no solo, em área de cana-de-açúcar. Revista Brasileira de Ciência do Solo , v. 37, n. 2, p. 1481-1491, 2013.), where the highest clay and silt contents were on the concave surface.
Spatial distributions of soil physical attributes in convex and concave landforms (0-5 cm depth).
Variability in relief conditions seems to be the main cause of spatial variability in soil physical attributes in the study area (Figure 4). Similar results were observed by Burak and Passos (2011BURAK, D. L.; PASSOS, R. R. Variabilidade espacial das frações texturas e da porosidade do solo em relevo ondulado sob cultivo de café conilon. Enciclopédia Biosfera, v.7, n. 12, p. 01-11, 2011.), who studied the variability of textural fractions and soil porosity in wavy relief, where regions with lower superficial flow and slope had higher clay contents and lower macroporosity values. Understanding the spatial distribution of soil properties due to morphologic variation and their effects on forest ecological succession is important for developing strategies aimed at the recovery of anthropogenic environments.
CONCLUSIONS
The distribution patterns of the soil physical and chemical attributes were different between the two landforms evaluated, which reflected the different effects of slope, altitude and relief characteristics.
Landform surface curvature is a driving factor for the spatial variability of soil physical attributes, where soil bulk density and silt content increased from the shoulder to the footslope in concave and convex landforms. However, an opposite pattern was observed for clay content in the convex landform.
All soil chemical attributes showed a random pattern of spatial variability in both landforms. However, the highest values of soil chemical properties were observed in the concave landform.
REFERENCES
- AQUINO, R. E.; COSTA CAMPOS, M. C.; OLIVEIRA, I. A.; SIQUEIRA, D. S.; RODRIGUES SOARES, M. D.; FREITAS, L. Técnicas geoestatísticas na avaliação de atributos químicos em Cambissolo com agrofloresta e cana-de-açúcar em Humaitá, Amazonas. Revista Brasileira de Ciências Agrárias, v. 10, n. 4, p. 544-552, 2015.
- ARTUR, A. G.; OLIVEIRA, D. P.; COSTA, M. C. G.; ROMERO, R. E.; SILVA, M. V. C.; FERREIRA, T. O. Variabilidade espacial dos atributos químicos do solo, associada ao microrrelevo. Revista Brasileira Engenharia Agrícola e Ambiental, v. 18, n. 2, p.141-149, 2014.
- BURAK, D. L.; PASSOS, R. R. Variabilidade espacial das frações texturas e da porosidade do solo em relevo ondulado sob cultivo de café conilon. Enciclopédia Biosfera, v.7, n. 12, p. 01-11, 2011.
- CAMBARDELLA, C. A.; MOORMAN, T. B.; NOVAK, J. M.; PARKIN, T. B.; KARLEN, D. L.; TURCO, R. F.; KONOPKA, A. E. Field scale variability of soil properties in central Iowa soils. Soil Science Society of America Journal, v. 58, n. 5, p.1501-1511, 1994.
- CAMPOS, M. C. C.; RIBEIRO, M. R.; SOUZA JUNIOR, V. S.; RIBEIRA FILHO, M. R.; OLIVEIRA, I. A. Interferência dos pedoambientes nos atributos do solo em uma topossequência da transição Campos/Floresta. Revista Ciência Agronômica, Fortaleza, v. 41, n. 4, p. 527-535, 2010.
- CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Spatial variability of physical atributes of an Alfisol under diferent hillslope curvatures. Revista Brasileira de Ciência do Solo, v. 34, n. 3, p. 617-30, 2010.
- CEDDIA, M. B.; VIEIRA, S. R.; VILLELA, A. L. O.; MOTA, L. S.; ANJOS, L. H. C.; CARVALHO, D. F. Topography and spatial variability of soil physical properties. Scientia Agricola, v. 66, n. 3, p. 338-352, 2009.
- BRASIL. Ministério do Meio Ambiente, Resolução Conama 006 de 4 de maio de 1994. Brasília, DF: Ministério do Meio Ambiente. Disponível em: <Disponível em: http://www.mma.gov.br/conama/legiano1.cfm/ >. Acessado em 20 de maio de 2012.
» http://www.mma.gov.br/conama/legiano1.cfm/ - EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA. Manual de métodos de análise de solo. Rio de Janeiro, EMBRAPA/SNLCS. 1997, 212p.
- FERREIRA-JÚNIOR, W. G.; SILVA, A. F.; SCHAEFER, C. E. G. R.; MEIRA NETO, J. A. A.; DIAS, A. S.; IGNÁCIO, M.; MEDEIROS, M. C. M. P. Influence of soils and topographic gradientes on tree species distribution in a brazilian Atlantic Tropical Semideciduos Forest. Edinburgh Journal of Botany, v. 64, n. 2, p. 137-157, 2007.
- GANDOLFI, S. História natural de uma Floresta Estacional Semidecidual no município de Campinas. 2000. 520p. Tese de Doutorado. Universidade Estadual de Campinas, Campinas.
- INSTITUTO BRASILEIRO DE GEOGRAFIA E ESTATÍSTICA. Manual Técnico da Vegetação Brasileira: Sistema fitogeográfico inventário das formações florestais e campestres, técnicas e manejo de coleções botânicas, procedimentos para mapeamentos. 2ª ed. Rio de Janeiro. 2012, 271p.
- JUNIOR O. , O. S.; SANTOS, A. C.; SANTOS, P. M.; ROCHA, J. M. L.; ALEXANDRINO, E. Distribuição espacial de atributos químicos do solo em áreas sob diferentes usos agrícolas. Revista Engenharia na Agricultura, v. 22, n. 5, p. 442-455, 2014.
- LIMA, J. S. S.; SOUZA, G. S.; SILVA, S. S. Amostragem e variabilidade espacial de atributos químicos do solo em área de vegetação natural em regeneração. Revista Árvore, Viçosa, v. 34, n. 1, p. 127-136, 2010.
- MENEZES, C. E. G. Integridade de paisagem, manejo e atributos do solo no Médio Vale do Paraíba do Sul, Pinheiral-RJ. 2008. 175p. Tese de Doutorado. Universidade Federal Rural do Rio de Janeiro, Seropédica, 2008.
- MONTANARI, R.; MARQUES JÚNIOR, J.; CAMPOS, M. C. C.; SOUZA, Z. M.; CAMARGO, Z. A. Caracterização mineralógica de Latossolos em diferentes feições do relevo na região de Jaboticabal, SP. Revista Ciência Agronômica , v. 41, n. 2, p. 191-199, 2010.
- MONTANARI, R.; MARQUES JÚNIOR.; PEREIRA, G. T.; SOUZA, Z. M. Forma da paisagem como critério para otimização amostral de Latossolos sob cultivo de cana-de-açúcar. Pesquisa Agropecuária Brasileira, v. 40, n. 1, p. 69-77, 2005.
- OLIVEIRA, D. P.; FERREIRA, T. O.; ROMERO, R, E.; FARIAS, P. R. S.; COSTA, M. C. G. Microrrelevo e a distribuição de frações granulométricas em Cambissolos de origem calcária. Revista Ciência Agronômica , v. 44, n. 4, p. 676-684, 2013.
- OLIVEIRA, I. R. D.; TEIXEIRA, D. D. B.; PANOSSO, A. R.; CAMARGO, L. A.; MARQUES JÚNIOR, J.; PEREIRA, G. T. Modelagem geoestatística das incertezas da distribuição espacial do fósforo disponível no solo, em área de cana-de-açúcar. Revista Brasileira de Ciência do Solo , v. 37, n. 2, p. 1481-1491, 2013.
- PINTO, S. I. C.; MARTINS, S. V.; BARROS, N. F.; DIAS, H. C. T.; KUNZ, S. H. Influence of environmental variables on the shrub and tree species distribuition in two Semideciduous Forest sites in Viçosa, Minas Gerais, Brazil. Revista de Biologia Tropical, v. 56, n. 3, p. 1557-1569, 2008.
- RASSOL, S. N.; GAIKWAD, S. W.; TALAT, M. A. Relationships between soil properties and slope segments of sallar wullarhama watershed in the Liddar Catchment of Jammu and Kashmir. Asian Journal of Engineering Research, v. 2, n. 2, p. 1-10, 2014.
- RESENDE, M.; CURI, N.; RESENDE, S.B.; CORRÊA, G.F.; KER, J.C. Pedologia: Base para distinção de ambientes. 6.ed. Lavras: UFLA, 2014. 378p.
- SANCHEZ, R. B.; MARQUES JÚNIOR, J.; SOUZA, Z. M.; PEREIRA, G. T.; MARTINS FILHO, M. V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, v. 68, n. 4, p. 1095-1103, 2009.
- SANTOS, G. L. Efeito da Pedoforma no Processo de Sucessão Secundária em Fragmentos Florestais na Região do Médio Vale do Paraíba do Sul, Pinheiral, RJ. 2014. 130p. Dissertação de Mestrado. Universidade Federal Rural do Rio de Janeiro, Seropédica, 2014.
- SIBALDELLI, R. N. R.; CARVALHO, J. D. F. C.; DE OLIVEIRA, M. C. N. Uso de geoestatística no estudo da variabilidade espacial da capacidade de troca de cátions do solo. Global Science and Technology, v. 8, n. 1, p. 141-156, 2015.
- SIRTOLI, A. E.; SILVEIRA, C. T.; SILVA, C. R.; MONTOVANI, L. E.; RIBEIRO, S. R. A.; OKA-FIORI, C. Atributos topográficos secundários no mapeamento de pedoformas. Geociências, v. 27, n. 1, p. 763-78, 2008.
- SOUZA, P. B.; LELIS, J. J.; SCHAEFER, C. E. G. R.; SOUZA, A. L.; MEIRA NETO, J. A. A. Distribution of tree species in a geomorphological and pedological gradient of Sub montane Semidecidual Seasonal Forest in the vicinty of Rio Doce state park, Minas Gerais. Revista Árvore , v. 30, n. 4, p. 707-718, 2012.
- SOUZA, G. S.; LIMA, J. S. S.; SILVA, A. S.; OLIVEIRA, R. B. Variabilidade espacial de atributos químicos em um Argissolo sob Pastagem. Acta Scientiarum Agronomy, v. 30, n. 4, p. 589-596, 2008.
- SOUZA, Z. M.; MARQUES JÚNIOR, J.; PEREIRA, G.T. Influência da pedoforma na variabilidade espacial de alguns atributos físicos e hídricos de um Latossolo sob o cultivo da cana-de-açúcar. Irriga, Piracicaba, v. 9, n.1, p.1-11, 2004.
- SOUZA, C. K,; MARQUES JÚNIOR, J.; MARTINS FILHO, M. V.; PEREIRA, G. T. Influência do relevo e erosão na variabilidade espacial de um Latossolo em Jaboticabal (SP). Revista Brasileira de Ciência do Solo , v. 27, n. 1, p. 1067-74, 2003.
- TROEH, F. R. Landform equations fitted to contour maps. American Journal of Science, 263:616-27, 1965.
- VIEIRA, S. R. Geoestatística em estudos de variabilidade espacial do solo. In: Novais, R. F.; Alvarez, VH & Schaefer, G.R.; editores. Tópicos em Ciência do Solo. Viçosa, MG: Sociedade Brasileira de Ciência do Solo, v. 1, n. 1, p. 1-54, 2000.
Datas de Publicação
-
Publicação nesta coleção
Oct-Dec 2016
Histórico
-
Recebido
05 Jun 2016 -
Aceito
07 Out 2016