Both empirical and process-simulation models are useful for evaluating the effects of management practices on environmental quality and crop yield. The use of these models is limited, however, because they need many soil property values as input. The first step towards modelling is the collection of input data. Soil properties can be highly variable spatially and temporally, and measuring them is time-consuming and expensive. Efficient methods, which consider the uncertainty and cost of measurements, for estimating soil hydraulic properties form the main thrust of this study. Hydraulic properties are affected by other soil physical, and chemical properties, therefore it is possible to develop empirical relations to predict them. This idea quantified is called a pedotransfer function. Such functions may be global or restricted to a country or region. The different classification of particle-size fractions used in Australia compared with other countries presents a problem for the immediate adoption of exotic pedotransfer functions. A database of Australian soil hydraulic properties has been compiled. Pedotransfer functions for estimating water-retention and saturated hydraulic conductivity from particle size and bulk density for Australian soil are presented. Different approaches for deriving hydraulic transfer functions have been presented and compared. Published pedotransfer functions were also evaluated, generally they provide a satisfactory estimation of water retention and saturated hydraulic conductivity depending on the spatial scale and accuracy of prediction. Several pedotransfer functions were developed in this study to predict water retention and hydraulic conductivity. The pedotransfer functions developed here may predict adequately in large areas but for site-specific applications local calibration is needed. There is much uncertainty in the input data, and consequently the transfer functions can produce varied outputs. Uncertainty analysis is therefore needed. A general approach to quantifying uncertainty is to use Monte Carlo methods. By sampling repeatedly from the assumed probability distributions of the input variables and evaluating the response of the model the statistical distribution of the outputs can be estimated. A modified Latin hypercube method is presented for sampling joint multivariate probability distributions. This method is applied to quantify the uncertainties in pedotransfer functions of soil hydraulic properties. Hydraulic properties predicted using pedotransfer functions developed in this study are also used in a field soil-water model to analyze the uncertainties in the prediction of dynamic soil-water regimes. The use of the disc permeameter in the field conventionally requires the placement of a layer of sand in order to provide good contact between the soil surface and disc supply membrane. The effect of sand on water infiltration into the soil and on the estimate of sorptivity was investigated. A numerical study and a field experiment on heavy clay were conducted. Placement of sand significantly increased the cumulative infiltration but showed small differences in the infiltration rate. Estimation of sorptivity based on the Philip's two term algebraic model using different methods was also examined. The field experiment revealed that the error in infiltration measurement was proportional to the cumulative infiltration curve. Infiltration without placement of sand was considerably smaller because of the poor contact between the disc and soil surface. An inverse method for predicting soil hydraulic parameters from disc permeameter data has been developed. A numerical study showed that the inverse method is quite robust in identifying the hydraulic parameters. However application to field data showed that the estimated water retention curve is generally smaller than the one obtained in laboratory measurements. Nevertheless the estimated near-saturated hydraulic conductivity matched the analytical solution quite well. Th author believes that the inverse method can give a reasonable estimate of soil hydraulic parameters. Some experimental and theoretical problems were identified and discussed. A formal analysis was carried out to evaluate the efficiency of the different methods in predicting water retention and hydraulic conductivity. The analysis identified the contribution of individual source of measurement errors to the overall uncertainty. For single measurements, the inverse disc-permeameter analysis is economically more efficient than using pedotransfer functions or measuring hydraulic properties in the laboratory. However, given the large amount of spatial variation of soil hydraulic properties it is perhaps not surprising that lots of cheap and imprecise measurements, e.g. by hand texturing, are more efficient than a few expensive precise ones.

The intensification of land-water management due to agriculture, forestry, and urbanization is a global phenomenon increasing the pressure on worlds water resources and threatening water security in North America. The Prairie Pothole Region of North America covers approximately 775,000 km2 and contains millions of wetlands that serve important hydrological and ecological functions. The unique hummocky topography and the variable effect of different processes contribute to high spatio-temporal variability in soil water, posing major challenges in hydrological studies. The objectives of this study were to a) examine the spatial pattern of soil water storage and its scale and location characteristics; and b) to identify its controls at multiple scales. Soil water content at 20 cm intervals down to 140 cm was measured along a transect extending over several knolldepression cycles in a hummocky landscape. High water storage in depressions and low water storage on the knolls created a spatial pattern that was inversely related to elevation. Spatial patterns were strongly similar within any given season (intra-season rank correlation coefficient as high as 0.99), moreso than between the same season over different years (inter-annual rank correlation coefficient as high as 0.97). Less similar spatial patterns were observed between different seasons (inter-season rank correlation coefficients as high as 0.90). While the intra-season and inter-annual spatial patterns were similar at scales >18 m, the inter-season spatial patterns were similar at much large scales (>72 m). This may be due to the variations in landform elements and micro-topography. The similarity at scales >72 m were present at any time and depth. However, small- and medium-scale spatial patterns changed with depth and with season due to a change in the hydrological processes. The relative dominance of a given set of processes operating both within a season and for the same season over different years yielded strong intra-season and inter-annual similarity at scales >18 m. Moreover, similarity was stronger with increasing depth, and was thought to be due to the dampening effect of overlying soil layers that are more dynamic. Similarity of spatial patterns over time helps to identify the location that best represents the field averaged soil water and improves sampling efficiency. Change in the similarity of scales of spatial pattern helps identify the change in sampling domain as controlled by hydrological processes. The scale information can be used to improve prediction for use in environmental management and modeling of different surface and subsurface hydrological processes. The similarity of spatial pattern between the surface and subsurface layers help make inferences on deep layer hydrological processes as well as groundwater dynamics from surface water measurements.

The intensification of land-water management due to agriculture, forestry, and urbanization is a global phenomenon increasing the pressure on worlds water resources and threatening water security in North America. The Prairie Pothole Region of North America covers approximately 775,000 km2 and contains millions of wetlands that serve important hydrological and ecological functions. The unique hummocky topography and the variable effect of different processes contribute to high spatio-temporal variability in soil water, posing major challenges in hydrological studies. The objectives of this study were to a) examine the spatial pattern of soil water storage and its scale and location characteristics; and b) to identify its controls at multiple scales. Soil water content at 20 cm intervals down to 140 cm was measured along a transect extending over several knolldepression cycles in a hummocky landscape. High water storage in depressions and low water storage on the knolls created a spatial pattern that was inversely related to elevation. Spatial patterns were strongly similar within any given season (intra-season rank correlation coefficient as high as 0.99), moreso than between the same season over different years (inter-annual rank correlation coefficient as high as 0.97). Less similar spatial patterns were observed between different seasons (inter-season rank correlation coefficients as high as 0.90). While the intra-season and inter-annual spatial patterns were similar at scales >18 m, the inter-season spatial patterns were similar at much large scales (>72 m). This may be due to the variations in landform elements and micro-topography. The similarity at scales >72 m were present at any time and depth. However, small- and medium-scale spatial patterns changed with depth and with season due to a change in the hydrological processes. The relative dominance of a given set of processes operating both within a season and for the same season over different years yielded strong intra-season and inter-annual similarity at scales >18 m. Moreover, similarity was stronger with increasing depth, and was thought to be due to the dampening effect of overlying soil layers that are more dynamic. Similarity of spatial patterns over time helps to identify the location that best represents the field averaged soil water and improves sampling efficiency. Change in the similarity of scales of spatial pattern helps identify the change in sampling domain as controlled by hydrological processes. The scale information can be used to improve prediction for use in environmental management and modeling of different surface and subsurface hydrological processes. The similarity of spatial pattern between the surface and subsurface layers help make inferences on deep layer hydrological processes as well as groundwater dynamics from surface water measurements.

This study examined the influence of different management practices under sugarcane production on soil chemical, biological and physical properties on a long-term (59yr) field experiment in KwaZulu-Natal. These management practices included conventional post-harvest burning of crop residues, with the tops either left on the soil surface or with tops removed, and green cane harvesting with the retention of crop residues on the soil surface as a trash blanket. Each of these treatments were either fertilized or did not receive fertilizer. The data collected was used to evaluate the effect of crop residue management on soil quality under sugarcane production. Soil organic matter content increased from 39 g C kg¯¹ soil, under conventional burning to 55 g C kg¯¹ soil in the surface 10 cm under green cane harvesting where crop residues are returned to the soil. It also resulted in greater recycling of nutrients and increases in exchangeable K and Ca and extractable P. Fertilizer application resulted in a build-up of soil nutrients in combination with trash retention. Fertilizer application induced increases in exchangeable P and also some accumulation in soil organic P. Exchangeable and non-exchangeable K concentrations were also higher in fertilized than non-fertilized treatments. However, nitrogen fertilizer application and, to a lesser degree, organic matter mineralization, resulted in soil acidification to a depth of 30 cm. Acidification in the fertilized treatments resulted in a concomitant increase in exchangeable acidity and exchangeable Al, due to the increase in H⁺ ions and solubilized Al species on exchange sites. Base cations moved into soil solution and were leaching to lower soil layers. The decrease in soil pH resulted in the surface charge conferred on the variable charge surfaces on soil colloids becoming less negative and as a result there was a decrease in ECEC. Acidification in fertilized treatment not only increased exchangeable Al but also the buffering reserve of non-exchangeable al; both that complexed with soil organic matter (CuCl₂- extractable) and that present as hydroxy - Al associated with mineral colloids (ammonium acetate - extractable). The increased organic matter content under trash retention resulted in an increase in ECEC. This enabled the soil to retain greater amounts of Ca²⁺, Mg²⁺ and K⁺ which were returned to the soil in the trash. Both residue retention and fertilizer application had a positive effect on the microbial biomass C and N and the microbial quotient increased from 0.39% to 0.86% as organic C increased from 39 g C kg¯¹ soil under burnt treatments to 55 g C kg¯¹ soil under trashed, fertilized treatments. This increase was associated with increased concentrations of labile organic material (K₂S0₄₋extractable) present as well as increased amounts of nutrients being cycled through the plant-soil system. The light fraction organic matter also increased with increasing returns of organic residues. However, the large active microbial biomass under the trashed, fertilized treatment resulted in an increased turnover rate of this fraction and consequently resulted in lower LF dry matter, C and N than in the unfertilized treatment. Aggregate stability increased with increasing amounts of organic material returned due to trash retention. Nevertheless, fertilized treatments induced a lower aggregate stability than unfertilized ones, despite the tendancy for the latter to have higher organic C and microbial biomass values. This was attributed to an increase in the proportion of exchangeable cations present in monovalent form (due to application of fertilizer K and leaching of Ca and Mg) favouring dispersion and a decline in aggregate stability. Green cane harvesting resulted in an increase in microbial activity (basal respiration, FDA hydrolytic activity, arginine ammonification rate and dehydrogenase activity) and in the activity of specific soil enzymes involved in turnover of C, N, P and S to a depth of 30 cm. Increased activities of these enzymes reflect a higher rate of turnover of C, N, P and S. The metabolic quotient decreased with increasing residue return, indicating a more metabolically efficient microbial community. Fertilizer application resulted in a variable effect on enzyme activity. Long-term fertilizer application resulted in an increase in invertase and acid phosphatase, a decrease in L-histidase and arylsulphatase and had little effect on protease and alkaline phosphatase. These variable effects were explained in terms of an interaction between fertilizer - induced increases in Corg and soil nutrient status and fertilizer - induced soil acidification. The size and activity of the soil microbial biomass was studied in the plant row and in the inter-row of a sugarcane field under burning or green cane harvesting. Soils were sampled to 30 cm depth in (i) the centre of the plant row, (ii) 30 cm out from the row centre and (iii) 60 cm out from the row centre (i.e. the middle of the inter-row area). Under burning, the only substantial input of organic matter to the soil was from root turnover in the row area where the root biomass was concentrated. As a consequence, the size (microbial biomass C) and activity (basal respiration of the soil microbial community were concentrated in the row. However, under green cane harvesting there was a large input of organic matter in the inter-row area in the form of the trash blanket itself and through turnover of crop roots that were concentrated in the surface 10 cm of the soil below the blanket. As a result, soil microbial activity was considerably higher in the interrow area under green cane harvesting than under burning. Phospholipids are essential membrane components of microorganisms and a good correlation was found between the total PLFA's extracted from soils and the microbial biomass C, indicating that phospholipids are an accurate measurement of living biomass. Mutivariate statistical analysis (PCA) was used to separate different PLFA profiles under burning versus trash retention and under different land uses (sugarcane, maize, annual and perennial pasture and undisturbed veld). Soil organic matter content contributed the greatest variance in the data along the first axis. That is increasing soil organic matter return not only increased the size of the microbial biomass, but also affected the composition of the microbial community. There was a shift in the different sub-fractions under different management practices. MUFA's are general biomarkers of Gram negative bacteria and were found to be a sensitive indicator of higher substrate availability (i.e they increased under green cane harvesting). Fungal biomarkers indicated an increased fungal biomass associated with surface application of residues. Soil physical conditions were considered to be a contributing factor to the shift in microbial community structure. Increased organic matter content improved soil physical conditions and preferentially stimulated the growth of aerobic microorganisms. In addition to this, the proportion of SATFA (gram positive bacteria) was found to increase in response to burning. This increase was attributed to the survival mechanisms of these microorganisms (i.e. endospore formation). It was found that the conversion from burning to trash management changes the composition of the soil microbial community. The effect of management practices on soil functional diversity was also evaluated using two methods (i.e. Biolog plates and substrate induced respiration (SIR)). Biolog plates are a selective technique that stimulate growth of a small proportion of the soil microbial community whereas the SIR technique measures the activity of the metabolically active microbial community in situ. As a result the SIR method separated treatments more effectively than Biolog plates (i.e. annually tilled treatments, permanent grassland sites and fertilized and unfertilized treatments). The quantity and the quality of organic C supply influenced the catabolic diversity. Conversion from burning to green cane harvesting greatly increased catabolic evenness and richness and therefore presumably also tended to increase the resilience of the soil to stress and disturbance particularly in relation to decomposition functions. It was concluded that conversion from preharvest burning to green cane harvesting results in an increase in soil organic matter content, an improvement in soil structure and soil nutrient status, an increase in the size, activity, taxonomic and functional diversity of the soil microbial community. The practice should therefore be promoted to the South African sugar industry. / Thesis (Ph.D.)-University of Natal, Pietermaritzburg, 2003.

There is evidence that soil compaction, through the use of mechanised equipment, causes detrimental effects to soil quality and reduces long-term productivity of soils. For economic reasons, farmers need to purchase larger, heavier machinery in order to cultivate larger areas under crops, resulting in larger forces on the soil. The severity of soil compaction is governed by various soil and vehicle properties and normally causes an increase in the soil's bulk density and a decrease in the air filled porosity. These changes in soil properties have negative effects on crop production and environmental sustainability. The aim of this study was to investigate and develop a model based decision support system for soil compaction management and research. Soil compaction occurs during the transfer of stresses from the tyre interface into the soil. Numerically, it has been modelled using both mechanistic and empirical models, which attempt to simulate the stress propagation and also sometimes the consequent damage to the soil. The SOCOMO soil compaction model is described and this model computes the stress at a point in the soil for any given horizontal and vertical stress distribution at the soil / tyre interface. It has been successfully used in the Netherlands and in Sweden to map the impact on the soil. The SOCOMO model was tested and verified at a forestry site in Richmond, KwaZulu- Natal. Relationships to determine bulk density were also tested and verified. The SOCOMO model performs satisfactory (RMSE = 47.9 kPa), although it tends to overestimate the pressures within the soil. This could be as a result of the high organic carbon content in the particular soil. Models predicting bulk density also performed satisfactory (RMSE = 69.9 kg.m" ), but resultant densities in the soil are generally underestimated. Future research is needed to find better relationships to estimate changes in dry bulk density and to test the model on a wider range of soils. If the model performs satisfactory it could provide a useful tool to determine the impact of soil compaction on crop yield. / Thesis (M.Sc.Eng.)-University of KwaZulu-Natal, 2006.

The Umzimvubu catchment on the eastern coastal escarpment of South Africa is sensitive to anthropogenic influences,with commercial and subsistence agriculture, irrigation, domestic and rural settlements and forestry compete for water use. An adequate supply of water to the region is seen as imperative in the light of the recent establishment of forest cultivation. In order to provide a sound assessment of the impacts of afforestation on the catchment, the subsurface hydrological processes of hillslopes on the Molteno sedimentary formations of the region must be clearly understood. Since the runoff hydrograph is, to a large degree, dependent on the subsurface processes, a number of models that simulate small catchment runoff have been developed. However, recent successful application of tracer techniques to hydrological modelling has shown that the subsurface processes are still not fully understood (Schultz, 1999), and whether or not the subsurface processes are modelled adequately is most often not verified, since there is a lack of relevant data. It is, therefore imperative that the subsurface component of these small catchment runoff models be improved. This can be achieved by first observing detailed subsurface water dynamics and assessing these against the catchment runoff response. In this dissertation, results from a detailed experiment that was initiated in a 1.5 km2 catchment in the northern East Cape Province are shown. Nests of automated tensiometers, groundwater level recorders and weather stations have been placed at critical points around the catchment, and these , together with soil hydraulic and physical characteristics are used to define and identify the dominant hillslope processes. Two crump weirs record runoff from these hillslopes. The results of this subsurface study highlight the dynamics of surface and subsurface water in the hillslope transects. It is evident that the subsurface processes are strongly influenced by the -bedrock topography as well as the soil characteristics, such as macropore flow and deep percolation. Using the monitored data and 2-D vadose zone modelling, the dominant hillslope processes have been defined and are used to aid in the selection of critical parameters to be used in estimating the catchment runoff. Results show that a clear understanding of the subsurface dynamics can lead to a realistic estimation of catchment scale runoff response. / Thesis (M.Sc.)-University of Natal, Pietermaritzburg, 2000.

Forest ecosystems undergo both natural and human induced disturbances. Depending on disturbance type, soil physical and chemical parameters show different response patterns during the recovery phase. An added level of complication is the ecological site types occurring throughout a forested area. The identification of indicators of soil fertility and the successful emulation of a natural disturbance regime were the scope of this research. / The research presented herein took place in the Manitoba Model Forest (MBMF), located in eastern Manitoba, where the natural disturbance regime is wildfire. Timber harvest strategies used in the area are designed to emulate a wildfire (5% retention of standing timber and extensive slash inputs) and clearcut harvesting. The objective of this study was to document changes in forest floor and soil properties prior to and following harvesting, and to compare these properties to those found in a small wildfire that burned in the MBMF in late summer 1998, which serves as a benchmark to the harvest. In the two study areas, both thin mineral soil (5--20 cm) and moderately deep mineral soil (20--100 cm) ecosite types were monitored at four dates over a two year period for soil physical and nutritional response patterns.

Both empirical and process-simulation models are useful for evaluating the effects of management practices on environmental quality and crop yield. The use of these models is limited, however, because they need many soil property values as input. The first step towards modelling is the collection of input data. Soil properties can be highly variable spatially and temporally, and measuring them is time-consuming and expensive. Efficient methods, which consider the uncertainty and cost of measurements, for estimating soil hydraulic properties form the main thrust of this study. Hydraulic properties are affected by other soil physical, and chemical properties, therefore it is possible to develop empirical relations to predict them. This idea quantified is called a pedotransfer function. Such functions may be global or restricted to a country or region. The different classification of particle-size fractions used in Australia compared with other countries presents a problem for the immediate adoption of exotic pedotransfer functions. A database of Australian soil hydraulic properties has been compiled. Pedotransfer functions for estimating water-retention and saturated hydraulic conductivity from particle size and bulk density for Australian soil are presented. Different approaches for deriving hydraulic transfer functions have been presented and compared. Published pedotransfer functions were also evaluated, generally they provide a satisfactory estimation of water retention and saturated hydraulic conductivity depending on the spatial scale and accuracy of prediction. Several pedotransfer functions were developed in this study to predict water retention and hydraulic conductivity. The pedotransfer functions developed here may predict adequately in large areas but for site-specific applications local calibration is needed. There is much uncertainty in the input data, and consequently the transfer functions can produce varied outputs. Uncertainty analysis is therefore needed. A general approach to quantifying uncertainty is to use Monte Carlo methods. By sampling repeatedly from the assumed probability distributions of the input variables and evaluating the response of the model the statistical distribution of the outputs can be estimated. A modified Latin hypercube method is presented for sampling joint multivariate probability distributions. This method is applied to quantify the uncertainties in pedotransfer functions of soil hydraulic properties. Hydraulic properties predicted using pedotransfer functions developed in this study are also used in a field soil-water model to analyze the uncertainties in the prediction of dynamic soil-water regimes. The use of the disc permeameter in the field conventionally requires the placement of a layer of sand in order to provide good contact between the soil surface and disc supply membrane. The effect of sand on water infiltration into the soil and on the estimate of sorptivity was investigated. A numerical study and a field experiment on heavy clay were conducted. Placement of sand significantly increased the cumulative infiltration but showed small differences in the infiltration rate. Estimation of sorptivity based on the Philip's two term algebraic model using different methods was also examined. The field experiment revealed that the error in infiltration measurement was proportional to the cumulative infiltration curve. Infiltration without placement of sand was considerably smaller because of the poor contact between the disc and soil surface. An inverse method for predicting soil hydraulic parameters from disc permeameter data has been developed. A numerical study showed that the inverse method is quite robust in identifying the hydraulic parameters. However application to field data showed that the estimated water retention curve is generally smaller than the one obtained in laboratory measurements. Nevertheless the estimated near-saturated hydraulic conductivity matched the analytical solution quite well. Th author believes that the inverse method can give a reasonable estimate of soil hydraulic parameters. Some experimental and theoretical problems were identified and discussed. A formal analysis was carried out to evaluate the efficiency of the different methods in predicting water retention and hydraulic conductivity. The analysis identified the contribution of individual source of measurement errors to the overall uncertainty. For single measurements, the inverse disc-permeameter analysis is economically more efficient than using pedotransfer functions or measuring hydraulic properties in the laboratory. However, given the large amount of spatial variation of soil hydraulic properties it is perhaps not surprising that lots of cheap and imprecise measurements, e.g. by hand texturing, are more efficient than a few expensive precise ones.

Thesis (MScFor)--University of Stellenbosch, 2007. / Bibliography. Also available via the Internet.

Thesis (M.Sc.) - University of KwaZulu-Natal, Pietermaritzburg, 2009. / Full text also available online. Scroll down for electronic link.