Semivariogram models for estimating fig fly population density throughout the year

The objective of this work was to select semivariogram models to estimate the population density of fig fly (Zaprionus indianus; Diptera: Drosophilidae) throughout the year, using ordinary kriging. Nineteen monitoring sites were demarcated in an area of 8,200 m2, cropped with six fruit tree species: persimmon, citrus, fig, guava, apple, and peach. During a 24 month period, 106 weekly evaluations were done in these sites. The average number of adult fig flies captured weekly per trap, during each month, was subjected to the circular, spherical, pentaspherical, exponential, Gaussian, rational quadratic, hole effect, K-Bessel, J-Bessel, and stable semivariogram models, using ordinary kriging interpolation. The models with the best fit were selected by cross-validation. Each data set (months) has a particular spatial dependence structure, which makes it necessary to define specific models of semivariograms in order to enhance the adjustment to the experimental semivariogram. Therefore, it was not possible to determine a standard semivariogram model; instead, six theoretical models were selected: circular, Gaussian, hole effect, K-Bessel, J-Bessel, and stable.

Amostragem e variabilidade espacial de atributos químicos do solo em área de vegetação natural em regeneração

The objectives of this work was to estimate the number of soil subsamples considering the classical statistics and geostatistics and determine the spatial variability of soil fertility attributes of an Ultisol, with clay texture, in an area of regenerating natural vegetation in Alegre - ES. Soil samples were collected in a depth of 0.0-0.2 m, at the crossing points of a regular grid, comprising a total of 64 points located at 10 m-intervals. The area presented low fertility soil. Considering a variation of 5% around the mean in the classic statistics, it is necessary a larger number of samples in relation to geostatistics. All the chemical attributes showed moderate to high spatial dependence, except for the effective cation exchange capacity (CECe), which showed pure nugget effect. The spherical semivariogram model gave the best fit to the data. Isoline maps allowed visualizing the differentiated spatial distribution of the contents of soil chemical attributes.

Variabilidade espacial de atributos físicos de um argissolo vermelho-amarelo sob pastagem e vegetação secundária em regeneração natural

This work was carried out with the objective of studying the spatial variability of the physical attributes of a Red-Yellow Ultisol under pasture and secondary vegetation in natural regeneration. Two areas were chosen in a hillside, with the soil sampling to the depth of 0-0.2 m, with the georeferenced points in a regular grid of 10x10 m, totalizing 64 points. In each point it was evaluated the total volume of porosity, macroporosity, microporosity, bulk density, soil penetration resistance and soil water content. The studied attributes in the pasture area present indicator of soil compaction for the animals' traffic, with moderate and strong structure of spatial dependence, except for the macroporosity and penetration resistance. In the area of secondary vegetation (VN) only the macroporosity does not present spatial dependence. The total volume of porosity and the bulk density present the same spatial standard in the area under pasture.

Variabilidade espacial da textura de um latossolo sob cultivo de citros

Relief influences soil texture variability, since it contributes to the time of exposition of the materials to weathering factors. Our work was carried out in the city of Gavião Peixoto (SP), with the objective of characterizing the spatial variability of texture of a dystrophic Red Latosol cultivated with citrus. The hillside was divided into three segments: top, stocking lean and inferior lean. Soil samples were collected in a grid with regular intervals of 50 m, at the depths of 0.0-0.2 m and 0.6-0.8 m, comprising a total of 332 points in an area of 83.5 ha. The data were submitted to descriptive and geostatistics analyses (semivariogram modeling and kriging maps). The spatial behavior of the texture of oxisols is directly related to the relief forms in this study, which controls the direction of surface and subsurface water flows. The concept of homogeneity of clay distribution in the Oxisol profile is a piece of information that can be adjusted by knowing the spatial pattern of this distribution in different relief forms.

Distribuição espacial de atributos de solo e de plantas na cultura do mamoeiro no norte capixaba

No Brasil, as lavouras de mamoeiro das planícies dos tabuleiros costeiros são as que melhor desenvolvem e aplicam tecnologias para a produção de mamão no mundo. O objetivo foi aplicar a estatística clássica e a geoestatística no mapeamento e na correlação da variabilidade espacial de atributos químicos e físicos de solo e de plantas de mamoeiro (Carica papaya L.) de uma lavoura comercial do norte capixaba cultivada em um Argissolo típico dos tabuleiros costeiros. O solo de textura arenosa de caráter coeso foi preparado convencionalmente e cultivado com mamoeiro variedade Golden THB. Após a sexagem, procederam-se as amostragens de solo, amostrado na projeção da copa (0-0,20 e 0,20-0,40 m) para a determinação dos atributos químicos e físicos, e de atributos biométricos das plantas em uma área de 1,2 ha (114 x 110 m) totalizando 129 pontos amostrais georreferenciados. Ao nono mês após o transplantio, registrou-se a altura da colheita dos primeiros frutos, o número e a massa dos frutos colhidos para estimativa da produtividade, amostrando três plantas por ponto amostral durante três meses. Os dados foram submetidos à análise estatística descritiva e à correlação de Pearson. A dependência espacial das variáveis foi analisada através da ferramenta geoestatística, com obtenção de semivariograma e os mapas de distribuição das variáveis. A maior parte dos atributos de solo e de plantas de mamoeiro apresenta dependência espacial e é mapeada adequadamente. Há correlação de dependência vertical para densidade do solo, argila, silte, resistência do solo à penetração na linha de plantio e na rua e volume total de poros. Dos atributos químicos não ocorre este comportamento apenas para K, Al e Sat K. As frações areia e argila foram os principais atributos a constituírem correlação com os demais. Há poucas correlações dos atributos do solo com os atributos biométricos e a produtividade do mamoeiro. Ocorre correlação positiva entre a produtividade inicial do mamoeiro com características biométricas ideias para as plantas de mamoeiro. A fertilidade e o preparo do solo são expressivos para o desenvolvimento do mamoeiro e para a variabilidade espacial dos atributos avaliados.

Geostatistical Mapping of Outfall Plume Dispersion Data gathered with an autonomous underwater vehicle

The main purpose of this study was to examine the applicability of geostatistical modeling to obtain valuable information for assessing the environmental impact of sewage outfall discharges. The data set used was obtained in a monitoring campaign to S. Jacinto outfall, located off the Portuguese west coast near Aveiro region, using an AUV. The Matheron’s classical estimator was used the compute the experimental semivariogram which was fitted to three theoretical models: spherical, exponential and gaussian. The cross-validation procedure suggested the best semivariogram model and ordinary kriging was used to obtain the predictions of salinity at unknown locations. The generated map shows clearly the plume dispersion in the studied area, indicating that the effluent does not reach the near by beaches. Our study suggests that an optimal design for the AUV sampling trajectory from a geostatistical prediction point of view, can help to compute more precise predictions and hence to quantify more accurately dilution. Moreover, since accurate measurements of plume’s dilution are rare, these studies might be very helpful in the future for validation of dispersion models.

Scaling of semivariograms and the kriging estimation of field-measured properties

Two methods were evaluated for scaling a set of semivariograms into a unified function for kriging estimation of field-measured properties. Scaling is performed using sample variances and sills of individual semivariograms as scale factors. Theoretical developments show that kriging weights are independent of the scaling factor which appears simply as a constant multiplying both sides of the kriging equations. The scaling techniques were applied to four sets of semivariograms representing spatial scales of 30 x 30 m to 600 x 900 km. Experimental semivariograms in each set successfully coalesced into a single curve by variances and sills of individual semivariograms. To evaluate the scaling techniques, kriged estimates derived from scaled semivariogram models were compared with those derived from unscaled models. Differences in kriged estimates of the order of 5% were found for the cases in which the scaling technique was not successful in coalescing the individual semivariograms, which also means that the spatial variability of these properties is different. The proposed scaling techniques enhance interpretation of semivariograms when a variety of measurements are made at the same location. They also reduce computational times for kriging estimations because kriging weights only need to be calculated for one variable. Weights remain unchanged for all other variables in the data set whose semivariograms are scaled.

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.

Spatial variability of physical attributes of an alfisol under different hillslope curvatures

The influence of relief forms has been studied by several authors and explains the variability in the soil attributes of a landscape. Soil physical attributes depend on relief forms, and their assessment is important in mechanized agricultural systems, such as of sugarcane. This study aimed to characterize the spatial variability in the physical soil attributes and their relationship to the hillslope curvatures in an Alfisol developed from sandstone and growing sugarcane. Grids of 100 x 100 m were delimited in a convex and a concave area. The grids had a regular spacing of 10 x 10 m, and the crossing points of this spacing determined a total of 121 georeferenced sampling points. Samples were collected to determine the physical attributes related to soil aggregates, porosity, bulk density, resistance to penetration and moisture within the 0-0.2 and 0.2-0.4 m depth. Statistical analyses, geostatistics and Student's t-tests were performed with the means of the areas. All attributes, except aggregates > 2 mm in the 0-0.2 m depth and macroporosity at both depths, showed significant differences between the hillslope curvatures. The convex area showed the highest values of the mean weighted diameter, mean geometric diameter, aggregates > 2 mm, 1-2 mm aggregates, total porosity and moisture and lower values of bulk density and resistance to penetration in both depth compared to the concave area. The number of soil attributes with greater spatial variability was higher in the concave area.

Pedotransfer functions related to spatial variability of water retention attributes for lowland soils

The estimation of non available soil variables through the knowledge of other related measured variables can be achieved through pedotransfer functions (PTF) mainly saving time and reducing cost. Great differences among soils, however, can yield non desirable results when applying this method. This study discusses the application of developed PTFs by several authors using a variety of soils of different characteristics, to evaluate soil water contents of two Brazilian lowland soils. Comparisons are made between PTF evaluated data and field measured data, using statistical and geostatistical tools, like mean error, root mean square error, semivariogram, cross-validation, and regression coefficient. The eight tested PTFs to evaluate gravimetric soil water contents (Ug) at the tensions of 33 kPa and 1,500 kPa presented a tendency to overestimate Ug 33 kPa and underestimate Ug1,500 kPa. The PTFs were ranked according to their performance and also with respect to their potential in describing the structure of the spatial variability of the set of measured values. Although none of the PTFs have changed the distribution pattern of the data, all resulted in mean and variance statistically different from those observed for all measured values. The PTFs that presented the best predictive values of Ug33 kPa and Ug1,500 kPa were not the same that had the best performance to reproduce the structure of spatial variability of these variables.

Spatial patterns of preconsolidation pressure and soil moisture along transects in two directions under coffee

Information on the spatial structure of soil physical and structural properties is needed to evaluate the soil quality. The purpose of this study was to investigate the spatial behavior of preconsolidation pressure and soil moisture in six transects, three selected along and three across coffee rows, at three different sites under different tillage management systems. The study was carried out on a farm, in Patrocinio, state of Minas Gerais, in the Southeast of Brazil (18 º 59 ' 15 '' S; 46 º 56 ' 47 '' W; 934 m asl). The soil type is a typic dystrophic Red Latosol (Acrustox) and consists of 780 g kg-1 clay; 110 g kg-1 silt and 110 g kg-1 sand, with an average slope of 3 %. Undisturbed soil cores were sampled at a depth of 0.10-0.13 m, at three different points within the coffee plantation: (a) from under the wheel track, where equipment used in farm operations passes; (b) in - between tracks and (c) under the coffee canopy. Six linear transects were established in the experimental area: three transects along and three across the coffee rows. This way, 161 samples were collected in the transect across the coffee rows, from the three locations, while 117 samples were collected in the direction along the row. The shortest sampling distance in the transect across the row was 4 m, and 0.5 m for the transect along the row. No clear patterns of the preconsolidation pressure values were observed in the 200 m transect. The results of the semivariograms for both variables indicated a high nugget value and short range for the studied parameters of all transects. A cyclic pattern of the parameters was observed for the across-rows transect. An inverse relationship between preconsolidation pressure and soil moisture was clearly observed in the samples from under the track, in both directions.

Mineralogy of the clay fraction of alfisols in two slope curvatures: IV - spatial correlation with physical properties

Although the influence of clay mineralogy on soil physical properties has been widely studied, spatial relationships between these features in Alfisols have rarely been examined. The purpose of this work was to relate the clay minerals and physical properties of an Alfisol of sandstone origin in two slope curvatures. The crystallographic properties such as mean crystallite size (MCS) and width at half height (WHH) of hematite, goethite, kaolinite and gibbsite; contents of hematite and goethite; aluminium substitution (AS) and specific surface area (SSA) of hematite and goethite; the goethite/(goethite+hematite) and kaolinite/(kaolinite+gibbsite) ratios; and the citrate/bicarbonate/dithionite extractable Fe (Fe d) were correlated with the soil physical properties through Pearson correlation coefficients and cross-semivariograms. The correlations found between aluminium substitution in goethite and the soil physical properties suggest that the degree of crystallinity of this mineral influences soil properties used as soil quality indicators. Thus, goethite with a high aluminium substitution resulted in large aggregate sizes and a high porosity, and also in a low bulk density and soil penetration resistance. The presence of highly crystalline gibbsite resulted in a high density and micropore content, as well as in smaller aggregates. Interpretation of the cross-semivariogram and classification of landscape compartments in terms of the spatial dependence pattern for the relief-dependent physical and mineralogical properties of the soil proved an effective supplementary method for assessing Pearson correlations between the soil physical and mineralogical properties.

Water storage variability in a vineyard soil in the southern highlands of Santa Catarina state

In the subtropical regions of southern Brazil, rainfall distribution is uneven, which results in temporal variability of soil water storage. For grapes, water is generally available in excess and water deficiency occurs only occasionally. Furthermore, on the Southern Plateau of Santa Catarina, there are differences in soil properties, which results in high spatial variability. These two factors affect the composition of wine grapes. Spatio-temporal analyses are therefore useful in the selection of cultural practices as well as of adequate soils for vineyards. In this way, well-suited areas can produce grapes with a more appropriate composition for the production of quality wines. The aim of this study was to evaluate the spatio-temporal variability of water storage in a Cambisol during the growth cycle of a Cabernet Sauvignon vineyard and its relation to selected soil properties. The experimental area consisted of a commercial 8-year-old vineyard in São Joaquim, Santa Catarina, Brazil. A sampling grid with five rows and seven points per row, spaced 12 m apart, was outlined on an area of 3,456 m². Soil samples were collected with an auger at these points, 0.30 m away from the grapevines, in the 0.00-0.30 m layer, to determine gravimetric soil moisture. Measurements were taken once a week from December 2008 to April 2009, and every two weeks from December 2009 to March 2010. In December 2008, undisturbed soil samples were collected to determine bulk density, macro- and microporosity, and disturbed samples were used to quantify particle size distribution and organic carbon content. Results were subjected to descriptive analysis and semivariogram analysis, calculating the mean relative difference and the Pearson correlation. The average water storage in a Cambisol under grapevine on ridges had variable spatial dependence, i.e., the lower the average water storage, the higher the range of spatial dependence. Water storage had a stable spatial pattern during the trial period, indicating that the points with lower water storage or points with higher water storage during a certain period maintain these conditions throughout the experimental period. The relative difference is a simple method to identify positions that represent the average soil water storage more adequately at any time for a given area.

USE OF SCALED SEMIVARIOGRAMS IN THE PLANNING SAMPLE OF SOIL PHYSICAL PROPERTIES IN SOUTHERN AMAZONAS, BRAZIL

There is a great lack of information from soil surveys in the southern part of the State of Amazonas, Brazil. The use of tools such as geostatistics may improve environmental planning, use and management. In this study, we aimed to use scaled semivariograms in sample design of soil physical properties of some environments in Amazonas. We selected five areas located in the south of the state of Amazonas, Brazil, with varied soil uses, such as forest, archaeological dark earth (ADE), pasture, sugarcane cropping, and agroforestry. Regular mesh grids were set up in these areas with 64 sample points spaced at 10 m from each other. At these points, we determined the particle size composition, soil resistance to penetration, moisture, soil bulk density and particle density, macroporosity, microporosity, total porosity, and aggregate stability in water at a depth of 0.00-0.20 m. Descriptive and geostatistical analyses were performed. The sample density requirements were lower in the pasture area but higher in the forest. We concluded that managed-environments had differences in their soil physical properties compared to the natural forest; notably, the soil in the ADE environment is physically improved in relation to the others. The physical properties evaluated showed a structure of spatial dependence with a slight variability of the forest compared to the others. The use of the range parameter of the semivariogram analysis proved to be effective in determining an ideal sample density.

Estudio de la distribución espacial de Cydia pomonella (L.) y Pandemis heparana (Denis & Schiffermüller) en Torregrossa (Lleida) mediante métodos geoestadísticos

El estudio de la distribución espacial de una especie por me'todos geoestadísticos se realiza mediante el conocimiento de la función semivariograma. Después de calcular el semivariograma se procede a la estimación de la variable regionalizada en cualquier punto de la zona de estudio. Esta estimación se realiza mediante técnicas de interpolación lineal llamadas «krigeado», en honor a Krige y Matheron, fundadores de la geoestadística. El «krigeado» se basa en la minimización de la varianza del error en cada punto de estudio, previamente localizado en el espacio por sus coordenadas de situación. Cydia pomonella (L.) y Pandemis heparana (Denis & Schiffermüller) son dos tortrícidos plaga del manzano y del peral. La estimación de sus poblaciones se realiza mediante trampas de feromona y es posible disponer de una amplia base de datos. El objetivo de este trabajo fue analizar la idoneidad de los métodos geoestadísticos para el estudio de poblaciones de insectos y aplicarlas al caso concreto de C. pomonella y P. heparana. Se utilizaron las capturas en 55 estaciones con trampas de feromonas (difusor de origen Wageningen) colocadas en parcelas comerciales de manzano y peral en 1996 y 1997 en el término municipal de Torregrossa (Lleida). La idoneidad de los métodos geoestadísticos quedó demostrada por el hecho de que la variable número acumulado de machos por trampa fue regionalizable. Una vez calculadas las funciones semivariograma para cada especie y año, se han dibujado los mapas de distribución mediante el uso de isolíneas. En el futuro, se plantea la posibilidad de ampliar la zona de estudio a toda la zona frutera de Lleida y analizar la influencia de variables independientes (climáticas ...) sobre la distribución espacial mediante métodos de «co-krigeado».