Different responses in bulk density and saturated hydraulic conductivity to soil deformation by logging machinery on a Ferralsol under native forest

Ferralsols have high structural stability, although structural degradation has been observed to result from forest to tillage or pasture conversion. An experimental series of forest skidder passes in an east Amazonian natural forest was performed for testing the effects of mechanical stress during selective logging operations on a clay-rich Ferralsol under both dry and wet soil conditions. Distinct ruts formed up to 25 cm depth only under wet conditions. After nine passes the initially very low surface bulk density of between 0.69 and 0.80 g cm(-3) increased to 1.05 g cm(-3) in the wet soil and 0.92 g cm(-3) in the dry soil. Saturated hydraulic conductivities, initially > 250 mm h(-1), declined to a minimum of around 10 mm h(-1) in the wet soil after the first pass, and in the dry soil more gradually after nine passes. The contrasting response of bulk density and saturated hydraulic conductivity is explained by exposure of subsoil material at the base of the ruts where macrostructure rapidly deteriorated under wet conditions. We attribute the resultant moderately high hydraulic conductivities to the formation of stable microaggregates with fine sand to coarse silt textures. We conclude that the topsoil macrostructure of Ferralsols is subject to similar deterioration to that of Luvisols in temperate zones. The stable microstructure prevents marked compaction and decrease in hydraulic conductivity under wetter and more plastic soil conditions. However, typical tropical storms may regularly exceed the infiltration capacity of the deformed soils. In the deeper ruts water may concentrate and cause surface run-off, even in gently sloping areas. To avoid soil erosion, logging operations in sloping areas should therefore be restricted to dry soil conditions when rut formation is minimal.

Viscous-Flow Approach to In Situ Infiltration and In Vitro Saturated Hydraulic Conductivity Determination

Infiltration is dominantly gravity driven, and a viscous-flow approach was developed. Laminar film flow equilibrates gravity with the viscous force and a constant flow velocity evolves during a period lasting 3/2 times the duration of a constant input rate, qS. Film thickness F and the specific contact area L of the film per unit soil volume are the key parameters. Sprinkler irrigation produced in situ time series of volumetric water contents, θ(z,t), as determined with TDR probes. The wetting front velocity v and the time series of the mobile water content, w(z,t) were deduced from θ(z,t). In vitro steady flow in a core of saturated soil provided volume flux density, q(z,t), and flow velocity, v, as determined from a heat front velocity. The F and L parameters of the in situ and the in vitro experiments were compared. The macropore-flow restriction states that, for a particular permeable medium, the specific contact area L must be independent from qS i.e., dL/dqS = 0. If true, then the relationship of qS ∝ v3/2 could scale a wide range of input rates 0 ≤ qS ≤ saturated hydraulic conductivity, Ksat, into a permeable medium, and kinematic-wave theory would become a versatile tool to deal with non-equilibrium flow. The viscous-flow approach is based on hydromechanical principles similar to Darcy’s law, but currently it is not suited to deduce flow properties from specified individual spatial structures of permeable media.

Evaluation of the Relationship Between Saturated Hydraulic Conductivity and Grain-Size Distribution of Fluvio-Glacial Deposits, Puget Lowland, Washington

The relationship between saturated hydraulic conductivity (Ks) and grain-size distribution was evaluated for 49 sites underlain by either glacially over consolidated or normally consolidated fluvio-glacial deposits in the Puget Lowland. A linear regression comprising pairs of grain-size analyses and pilot infiltration tests predicts Ks with a 1 sigma uncertainty of a factor of about 3.5 with 70% of the population variance accounted for. The correlation coefficient R^2 of about 0.90 shows that there is a strong correlation between the grain-size distribution and Ks. In contrast, a widely applied analysis proposed by Massmann (2003) explains only 20% of the population variance for normally consolidated materials with an R^2 of only 0.15. That analysis entirely fails to explain the population variance for over consolidated materials. The method developed in this study is recommended for determination of Ks for fluvio-glacial deposits of the Puget Lowland.

Effect of cultivation on soil C contents and saturated hydraulic conductivity

The irrigation of pasture with saline, Na-contaminated industrial wastewater typically results in an increase in soil ESP. From current knowledge (derived largely from cultivated agricultural soils), although these sodic soils are likely to remain stable whilst irrigated with effluent (due to the effluent’s large electrolyte concentration), during rainfall periods of low electrolyte concentration these soils would be expected to disperse. However, effluent irrigated pasture soils have been observed to maintain their structure even during intense rainfall events. Three soil types were collected (Sodosol, Vertosol and Dermosol), each with a cultivated/non-cultivated pair. The soils were equilibrated with various SAR solutions and then leached with deionised water to allow the measurement of saturated hydraulic conductivity (Ksat). At low SARs, Ksat tended to be greater in non-cultivated than cultivated soils and is attributable to a loss of structure associated with cultivation. In addition, as SAR increased, the reduction in relative Ksat tended to be significantly greater in cultivated than non-cultivated soils. The relatively rapid saturated hydraulic conductivity in the non-cultivated soils at large SARs is due to a greater aggregate stability due to greater soil C content. For the sustainable disposal of saline effluent, it is therefore necessary to ensure that soils remain undisturbed and preferably under pasture, thus maximising soil structural stability and hydraulic conductivity.

Physical Characteristics and Estimated Hydraulic Conductivity from Grain-Size Distribution for Vashon Till, South-Central Puget Lowland, Washington

Population growth, urban development, and increased commercial and industrial activity in the south-central Puget Lowlands of Washington State has led to an increased demand for groundwater. The Vashon till is a glacially consolidated, low-permeability unit comprising unstratified clay, silt, cobbles and boulders with ubiquitous coarse-grained lenses and is an extensive surficial unit throughout the south-central Puget Lowland. Thus, understanding the physical and hydrological characteristics – specifically, the hydraulic conductivity – of this unit is a necessary component of a groundwater model. This study provides (1) a record of the physical characteristics of Vashon till deposits within the study area; and (2) an estimate of the highest, lowest, and average value of saturated hydraulic conductivity based on the grain-size distribution of Vashon till samples collected from six field sites in the Puyallup River Watershed. Analysis shows that the average moisture content ranges between about 1 and 6%, average dry bulk density is about 2.20 g/cm3, and average porosity is about 17%. Grain-size distributions show that half of the samples analyzed are well graded, while the other half is poorly graded. Grain-size distributions also show an average d10 value of about 0.20 mm, and average ff values ≤ 16%, which are key values in estimating the saturated hydraulic conductivity of over-consolidated glacial deposits. Based on these observed values, the estimates of hydraulic conductivity range from a minimum of 0.02 m/d to a maximum of 1.38 m/d in within the general Vashon till.

Spatial variability of satured soil hydraulic conductivity in the region of Araguaia River - Brazil.

This study evaluates the spatial variability of saturated hydraulic conductivity in the soil in an area of 51,850 ha at the headwaters of the Araguaia River MT/GO. This area is highly vulnerable because it is a location of recharging through natural water infiltration of the Guarani Aquifer System and an area of intense increases in agriculture since its adoption by growers in the last 30 years. Soil samples were collected at 383 points, geographically located by GPS. The samples were collected from depths of 0 - 20 cm and 60 - 80 cm. Exploratory statistics and box-plot were used in the descriptive analysis and semivariogram were constructed to determine the spatial model. The exploratory analysis showed that the mean hydraulic conductivity in the superficial layer was less than at the level of 60-80 cm; however, the greatest variability evaluated with a coefficient of variation also was from this layer. Data tended towards a normal distribution. These results can be explained by the greater soil compaction in the superficial layer. The semivariogram models, adjusted for the two layers, were exponential and demonstrated moderate and strong dependence, with ranges of 5000 and 3000 utm respectively. It was concluded that soil use is influencing the spatial distribution model of the hydraulic conductivity in the region.

Hydraulic conductivity increases in a sodic clay soil in response to gypsum applications impacts of bulk density and cation exchange

Amelioration of sodic soils is commonly achieved by applying gypsum, which increases soil hydraulic conductivity by altering soil chemistry. The magnitude of hydraulic conductivity increases expected in response to gypsum applications depends on soil properties including clay content, clay mineralogy, and bulk density. The soil analyzed in this study was a kaolinite rich sodic clay soil from an irrigated area of the Lower Burdekin coastal floodplain in tropical North Queensland, Australia. The impact of gypsum amelioration was investigated by continuously leaching soil columns with a saturated gypsum solution, until the hydraulic conductivity and leachate chemistry stabilized. Extended leaching enabled the full impacts of electrolyte effects and cation exchange to be determined. For the columns packed to 1.4 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.41 ± 0.06 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 4.3 ± 2.12 mEq/100 g, and hydraulic conductivity increased to 0.15 ± 0.04 cm/d. For the columns packed to 1.3 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.51 ± 0.03 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 0.55 ± 0.36 mEq/100 g, and hydraulic conductivity increased to 0.96 ± 0.53 cm/d. The results of this study highlight that both sodium and magnesium need to be taken into account when determining the suitability of water quality for irrigation of sodic soils and that soil bulk density plays a major role in controlling the extent of reclamation that can be achieved using gypsum applications.

On the effective hydraulic conductivity and macrodispersivity for density-dependent groundwater flow

In this paper, semi-analytical expressions of the effective hydraulic conductivity ( KE) and macrodispersivity ( αE) for 3D steady-state density-dependent groundwater flow are derived using a stationary spectral method. Based on the derived expressions, we present the dependence of KE and αE on the density of fluid under different dispersivity and spatial correlation scale of hydraulic conductivity. The results show that the horizontal KE and αE are not affected by density-induced flow. However, due to gravitational instability of the fluid induced by density contrasts, both vertical KE and αE are found to be reduced slightly when the density factor ( γ ) is less than 0.01, whereas significant decreases occur when γ exceeds 0.01. Of note, the variation of KE and αE is more significant when local dispersivity is small and the correlation scale of hydraulic conductivity is large.

Determination of vertical hydraulic conductivity of aquitards in a multilayered leaky system using water-level signals in adjacent aquifers

This paper presents a methodology for determining the vertical hydraulic conductivity (Kv) of an aquitard, in a multilayered leaky system, based on the harmonic analysis of arbitrary water-level fluctuations in aquifers. As a result, Kv of the aquitard is expressed as a function of the phase-shift of water-level signals measured in the two adjacent aquifers. Based on this expression, we propose a robust method to calculate Kv by employing linear regression analysis of logarithm transformed frequencies and phases. The frequencies, where the Kv are calculated, are identified by coherence analysis. The proposed methods are validated by a synthetic case study and are then applied to the Westbourne and Birkhead aquitards, which form part of a five-layered leaky system in the Eromanga Basin, Australia.

Laboratory measurement of hydraulic conductivity functions of two unsaturated sandy soils during drying and wetting processes

The importance of applying unsaturated soil mechanics to geotechnical engineering design has been well understood. However, the consumption of time and the necessity for a specific laboratory testing apparatus when measuring unsaturated soil properties have limited the application of unsaturated soil mechanics theories in practice. Although methods for predicting unsaturated soil properties have been developed, the verification of these methods for a wide range of soil types is required in order to increase the confidence of practicing engineers in using these methods. In this study, a new permeameter was developed to measure the hydraulic conductivity of unsaturated soils using the steady-state method and directly measured suction (negative pore-water pressure) values. The apparatus is instrumented with two tensiometers for the direct measurement of suction during the tests. The apparatus can be used to obtain the hydraulic conductivity function of sandy soil over a low suction range (0-10 kPa). Firstly, the repeatability of the unsaturated hydraulic conductivity measurement, using the new permeameter, was verified by conducting tests on two identical sandy soil specimens and obtaining similar results. The hydraulic conductivity functions of the two sandy soils were then measured during the drying and wetting processes of the soils. A significant hysteresis was observed when the hydraulic conductivity was plotted against the suction. However, the hysteresis effects were not apparent when the conductivity was plotted against the volumetric water content. Furthermore, the measured unsaturated hydraulic conductivity functions were compared with predictions using three different predictive methods that are widely incorporated into numerical software. The results suggest that these predictive methods are capable of capturing the measured behavior with reasonable agreement.

Analysis and modelling of the hydraulic conductivity in aquitards : application to the Galilee Basin and the Great Artesian Basin, Australia

In large sedimentary basins with layers of different rocks, the groundwater flow between aquifers depends on the hydraulic conductivity (K) of the separating low-permeable rocks, or aquitards. Three methods were developed to evaluate K in aquitards for areas with limited field data: • Coherence and harmonic analysis: estimates the regional-scale K based on water-level fluctuations in adjacent aquifers. • Cokriging and Bayes' rule: infers K from downhole geophysical logs. • Fluvial process model: reproduces the lithology architecture of sediment formations which can be converted to K. These proposed methods enable good estimates of K and better planning of further drillholes.

Hydraulic conductivity fields: Gaussian or not?

Hydraulic conductivity (K) fields are used to parameterize groundwater flow and transport models. Numerical simulations require a detailed representation of the K field, synthesized to interpolate between available data. Several recent studies introduced high-resolution K data (HRK) at the Macro Dispersion Experiment (MADE) site, and used ground-penetrating radar (GPR) to delineate the main structural features of the aquifer. This paper describes a statistical analysis of these data, and the implications for K field modeling in alluvial aquifers. Two striking observations have emerged from this analysis. The first is that a simple fractional difference filter can have a profound effect on data histograms, organizing non-Gaussian ln K data into a coherent distribution. The second is that using GPR facies allows us to reproduce the significantly non-Gaussian shape seen in real HRK data profiles, using a simulated Gaussian ln K field in each facies. This illuminates a current controversy in the literature, between those who favor Gaussian ln K models, and those who observe non-Gaussian ln K fields. Both camps are correct, but at different scales.

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.