Extraction of soil water by plants: development and validation of a model

A quantitative model of water movement within the immediate vicinity of an individual root is developed and results of an experiment to validate the model are presented. The model is based on the assumption that the amount of water transpired by a plant in a certain period is replaced by an equal volume entering its root system during the same time. The model is based on the Darcy-Buckingham equation to calculate the soil water matric potential at any distance from a plant root as a function of parameters related to crop, soil and atmospheric conditions. The model output is compared against measurements of soil water depletion by rice roots monitored using γ-beam attenuation in a greenhouse of the Escola Superior de Agricultura "Luiz de Queiroz"/Universidade de São Paulo(ESALQ/USP) in Piracicaba, State of São Paulo, Brazil, in 1993. The experimental results are in agreement with the output from the model. Model simulations show that a single plant root is able to withdraw water from more than 0.1 m away within a few days. We therefore can assume that root distribution is a less important factor for soil water extraction efficiency.

Tropical terrestrial model ecosystems for evaluation of soil fauna and leaf litter quality effects on litter consumption, soil microbial biomass and plant growth

The aim of this work was to evaluate whether terrestrial model ecosystems (TMEs) are a useful tool for the study of the effects of litter quality, soil invertebrates and mineral fertilizer on litter decomposition and plant growth under controlled conditions in the tropics. Forty-eight intact soil cores (17.5-cm diameter, 30-cm length) were taken out from an abandoned rubber plantation on Ferralsol soil (Latossolo Amarelo) in Central Amazonia, Brazil, and kept at 28ºC in the laboratory during four months. Leaf litter of either Hevea pauciflora (rubber tree), Flemingia macrophylla (a shrubby legume) or Brachiaria decumbens (a pasture grass) was put on top of each TME. Five specimens of either Pontoscolex corethrurus or Eisenia fetida (earthworms), Porcellionides pruinosus or Circoniscus ornatus (woodlice), and Trigoniulus corallinus (millipedes) were then added to the TMEs. Leaf litter type significantly affected litter consumption, soil microbial biomass and nitrate concentration in the leachate of all TMEs, but had no measurable effect on the shoot biomass of rice seedlings planted in top soil taken from the TMEs. Feeding rates measured with bait lamina were significantly higher in TMEs with the earthworm P. corethrurus and the woodlouse C. ornatus. TMEs are an appropriate tool to assess trophic interactions in tropical soil ecossistems under controlled laboratory conditions.

Temporal variability of soil CO2 emission after conventional and reduced tillage described by an exponential decay in time model

To study Assessing the impact of tillage practices on soil carbon losses dependents it is necessary to describe the temporal variability of soil CO2 emission after tillage. It has been argued that large amounts of CO2 emitted after tillage may serve as an indicator for longer-term changes in soil carbon stocks. Here we present a two-step function model based on soil temperature and soil moisture including an exponential decay in time component that is efficient in fitting intermediate-term emission after disk plow followed by a leveling harrow (conventional), and chisel plow coupled with a roller for clod breaking (reduced) tillage. Emission after reduced tillage was described using a non-linear estimator with determination coefficient (R²) as high as 0.98. Results indicate that when emission after tillage is addressed it is important to consider an exponential decay in time in order to predict the impact of tillage in short-term emissions.

Mathematical model to estimate of the deterioration of wooden poles in contact with soil used in rural areas

In São Paulo State, mainly in rural areas, the utilization of wooden poles is observed for different purposes. In this context, wood in contact with the ground presents faster deterioration, which is generally associated to environmental factors and, especially to the presence of fungi and insects. With the use of mathematical models, the useful life of wooden structures can be predicted by obtaining "climatic indexes" to indicate, comparatively among the areas studied, which have more or less tendency to fungi and insects attacks. In this work, by using climatological data of several cities at São Paulo State, a simplified mathematical model was obtained to measure the aggressiveness of the wood in contact with the soil.

Comparison of maps of spatial variability of soil resistance to penetration constructed with and without covariables using a spatial linear model

A study about the spatial variability of data of soil resistance to penetration (RSP) was conducted at layers 0.0-0.1 m, 0.1-0.2 m and 0.2-0.3 m depth, using the statistical methods in univariate forms, i.e., using traditional geostatistics, forming thematic maps by ordinary kriging for each layer of the study. It was analyzed the RSP in layer 0.2-0.3 m depth through a spatial linear model (SLM), which considered the layers 0.0-0.1 m and 0.1-0.2 m in depth as covariable, obtaining an estimation model and a thematic map by universal kriging. The thematic maps of the RSP at layer 0.2-0.3 m depth, constructed by both methods, were compared using measures of accuracy obtained from the construction of the matrix of errors and confusion matrix. There are similarities between the thematic maps. All maps showed that the RSP is higher in the north region.

Dose-response curve to soil applied herbicides and susceptibility evaluation of different amaranthus species using model identity

Greenhouse studies were conducted in 2008-2009 with the objective of adjusting dose-response curves of the main soil-applied herbicides currently used in cotton for the control of Amaranthus viridis, A. hybridus, A. spinosus, A. lividus, as well as comparing susceptibility among different species, using the identity test models. Thirty six individual experiments were simultaneously carried out in greenhouse, in a sandy clay loam soil (21% clay, 2.36% OM) combining increasing doses of the herbicides alachlor, clomazone, diuron, oxyfluorfen, pendimethalin, prometryn, S-metolachlor, and trifluralin applied to each species. Dose-response curves were adjusted for visual weed control at 28 days after herbicide application and doses required for 80% (C80) and 95% (C95) control were calculated. All herbicides, except clomazone and trifluralin, provided efficient control of most Amaranthus species, but substantial differences in susceptibility to herbicides were found. In general, A. lividus was the least sensitive species, whereas A. spinosus demonstrated the highest sensitivity to herbicides. Alachlor, diuron, oxyfluorfen, pendimethalin, S-metolachlor, and prometryn are efficient alternatives to control Amaranthus spp. in a range of doses that are currently lower than those recommended to cotton.

Evaluation of a linear spectral mixture model and vegetation indices (NDVI and EVI) in a study of schistosomiasis mansoni and Biomphalaria glabrata distribution in the state of Minas Gerais, Brazil

This paper analyses the associations between Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) on the prevalence of schistosomiasis and the presence of Biomphalaria glabrata in the state of Minas Gerais (MG), Brazil. Additionally, vegetation, soil and shade fraction images were created using a Linear Spectral Mixture Model (LSMM) from the blue, red and infrared channels of the Moderate Resolution Imaging Spectroradiometer spaceborne sensor and the relationship between these images and the prevalence of schistosomiasis and the presence of B. glabrata was analysed. First, we found a high correlation between the vegetation fraction image and EVI and second, a high correlation between soil fraction image and NDVI. The results also indicate that there was a positive correlation between prevalence and the vegetation fraction image (July 2002), a negative correlation between prevalence and the soil fraction image (July 2002) and a positive correlation between B. glabrata and the shade fraction image (July 2002). This paper demonstrates that the LSMM variables can be used as a substitute for the standard vegetation indices (EVI and NDVI) to determine and delimit risk areas for B. glabrata and schistosomiasis in MG, which can be used to improve the allocation of resources for disease control.

Stratigraphical discontinuities, tropical landscape evolution and soil distribution relationships in a case study in SE-Brazil

On a regional summit surface in the county of Piracicaba (SP) within the Peripheric Depression of São Paulo, formed of discontinued flattened tops, there is an abrupt transition between a Typic Hapludox and a Kandiudalfic Eutrudox, together with two stoneline layers. Using stratigraphical, mineralogical, and cartographic studies, this transition and the soil distribution of this surface were studied, correlating them with the different parent materials and the morphoclimatic model of landscape evolution in Southeastern Brazil. The Typic Hapludox was formed on a sandy Cenozoic deposit (Q) that overlies a pellitic deposit of the Iratí formation (Pi), representing a regional erosive discordance. Westwards to the Piracicaba River, this sequence is interrupted by a diabase sill overlain by a red clayey material which gave origin to the Kandiudalfic Eutrudox. Two post-Permian depositional events were identified by the two stonelines and stratigraphical discontinuities. The first event generated the deposition of a sandy sediment in the form of levelled alveoluses on regional barriers, most of these formed by dikes and diabase sills, probably during a drier phase. The second depositional event, leading to the deposition of the red clay was probably the dissection of the previously formed pediplane during a humid climate, followed by another pedimentation process during a later, drier period.

Reconsolidation of the soil surface after tillage discontinuity, with and without cultivation, related to erosion and its prediction with RUSLE

Site-specific regression coefficient values are essential for erosion prediction with empirical models. With the objective to investigate the surface-soilconsolidation factor, Cf, linked to the RUSLE's prior-land-use subfactor, PLU, an erosion experiment using simulated rainfall on a 0.075 m m-1 slope, sandy loam Paleudult soil, was conducted at the Agriculture Experimental Station of the Federal University of Rio Grande do Sul (EEA/UFRGS), in Eldorado do Sul, State of Rio Grande do Sul, Brazil. Firstly, a row-cropped area was excluded from cultivation (March 1995), the existing crop residue removed from the field, and the soil kept clean-tilled the rest of the year (to get a degraded soil condition for the intended purpose of this research). The soil was then conventional-tilled for the last time (except for a standard plot which was kept continuously cleantilled for comparison purposes), in January 1996, and the following treatments were established and evaluated for soil reconsolidation and soil erosion until May 1998, on duplicated 3.5 x 11.0 m erosion plots: (a) fresh-tilled soil, continuously in clean-tilled fallow (unit plot); (b) reconsolidating soil without cultivation; and (c) reconsolidating soil with cultivation (a crop sequence of three corn- and two black oats cycles, continuously in no-till, removing the crop residues after each harvest for rainfall application and redistributing them on the site after that). Simulated rainfall was applied with a Swanson's type, rotating-boom rainfall simulator, at 63.5 mm h-1 intensity and 90 min duration, six times during the two-and-half years of experimental period (at the beginning of the study and after each crop harvest, with the soil in the unit plot being retilled before each rainfall test). The soil-surface-consolidation factor, Cf, was calculated by dividing soil loss values from the reconsolidating soil treatments by the average value from the fresh-tilled soil treatment (unit plot). Non-linear regression was used to fit the Cf = e b.t model through the calculated Cf-data, where t is time in days since last tillage. Values for b were -0.0020 for the reconsolidating soil without cultivation and -0.0031 for the one with cultivation, yielding Cf-values equal to 0.16 and 0.06, respectively, after two-and-half years of tillage discontinuation, compared to 1.0 for fresh-tilled soil. These estimated Cf-values correspond, respectively, to soil loss reductions of 84 and 94 %, in relation to soil loss from the fresh-tilled soil, showing that the soil surface reconsolidated intenser with cultivation than without it. Two distinct treatmentinherent soil surface conditions probably influenced the rapid decay-rate of Cf values in this study, but, as a matter of a fact, they were part of the real environmental field conditions. Cf-factor curves presented in this paper are therefore useful for predicting erosion with RUSLE, but their application is restricted to situations where both soil type and particular soil surface condition are similar to the ones investigate in this study.

Sampling and TDR probe insertion in the determination of the volumetric soil water content

Volumetric soil water content (theta) can be evaluated in the field by direct or indirect methods. Among the direct, the gravimetric method is regarded as highly reliable and thus often preferred. Its main disadvantages are that sampling and laboratory procedures are labor intensive, and that the method is destructive, which makes resampling of a same point impossible. Recently, the time domain reflectometry (TDR) technique has become a widely used indirect, non-destructive method to evaluate theta. In this study, evaluations of the apparent dielectric number of soils (epsilon) and samplings for the gravimetrical determination of the volumetric soil water content (thetaGrav) were carried out at four sites of a Xanthic Ferralsol in Manaus - Brazil. With the obtained epsilon values, theta was estimated using empirical equations (thetaTDR), and compared with thetaGrav derived from disturbed and undisturbed samples. The main objective of this study was the comparison of thetaTDR estimates of horizontally as well as vertically inserted probes with the thetaGrav values determined by disturbed and undisturbed samples. Results showed that thetaTDR estimates of vertically inserted probes and the average of horizontally measured layers were only slightly and insignificantly different. However, significant differences were found between the thetaTDR estimates of different equations and between disturbed and undisturbed samples in the thetaGrav determinations. The use of the theoretical Knight et al. model, which permits an evaluation of the soil volume assessed by TDR probes, is also discussed. It was concluded that the TDR technique, when properly calibrated, permits in situ, nondestructive measurements of q in Xanthic Ferralsols of similar accuracy as the gravimetric method.

State-space approach to evaluate the relation between soil physical and chemical properties

The state-space approach is used to evaluate the relation between soil physical and chemical properties in an area cultivated with sugarcane. The experiment was carried out on a Rhodic Kandiudalf in Piracicaba, State of São Paulo, Brazil. Sugarcane was planted on an area of 0.21 ha i.e., in 15 rows 100 m long, spaced 1.4 m. Soil water content, soil organic matter, clay content and aggregate stability were sampled along a transect of 84 points, meter by meter. The state-space approach is used to evaluate how the soil water content is affected by itself and by soil organic matter, clay content, and aggregate stability of neighboring locations, in different combinations, aiming to contribute to a better understanding of the relation among these variables in the soil. Results show that soil water contents were successfully estimated by this approach. Best performances were found when the estimate of soil water content at locations i was related to soil water content, clay content and aggregate stability at locations i-1. Results also indicate that this state-space model using all series describes the soil water content better than any equivalent multiple regression equation.

Applicability of a distributed watershed pollution model in a data-poor environment in Santa Catarina State, Brazil

Intensification of agricultural production without a sound management and regulations can lead to severe environmental problems, as in Western Santa Catarina State, Brazil, where intensive swine production has caused large accumulations of manure and consequently water pollution. Natural resource scientists are asked by decision-makers for advice on management and regulatory decisions. Distributed environmental models are useful tools, since they can be used to explore consequences of various management practices. However, in many areas of the world, quantitative data for model calibration and validation are lacking. The data-intensive distributed environmental model AgNPS was applied in a data-poor environment, the upper catchment (2,520 ha) of the Ariranhazinho River, near the city of Seara, in Santa Catarina State. Steps included data preparation, cell size selection, sensitivity analysis, model calibration and application to different management scenarios. The model was calibrated based on a best guess for model parameters and on a pragmatic sensitivity analysis. The parameters were adjusted to match model outputs (runoff volume, peak runoff rate and sediment concentration) closely with the sparse observed data. A modelling grid cell resolution of 150 m adduced appropriate and computer-fit results. The rainfall runoff response of the AgNPS model was calibrated using three separate rainfall ranges (< 25, 25-60, > 60 mm). Predicted sediment concentrations were consistently six to ten times higher than observed, probably due to sediment trapping along vegetated channel banks. Predicted N and P concentrations in stream water ranged from just below to well above regulatory norms. Expert knowledge of the area, in addition to experience reported in the literature, was able to compensate in part for limited calibration data. Several scenarios (actual, recommended and excessive manure applications, and point source pollution from swine operations) could be compared by the model, using a relative ranking rather than quantitative predictions.

Assessment of the soil compaction of two ultisols caused by logging operations

The impact of wood loads on bulk density and preconsolidation pressure and of harvester and forwarder traffic on rut depth, bulk density and preconsolidation pressure of two Ultisols were examined in this study. Our objective was to quantify the threshold beyond which significant soil compaction and rutting would occur. This study was carried out in the county of Eunápolis, state of Bahia, Brazil, (16 º 23 ' 17 '' S and 39 º 10 ' 06 '' W; altitude 80 m asl) in two Ultisols (PAd2 and PAd3) with different texture classes, in experimental areas with eucalypt plantation. The study involved measurements at the wood load site and machine driving at specific locations in the forest during logging operations. The treatments consisted of one harvester pass and, 8, 16 and 40 passes of a fully loaded forwarder. Thresholds were established based on the rut depth and percentage of preconsolidation pressure values in the region of additional soil compaction defined in the bearing capacity model. The percentage of soil samples with values of preconsolidation pressure in the region of additional soil compaction indicated a greater susceptibility of PAd3 than of PAd2 to soil compaction. The threshold levels established here based on preconsolidation pressure and rut depth indicated that no more than eight forwarder passes should be allowed in loading operations in order to minimize soil compaction.

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.

The combined use of reflectance, emissivity and elevation Aster/Terra data for tropical soil studies

Reflectance, emissivity and elevation data of the sensor ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer)/Terra were used to characterize soil composition variations according to the toposequence position. Normalized data of SWIR (shortwave infrared) reflectance and TIR (thermal infrared) emissivity, coupled to a soil-fraction image from a spectral mixture model, were evaluated to separate bare soils from nonphotosynthetic vegetation. Regression relationships of some soil properties with reflectance and emissivity data were then applied on the exposed soil pixels. The resulting estimated values were plotted on the ASTER-derived digital elevation model. Results showed that the SWIR bands 5 and 6 and the TIR bands 10 and 14 measured the clay mineral absorption band and the quartz emissivity feature, respectively. These bands improved also the discrimination between nonphotosynthetic vegetation and soils. Despite the differences in pixel size and field sampling size, some soil properties were correlated with reflectance (R² of 0.65 for Al2O3 in band 6; 0.61 for Fe2O3 in band 3) and emissivity (R² of 0.65 for total sand fraction in the 10/14 band ratio). The combined use of reflectance, emissivity and elevation data revealed variations in soil composition with topography in specific parts of the landscape. From higher to lower slope positions, a general decrease in Al2O3 and increase in total sand fraction was observed, due to the prevalence of Rhodic Acrustox at the top and its gradual transition to Typic Acrustox at the bottom.