ABSTRACT Pedosphere-hydrosphere interface accounts for the association between soil hydrology and landscape, represented by topographic and Remote Sensing data support and integration. This study aimed to analyze different statistical radiometric and spectral data selection methods and dimensionality reduce environment-related data to support the classification of soil physical-hydric properties, such as soil basic infiltration rate (bir) and saturated hydraulic conductivity (Ksat); as well as to act in data mining processes applied to hydropedological properties digital mapping. Accordingly, research integrated information from Visible to Infrared (VIS-IR) spectral indices and Sentinel’s 2A mission Multispectral Instrument (MSI) sensor bands, terrain numerical modeling and aerogeophysics set to model soil-water content in two soil layers (0.00-0.20 m and 0.20-0.40 m). Pre-processed data were subjected to statistical analysis (multivariate and hypothesis tests); subsequently, the methods were applied (variation inflation factor - VIF, Stepwise Akaike information criterion – Stepwise AIC, and recursive feature elimination - RFE) to mine covariates used for Random Forest modeling. Based on the results, there were distinctions and singularities in spectral and radiometric data selection for each adopted method; the importance degree, and contribution of each one to soil physical-hydric properties have varied. According to the applied statistical metrics and decision-making criteria (highest R2 and lowest RMSE / MAE), the chosen methods were RFE (0.00-0.20 m layers) and Stepwise AIC (0.20-0.40 m layers) - both concerned with the assessed variables (bir and Ksat). This approach captured the importance of environmental variables and highlighted their potential use in hydropedological digital mapping at Guapi-Macacu watershed.

Abstract Although the conservation of tropical biodiversity depends on protected areas, there is still a very large ‘gap’ of knowledge on the flora of Brazilian reserves, especially in the Northeast region of Brazil. Field and herbarium surveys of the phanerogamic flora of the Ubajara National Park, located on the Brazilian Northeast, were made and analyses on phenology and dispersal syndromes were performed. 418 taxa (213 trees and shrubs, 100 terrestrial herbs, 68 climbing plants, 33 sub-shrubs, two epiphytes, one hemiparasite and one aquatic herb) were recorded. The most representative families were: Fabaceae, Malvaceae, Asteraceae, Rubiaceae and Euphorbiaceae. The annual flowering / fruiting peak hypothesis was not fully confirmed, therefore, the forest may be an important food resource for the fauna all year long (especially in the moister region). Zoochory was the predominant dispersal syndrome in the moister area, whereas, autochory and anemochory together, predominated in the drier area.

ABSTRACT Notwithstanding the importance of soil surveys, advances in digital soil mapping have mainly focused on mapping soil attributes or properties rather than developing digital maps of soil units or soil classes. The purpose of this research was to develop digital soil unit maps based on primary soil data collection in areas without previously collected soil information. The covariate variability, the random effect across the data subset and the map outputs were the focuses of this study. We used five datasets with four models (Random Forest - RF, Gradient Boosted Machine - GBM, C5.0, and multinomial log-linear model - MLR). The covariates were grouped into five datasets, where four were grouped by Region Of Interest per Class (ROIC) and one was not grouped by ROIC. To evaluate the random effect to split the dataset, we ran each model 50 times and observed the overall accuracy (OA) and kappa index, and uncertainty, majority and variety maps. The OA of Dataset01 to 04 was lower than to Dataset05 accuracy. However, map outputs of RF and GBM for Dataset01 and Dataset05 had the same majority prediction. It seems that RF and GBM produce consistent results in map outputs according to this methodology and pedologist expertise. To evaluate the uncertainty and the consistency of soil unit prediction, we used the majority maps process. Random Forest, similar to GBM, presented the best results. The increase in the number of covariates was not a guarantee of improvement in the OA or in the quality of the map output. Geographic position and distance raster did not improve the map output according to expert evaluation. Because the variance between the ROICs, when the training and validation datasets were split based on it, the subsets are quite different in relation to the covariates, and this is the reason for the worse results of this model, comparing with the Dataset05. On the other hand, when considering one complete dataset not based on ROICs, the variance of training and validation subsets is lower and produced more accurate parameters of quality.

ABSTRACT: Users of soil survey products are mostly interested in understanding how soil properties vary in space and time. The aim of digital soil mapping (DSM) is to represent the spatial variability of soil properties quantitatively to support decision-making. The goal of this study is to evaluate DSM techniques (Regression Trees - RT and Multiple Linear Regressions - MLR) and the ability of these tools to predict mineral fraction content under a wide variability of landscapes. The study site was the entire Guapi-Macacu watershed (1,250.78 km2) in the state of Rio de Janeiro in the Southeast region of Brazil. Terrain attributes and remote sensing data (with 30 m of spatial resolution) were used to represent landscape co-variables selected as an input in predictive models in order to develop the explanatory variables. The selection of sampling sites was based on the Latin Hypercube algorithm. A representative set of one hundred points with feasible field access was chosen. Different input databases were tested for prediction of mineral fraction content (harmonized and original data). The Spline algorithm was used to harmonize data according to the GlobalSoil. Net consortium standards. The results showed better performance from the RT models, using input from an average of six covariates; the simplest MLR model used twice as many input variables, creating more complex models without gaining precision. Furthermore, better R2 values were obtained using RT models, irrespective of harmonization of soil data. The harmonized dataset from the 0.00-0.05 and 0.05-0.15 m layers, in general, presented better results for the clay and silt, with R2 values of 0.52 (0.00-0.05 m) and 0.69 (0.05-0.15 m), respectively. Prediction of sand content showed better results when the original depth data was used as an input, although all regression tree models had R2 values greater than 0.52. The RT models provided a better statistical index than MLR for all predicted properties; however, the variance between models suggests similarity of performance. Regarding harmonization of soil data, both input databases (harmonized or not) can be used to predict soil properties, since the variance of model performance was low and generalization of the soil maps showed similar trends. The products obtained from the digital soil mapping approach make it possible to integrate the factor of uncertainties, providing easier interpretation for soil management and land use decisions.

ABSTRACT: Planning sustainable use of land resources requires reliable information about spatial distribution of soil physical and chemical properties related to environmental processes and ecosystemic functions. In this context, cation exchange capacity (CEC) is a fundamental soil quality indicator; however, it takes money and time to obtain this data. Although many studies have been conducted to spatially quantify soil properties on various scales and in different environments, not much is known about interactions between soil properties and environmental covariates in the Brazilian semiarid region. The goal of this study was to evaluate the efficiency of random forest and cokriging models applied to predict CEC in the Brazilian semiarid region. The covariates used to predict CEC consist of images from Landsat 5 TM and a legacy soil map (scale 1:10,000). The sample set comprises 499 samples from the topsoil layer (0.00-0.20 m), where 375 samples were used in training processes and 124 as validation samples. The cokriging model (R2 = 0.57 and RMSE = 7.22 cmolc kg-1) performed better in predicting CEC than the random forest model (R2 = 0.47 and RMSE = 7.89 cmolc kg-1). The approach used showed potential for estimating CEC content in the Brazilian semiarid region by using covariates obtained from orbital remote sensing and the legacy soil map.

ABSTRACT: The pattern of variation in soil and landform properties in relation to environmental covariates are closely related to soil type distribution. The aim of this study was to apply digital soil mapping techniques to analysis of the pattern of soil property variation in relation to environmental covariates under periglacial conditions at Keller Peninsula, Maritime Antarctica. We considered the hypothesis that covariates normally used for environmental correlation elsewhere can be adequately employed in periglacial areas in Maritime Antarctica. For that purpose, 138 soil samples from 47 soil sites were collected for analysis of soil chemical and physical properties. We tested the correlation between soil properties (clay, potassium, sand, organic carbon, and pH) and environmental covariates. The environmental covariates selected were correlated with soil properties according to the terrain attributes of the digital elevation model (DEM). The models evaluated were linear regression, ordinary kriging, and regression kriging. The best performance was obtained using normalized height as a covariate, with an R2 of 0.59 for sand. In contrast, the lowest R2 of 0.15 was obtained for organic carbon, also using the regression kriging method. Overall, results indicate that, despite the predominant periglacial conditions, the environmental covariates normally used for digital terrain mapping of soil properties worldwide can be successfully employed for understanding the main variations in soil properties and soil-forming factors in this region.

Resumo O objetivo deste trabalho foi avaliar ferramentas pedométricas para caracterização da composição granulométrica (areia, silte e argila) e comparação de perfis de solos hidromórficos. O estudo foi realizado na bacia hidrográfica do rio Guapi-Macacu, RJ. Foram aplicados o algoritmo "slice-wise" para o fatiamento do perfil em intervalos de 1 cm, e a função "spline" para harmonizar o conjunto de dados segundo seis intervalos de profundidade predefinidos. A análise do conjunto de perfis, por meio das funções de profundidade (slice-wise e spline) e de dendrograma de dissimilaridade, revelou que Planossolos e Gleissolos são relativamente rasos e apresentam camadas arenosas nos horizontes superficiais e teores de argila e silte irregulares em subsuperficie. Essas duas ordens ocorrem na área da bacia hidrográfica geralmente em planícies aluviais, em vales e em sopés de elevação. Os Neossolos Flúvicos apresentaram maiores quantidades de argila e silte em superfície, as quais decrescem em profundidade, e são mais profundos e ocorrem nos terraços aluviais. A abordagem permitiu descrever a variabilidade da textura dos solos hidromórficos em perfil e realizar a análise de similaridade entre as ordens, o que pode apoiar a distinção de fases em unidades de mapeamento e a identificação de critérios quantitativos para a classificação dos solos.

Abstract The objective of this work was to evaluate pedometric tools to characterize the particle size distribution (sand, silt, and clay) and to compare the profiles of hydromorphic soils. The study was carried out in the Guapi-Macacu watershed, in the state of Rio de Janeiro, Brazil. The slice-wise algorithm was applied to slice the profile data at 1-cm intervals, and the equal-area spline function for database harmonization, according to six predefined depth intervals. The analysis of the soil profile collection, using the soil-depth functions (slice-wise and spline) and the similarity analysis, revealed that Alfisols and Entisols are relatively shallow and have coarse texture in the topsoil layer and a more irregular distribution of clay and silt along the soil profile. These two orders usually occur in the watershed area related to floodplains, valleys, and footslopes. Fluvisols had higher amounts of clay and silt in the topsoil layer, which decrease with soil depth, and they are deeper and occur in alluvial terraces. The approach allowed to characterize the variability of texture along hydromorphical soil profiles and to perform a similarity analysis between soil orders, which can support the differentiation of soil mapping units and the identification of quantitative criteria for soil classification.

Resumo O objetivo deste trabalho foi avaliar a influência da resolução espacial do modelo digital de elevação e da eficiência de modelos Random Forest sobre a predição dos teores de areia, argila e carbono orgânico, com uso de número reduzido de amostras. O trabalho foi realizado em área de Cerrado com diversidade litológica, no Estado do Mato Grosso do Sul, tendo-se utilizado atributos morfométricos, dados do sensor TM Landsat 5 e litologia como covariáveis preditoras. Dados da camada superficial (0,0-0,2 m) de 175 perfis de solos (0,009 perfis km-2) e de 26 covariáveis preditoras foram utilizados com resolução espacial de 30 (conjunto 1) e 90 m (conjunto 2). A análise realizada pelo Random Forest mostrou que as covariáveis de nível de base do canal de drenagem, da elevação e da litologia foram as mais importantes para explicar a variabilidade. A validação dos modelos apresentou similaridade entre os conjuntos quanto à predição de areia, argila e carbono orgânico, o que explica os seguintes valores de variabilidade espacial, respectivamente: 44, 40 e 33%, para a resolução de 30 m; e de 45, 46 e 33%, para a resolução de 90 m. A resolução espacial das covariáveis preditoras tem pouca influência sobre a predição dos atributos, e a abordagem por Random Forest apresenta potencial de utilização para estimar atributos do solo.

Abstract The objective of this work was to evaluate the effect of the digital elevation model spatial resolution and of the efficiency of Random Forest models on the prediction of sand, clay, and organic carbon contents, using few soil samples. The study was carried out in a Cerrado area with lithological diversity, in the state of Mato Grosso do Sul, Brazil, using morphometric attributes, TM Landsat 5 sensor data, and lithology as predictive covariates. The surface layer data (0.0-0.2 m) of 175 soil profiles (0,009 profiles km-2) and of 26 predictor covariates were used with 30 (set 1) and 90-m (set 2) spatial resolutions. The performed analysis by Random Forest models showed that channel base level, elevation, and lithology were the most important ones to explain the variability. The validation of the models showed similarity among sets for the prediction of sand, clay, and organic carbon contents, which explains the following values of spatial variability, respectively: 44, 40, and 33%, for the spatial resolution of 30 m; and 45, 46, and 33%, for the spatial resolution of 90 m. The spatial resolution of the predictive covariates has little effect on attribute predictions, and the Random Forest approach has potential use for estimating soil properties.

Resumo O objetivo deste trabalho foi o desenvolvimento de modelos com diferentes conjuntos de dados, para estimar a densidade de solos de regiões tropicais montanhosas, a partir de atributos de solos comumente encontrados nas análises de perfis de solos descritos nos levantamentos regionais. O conjunto total de dados compõe-se de 163 amostras e foi dividido em seis grupamentos, dos quais três com 73 amostras, com o máximo de 32 covariáveis, e três com 163 amostras, com o máximo de 18 covariáveis. Testaram-se modelos de regressão linear múltipla (RLM) e randomForest (RF). A menor incerteza entre os modelos foi alcançada pelo RLM2, com R2 de 0,56, 13 covariáveis e 73 amostras. Nos grupamentos com 163 amostras, os melhores modelos foram os RF, com R2 médio de 0,48. A raiz quadrada da média do erro ao quadrado variou entre 0,09 e 0,14. As covariáveis mais importantes no modelo RF foram: carbono orgânico, hidrogênio, areia fina e grossa, saturação por bases e capacidade de troca catiônica. Pelo método "stepwise regression", as variáveis mais importantes foram: a relação silte/argila; areia grossa e fina; carbono orgânico; saturação por bases; e potássio.

Abstract The objective of this work was the development of models with different sets of data for estimating soil bulk density in tropical mountainous regions, from soil attributes commonly found in the analyses of soil profiles described in regional surveys. The complete dataset is composed of 163 samples and it was divided into six groups, of which three groups have 73 samples and the maximum of 32 covariables, and three have 163 samples and the maximum of 18 covariables. The linear regression (RLM) and randomForest (RF) models were tested. The lowest uncertainty between the models was achieved by RLM2, with R2 of 0.56, 13 covariables, and 73 samples. Considering the groups with 163 samples, the best models were the RFs with mean R2 of 0.48. The root mean squared error ranged between 0.09 and 0.14. The most important covariables in the RF model were: organic carbon, hydrogen, fine and coarse sand, base saturation, and cation exchange capacity. By the stepwise backward regression, the main covariables were: silt and clay relation; fine and coarse sand; organic carbon; base saturation; and potassium.

Resumo:O objetivo deste trabalho foi estabelecer frações ótimas da dosagem total de fósforo, a serem aplicadas no pré-plantio e na fertirrigação do tomateiro de crescimento determinado, bem como determinar o efeito das diferentes formas de aplicação do nutriente sobre sua distribuição espacial no solo. Os tratamentos consistiram de cinco esquemas de aplicação do fertilizante: 100, 75, 50, 25 e 0% da dose total de P, dividida entre o pré-plantio (superfosfato triplo incorporado ao solo) e a fertirrigação (ácido fosfórico). O experimento foi realizado em Latossolo Vermelho argiloso com duas classes de disponibilidade de P - 7,5 e 32,9 mg dm-3 (Mehlich-1). Os tratamentos não apresentaram efeito sobre as variáveis de produção do tomateiro cultivado no solo com alta disponibilidade de P. Entretanto, no solo com menor disponibilidade do nutriente, a fertirrigação com P apresentou efeito quadrático sobre a produtividade de frutos, a receita líquida e a taxa líquida de retorno. Os melhores resultados foram obtidos com a aplicação de metade da dose total de P pela fertirrigação. A fertirrigação possibilitou movimento vertical do nutriente no solo para até 30 cm de profundidade, e movimento horizontal de até 15 cm.

Abstract:The objective of this work was to establish the optimal fractions of the total amount of phosphorus to be applied in preplanting and in drip fertigation of determinate-growth tomato, as well as to determine the effect of the different forms of nutrient application on their soil spatial distribution. Treatments consisted of five schemes of fertilizer application: 100, 75, 50, 25, and 0% of P total dose divided between the pre-planting (triple superphosphate in band placement) and the fertigation (phosphoric acid). The experiment was carried out in a clayed Oxisol with two P availability classes - 7.5 and 32.9 mg dm-3 (Mehlich-1). The treatments showed no effect on tomato production variables in the soil with high P availability. However, in the soil with smaller availability of the nutrient, P fertigation showed a quadratic effect on fruit yield, net income, and net rate of return. The best results were obtained with half of the P total dose applied by fertigation. The fertigation allowed P movement in the soil up to 30 cm deep, and horizontal movement up to 15 cm.

A informação de solo é necessária para o gerenciamento do ambiente agrícola. O objetivo deste trabalho foi aplicar redes neurais artificiais (RNAs) para a predição de classes de solos, utilizando como fonte de dados produtos de sensores remotos orbitais, atributos do terreno derivados de um modelo digital de elevação e informação da geologia local, visando avaliar a utilização dessa abordagem no mapeamento digital de solos, em área com elevado grau de diversidade litológica na Serra do Mar. O simulador de redes neurais utilizado foi o JavaNNS e o algoritmo de aprendizado, o backpropagation . Para a predição das classes de solos, testaram-se diferentes combinações entre as variáveis discriminantes selecionadas: elevação, declividade, aspecto, curvatura, plano de curvatura, perfil de curvatura, índice topográfico, radiação solar, fator topográfico LS, informações da geologia local e índices minerais de argila, óxidos de ferro e vegetação por diferença normalizada (NDVI), derivados de uma imagem do sensor ETM+ do LANDSAT 7. Dos conjuntos testados, os melhores resultados foram obtidos com todas as variáveis discriminantes associadas às informações de geologia, alcançando exatidão global entre 93,2 e 95,6 % e índice Kappa entre 0,924 e 0,951 (conjunto 13). Excluindo a variável perfil de curvatura (conjunto 12), a exatidão global alcançada oscilou entre 93,9 e 95,4 % e o Kappa entre 0,932 e 0,948. Os mapas inferidos pelo classificador por redes neurais evidenciaram coerência e semelhança com o mapa de solos convencional, produzido para área de estudo, apresentando porém, mais detalhes espaciais. Os resultados apresentaram o potencial de utilização de RNAs na predição de classes de solos de áreas montanhosas com diversidade litológica.

Soil information is needed for managing the agricultural environment. The aim of this study was to apply artificial neural networks (ANNs) for the prediction of soil classes using orbital remote sensing products, terrain attributes derived from a digital elevation model and local geology information as data sources. This approach to digital soil mapping was evaluated in an area with a high degree of lithologic diversity in the Serra do Mar. The neural network simulator used in this study was JavaNNS and the backpropagation learning algorithm. For soil class prediction, different combinations of the selected discriminant variables were tested: elevation, declivity, aspect, curvature, curvature plan, curvature profile, topographic index, solar radiation, LS topographic factor, local geology information, and clay mineral indices, iron oxides and the normalized difference vegetation index (NDVI) derived from an image of a Landsat-7 Enhanced Thematic Mapper Plus (ETM+) sensor. With the tested sets, best results were obtained when all discriminant variables were associated with geological information (overall accuracy 93.2 - 95.6 %, Kappa index 0.924 - 0.951, for set 13). Excluding the variable profile curvature (set 12), overall accuracy ranged from 93.9 to 95.4 % and the Kappa index from 0.932 to 0.948. The maps based on the neural network classifier were consistent and similar to conventional soil maps drawn for the study area, although with more spatial details. The results show the potential of ANNs for soil class prediction in mountainous areas with lithological diversity.

As propriedades dos solos têm grande impacto sobre aspectos econômicos e ambientais da produção agropecuária. As relações quantitativas entre as propriedades dos solos e os fatores que condicionam sua variabilidade são a base do mapeamento digital de solos. Os modelos preditivos de propriedades dos solos avaliados neste trabalho são os estatísticos (Regressão Linear Múltipla-RLM) e geoestatísticos (krigagem ordinária e cokrigagem). Este estudo foi desenvolvido para o município de Bom Jardim, RJ, e usou um banco de dados de solos com 208 pontos amostrais. Foram avaliados modelos preditivos para as frações areia, silte e argila, pH em água e carbono orgânico para seis profundidades, de acordo com as especificações do consórcio de mapeamento digital de solos em nível global (GlobalSoilMap). Utilizaram-se covariáveis preditoras contínuas e categóricas, estas últimas para avaliar suas contribuições ao modelo. Apenas as covariáveis ambientais elevação, aspecto, índice de potência de fluxo (SPI), índice de umidade (SWI), índice de vegetação por diferença normalizada (NDVI) e relação entre bandas b3/b2 apresentaram correlação significativa com as propriedades do solo. Os modelos preditivos tiveram em média coeficiente de determinação de 0,21. Os modelos preditivos que apresentaram os melhores resultados foram os geoestatísticos, com o maior coeficiente de determinação 0,43 associado à propriedade areia entre 60 e 100 cm de profundidade. A utilização de conjunto de dados de solos esparsos para mapeamento digital de propriedades de solos pode explicar apenas uma parte da variação espacial dessas propriedades. Os resultados podem estar relacionados à densidade de amostragem, à quantidade e qualidade das covariáveis ambientais usadas e aos modelos preditivos utilizados.

Soil properties have an enormous impact on economic and environmental aspects of agricultural production. Quantitative relationships between soil properties and the factors that influence their variability are the basis of digital soil mapping. The predictive models of soil properties evaluated in this work are statistical (multiple linear regression-MLR) and geostatistical (ordinary kriging and co-kriging). The study was conducted in the municipality of Bom Jardim, RJ, using a soil database with 208 sampling points. Predictive models were evaluated for sand, silt and clay fractions, pH in water and organic carbon at six depths according to the specifications of the consortium of digital soil mapping at the global level (GlobalSoilMap). Continuous covariates and categorical predictors were used and their contributions to the model assessed. Only the environmental covariates elevation, aspect, stream power index (SPI), soil wetness index (SWI), normalized difference vegetation index (NDVI), and b3/b2 band ratio were significantly correlated with soil properties. The predictive models had a mean coefficient of determination of 0.21. Best results were obtained with the geostatistical predictive models, where the highest coefficient of determination 0.43 was associated with sand properties between 60 to 100 cm deep. The use of a sparse data set of soil properties for digital mapping can explain only part of the spatial variation of these properties. The results may be related to the sampling density and the quantity and quality of the environmental covariates and predictive models used.

A amostragem é uma das etapas mais importantes dos levantamentos de solos. No entanto, os esquemas de amostragem utilizados nos levantamentos convencionais têm se evidenciado inadequados para o mapeamento digital de solos, pois podem comprometer os resultados e, além disso, não possibilitam a realização de análises estatísticas. Este estudo teve por objetivo avaliar o método de amostragem do hipercubo latino condicionado (cLHS, sigla em inglês), na presença de covariáveis ambientais (elevação, declividade, curvatura e mapa de uso e cobertura do solo), em comparação com a amostragem aleatória, na alocação de 100 pontos amostrais, buscando maior representatividade das características ambientais da bacia do rio Guapi-Macacu. O desempenho dos métodos foi avaliado pela análise qualitativa dos histogramas de frequência e das análises estatísticas pelos testes F, T de Student e Kolmogorov-Smirnov (K-S), para cada covariável. Os resultados apresentaram que os pontos selecionados pelo método cLHS possuíam distribuição geográfica mais adequada do que aqueles obtidos pela amostragem aleatória. Além disso, o método cLHS preservou mais a distribuição de frequência das covariáveis contínuas do que a amostragem aleatória; para covariável categórica uso e cobertura do solo os métodos foram equivalentes. Os testes estatísticos confirmaram o melhor desempenho do método cLHS, cujas amostras não diferiram estatisticamente da bacia. Entretanto, a amostragem aleatória apresentou diferença estatística para com a bacia, para todas as covariáveis contínuas para pelo menos um dos testes utilizados. Assim, o método cLHS pode ser considerado como um método satisfatório para seleção de locais de amostragem em áreas heterogêneas similares as deste estudo, visando a utilização no mapeamento digital de solos.

Soil sampling is one of the most important stages in soil surveys. However, in conventional soil surveys, sampling schemes have proven to be inadequate for digital soil mapping because they may affect the results obtained and, in addition, they do not allow statistical estimations to be performed. This study attempts to evaluate the conditioned Latin Hypercube Sampling (cLHS) method with four environmental covariates (elevation, slope, curvature and land cover/land use map), in comparison to random sampling, in allocating 100 sample points, seeking a more representative description of the environmental characteristics of the Macacu watershed area. The performance of the methods was evaluated by qualitative assessment of histograms of distribution, and statistical analyses by the F test, Student T Test and Kolmogorov-Smirnov test (K-S) for each covariable. The sampling points selected by the cLHS method showed better geographic distribution than those obtained by random sampling. Furthermore, the cLHS method better maintained the frequency distribution of the continuous covariables than random sampling did, and the categorical covariable of land cover/land use showed no significance differences between the cLHS and random methods. Statistical analysis confirmed the better performance of the cLHS method compared to the random method, and the cLHS samples showed no significance differences from the watershed data. In contrast, the random samples showed significance differences from the watershed data for all the continuous covariables for at least one of the statistical tests used. Thus, the cLHS method is a satisfactory method for selecting points for soil sampling in heterogeneous areas, like those of the present study, for the purpose of use in digital soil mapping.

O trabalho teve como objetivo apresentar uma proposta metodológica para o zoneamento agroecológico do Estado do Mato Grosso do Sul, tendo como estudo de caso o município de Bandeirantes. Foram utilizadas geotecnologias, tais como: imagens de satélite (Landast 5), sistema de posicionamento global, sistemas de informação geográfica, cartografia digital e software de processamento digital de imagens, visando construir parâmetros ambientais para uma análise integrada de dados, associados com atributos do terreno derivados de um modelo digital de elevação, dados de propriedades dos solos e de clima. Estas geotecnologias foram empregadas para calcular e espacializar o potencial natural de erosão das terras, a fertilidade natural, drenagem interna e a capacidade de retenção de umidade dos solos. A análise integrada dos parâmetros ambientais, do uso atual das terras e da espacialização da legislação ambiental, permitiu estratificar o município em diferentes zonas agroecológicas recomendadas para o uso agrícola; para o uso com pastagens; para a conservação dos recursos naturais e para a recuperação ambiental. No município de Bandeirantes as zonas agroecológicas recomendadas para o uso com lavouras somam 1.150 km2 (37%) e para pastagens 1.658 km2 (53%), As recomendadas para conservação dos recursos naturais e para recuperação ambiental correspondem a 310 km2, isto é, 10% das terras do município. Pôde-se concluir que a análise integrada de dados ambientais associados às geotecnologias possibilitou a avaliação quantitativa e qualitativa permitindo à estratificação do município em diferentes zonas agroecológicas com indicação de áreas passíveis de exploração agrícola sustentável

This study aimed to present the methodological approach of the agroecological zoning of Mato Grosso do Sul State. The approach incorporated quantitative and qualitative assessment, with geotechnologies (satellite image–landsat 5, global positioning system, geographic information system, digital mapping, software of digital image processing) adoption to perform an integrated analysis of environmental data. Primary soils information associated with terrain attributes and climate characteristics were used with geotechnologies to evaluation of the natural erosion potential, natural fertility, internal drainage and moisture soil retention capacity of soil. The land use, and land cover and the environmental laws were observed. The integrated analysis of these data allowed environmental stratification in different agroecological zones. These zones are defined to agricultural use; pastures; conservation of natural resources or ecological restoration. In Bandeirantes case study, the agroecological zones recommended for agriculture sum 1,150 km2. The evaluation for conservation of natural resources and ecological results in 310 km2, and for the areas recommended for pastures sum 1,650 km2 or 53% of Bandeirantes municipality. This study allows conclude that the integrated analysis of environmental data with application of geotechnologies can subsidize a quantitative and qualitative assessment, stratifying the lands in different agroecological zones, each one with recommendations for sustainable use, and contributing as a powerful tool for spatial planning.

A topografia do terreno controla muitos processos hidrológicos, geomorfológicos e pedológicos, influenciando a radiação solar, precipitação pluvial, escoamento superficial, evaporação, regime de umidade do solo e tipo de vegetação, entre outros. Desta maneira, atributos topográficos derivados de um modelo digital de elevação foram utilizados para explicar a diferenciação de solos com horizonte B textural da região de ''Mar de Morros'', no noroeste do Estado do Rio de Janeiro. Foram estudados oito perfis de solo (quatro Argissolos Vermelhos e quatro Argissolos Vermelho-Amarelos), desenvolvidos de granulitos noríticos nos municípios de São José de Ubá e Itaperuna. Os dados dos atributos topográficos em cada perfil de solo (elevação, declividade, aspecto, plano e perfil de curvatura, radiação solar e índice topográfico combinado) foram extraídos, utilizando-se uma janela de 7 x 7 pixels, com os perfis ocupando o centro dessa janela. Os Argissolos que ocorrem nas encostas noroeste e nordeste recebem maior quantidade de radiação solar (mais quentes e secas) e por isto são vermelhos, eutróficos e não apresentam horizonte B latossólico, enquanto os solos das encostas sudeste (relativamente mais frias e úmidas) são vermelho-amarelados, distróficos e com horizonte B latossólico. A radiação solar, a declividade e o aspecto são os atributos topográficos relacionados com as diferenças entre os perfis de solos estudados. Os atributos topográficos analisados pelas técnicas multivariadas apresentaram-se efetivos para a distinção entre os Argissolos Vermelho-Amarelos distróficos e os Argissolos Vermelhos eutróficos, exceto o P7 (PVe). Estudos relacionando a influência da radiação solar na pedogênese são escassos em condições tropicais e, embora menos atuante em regiões tropicais, foi a radiação solar que determinou diretamente a diferenciação dos solos da área estudada.

The topographic conditions controls many hydrological, geomorphological and pedological processes, affectting the solar radiation, pluvial precipitation, superficial water flow, evaporation, soil moisture and vegetation types, among others. Topographic properties derived from a digital elevation model were used to explain the differentiation between soils with an argillic horizon in the region of ''Mar de Morros'' in the northwestern State of Rio de Janeiro. Eight soil profiles (four Red Argisols and four Red-Yellow Argisols) developed from noritic granulites, in the municipalities of São José de Ubá and Itaperuna, were studied. The topographic data of each soil profile (elevation, slope, aspect, curvature plan and profile, solar radiation and combined topographic index) were extracted using a 7 x 7 pixel window, with the profiles in the center. Argisols on slopes facing northwest and northeast receive most solar radiation (warmer and drier) and for this reason are red, eutrophic and have no oxic horizon, whereas soils on slopes facing southeast (relatively cold and wet) are yellowish-red, dystrophic and have an oxic horizon. The solar radiation, slope and aspect are the topographic attributes that differ between the soils profiles. The topographic properties analyzed by multivariate techniques were effective to distinguish eutrophic Red Argisols from dystrophic Red-Yellow Argisols, except for P7 (PVe). Studies on the influence of solar radiation on pedogenesis under tropical conditions are scarce, and although less active in tropical regions directly determined the soil differentiation of the studied area.