Microtopography,water storage and flow patterns in a fine‐textured soil under agricultural use

Agricultural use of soils implies tillage and often compaction and therefore influences processes on soil surface and affects infiltration of water into the subsoil. Although many studies on soil surface processes or flow patterns in soils exist, works relating both are rare in literature. We did two tracer experiments with Brilliant Blue FCF on a tilled and compacted plot and a non‐tilled one to investigate water storage on the soil surface during simulated rainfall and changes of soil microtopography, to analyse the associated flow patterns in the soil and to relate both to tillage and compaction. Our results show that storage was larger on the tilled and compacted plot than on the non‐tilled one. After tillage, transport processes above the plough pan were partly disconnected from those underneath because macropores were disrupted and buried by the tillage operation. However, preferential flow along cracks occurred on both plots and the macropores buried below the tillage pan still functioned as preferential flow paths. Therefore, we conclude that the studied soil is susceptible to deep vertical solute propagation at dry conditions when cracks are open, irrespective of tillage and compaction. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

Fingered flow: The creator of sand columns in dune and beach sands

Sand columns, sand cones, sand mushrooms and other striking sand forms are frequently observed in the Dutch and German beach and dune sands. This paper aims to clarify the mechanism of sand column formation. Recently it has become evident that homogeneous beach and dune sands often become irregularly wetted by infiltrating rainwater. In otherwise dry sandy soils, wet preferential flow paths (‘fingers’) may develop. At two test sites the volumetric soil moisture content varied between 0·2 and 12·0 per cent. The wet fingers represent the premature state of sand columns. When the dry sand in between these fingers is blown away by the wind, the more resistant wet sand of the fingers will remain in its place and appear as sand columns at the surface. As a result of wind and erosive sand drifts, striking sand forms may be formed.     

The effect of alternating different water qualities on accumulation and leaching of solutes in a Mediterranean cracking soil

The relevance of bypass flow on water flow, solute or pesticide transport is becoming increasingly recognized. Recent investigations proved that soil salinization may be influenced by bypass flow, i.e. the rapid transport of water and solutes via macropores and/or shrinkage cracks to subsoil and groundwater. This paper explores the role of bypass flow in the process of accumulation and leaching of solutes, as well as of sodium, in a Mediterranean cracking soil irrigated with saline/sodic waters. The results of bypass flow experiments performed on undisturbed soil cores showed that leaching of solutes occurred in concomitance with bypass fluxes when a low salinity solution was alternated with a high salinity solution. Exchange of solutes between the incoming solution and the soil matrix occurred during the bypass flow events at the contact surfaces (cracks walls) between the solution and the soil matrix and where cracks terminated in the soil samples. Concomitant exchanges of sodium were indicated by measurements performed in the effluent solution during the bypass flow measurements. The amount of Sodium released from the soil during the bypass flow events, as well as that of the soluble salts leached from the soil, were found to depend on the degree of soil cracking. These results indicated that:

• 1 in management of irrigation in cracking soils, under the occurrence of bypass fluxes, alternating a low salinity/sodicity water with a high salinity/sodicity solution can be effective for preventing salinization and sodification:
• 2 greater efficiency of removal of sodium/soluble salts can be obtained if application of the leaching solution is performed when the soil is at a considerable degree of cracking.

Copyright © 2002 John Wiley & Sons, Ltd.     

Characterizing heterogeneous soil water flow and solute transport using information measures

Heterogeneous water flow and solute transport in soils are an important phenomenon and difficult to be characterized. The objectives of this study were to investigate the heterogeneity of solute transport related to heterogeneous soil water flow using dye infiltration experiments, and to characterize heterogeneous water flow and solute transport in soils using the information theory. Field experiments of dye infiltration were performed in four plots. Various information measures were applied to characterize information content and complexity of water flow and solute transport in soils. Information contents and complexities of the maximum and apparent infiltration depths, and the mean and standard deviation of concentrations in the vertical direction of the plots were calculated. More heterogeneous processes of soil water flow and transport result in higher information/complexity values. The probability distributions of mean concentration were similar to those of the corresponding apparent infiltration depths for the plots, indicating that heterogeneity of dye concentrations was closely related to that of soil water flow. However, the range of information entropy and complexity of the water flow sequences was much narrower than that of the sequences of the concentrations. The results suggested that the transport processes were more heterogeneous than the water flow processes. Compared with the probability distributions of flow parameters, the information measures appeared to be a more versatile tool to describe flow and transport heterogeneities in soils.     

Solute transport and water content measurements in clay soils using time domain reflectometry

Abstract

Clayey and saline soils have been shown to be problematic for time domain reflectometry (TDR) measurements. This study presents some of these problems and discusses solutions to them. Thirteen solute transport experiments were carried out in three undisturbed soil columns of swelling clay soil from Tunisia, labelled S1, S2, and S3 respectively. The columns were collected at three different physiographical regions within a catchment. Water fluxes ranged from 1.2 to 7.2 cm day^?1. The large solute transport heterogeneity and large tailing indicated that preferential flow was most pronounced in S1. The preferential flow took place in voids between structural elements and in wormholes. In S3, preferential flow was also evident, but not to the same extent as in S1. In S2, the solute transport was more uniform with little preferential flow. The heterogeneity of the solute transport increased with the water flux in S1 and to a smaller extent in S3, whereas it remained constant in S2. In a previous dye experiment in the field, preferential flow in cracks was observed at those sites where S1 and S3 were collected. In the column experiments, preferential flow in these cracks was less due to the higher initial water content compared to the dye experiments, indicating that the desiccation cracks were closed by the swelling clay.     

Differences in preferential flow with antecedent moisture conditions and soil texture: Implications for subsurface P transport

Preferential flowpaths transport phosphorus (P) to agricultural tile drains. However, if and to what extent this may vary with soil texture, moisture conditions, and P placement is poorly understood. This study investigated (a) interactions between soil texture, antecedent moisture conditions, and the relative contributions of matrix and preferential flow and (b) associated P distributions through the soil profile when fertilizers were applied to the surface or subsurface. Brilliant blue dye was used to stain subsurface flowpaths in clay and silt loam plots during simulated rainfall events under wet and dry conditions. Fertilizer P was applied to the surface or via subsurface placement to plots of different soil texture and moisture condition. Photographs of dye stains were analysed to classify the flow patterns as matrix dominated or macropore dominated, and soils within plots were analysed for their water‐extractable P (WEP) content. Preferential flow occurred under all soil texture and moisture conditions. Dye penetrated deeper into clay soils via macropores and had lower interaction with the soil matrix, compared with silt loam soil. Moisture conditions influenced preferential flowpaths in clay, with dry clay having deeper infiltration (92 ± 7.6 cm) and less dye–matrix interaction than wet clay (77 ± 4.7 cm). Depth of staining did not differ between wet (56 ± 7.2 cm) and dry (50 ± 6.6 cm) silt loam, nor did dominant flowpaths. WEP distribution in the top 10 cm of the soil profile differed with fertilizer placement, but no differences in soil WEP were observed at depth. These results demonstrate that large rainfall events following drought conditions in clay soil may be prone to rapid P transport to tile drains due to increased preferential flow, whereas flow in silt loams is less affected by antecedent moisture. Subsurface placement of fertilizer may minimize the risk of subsurface P transport, particularily in clay.     

Modelling evaporation processes in soils from the Huasco salt flat basin,Chile

The need to understand and simulate hydrological phenomena and their interactions, and the impact of anthropogenic and climate changes on natural environments have promoted the study of evaporation from bare soils in arid climates. In closed Altiplano basins, such as those encountered in arid and hyper arid basins in northern Chile, evaporation from shallow groundwater is the main source of aquifer depletion, and thus, its study is crucial for water resources management. The objective of this work is to understand the mechanisms of evaporation in saline soils with shallow water tables, in order to better quantify evaporation fluxes and improve our understanding of the water balance in these regions. To achieve this objective, a model that couples fluid flow with heat transfer was developed and calibrated using column experiments with saline soils from the Huasco salt flat basin, Chile. The model enables determination of both liquid and water vapour fluxes, as well as the location of the evaporation front. Experimental results showed that salt transport inside the soil profile modified the water retention curve, highlighting the importance of including salt transport when modelling the evaporation processes in these soils. Indeed, model simulations only agreed with the experimental data when the effect of salt transport on water retention curves was taken into account. Model results also showed that the evaporation front is closer to the soil surface as the water table depth reduces. Therefore, the model allows determining the groundwater level depth that results in disconnection of liquid fluxes in the vadose zone. A sensitivity analysis allowed understanding the effect of water‐flux enhancements mechanisms on soil evaporation. The results presented in this study are important as they allow quantifying the evaporation that occurs in bare soils from Altiplano basins, which is typically the main water discharge in these closed basins. Copyright © 2016 John Wiley & Sons, Ltd.     

Comparative performance of soil water balance models in computing semi-arid aquifer recharge

Abstract

Estimating groundwater recharge is essential to ensure the sustainable use of groundwater resources, particularly in arid and semi-arid regions. Soil water balances have been frequently advocated as valuable tools to estimate groundwater recharge. This article compares the performance of three soil water balance models (Hydrobal, Visual Balan v2.0 and Thornthwaite) in the Ventós-Castellar aquifer, Spain. The models were used to simulate wet and dry years. Recharge estimates were transformed into water table fluctuations by means of a lumped groundwater model. These, in turn, were calibrated against piezometric data. Overall, the Hydrobal model shows the best fit between observed and calculated levels (r² = 0.84), highlighting the role of soil moisture and vegetation in recharge processes.

Editor D. Koutsoyiannis; Associate editor X. Chen

Citation Touhami, I., et al., 2014. Comparative performance of soil water balance models in computing semi-arid aquifer recharge. Hydrological Sciences Journal, 59 (1), 193–203.     

The distribution of deuterium and O in dry soils : 1. Theory

A mathematical model is developed describing the shape of H₂¹⁸O and HDO depth profiles which result from evaporation of water from dry soil under quasi-steady state conditions. Typically, isotope concentrations rise from a minimum at the soil surface to a maximum a short distance beneath the surface, and then decrease approximately exponentially to constant concentrations at depth. The model predicts that for isothermal conditions, the slope of the relationship between ¹⁸O and deuterium δ-values of samples of the soil water will be 30% lower for a dry soil than for a wet soil evaporating under the same conditions. It is concluded that low slopes should be indicative of soil water or groundwater recharged under arid or semi-arid conditions. Using the shape of the ¹⁸O and deuterium profiles, three independent methods for estimating evaporation for dry soils are developed. When water loss occurs by both transpiration and evaporation, the slope of the ¹⁸O-D relationship should be slightly lower than that for a site where water loss occurs by evaporation alone.     

VARIABILITY OF MELTWATER AND SOLUTE FLUXES FROM HOMOGENEOUS MELTING SNOW AT THE LABORATORY SCALE

Four experiments were performed to examine the relationship between the meltwater flow field and ion release from melting snow. A 0.4 m³ volume of snow was placed in a Plexiglass box and melted from above using a heating plate. The meltwater and solute fluxes issuing from the bottom of the snow were monitored. In experiments with NaCl tracer added to the snow, the solute concentrations were generally lower in the flow fingers than in the background wetting front. Dye tracer experiments revealed contemporaneous areas of concentrated dye and dilute meltwater in flow fingers. This suggests that the meltwater in flow fingers is diluted by low concentration water from the top of the snowpack. Flow fingers contribute more meltwater flux primarily because the flow is maintained for a longer period of time than in the non-finger areas; however, the relative contribution of flow fingers to solute flux was apparently not as great as that of the background wetting front because of dilution of solute in the flow finger areas.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

Groundwater modelling with limited data sets: the Chari–Logone area (Lake Chad Basin,Chad)

One of the most important issues for water resource management is developing strategies for groundwater modelling that are adaptable to data scarcity. These strategies are particularly important in arid and semi‐arid areas where access to data is poor and data collection is difficult, such as the Lake Chad Basin in Africa. In the present study, we establish a numerical groundwater flow model and evaluate the effects of dry and wet periods on groundwater recharge in the Chari–Logone area (96 000 km²) of the Lake Chad Basin. Boundary conditions, flow direction, sources, and sinks for the Chari–Logone local model were obtained by revising and remodelling the Lake Chad Basin regional hydrogeological model (508 400 km²) developed by the BRGM (Bureau de Recherches Géologiques et Minières) in the 1990s. The simulated aquifer water level showed good agreement with observed levels. Aquifer recharge is primarily determined by river–aquifer interactions and mostly occurs in the southern section of the study area. In wet years, groundwater recharge also occurs in the N'Djamena area. The approach we adopted provided relevant results and was useful as an initial step in more detailed modelling of the area. It also proved to be a useful method for groundwater modelling in large semi‐arid and arid regions where available data are scarce. Copyright © 2013 John Wiley & Sons, Ltd.     

Using time‐domain reflectometry to characterize shallow solute transport in an oak woodland hillslope in northern California,USA

The natural heterogeneity of water and solute movement in hillslope soils makes it difficult to accurately characterize the transport of surface‐applied pollutants without first gathering spatially distributed hydrological data. This study examined the application of time‐domain reflectometry (TDR) to measure solute transport in hillslopes. Three different plot designs were used to examine the transport of a conservative tracer in the first 50 cm of a moderately sloping soil. In the first plot, which was designed to examine spatial variability in vertical transport in a 1·2 m² plot, a single probe per meter was found to adequately characterize vertical solute travel times. In addition, a dye and excavation study in this plot revealed lateral preferential flow in small macropores and a transport pattern where solute is focused vertically into preferential flow pathways. The bypass flow delivers solute deeper in the soil, where lateral flow occurs. The second plot, designed to capture both vertical and lateral flow, provided additional evidence confirming the flow patterns identified in the excavation of the first plot. The third plot was designed to examine lateral flow and once again preferential flow of the tracer was observed. In one instance rapid solute transport in this plot was estimated to occur in as little as 3% of the available pore space. Finally, it was demonstrated that the soil anisotropy, although partially responsible for lateral subsurface transport, may also homogenize the transport response across the hillslope by decreasing vertical solute spreading. Copyright © 2002 John Wiley & Sons, Ltd.     

Comparison of soil salinity and solute transport for different cultivated soil types in northeastern Egypt / Comparaison de la salinité du sol et du transport de solutés pour différents types de sols cultivés du nord-est de l'Egypte

Abstract

The soil salinity distribution and solute transport properties of three different soil types were investigated and compared within a project area in northeastern Egypt. For this purpose, dye tracer experiments and salinity sampling were carried out. The resulting salinity maps showed that the soil salinity in the cultivated western site of the project area is 8–10 times higher than that in the cultivated eastern site. However, the cultivated soil displayed significantly lower salinity with higher uniformity as compared to the uncultivated soil. The preferential flow phenomenon was less apparent in the cultivated soil. This is mainly due to tillage which disrupts the structure of the soil so that deep cracks are no longer connected to the soil surface. This reduces the risk for groundwater contamination through preferential flow. The study showed that careful and continuous monitoring of the salinity status is needed now and in the future.     

Soil moisture states,lateral flow,and streamflow generation in a semi‐arid,snowmelt‐driven catchment

Hydraulic connectivity on hillslopes and the existence of preferred soil moisture states in a catchment have important controls on runoff generation. In this study we investigate the relationships between soil moisture patterns, lateral hillslope flow, and streamflow generation in a semi‐arid, snowmelt‐driven catchment. We identify five soil moisture conditions that occur during a year and present a conceptual model based on field studies and computer simulations of how streamflow is generated with respect to the soil moisture conditions. The five soil moisture conditions are (1) a summer dry period, (2) a transitional fall wetting period, (3) a winter wet, low‐flux period, (4) a spring wet, high‐flux period, and (5) a transitional late‐spring drying period. Transitions between the periods are driven by changes in the water balance between rain, snow, snowmelt and evapotranspiration. Low rates of water input to the soil during the winter allow dry soil regions to persist at the soil–bedrock interface, which act as barriers to lateral flow. Once the dry‐soil flow barriers are wetted, whole‐slope hydraulic connectivity is established, lateral flow can occur, and upland soils are in direct connection with the near‐stream soil moisture. This whole‐slope connectivity can alter near‐stream hydraulics and modify the delivery of water, pressure, and solutes to the stream. Copyright © 2005 John Wiley & Sons, Ltd.     

Infiltration-recharge through wadi beds in arid regions

Abstract

Groundwater recharge in arid regions is intermittent and usually occurs as a result of flood flow transmission losses in dry wadi channels. Hydrograph characteristics play a dominant role in determining the amount of channel abstraction in relation to the width of the wetted perimeter and the time of inundation, and the subsequent groundwater recharge. Large variations in the magnitude of channel losses result mainly from the diversity in inflow volumes. The magnitude of groundwater recharge in relation to bed transmission losses is dependent on flood volume and duration, soil moisture content and physical soil profile characteristics. Runoff volume and duration are the dominant factors influencing the cumulative infiltrated volume and recharge to shallow water tables. Taking into consideration the influence of various hydrological and channel characteristics, several regression equations are suggested to estimate the transmission losses from a wadi bed and the groundwater recharge.     

The influence of preferential flow on hillslope hydrology in a semi‐arid watershed (in the Spanish Dehesas)

Preferential flow is known to influence hillslope hydrology in many areas around the world. Most research on preferential flow has been performed in temperate regions. Preferential infiltration has also been found in semi‐arid regions, but its impact on the hydrology of these regions is poorly known. The aim of this study is to describe and quantify the influence of preferential flow on the hillslope hydrology from small scale (infiltration) to large scale (subsurface stormflow) in a semi‐arid Dehesa landscape. Precipitation, soil moisture content, piezometric water level and discharge data were used to analyse the hydrological functioning of a catchment in Spain. Variability of soil moisture content during the transition from dry to wet season (September to November) within horizontal soil layers leads to the conclusion that there is preferential infiltration into the soils. When the rainfall intensity is high, a water level rapidly builds up in the piezometer pipes in the area, sometimes even reaching soil surface. This water level also drops back to bedrock within a few hours (under dry catchment conditions) to days (under wet catchment conditions). As the soil matrix is not necessarily wet while this water layer is built up, it is thought to be a transient water table in large connected pores which drain partly to the matrix, partly fill up bedrock irregularities and partly drain through subsurface flow to the channels. When the soil matrix becomes wetter the loss of water from macropores to the matrix and bedrock decreases and subsurface stormflow increases. It may be concluded that the hillslope hydrological system consists of a fine matrix domain and a macropore domain, which have their own flow characteristics but which also interact, depending on the soil matrix and macropore moisture contents. The macropore flow can result in subsurface flow, ranging from 13% contribution to total discharge for a large event of high intensity rainfall or high discharge to 80% of total discharge for a small event with low intensity rainfall or low discharge. During large events the fraction of subsurface stormflow in the discharge is suppressed by the large amount of surface runoff. Copyright © 2008 John Wiley & Sons, Ltd.     

On Evaluating Characteristics of the Solute Transport in the Arid Vadose Zone

The transport of bromide (Br) under matric heads of 0, ?2, ?5, and ?10 cm using undisturbed soil columns was investigated for understanding the solute transport in arid soils. Undisturbed soil cores were collected at ground surface, directly below where tension infiltrometer measurements were made in the Amargosa Desert, Nevada, United States. Laboratory experiments were conducted by introducing water containing Br tracer into a soil column maintained at steady‐state conditions. The observed data of breakthrough curves (BTC) were well fitted to an one‐region model, except for the cores at saturation, and a core at the matric head of ?5 cm, from which the observed data were better fitted to a two‐region model. Fitted pore water velocities with the one‐region model ranged from 1.2 to 56.6 cm/h, and fitted dispersion coefficients (D) ranged from 2.2 to 100 cm²/h. Results for the core analyzed with the two‐region model indicated that D ranged from 27.6 to 70.9 cm²/h at saturation, and 25.7 cm²/h at the matric head of ?5 cm; fraction of mobile water (β) ranged from 0.18 to 0.65, and mass transfer coefficient (ω) ranged from 0.006 to 0.03. In summary, the water fluxes and Br dispersion coefficients at investigated matric heads were very high due to the coarseness of the soils and possibly due to preferential flow pathways. These high water fluxes and Br dispersion coefficients would lead to a higher risk of deeper leaching accumulating nitrate nitrogen to the groundwater, and have significant effects on the desert ecosystem.