Sampling strategies for tropical forest nutrient cycling studies: a case study in São Paulo, Brazil

The precise sampling of soil, biological or micro climatic attributes in tropical forests, which are characterized by a high diversity of species and complex spatial variability, is a difficult task. We found few basic studies to guide sampling procedures. The objective of this study was to define a sampling strategy and data analysis for some parameters frequently used in nutrient cycling studies, i. e., litter amount, total nutrient amounts in litter and its composition (Ca, Mg, Κ, Ν and P), and soil attributes at three depths (organic matter, Ρ content, cation exchange capacity and base saturation). A natural remnant forest in the West of São Paulo State (Brazil) was selected as study area and samples were collected in July, 1989. The total amount of litter and its total nutrient amounts had a high spatial independent variance. Conversely, the variance of litter composition was lower and the spatial dependency was peculiar to each nutrient. The sampling strategy for the estimation of litter amounts and the amount of nutrient in litter should be different than the sampling strategy for nutrient composition. For the estimation of litter amounts and the amount of nutrients in litter (related to quantity) a large number of randomly distributed determinations are needed. Otherwise, for the estimation of litter nutrient composition (related to quality) a smaller amount of spatially located samples should be analyzed. The determination of sampling for soil attributes differed according to the depth. Overall, surface samples (0-5 cm) showed high short distance spatial dependent variance, whereas, subsurface samples exhibited spatial dependency in longer distances. Short transects with sampling interval of 5-10 m are recommended for surface sampling. Subsurface samples must also be spatially located, but with transects or grids with longer distances between sampling points over the entire area. Composite soil samples would not provide a complete understanding of the relation between soil properties and surface dynamic processes or landscape aspects. Precise distribution of Ρ was difficult to estimate.

Variability of apparently homogeneous soilscapes in São Paulo state, Brazil: I. spatial analysis

The spatial variability of strongly weathered soils under sugarcane and soybean/wheat rotation was quantitatively assessed on 33 fields in two regions in São Paulo State, Brazil: Araras (15 fields with sugarcane) and Assis (11 fields with sugarcane and seven fields with soybean/wheat rotation). Statistical methods used were: nested analysis of variance (for 11 fields), semivariance analysis and analysis of variance within and between fields. Spatial levels from 50 m to several km were analyzed. Results are discussed with reference to a previously published study carried out in the surroundings of Passo Fundo (RS). Similar variability patterns were found for clay content, organic C content and cation exchange capacity. The fields studied are quite homogeneous with respect to these relatively stable soil characteristics. Spatial variability of other characteristics (resin extractable P, pH, base- and Al-saturation and also soil colour), varies with region and, or land use management. Soil management for sugarcane seems to have induced modifications to greater depths than for soybean/wheat rotation. Surface layers of soils under soybean/wheat present relatively little variation, apparently as a result of very intensive soil management. The major part of within-field variation occurs at short distances (< 50 m) in all study areas. Hence, little extra information would be gained by increasing sampling density from, say, 1/km² to 1/50 m². For many purposes, the soils in the study regions can be mapped with the same observation density, but residual variance will not be the same in all areas. Bulk sampling may help to reveal spatial patterns between 50 and 1.000 m.

Variability of apparently homogeneous soilscapes in São Paulo state, Brazil: II. quality of soil maps

The quality of semi-detailed (scale 1:100.000) soil maps and the utility of a taxonomically based legend were assessed by studying 33 apparently homogeneous fields with strongly weathered soils in two regions in São Paulo State: Araras and Assis. An independent data set of 395 auger sites was used to determine purity of soil mapping units and analysis of variance within and between mapping units and soil classification units. Twenty three soil profiles were studied in detail. The studied soil maps have a high purity for some legend criteria, such as B horizon type (> 90%) and soil texture class (> 80%). The purity for the "trophic character" (eutrophic, dystrophic, allic) was only 55% in Assis. It was 88% in Araras, where many soil units had been mapped as associations. In both regions, the base status of clay-textured soils was generally better than suggested by the maps. Analysis of variance showed that mapping was successful for "durable" soil characteristics such as clay content (> 80% of variance explained) and cation exchange capacity (≥ 50% of variance explained) of 0-20 and 60-80 cm layers. For soil characteristics that are easily modified by management, such as base saturation of the 0-20 cm layer, the maps had explained very little (< 15%) of the total variance in the study areas. Intermediate results were obtained for base saturation of the 60-80 cm layer (56% in Assis; 42% in Araras). Variance explained by taxonomic groupings that formed the basis for the legend of the soil maps was similar to, often even smaller than, variance explained by mapping units. The conclusion is that map boundaries have been very carefully located, but descriptions of mapping units could be improved. In future mappings, this could possibly be done at low cost by (a) bulk sampling to remove short range variation and enhance visualization of spatial patterns at distances > 100 m; (b) taking advantage of correlations between easily measured soil characteristics and chemical soil properties and, (c) unbending the link between legend criteria and a taxonomic system. The maps are well suited to obtain an impression of land suitability for high-input farming. Additional field work and data on former land use/management are necessary for the evaluation of chemical properties of surface horizons.

Rice grown in nutrient solution with doses of manganese and silicon

Although silicon is not recognized as a nutrient, it may benefit rice plants and may alleviate the Mn toxicity in some plant species. The dry matter yield (root, leaf, sheaths and leaf blade) and plant architecture (angle of leaf insertion and leaf arc) were evaluated in rice plants grown in nutrient solutions with three Mn doses, with and without Si addition. The treatments were arranged in a 2 x 3 factorial [with and without (2 mmol L-1) Si; three Mn doses (0.5; 2.5 and 10 µmol L-1)], in a randomized block design with 4 replications. The experimental unit was a 4 L plastic vase with 4 rice (Metica-1 cultivar) plants. Thirty nine days after keeping the seedlings in the nutrient solution the plant dry matter yield was determined; the angle of leaf insertion in the sheath and the leaf arc were measured; and the Si and Mn concentrations in roots, sheaths and leaves were determined. The analysis of variance (F test at 5 and 1 % levels) and the regression analysis (for testing plant response to Mn with the Si treatments) were performed. The Si added to the nutrient solution increased the dry matter yield of roots, sheaths and leaf blades and also decreased the angle of leaf blade insertion into the sheath and the foliar arc in the rice plant. Additionally, it ameliorated the rice plant architecture which allowed an increase in the dry matter yield. Similarly, the addition of Mn to the solution improved the architecture of the rice plants with gain in dry matter yield. As Si was added to the nutrient solution, the concentration of Mn in leaves decreased and in roots increased thus alleviating the toxic effects of Mn on the plants.

Enzymatic activity measured by microcalorimetry in soil amended with organic residues

Enzymatic activity is an important property for soil quality evaluation. Two sequences of experiments were carried out in order to evaluate the enzymatic activity in a soil (Rhodic Eutrudox) amended with cattle manure, earthworm casts, or sewage sludges from the municipalities of Barueri and Franca. The activity of commercial enzymes was measured by microcalorimetry in the same soil samples after sterilization. In the first experiment, the enzyme activities of cellulase, protease, and urease were determined in the soil samples during a three month period. In the second sequence of experiments, the thermal effect of the commercial enzymes cellulase, protease, and urease on sterilized soil samples under the same tretaments was monitored for a period of 46 days. The experimental design was randomized and arranged as factorial scheme in five treatments x seven samplings with five replications. The treatment effects were statistically evaluated by one-way analysis of variance. Tukey´s test was used to compare means at p < 0.05. The presence of different sources of organic residues increased the enzymatic activity in the sampling period. Cattle manure induced the highest enzymatic activity, followed by municipal sewage sludge, whereas earthworm casts induced the lowest activity, but differed from control treatment. The thermal effect on the enzyme activity of commercial cellulase, protease, and urease showed a variety of time peaks. These values probably oscillated due to soil physical-chemical factors affecting the enzyme activity on the residues.

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.

Relationships between microbial activity and soil physical and chemical properties in native and reforested Araucaria angustifolia forests in the state of São Paulo, Brazil

Araucaria angustifolia (Bert.) O. Kuntze is the main component of the Mixed Ombrophilous forest and, in the State of São Paulo, it is associated with a high diversity of soil organisms, essential for the maintenance of soil quality, making the conservation of this ecosystem a major and pressing challenge. The objective of this study was to identify the physical and chemical properties that are most closely correlated with dehydrogenase enzyme activity, basal respiration and microbial biomass under native (NF) and replanted (RF) Araucaria angustifolia forests in three regions of the state of São Paulo, in winter and summer. The main differentiating factors between the areas were also determined. Each forest was represented by three true replications; at each site, from around the araucaria trees, 15 soil samples (0-20 cm) were collected to evaluate the soil physical, chemical and microbiological properties. At the same points, forest litter was sampled to assess mass and chemical properties. The following microbiological properties were evaluated: microbial biomass carbon (MBC), basal respiration (CO2-C), metabolic quotient (Q: CO2), dehydrogenase enzyme activity (DHA) as well as the physical properties (moisture, bulk density, macroporosity and total porosity), soil chemical properties [pH, organic carbon (org-C), P, Ca, K, Mg, Al, H+Al], litter dry mass, and C, N and S contents. The data were subjected to analysis of variance (TWO-WAY: ANOVA). A Canonical Discriminant Analysis (CDA) and a Canonical Correlation Analysis (CCA) were also performed. In the soil under NF, the values of K, P, soil macroporosity, and litter dry mass were higher and Q: CO2 and DHA lower, regardless of the sampling period, and DHA was lower in winter. In the RF areas, the levels of moisture, porosity and Q: CO2 were higher in both sampling periods, and DHA was higher in winter. The MBC was only higher under NF in the summer, while the litter contents of C, N and S were greater in winter. In winter, CCA showed a high correlation of DHA with CO2-C, pH and H+Al, while in the summer org-C, moisture, Mg, pH and litter C were more associated with DHA and CO2-C. The CDA indicated H+Al, available P, total porosity, litter S content, and soil moisture as the most discriminating variables between NF and RF, but moisture was the most relevant, in both seasons and CO2-C only in winter. The combined analysis of CCA and CDA showed that the contribution of the microbiological variables to a differentiation of the areas was small at both samplings, which may indicate that the period after reforestation was long enough to allow an almost complete recovery of the microbial activity.

Arsenic toxicity in Acacia mangium willd. and mimosa Caesalpiniaefolia benth. seedlings

Acacia mangium and Mimosa caesalpiniaefolia are fast-growing woody fabaceous species that might be suitable for phytoremediation of arsenic (As)-contaminated sites. To date, few studies on their tolerance to As toxicity have been published. Therefore, this study assessed As toxicity symptoms in A. mangium and M. caesalpiniaefolia seedlings under As stress in a greenhouse. Seedlings of Acacia mangium and M. caesalpiniaefolia were grown for 120 d in an Oxisol-sand mixture with 0, 50, 100, 200, and 400 mg kg-1 As, in four replications in four randomized blocks. The plants were assessed for visible toxicity symptoms, dry matter production, shoot/root ratio, root anatomy and As uptake. Analyses of variance and regression showed that the growth of A. mangium and M. caesalpiniaefolia was severely hindered by As, with a reduction in dry matter production of more than 80 % at the highest As rate. The root/shoot ratio increased with increasing As rates. At a rate of 400 mg kg-1 As, whitish chlorosis appeared on Mimosa caesalpiniaefolia seedlings. The root anatomy of both species was altered, resulting in cell collapse, death of root buds and accumulation of phenolic compounds. Arsenic concentration was several times greater in roots than in shoots, with more than 150 and 350 mg kg-1 in M. caesalpiniaefolia and A. mangium roots, respectively. These species could be suitable for phytostabilization of As-contaminated sites, but growth-stimulating measures should be used.

Variability of root traits in common bean genotypes at different levels of phosphorus supply and ontogenetic stages

Selection of common bean (Phaseolus vulgaris L.) cultivars with enhanced root growth would be a strategy for increasing P uptake and grain yield in tropical soils, but the strong plasticity of root traits may compromise their inclusion in breeding programs. The aim of this study was to evaluate the magnitude of the genotypic variability of root traits in common bean plants at two ontogenetic stages and two soil P levels. Twenty-four common bean genotypes, comprising the four growth habits that exist in the species and two wild genotypes, were grown in 4 kg pots at two levels of applied P (20 and 80 mg kg-1) and harvested at the stages of pod setting and early pod filling. Root area and root length were measured by digital image analysis. Significant genotype × P level and genotype × harvest interactions in analysis of variance indicate that the genotypic variation of root traits depended on soil nutrient availability and the stage at which evaluation was made. Genotypes differed for taproot mass, basal and lateral root mass, root area and root length at both P levels and growth stages; differences in specific root area and length were small. Genotypes with growth habits II (upright indeterminate) and III (prostrate indeterminate) showed better adaptation to limited P supply than genotypes of groups I (determinate) and IV (indeterminate climbing). Between the two harvests, genotypes of groups II and III increased the mass of basal and lateral roots by 40 and 50 %, respectively, whereas genotypes of groups I and IV by only 7 and 19 %. Values of the genotypic coefficient of determination, which estimates the proportion of phenotypic variance resulting from genetic effects, were higher at early pod filling than at pod setting. Correlations between shoot mass and root mass, which could indicate indirect selection of root systems via aboveground biomass, were higher at early pod filling than at pod setting. The results indicate that selection for root traits in common bean genotypes should preferentially be performed at the early pod-filling stage.

Can Soil Penetration Resistance and Bulk Density Be Determined in a Single Undisturbed Sample?

Soil quality indicators such as penetration resistance (PR) and bulk density (BD) are traditionally determined in a single undisturbed soil sample. The aim of this study was to assess the effect of PR measurements of undisturbed samples on the determination of BD in the same sample of two soils differing in clay contents. To this end, samples were collected from the 0.00-0.10 and 0.10-0.20 m layers of two soils of clayey and very clayey texture. Volumetric rings were used to collect a total of 120 undisturbed soil samples from each soil layer that were divided into two subsets containing 60 units each. One sample set, designated “perforated samples”, was used to determine PR and BD in the same undisturbed sample; the other, named “intact samples”, was used to determine BD only. Bulk density values for perforated and intact samples were compared by analysis of variance, using a completely randomized experimental design. Means were compared by the t-test at 5 %. The BD values for the clayey soil were similar in perforated and intact samples from the two layers. However, BD of the very clayey soil was lower in the perforated than in the intact samples at both depths. Therefore, PR and BD in clayey soils can be accurately determined in the same undisturbed sample whereas in very clayey soils, different samples are required for this purpose.

Diallel analysis for grain yield and mineral absorption rate of soybeans grown in acid brazilian savannah soil

High available aluminium and low levels of calcium below the ploughed zone of the soil are limiting factors for agricultural sustainability in the Brazilian Cerrados (Savannahs). The mineral stresses compound with dry spells effect by preventing deep root growth of cultivated plants and causes yield instability. The mode of inheritance for grain yield and mineral absorption ratio of a diallel cross in soybeans [Glycine max (L.) Merrill] grown in high and low Al areas was identified. Differences among the genotypes for grain yield were more evident in the high Al, by grouping tolerant and non-tolerant genotypes for their respective arrays in the hybrids. A large proportion of genetic variance was additive for grain yield and mineral absorption ratio in both environments. High heritability values suggest that soybeans can be improved by crosses among Al-tolerant genotypes, using modified pedigree, early generation and recurrent selection schemes.

Phenotypic stability and genetic gains in six-year girth growth of Hevea clones

Rubber tree [Hevea brasiliensis (Willd. ex Adr. de Juss.) Müell. Arg.] budgrafts of seven clones were evaluated on five contrasting sites in the plateau region of the São Paulo State, Brazil. The objective of this work was to study the phenotypic stability for girth growth. The experimental design was a randomized block design with three replications and seven treatments. Analysis of variance of girth at six-year plant growth indicated a highly significant clone x site interaction. Only linear sites and clone x site components of clone x year interaction were significant, indicating that the performance of clones over sites for this trait could be predicted. The clones GT 1 and PB 235 showed the greatest stability in relation to girth growth, with foreseen responses to change, introduced in the sites. The clones PB 235 and IAN 873 showed significative difference in relation to regression coefficient, representing clones with specific adaptability on favorable and unfavorable sites respectively. The clone GT 1 became the most promissory one in the study of stability and adaptability even showing low girth growth. Expected genetic gains from planting sites, along with estimates of clonal variance and repeatability of clonal means are generally greatest or close to the greatest when selection is done at the same site.

Analysis of variance of primary data on plant growth analysis

Plant growth analysis presents difficulties related to statistical comparison of growth rates, and the analysis of variance of primary data could guide the interpretation of results. The objective of this work was to evaluate the analysis of variance of data from distinct harvests of an experiment, focusing especially on the homogeneity of variances and the choice of an adequate ANOVA model. Data from five experiments covering different crops and growth conditions were used. From the total number of variables, 19% were originally homoscedastic, 60% became homoscedastic after logarithmic transformation, and 21% remained heteroscedastic after transformation. Data transformation did not affect the F test in one experiment, whereas in the other experiments transformation modified the F test usually reducing the number of significant effects. Even when transformation has not altered the F test, mean comparisons led to divergent interpretations. The mixed ANOVA model, considering harvest as a random effect, reduced the number of significant effects of every factor which had the F test modified by this model. Examples illustrated that analysis of variance of primary variables provides a tool for identifying significant differences in growth rates. The analysis of variance imposes restrictions to experimental design thereby eliminating some advantages of the functional growth analysis.

Genotype and environment interaction on soybean yield in Mato Grosso State, Brazil

The objective of this work was to investigate the genotype-environment interaction in Mato Grosso State, MT. The relative importance of locations, years, sowing dates and cultivars and their interactions was analyzed from data collected in regional yield trials performed in a randomized complete block design with four replications, from 1994-1995 through 1999-2000, in nine locations and two sowing dates. Individual and pooled analyses of variance over years and locations were performed. Complementary analyses of variances partitioned MT State in two main and five smaller regions, respectively: North and South of Cuiabá; and MT-South-A (Pedra Preta area), MT-South-B (Rondonópolis and Itiquira), MT-East (Primavera do Leste and Campo Verde), MT-Central (Nova Mutum, Lucas do Rio Verde and Sorriso) and MT-Parecis (Campo Novo dos Parecis and Sapezal). Locations are relatively more important than years for yield testing soybeans in the MT State, therefore, investment should be made in increasing locations rather than years to improve experimental precision. Partitioning the MT State into regions has little impact on soybean yield testing results and, consequently, on the efficiency of the soybean breeding program in the State. Breeding genotypes with broad adaptation for the MT State is an efficient strategy for cultivar development.

Sample size for full-sib family evaluation in sugarcane

The objective of this study was to determine the minimum number of plants per plot that must be sampled in experiments with sugarcane (Saccharum officinarum) full-sib families in order to provide an effective estimation of genetic and phenotypic parameters of yield-related traits. The data were collected in a randomized complete block design with 18 sugarcane full-sib families and 6 replicates, with 20 plants per plot. The sample size was determined using resampling techniques with replacement, followed by an estimation of genetic and phenotypic parameters. Sample-size estimates varied according to the evaluated parameter and trait. The resampling method permits an efficient comparison of the sample-size effects on the estimation of genetic and phenotypic parameters. A sample of 16 plants per plot, or 96 individuals per family, was sufficient to obtain good estimates for all traits considered of all the characters evaluated. However, for Brix, if sample separation by trait were possible, ten plants per plot would give an efficient estimate for most of the characters evaluated.