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.

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.

Geostical analysis of spatial dynamics of Citrus Variegated Chlorossis

The spatial dynamics of Citrus Variegated Chlorosis (CVC) was studied in a five-year old commercial orchard of 'Valencia' sweet orange (Citrus sp.) trees, located in the northern region of the state of São Paulo, Brazil. One thousand trees were assessed in 25 rows of 40 trees, planted at 8 x 5 m spacing. Disease incidence data were taken beginning in March 1994 and ending in January 1996, at intervals of four to five months. Disease aggregation was observed through the dispersion index analysis (Ib), which was calculated by dividing the area into quadrants. CVC spatial dynamics was examined using semivariogram analysis, which revealed that the disease was aggregated in the field forming foci of 10 to 14 m. For each well-fitted model, a kriging map was created to better visualize the distribution of the disease. The spherical model was the best fit for the data in this study. Kriging maps also revealed that the incidence of CVC increased in periods during which the trees underwent vegetative growth, coinciding with greater expected occurrence of insect vectors of the bacterium in the field.

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 air dry bulb temperature and black globe humidity index in a broiler house during the heating phase

The air dry-bulb temperature (t db),as well as the black globe humidity index (BGHI), exert great influence on the development of broiler chickens during their heating phase. Therefore, the aim of this study was to analyze the structure and the magnitude of the t db and BGHI spatial variability, using geostatistics tools such as semivariogram analysis and also producing kriging maps. The experiment was conducted in the west mesoregion of the states of Minas Gerais in 2010, in a commercial broiler house with heating system consisting of two furnaces that heat the air indirectly, in the firsts 14 days of the birds' life. The data were registered at intervals of five minutes in the period from 8 a.m. to 10 a.m. The variables were evaluated by variograms fitted by residual maximum likelihood (REML) testing the Spherical and Exponential models. Kriging maps were generated based on the best model used to fit the variogram. It was possible to characterize the variability of the t db and BGHI, which allowed observing the spatial dependence by using geostatistics techniques. In addition, the use of geostatistics and distribution maps made possible to identify problems in the heating system in regions inside the broiler house that may harm the development of chicks.

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.