HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Ross, P. J. Articles by Bristow, K. L. Search for Related Content PubMed Articles by Ross, P. J. Articles by Bristow, K. L. Agricola Articles by Ross, P. J. Articles by Bristow, K. L.

Simulating Water Movement in Layered and Gradational Soils Using the Kirchhoff Transform

Div. of Soils, CSIRO, Private Mail Bag, PO Aitkenvale, Townsville, Qld. 4814, Australia

*Corresponding author.

ABSTRACT

Changes in soil hydraulic properties with depth are common in field soils and need to be accounted for in simulations of field soil-water balances. We show that this is possible using the Kirchhoff transform (K transform) when solving Richards' equation numerically, even though the K transform is not continuous across boundaries between layers with different hydraulic properties. The solution is achieved by dividing the soil into depth elements, using the K transform to compute flow within each element, and then using matric potential (which is continuous across boundaries) to couple the elements. Combination of the K transform with iterative, implicit, mass-conserving methods of solving the flow equations results in efficient solutions that take no more than a few seconds on a personal computer to simulate infiltration. This is achieved with errors of only a few percent in the amount and distribution of soil water. Illustrative simulations using hydraulic properties that represent different soil profiles are given. These include simulation of infiltration into a layered profile consisting of tilled soil with a crust overlying an undisturbed soil, and gradational profiles with hydraulic properties of a sand varying linearly with depth to those of a clay, and conversely. Development and dissipation of surface ponding and perched water tables during simulation runs presented no difficulties. Nor did use of the Campbell formulations for the soil hydraulic properties, despite the discontinuity at the wet end.

Received for publication September 25, 1989.

This article has been cited by other articles:

				P. J. Ross Modeling Soil Water and Solute Transport--Fast, Simplified Numerical Solutions Agron. J., November 1, 2003; 95(6): 1352 - 1361. [Abstract] [Full Text] [PDF]

				P. J. Gregory, L. P. Simmonds, and C. J. Pilbeam Soil Type, Climatic Regime, and the Response of Water Use Efficiency to Crop Management Agron. J., September 1, 2000; 92(5): 814 - 820. [Abstract] [Full Text]

				J.-P. Vandervaere, M. Vauclin, and D. E. Elrick Transient Flow from Tension Infiltrometers: I. The Two-Parameter Equation Soil Sci. Soc. Am. J., July 1, 2000; 64(4): 1263 - 1272. [Abstract] [Full Text]