Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Structural changes in latosols of the cerrado region: II - soil compressive behavior and modeling of additional compaction

Currently in Brazil, as in other parts of the world, the concern is great with the increase of degraded agricultural soil, which is mostly related to the occurrence of soil compaction. Although soil texture is recognized as a very important component in the soil compressive behaviors, there are few studies that quantify its influence on the structural changes of Latosols in the Brazilian Cerrado region. This study aimed to evaluate structural changes and the compressive behavior of Latosols in Rio Verde, Goiás, through the modeling of additional soil compaction. The study was carried out using five Latosols with very different textures, under different soil compaction levels. Water retention and soil compression curves, and bearing capacity models were determined from undisturbed samples collected on the B horizons. Results indicated that clayey and very clayey Latosols were more susceptible to compression than medium-textured soils. Soil compression curves at density values associate with edaphic functions were used to determine the beneficial pressure (σ b) , i.e., pressure with optimal water retention, and critical pressure (σcrMAC), i.e., pressure with macroporosity below critical levels. These pressure values were higher than the preconsolidation pressure (σp), and therefore characterized as additional compaction. Based on the compressive behavior of these Latosols, it can be concluded that the combined preconsolidation pressure, beneficial pressure and critical pressure allow a better understanding of compression processes of Latosols.

Modeling changes in organic carbon stocks for distinct soils in southeastern brazil after four eucalyptus rotations using the century model

Soil organic matter (SOM) plays an important role in carbon (C) cycle and soil quality. Considering the complexity of factors that control SOM cycling and the long time it usually takes to observe changes in SOM stocks, modeling constitutes a very important tool to understand SOM cycling in forest soils. The following hypotheses were tested: (i) soil organic carbon (SOC) stocks would be higher after several rotations of eucalyptus than in low-productivity pastures; (ii) SOC values simulated by the Century model would describe the data better than the mean of observations. So, the aims of the current study were: (i) to evaluate the SOM dynamics using the Century model to simulate the changes of C stocks for two eucalyptus chronosequences in the Rio Doce Valley, Minas Gerais State, Brazil; and (ii) to compare the C stocks simulated by Century with the C stocks measured in soils of different Orders and regions of the Rio Doce Valley growing eucalyptus. In Belo Oriente (BO), short-rotation eucalyptus plantations had been cultivated for 4.0; 13.0, 22.0, 32.0 and 34.0 years, at a lower elevation and in a warmer climate, while in Virginópolis (VG), these time periods were 8.0, 19.0 and 33.0 years, at a higher elevation and in a milder climate. Soil samples were collected from the 0-20 cm layer to estimate C stocks. Results indicate that the C stocks simulated by the Century model decreased after 37 years of poorly managed pastures in areas previously covered by native forest in the regions of BO and VG. The substitution of poorly managed pastures by eucalyptus in the early 1970´s led to an average increase of C of 0.28 and 0.42 t ha-1 year-1 in BO and VG, respectively. The measured C stocks under eucalyptus in distinct soil Orders and independent regions with variable edapho-climate conditions were not far from the values estimated by the Century model (root mean square error - RMSE = 20.9; model efficiency - EF = 0.29) despite the opposite result obtained with the statistical procedure to test the identity of analytical methods. Only for lower soil C stocks, the model over-estimated the C stock in the 0-20 cm layer. Thus, the Century model is highly promising to detect changes in C stocks in distinct soil orders under eucalyptus, as well as to indicate the impact of harvest residue management on SOM in future rotations.

Determination of hydraulic properties of a tropical soil of Hawaii using column experiments and inverse modeling

A method for determining soil hydraulic properties of a weathered tropical soil (Oxisol) using a medium-sized column with undisturbed soil is presented. The method was used to determine fitting parameters of the water retention curve and hydraulic conductivity functions of a soil column in support of a pesticide leaching study. The soil column was extracted from a continuously-used research plot in Central Oahu (Hawaii, USA) and its internal structure was examined by computed tomography. The experiment was based on tension infiltration into the soil column with free outflow at the lower end. Water flow through the soil core was mathematically modeled using a computer code that numerically solves the one-dimensional Richards equation. Measured soil hydraulic parameters were used for direct simulation, and the retention and soil hydraulic parameters were estimated by inverse modeling. The inverse modeling produced very good agreement between model outputs and measured flux and pressure head data for the relatively homogeneous column. The moisture content at a given pressure from the retention curve measured directly in small soil samples was lower than that obtained through parameter optimization based on experiments using a medium-sized undisturbed soil column.

Correction of resistance to penetration by pedofunctions and a reference soil water content

The soil penetration resistance is an important indicator of soil compaction and is strongly influenced by soil water content. The objective of this study was to develop mathematical models to normalize soil penetration resistance (SPR), using a reference value of gravimetric soil water content (U). For this purpose, SPR was determined with an impact penetrometer, in an experiment on a Dystroferric Red Latossol (Rhodic Eutrudox), at six levels of soil compaction, induced by mechanical chiseling and additional compaction by the traffic of a harvester (four, eight, 10, and 20 passes); in addition to a control treatment under no-tillage, without chiseling or additional compaction. To broaden the range of U values, SPR was evaluated in different periods. Undisturbed soil cores were sampled to quantify the soil bulk density (BD). Pedotransfer functions were generated correlating the values of U and BD to the SPR values. By these functions, the SPR was adequately corrected for all U and BD data ranges. The method requires only SPR and U as input variables in the models. However, different pedofunctions are needed according to the soil layer evaluated. After adjusting the pedotransfer functions, the differences in the soil compaction levels among the treatments, previously masked by variations of U, became detectable.

Fuzzy logic applied to the modeling of water dynamics in an Oxisol in northeastern Brazil

Modeling of water movement in non-saturated soil usually requires a large number of parameters and variables, such as initial soil water content, saturated water content and saturated hydraulic conductivity, which can be assessed relatively easily. Dimensional flow of water in the soil is usually modeled by a nonlinear partial differential equation, known as the Richards equation. Since this equation cannot be solved analytically in certain cases, one way to approach its solution is by numerical algorithms. The success of numerical models in describing the dynamics of water in the soil is closely related to the accuracy with which the water-physical parameters are determined. That has been a big challenge in the use of numerical models because these parameters are generally difficult to determine since they present great spatial variability in the soil. Therefore, it is necessary to develop and use methods that properly incorporate the uncertainties inherent to water displacement in soils. In this paper, a model based on fuzzy logic is used as an alternative to describe water flow in the vadose zone. This fuzzy model was developed to simulate the displacement of water in a non-vegetated crop soil during the period called the emergency phase. The principle of this model consists of a Mamdani fuzzy rule-based system in which the rules are based on the moisture content of adjacent soil layers. The performances of the results modeled by the fuzzy system were evaluated by the evolution of moisture profiles over time as compared to those obtained in the field. The results obtained through use of the fuzzy model provided satisfactory reproduction of soil moisture profiles.

Quantification of permanent and variable charges in reference soils of the State of Pernambuco, Brazil

The electrical charges in soil particles are divided into structural or permanent charges and variable charges. Permanent charges develop on the soil particle surface by isomorphic substitution. Variable charges arise from dissociation and association of protons (H+), protonation or deprotonation, and specific adsorption of cations and anions. The aim of this study was to quantify the permanent charges and variable charges of Reference Soils of the State of Pernambuco, Brazil. To do so, 24 subsurface profiles from different regions (nine in the Zona da Mata, eight in the Agreste, and seven in the Sertão) were sampled, representing approximately 80 % of the total area of the state. Measurements were performed using cesium chloride solution. Determination was made of the permanent charges and the charges in regard to the hydroxyl functional groups through selective ion exchange of Cs+ by Li+ and Cs+ by NH4+, respectively. All the soils analyzed exhibited variable cation exchange capacity, with proportions from 0.16 to 0.60 and an average of 0.40 when related to total cation exchange capacity.

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.

Background and Reference Values of Metals in Soils from Paraíba State, Brazil

ABSTRACT Soil contamination by heavy metals threatens ecosystems and human health. Environmental monitoring bodies need reference values for these contaminants to assess the impacts of anthropogenic activities on soil contamination. Quality reference values (QRVs) reflect the natural concentrations of heavy metals in soils without anthropic interference and must be regionally established. The aim of this study was to determine the natural concentrations and quality reference values for the metals Ag, Ba, Cd, Co, Cu, Cr, Mo, Ni, Pb, Sb and Zn in soils of Paraíba state, Brazil. Soil samples were collected from 94 locations across the state in areas of native vegetation or with minimal anthropic interference. The quality reference values (QRVs) were (mg kg-1): Ag (<0.53), Ba (117.41), Cd (0.08), Co (13.14), Cu (20.82), Cr (48.35), Mo (0.43), Ni (14.44), Sb (0.61), Pb (14.62) and Zn (33.65). Principal component analysis grouped the metals Cd, Cr, Cu, Ni, Pb and Sb (PC1); Ag (PC2); and Ba, Co, Fe, Mn and Zn (PC3). These values were made official by Paraíba state through Normativa Resolution 3602/2014.

Identification of reference genes for expression analysis by real-time quantitative PCR in drought-stressed soybean

The objective of this work was to validate, by quantitative PCR in real time (RT-qPCR), genes to be used as reference in studies of gene expression in soybean in drought-stressed trials. Four genes commonly used in soybean were evaluated: Gmβ-actin, GmGAPDH, GmLectin and GmRNAr18S. Total RNA was extracted from six samples: three from roots in a hydroponic system with different drought intensities (0, 25, 50, 75 and 100 minutes of water stress), and three from leaves of plants grown in sand with different soil moistures (15, 5 and 2.5% gravimetric humidity). The raw cycle threshold (Ct) data were analyzed, and the efficiency of each primer was calculated for an overall analysis of the Ct range among the different samples. The GeNorm application was used to evaluate the best reference gene, according to its stability. The GmGAPDH was the least stable gene, with the highest mean values of expression stability (M), and the most stable genes, with the lowest M values, were the Gmβ-actin and GmRNAr18S, when both root and leaves samples were tested. These genes can be used in RT-qPCR as reference gene for expression analysis.

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.

Suggestion of a national diagnostic reference level for 18F-FDG/PET scans in adult cancer patients in Brazil

Objective To suggest a national value for the diagnostic reference level (DRL) in terms of activity in MBq.kg–1, for nuclear medicine procedures with fluorodeoxyglucose (18F-FDG) in whole body positron emission tomography (PET) scans of adult patients. Materials and Methods A survey on values of 18F-FDG activity administered in Brazilian clinics was undertaken by means of a questionnaire including questions about number and manufacturer of the installed equipment, model and detector type. The suggested DRL value was based on the calculation of the third quartile of the activity values distribution reported by the clinics. Results Among the surveyed Brazilian clinics, 58% responded completely or partially the questionnaire; and the results demonstrated variation of up to 100% in the reported radiopharmaceutical activity. The suggested DRL for 18F-FDG/PET activity was 5.54 MBq.kg–1 (0.149 mCi.kg–1). Conclusion The present study has demonstrated the lack of standardization in administered radiopharmaceutical activities for PET procedures in Brazil, corroborating the necessity of an official DRL value to be adopted in the country. The suggested DLR value demonstrates that there is room for optimization of the procedures and 18F-FDG/PET activities administered in Brazilian clinics to reduce the doses delivered to patients. It is important to highlight that this value should be continually revised and optimized at least every five years.

Complete internal audit of a mammography service in a reference institution for breast imaging

Objective Undertaking of a complete audit of the service of mammography, as recommended by BI-RADS®, in a private reference institution for breast cancer diagnosis in the city of São Paulo, SP, Brazil, and comparison of results with those recommended by the literature. Materials and Methods Retrospective, analytical and cross-sectional study including 8,000 patients submitted to mammography in the period between April 2010 and March 2011, whose results were subjected to an internal audit. The patients were followed-up until December 2012. Results The radiological classification of 7,249 screening mammograms, according to BI-RADS, was the following: category 0 (1.43%), 1 (7.82%), 2 (80.76%), 3 (8.35%), 4 (1.46%), 5 (0.15%) and 6 (0.03%). The breast cancer detection ratio was 4.8 cases per 1,000 mammograms. Ductal carcinoma in situ was found in 22.8% of cases. Positive predictive values for categories 3, 4 and 5 were 1.3%, 41.3% and 100%, respectively. In the present study, the sensitivity of the method was 97.1% and specificity, 97.4%. Conclusion The complete internal audit of a service of mammography is essential to evaluate the quality of such service, which reflects on an early breast cancer detection and reduction of mortality rates.