Application of two methods to determine hydraulic conductivity with disc permeameters on sloping land

The purpose of the present study was to compare two methods for estimating the hydraulic conductivity near saturation with disc permeameters, and to identify their merits when studying runoff on sloping land. The soil's hydraulic conductivity (K) was measured with disc permeameters at a sequence of nominal pressures (ψ) in three blocks with average slope gradients of 11.0% (two occasions), 21.5% (two occasions) and 29.3% (one occasion), respectively, within a sloping area, 40 m wide and 100 m long. Two different methods were used. In the first, the split‐location method, the permeameter was moved to an adjacent spot after measurement at each applied ψ. The estimate of K(ψ) was based on measured sorptivity, steady‐state volumetric flow, initial volumetric water content and the volumetric water content at the applied ψ. In the second method, the one‐location method, the permeameter was not moved during the measurements at each sequence of applied ψ and the estimate of K(ψ) was based only on steady‐state volumetric flow and piece‐wise application of the exponential relation between K and ψ. The latter method generally gave smaller estimates of K than the former on the gentle slopes. These differences were smaller or negligible on steeper slopes. The slope gradient and the conditions in the uppermost soil layers had a definite influence on the values of K obtained. The one‐location method is recommended in studies in which disturbance of the soil surface must be kept to a minimum, as is the case in experimental plots.     

Response of soil structure and hydraulic conductivity to reduced tillage and animal manure in a temperate loamy soil

We studied the combined effects of reduced tillage and animal manure on soil structure and hydraulic conductivity (K) in the 2–10 and 12–20 cm layers in a loamy soil. The study was performed at the end of a 7‐yr field trial and included three tillage treatments (mouldboard ploughing until 25 cm depth: MP, shallow tillage until 12 cm depth: ST, no‐till: NT) and two fertilizer application treatments (mineral or poultry manure). Soil structure was assessed through bulk density (ρ_b), micromorphological and macropore‐space characteristics. K was measured in situ at ?0.6, ?0.2 and ?0.05 kPa. Untilled layers had a vermicular microstructure resulting from earthworm activity, whereas tilled layers displayed a mixture of crumb and channel microstructures. Untilled layers had the highest ρ_b and twice as much lower total macroporosity area (pores > 240 μm in equivalent diameter) than tilled layers, reflected by the smallest area of macropores 310–2000 μm in diameter and the smallest area of large complex macropores. K under untilled layers was 12–62% lower than that under tilled layers, but differences were statistically significant only at ?0.05 kPa in the 2–10 cm. No significant interaction between tillage and nutrient application treatments was detected for all properties. Compared with mineral fertilizer, poultry manure resulted in a similar ρ_b but 20% greater total macroporosity area and 30% higher K at ?0.2 kPa. Overall, the sensitivity of soil structure and K to poultry manure were relatively small compared with tillage. We suggest that cultivation practices other than animal manure application are needed to improve physical properties under reduced tillage.     

Determination of unsaturated soil hydraulic properties at different applied tensions and water qualities

Irrigation with low-quality water may change soil hydraulic properties due to excessive electrical conductivity (EC_w) and sodium adsorption ratio (SAR_w). Field experiments were conducted to determine the effects of water quality (EC_w of 0.5–20 dS m^?1 and SAR_w of 0.5–40 mol^0.5 l^?0.5) on the hydraulic properties of a sandy clay loam soil (containing ～421 g gravel kg^?1 soil) at applied tensions of 0–0.2 m. The mean unsaturated hydraulic conductivity [K(ψ)], sorptive number (α) and sorptivity coefficient (S) varied with change in EC_w and SAR_w as quadratic or power equations, whereas macroscopic capillary length, λ, varied as quadratic or logarithmic equations. The maximum value of K(ψ) was obtained with a EC_w/SAR_w of 10 dS m^?1/20 mol^0.5 l^?0.5 at tensions of 0.2 and 0.15 m, and with 10 dS m^?1/10 mol^0.5 l^?0.5 at other tensions. Changes in K(ψ) due to the application of EC_w and SAR_w decreased as applied tension increased. Analysis indicated that 13.7 and 86.3% of water flow corresponded to soil pore diameters <1.5 and >1.5 μm, respectively, confirming that macropores are dominant in the studied soil. The findings indicated that use of saline waters with an EC of <10 dS m^?1 can improve soil hydraulic properties in such soils. Irrigation waters with SAR_w < 20 mol^0.5 l^?0.5 may not adversely affect hydraulic attributes at early time; although higher SAR_w may negatively affect them.     

Hydraulic conductivity and porosity under conventional and no-tillage and the effect of three species of cover crop in northern France

We studied soil hydraulic conductivity (K) and porosity in five combinations of soil tillage and cover crop management systems. Treatments were winter wheat (Triticum aestivum L.) grown on a conventionally tilled soil (CT), on a no‐till soil (NT), and on an NT with three different cover crops: red fescue (Festuca rubra L.; Fr), bird's‐foot‐trefoil (Lotus corniculatus L.; Lc) and alfalfa (Medicago sativa L.; Ms). Measurements were made on a loamy soil in Grignon, France, in November 2004, May 2005 and October 2005. K and mean size of hydraulically active pores were measured in situ at three water potentials (?0.6, ?0.2 and ?0.05 kPa) at the soil surface and at 10 cm depth. In November 2004 and May 2005, pore space was described using 2D image analysis of pores on undisturbed soil samples in the 0–10 cm layer and in the 10–20 cm layer. The major differences were caused by soil tillage that created two heterogeneous soil layers and increased K in the 0–10 cm layer relative to NT. The effects of cover crop on K and porosity were not affected by the root type: there were no major differences between the grass cover crop (fibrous‐root type) and the leguminous ones (tap‐root type). However, we recorded larger functional pores and more tubules in the no‐till treatments with a cover crop, compared with the no‐till treatment without cover crop; this was probably the result of root activity. Although these changes generally did not result in larger values of K, they participated in the maintenance of soil structure and K over time.     

Prediction of hydraulic conductivity and sorptivity in soils at steady-state infiltration

The purpose of this study was (1) to find a matching factor (u) between infiltration rate and hydraulic conductivity during steady-state infiltration, and (2) to propose equations based on infiltration and soil moisture-retention functions for prediction of the hydraulic conductivity K(θ) within the rapidly (non-capillary) drainable pores (RDP) and capillary-matrix pores of soils. The K(θ) of capillary pores was divided into K(θ)_SDP, K(θ)_WHP and K(θ)_FCP within slowly drainable pores (SDP), water-holding pores (WHP) and fine capillary pores (FCP), respectively. Five soil profiles of calcareous sandy loam, alluvial saline and non-saline clay, located at the Nile Delta, were used to apply the proposed equations. The highest and the lowest values of K(θ)_RDP were observed in calcareous and saline clay soil profiles, respectively. Values of K(θ)_RDP remained higher than those for capillary pores in the studied soils. The predicted values of K(θ) in capillary and non-capillary pores classes were in the expected range for unsaturated hydraulic conductivity. Water sorptivity (S) was determined at initial unsaturated soil water conditions and calculated at steady-state infiltration (S _w) using a derived equation. There was a decrease in S with an increase in soil water content; i.e. at steady-state infiltration, S decreased by 35–40% in calcareous soils and by 45–60% in alluvial clay soils. The parameter values of u and S _w tended to be uniform in calcareous soils, but nonuniform in saline and non-saline clay soils.     

Effect of long‐term swine‐manure application on soil hydraulic properties and heavy metal behaviour

This study evaluated the effect of 13 years of swine‐manure application on the changes in soil hydraulic properties, and as associated physicochemical properties, with a focus on heavy metal mobility. Various soil hydraulic properties were measured, including soil water retention (SWR), saturated field hydraulic conductivity (K_fs) and unsaturated field hydraulic conductivity (K_funsat) using a disc infiltrometer. Heavy metal mobility was evaluated with a sequential extraction procedure. At 0–30 cm soil depth in the heavily manured plot (SMhigh plot), SWR at 0 to ?100 kPa was significantly larger than in plots amended with a standard amount of manure (SMstd plot) or with chemical fertilizer (CF plot). K_fs and K_funsat values in both manure‐amended plots were less than in the CF plot under dry soil conditions but greater than those of the CF plot under wet soil conditions. Furthermore, K_fs and K_funsat did not necessarily increase with manure application rates. On the other hand, high‐mobility metal fractions, such as the exchangeable fraction of Zn, and the CH₃CO₂Na‐extractable fraction of Zn and Mn, and the metal–organic complex fractions of Zn, Cu and Mn, increased with the greater manure application rate. In addition, low‐mobility metal fractions, the organically bound fractions of Zn, Cu and Mn in the high SM plot and the easily reducible metal oxide fraction of Mn in both manure‐amended plots were probably affected and released into high‐mobility fractions. This indicated that manure application changed the soil redox conditions by improving the soil structure, depending on the water content of soil pores. Despite the reduction of K_fs and K_funsat by heavy manure application, the transport of high‐mobility metal fractions with either surface water flow or infiltration water flow could be controlled by soil water content at the beginning of a rain or irrigation event.     

Artificial neural networks for predicting unsaturated soil hydraulic characteristics at different applied tensions

Direct measurement of unsaturated hydraulic parameters that are crucial inputs in any modeling of water flow or solute transport through the vadose zone is costly and time-consuming. Therefore, indirect methods like artificial neural networks (ANNs) can be used to estimate these parameters. Different ANNs conditions [two training algorithms (Trainlm and Traingd), two transfer functions (Tansig and Logsig), and different combinations of the input variables sand, silt, clay, bulk density (BD), soil organic matter (SOM) and initial (?_i) and saturated (?_s) volumetric water content] were used to predict the unsaturated hydraulic conductivity, K _ψ [at six applied tensions (ψ) of 0–0.2 m] and sorptive number, α (at five ψ values of 0.03–0.2 m) that correspond to 138 soil samples from two neighboring study areas located in the Agricultural College, Shiraz University, Islamic Republic of Iran. A four-layer ANNs with three and four nodes in the hidden layers, performed the best in predicting K _ψ and α. Traingd produced the best predictions over all input variables when Tansig and Logsig transfer functions were used for K _ψ and α at different values of ψ, respectively. ‘Silt + clay + sand + BD + SOM’ was the most basic influential input variable for predictions of K _ψ and α at almost all values of ψ. The accuracy of ANNs-predicted K _ψ decreased with decreasing ψ from 0.2 to 0.1 m, followed by an increase to a ψ value of 0 m; a specific relationship was not observed for α. Therefore, ANNs can be used to predict K _ψ with greater confidence at moderate values of ψ than at lower or higher values. The normalized root mean square error, NRMSE, confirmed that ANNs predictions for K _ψ were more accurate than predictions for α. Because reliable predictions were obtained for K _ψ, and to a lesser extent for α, it is recommended that such intelligence models are used to predict these vital soil hydraulic attributes.     

Estimation of saturated hydraulic conductivity from double‐ring infiltrometer measurements

This research aims to determine soil vertical saturated hydraulic conductivity (K_s) in situ from the measured steady infiltration rate (I), initial soil properties and double‐ring infiltrometer (DRI) test data. Characterizing the effects of these variables on the measured steady infiltration rate will enable more accurate prediction of K_s. We measured the effects of the ring diameter, head of ponding, ring depth, initial effective saturation and soil macroscopic capillary length on measured steady infiltration rates. We did this by simulating 864 DRI tests with the finite element program HYDRUS‐2D and by conducting 39 full‐scale in situ DRI tests, 30 Mini‐Disk infiltrometer experiments and four Guelph Permeameter tests. The M5′ model trees and genetic programming (GP) methods were applied to the data to establish formulae to predict the K_s of sandy to sandy‐clay soils. The nine field DRI tests were used to verify the computer models. We determined the accuracy of the methods with 30% of the simulated DRI data to compare I/K_S values of the finite element models with estimates from the suggested formulae. We also used the suggested formulae to predict the K_s values of 30 field DRI experiments and compared them with values measured by Guelph Permeameter tests. Compared with the GP method, the M5′ model was better at predicting K_S, with a correlation coefficient of 0.862 and root mean square error (RMSE) of 0.282 cm s^?1. In addition, the latter method estimated K_svalues of the field experiments more accurately, with an RMSE of 0.00346 cm s^?1.     

Comparison of soil water chemistry and sample size requirements for pan vs tension lysimeters

A study was conducted to compare soil leachate chemistry and determine sample size requirements for tension vs pan (zero-tension) lysimeters. Analyses were performed on an annual and seasonal basis for one year of data collected at Pea Vine Hill, a forested site in southwestern Pennsylvania. On an annual basis, SO₄ ^?2, Ca⁺², Mg⁺², Mn⁺², K⁺ and specific conductance were significantly higher in tension lysimeter samples but no chemical species were significantly higher in pan lysimeters. Seasonal comparisons indicated chemical differences between lysimeter types were variable with more significant deviations present during wet periods. Nearly all significant seasonal differences were comprised of higher concentrations in tension compared to pan lysimeters. Disparities in leachate chemistry between lysimeter types were ascribed to different sources of water collected by the instruments especially during wet periods. Sample size requirements were calculated for two biweekly periods for each lysimeter type at three confidence levels. Based upon calculated sample demands, pan lysimeter soil leachate chemistry could be characterized with fewer samples than tension lysimeters. Less than .30 samples were generally necessary for pan B-horizon lysimeters at the 70% confidence level. Sample requirements were usually unreasonable at higher confidence levels.     

Tillage effects on shear strength and bulk density of soil aggregates

An experiment was conducted to determine the effect of four tillage systems (moldboard plow, chisel plow, Paraplow and no-till) on soil aggregate shear strength and bulk density. Two soils, a Canisteo clay loam (fine-loamy, mixed (calcareous), mesic, Typic Haplaquoll) and a Haig silt loam (fine, montmorillonitic, mesic, Typic Argiaquoll) were used in this study. Soil samples were collected from the 0.075–0.15-m-depth increment in 1983 and the 0.075–0.15- and 0.225–0.30-m-depth increments in 1985. Shear strength of soil aggregates 0.02–0.03 m in diameter was measured by a fall-cone penetrometer and bulk density of the same aggregates was measured by gamma-ray attenuation. Aggregates were tested at soil water matric potentials (ψ_m) of −0.2, −1.1 and −4.0 kPa in 1983 and at ψ_m of −0.2, −1.1, −4.0 and −7.9 kPa in 1985. Tillage for the 1983 growing season was conducted under very wet conditions, whereas tillage for the 1985 growing season was conducted under much drier conditions. Samples collected in 1983 showed little tillage effect on shear strength or bulk density. In 1985, tillage had an effect on shear strength and bulk density for the Haig soil, but not for the Canisteo soil. Much of the tillage effect on soil aggregate shear strength could be explained by tillage-induced changes in the aggregate bulk density. As bulk density decreased, soil aggregate shear strength decreased.Sampling depth had no effect on soil aggregate shear strength or bulk density. Matric potential had an effect on soil aggregate shear strength and bulk density. As matric potential decreased, both shear strength and bulk density increased.     

Surface soil hydraulic properties in four soil series under different land uses and their temporal changes

The concepts of “genoform” and “phenoform” distinguish the genetically-defined soil series and the variation of soil properties resulted from different land uses and management practices. With the repeated field measurements over time, we attempted to understand the difference of soil hydraulic properties among different land uses for a given soil series, and their temporal dynamics. Four soil series (Glenelg, Hagerstown, Joanna, and Morrison) in Pennsylvania with contrasting textures, structures, and parent materials were investigated. Within each soil series, four common land uses (woodland, cropland, pasture, and urban) were examined. At each site of soil series–land use combination, field-saturated and near-saturated hydraulic conductivities, K(ψ), were measured at the soil surface using standard tension infiltrometers at water supply potentials (ψ) of − 0.12, − .06, − 0.03, − 0.02, − 0.01, and 0 m. Surface infiltration measurements were repeated at each site in May and October from 2004 to 2006. The analysis of variance indicated that the measurement time (May or October) had the greatest impact on all measured hydraulic conductivities (p < 0.001), followed by the land use (p < 0.05 for K_ψ = 0 and K_{ψ = − 0.06}) and soil series (p < 0.06 for K_{ψ = − 0.01} to K_{ψ = − 0.03}). The interactions between the time and land use and between the soil series and land use were statistically significant for K_ψ = 0 and K_{ψ = − 0.01.} When separated by the measurement time, land use showed greater impacts in October than in May, while soil series had greater impacts in May than in October. Among the four land uses, woodland showed less obvious temporal change compared to the other three land uses because of less human-induced impacts and more consistent ground cover. Other three land uses generally showed a higher hydraulic conductivity in May than in October due to the drier initial soil moisture condition and related management practices in the spring that gave rise to more significant macropore flow. The results suggested that the initial soil moisture is an important variable that drives the temporal variation of the surface soil hydraulic properties.     

The role of organic acids in the soil solution chemistry of a podzolized soil

The soil solution chemistry of a podzolized soil in the north of Sweden was monitored for four years using percolation lysimeters. Weak organic acids were a major constituent of the soil solution and are important because of their ability to form complexes with aluminium. Dissolved organics leached from the mor layer enhance the weathering rate in the eluvial horizon by forming complexes with aluminium, especially during the autumn when the leaching of dissolved organics was greatest. The weak organic acids were titrated and their pK_a values were evaluated. Aluminium was speciated with an ion-exchange method and by applying equilibrium calculations. Formation constants for the organic aluminium complexes were calculated to be log K_Along=5.42±0.32 m ^?1 (n=13) in spring and summer and log K_Alorg=4.87±0.14 m ^?1 (n=6) in autumn. Equilibria of Al³⁺ with solid phases were also examined using solubility constants. Percolation lysimeters below undisturbed and cut-off mor layers were compared.     

Quantity‐intensity parameters of potassium in relation to uptake by guinea‐corn in representative soils of the ecological zones of Nigeria

Abstract

Quantity‐intensity (Q‐I) relation studies were often used to supplement information obtained from conventional soil tests for the estimation of potassium (K) needs of crops. With a view to ascertaining the reliability of the Q‐I relation parameters for comprehensive characterization of K dynamics in typical Nigerian soils, K values derived from Q‐I isotherms were related to neutral normal ammonium acetate (1 N NH₄OAc, pH 7.0) (exchangeable) K, other soil K forms [non‐exchangeable (K_ne), exchange (K_e), mineral K (K _m ), and solution K (K_s)] and the K uptake by Guinea‐corn (Sorghum bicolor, var. LS 187) subjected to weekly cuts in Neubauer cultivation vessels. Most of the soil K (about 98%) was in the form of soil minerals while less than 1% was plant available whereas about 1% was trapped within the interlattice layers of the clay minerals (as fixed K or K_ne). Mineral K (K_m) content was closely related to total K (K_t), but not to the other forms, K_ne, K_e, and K_s. A close relationship was noted between the two components of labile K (K_e and K_s). Except for % K saturation, the relationships between the K measurements with plant response were poor. The results of these investigations clearly demonstrate that the Q‐I relation could not adequately characterize the K dynamics of these tropical soils.     

Ability of Forage Grasses Exposed to Atrazine and Isoxaflutole to Reduce Nutrient Levels in Soils and Shallow Groundwater

Abstract

Successful implementation of vegetative buffers requires inclusion of plant species that facilitate rapid dissipation of deposited contaminants before they have a chance to be transported in surface runoff or to shallow groundwater. Thirty‐six field lysimeters with six different ground covers [bare ground, orchardgrass (Dactylis glomerata L.), tall fescue (Festuca arundinacea Schreb.), smooth bromegrass (Bromus inermis Leyss.), timothy (Phleum pratense L.), and switchgrass (Panicum virgatum L.)] were established to evaluate the ability of grasses to reduce nutrient levels in soils and shallow groundwater. Nitrate (NO₃ ^?) and orthophosphate (PO₄ ^3?) were uniformly applied to each lysimeter. In addition, half of the lysimeters received an application of atrazine, and the other half received isoxaflutole (Balance?) at levels indicative of surface runoff from cropland. The leachate from each lysimeter was collected after major rainfall events during a 25‐day period, and soil was collected from each lysimeter at the end of the 25‐day period. Water samples were analyzed for NO₃‐N and PO₄‐P, and soil samples were analyzed for NO₃‐N. Grass treatments reduced NO₃‐N levels in leachate by 74.5 to 99.7% compared to the bare ground control, but timothy was significantly less effective at reducing NO₃‐N leaching than the other grasses. Grass treatments reduced residual soil NO₃‐N levels by 40.9 to 91.2% compared to the control, with tall fescue, smooth bromegrass, and switchgrass having the lowest residual levels. Switchgrass decreased PO₄‐P leaching to the greatest extent, reducing it by 60.0 to 74.2% compared to the control. The ability of the forage grasses to reduce nutrient levels in soil or shallow groundwater were not significant between herbicide treatments. Quantification of microbial NO₃ ^? dissipation rates in soil suggested that denitrification was greatest in switchgrass, smooth bromegrass, and tall fescue treatments. The overall performance of these three grasses indicated that they are the most suitable for use in vegetative buffers because of their superior ability to dissipate soil NO₃ ^? and reduce nutrient transport to shallow groundwater.     

Water and Bromide Recovery in Wick and Pan Lysimeters under Conventional and Zero Tillage

Abstract

Quantifying in situ solute transport through soils and the landscape has been widely acknowledged as important and yet challenging. The objective of this study was to evaluate water and bromide movement in no‐tilled (NT) and conventionally tilled (CT) corn using two different types of in situ lysimeters—pan and capillary wick—for single rainfall events. Four zero‐tension pan and four capillary‐wick lysimeters were installed 1.2 m deep on opposite sides of four soil pits. Two were under NT corn, and two were under CT corn. Bromide (Br) was either surface applied or applied with an initial 25 to 27 mm of irrigation (33 to 34 g Br m²). A total of 120 to 147 mm of irrigation was applied continuously at 8.8 mm h^?1. Leachate was collected on 15 min intervals for 24 h and on greater intervals for up to 350 h. Lysimeter discharge and Br concentration were determined for each interval. After drainage began and until rainfall was discontinued, the water drainage rate was, on average, greater in NT (7.2 mm h^?1) than in CT (5.6 mm h^?1) based on results from the pan lysimeters. By contrast, the water drainage rate for the wick lysimeters was, on average, greater in CT (7.3 mm h^?1) than in NT (3.0 mm h^?1). The wick lysimeter appears to have behaved as a sink under the CT conditions, likely representing water flow in smaller channels. Under NT conditions, greater discharge observed with the pan lysimeter implicates the response from larger channels as the conduit for water flow. Flow‐weighted mean Br concentration was less when Br was applied on the soil surface (17.9 mg L^?1) than when Br was applied with the irrigation water (50.6 mg L^?1). Implications from preferential flow studies are often determined based on a single method of evaluation for solute transport, which are likely subject to the limitation of the method used. This study illustrates that contrary to the conventional understanding about preferential flow in NT, water flow and Br transport to the 1.2-m depth was as great as or greater with CT than with NT based on the results from the wick lysimeters for single rainfall events.     

Nutrient leaching from clay soil monoliths with variable past manure inputs

Overall, arable soils in Sweden are currently generally close to phosphorus (P) balance, but excessive P accumulation has occurred on animal fur farms, i.e., those rearing mink (Mustela vison) and foxes (Alopex lagopus and Vulpes vulpes). Manure P from these farms has sometimes regarded as sparingly soluble. Laboratory lysimeter topsoil trials with simulated rain demonstrated that potential leaching of P in dissolved reactive form (DRP) can be very high, even for heavy clay (50%–65%) soils. The Swedish/Norwegian soil test P‐AL (soil P extracted with acid NH₄ lactate, AL) proved useful as a potential indicator of DRP leaching risk (regression coefficient [R²] = 0.89) from fur farms. The upper 5‐cm soil layer, with 190% higher (median) soil P status than the 5–20 cm layer, was the major source of potential DRP leaching through soil columns at the site, despite having been under grass or green fallow for the past 8 y. In percolate from topsoil lysimeters, DRP concentration increased by 0.29 mg L^–1 after the long‐term manure application but only by 0.14 mg L^–1 after the single slurry application when compared to no addition of slurry. Therefore, the build‐up to a high soil P status due to the long‐term application of mink manure was more important than a single application of pig slurry at a rate corresponding to 22 kg P ha^–1 with respect to soil leachate DRP losses in this lysimeter study. The study stresses the importance of precision farming, in which the amount of slurry‐P applied is based on testing the already existing soil P content.     

Spatial variability of physico-chemical properties and hydraulic characteristics of a gravelly calcareous soil

Due to the existence of gravelly soils and the lack of sufficient research on such soils, this study was carried out on a gravelly calcareous soil. Selected physico–chemical and hydraulic soil attributes were determined at 69 points on a nested-sampling design. Hydraulic characteristics including unsaturated hydraulic conductivity (K _ψ) and sorptivity coefficient (S) at six applied tensions of 0 to 0.2 m, and sorptive number (α) and macroscopic capillary length (λ) at five applied tensions of 0.03 to 0.2 m were determined using a tension disc infiltrometer. Hydraulic and chemical soil attributes were the highest and the lowest variants, respectively. The maximum and minimum values for the coefficient of variation (CV) in all the measured physico-chemical and hydraulic soil attributes were obtained for α_0.2 and soil pH, respectively. Minimum, maximum, mean and variance values of K _ψ, S and α decreased as applied tension increased. Although the pattern was reversed for λ. The majority of soil attributes showed the spatial structure with dominant spherical and exponential models for physico-chemical and hydraulic attributes, respectively. Range values of semi-variograms were obtained between 4.6 m (for α_0.03) and 211 m (for clay, gravel content and soluble Mg). In general, range values were 99.60, 82.05 and 40.2 m for physical, chemical and hydraulic soil attributes, respectively, indicating that the physical soil attributes influenced neighboring values over greater distances than the other soil attributes. This enables soil scientists to use measured soil physical data over greater distances to estimate attributes in the unsampled locations.     

Spatio-temporal variation of anisotropy of saturated hydraulic conductivity in a tilled sandy loam soil

Knowledge on anisotropy of saturated hydraulic conductivity can improve the understanding of transport phenomena in soil. We hypothesized that saturated hydraulic conductivity (K_s) in the upper part of the root zone of an agricultural sandy loam soil was anisotropic at different soil depths and times after tillage. K_s was measured on undisturbed 100 cm³ core samples taken in the horizontal and vertical directions in up to four soil layers (Surf: surface layer (0–5 cm); Top: topsoil (10–15 cm); Trans: transition layer between topsoil and subsoil; Sub: subsoil (40–60 cm)) 1, 8 and 32 months, respectively, after mouldboard ploughing and drilling. The ratio between estimated geometric mean values for K_s in the vertical and the horizontal directions (K_ms,v/K_ms,h) was used to test the hypotheses. A total of 669 soil samples were analysed.K_ms,v/K_ms,h varied with time after tillage and between soil layers. One month after ploughing, K_ms,v/K_ms,h was <0.23 (P = 0.975) in the Trans layer with an average value of 0.084, i.e. K_ms,h was 12 times larger than K_ms,v. Anisotropy could not be documented in this layer 8 or 32 months after ploughing, i.e. K_ms,v/K_ms,h was not significantly different from 1.0. For the Surf and Top layers 32 months after ploughing, K_ms,v/K_ms,h was in the intervals 1.4–50 and 3.1–77, respectively, (P = 0.95) with average values of 8.4 and 15, respectively. Thus, K_ms,v was 8.4 respectively 15 times larger than K_ms,h in the two layers. Anisotropy was not found in these layers 1 or 8 months after tillage. Strong anisotropy was found in the Sub layer with K_ms,v/K_ms,h averaging to 14 and 32, respectively, 8 and 32 months after tillage. K_ms,v and K_ms,h generally decreased with time in the Surf, Top and Trans layers, except in the vertical direction in the Top layer between 8 and 32 months after ploughing, and in the Trans layer between 1 and 8 months after ploughing. Overall, the geometric means of K_s varied between 10^−4.0 and 10^−7.1 m s⁻¹.The results may reflect systematic effects of soil settlement and drying/wetting phenomena coupled with biological activity and the existence of stable, vertically oriented biopores in the subsoil. It appears to be necessary to consider anisotropy of K_s and its variation in the analysis and modelling of water flow and chemical transport in agricultural soils, particularly to explain heterogeneous flow phenomena at the plot and field scales.     

Simulation of water flow and solute transport in free-drainage lysimeters and field soils with heterogeneous structures

Lysimeters are valuable for studying the fate and transport of chemicals in soil. Large‐scale field lysimeters are used to assess pesticide behaviour and radionuclide transport, and are assumed to represent natural field conditions better than laboratory columns. Field lysimeters are usually characterized by a free‐draining lower boundary. As a result, the hydraulic gradient is disrupted, and leachate cannot be collected until the bottom of the lysimeter becomes saturated. We compared heterogeneously structured, free‐drainage lysimeters and field soils with respect to water flow and solute transport. Numerical simulations were carried out in a two‐dimensional heterogeneous sandy soil under unsaturated water flow conditions with the CHAIN_2D code. Three different soil structures (isotropic, horizontal, and vertical) were generated, and Miller–Miller similitude was used to scale the hydraulic properties of the soil. The results showed that ponding occurs at the bottom of the lysimeter for the three soil structures and that it occurred faster and was more pronounced with the vertical structure (preferential flow effect). Breakthrough curves of a conservative solute (bromide) showed that solutes are moving faster in the field than in the lysimeters. Fewer differences between lysimeters and field soils were found with the horizontal soil structure than with the isotropic and vertical structures.