Hydraulic conductivity and water retention in leptosols-regosols and saprolite derived from sandstone, Brazil

Leptosols and Regosols are soils with a series of restrictions for use, mainly related to the effective depth, which have been poorly studied in Brazil. These soils, when derived from sedimentary rocks should be treated with particular care to avoid environmental damage such as aquifer contamination. The purpose of this study was to verify the behavior of hydraulic conductivity and water retention capacity in profiles of Leptosols and Regosols derived from sandstone of the Caturrita formation in Rio Grande do Sul state. The morphology, particle size distribution, porosity, soil density (Ds), saturated hydraulic conductivity (Ks), basic water infiltration in the field (BI) and water retention were determined in soil and saprolite samples of six soil profiles. High Ds, low macroporosity and high microporosity were observed in the profiles, resulting in a low Ks and BI, even under conditions of sandy texture and a highly fractured saprolite layer. The variation coefficients of data of Ks and BI were high among the studied profiles and between replications of a same profile. Water retention of the studied soils was higher in Cr layers than in the A horizons and the volume of plant-available water greater and variable among A horizons and Cr layers.

Hydraulic conductivity and soil-sewage sludge interactions

One of the main problems faced by humanity is pollution caused by residues resulting from the production and use of goods, e.g, sewage sludge. Among the various alternatives for its disposal, the agricultural use seems promising. The purpose of this study was to evaluate the hydraulic conductivity and interaction of soil with sandy-silty texture, classified as Spodosols, from the Experimental Station Itapirema - IPA, in Goiana, state of Pernambuco, in mixtures with sewage sludge from the Mangueira Sewage Treatment Station, in the city of Recife, Pernambuco at rates of 25, 50 and 75 Mg ha-1. Tests were conducted to let water percolate the natural saturated soil and soil-sludge mixtures to characterize their physical, chemical, and microstructural properties as well as hydraulic conductivity. Statistical data analysis showed that the presence of sewage sludge in soils leads to an increase of the < 0.005 mm fraction, reduction in real specific weight and variation in optimum moisture content from 11.60 to 12.90 % and apparent specific dry weight from 17.10 and 17.50 kN m-3. In the sludge-soil mixture, the quartz grains were covered by sludge and filling of the empty soil macropores between grains. There were changes in the chemical characteristics of soil and effluent due to sewage sludge addition and a small decrease in hydraulic conductivity. The results indicate the possibility that soil acidity influenced the concentrations of the elements found in the leachate, showing higher levels at higher sludge doses. It can be concluded that the leaching degree of potentially toxic elements from the sewage sludge treatments does not harm the environment.

Assessment of pedotransfer functions for estimating soil water retention curves for the amazon region

Knowledge of the soil water retention curve (SWRC) is essential for understanding and modeling hydraulic processes in the soil. However, direct determination of the SWRC is time consuming and costly. In addition, it requires a large number of samples, due to the high spatial and temporal variability of soil hydraulic properties. An alternative is the use of models, called pedotransfer functions (PTFs), which estimate the SWRC from easy-to-measure properties. The aim of this paper was to test the accuracy of 16 point or parametric PTFs reported in the literature on different soils from the south and southeast of the State of Pará, Brazil. The PTFs tested were proposed by Pidgeon (1972), Lal (1979), Aina & Periaswamy (1985), Arruda et al. (1987), Dijkerman (1988), Vereecken et al. (1989), Batjes (1996), van den Berg et al. (1997), Tomasella et al. (2000), Hodnett & Tomasella (2002), Oliveira et al. (2002), and Barros (2010). We used a database that includes soil texture (sand, silt, and clay), bulk density, soil organic carbon, soil pH, cation exchange capacity, and the SWRC. Most of the PTFs tested did not show good performance in estimating the SWRC. The parametric PTFs, however, performed better than the point PTFs in assessing the SWRC in the tested region. Among the parametric PTFs, those proposed by Tomasella et al. (2000) achieved the best accuracy in estimating the empirical parameters of the van Genuchten (1980) model, especially when tested in the top soil layer.

Stabilization of soil hydraulic properties under a long term no-till system

The area under the no-tillage system (NT) has been increasing over the last few years. Some authors indicate that stabilization of soil physical properties is reached after some years under NT while other authors debate this. The objective of this study was to determine the effect of the last crop in the rotation sequence (1st year: maize, 2nd year: soybean, 3rd year: wheat/soybean) on soil pore configuration and hydraulic properties in two different soils (site 1: loam, site 2: sandy loam) from the Argentinean Pampas region under long-term NT treatments in order to determine if stabilization of soil physical properties is reached apart from a specific time in the crop sequence. In addition, we compared two procedures for evaluating water-conducting macroporosities, and evaluated the efficiency of the pedotransfer function ROSETTA in estimating the parameters of the van Genuchten-Mualem (VGM) model in these soils. Soil pore configuration and hydraulic properties were not stable and changed according to the crop sequence and the last crop grown in both sites. For both sites, saturated hydraulic conductivity, K0, water-conducting macroporosity, εma, and flow-weighted mean pore radius, R0ma, increased from the 1st to the 2nd year of the crop sequence, and this was attributed to the creation of water-conducting macropores by the maize roots. The VGM model adequately described the water retention curve (WRC) for these soils, but not the hydraulic conductivity (K) vs tension (h) curve. The ROSETTA function failed in the estimation of these parameters. In summary, mean values of K0 ranged from 0.74 to 3.88 cm h-1. In studies on NT effects on soil physical properties, the crop effect must be considered.

Physical-hydraulic properties of a sandy loam typic paleudalf soil under organic cultivation of 'montenegrina' mandarin (Citrus deliciosa Tenore)¹

Citrus plants are the most important fruit species in the world, with emphasis to oranges, mandarins and lemons. In Rio Grande do Sul, Brazil, most fruit production is found on small properties under organic cultivation. Soil compaction is one of the factors limiting production and due to the fixed row placement of this crop, compaction can arise in various manners in the interrows of the orchard. The aim of this study was to evaluate soil physical properties and water infiltration capacity in response to interrow management in an orchard of mandarin (Citrus deliciosa Tenore 'Montenegrina') under organic cultivation. Interrow management was performed through harrowing, logs in em "V", mowing, and cutting/knocking down plants with a knife roller. Soil physical properties were evaluated in the wheel tracks of the tractor (WT), between the wheel tracks (BWT), and in the area under the line projection of the canopy (CLP), with undisturbed soil samples collected in the 0.00-0.15, 0.15-0.30, 0.30-0.45, and 0.45-0.60 m layers, with four replicates. The soil water infiltration test was performed using the concentric cylinder method, with a maximum time of 90 min for each test. In general, soil analysis showed a variation in the physical-hydraulic properties of the Argissolo Vermelho-Amarelo distrófico arênico (sandy loam Typic Paleudalf) in the three sampling sites in all layers, regardless of the management procedure in the interrows. Machinery traffic leads to heterogeneity in the soil physical-hydraulic properties in the interrows of the orchard. Soil porosity and bulk density are affected especially in the wheel tracks of the tractor (WT), which causes a reduction in the constant rate of infiltration and in the accumulated infiltration of water in this sampling site. The use of the disk harrow and mower leads to greater harmful effects on the soil, which can interfere with mandarin production.

MODELING OF SOIL LOAD-BEARING CAPACITY AS A FUNCTION OF SOIL MECHANICAL RESISTANCE TO PENETRATION

Estimation of soil load-bearing capacity from mathematical models that relate preconsolidation pressure (σp) to mechanical resistance to penetration (PR) and gravimetric soil water content (U) is important for defining strategies to prevent compaction of agricultural soils. Our objective was therefore to model the σp and compression index (CI) according to the PR (with an impact penetrometer in the field and a static penetrometer inserted at a constant rate in the laboratory) and U in a Rhodic Eutrudox. The experiment consisted of six treatments: no-tillage system (NT); NT with chiseling; and NT with additional compaction by combine traffic (passing 4, 8, 10, and 20 times). Soil bulk density, total porosity, PR (in field and laboratory measurements), U, σp, and CI values were determined in the 5.5-10.5 cm and 13.5-18.5 cm layers. Preconsolidation pressure (σp) and CI were modeled according to PR in different U. The σp increased and the CI decreased linearly with increases in the PR values. The correlations between σp and PR and PR and CI are influenced by U. From these correlations, the soil load-bearing capacity and compaction susceptibility can be estimated by PR readings evaluated in different U.

Surface Complexation Modeling in Variable Charge Soils: Charge Characterization by Potentiometric Titration

ABSTRACT Intrinsic equilibrium constants of 17 representative Brazilian Oxisols were estimated from potentiometric titration measuring the adsorption of H+ and OH− on amphoteric surfaces in suspensions of varying ionic strength. Equilibrium constants were fitted to two surface complexation models: diffuse layer and constant capacitance. The former was fitted by calculating total site concentration from curve fitting estimates and pH-extrapolation of the intrinsic equilibrium constants to the PZNPC (hand calculation), considering one and two reactive sites, and by the FITEQL software. The latter was fitted only by FITEQL, with one reactive site. Soil chemical and physical properties were correlated to the intrinsic equilibrium constants. Both surface complexation models satisfactorily fit our experimental data, but for results at low ionic strength, optimization did not converge in FITEQL. Data were incorporated in Visual MINTEQ and they provide a modeling system that can predict protonation-dissociation reactions in the soil surface under changing environmental conditions.

Surface Complexation Modeling in Variable Charge Soils: Prediction of Cadmium Adsorption

ABSTRACT Intrinsic equilibrium constants for 22 representative Brazilian Oxisols were estimated from a cadmium adsorption experiment. Equilibrium constants were fitted to two surface complexation models: diffuse layer and constant capacitance. Intrinsic equilibrium constants were optimized by FITEQL and by hand calculation using Visual MINTEQ in sweep mode, and Excel spreadsheets. Data from both models were incorporated into Visual MINTEQ. Constants estimated by FITEQL and incorporated in Visual MINTEQ software failed to predict observed data accurately. However, FITEQL raw output data rendered good results when predicted values were directly compared with observed values, instead of incorporating the estimated constants into Visual MINTEQ. Intrinsic equilibrium constants optimized by hand calculation and incorporated in Visual MINTEQ reliably predicted Cd adsorption reactions on soil surfaces under changing environmental conditions.

Modeling growth and yield of loblolly pinestands under intensive management

Abstract:The objective of this work was to develop and validate a prognosis system for volume yield and basal area of intensively managed loblolly pine (Pinus taeda) stands, using stand and diameter class models compatible in basal area estimates. The data used in the study were obtained from plantations located in northern Uruguay. For model validation without data loss, a three-phase validation scheme was applied: first, the equations were fitted without the validation database; then, model validation was carried out; and, finally, the database was regrouped to recalibrate the parameter values. After the validation and final parameterization of the models, a simulation of the first commercial thinning was carried out. The developed prognosis system was precise and accurate in estimating basal area production per hectare or per diameter classes. There was compatibility in basal area estimates between diameter class and whole stand models, with a mean difference of -0.01 m2ha-1. The validation scheme applied is logic and consistent, since information on the accuracy and precision of the models is obtained without the loss of any information in the estimation of the models' parameters.

Phase diagram and complex patterns in the modeling of the bromate-oxalic acid-Ce-acetone oscillating reaction

Simulations have been carried out on the bromate - oxalic acid - Ce(IV) - acetone oscillating reaction, under flow conditions, using Field and Boyd's model (J. Phys. Chem. 1985, 89, 3707). Many different complex dynamic behaviors were found, including simple periodic oscillations, complex periodic oscillations, quasiperiodicity and chaos. Some of these complex oscillations can be understood as belonging to a Farey sequence. The many different behaviors were systematized in a phase diagram which shows that some regions of complex patterns were nested with one inside the other. The existence of almost all known dynamic behavior for this system allows the suggestion that it can be used as a model for some very complex phenomena that occur in biological systems.

Supercritical CO2 recovery of caffeine from green coffee oil: new experimental solubility data and modeling

The caffeine solubility in supercritical CO2 was studied by assessing the effects of pressure and temperature on the extraction of green coffee oil (GCO). The Peng-Robinson¹ equation of state was used to correlate the solubility of caffeine with a thermodynamic model and two mixing rules were evaluated: the classical mixing rule of van der Waals with two adjustable parameters (PR-VDW) and a density dependent one, proposed by Mohamed and Holder² with two (PR-MH, two parameters adjusted to the attractive term) and three (PR-MH3 two parameters adjusted to the attractive and one to the repulsive term) adjustable parameters. The best results were obtained with the mixing rule of Mohamed and Holder² with three parameters.

Temperature measurements by oh lif and chemiluminescence kinetic modeling for ethanol flames

OH LIF-thermometry was applied to premixed ethanol flames at atmospheric pressure in a burner for three flame conditions. Flame temperatures were simulated from energy equation with PREMIX code of CHEMKIN software package for comparison. A kinetic modeling based on a model validated through chemiluminescence measurements and on a set of reactions for nitrogen chemistry was evaluated. Marinov's mechanism was also tested. Sensitivity analysis was performed for fuel-rich flame condition with Φ = 1.34. Simulated temperatures from both reaction mechanisms evaluated were higher than experimental values. However, the proposed kinetic modeling resulted in temperature profiles qualitatively very close to the experimental.

QSAR modeling: um novo pacote computacional open source para gerar e validar modelos QSAR

QSAR modeling is a novel computer program developed to generate and validate QSAR or QSPR (quantitative structure- activity or property relationships) models. With QSAR modeling, users can build partial least squares (PLS) regression models, perform variable selection with the ordered predictors selection (OPS) algorithm, and validate models by using y-randomization and leave-N-out cross validation. An additional new feature is outlier detection carried out by simultaneous comparison of sample leverage with the respective Studentized residuals. The program was developed using Java version 6, and runs on any operating system that supports Java Runtime Environment version 6. The use of the program is illustrated. This program is available for download at lqta.iqm.unicamp.br.

Kinetic modeling of the biosorption of Cd2+ ions from aqueous solutions onto Eichhornia crassipes roots using potentiometry: low-cost alternative to conventional methods

This work presents the use of potentiometric measurements for kinetic studies of biosorption of Cd2+ ions from aqueous solutions on Eichhornia crassipes roots. The open circuit potential of the Cd/Cd2+ electrode of the first kind was measured during the bioadsorption process. The amount of Cd2+ ions accumulated was determined in real time. The data were fit to different models, with the pseudo-second-order model proving to be the best in describing the data. The advantages and limitations of the methodology proposed relative to the traditional method are discussed.

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.