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).

An improved method for measuring soil microbial activity by gas phase flow injection analysis

The rate of carbon dioxide production is commonly used as a measure of microbial activity in the soil. The traditional method of CO2 determination involves trapping CO2 in an alkali solution and then determining CO2 concentration indirectly by titration of the remaining alkali in the solution. This method is still commonly employed in laboratories throughout the world due to its relative simplicity and the fact that it does not require expensive, specific equipment. However, there are several drawbacks: the method is time-consuming, requires large amounts of chemicals and the consistency of results depends on the operator's skills. With this in mind, an improved method was developed to analyze CO2 captured in alkali traps, which is cheap and relatively simple, with a substantially shorter sample handling time and reproducibility equivalent to the traditional titration method. A comparison of the concentration values determined by gas phase flow injection analysis (GPFIA) and titration showed no significant difference (p > 0.05), but GPFIA has the advantage that only a tenth of the sample volume of the titration method is required. The GPFIA system does not require the purchase of new, costly equipment but the device was constructed from items commonly found in laboratories, with suggestions for alternative configurations for other detection units. Furthermore, GPFIA for CO2 analysis can be equally applied to samples obtained from either the headspace of microcosms or from a sampling chamber that allows CO2 to be released from alkali trapping solutions. The optimised GPFIA method was applied to analyse CO2 released from degrading hydrocarbons from a site contaminated by diesel spillage.

SPATIAL UNCERTAINTY OF NUTRIENT LOSS BY EROSION IN SUGARCANE HARVESTING SCENARIOS

The assessment of spatial uncertainty in the prediction of nutrient losses by erosion associated with landscape models is an important tool for soil conservation planning. The purpose of this study was to evaluate the spatial and local uncertainty in predicting depletion rates of soil nutrients (P, K, Ca, and Mg) by soil erosion from green and burnt sugarcane harvesting scenarios, using sequential Gaussian simulation (SGS). A regular grid with equidistant intervals of 50 m (626 points) was established in the 200-ha study area, in Tabapuã, São Paulo, Brazil. The rate of soil depletion (SD) was calculated from the relation between the nutrient concentration in the sediments and the chemical properties in the original soil for all grid points. The data were subjected to descriptive statistical and geostatistical analysis. The mean SD rate for all nutrients was higher in the slash-and-burn than the green cane harvest scenario (Student’s t-test, p<0.05). In both scenarios, nutrient loss followed the order: Ca>Mg>K>P. The SD rate was highest in areas with greater slope. Lower uncertainties were associated to the areas with higher SD and steeper slopes. Spatial uncertainties were highest for areas of transition between concave and convex landforms.