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.

Multivariate analysis and geostatistics of the fertility of a humic rhodic hapludox under coffee cultivation

The spatial variability of soil and plant properties exerts great influence on the yeld of agricultural crops. This study analyzed the spatial variability of the fertility of a Humic Rhodic Hapludox with Arabic coffee, using principal component analysis, cluster analysis and geostatistics in combination. The experiment was carried out in an area under Coffea arabica L., variety Catucai 20/15 - 479. The soil was sampled at a depth 0.20 m, at 50 points of a sampling grid. The following chemical properties were determined: P, K+, Ca2+, Mg2+, Na+, S, Al3+, pH, H + Al, SB, t, T, V, m, OM, Na saturation index (SSI), remaining phosphorus (P-rem), and micronutrients (Zn, Fe, Mn, Cu and B). The data were analyzed with descriptive statistics, followed by principal component and cluster analyses. Geostatistics were used to check and quantify the degree of spatial dependence of properties, represented by principal components. The principal component analysis allowed a dimensional reduction of the problem, providing interpretable components, with little information loss. Despite the characteristic information loss of principal component analysis, the combination of this technique with geostatistical analysis was efficient for the quantification and determination of the structure of spatial dependence of soil fertility. In general, the availability of soil mineral nutrients was low and the levels of acidity and exchangeable Al were high.

N2O emissions from a cultivated mollisol: optimal time of day for sampling and the role of soil temperature

The correct use of closed field chambers to determine N2O emissions requires defining the time of day that best represents the daily mean N2O flux. A short-term field experiment was carried out on a Mollisol soil, on which annual crops were grown under no-till management in the Pampa Ondulada of Argentina. The N2O emission rates were measured every 3 h for three consecutive days. Fluxes ranged from 62.58 to 145.99 ∝g N-N2O m-2 h-1 (average of five field chambers) and were negatively related (R² = 0.34, p < 0.01) to topsoil temperature (14 - 20 ºC). N2O emission rates measured between 9:00 and 12:00 am presented a high relationship to daily mean N2O flux (R² = 0.87, p < 0.01), showing that, in the study region, sampling in the mornings is preferable for GHG.

Optimization of the sampling scheme for maps of physical and chemical properties estimated by kriging

The sampling scheme is essential in the investigation of the spatial variability of soil properties in Soil Science studies. The high costs of sampling schemes optimized with additional sampling points for each physical and chemical soil property, prevent their use in precision agriculture. The purpose of this study was to obtain an optimal sampling scheme for physical and chemical property sets and investigate its effect on the quality of soil sampling. Soil was sampled on a 42-ha area, with 206 geo-referenced points arranged in a regular grid spaced 50 m from each other, in a depth range of 0.00-0.20 m. In order to obtain an optimal sampling scheme for every physical and chemical property, a sample grid, a medium-scale variogram and the extended Spatial Simulated Annealing (SSA) method were used to minimize kriging variance. The optimization procedure was validated by constructing maps of relative improvement comparing the sample configuration before and after the process. A greater concentration of recommended points in specific areas (NW-SE direction) was observed, which also reflects a greater estimate variance at these locations. The addition of optimal samples, for specific regions, increased the accuracy up to 2 % for chemical and 1 % for physical properties. The use of a sample grid and medium-scale variogram, as previous information for the conception of additional sampling schemes, was very promising to determine the locations of these additional points for all physical and chemical soil properties, enhancing the accuracy of kriging estimates of the physical-chemical properties.

Sampling Position under No-Tillage System Affects the Results of Soil Physical Properties

ABSTRACT Understanding the spatial behavior of soil physical properties under no-tillage system (NT) is required for the adoption and maintenance of a sustainable soil management system. The aims of this study were to quantify soil bulk density (BD), porosity in the soil macropore domain (PORp) and in the soil matrix domain (PORm), air capacity in the soil matrix (ACm), field capacity (FC), and soil water storage capacity (FC/TP) in the row (R), interrow (IR), and intermediate position between R and IR (designated IP) in the 0.0-0.10 and 0.10-0.20 m soil layers under NT; and to verify if these soil properties have systematic variation in sampling positions related to rows and interrows of corn. Soil sampling was carried out in transect perpendicular to the corn rows in which 40 sampling points were selected at each position (R, IR, IP) and in each soil layer, obtaining undisturbed samples to determine the aforementioned soil physical properties. The influence of sampling position on systematic variation of soil physical properties was evaluated by spectral analysis. In the 0.0-0.1 m layer, tilling the crop rows at the time of planting led to differences in BD, PORp, ACm, FC and FC/TP only in the R position. In the R position, the FC/TP ratio was considered close to ideal (0.66), indicating good water and air availability at this sampling position. The R position also showed BD values lower than the critical bulk density that restricts root growth, suggesting good soil physical conditions for seed germination and plant establishment. Spectral analysis indicated that there was systematic variation in soil physical properties evaluated in the 0.0-0.1 m layer, except for PORm. These results indicated that the soil physical properties evaluated in the 0.0-0.1 m layer were associated with soil position in the rows and interrows of corn. Thus, proper assessment of soil physical properties under NT must take into consideration the sampling positions and previous location of crop rows and interrows.

Towards a universal sampling protocol for soil biotas in the humid tropics

This paper reviews the methods for the inventory of below-ground biotas in the humid tropics, to document the (hypothesized) loss of soil biodiversity associated with deforestation and agricultural intensification at forest margins. The biotas were grouped into eight categories, each of which corresponded to a major functional group considered important or essential to soil function. An accurate inventory of soil organisms can assist in ecosystem management and help sustain agricultural production. The advantages and disadvantages of transect-based and grid-based sampling methods are discussed, illustrated by published protocols ranging from the original "TSBF transect", through versions developed for the alternatives to Slash-and-Burn Project (ASB) to the final schemes (with variants) adopted by the Conservation and Sustainable Management of Below-ground Biodiversity Project (CSM-BGBD). Consideration is given to the place and importance of replication in below-ground biological sampling and it is argued that the new sampling protocols are inclusive, i.e. designed to sample all eight biotic groups in the same field exercise; spatially scaled, i.e. provide biodiversity data at site, locality, landscape and regional levels, and link the data to land use and land cover; and statistically robust, as shown by a partial randomization of plot locations for sampling.

Prediction of soil orders with high spatial resolution: response of different classifiers to sampling density

The objective of this work was to evaluate sampling density on the prediction accuracy of soil orders, with high spatial resolution, in a viticultural zone of Serra Gaúcha, Southern Brazil. A digital elevation model (DEM), a cartographic base, a conventional soil map, and the Idrisi software were used. Seven predictor variables were calculated and read along with soil classes in randomly distributed points, with sampling densities of 0.5, 1, 1.5, 2, and 4 points per hectare. Data were used to train a decision tree (Gini) and three artificial neural networks: adaptive resonance theory, fuzzy ARTMap; self‑organizing map, SOM; and multi‑layer perceptron, MLP. Estimated maps were compared with the conventional soil map to calculate omission and commission errors, overall accuracy, and quantity and allocation disagreement. The decision tree was less sensitive to sampling density and had the highest accuracy and consistence. The SOM was the less sensitive and most consistent network. The MLP had a critical minimum and showed high inconsistency, whereas fuzzy ARTMap was more sensitive and less accurate. Results indicate that sampling densities used in conventional soil surveys can serve as a reference to predict soil orders in Serra Gaúcha.

Comparação de duas metodologias de amostragem atmosférica com ferramenta estatística não paramétrica

In atmospheric aerosol sampling, it is inevitable that the air that carries particles is in motion, as a result of both externally driven wind and the sucking action of the sampler itself. High or low air flow sampling speeds may lead to significant particle size bias. The objective of this work is the validation of measurements enabling the comparison of species concentration from both air flow sampling techniques. The presence of several outliers and increase of residuals with concentration becomes obvious, requiring non-parametric methods, recommended for the handling of data which may not be normally distributed. This way, conversion factors are obtained for each of the various species under study using Kendall regression.

Use data from permanent sampling points in growth and yield modeling

This study was carried to evaluate the efficiency of the Bitterlich method in growth and yield modeling of the even-aged Eucalyptus stands. 25 plots were setup in Eucalyptus grandis cropped under a high bole system in the Central Western Region of Minas Gerais, Brazil. The sampling points were setup in the center of each plot. The data of four annual mesurements were colleted and used to adjust the three model types using the age, the site index and the basal area as independent variables. The growths models were fitted for volume and mass of trees. The efficiency of the Bitterlich method was confirmed for generating the data for growth and yield modeling.

Accuracy and efficiency evaluation of point-centered quarter method variations for vegetation sampling in an araucaria forest

In order to verify Point-Centered Quarter Method (PCQM) accuracy and efficiency, using different numbers of individuals by per sampled area, in 28 quarter points in an Araucaria forest, southern Paraná, Brazil. Three variations of the PCQM were used for comparison associated to the number of sampled individual trees: standard PCQM (SD-PCQM), with four sampled individuals by point (one in each quarter), second measured (VAR1-PCQM), with eight sampled individuals by point (two in each quarter), and third measuring (VAR2-PCQM), with 16 sampled individuals by points (four in each quarter). Thirty-one species of trees were recorded by the SD-PCQM method, 48 by VAR1-PCQM and 60 by VAR2-PCQM. The level of exhaustiveness of the vegetation census and diversity index showed an increasing number of individuals considered by quadrant, indicating that VAR2-PCQM was the most accurate and efficient method when compared with VAR1-PCQM and SD-PCQM.

Determination of management zones from normalized and standardized equivalent produtivity maps in the soybean culture

Through the site-specific management, the precision agriculture brings new techniques for the agricultural sector, as well as a larger detailing of the used methods and increase of the global efficiency of the system. The objective of this work was to analyze two techniques for definition of management zones using soybean yield maps, in a productive area handled with localized fertilization and other with conventional fertilization. The sampling area has 1.74 ha, with 128 plots with site-specific fertilization and 128 plots with conventional fertilization. The productivity data were normalized by two techniques (normalized and standardized equivalent productivity), being later classified in management zones. It can be concluded that the two methods of management zones definition had revealed to be efficient, presenting similarities in the data disposal. Due to the fact that the equivalent standardized productivity uses standard score, it contemplates a better statistics justification.

Soil sampling intensity and spatial distribution pattern of soils attributes and corn yield in no-tillage system

Taking into account that the sampling intensity of soil attributes is a determining factor for applying of concepts of precision agriculture, this study aims to determine the spatial distribution pattern of soil attributes and corn yield at four soil sampling intensities and verify how sampling intensity affects cause-effect relationship between soil attributes and corn yield. A 100-referenced point sample grid was imposed on the experimental site. Thus, each sampling cell encompassed an area of 45 m² and was composed of five 10-m long crop rows, where referenced points were considered the center of the cell. Samples were taken from at 0 to 0.1 m and 0.1 to 0.2 m depths. Soil chemical attributes and clay content were evaluated. Sampling intensities were established by initial 100-point sampling, resulting data sets of 100; 75; 50 and 25 points. The data were submitted to descriptive statistical and geostatistics analyses. The best sampling intensity to know the spatial distribution pattern was dependent on the soil attribute being studied. The attributes P and K+ content showed higher spatial variability; while the clay content, Ca2+, Mg2+ and base saturation values (V) showed lesser spatial variability. The spatial distribution pattern of clay content and V at the 100-point sampling were the ones which best explained the spatial distribution pattern of corn yield.

Influence of the sample area in the variability of losses in the mechanical harvesting of soybeans

The mechanical harvesting is an important stage in the production process of soybeans and, in this process; the loss of a significant number of grains is common. Despite the existence of mechanisms to monitor these losses, it is still essential to use sampling methods to quantify them. Assuming that the size of the sample area affects the reliability and variability between samples in quantifying losses, this paper aimed to analyze the variability and feasibility of using different sizes of sample area (1, 2 and 3 m²) in quantifying losses in the mechanical harvesting of soybeans. Were sampled 36 sites and the cutting losses, losses by other mechanisms of the combine and total losses were evaluated, as well as the water content in seeds, straw distribution and crop productivity. Data were subjected to statistical analysis (descriptive statistics and analysis of variance) and Statistical Control Process (SCP). The coefficients of variation were similar for the three frames available. Combine losses showed stable behavior, whereas cutting losses and total losses showed unstable behavior. The frame size did not affect the quantification and variability of losses in the mechanical harvesting of soybeans, thus a frame of 1 m² can be used for determining losses.

Spatial statistics applied to soybean production data from Paraná State for 2003-04 to 2009-10 crop-years

In the current study, we performed a soybean production spatial distribution analysis in Paraná State. Seven crop-year data, from 2003-04 to 2009-10, obtained from the Paraná Department of Agriculture and Supply (SEAB) were used to develop a Boxmap for each crop-year, show soybean production throughout this time interval. Moran's index was used to measure spatial autocorrelation among municipalities at an aggregate level, while LISA index local correlation. For each index, different contiguity matrix and order were used and there was a significance level study. As a result, we have showed spatial relationship among cities regarding the production, which allowed the indication of high and low production clusters. Finally, identifying main soybean-producing cities, what may provide supply chain members with information to strengthen the crop production in Paraná.

Comparison of two methods of tear sampling for protein quantification by Bradford method

The aim of this study was to compare two methods of tear sampling for protein quantification. Tear samples were collected from 29 healthy dogs (58 eyes) using Schirmer tear test (STT) strip and microcapillary tubes. The samples were frozen at -80ºC and analyzed by the Bradford method. Results were analyzed by Student's t test. The average protein concentration and standard deviation from tears collected with microcapillary tube were 4.45mg/mL ±0.35 and 4,52mg/mL ±0.29 for right and left eyes respectively. The average protein concentration and standard deviation from tears collected with Schirmer Tear Test (STT) strip were and 54.5mg/mL ±0.63 and 54.15mg/mL ±0.65 to right and left eyes respectively. Statistically significant differences (p<0.001) were found between the methods. In the conditions in which this study was conducted, the average protein concentration obtained with the Bradford test from tear samples obtained by Schirmer Tear Test (STT) strip showed values higher than those obtained with microcapillary tube. It is important that concentration of tear protein pattern values should be analyzed according the method used to collect tear samples.