Heterogeneous Sandy Soil Research Articles

Numerical investigation of streamtube approach to model flow in heterogeneous unsaturated sandy soils

The streamtube approximation has been commonly used for upscaling of flow and transport in spatially heterogeneous unsaturated soils. This study investigates the prediction ability of streamtube approach to model flow in heterogeneous unsaturated sandy soils under a rainfall event at a plot scale. Towards this aim, a numerical model for solving 3-D Richards’s equation was developed and verified for 1-D and 3-D cases. Turning band method was used to generate the random field of soil hydraulic parameters. The 3-D soil field heterogeneity was represented by considering variation in both van Genuchten parameter (α) and Ks. For streamtube models, two cases of heterogeneity were investigated: i) variability in Ks only, and ii) variability in both Ks and α. Both the parameters were assumed to follow a log-normal distribution. For different boundary conditions at the soil surface and varying soil heterogeneity, the average Darcy's flux and mean field-scale soil moisture obtained using streamtube approach were compared with 3-D numerical solutions. The results identify that for high recharge rate at the soil surface, streamtube model is a good approximation of the 3-D heterogeneous soil field, and the average difference in mean soil moisture obtained from 3-D model and streamtube based 1-D model was below 10% when variability in both Ks and α was considered for all timesteps. The inclusion of variability in α leads to improvement in prediction capability of streamtube model. Further, regression equations were developed that can be utilized to predict the mean field soil moisture for heterogeneous 3-D sandy soil fields given the mean field soil moisture obtained from streamtube approach for both the cases. The regressions equations follow exponential model for case-1 and linear model for case-2 with R2 values 0.80 and 0.99, respectively. The developed equations can be reliably used to predict mean field soil moisture for heterogeneous sandy soil fields.

Changes in functional diversity of the soil microbial community in a heterogeneous sandy soil after long-term fertilization with cattle manure and mineral fertilizer

The effects of cattle without and with biodynamic preparations on functional diversity of the soil microbial community were investigated in comparison with mineral fertilization (+ straw incorporation) at low and high application rates to a sandy soil of the long-term fertilization trial in Darmstadt, Germany. The multi-SIR method was used for investigating the respiratory response to 17 individual substrates in surface arable soils (0–5cm). Multivariate analysis made it possible to significantly separate mineral fertilizer and manure treatments, mainly on the basis of carbohydrates and amino acids. Correlation analysis of the extracted discriminant function showed that the changes in the catabolic profiles of soil microorganisms between the mineral fertilizer and manure treatments were partly caused by differences in soil pH and soil organic C content. The functional diversity of soil microorganisms caused by differences in the community level physiological profile decreased as a result of long-term mineral fertilization with straw incorporation in comparison with farmyard manure application.

Analyzing the efficiency of soil amendments and irrigation for plant production on heterogeneous sandy soils under greenhouse conditions

AbstractVariability in soil properties is a complication for fertilization, irrigation, and amendment application. However, only limited progress has been made in managing soil variability for uniform productivity and increased water‐use efficiency. This study was designed to ameliorate the poor‐productivity areas of the variable sandy soils in Florida citrus groves by using frequent small irrigations and applying organic and inorganic soil amendments. Two greenhouse experiments were set up with sorghum and radish as bioassay crops in a randomized complete block design (RCBD). The factors studied were two soil‐productivity classes (very poor and very good), two water contents (50% and 100% of field capacity), two amendments (phosphatic clay and Fe humate), and two amendment rates (10 and 25 g kg–1 for sorghum and 50 and 100 g kg–1 for radish). Amendments applied at 50 and 100 g kg–1 increased the water‐holding capacity (WHC) of poor soil by 2‐ to 6‐fold, respectively. The lower rates (10 and 25 g kg–1) of amendments were not effective in enhancing sorghum growth. The higher rates (50 and 100 g kg–1) doubled the radish growth as compared to the control. The results indicate that rates greater than 50 g kg–1 of both amendments were effective in improving water retention and increasing productivity. Irrigation treatment of 100% of field capacity (FC) increased the sorghum and radish growth by about 2‐fold as compared with the 50%–water content treatment. The results suggest that the root‐zone water content should be maintained near FC by frequent small irrigations to enhance water availability in excessively drained sandy soils. In addition, application of soil amendments in the root zone can enhance the water retention of these soils. Furthermore, managing variable sandy soils with WHC‐based irrigation can increase water uptake and crop production in the poor areas of the grove.

Pedotransfer functions in soil electrical resistivity estimation

Surface electrical resistivity tomography (ERT) is recognized as a powerful non-invasive soil survey and monitoring method. Relationships between electrical resistivity (ER) and soil water contents that are needed to infer the spatial distribution of soil moisture from the ERT results, are known to reflect soil properties. However, pedotransfer relationships that relate ER, water content, and basic soil properties other than clay content have not been developed. The objective of this work was to evaluate the improvement in the estimates of ER from soil water contents if soil basic properties are included in the pedotransfer equations. Soil — coarse-loamy, siliceous, mesic Typic (or Aquic) Hapludult — was sampled in 12 locations at ten depths across a 30 × 30-m experimental site. The resistivity was measured in disturbed samples in four-electrode cells. Water content, bulk density, texture, organic carbon, and pH were measured in these samples. The coefficient and the exponent in the power law-dependence of ER on volumetric water content were approximated with linear regressions that included any number of measured soil basic properties as inputs. A total of 1024 regressions were compared using the Akaike information criterion, and the best model was found that provided the best approximation with minimum number of regression terms. Six soil textural classes were encountered in the samples. When no pedotransfer relationships were included, the R 2 of the log-log relationships between ER and water content increased from 0.45 at 10 cm to 0.91 at 190 cm. Splitting the dataset into topsoil and subsoil subsets resulted in favorable pedotransfer relationships with R 2 about 0.85 both in topsoil and subsoil. Bulk density, clay content, sand content and pH appeared to be influential soil properties. Overall, accuracy of relationships between soil electrical resistivity and volumetric water content at ten depths in highly heterogeneous sandy soil could be substantially improved if pedotransfer relationships were established. Future work is needed on the applicability of data obtained on repacked samples to undisturbed soils and on the selection of input variables and evaluation methods for pedotransfer electrical resistivity pedotransfer functions.

A spatial approach to soil‐ecological experimentation at landscape scale

AbstractThe upscaling of soil‐ecological processes to larger landscape units represents a special challenge to soil ecology. Results from micro‐ or mesocosms cannot easily be transferred to other scales because effects are often scale‐dependent. In this context, field experiments which take into account the heterogeneity of the landscape may be promising. Therefore, we carried out an experiment based on a transect study in the agrolandscape of NE Germany on heterogeneous sandy soil in which the feeding activity of the soil‐organism community was assessed by means of the bait‐lamina test at each of the 101 transect locations. At every 4th position, prior to the measurement the soil biota were stimulated by a treatment consisting of adding easily available C and water to the soil. Our aim was to test whether this kind of spatial approach enables to separate effects induced by treatments from landscape effects. The results showed a highly variable feeding activity along the transect after 4 weeks. Despite this variability, a basic trend could be identified which was related to a landscape factor, i.e., the relief. On upper‐slope positions, the feeding activity tended to be less in comparison to positions down‐slope. At every 4th position of the transect, the stimulating effect of the substrate and water addition could be clearly detected and quantified with spectral and cross‐spectral analysis. It is concluded that effects of treatments in heterogeneous landscapes may be distinguished from site effects when the signal‐to‐noise ratio is high and soil and treatment effects on the variable of interest are sufficiently different from one another. In a heterogeneous landscape with gradients of site properties, a treatment based on the frequency domain and applied in regular intervals can be distinguished with spectral analysis techniques.

Study of water tension differences in heterogeneous sandy soils using surface ERT

Herbaceous vegetation in the Sahel grows almost exclusively on sandy soils which preferentially retain water through infiltration and storage. The hydrological functioning of these sandy soils during rain cycles is unknown. One way to tackle this issue is to spatialize variations in water content but these are difficult to measure in the vadose zone. We investigated the use of Electrical Resistivity Tomography (ERT) as a technique for spatializing resistivity in a non-destructive manner in order to improve our knowledge of relevant hydrological processes. To achieve this, two approaches were examined. First, we focused on a possible link between water tension (which is much easier to measure in the field by point measurements than water content), and resistivity (spatialized with ERT). Second, because ERT is affected by solution non-uniqueness and reconstruction smoothing, we improved the accuracy of ERT inversion by comparing calculated solutions with in-situ resistivity measurements. We studied a natural microdune during a controlled field experiment with artificial sprinkling which reproduced typical rainfall cycles. We recorded temperature, water tension and resistivity within the microdune and applied surface ERT before and after the 3 rainfall cycles. Soil samples were collected after the experiment to determine soil physical characteristics. An experimental relationship between water tension and water content was also investigated. Our results showed that the raw relationship between calculated ERT resistivity and water tension measurements in sand is highly scattered because of significant spatial variations in porosity. An improved correlation was achieved by using resistivity ratio and water tension differences. The slope of the relationship depends on the soil solution conductivity, as predicted by Archie's law when salted water was used for the rain simulation. We found that determining the variations in electrical resistivity is a sensitive method for spatializing the differences in water tension which are directly linked with infiltration and evaporation/drainage processes in the vadose zone. However, three factors complicate the use of this approach. Firstly, the relation between water tension and water content is generally non-linear and dependent on the water content range. This could limit the use of our site-specific relations for spatializing water content with ERT through tension. Secondly, to achieve the necessary optimization of ERT inversion, we used destructive resistivity measurements in the soil, which renders ERT less attractive. Thirdly, we found that the calculated resistivity is not always accurate because of the smoothing involved in surface ERT data inversion. We conclude that further developments are needed into ERT image reconstruction before water tension (and water content) can be spatialized in heterogeneous sandy soils with the accuracy needed to routinely study their hydrological functioning.

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

SummaryLysimeters 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.