Distribuição espacial e raio de agregação de cancro cítrico definidos por geoestatística

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

Geoestatística na avaliação dos atributos físicos em latossolo sob floresta nativa e pastagem na Região de Manicoré, Amazonas

The use of geostatistical techniques allows detection of the existence of dependence and the spatial distribution of soil properties, thus constituting an important tool in the analysis and detailed description of the behavior of soil physical properties. The aim of the present study was to use geostatistics in assessment of physical properties in a Latossolo (Oxisol) dystrophic under native forest and pasture in the Amazon region of Manicore. Grids with of 70 x 70 m were established in native forest and pasture, and points were marked in these grids spaced at every 10 m, for a total of 64 points. These points were then georeferenced and in each one, soil samples (128) were collected at the depths of 0.00-0.20 and 0.40-0.60 m for a survey of their physical properties. These grids are parallel at a distance of 100 m from one another. The following determinations were made: texture, bulk density and particle density, macroporosity, microporosity, total porosity and aggregate stability in water. After tabulating the data, descriptive statistical analysis and geostatistical analysis were performed. The pasture had a slight variation in its physical properties in relation to native forest, with a high coefficient of variation and weak spatial dependence. The scaled semivariograms were able to satisfactorily reproduce the spatial behavior of the properties in the same pattern as the individual semivariograms, and the use of the parameter range of the semivariogram was efficient for determining the optimal sampling density for the environments under study. The geostatistical results indicate that the removal of native forest for establishing pasture altered the natural variability of the physical properties.

Hot spots, hot moments, and spatio-temporal controls on soil CO2 efflux in a water-limited ecosystem

Soil CO2 efflux is the primary source of CO2 emissions from terrestrial ecosystems to the atmosphere. The rates of this flux vary in time and space producing hot moments (sudden temporal high fluxes) and hot spots (spatially defined high fluxes), but these high reaction rates are rarely studied in conjunction with each other. We studied temporal and spatial variation of soil CO2 efflux in a water-limited Mediterranean ecosystem in Baja California, Mexico. Soil CO2 efflux increased 522% during a hot moment after rewetting of soils following dry summer months. Monthly precipitation was the primary driver of the seasonal trend of soil CO2 efflux (including the hot moment) and through changes in soil volumetric water content (VWC) it influenced the relationship between CO2 efflux and soil temperature. Geostatistical analyses showed that the spatial dependence of soil CO2 efflux changed between two contrasting seasons (dry and wet). During the dry season high soil VWC was associated with high soil CO2 efflux, and during the wet season the emergence of a hot spot of soil CO2 efflux was associated with higher root biomass and leaf area index. These results suggest that sampling designs should accommodate for changes in spatial dependence of measured variables. The spatio-temporal relationships identified in this study are arguably different from temperate ecosystems where the majority of soil CO2 efflux research has been done. This study provides evidence of the complexity of the mechanisms controlling the spatio-temporal variability of soil CO2 efflux in water-limited ecosystems. (C) 2014 Elsevier Ltd. All rights reserved.

Distribuição espacial de Enneothrips flavens, Moulton 1941 (Thysanoptera: Thripidae), na cultura do amendoim, utilizando geoestatística

Pós-graduação em Agronomia (Entomologia Agrícola) - FCAV

Re-analysis of pluvial precipitation in southern Brazil

This study examined precipitation in southern Brazil based on a data set provided by the Brazilian National Water Agency, covering the period from 1976 to 2010. Data were homogenized using the R software and the Climatol subroutine, which allow completing missing data. Isohyets were drawn using the Geostatistics software to obtain a semivariogram for each analysis. There was a remarkable interannual variability in this region, with positive anomalies in the warm phase (El Nino) and negative anomalies in the cold phase (La Nina) of ENSO. Also, the responses of this variability were not uniform in the entire region, since there was variability from year to year and from event to event.

Spatial intra-annual variability of precipitation based on geostatistics. a case study for the Paraiba do Sul Basin, southeastern Brazil

The emphasis in this research is to evaluate the spatial distribution of the precipitation using a geostatistics approach. Seasonal time scales records considering DJF, MAM, JJA e SON periods performed the analysis. Procedures to evaluate the variogram selection and to produce kriging maps were performed in a GIS environment (ArcGIS®). The results showed that kriging method was very suitable to detect both large changes in the whole area as those local small and subtle changes. Kriging demonstrated be a powerful statistical interpolation method that might be very useful in regions with great complexity in climatology and geomorphology.

SPATIAL VARIABILITY OF THE RAINFALL EROSIVE POTENTIAL IN THE STATE OF MATO GROSSO DO SUL, BRAZIL

Information about rainfall erosivity is important during soil and water conservation planning. Thus, the spatial variability of rainfall erosivity of the state Mato Grosso do Sul was analyzed using ordinary kriging interpolation. For this, three pluviograph stations were used to obtain the regression equations between the erosivity index and the rainfall coefficient EI30. The equations obtained were applied to 109 pluviometric stations, resulting in EI30 values. These values were analyzed from geostatistical technique, which can be divided into: descriptive statistics, adjust to semivariogram, cross-validation process and implementation of ordinary kriging to generate the erosivity map. Highest erosivity values were found in central and northeast regions of the State, while the lowest values were observed in the southern region. In addition, high annual precipitation values not necessarily produce higher erosivity values.

Spatial variability of soil fertility and its relationship with seed physiological potential in a soybean production area

The correlation of soil fertility x seed physiological potential is very important in the area of seed technology but results published with that theme are contradictory. For this reason, this study to evaluate the correlations between soil chemical properties and physiological potential of soybean seeds. On georeferenced points, both soil and seeds were sampled for analysis of soil fertility and seed physiological potential. Data were assessed by the following analyses: descriptive statistics; Pearson's linear correlation; and geostatistics. The adjusted parameters of the semivariograms were used to produce maps of spatial distribution for each variable. Organic matter content, Mn and Cu showed significant effects on seed germination. Most variables studied presented moderate to high spatial dependence. Germination and accelerated aging of seeds, and P, Ca, Mg, Mn, Cu and Zn showed a better fit to spherical semivariogram: organic matter, pH and K had a better fit to Gaussian model; and V% and Fe showed a better fit to the linear model. The values for range of spatial dependence varied from 89.9 m for P until 651.4 m for Fe. These values should be considered when new samples are collected for assessing soil fertility in this production area.

Spatial Patterns of Intraseasonal Variability of Chlorophyll and Sea Surface Temperature in the California Current

Six years of daily satellite data are used to quantify and map intraseasonal variability of chlorophyll and sea surface temperature (SST) in the California Current. We define intraseasonal variability as temporal variation remaining after removal of interannual variability and stationary seasonal cycles. Semivariograms are used to quantify the temporal structure of residual time series. Empirical orthogonal function (EOF) analyses of semivariograms calculated across the region isolate dominant scales and corresponding spatial patterns of intraseasonal variability. The mode 1 EOFs for both chlorophyll and SST semivariograms indicate a dominant timescale of similar to 60 days. Spatial amplitudes and patterns of intraseasonal variance derived from mode 1 suggest dominant forcing of intraseasonal variability through distortion of large scale chlorophyll and SST gradients by mesoscale circulation. Intraseasonal SST variance is greatest off southern Baja and along southern Oregon and northern California. Chlorophyll variance is greatest over the shelf and slope, with elevated values closely confined to the Baja shelf and extending farthest from shore off California and the Pacific Northwest. Intraseasonal contributions to total SST variability are strongest near upwelling centers off southern Oregon and northern California, where seasonal contributions are weak. Intraseasonal variability accounts for the majority of total chlorophyll variance in most inshore areas save for southern Baja, where seasonal cycles dominate. Contributions of higher EOF modes to semivariogram structure indicate the degree to which intraseasonal variability is shifted to shorter timescales in certain areas. Comparisons of satellite-derived SST semivariograms to those calculated from co-located and concurrent buoy SST time series show similar features.

The Discontinuous Nature of Kriging Interpolation for Digital Terrain Modeling

Kriging is a widely employed method for interpolating and estimating elevations from digital elevation data. Its place of prominence is due to its elegant theoretical foundation and its convenient practical implementation. From an interpolation point of view, kriging is equivalent to a thin-plate spline and is one species among the many in the genus of weighted inverse distance methods, albeit with attractive properties. However, from a statistical point of view, kriging is a best linear unbiased estimator and, consequently, has a place of distinction among all spatial estimators because any other linear estimator that performs as well as kriging (in the least squares sense) must be equivalent to kriging, assuming that the parameters of the semivariogram are known. Therefore, kriging is often held to be the gold standard of digital terrain model elevation estimation. However, I prove that, when used with local support, kriging creates discontinuous digital terrain models, which is to say, surfaces with “rips” and “tears” throughout them. This result is general; it is true for ordinary kriging, kriging with a trend, and other forms. A U.S. Geological Survey (USGS) digital elevation model was analyzed to characterize the distribution of the discontinuities. I show that the magnitude of the discontinuity does not depend on surface gradient but is strongly dependent on the size of the kriging neighborhood.

Soil fragmentation study applying different tillage systems

Runoff generation depends on rainfall, infiltration, interception, and surface depressional storage. Surface depressional storage depends on surface microtopography, usually quantified trough soil surface roughness (SSR). SSR is subject to spatial and temporal changes that create a high variability. In an agricultural environment, tillage operations produce abrupt changes in roughness. Subsequent rainfall gradually decreases roughness. Beside it, local variation in soil properties and hydrology cause its SSR to vary spatially at different scales. The methods commonly used to measure it involve collecting point elevations in regular grids using laser profilers or scanners, digital close range stereo-photogrammetry and terrestrial laser scanning or LIDAR systems. In this case, a laser-scanning instrument was used to obtain representative digital elevation models (DEMs) at a grid resolution of 7.2x7.2mm that cover an area of 0.9x0.9m. The DEMs were obtained from two study sites with different soils. The first study site was an experimental field on which five conventional tillage methods were applied. The second study site was a large olive orchard with trees planted at 7.5x5.0m and bare soils between rows. Here, three tillage treatments were applied. In this work we have evaluated the spatial variability of SSR at several scales studying differences in height calculated from points separated by incremental distances h were raised to power values q (from 0 to 4 in steps of 0.1). The q = 2 data were studied as a semivariogram model. The logarithm of average differences plotted vs. log h were characterized by their slope, ?(q). Structure functions [?(q) vs. q] were fitted showing that data had nonlinear structure functions typical of multiscale phenomena. Comparisson of the two types of soil in their respective structure functions are shown.

Spatial variability patterns of some Vertisol properties at a field scale using standardized data

Spatial variability of Vertisol properties is relevant for identifying those zones with physical degradation. In this sense, one has to face the problem of identifying the origin and distribution of spatial variability patterns. The objectives of the present work were (i) to quantify the spatial structure of different physical properties collected from a Vertisol, (ii) to search for potential correlations between different spatial patterns and (iii) to identify relevant components through multivariate spatial analysis. The study was conducted on a Vertisol (Typic Hapludert) dedicated to sugarcane (Saccharum officinarum L.) production during the last sixty years. We used six soil properties collected from a squared grid (225 points) (penetrometer resistance (PR), total porosity, fragmentation dimension (Df), vertical electrical conductivity (ECv), horizontal electrical conductivity (ECh) and soil water content (WC)). All the original data sets were z-transformed before geostatistical analysis. Three different types of semivariogram models were necessary for fitting individual experimental semivariograms. This suggests the different natures of spatial variability patterns. Soil water content rendered the largest nugget effect (C0 = 0.933) while soil total porosity showed the largest range of spatial correlation (A = 43.92 m). The bivariate geostatistical analysis also rendered significant cross-semivariance between different paired soil properties. However, four different semivariogram models were required in that case. This indicates an underlying co-regionalization between different soil properties, which is of interest for delineating management zones within sugarcane fields. Cross-semivariograms showed larger correlation ranges than individual, univariate, semivariograms (A ≥ 29 m). All the findings were supported by multivariate spatial analysis, which showed the influence of soil tillage operations, harvesting machinery and irrigation water distribution on the status of the investigated area.

Análisis de la distribución espacial de los parámetros acústicos en recintos

El principal objetivo de este estudio es la evaluación de la distribución espacial de los parámetros acústicos en un recinto a través de la aplicación de técnicas geoestadísticas como el método Kriging. Mediante el uso de la herramienta de análisis espacial ArcMap, perteneciente a la plataforma ArcGIS, se ha analizado el comportamiento acústico del Salón de Actos común a la Escuela Técnica Superior de Ingeniería y Sistemas de Telecomunicación (ETSIST) y la Escuela Técnica Superior de Ingeniería de Sistemas Informáticos (ETSISI), ambas situadas en el Campus Sur de la Universidad Politécnica de Madrid. Se han realizado mediciones in-situ del recinto no ocupado utilizando la herramienta de medición DIRAC Room Acoustics y el método de la respuesta impulsiva integrada, extrayéndose los parámetros acústicos de tiempo de reverberación (RT), tiempo de reverberación inicial (EDT), fuerza sonora relativa (Grel), claridad (C80), tiempo central (Ts), definición (D50) e Índice de Transmisión Rápida de la Palabra (RASTI). Se ha analizado la adecuación de los valores observados de Grel, C80 y Ts al modelo teórico de Barron y estimado, mediante el método Kriging Ordinario, el conjunto de parámetros medidos en el recinto, obteniéndose los semivariograma y mapas de estimación correspondientes. Además, se ha evaluado la calidad de la estimación en base a un número de puntos de medición reducido. A la vista de los resultados obtenidos, en general, el método Kriging puede considerarse un buen interpolador de los parámetros acústicos en un recinto, observándose que los parámetros que evalúan relaciones energéticas, especialmente la fuerza sonora relativa (Grel) proporcionan mejores estimaciones en comparación con aquellos relacionados con la reverberación y la inteligibilidad del habla. El coeficiente de determinación (R2) constituye una medida útil para evaluar la precisión de la estimación. Además, la entropía de los datos observados puede ser un buen indicador a priori de la precisión de la estimación. Asimismo, se ha demostrado que, basándose en un reducido número de puntos de medición, es posible obtener una estimación precisa de los parámetros acústicos de fuerza sonora relativa (Grel) y tiempo central (Ts). ABSTRACT. This project aims to evaluate the feasibility of using geostatistical techniques such as Kriging on the analysis of the spatial distribution of the acoustic parameters in rooms. The acoustic behaviour of the Assembly Hall of the ETSIST and ETSISI (Universidad Politécnica de Madrid) is investigated using ArcMap, which is the main component of ArcGIS suite of geospatial processing programs. For this purpose, in-situ acoustic measurements are carried out in the unoccupied room using DIRAC Room Acoustics software. The following acoustic parameters are measured by means of the integrated impulse response method for further examination: Reverberation Time (RT), Early Decay Time (EDT), Relative Strength (Grel), Clarity (C80), Centre Time (Ts), Definition (D50) and Rapid Speech Transmission Index (RASTI). Goodness-of-fit of measured Grel, C80 and Ts values to Barron’s theory is determined and Ordinary Kriging is applied to all the measured parameters in order to calculate the semivariogram and prediction surfaces. The prediction performance is also analysed when significantly fewer receiver positions are used for the prediction. The experimental results obtained lead to conclude that Kriging can be successfully applied to room acoustics. Energy\based acoustic parameters can be estimated with higher accuracy compared to those related to reverberation and speech intelligibility. Coefficient of determination (R2) is a reliable statistic for assessing the prediction accuracy, for which measured data entropy can also be a good a priori indicator. Furthermore, based on fewer receiver positions, it is demonstrated that accurate predictions of Grel and Ts can be achieved.

Linking riparian vegetation spatial structure in Australian tropical savannas to ecosystem health indicators: semi-variogram analysis of high spatial resolution satellite imagery

Government agencies responsible for riparian environments are assessing the combined utility of field survey and remote sensing for mapping and monitoring indicators of riparian zone health. The objective of this work was to determine if the structural attributes of savanna riparian zones in northern Australia can be detected from commercially available remotely sensed image data. Two QuickBird images and coincident field data covering sections of the Daly River and the South Alligator River - Barramundie Creek in the Northern Territory were used. Semi-variograms were calculated to determine the characteristic spatial scales of riparian zone features, both vegetative and landform. Interpretation of semi-variograms showed that structural dimensions of riparian environments could be detected and estimated from the QuickBird image data. The results also show that selecting the correct spatial resolution and spectral bands is essential to maximize the accuracy of mapping spatial characteristics of savanna riparian features. The distribution of foliage projective cover of riparian vegetation affected spectral reflectance variations in individual spectral bands differently. Pan-sharpened image data enabled small-scale information extraction (< 6 m) on riparian zone structural parameters. The semi-variogram analysis results provide the basis for an inversion approach using high spatial resolution satellite image data to map indicators of savanna riparian zone health.

Efeitos diretos e indiretos de parâmetros dos semivariogramas referentes aos atributos químicos do solo sobre a produtividade de grãos de soja

In order to add value to soybens crops, and hence the marketing, medium and large producers have been using precision agriculture techniques (PA), as the Remote Sensing, Geographic Information Systems (GIS) and positioning satellite, to assist the management of crops. Thus, given the economic relevance of that culture to the southwest of Paraná State and Brazil, scientific studies to increase their productivity and profitability are of main importance. The objective of this study was to evaluate the correlation between the chemical soil properties and soybean yield for each estimated parameter of semivariogram (range, nugget and level effect), and the deployment of these correlations in direct and indirect effects, aiming to improve the mapping process of spatial variability of soil chemical properties for use in PA. The hypothesis is that not all attributes of soil used to estimate the semivariogram parameters has a direct effect on productivity, and that even in groups of plants within a larger area it is possible to estimate the parameters of the semivariograms. The experiment was conducted in a commercial area of 19.7 ha, located in the city of Pato Branco - PR, central geographic coordinates 26º 11 '35 "South, 52 43' 05" West longitude, and average altitude of 780 m. The area is planted with soybeans for over 30 years, currently being adopted to cultivate Brasmax Target RR - Don Mario 5.9i, with row spacing of 0.50 m and 13 plants m-1, totaling 260,000 plants ha-1. For georeferencing of the area of study and sampling points was used a couple of topographic ProMarkTM3 receptors, making a relative positioning to obtain the georeferenced coordinates. To collect data (chemical analyzes of soil and crop yield) were sampled 10 blocks in the experimental area, each with an area of 20 m2 (20 meters long x 1 meter wide) containing two spaced adjacent rows of 0.5 m. Each block was divided into 20 portions of 1 m2, and from each were collected four subsamples at a distance of 0.5 m in relation to the lines of blocks, making up a sample depth for 0-10 cm a sample to 10-20 cm for each plot, totaling 200 samples for each depth. The soybean crop was performed on the blocks depending on maturity, and in each block was considered a bundle at each meter. In the data analysis, it was performed a diagnosis of multicollinearity, and subsequently a path analysis of the main variables according to the explanatory variables (range of chemical attributes: pH, K, P, Ca, etc.). The results obtained by the path analysis of the parameters of the semivariogram of soil chemical properties, indicated that only the Fe, Mg, Mn, organic matter (OM), P and Saturation by bases (SB) exerted direct and indirect effects on soybean productivity, although they have not presented spatial variability, indicating that the distribution of blocks in the area was unable to identify the spatial dependence of these elements, making it impossible to draw up maps of the chemical attributes for use in PA.