Implications of the variability in soil penetration resistance for statistical analysis

Soil penetration resistance (PR) is a measure of soil compaction closely related to soil structure and plant growth. However, the variability in PR hampers the statistical analyses. This study aimed to evaluate the variability of soil PR on the efficiency of parametric and nonparametric analyses in indentifying significant effects of soil compaction and to classify the coefficient of variation of PR into low, medium, high and very high. On six dates, the PR of a typical dystrophic Red Ultisol under continuous no-tillage for 16 years was measured. Three tillage and/or traffic conditions were established with the application of: (i) no chiseling or additional traffic, (ii) additional compaction, and (iii) chiseling. On each date, the nineteen PR data (measured at every 1.5 cm to a depth of 28.5 cm) were grouped in layers with different thickness. In each layer, the treatment effects were evaluated by variance (ANOVA) and Kruskal-Wallis analyses in a completely randomized design, and the coefficients of variation of all analyses were classified (low, intermediate, high and very high). The ANOVA performed better in discriminating the compaction effects, but the rejection rate of null hypothesis decreased from 100 to 80 % when the coefficient of variation increased from 15 to 26 %. The values of 15 and 26 % were the thresholds separating the low/intermediate and the high/very high coefficient variation classes of PR in this Ultisol.

Heavy rainfall equations for Santa Catarina, Brazil

Knowledge of intensity-duration-frequency (IDF) relationships of rainfall events is extremely important to determine the dimensions of surface drainage structures and soil erosion control. The purpose of this study was to obtain IDF equations of 13 rain gauge stations in the state of Santa Catarina in Brazil: Chapecó, Urussanga, Campos Novos, Florianópolis, Lages, Caçador, Itajaí, Itá, Ponte Serrada, Porto União, Videira, Laguna and São Joaquim. The daily rainfall data charts of each station were digitized and then the annual maximum rainfall series were determined for durations ranging from 5 to 1440 min. Based on these, with the Gumbel-Chow distribution, the maximum rainfall was estimated for durations ranging from 5 min to 24 h, considering return periods of 2, 5, 10, 20, 25, 50, and 100 years,. Data agreement with the Gumbel-Chow model was verified by the Kolmogorov-Smirnov test, at 5 % significance level. For each rain gauge station, two IDF equations of rainfall events were adjusted, one for durations from 5 to 120 min and the other from 120 to 1440 min. The results show a high variability in maximum intensity of rainfall events among the studied stations. Highest values of coefficients of variation in the annual maximum series of rainfall were observed for durations of over 600 min at the stations of the coastal region of Santa Catarina.

Modelling of soil penetration resistance for an oxisol under no-tillage

Soil penetration resistance is an important property that affects root growth and elongation and water movement in the soil. Since no-till systems tend to increase organic matter in the soil, the purpose of this study was to evaluate the efficiency with which soil penetration resistance is estimated using a proposed model based on moisture content, density and organic matter content in an Oxisol containing 665, 221 and 114 g kg-1 of clay, silt and sand respectively under annual no-till cropping, located in Londrina, Paraná State, Brazil. Penetration resistance was evaluated at random locations continually from May 2008 to February 2011, using an impact penetrometer to obtain a total of 960 replications. For the measurements, soil was sampled at depths of 0 to 20 cm to determine gravimetric moisture (G), bulk density (D) and organic matter content (M). The penetration resistance curve (PR) was adjusted using two non-linear models (PR = a Db Gc and PR' = a Db Gc Md), where a, b, c and d are coefficients of the adjusted model. It was found that the model that included M was the most efficient for estimating PR, explaining 91 % of PR variability, compared to 82 % of the other model.

Duration-frequency relationships of heavy rainfall in Santa Catarina, Brazil

The purpose of this study was to adjust equations that establish relationships between rainfall events with different duration and data from weather stations in the state of Santa Catarina, Brazil. In this study, the relationships between different duration heavy rainfalls from 13 weather stations of Santa Catarina were analyzed. From series of maximum annual rainfalls, and using the Gumbel-Chow distribution, the maximum rainfall for durations between 5 min and 24 h were estimated considering return periods from 2 to 100 years. The data fit to the Gumbel-Chow model was verified by the Kolmogorov-Smirnov test at 5 % significance. The coefficients of Bell's equation were adjusted to estimate the relationship between rainfall duration t (min) and the return period T (y) in relation to the maximum rainfall with a duration of 1 hour and a 10 year return period. Likewise, the coefficients of Bell's equation were adjusted based on the maximum rainfall with a duration of 1 day and a 10 year return period. The results showed that these relationships are viable to estimate short-duration rainfall events at locations where there are no rainfall records.

Performance of a reciprocal shaker in mechanical dispersion of soil samples for particle-size analysis

The dispersion of the samples in soil particle-size analysis is a fundamental step, which is commonly achieved with a combination of chemical agents and mechanical agitation. The purpose of this study was to evaluate the efficiency of a low-speed reciprocal shaker for the mechanical dispersion of soil samples of different textural classes. The particle size of 61 soil samples was analyzed in four replications, using the pipette method to determine the clay fraction and sieving to determine coarse, fine and total sand fractions. The silt content was obtained by difference. To evaluate the performance, the results of the reciprocal shaker (RSh) were compared with data of the same soil samples available in reports of the Proficiency testing for Soil Analysis Laboratories of the Agronomic Institute of Campinas (Prolab/IAC). The accuracy was analyzed based on the maximum and minimum values defining the confidence intervals for the particle-size fractions of each soil sample. Graphical indicators were also used for data comparison, based on dispersion and linear adjustment. The descriptive statistics indicated predominantly low variability in more than 90 % of the results for sand, medium-textured and clay samples, and for 68 % of the results for heavy clay samples, indicating satisfactory repeatability of measurements with the RSh. Medium variability was frequently associated with silt, followed by the fine sand fraction. The sensitivity analyses indicated an accuracy of 100 % for the three main separates (total sand, silt and clay), in all 52 samples of the textural classes heavy clay, clay and medium. For the nine sand soil samples, the average accuracy was 85.2 %; highest deviations were observed for the silt fraction. In relation to the linear adjustments, the correlation coefficients of 0.93 (silt) or > 0.93 (total sand and clay), as well as the differences between the angular coefficients and the unit < 0.16, indicated a high correlation between the reference data (Prolab/IAC) and results obtained with the RSh. In conclusion, the mechanical dispersion by the reciprocal shaker of soil samples of different textural classes was satisfactory. The results allowed recommending the use of the equipment at low agitation for particle size- analysis. The advantages of this Brazilian apparatus are its low cost, the possibility to simultaneously analyze a great number of samples using ordinary, easily replaceable glass or plastic bottles.

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.

Using numerical classification of profiles based on Vis-NIR spectra to distinguish soils from the Piracicaba Region, Brazil

Considering that information from soil reflectance spectra is underutilized in soil classification, this paper aimed to evaluate the relationship of soil physical, chemical properties and their spectra, to identify spectral patterns for soil classes, evaluate the use of numerical classification of profiles combined with spectral data for soil classification. We studied 20 soil profiles from the municipality of Piracicaba, State of São Paulo, Brazil, which were morphologically described and classified up to the 3rd category level of the Brazilian Soil Classification System (SiBCS). Subsequently, soil samples were collected from pedogenetic horizons and subjected to soil particle size and chemical analyses. Their Vis-NIR spectra were measured, followed by principal component analysis. Pearson's linear correlation coefficients were determined among the four principal components and the following soil properties: pH, organic matter, P, K, Ca, Mg, Al, CEC, base saturation, and Al saturation. We also carried out interpretation of the first three principal components and their relationships with soil classes defined by SiBCS. In addition, numerical classification of the profiles based on the OSACA algorithm was performed using spectral data as a basis. We determined the Normalized Mutual Information (NMI) and Uncertainty Coefficient (U). These coefficients represent the similarity between the numerical classification and the soil classes from SiBCS. Pearson's correlation coefficients were significant for the principal components when compared to sand, clay, Al content and soil color. Visual analysis of the principal component scores showed differences in the spectral behavior of the soil classes, mainly among Argissolos and the others soils. The NMI and U similarity coefficients showed values of 0.74 and 0.64, respectively, suggesting good similarity between the numerical and SiBCS classes. For example, numerical classification correctly distinguished Argissolos from Latossolos and Nitossolos. However, this mathematical technique was not able to distinguish Latossolos from Nitossolos Vermelho férricos, but the Cambissolos were well differentiated from other soil classes. The numerical technique proved to be effective and applicable to the soil classification process.

Variability of soil fertility properties in areas planted to sugarcane in the State of Goias, Brazil

Soil sampling should provide an accurate representation of a given area so that recommendations for amendments of soil acidity, fertilization and soil conservation may be drafted to increase yield and improve the use of inputs. The aim of this study was to evaluate the variability of soil fertility properties of Oxisols in areas planted to sugarcane in the State of Goias, Brazil. Two areas of approximately 8,100 m² each were selected, representing two fields of the Goiasa sugarcane mill in Goiatuba. The sugarcane crop had a row spacing of 1.5 m and subsamples were taken from 49 points in the row and 49 between the row with a Dutch auger at depths of 0.0-0.2 and 0.2-0.4 m, for a total of 196 subsamples for each area. The samples were individually subjected to chemical analyses of soil fertility (pH in CaCl2, potential acidity, organic matter, P, K, Ca and Mg) and particle size analysis. The number of subsamples required to compose a sample within the acceptable ranges of error of 5, 10, 20 and 40 % of each property were computed from the coefficients of variation and the Student t-value for 95 % confidence. The soil properties under analysis exhibited different variabilities: high (P and K), medium (potential acidity, Ca and Mg) and low (pH, organic matter and clay content). Most of the properties analyzed showed an error of less than 20 % for a group of 20 subsamples, except for P and K, which were capable of showing an error greater than 40 % around the mean. The extreme variability in phosphorus, particularly at the depth of 0.2-0.4 m, attributed to banded application of high rates of P fertilizers at planting, places limitations on assessment of its availability due to the high number of subsamples required for a composite sample.

Multivariate analysis applied to the study of the relationship between soil and plant properties in a peach orchard

In the State of Rio Grande do Sul, the municipality of Pelotas is responsible for 90 % of peach production due to its suitable climate and soil conditions. However, there is the need for new studies that aim at improved fruit quality and increased yield. The aim of this study was to evaluate the relationship that exists between soil physical properties and properties in the peach plant in the years 2010 and 2011 by the technique of multivariate canonical correlation. The experiment was conducted in a peach orchard located in the municipality of Morro Redondo, RS, Brazil, where an experimental grid of 101 plants was established. In a trench dug beside each one of the 101 plants, soil samples were collected to determine silt, clay, and sand contents, soil density, total porosity, macroporosity, microporosity, and volumetric water content in the 0.00-0.10 and 0.10-0.20 m layers, as well as the depth of the A horizon. In each plant and in each year, the following properties were assessed: trunk diameter, fruit size and number of fruits per plant, average weight of the fruit per plant, fruit pulp firmness, Brix content, and yield from the orchard. Exploratory analysis of the data was undertaken by descriptive statistics, and the relationships between the physical properties of the soil and of the plant were assessed by canonical correlation analysis. The results showed that the clay and microporosity variables were those that exhibited the highest coefficients of canonical cross-loading with the plant properties in the soil layers assessed, and that the variable of mean weight of the fruit per plant was that which had the highest coefficients of canonical loading within the plant group for the two years assessed.

Impacts of deforestation on water balance components of a watershed on the Brazilian East Coast

The Brazilian East coast was intensely affected by deforestation, which drastically cut back the original biome. The possible impacts of this process on water resources are still unknown. The purpose of this study was an evaluation of the impacts of deforestation on the main water balance components of the Galo creek watershed, in the State of Espírito Santo, on the East coast of Brazil. Considering the real conditions of the watershed, the SWAT model was calibrated with data from 1997 to 2000 and validated for the period between 2001 and 2003. The calibration and validation processes were evaluated by the Nash-Sutcliffe efficiency coefficient and by the statistical parameters (determination coefficient, slope coefficient and F test) of the regression model adjusted for estimated and measured flow data. After calibration and validation of the model, new simulations were carried out for three different land use scenarios: a scenario in compliance with the law (C1), assuming the preservation of PPAs (permanent preservation areas); an optimistic scenario (C2), which considers the watershed to be almost entirely covered by native vegetation; and a pessimistic scenario (C3), in which the watershed would be almost entirely covered by pasture. The scenarios C1, C2 and C3 represent a soil cover of native forest of 76, 97 and 0 %, respectively. The results were compared with the simulation, considering the real scenario (C0) with 54 % forest cover. The Nash-Sutcliffe coefficients were 0.65 and 0.70 for calibration and validation, respectively, indicating satisfactory results in the flow simulation. A mean reduction of 10 % of the native forest cover would cause a mean annual increase of approximately 11.5 mm in total runoff at the watershed outlet. Reforestation would ensure minimum flows in the dry period and regulate the maximum flow of the main watercourse of the watershed.

Public-domain software for root image analysis

In the search for high efficiency in root studies, computational systems have been developed to analyze digital images. ImageJ and Safira are public-domain systems that may be used for image analysis of washed roots. However, differences in root properties measured using ImageJ and Safira are supposed. This study compared values of root length and surface area obtained with public-domain systems with values obtained by a reference method. Root samples were collected in a banana plantation in an area of a shallower Typic Carbonatic Haplic Cambisol (CXk), and an area of a deeper Typic Haplic Ta Eutrophic Cambisol (CXve), at six depths in five replications. Root images were digitized and the systems ImageJ and Safira used to determine root length and surface area. The line-intersect method modified by Tennant was used as reference; values of root length and surface area measured with the different systems were analyzed by Pearson's correlation coefficient and compared by the confidence interval and t-test. Both systems ImageJ and Safira had positive correlation coefficients with the reference method for root length and surface area data in CXk and CXve. The correlation coefficient ranged from 0.54 to 0.80, with lowest value observed for ImageJ in the measurement of surface area of roots sampled in CXve. The IC (95 %) revealed that root length measurements with Safira did not differ from that with the reference method in CXk (-77.3 to 244.0 mm). Regarding surface area measurements, Safira did not differ from the reference method for samples collected in CXk (-530.6 to 565.8 mm²) as well as in CXve (-4231 to 612.1 mm²). However, measurements with ImageJ were different from those obtained by the reference method, underestimating length and surface area in samples collected in CXk and CXve. Both ImageJ and Safira allow an identification of increases or decreases in root length and surface area. However, Safira results for root length and surface area are closer to the results obtained with the reference method.

Transfer of cadmium and barium from soil to crops grown in tropical soils

Phytotoxicity and transfer of potentially toxic elements, such as cadmium (Cd) or barium (Ba), depend on the availability of these elements in soils and on the plant species exposed to them. With this study, we aimed to evaluate the effect of Cd and Ba application rates on yields of pea (Pisum sativum L.), sorghum (Sorghum bicolor L.), soybean (Glycine max L.), and maize (Zea mays L.) grown under greenhouse conditions in an Oxisol and an Entisol with contrasting physical and chemical properties, and to correlate the amount taken up by plants with extractants commonly used in routine soil analysis, along with transfer coefficients (Bioconcentration Factor and Transfer Factor) in different parts of the plants. Plants were harvested at flowering stage and measured for yield and Cd or Ba concentrations in leaves, stems, and roots. The amount of Cd accumulated in the plants was satisfactorily evaluated by both DTPA and Mehlich-3 (M-3). Mehlich-3 did not relate to Ba accumulated in plants, suggesting it should not be used to predict Ba availability. The transfer coefficients were specific to soils and plants and are therefore not recommended for direct use in risk assessment models without taking soil properties and group of plants into account.

Analysis of experiments in square lattice with emphasis on variance components. i. Individual analysis

This paper focused on four alternatives of analysis of experiments in square lattice as far as the estimation of variance components and some genetic parameters are concerned: 1) intra-block analysis with adjusted treatment and blocks within unadjusted repetitions; 2) lattice analysis as complete randomized blocks; 3) intrablock analysis with unadjusted treatment and blocks within adjusted repetitions; 4) lattice analysis as complete randomized blocks, by utilizing the adjusted means of treatments, obtained from the analysis with recovery of interblock information, having as mean square of the error the mean effective variance of this same analysis with recovery of inter-block information. For the four alternatives of analysis, the estimators and estimates were obtained for the variance components and heritability coefficients. The classification of material was also studied. The present study suggests that for each experiment and depending of the objectives of the analysis, one should observe which alternative of analysis is preferable, mainly in cases where a negative estimate is obtained for the variance component due to effects of blocks within adjusted repetitions.