Variabilidade espacial de atributos físicos em um Latossolo Vermelho distrófico sob semeadura direta em Selvíria (MS)

Em uma paisagem natural, os solos apresentam ampla variação de propriedades morfológicas, físicas, químicas e mineralógicas, resultante da superposição dos diversos fatores de formação envolvidos. No caso de uma paisagem cultivada, existem outras fontes de heterogeneidade no solo devidas ao manejo exercido pelo homem. O conhecimento dessa variação é importante para o levantamento e classificação dos solos, desenvolvimento de esquemas de amostragem e definições de práticas de manejos. O presente trabalho foi realizado em agosto e setembro de 1999, na Fazenda Experimental de Ensino e Pesquisa da Faculdade de Engenharia de Ilha Solteira/UNESP, localizada em Selvíria (MS), com o objetivo de estudar a variabilidade espacial de alguns atributos físicos de um Latossolo Vermelho distrófico, cultivado no sistema de semeadura direta. A amostragem do solo foi realizada de acordo com um desenho apropriado para a análise geoestatística na forma de uma malha em uma área de 150 m de comprimento, na direção sul, e 30 m de largura, na direção norte, totalizando 103 pontos eqüidistantes de 10 em 10 m e 39 pontos eqüidistantes de 1 m distribuídos aleatoriamente pela malha. De cada ponto definido pela malha, retiraram-se amostras nas profundidades de 0,00-0,05 m e 0,15-0,20 m. Os dados foram avaliados primeiramente por uma análise estatística exploratória, calculando-se a média, distribuição de freqüência, variância, coeficiente de variação, coeficiente de assimetria e coeficiente de curtose. Posteriormente, a dependência espacial foi verificada por meio de semivariogramas. Os atributos microporosidade, porosidade total e densidade do solo seguiram a distribuição normal, enquanto a macroporosidade e a resistência à penetração seguiram a distribuição lognormal. A umidade do solo apresentou uma distribuição mais próxima da lognormal. Os maiores coeficientes de variação foram observados para a macroporosidade e resistência à penetração, tendo as demais variáveis apresentado coeficiente de variação abaixo de 10%. Foi observada uma dependência espacial moderada para todas as variáveis estudadas. O alcance da dependência espacial variou de 8,36 m (umidade do solo) a 58,80 m (resistência à penetração).

Evaluation of occupational exposure: comparison of biological and environmental variabilities using physiologically based toxicokinetic modeling

PURPOSE: Few studies compare the variabilities that characterize environmental (EM) and biological monitoring (BM) data. Indeed, comparing their respective variabilities can help to identify the best strategy for evaluating occupational exposure. The objective of this study is to quantify the biological variability associated with 18 bio-indicators currently used in work environments. METHOD: Intra-individual (BV(intra)), inter-individual (BV(inter)), and total biological variability (BV(total)) were quantified using validated physiologically based toxicokinetic (PBTK) models coupled with Monte Carlo simulations. Two environmental exposure profiles with different levels of variability were considered (GSD of 1.5 and 2.0). RESULTS: PBTK models coupled with Monte Carlo simulations were successfully used to predict the biological variability of biological exposure indicators. The predicted values follow a lognormal distribution, characterized by GSD ranging from 1.1 to 2.3. Our results show that there is a link between biological variability and the half-life of bio-indicators, since BV(intra) and BV(total) both decrease as the biological indicator half-lives increase. BV(intra) is always lower than the variability in the air concentrations. On an individual basis, this means that the variability associated with the measurement of biological indicators is always lower than the variability characterizing airborne levels of contaminants. For a group of workers, BM is less variable than EM for bio-indicators with half-lives longer than 10-15 h. CONCLUSION: The variability data obtained in the present study can be useful in the development of BM strategies for exposure assessment and can be used to calculate the number of samples required for guiding industrial hygienists or medical doctors in decision-making.

Soil water dynamics related to the degree of compaction of two brazilian oxisols under no-tillage

Soil water properties are related to crop growth and environmental aspects and are influenced by the degree of soil compaction. The objective of this study was to determine the water infiltration and hydraulic conductivity of saturated soil under field conditions in terms of the compaction degree of two Oxisols under a no-tillage (NT). Two commercial fields were studied in the state of Rio Grande do Sul, Brazil: one a Haplortox after 14 years under NT; the other a Hapludox after seven years under NT. Maps (50 x 30 m) of the levels of mechanical penetration resistance (PR) were drawn based on the kriging method, differentiating three compaction degrees (CD): high, intermediate and low. In each CD area, the infiltration rate (initial and steady-state) and cumulative water infiltration were measured using concentric rings, with six replications, and the saturated hydraulic conductivity (K(θs)) was determined using the Guelph permeameter. Statistical evaluation was performed based on a randomized design, using the least significant difference (LSD) test and regression analysis. The steady-state infiltration rate was not influenced by the compaction degree, with mean values of 3 and 0.39 cm h-1 in the Haplortox and the Hapludox, respectively. In the Haplortox, saturated soil hydraulic conductivity was 26.76 cm h-1 at a low CD and 9.18 cm h-1 at a high CD, whereas in the Hapludox, this value was 5.16 cm h-1 and 1.19 cm h-1 for the low and high CD, respectively. The compaction degree did not affect the initial and steady-state water infiltration rate, nor the cumulative water infiltration for either soil type, although the values were higher for the Haplortox than the Hapludox.

Spatial variability of soil chemical properties after coffee tree removal

Assessing the spatial variability of soil chemical properties has become an important aspect of soil management strategies with a view to higher crop yields with minimal environmental degradation. This study was carried out at the Centro Experimental of the Instituto Agronomico, in Campinas, São Paulo, Brazil. The aim was to characterize the spatial variability of chemical properties of a Rhodic Hapludox on a recently bulldozer-cleaned area after over 30 years of coffee cultivation. Soil samples were collected in a 20 x 20 m grid with 36 sampling points across a 1 ha area in the layers 0.0-0.2 and 0.2-0.4 m to measure the following chemical properties: pH, organic matter, K+, P, Ca2+, Mg2+, potential acidity, NH4-N, and NO3-N. Descriptive statistics were applied to assess the central tendency and dispersion moments. Geostatistical methods were applied to evaluate and to model the spatial variability of variables by calculating semivariograms and kriging interpolation. Spatial dependence patterns defined by spherical model adjusted semivariograms were made for all cited soil properties. Moderate to strong degrees of spatial dependence were found between 31 and 60 m. It was still possible to map soil spatial variability properties in the layers 0-20 cm and 20-40 cm after plant removal with bulldozers.

Diagnostic techniques applied in geostatistics for agricultural data analysis

The structural modeling of spatial dependence, using a geostatistical approach, is an indispensable tool to determine parameters that define this structure, applied on interpolation of values at unsampled points by kriging techniques. However, the estimation of parameters can be greatly affected by the presence of atypical observations in sampled data. The purpose of this study was to use diagnostic techniques in Gaussian spatial linear models in geostatistics to evaluate the sensitivity of maximum likelihood and restrict maximum likelihood estimators to small perturbations in these data. For this purpose, studies with simulated and experimental data were conducted. Results with simulated data showed that the diagnostic techniques were efficient to identify the perturbation in data. The results with real data indicated that atypical values among the sampled data may have a strong influence on thematic maps, thus changing the spatial dependence structure. The application of diagnostic techniques should be part of any geostatistical analysis, to ensure a better quality of the information from thematic maps.

Spatial and temporal variability of crop yield and some Rhodic Hapludox properties under no-tillage

Soil properties are closely related with crop production and spite of the measures implemented, spatial variation has been repeatedly observed and described. Identifying and describing spatial variations of soil properties and their effects on crop yield can be a powerful decision-making tool in specific land management systems. The objective of this research was to characterize the spatial and temporal variations in crop yield and chemical and physical properties of a Rhodic Hapludox soil under no-tillage. The studied area of 3.42 ha had been cultivated since 1985 under no-tillage crop rotation in summer and winter. Yield and soil property were sampled in a regular 10 x 10 m grid, with 302 sample points. Yields of several crops were analyzed (soybean, maize, triticale, hyacinth bean and castor bean) as well as soil chemical (pH, Soil Organic Matter (SOM), P, Ca2+, Mg2+, H + Al, B, Fe, Mn, Zn, CEC, sum of bases (SB), and base saturation (V %)) and soil physical properties (saturated hydraulic conductivity, texture, density, total porosity, and mechanical penetration resistance). Data were analyzed using geostatistical analysis procedures and maps based on interpolation by kriging. Great variation in crop yields was observed in the years evaluated. The yield values in the Northern region of the study area were high in some years. Crop yields and some physical and soil chemical properties were spatially correlated.

Local influence for spatial analysis of soil physical properties and soybean yield using student's t-distribution

The modeling and estimation of the parameters that define the spatial dependence structure of a regionalized variable by geostatistical methods are fundamental, since these parameters, underlying the kriging of unsampled points, allow the construction of thematic maps. One or more atypical observations in the sample data can affect the estimation of these parameters. Thus, the assessment of the combined influence of these observations by the analysis of Local Influence is essential. The purpose of this paper was to propose local influence analysis methods for the regionalized variable, given that it has n-variate Student's t-distribution, and compare it with the analysis of local influence when the same regionalized variable has n-variate normal distribution. These local influence analysis methods were applied to soil physical properties and soybean yield data of an experiment carried out in a 56.68 ha commercial field in western Paraná, Brazil. Results showed that influential values are efficiently determined with n-variate Student's t-distribution.

Uncertainties in the prediction of spatial variability of soil CO2 emissions and related properties

The soil CO2 emission has high spatial variability because it depends strongly on soil properties. The purpose of this study was to (i) characterize the spatial variability of soil respiration and related properties, (ii) evaluate the accuracy of results of the ordinary kriging method and sequential Gaussian simulation, and (iii) evaluate the uncertainty in predicting the spatial variability of soil CO2 emission and other properties using sequential Gaussian simulations. The study was conducted in a sugarcane area, using a regular sampling grid with 141 points, where soil CO2 emission, soil temperature, air-filled pore space, soil organic matter and soil bulk density were evaluated. All variables showed spatial dependence structure. The soil CO2 emission was positively correlated with organic matter (r = 0.25, p < 0.05) and air-filled pore space (r = 0.27, p < 0.01) and negatively with soil bulk density (r = -0.41, p < 0.01). However, when the estimated spatial values were considered, the air-filled pore space was the variable mainly responsible for the spatial characteristics of soil respiration, with a correlation of 0.26 (p < 0.01). For all variables, individual simulations represented the cumulative distribution functions and variograms better than ordinary kriging and E-type estimates. The greatest uncertainties in predicting soil CO2 emission were associated with areas with the highest estimated values, which produced estimates from 0.18 to 1.85 t CO2 ha-1, according to the different scenarios considered. The knowledge of the uncertainties generated by the different scenarios can be used in inventories of greenhouse gases, to provide conservative estimates of the potential emission of these gases.

Sugarcane trash management assessed by the interaction of yield with soil properties

Currently, sugarcane plays an important global role, particularly with a view to alternative energy sources. Thus, in a sugarcane field of the mill Vale do Paraná S/A Álcool e Açúcar, Rubineia, São Paulo State, managed under two green cane harvest systems (cane trash left on and cane trash removed from the soil), Pearson and spatial correlations between the sugarcane yield (variety RB855035 in the third cut) and soil physical and chemical properties were studied to identify the property best correlated with stalk yield and the best harvest method. For this purpose, two geostatistical grids (121 sampling points on 1.30 ha) were installed on a eutrophic Red Argisol (homogeneous slope of 0.065 m m-1), in 2011, to determine the properties: stalk yield and sugarcane plant population, and soil resistance to penetration, gravimetric moisture, bulk density, and carbon stock, in the layers 0-0.20 and 0.20-0.40 m. The data were analyzed by descriptive, linear correlation and geostatistical analysis. In both treatments, the property stand density was best correlated with sugarcane yield (r = 0.725 in the trash mulching treatment - TM and r = 0.769 in the trash removal treatment - TR). However, in relation to the soil properties, bulk density (0-0.20 m) was best correlated (r = 0.305 in TM, r = 0.211 in TR). Similarly, from the spatial point of view, stand density was the property that best explained the sugarcane yield. However, in the TM treatment the density (0.20-0.40 m) was the only soil property spatially correlated with stalk yield. The carbon stock in the soil of the TM was 11.5 % higher than in the TR treatment. Results of the TM treatment were best, also with regard to soil management and conservation.

Impacts of land leveling on lowland soil physical properties

The practice of land leveling alters the soil surface to create a uniform slope to improve land conditions for the application of all agricultural practices. The aims of this study were to evaluate the impacts of land leveling through the magnitudes, variances and spatial distributions of selected soil physical properties of a lowland area in the State of Rio Grande do Sul, Brazil; the relationships between the magnitude of cuts and/or fills and soil physical properties after the leveling process; and evaluation of the effect of leveling on the spatial distribution of the top of the B horizon in relation to the soil surface. In the 0-0.20 m layer, a 100-point geo-referenced grid covering two taxonomic soil classes was used in assessment of the following soil properties: soil particle density (Pd) and bulk density (Bd); total porosity (Tp), macroporosity (Macro) and microporosity (Micro); available water capacity (AWC); sand, silt, clay, and dispersed clay in water (Disp clay) contents; electrical conductivity (EC); and weighted average diameter of aggregates (WAD). Soil depth to the top of the B horizon was also measured before leveling. The overall effect of leveling on selected soil physical properties was evaluated by paired "t" tests. The effect on the variability of each property was evaluated through the homogeneity of variance test. The thematic maps constructed by kriging or by the inverse of the square of the distances were visually analyzed to evaluate the effect of leveling on the spatial distribution of the properties and of the top of the B horizon in relation to the soil surface. Linear regression models were fitted with the aim of evaluating the relationship between soil properties and the magnitude of cuts and fills. Leveling altered the mean value of several soil properties and the agronomic effect was negative. The mean values of Bd and Disp clay increased and Tp, Macro and Micro, WAD, AWC and EC decreased. Spatial distributions of all soil physical properties changed as a result of leveling and its effect on all soil physical properties occurred in the whole area and not specifically in the cutting or filling areas. In future designs of leveling, we recommend overlaying a cut/fill map on the map of soil depth to the top of the B horizon in order to minimize areas with shallow surface soil after leveling.

Evaluation of statistical and geostatistical models of digital soil properties mapping in tropical mountain regions

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.

SPATIAL CORRELATION BETWEEN PHYSICAL PROPERTIES OF SOIL AND WEEDS IN TWO MANAGEMENT SYSTEMS

The spatial correlation between soil properties and weeds is relevant in agronomic and environmental terms. The analysis of this correlation is crucial for the interpretation of its meaning, for influencing factors such as dispersal mechanisms, seed production and survival, and the range of influence of soil management techniques. This study aimed to evaluate the spatial correlation between the physical properties of soil and weeds in no-tillage (NT) and conventional tillage (CT) systems. The following physical properties of soil and weeds were analyzed: soil bulk density, macroporosity, microporosity, total porosity, aeration capacity of soil matrix, soil water content at field capacity, weed shoot biomass, weed density, Commelina benghalensis density, and Bidens pilosa density. Generally, the ranges of the spatial correlations were higher in NT than in CT. The cross-variograms showed that many variables have a structure of combined spatial variation and can therefore be mapped from one another by co-kriging. This combined variation also allows inferences about the physical and biological meanings of the study variables. Results also showed that soil management systems influence the spatial dependence structure significantly.

ASPECTS OF THE SILVOPASTORAL SYSTEM CORRELATED WITH PROPERTIES OF A TYPIC QUARTZIPSAMMENT (ENTISOL) IN MATO GROSSO DO SUL, BRAZIL

In Brazil, grazing mismanagement may lead to soil and pasture degradation. To impede this process, integrated cropping systems such as silvopasture have been an effective alternative, allied with precision agriculture based on soil mapping for site-specific management. In this study, we aimed to define the soil property that best sheds light on the variability of eucalyptus and forage yield. The experiment was conducted in the 2011/12 crop year in Ribas do Rio Pardo, Mato Grosso do Sul State, Brazil. We analyzed linear and spatial correlations between eucalyptus traits and physical properties of a Typic Quartzipsamment at two depths (0.00-0.10 and 0.10-0.20 m). For that purpose, we set up a geostatistical grid for collection at 72 points. Gravimetric moisture in the 0.00-0.10 m layer is an important index of soil physical quality, showing correlation to eucalyptus circumference at breast height (CBH) in a Typic Quartzipsamment. With an increase in resistance to penetration in the soil surface layer, there is an increase in eucalyptus height and in neutral detergent fiber content in the forage crop. From a spatial point of view, the height of eucalyptus and the neutral detergent fiber of forage can be estimated by co-kriging analysis with soil resistance to penetration. Resistance to penetration values above 2.3 MPa indicated higher yielding sites.

Monitoring network optimisation for spatial data classification using support vector machines

The paper presents a novel method for monitoring network optimisation, based on a recent machine learning technique known as support vector machine. It is problem-oriented in the sense that it directly answers the question of whether the advised spatial location is important for the classification model. The method can be used to increase the accuracy of classification models by taking a small number of additional measurements. Traditionally, network optimisation is performed by means of the analysis of the kriging variances. The comparison of the method with the traditional approach is presented on a real case study with climate data.

A produtividade dos autores sobre a Lei de Lotka

Estuda a produção estratificada dos autores produtores da literatura sobre a Lei de Lotka de 1922 a 2003 e analisa essa produtividade através dos modelos Poisson lognormal e Gauss Poisson inversa generalizada. Para tanto, faz uso dos três tipos de contagem da literatura produzida: contagem direta, contagem completa e contagem fracionada. Os dados da pesquisa são avaliados usando o teste qui-quadrado ao 0.05 nível de significância. Ambos os modelos ajustam-se muito bem à distribuição da literatura produzida, porém a distribuição Poisson Gauss inversa generalizada produz um chi-quadrado menor e prediz melhor o total de autores do que a distribuição Poisson Lognormal.