Field hydraulic properties of an Alfisol under various fallow systems in southwestern Nigeria

The effects of various fallow management systems and cropping intensities on water infiltration were measured on an Alfisol at Ibadan in southwestern Nigeria. The objective was to determine the influence of the land use systems (a combination of crop-fallow sequences and intercropping types) on soil hydraulic properties obtained by disc permeameter and double-ring infiltration measurements. The experiment was established in 1989 as a split-plot design with four replications. The main plots were natural fallow, planted Pueraria phaseoloides and planted Leucaena leucocephala. The subplots were 1 year of maize/cassava intercrop followed by 3-year fallow (25% cropping intensity), or 2-year fallow (33% cropping intensity), or 1-year fallow (50% cropping intensity), or no fallow period (100% cropping intensity). Water infiltration rates and sorptivities were measured under saturated and unsaturated flow. Irrespective of land use, infiltration rates at the soil surface (121-324 cm h(-1)) were greater than those measured at 30 cm depth (55-144 cm h(-1)). This indicated that fewer large pores were present below 30 cm depth compared with 0-30 cm, depth. Despite some temporal variation, sorptivities with the highest mean value of 93.5 cm h(-1/2) increased as the cropping intensity decreased, suggesting a more continuous macropore system under less intensive land use systems. This was most likely due to continuous biopores created by perennial vegetation under long fallow systems. Intercropped maize and cassava yields also increased as cropping intensity decreased. The weak relationship between crop yields and hydraulic conductivity/infiltration rates suggests that the rates were not limiting.

Hydraulic properties of layered soils influence survival of Rhodes grass (Chloris gayana Kunth.) during water stress

Survival of vegetation on soil-capped mining wastes is often impaired during dry seasons due to the limited amount of water stored in the shallow soil capping. Growth and survival of Rhodes grass (Chloris gayana) during soil drying on various layered capping sequences constructed of combinations of topsoil, subsoil, seawater-neutralised residue sand and low grade bauxite was determined in a glasshouse. The aim was to describe the survival of Rhodes grass in terms of plant and soil water relationships. The soil water characteristic curve and soil texture analysis was a good predictor of plant survival. The combination of soil with a high water holding capacity and low soil water diffusivity (e.g. subsoil with high clay contents) with soil having a high water holding capacity and high diffusivity (e.g. residue sand) gave best survival during drying down (up to 88 days without water), whereas topsoil and low grade bauxite were unsuitable (plants died within 18-39 days). Clayey soil improved plant survival by triggering a water stress response during peak evaporative water demand once residue sand dried down and its diffusivity fell below a critical range. Thus, for revegetation in seasonally dry climates, soil capping should combine one soil with low diffusivity and one or more soils with high total water holding capacity and high diffusivity.

Analysis of soil wetting and solute transport in subsurface trickle irrigation

The increased use of trickle or drip irrigation is seen as one way of helping to improve the sustainability of irrigation systems around the world. However, soil water and solute transport properties and soil profile characteristics are often not adequately incorporated in the design and management of trickle systems. In this paper, we describe results of a simulation study designed to highlight the impacts of soil properties on water and solute transport from buried trickle emitters. The analysis addresses the influence of soil hydraulic properties, soil layering, trickle discharge rate, irrigation frequency, and timing of nutrient application on wetting patterns and solute distribution. We show that (1) trickle irrigation can improve plant water availability in medium and low permeability fine-textured soils, providing that design and management are adapted to account for their soil hydraulic properties, (2) in highly permeable coarse-textured soils, water and nutrients move quickly downwards from the emitter, making it difficult to wet the near surface zone if emitters are buried too deep, and (3) changing the fertigation strategy for highly permeable coarse-textured soils to apply nutrients at the beginning of an irrigation cycle can maintain larger amounts of nutrient near to and above the emitter, thereby making them less susceptible to leaching losses. The results demonstrate the need to account for differences in soil hydraulic properties and solute transport when designing irrigation and fertigation management strategies. Failure to do this will result in inefficient systems and lost opportunities for reducing the negative environmental impacts of irrigation.

Prediction of steady-state flux through variably saturated zones within a septic absorption trench

This paper describes effluent flow dynamics within a septic absorption system and the prediction of flow through the biomat and sub-biomat zone. Using soil hydraulic properties in a one dimensional model we demonstrate how soil hydraulic properties interact with biomat resistances to determine long-term acceptance rate (LTAR). The LTAR is a key parameter used in the Australian and New Zealand Standard AS1547:2000 to calculate the area of trench required to ensure trenches are not overloaded. Results show that several orders of magnitude variation in saturated hydraulic conductivity (Ks) collapse to a one order of magnitude variation in LTAR. These results are calculated from a model using basic flow theory, allowing LTAR to be estimated for any combination of biomat resistance and soil hydraulic properties. To increase the reliability of prediction of septic trench hydrology, HYDRUS 2D was used to model two dimensional flow. For more permeable soils, the exfiltration zone above sidewall biomat growth is shown to be a key pathway for excess effluent flow.

A technique to assess small-scale heterogeneity of chemical properties in soil aggregates

A new method is presented which allows the separation of the soil aggregate exterior from the aggregate core. The method employs a combination of aggregate freezing with rapid separation of aggregate exteriors using ultrasonic energy. The factors influencing the thickness of the removed aggregate surface layer include water content of the aggregate prior to freezing, temperature difference between that of the frozen aggregate and that of the liquid it is submerged in during sonification, sonification time and energy, and the type of the immersion liquid. The success of the method and the thickness of the removed aggregate surface were examined using barium ( Ba2+) as a tracer. Barium ( as BaCl2) is rapidly absorbed by soil and is present at only very low levels in natural soils. Surface layers of 0.2 - 0.4 cm thickness were successfully removed from aggregates of 1 - 4 cm diameter. Two examples are given from soils in northern NSW to demonstrate the occurrence of small- scale heterogeneity in soil chemical properties. Compared with the surface fraction, a 4 - 7% higher calcium concentration was found in the core fraction of a clay loam soil ( Dermosol). Conversely, on a cracking clay soil ( Vertosol), atrazine concentration was around 15 times greater in the aggregate surface fractions compared with core fractions. Compared with the traditional estimation of soil chemical properties on homogenised bulk soil samples, it is suggested that separate analysis of aggregate surface and core fractions could provide useful additional information on the relationships between soil properties and environmental responses.

Heterogeneity of physico-chemical properties in structured soils and its consequences

Structured soils are characterized by the presence of inter- and intra-aggregate pore systems and aggregates, which show varying chemical, physical, and biological properties depending on the aggregate type and land use system. How far these aspects also affect the ion exchange processes and to what extent the interaction between the carbon distribution and kind of organic substances affect the internal soil strength as well as hydraulic properties like wettability are still under discussion. Thus, the objective of this research was to clarify the effect of soil aggregation on physical and chemical properties of structured soils at two scales: homogenized material and single aggregates. Data obtained by sequentially peeling off soil aggregates layers revealed gradients in the chemical composition from the aggregate surface to the aggregate core. In aggregates from long term untreated forest soils we found lower amounts of carbon in the external layer, while in arable soils the differentiation was not pronounced. However, soil aggregates originating from these sites exhibited a higher concentration of microbial activity in the outer aggregate layer and declined towards the interior. Furthermore, soil depth and the vegetation type affected the wettability. Aggregate strength depended. on water suction and differences in tillage treatments.

Efffluent flux prediction in variably saturated soil zones within a septic tank-soil absorption trench

The treatment and hydraulic mechanisms in a septic tank-soil absorption system ( SAS) are highly influenced by the clogging layer or biomat zone which develops on bottom and lower sidewall surfaces within the trench. Flow rates through the biomat and sub-biomat zones are governed largely by the biomat hydraulic properties (resistance and hydraulic conductivity) and the unsaturated hydraulic conductivity of the underlying soil. One- and 2-dimensional models were used to investigate the relative importance of sidewall and vertical flow rates and pathways in SAS. Results of 1-dimensional modelling show that several orders of magnitude variation in saturated hydraulic conductivity (Ks) reduce to a 1 order of magnitude variation in long-term flow rates. To increase the reliability of prediction of septic trench hydrology, HYDRUS-2D was used to model 2-dimensional flow. In the permeable soils, under high trench loading, effluent preferentially flowed in the upper region of the trench where no resistant biomat was present (the exfiltration zone). By comparison, flow was more evenly partitioned between the biomat zones and the exfiltration zones of the low permeability soil. An increase in effluent infiltration corresponded with a greater availability of exfiltration zone, rather than a lower resistance of biomat. Results of modelling simulations demonstrated the important role that a permeable A horizon may play in limiting surface surcharge of effluent under high trench hydraulic loading.

Long-term flow rates and biomat zone hydrology in soil columns receiving septic tank effluent

Soil absorption systems (SAS) are used commonly to treat and disperse septic tank effluent (STE). SAS can hydraulically fail as a result of the low permeable biomat zone that develops on the infiltrative surface. The objectives of this experiment were to compare the hydraulic properties of biomats grown in soils of different textures, to investigate the long-term acceptance rates (LTAR) from prolonged application of STE, and to assess if soils were of major importance in determining LTAR. The STE was applied to repacked sand, Oxisol and Vertisol soil columns over a period of 16 months, at equivalent hydraulic loading rates of 50, 35 and 8 L/m(2)/d, respectively Infiltration rates, soil matric potentials, and biomat hydraulic properties were measured either directly from the soil columns or calculated using established soil physics theory. Biomats 1 to 2 cm thick developed in all soils columns with hydraulic resistances of 27 to 39 d. These biomats reduced a 4 order of magnitude variation in saturated hydraulic conductivity (K.) between the soils to a one order of magnitude variation in LTAR. A relationship between biomat resistance and organic loading rate was observed in all soils. Saturated hydraulic conductivity influenced the rate and extent of biomat development. However, once the biomat was established, the LTAR was governed by the resistance of the biomat and the sub-biomat soil unsaturated flow regime induced by the biomat. Results show that whilst initial soil K. is likely to be important in the establishment of the biomat zone in a trench, LTAR is determined by the biomat resistance and the unsaturated soil hydraulic conductivity, not the K, of a soil. The results call into question the commonly used approach of basing the LTAR, and ultimately trench length in SAS, on the initial K, of soils. (c) 2006 Elsevier Ltd. All rights reserved.

A mathematical model for estimating the extent of solute- and water-flux heterogeneity in multiple sample percolation experiments

Multiple sampling is widely used in vadose zone percolation experiments to investigate the extent in which soil structure heterogeneities influence the spatial and temporal distributions of water and solutes. In this note, a simple, robust, mathematical model, based on the beta-statistical distribution, is proposed as a method of quantifying the magnitude of heterogeneity in such experiments. The model relies on fitting two parameters, alpha and zeta to the cumulative elution curves generated in multiple-sample percolation experiments. The model does not require knowledge of the soil structure. A homogeneous or uniform distribution of a solute and/or soil-water is indicated by alpha = zeta = 1, Using these parameters, a heterogeneity index (HI) is defined as root 3 times the ratio of the standard deviation and mean. Uniform or homogeneous flow of water or solutes is indicated by HI = 1 and heterogeneity is indicated by HI > 1. A large value for this index may indicate preferential flow. The heterogeneity index relies only on knowledge of the elution curves generated from multiple sample percolation experiments and is, therefore, easily calculated. The index may also be used to describe and compare the differences in solute and soil-water percolation from different experiments. The use of this index is discussed for several different leaching experiments. (C) 1999 Elsevier Science B.V. All rights reserved.

Linking physiological processes with mangrove forest structure: phosphorus deficiency limits canopy development, hydraulic conductivity and photosynthetic carbon gain in dwarf Rhizophora mangle

Spatial gradients in mangrove tree height in barrier islands of Belize are associated with nutrient deficiency and sustained flooding in the absence of a salinity gradient. While nutrient deficiency is likely to affect many parameters, here we show that addition of phosphorus (P) to dwarf mangroves stimulated increases in diameters of xylem vessels, area of conductive xylem tissue and leaf area index (LAI) of the canopy. These changes in structure were consistent with related changes in function, as addition of P also increased hydraulic conductivity (K-s), stomatal conductance and photosynthetic assimilation rates to the same levels measured in taller trees fringing the seaward margin of the mangrove. Increased xylem vessel size and corresponding enhancements in stern hydraulic conductivity in P fertilized dwarf trees came at the cost of enhanced midday loss of hydraulic conductivity and was associated with decreased assimilation rates in the afternoon. Analysis of trait plasticity identifies hydraulic properties of trees as more plastic than those of leaf structural and physiological characteristics, implying that hydraulic properties are key in controlling growth in mangroves. Alleviation of P deficiency, which released trees from hydraulic limitations, reduced the structural and functional distinctions between dwarf and taller fringing tree forms of Rhizophora mangle.

Facilitating soil-and plant-water research through the use of TDR

The development of TDR for measurement of soil water content and electrical conductivity has resulted in a large shift in measurement methods for a breadth of soil and hydrological characterization efforts. TDR has also opened new possibilities for soil and plant research. Five examples show how TDR has enhanced our ability to conduct our soil- and plant-water research. (i) Oxygen is necessary for healthy root growth and plant development but quantitative evaluation of the factors controlling oxygen supply in soil depends on knowledge of the soil water content by TDR. With water content information we have modeled successfully some impact of tillage methods on oxygen supply to roots and their growth response. (ii) For field assessment of soil mechanical properties influencing crop growth, water content capability was added to two portable soil strength measuring devices; (a) A TDT (Time Domain Transmittivity)-equipped soil cone penetrometer was used to evaluate seasonal soil strengthwater content relationships. In conventional tillage systems the relationships are dynamic and achieve the more stable no-tillage relationships only relatively late in each growing season; (b) A small TDR transmission line was added to a modified sheargraph that allowed shear strength and water content to be measured simultaneously on the same sample. In addition, the conventional graphing procedure for data acquisition was converted to datalogging using strain gauges. Data acquisition rate was improved by more than a factor of three with improved data quality. (iii) How do drought tolerant plants maintain leaf water content? Non-destructive measurement of TDR water content using a flat serpentine triple wire transmission line replaces more lengthy procedures of measuring relative water content. Two challenges remain: drought-stressed leaves alter salt content, changing electrical conductivity, and drought induced changes in leaf morphology affect TDR measurements. (iv) Remote radar signals are reflected from within the first 2 cm of soil. Appropriate calibration of radar imaging for soil water content can be achieved by a parallel pair of blades separated by 8 cm, reaching 1.7 cm into soil and forming a 20 cm TDR transmission line. The correlation between apparent relative permittivity from TDR and synthetic aperture radar (SAR) backscatter coefficient was 0.57 from an airborne flyover. These five examples highlight the diversity in the application of TDR in soil and plant research.

Traffic and residue cover effects on infiltration

Wheel traffic can lead to compaction and degradation of soil physical properties. This study, as part of a study of controlled traffic farming, assessed the impact of compaction from wheel traffic on soil that had not been trafficked for 5 years. A tractor of 40 kN rear axle weight was used to apply traffic at varying wheelslip on a clay soil with varying residue cover to simulate effects of traffic typical of grain production operations in the northern Australian grain belt. A rainfall simulator was used to determine infiltration characteristics. Wheel traffic significantly reduced time to ponding, steady infiltration rate, and total infiltration compared with non-wheeled soil, with or without residue cover. Non-wheeled soil had 4-5 times greater steady infiltration rate than wheeled soil, irrespective of residue cover. Wheelslip greater than 10% further reduced steady infiltration rate and total infiltration compared with that measured for self-propulsion wheeling (3% wheelslip) under residue-protected conditions. Where there was no compaction from wheel traffic, residue cover had a greater effect on infiltration capacity, with steady infiltration rate increasing proportionally with residue cover (R-2 = 0.98). Residue cover, however, had much less effect on infiltration when wheeling was imposed. These results demonstrated that the infiltration rate for the non-wheeled soil under a controlled traffic zero-till system was similar to that of virgin soil. However, when the soil was wheeled by a medium tractor wheel, infiltration rate was reduced to that of long-term cropped soil. These results suggest that wheel traffic, rather than tillage and cropping, might be the major factor governing infiltration. The exclusion of wheel traffic under a controlled traffic farming system, combined with conservation tillage, provides a way to enhance the sustainability of cropping this soil for improved infiltration, increased plant-available water, and reduced runoff-driven soil erosion.

Anomalous water absorption in porous materials

The absorption of fluid by unsaturated, rigid porous materials may be characterized by the sorptivity. This is a simple parameter to determine and is increasingly being used as a measure of a material's resistance to exposure to fluids (especially moisture and reactive solutes) in aggressive environments. The complete isothermal absorption process is described by a nonlinear diffusion equation, with the hydraulic diffusivity being a strongly nonlinear function of the degree of saturation of the material. This diffusivity can be estimated from the sorptivity test. In a typical test the cumulative absorption is proportional to the square root of time. However, a number of researchers have observed deviation from this behaviour when the infiltrating fluid is water and there is some potential for chemo-mechanical interaction with the material. In that case the current interpretation of the test and estimation of the hydraulic diffusivity is no longer appropriate. Kuntz and Lavallee (2001) discuss the anomalous behaviour and propose a non-Darcian model as a more appropriate physical description. We present an alternative Darcian explanation and theory that retrieves the earlier advantages of the simple sorptivity test in providing parametric information about the material's hydraulic properties and allowing simple predictive formulae for the wetting profile to be generated.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.