Variability of apparently homogeneous soilscapes in São Paulo state, Brazil: II. quality of soil maps

The quality of semi-detailed (scale 1:100.000) soil maps and the utility of a taxonomically based legend were assessed by studying 33 apparently homogeneous fields with strongly weathered soils in two regions in São Paulo State: Araras and Assis. An independent data set of 395 auger sites was used to determine purity of soil mapping units and analysis of variance within and between mapping units and soil classification units. Twenty three soil profiles were studied in detail. The studied soil maps have a high purity for some legend criteria, such as B horizon type (> 90%) and soil texture class (> 80%). The purity for the "trophic character" (eutrophic, dystrophic, allic) was only 55% in Assis. It was 88% in Araras, where many soil units had been mapped as associations. In both regions, the base status of clay-textured soils was generally better than suggested by the maps. Analysis of variance showed that mapping was successful for "durable" soil characteristics such as clay content (> 80% of variance explained) and cation exchange capacity (≥ 50% of variance explained) of 0-20 and 60-80 cm layers. For soil characteristics that are easily modified by management, such as base saturation of the 0-20 cm layer, the maps had explained very little (< 15%) of the total variance in the study areas. Intermediate results were obtained for base saturation of the 60-80 cm layer (56% in Assis; 42% in Araras). Variance explained by taxonomic groupings that formed the basis for the legend of the soil maps was similar to, often even smaller than, variance explained by mapping units. The conclusion is that map boundaries have been very carefully located, but descriptions of mapping units could be improved. In future mappings, this could possibly be done at low cost by (a) bulk sampling to remove short range variation and enhance visualization of spatial patterns at distances > 100 m; (b) taking advantage of correlations between easily measured soil characteristics and chemical soil properties and, (c) unbending the link between legend criteria and a taxonomic system. The maps are well suited to obtain an impression of land suitability for high-input farming. Additional field work and data on former land use/management are necessary for the evaluation of chemical properties of surface horizons.

Variability of apparently homogeneous soilscapes in São Paulo state, Brazil: I. spatial analysis

The spatial variability of strongly weathered soils under sugarcane and soybean/wheat rotation was quantitatively assessed on 33 fields in two regions in São Paulo State, Brazil: Araras (15 fields with sugarcane) and Assis (11 fields with sugarcane and seven fields with soybean/wheat rotation). Statistical methods used were: nested analysis of variance (for 11 fields), semivariance analysis and analysis of variance within and between fields. Spatial levels from 50 m to several km were analyzed. Results are discussed with reference to a previously published study carried out in the surroundings of Passo Fundo (RS). Similar variability patterns were found for clay content, organic C content and cation exchange capacity. The fields studied are quite homogeneous with respect to these relatively stable soil characteristics. Spatial variability of other characteristics (resin extractable P, pH, base- and Al-saturation and also soil colour), varies with region and, or land use management. Soil management for sugarcane seems to have induced modifications to greater depths than for soybean/wheat rotation. Surface layers of soils under soybean/wheat present relatively little variation, apparently as a result of very intensive soil management. The major part of within-field variation occurs at short distances (< 50 m) in all study areas. Hence, little extra information would be gained by increasing sampling density from, say, 1/km² to 1/50 m². For many purposes, the soils in the study regions can be mapped with the same observation density, but residual variance will not be the same in all areas. Bulk sampling may help to reveal spatial patterns between 50 and 1.000 m.

Strategy of specification of management areas: rice grain yield as related to soil fertility

It is well-known nowadays that soil variability can influence crop yields. Therefore, to determine specific areas of soil management, we studied the Pearson and spatial correlations of rice grain yield with organic matter content and pH of an Oxisol (Typic Acrustox) under no- tillage, in the 2009/10 growing season, in Selvíria, State of Mato Grosso do Sul, in the Brazilian Cerrado (longitude 51º24' 21'' W, latitude 20º20' 56'' S). The upland rice cultivar IAC 202 was used as test plant. A geostatistical grid was installed for soil and plant data collection, with 120 sampling points in an area of 3.0 ha with a homogeneous slope of 0.055 m m-1. The properties rice grain yield and organic matter content, pH and potential acidity and aluminum content were analyzed in the 0-0.10 and 0.10-0.20 m soil layers. Spatially, two specific areas of agricultural land management were discriminated, differing in the value of organic matter and rice grain yield, respectively with fertilization at variable rates in the second zone, a substantial increase in agricultural productivity can be obtained. The organic matter content was confirmed as a good indicator of soil quality, when spatially correlated with rice grain yield.

Management zones in agriculture acording to the soil and landscape variables

The study of spatial variability of soil and plants attributes, or precision agriculture, a technique that aims the rational use of natural resources, is expanding commercially in Brazil. Nevertheless, there is a lack of mathematical analysis that supports the correlation of these independent variables and their interactions with the productivity, identifying scientific standards technologically applicable. The aim of this study was to identify patterns of soil variability according to the eleven physical and seven chemical indicators in an agricultural area. It was used two multivariate techniques: the hierarchical cluster analysis (HCA) and the principal component analysis (PCA). According to the HCA, the area was divided into five management zones: zone 1 with 2.87ha, zone 2 with 0.8ha, zone 3 with 1.84ha, zone 4 with 1.33ha and zone 5 with 2.76ha. By the PCA, it was identified the most important variables within each zone: V% for the zone 1, CTC in the zone 2, levels of H+Al in the zone 4 and sand content and altitude in the zone 5. The zone 3 was classified as an intermediate zone with characteristics of all others. According to the results it is concluded that it is possible to separate into groups (management zones) samples with the same patterns of variability by the multivariate statistical techniques.

Strategy of specification of management areas: rice grain yield as related to soil fertility

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Cross-hole and seismic CPT tests in a tropical soil site

A large portion of Brazil is covered with tropical soils but literature about dynamic parameters of these soils is very limited. SCPT and cross-hole tests were carried out at an experimental research site inland in the state of São Paulo, Brazil. Shear wave velocities (VS) determined based on both tests are presented and compared. A good agreement was observed between both test results and the differences can be associated with soil variability, which was very sensitive to CPT tests. It was also observed that Go/q c ratio determined based on SCPT appears to be an interesting technique to help identify tropical soils. Copyright ASCE 2006.

Digital soil mapping: evaluation of sampling systems for soil surveys and refinement of soil maps at lower cost using legacy data

In soil surveys, several sampling systems can be used to define the most representative sites for sample collection and description of soil profiles. In recent years, the conditioned Latin hypercube sampling system has gained prominence for soil surveys. In Brazil, most of the soil maps are at small scales and in paper format, which hinders their refinement. The objectives of this work include: (i) to compare two sampling systems by conditioned Latin hypercube to map soil classes and soil properties; (II) to retrieve information from a detailed scale soil map of a pilot watershed for its refinement, comparing two data mining tools, and validation of the new soil map; and (III) to create and validate a soil map of a much larger and similar area from the extrapolation of information extracted from the existing soil map. Two sampling systems were created by conditioned Latin hypercube and by the cost-constrained conditioned Latin hypercube. At each prospection place, soil classification and measurement of the A horizon thickness were performed. Maps were generated and validated for each sampling system, comparing the efficiency of these methods. The conditioned Latin hypercube captured greater variability of soils and properties than the cost-constrained conditioned Latin hypercube, despite the former provided greater difficulty in field work. The conditioned Latin hypercube can capture greater soil variability and the cost-constrained conditioned Latin hypercube presents great potential for use in soil surveys, especially in areas of difficult access. From an existing detailed scale soil map of a pilot watershed, topographical information for each soil class was extracted from a Digital Elevation Model and its derivatives, by two data mining tools. Maps were generated using each tool. The more accurate of these tools was used for extrapolation of soil information for a much larger and similar area and the generated map was validated. It was possible to retrieve the existing soil map information and apply it on a larger area containing similar soil forming factors, at much low financial cost. The KnowledgeMiner tool for data mining, and ArcSIE, used to create the soil map, presented better results and enabled the use of existing soil map to extract soil information and its application in similar larger areas at reduced costs, which is especially important in development countries with limited financial resources for such activities, such as Brazil.

Soil organic carbon storage and soil properties for 50 soil profiles in the Lena River Delta including land form description and map

To project the future development of the soil organic carbon (SOC) storage in permafrost environments, the spatial and vertical distribution of key soil properties and their landscape controls needs to be understood. This article reports findings from the Arctic Lena River Delta where we sampled 50 soil pedons. These were classified according to the U.S.D.A. Soil Taxonomy and fall mostly into the Gelisol soil order used for permafrost-affected soils. Soil profiles have been sampled for the active layer (mean depth 58±10 cm) and the upper permafrost to one meter depth. We analyze SOC stocks and key soil properties, i.e. C%, N%, C/N, bulk density, visible ice and water content. These are compared for different landscape groupings of pedons according to geomorphology, soil and land cover and for different vertical depth increments. High vertical resolution plots are used to understand soil development. These show that SOC storage can be highly variable with depth. We recommend the treatment of permafrost-affected soils according to subdivisions into: the surface organic layer, mineral subsoil in the active layer, organic enriched cryoturbated or buried horizons and the mineral subsoil in the permafrost. The major geomorphological units of a subregion of the Lena River Delta were mapped with a land form classification using a data-fusion approach of optical satellite imagery and digital elevation data to upscale SOC storage. Landscape mean SOC storage is estimated to 19.2±2.0 kg C/m**2. Our results show that the geomorphological setting explains more soil variability than soil taxonomy classes or vegetation cover. The soils from the oldest, Pleistocene aged, unit of the delta store the highest amount of SOC per m**2 followed by the Holocene river terrace. The Pleistocene terrace affected by thermal-degradation, the recent floodplain and bare alluvial sediments store considerably less SOC in descending order.

DESARROLLO DE APLICACIONES ESTADÍSTICAS PARA LA AGRICULTURA DE PRECISIÓN. STATISTICAL APPLICATIONS FOR PRECISION AGRICULTURE

A partir de las últimas décadas se ha impulsado el desarrollo y la utilización de los Sistemas de Información Geográficos (SIG) y los Sistemas de Posicionamiento Satelital (GPS) orientados a mejorar la eficiencia productiva de distintos sistemas de cultivos extensivos en términos agronómicos, económicos y ambientales. Estas nuevas tecnologías permiten medir variabilidad espacial de propiedades del sitio como conductividad eléctrica aparente y otros atributos del terreno así como el efecto de las mismas sobre la distribución espacial de los rendimientos. Luego, es posible aplicar el manejo sitio-específico en los lotes para mejorar la eficiencia en el uso de los insumos agroquímicos, la protección del medio ambiente y la sustentabilidad de la vida rural. En la actualidad, existe una oferta amplia de recursos tecnológicos propios de la agricultura de precisión para capturar variación espacial a través de los sitios dentro del terreno. El óptimo uso del gran volumen de datos derivado de maquinarias de agricultura de precisión depende fuertemente de las capacidades para explorar la información relativa a las complejas interacciones que subyacen los resultados productivos. La covariación espacial de las propiedades del sitio y el rendimiento de los cultivos ha sido estudiada a través de modelos geoestadísticos clásicos que se basan en la teoría de variables regionalizadas. Nuevos desarrollos de modelos estadísticos contemporáneos, entre los que se destacan los modelos lineales mixtos, constituyen herramientas prometedoras para el tratamiento de datos correlacionados espacialmente. Más aún, debido a la naturaleza multivariada de las múltiples variables registradas en cada sitio, las técnicas de análisis multivariado podrían aportar valiosa información para la visualización y explotación de datos georreferenciados. La comprensión de las bases agronómicas de las complejas interacciones que se producen a la escala de lotes en producción, es hoy posible con el uso de éstas nuevas tecnologías. Los objetivos del presente proyecto son: (l) desarrollar estrategias metodológicas basadas en la complementación de técnicas de análisis multivariados y geoestadísticas, para la clasificación de sitios intralotes y el estudio de interdependencias entre variables de sitio y rendimiento; (ll) proponer modelos mixtos alternativos, basados en funciones de correlación espacial de los términos de error que permitan explorar patrones de correlación espacial de los rendimientos intralotes y las propiedades del suelo en los sitios delimitados. From the last decades the use and development of Geographical Information Systems (GIS) and Satellite Positioning Systems (GPS) is highly promoted in cropping systems. Such technologies allow measuring spatial variability of site properties including electrical conductivity and others soil features as well as their impact on the spatial variability of yields. Therefore, site-specific management could be applied to improve the efficiency in the use of agrochemicals, the environmental protection, and the sustainability of the rural life. Currently, there is a wide offer of technological resources to capture spatial variation across sites within field. However, the optimum use of data coming from the precision agriculture machineries strongly depends on the capabilities to explore the information about the complex interactions underlying the productive outputs. The covariation between spatial soil properties and yields from georeferenced data has been treated in a graphical manner or with standard geostatistical approaches. New statistical modeling capabilities from the Mixed Linear Model framework are promising to deal with correlated data such those produced by the precision agriculture. Moreover, rescuing the multivariate nature of the multiple data collected at each site, several multivariate statistical approaches could be crucial tools for data analysis with georeferenced data. Understanding the basis of complex interactions at the scale of production field is now within reach the use of these new techniques. Our main objectives are: (1) to develop new statistical strategies, based on the complementarities of geostatistics and multivariate methods, useful to classify sites within field grown with grain crops and analyze the interrelationships of several soil and yield variables, (2) to propose mixed linear models to predict yield according spatial soil variability and to build contour maps to promote a more sustainable agriculture.

Multiple Geotechnological Tools Applied to Digital Mapping of Tropical Soils

ABSTRACT In recent years, geotechnologies as remote and proximal sensing and attributes derived from digital terrain elevation models indicated to be very useful for the description of soil variability. However, these information sources are rarely used together. Therefore, a methodology for assessing and specialize soil classes using the information obtained from remote/proximal sensing, GIS and technical knowledge has been applied and evaluated. Two areas of study, in the State of São Paulo, Brazil, totaling approximately 28.000 ha were used for this work. First, in an area (area 1), conventional pedological mapping was done and from the soil classes found patterns were obtained with the following information: a) spectral information (forms of features and absorption intensity of spectral curves with 350 wavelengths -2,500 nm) of soil samples collected at specific points in the area (according to each soil type); b) obtaining equations for determining chemical and physical properties of the soil from the relationship between the results obtained in the laboratory by the conventional method, the levels of chemical and physical attributes with the spectral data; c) supervised classification of Landsat TM 5 images, in order to detect changes in the size of the soil particles (soil texture); d) relationship between classes relief soils and attributes. Subsequently, the obtained patterns were applied in area 2 obtain pedological classification of soils, but in GIS (ArcGIS). Finally, we developed a conventional pedological mapping in area 2 to which was compared with a digital map, ie the one obtained only with pre certain standards. The proposed methodology had a 79 % accuracy in the first categorical level of Soil Classification System, 60 % accuracy in the second category level and became less useful in the categorical level 3 (37 % accuracy).

Mapeamento das formas do relevo para estimativa de custos de fertilização em cana-de-açúcar

A modelagem matemática associada ao conhecimento da variabilidade dos atributos do solo e mapeamento das formas do relevo pode auxiliar no manejo da fertilidade do solo em usinas sucroalcooleiras. O presente trabalho teve como objetivo avaliar o uso da geoestatística e da modelagem matemática na estimativa de custos de fertilização, em diferentes formas do relevo. em uma área de 200 ha, foram identificadas duas formas de relevo, uma côncava e outra convexa, sendo os solos coletados nos pontos de cruzamento de uma malha, com intervalos regulares de 50 m, perfazendo um total de 623 pontos. As amostras foram submetidas a análises químicas, e, posteriormente, os dados foram avaliados por meio da estatística descritiva, geoestatística e modelagem matemática. Os resultados mostraram que, quando as formas do relevo são incorporadas às análises geoestatística e de modelagem matemática, ocorre aumento na eficiência de aplicação do calcário, fósforo e potássio no solo.

Superfícies geomórficas e atributos de Latossolos em uma seqüência Arenítico-Basáltica da região de Jaboticabal (SP)

Em um segmento de vertente com substrato de arenito em contato com basalto, regionalmente muito freqüente, pretendeu-se não só relacionar as superfícies geomórficas com os atributos físicos, químicos e mineralógicos dos Latossolos nelas encontrados, mas também testar métodos geoestatísticos para localização de limites dessas superfícies. Usando critérios geomorfológicos, três superfícies foram identificadas e topograficamente caracterizadas. Os solos foram amostrados, a intervalos regulares de 25 m, na profundidade de 0,6 a 0,8 m (topo do horizonte B), em uma transeção de 1.700 m perfazendo 109 pontos. Nas amostras, foram analisados: densidade de partículas, granulometria, CTC do solo, CTC da argila, Fe total da argila (ataque por H2SO4) e óxidos de Fe livres (por dissolução seletiva). A fração argila desferrificada foi analisada por difração de raios X. Com base na estratigrafia e variações do relevo local, foram identificadas e diferenciadas, no campo, três superfícies geomórficas. Analisaram-se também o perfil altimétrico e o modelo de elevação digital do terreno. Observou-se que as três diferentes superfícies estão bem relacionadas com os atributos físicos, químicos e mineralógicos dos seus respectivos solos. Na parte inferior desta vertente, superfície mais recente e sobre basalto, em Latossolo Vermelho eutroférrico típico, foram encontradas as maiores variabilidades da declividade, da argila e de Fe. As variações da inclinação do terreno, quando analisadas sistematicamente pelo split moving windows dissimilarity analysis (análise estatística de dissimilaridade, em segmentos móveis), mostraram que este método estatístico pode ser usado para ajudar a localizar os limites entre superfícies geomórficas. As variações dos solos da transeção, e arredores, mostraram-se relacionadas com idade, inclinação do terreno e litologia. O trabalho geomórfico detalhado forneceu importantes informações para subsidiar os trabalhos de levantamento de solos e de pedogênese.

Variabilidade espacial de fatores de erosão em Latossolo Vermelho eutroférrico sob cultivo de cana-de-açúcar

Visando a avaliar a variabilidade espacial de fatores de erosão em Latossolo Vermelho eutroférrico, foram obtidas amostras do solo em intervalos regulares de 50 m, em forma de grid, totalizando 206 pontos de amostragem. Foram coletadas amostras nas profundidades de 0,0-0,2 m para a determinação da composição granulométrica e do conteúdo de matéria orgânica. Os fatores de erosão locais, como erosividade (R), erodibilidade (K), relevo (LS), perda de solo (A), potencial natural de erosão (PNE), risco de erosão (RE) e expectativa de erosão (EE), foram avaliados. A variabilidade do solo medida pelo coeficiente de variação registrou-se média para K, alta para o RE e EE e muito alta para A, LS e PNE. As variáveis estudadas apresentaram estrutura de dependência espacial com grau moderado para as variáveis K, A, PNE e RE, e forte para o LS e EE. Mapas obtidos por krigagem foram apresentados para descrição dos padrões de distribuição dos fatores de erosão na paisagem.

Correlação da produtividade vegetal em função de atributos físico-químicos do solo na integração lavoura-pecuária

Pós-graduação em Agronomia - FEIS

Relationships of Photosynthetic Photon Flux Density, air Temperature and Humidity with Tomato Leaf Diffusive Conductance and Temperature

The objective was to study the leaf temperature (LT) and leaf diffusive vapor conductance (gs) responses to temperature, humidity and incident flux density of photosynthetically active photons (PPFD) of tomato plants grown without water restriction in a plastic greenhouse in Santa Maria, RS, Brazil. The plants were grown in substrate and irrigated daily. The gs was measured using a steady-state null-balance porometer on the abaxial face of the leaves during the daytime. Both leaf surfaces were measured in one day. The PPFD and LT were measured using the porometer. Leaf temperature was determined using an infrared thermometer, and air temperature and humidity were measured using a thermohygrograph. The leaves on the upper layer of the plants had higher gs than the lower layer. The relationship between the gs and PPFD was different for the two layers in the plants. A consistent relationship between the gs and atmospheric water demand was observed only in the lower layer. The LT tended to be lower than the air temperature. The mean value for the gs was 2.88 times higher on the abaxial than adaxial leaf surface.