Estimation of soil hydraulic properties and their uncertainty: comparison between laboratory and field experiment

Often the soil hydraulic parameters are obtained by the inversion of measured data (e.g. soil moisture, pressure head, and cumulative infiltration, etc.). However, the inverse problem in unsaturated zone is ill-posed due to various reasons, and hence the parameters become non-unique. The presence of multiple soil layers brings the additional complexities in the inverse modelling. The generalized likelihood uncertainty estimate (GLUE) is a useful approach to estimate the parameters and their uncertainty when dealing with soil moisture dynamics which is a highly non-linear problem. Because the estimated parameters depend on the modelling scale, inverse modelling carried out on laboratory data and field data may provide independent estimates. The objective of this paper is to compare the parameters and their uncertainty estimated through experiments in the laboratory and in the field and to assess which of the soil hydraulic parameters are independent of the experiment. The first two layers in the field site are characterized by Loamy sand and Loamy. The mean soil moisture and pressure head at three depths are measured with an interval of half hour for a period of 1 week using the evaporation method for the laboratory experiment, whereas soil moisture at three different depths (60, 110, and 200 cm) is measured with an interval of 1 h for 2 years for the field experiment. A one-dimensional soil moisture model on the basis of the finite difference method was used. The calibration and validation are approximately for 1 year each. The model performance was found to be good with root mean square error (RMSE) varying from 2 to 4 cm(3) cm(-3). It is found from the two experiments that mean and uncertainty in the saturated soil moisture (theta(s)) and shape parameter (n) of van Genuchten equations are similar for both the soil types. Copyright (C) 2010 John Wiley & Sons, Ltd.

Developing a reliable strategy to infer the effective soil hydraulic properties from field evaporation experiments for agro-hydrological models

The Richards equation has been widely used for simulating soil water movement. However, the take-up of agro-hydrological models using the basic theory of soil water flow for optimizing irrigation, fertilizer and pesticide practices is still low. This is partly due to the difficulties in obtaining accurate values for soil hydraulic properties at a field scale. Here, we use an inverse technique to deduce the effective soil hydraulic properties, based on measuring the changes in the distribution of soil water with depth in a fallow field over a long period, subject to natural rainfall and evaporation using a robust micro Genetic Algorithm. A new optimized function was constructed from the soil water contents at different depths, and the soil water at field capacity. The deduced soil water retention curve was approximately parallel but higher than that derived from published pedo-tranfer functions for a given soil pressure head. The water contents calculated from the deduced soil hydraulic properties were in good agreement with the measured values. The reliability of the deduced soil hydraulic properties was tested in reproducing data measured from an independent experiment on the same soil cropped with leek. The calculation of root water uptake took account for both soil water potential and root density distribution. Results show that the predictions of soil water contents at various depths agree fairly well with the measurements, indicating that the inverse analysis is an effective and reliable approach to estimate soil hydraulic properties, and thus permits the simulation of soil water dynamics in both cropped and fallow soils in the field accurately. (C) 2009 Elsevier B.V. All rights reserved.

Comparison among monitoring strategies to assess water flow dynamic and soil hydraulic properties in agricultural soils

Irrigated agriculture is usually performed in semi-arid regions despite scarcity of water resources. Therefore, optimal irrigation management by monitoring the soil is essential, and assessing soil hydraulic properties and water flow dynamics is presented as a first measure. For this purpose, the control of volumetric water content, θ, and pressure head, h, is required. This study adopted two types of monitoring strategies in the same experimental plot to control θ and h in the vadose zone: i) non-automatic and more time-consuming; ii) automatic connected to a datalogger. Water flux was modelled with Hydrus-1D using the data collected from both acquisition strategies independently (3820 daily values for the automatic; less than 1000 for the non-automatic). Goodness-of-fit results reported a better adjustment in case of automatic sensors. Both model outputs adequately predicted the general trend of θ and h, but with slight differences in computed annual drainage (711 mm and 774 mm). Soil hydraulic properties were inversely estimated from both data acquisition systems. Major differences were obtained in the saturated volumetric water content, θs, and the n and α van Genuchten model shape parameters. Saturated hydraulic conductivity, Ks, shown lower variability with a coefficient of variation range from 0.13 to 0.24 for the soil layers defined. Soil hydraulic properties were better assessed through automatic data acquisition as data variability was lower and accuracy was higher.

Applying HYDRUS to flow in a sodic clay soil with solution composition–dependent hydraulic conductivity

Management of sodic soils under irrigation often requires application of chemical ameliorants to improve permeability combined with leaching of excess salts. Modeling irrigation, soil treatments, and leaching in these sodic soils requires a model that can adequately represent the physical and chemical changes in the soil associated with the amelioration process. While there are a number of models that simulate reactive solute transport, UNSATCHEM and HYDRUS-1D are currently the only models that also include an ability to simulate the impacts of soil chemistry on hydraulic conductivity. Previous researchers have successfully applied these models to simulate amelioration experiments on a sodic loam soil. To further gauge their applicability, we extended the previous work by comparing HYDRUS simulations of sodic soil amelioration with the results from recently published laboratory experiments on a more reactive, repacked sodic clay soil. The general trends observed in the laboratory experiments were able to be simulated using HYDRUS. Differences between measured and simulated results were attributed to the limited flexibility of the function that represents chemistry-dependent hydraulic conductivity in HYDRUS. While improvements in the function could be made, the present work indicates that HYDRUS-UNSATCHEM captures the key changes in soil hydraulic properties that occur during sodic clay soil amelioration and thus extends the findings of previous researchers studying sodic loams.

Impacts of sodic soil amelioration on hydraulic conductivity and deep drainage in the Lower Burdekin

An understanding of the influence of soil chemistry on soil hydraulic properties is of critical importance for the management of sodic soils under irrigation. The hydraulic conductivity of sodic soils has been shown to be affected by properties of the applied solution including pH (Suarez et al. 1984), sodicity and salt concentration (McNeal and Coleman 1966). The changes in soil hydraulic conductivity are the result of changes in the spacing between clay layers in response to changes in soil solution chemistry. While the importance o f soil chemistry in controlling hydraulic conductivity is known, the exact impacts of sodic soil amelioration on hydraulic conductivity and deep drainage at a given location are difficult to predict. This is because the relationships between soil chemical factors and hydraulic conductivity are soil specific and because local site specific factors also need to be considered to determine the actual impacts on deep drainage rates.

Cement-fly ash stabilisation/solidification of contaminated soil: Performance properties and initiation of operating envelopes

This study was aimed at evaluating the mechanical and pH-dependent leaching performance of a mixed contaminated soil treated with a mixture of Portland cement (CEMI) and pulverised fuel ash (PFA). It also sought to develop operating envelopes, which define the range(s) of operating variables that result in acceptable performance. A real site soil with low contaminant concentrations, spiked with 3000mg/kg each of Cd, Cu, Pb, Ni and Zn, and 10,000mg/kg of diesel, was treated with one part CEMI and four parts PFA (CEMI:PFA=1:4) using different binder and water contents. The performance was assessed over time using unconfined compressive strength (UCS), hydraulic conductivity, acid neutralisation capacity (ANC) and pH-dependent leachability of contaminants. With binder dosages ranging from 5% to 20% and water contents ranging from 14% to 21% dry weight, the 28-day UCS was up to 500kPa and hydraulic conductivity was around 10-8m/s. With leachant pH extremes of 7.2 and 0.85, leachability of the contaminants was in the range: 0.02-3500mg/kg for Cd, 0.35-1550mg/kg for Cu, 0.03-92mg/kg for Pb, 0.01-3300mg/kg for Ni, 0.02-4010mg/kg for Zn, and 7-4884mg/kg for total petroleum hydrocarbons (TPHs), over time. Design charts were produced from the results of the study, which show the water and/or binder proportions that could be used to achieve relevant performance criteria. The charts would be useful for the scale-up and design of stabilisation/solidification (S/S) treatment of similar soil types impacted with the same types of contaminants. © 2013 Elsevier Ltd.

Linking watershed terrain and hydrology to soil chemical properties, microbial communities and impacts on soil organic C in a humid mid-latitude forested watershed

Understanding the response of humid mid-latitude forests to changes in precipitation, temperature, nutrient cycling, and disturbance is critical to improving our predictive understanding of changes in the surface-subsurface energy balance due to climate change. Mechanistic understanding of the effects of long-term and transient moisture conditions are needed to quantify
linkages between changing redox conditions, microbial activity, and soil mineral and nutrient interactions on C cycling and greenhouse gas releases. To illuminate relationships between the soil chemistry, microbial communities and organic C we established transects across hydraulic and topographic gradients in a small watershed with transient moisture conditions. Valley bottoms tend to be more frequently saturated than ridge tops and side slopes which generally are only saturated when shallow storm flow zones are active. Fifty shallow (~36”) soil cores were collected during timeframes representative of low CO2, soil winter conditions and high CO2, soil summer conditions. Cores were subdivided into 240 samples based on pedology and analyses of the geochemical (moisture content, metals, pH, Fe species, N, C, CEC, AEC) and microbial (16S rRNA gene
amplification with Illumina MiSeq sequencing) characteristics were conducted and correlated to watershed terrain and hydrology. To associate microbial metabolic activity with greenhouse gas emissions we installed 17 soil gas probes, collected gas samples for 16 months and analyzed them for CO2 and other fixed and greenhouse gasses. Parallel to the experimental efforts our data is being used to support hydrobiogeochemical process modeling by coupling the Community Land Model (CLM) with a subsurface process model (PFLOTRAN) to simulate processes and interactions from the molecular to watershed scales. Including above ground processes (biogeophysics, hydrology, and vegetation dynamics), CLM provides mechanistic water, energy, and organic matter inputs to the surface/subsurface models, in which coupled biogeochemical reaction
networks are used to improve the representation of below-ground processes. Preliminary results suggest that inclusion of above ground processes from CLM greatly improves the prediction of moisture response and water cycle at the watershed scale.

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.

Impacts of sodic soil amelioration on deep drainage

Groundwater tables are rising beneath irrigated fields in some areas of the Lower Burdekin in North Queensland, Australia. The soils where this occurs are predominantly sodic clay soils with low hydraulic conductivities. Many of these soils have been treated by applying gypsum or by increasing the salinity of irrigation water by mixing saline groundwater with fresh river water. While the purpose of these treatments is to increase infiltration into the surface soils and improve productivity of the root zone, it is thought that the treatments may have altered the soil hydraulic properties well below the root zone leading to increased groundwater recharge and rising water tables. In this paper we discuss the use of column experiments and HYDRUS modelling, with major ion reaction and transport and soil water chemistry-dependent hydraulic conductivity, to assess the likely depth, magnitude and timing of the impacts of surface soil amelioration on soil hydraulic properties below the root zone and hence groundwater recharge. In the experiments, columns of sodic clays from the Lower Burdekin were leached for extended periods of time with either gypsum solutions or mixed cation salt solutions and change s in hydraulic conductivity were measured. Leaching with a gypsum solution for an extended time period, until the flow rate stabilised, resulted in an approximately twenty fold increase in the hydraulic conductivity when compared with a low salinity, mixed cation solution. HYDRUS modelling was used to high light the role of those factors which might influence the impacts of soil treatment, particularly at depth, including the large amounts of rain during the relatively short wet season and the presence of thick low permeability clay layers.

Lantana camara L. (Verbenaceae) invasion effects on soil physicochemical properties

Lantana camara is a recognized weed of worldwide significance due to its extensive distribution and its impacts on primary industries and nature conservation. However, quantitative data on the impact of the weed on soil ecosystem properties are scanty, especially in SE Australia, despite the pervasive presence of the weed along its coastal and inland regions. Consequently, mineral soils for physicochemical analyses were collected beneath and away from L. camara infestations in four sites west of Brisbane, SE Australia. These sites (hoop pine plantation, cattle farm, and two eucalyptus forests with occasional grazing and a fire regime, respectively) vary in landscape and land-use types. Significant site effect was more frequently observed than effect due to invasion status. Nonetheless, after controlling for site differences, ~50% of the 23 soil traits examined differed significantly between infested and non-infested soils. Moisture, pH, Ca, total and organic C, and total N (but not exchangeable N in form of NO3-) were significantly elevated, while sodium, chloride, copper, iron, sulfur, and manganese, many of which can be toxic to plant growth if present in excess levels, were present at lower levels in soils supporting L. camara compared to soils lacking the weed. These results indicate that L. camara can improve soil fertility and influence nutrient cycling, making the substratum ideal for its own growth and might explain the ability of the weed to outcompete other species, especially native ones.

Cat's claw creeper vine, Macfadyena unguis-cati (Bignoniaceae), invasion impacts: comparative leaf nutrient content and effects on soil physicochemical properties

Macfadyena unguis-cati (L.) Gentry (Bignoniaceae) is a major environmental weed in coastal Queensland, Australia. There is a lack of quantitative data on its leaf chemistry and its impact on soil properties. Soils from infested vs uninfested areas, and leaves of M. unguis-cati and three co-occurring vine species (one exotic, two native) were collected at six sites (riparian and non-riparian) in south-eastern Queensland. Effects of invasion status, species, site and habitat type were examined using univariate and multivariate analyses. Habitat type had a greater effect on soil nutrients than on leaf chemistry. Invasion effect of M. unguis-cati on soil chemistry was more pronounced in non-riparian than in riparian habitat. Significantly higher values were obtained in M. unguis-cati infested (vs. uninfested) soils for ~50% of traits. Leaf ion concentrations differed significantly between exotic and native vines. Observed higher leaf-nutrient load (especially nitrogen, phosphorus and potassium) in exotic plants aligns with the preference of invasive plant species for disturbed habitats with higher nutrient input. Higher load of trace elements (aluminium, boron, cadmium and iron) in its leaves suggests that cycling of heavy-metal ions, many of which are potentially toxic at excess level, could be accelerated in soils of M. unguis-cati-invaded landscape. Although inferences from the present study are based on correlative data, the consistency of the patterns across many sites suggests that M. unguis-cati may improve soil fertility and influence nutrient cycling, perhaps through legacy effects of its own litter input.

Soil hydrological properties and conditions, site preparation, and the long-term performance of planted Scots pine (Pinus sylvestris L.) on upland forest sites in Finnish Lapland

Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) Karst.) forests dominate in Finnish Lapland. The need to study the effect of both soil factors and site preparation on the performance of planted Scots pine has increased due to the problems encountered in reforestation, especially on mesic and moist, formerly spruce-dominated sites. The present thesis examines soil hydrological properties and conditions, and effect of site preparation on them on 10 pine- and 10 spruce-dominated upland forest sites. Finally, the effects of both the site preparation and reforestation methods, and soil hydrology on the long-term performance of planted Scots pine are summarized. The results showed that pine and spruce sites differ significantly in their soil physical properties. Under field capacity or wetter soil moisture conditions, planted pines presumably suffer from excessive soil water and poor soil aeration on most of the originally spruce sites, but not on the pine sites. The results also suggested that site preparation affects the soil-water regime and thus prerequisites for forest growth over two decades after site preparation. High variation in the survival and mean height of planted pine was found. The study suggested that on spruce sites, pine survival is the lowest on sites that dry out slowly after rainfall events, and that height growth is the fastest on soils that reach favourable aeration conditions for root growth soon after saturation, and/or where the average air-filled porosity near field capacity is large enough for good root growth. Survival, but not mean height can be enhanced by employing intensive site preparation methods on spruce sites. On coarser-textured pine sites, site preparation methods don t affect survival, but methods affecting soil fertility, such as prescribed burning and ploughing, seem to enhance the height growth of planted Scots pines over several decades. The use of soil water content in situ as the sole criterion for sites suitable for pine reforestation was tested and found to be a relatively uncertain parameter. The thesis identified new potential soil variables, which should be tested using other data in the future.