scholarly journals Interplay between Fingering Instabilities and Initial Soil Moisture in Solute Transport through the Vadose Zone

Water ◽  
2020 ◽  
Vol 12 (3) ◽  
pp. 917
Author(s):  
Luis Cueto-Felgueroso ◽  
María José Suarez-Navarro ◽  
Xiaojing Fu ◽  
Ruben Juanes
Keyword(s):  
Spatial Distribution ◽  
Solute Transport ◽  
Soil Water ◽  
Water Flow ◽  
Vadose Zone ◽  
Preferential Flow ◽  
Solute Diffusion ◽  
Wetting Front ◽  
Initial Soil ◽  
The Impact

Modeling water flow and solute transport in the vadose zone is essential to understanding the fate of soil pollutants and their travel times towards groundwater bodies. It also helps design better irrigation strategies to control solute concentrations and fluxes in semiarid and arid regions. Heterogeneity, soil texture and wetting front instabilities determine the flow patterns and solute transport mechanisms in dry soils. When water is already present in the soil, the flow of an infiltration pulse depends on the spatial distribution of soil water and on its mobility. We present numerical simulations of passive solute transport during unstable infiltration of water into sandy soils that are prone to wetting front instability. We study the impact of the initial soil state, in terms of spatial distribution of water content, on the infiltration of a solute-rich water pulse. We generate random fields of initial moisture content with spatial structure, through multigaussian fields with prescribed correlation lengths. We characterize the patterns of water flow and solute transport, as well as the mass fluxes through the soil column. Our results indicate a strong interplay between preferential flow and channeling due to fingering and the spatial distribution of soil water at the beginning of infiltration. Fingering and initial water saturation fields have a strong effect on solute diffusion and dilution into the ambient water during infiltration, suggesting an effective separation between mobile and inmobile transport domains that are controlled by the preferential flow paths due to fingering.

Agricultural hillslope soil heterogeneity challenges: first experimental results from SUPREHILL critical zone observatory

2021 ◽  
Author(s):  
Vedran Krevh ◽  
Jasmina Defterdarović ◽  
Lana Filipović ◽  
Zoran Kovač ◽  
Steffen Beck-Broichsitter ◽  
...  
Keyword(s):  
Numerical Modeling ◽  
Soil Moisture ◽  
Solute Transport ◽  
Soil Water ◽  
Water Flow ◽  
Preferential Flow ◽  
Soil Heterogeneity ◽  
Local Scale ◽  
Soil Water Flow ◽  
Critical Zone Observatory

<p>SUPREHILL is a new (2020) and first Croatian critical zone observatory (CZO), focused on local scale agricultural e.g., vineyard hillslope processes. The experimental setup includes an extensive sensor-based network accompanied by weighing lysimeters and instruments for surface and subsurface hydrology measurement. The field measurements are supported by novel laboratory and numerical quantification methods for the determination of water flow and solute transport. This combined approach will allow the research team to accurately determine soil water balance components (soil water flow, preferential flow/transport pathways, surface runoff, evapotranspiration), the temporal origin of water in hillslope hydrology (isotopes), transport of agrochemicals, and to calibrate and validate numerical modeling procedures for describing and predicting soil water flow and solute transport. First results from sensors indicate increased soil moisture on the hilltop, which is supported by precipitation data from rain gauges and weighing lysimeters. The presence of a compacted soil horizon and compacted inter-row parts (due to trafficking) of the vineyard seems to be highly relevant in regulating water dynamics. Wick lysimeters confirm the sensor soil moisture data, while showing a significant difference in its repetitions which suggests a possibility of a preferential flow imposed by local scale soil heterogeneity. Measured values of surface and subsurface runoff suggest a crucial role of these processes in the hillslope hydrology, while slope and structure dynamics additionally influence soil hydraulic properties. We are confident that the CZO will give us new insights in the landscape heterogeneity and substantially increase our understanding regarding preferential flow and nonlinear solute transport, with results directly applicable in agricultural (sloped agricultural soil management) and environmental (soil and water) systems. Challenges remain in characterizing local scale soil heterogeneity, dynamic properties quantification and scaling issues for which we will rely on combining CZO focused measurements and numerical modeling after substantial data is collected.</p>

Influence of the shape of inter-horizon boundary and size of soil tongues on preferential flow under shallow groundwater conditions: A simulation study

2011 ◽  
Vol 91 (2) ◽  
pp. 211-221 ◽  
Author(s):  
Priyantha B. Kulasekera ◽  
Gary W. Parkin
Keyword(s):  
Solute Transport ◽  
Vadose Zone ◽  
Simulation Study ◽  
Soil Profile ◽  
Preferential Flow ◽  
Shallow Groundwater ◽  
Soil Horizon ◽  
Tongue Length ◽  
The Impact ◽  
Water And Solute Transport

Kulasekera, P. B. and Parkin, G. W. 2011. Influence of the shape of inter-horizon boundary and size of soil tongues on preferential flow under shallow groundwater conditions: A simulation study. Can. J. Soil Sci. 91: 211–221. Detailed studies of the impact of soil tongues at soil horizon interfaces are very important in understanding preferential flow processes through layered soils and in improving the accuracy of models predicting water and solute transport through the vadose zone. The implication of having soil tongues of different shapes and sizes created at the soil horizon interface on solute transport through a layered soil horizon was studied by simulating water and solute transport using the VS2DI model. This 2-D simulation study reconfirmed that soil tongues facilitate preferential flow, and the level of activeness of tongues may depend on the number of soil tongues, their spacing and distribution. Also, the size of the soil tongues (length and diameter at the interface between the soil horizons) and their shape influence the rate of preferential flow. Increasing tongue length consistently resulted in an increase in solute velocity across the entire soil profile regardless of the tongue shape; for example, a soil tongue of 0.25 m length increased solute velocity by about 1.5 times over a soil profile without tongues, but this increase might be different for soil types and groundwater conditions other than those considered in this study. Narrowing of tongues increased solute velocity, whereas increasing the number of tongues in a wider soil profile decreased the solute-front's velocity. As tongue length increased, the area containing solutes at prescribed elapsed times decreased. An implication of this study is that soil horizon tongue shape and spacing reduce pollutant residence times, hence inter-horizon boundary morphology should be considered when modelling transport through the vadose zone. As well, since the solute velocity behaviours of a triangular- and a wider rectangular-shaped tongue were nearly identical, simply measuring solute velocity in the field will reveal little information on the shape of a soil tongue.

scholarly journals Simulating preferential soil water flow and tracer transport using the Lagrangian Soil Water and Solute Transport Model

2019 ◽  
Vol 23 (10) ◽  
pp. 4249-4267 ◽  
Author(s):  
Alexander Sternagel ◽  
Ralf Loritz ◽  
Wolfgang Wilcke ◽  
Erwin Zehe
Keyword(s):  
Solute Transport ◽  
Soil Water ◽  
Water Flow ◽  
Preferential Flow ◽  
Tracer Transport ◽  
Model Concept ◽  
Soil Matrix ◽  
Matrix Domain ◽  
Flow Domain ◽  
Water And Solute Transport

Abstract. We propose an alternative model concept to represent rainfall-driven soil water dynamics and especially preferential water flow and solute transport in the vadose zone. Our LAST-Model (Lagrangian Soil Water and Solute Transport) is based on a Lagrangian perspective of the movement of water particles (Zehe and Jackisch, 2016) carrying a solute mass through the subsurface which is separated into a soil matrix domain and a preferential flow domain. The preferential flow domain relies on observable field data like the average number of macropores of a given diameter, their hydraulic properties and their vertical length distribution. These data may be derived either from field observations or by inverse modelling using tracer data. Parameterization of the soil matrix domain requires soil hydraulic functions which determine the parameters of the water particle movement and particularly the distribution of flow velocities in different pore sizes. Infiltration into the matrix and the macropores depends on their respective moisture state, and subsequently macropores are gradually filled. Macropores and matrix interact through diffusive mixing of water and solutes between the two flow domains, which again depends on their water content and matric potential at the considered depths. The LAST-Model is evaluated using tracer profiles and macropore data obtained at four different study sites in the Weiherbach catchment in southern Germany and additionally compared against simulations using HYDRUS 1-D as a benchmark model. While both models show qual performance at two matrix-flow-dominated sites, simulations with LAST are in better accordance with the fingerprints of preferential flow at the two other sites compared to HYDRUS 1-D. These findings generally corroborate the feasibility of the model concept and particularly the implemented representation of macropore flow and macropore–matrix exchange. We thus conclude that the LAST-Model approach provides a useful and alternative framework for (a) simulating rainfall-driven soil water and solute dynamics and fingerprints of preferential flow as well as (b) linking model approaches and field experiments. We also suggest that the Lagrangian perspective offers promising opportunities to quantify water ages and to evaluate travel and residence times of water and solutes by a simple age tagging of particles entering and leaving the model domain.

Pines influence hydrophysical parameters and water flow in a sandy soil

Biologia ◽  
2013 ◽  
Vol 68 (6) ◽  
Author(s):  
Ľubomír Lichner ◽  
Jozef Capuliak ◽  
Natalia Zhukova ◽  
Ladislav Holko ◽  
Henryk Czachor ◽  
...  
Keyword(s):  
Hydraulic Conductivity ◽  
Forest Soil ◽  
Soil Water ◽  
Water Flow ◽  
Sandy Soil ◽  
Preferential Flow ◽  
Water Repellency ◽  
Water Repellent ◽  
Wetting Front ◽  
Soil Water Repellency

AbstractPines, used for sand dune stabilization, can influence the hydrophysical parameters and water flow in an aeolian sandy soil considerably, mainly due to soil water repellency. Two sites, separated by distance of about 20 m, formed the basis of our study. A control soil (“Pure sand“) with limited impact of vegetation or organic matter was formed at 50 cm depth beneath a forest glade area. This was compared to a “Forest soil” in a 30-year old Scots pine (Pinus sylvestris) forest. Most of the hydrophysical parameters were substantially different between the two soil surfaces. The forest soil was substantially more water repellent and had two-times the degree of preferential flow compared to pure sand. Water and ethanol sorptivities, hydraulic conductivity, and saturated hydraulic conductivity were 1%, 84%, 2% and 26% those of the pure sand, respectively. The change in soil hydrophysical parameters due to soil water repellency resulted in preferential flow in the forest soil, emerging during a simulated heavy rain following a long hot, dry period. The wetting front established in pure sand exhibited a form typical of that for stable flow. Such a shape of the wetting front can be expected in the forest soil in spring, when soil water repellency is alleviated substantially.

scholarly journals Simulating preferential soil water flow and tracer transport using the Lagrangian Soil Water and Solute Transport Model

2019 ◽  
Author(s):  
Alexander Sternagel ◽  
Ralf Loritz ◽  
Wolfgang Wilcke ◽  
Erwin Zehe
Keyword(s):  
Solute Transport ◽  
Soil Water ◽  
Water Flow ◽  
Field Experiments ◽  
Preferential Flow ◽  
Sensitivity Analyses ◽  
Tracer Transport ◽  
Soil Matrix ◽  
Soil Water Flow ◽  
Water And Solute Transport

Abstract. We propose an alternative model to overcome these weaknesses of the Darcy-Richards approach and to simulate preferential soil water flow and tracer transport in macroporous soils. Our LAST-Model (Lagrangian Soil Water and Solute Transport) relies on a Lagrangian perspective on the movement of water particles carrying a solute mass through the soil matrix and macropores. We advance the model of Zehe and Jackisch (2016) by two main extensions: a) a new routine for solute transport within the soil matrix and b) the implementation of an additional 1-D preferential flow domain which simulates flow and transport in a population of macropores. Infiltration into the matrix and the macropores depends on their moisture state and subsequently macropores are gradually filled. Macropores and matrix interact through diffusive mixing of water and solutes between the two domains which depends on their water content and matric potential at the considered depths. The LAST-Model is then evaluated with sensitivity analyses and tested against tracer field experiments at three different sites. The results show the internal and physical validity of the model and the robustness of our solute transport and the linear mixing approach. Further, the model is able to simulate preferential flow through macropores and to depict well the observed 1-D solute mass profile of a tracer experiment with a high computational efficiency and short simulation times.

Simulating preferential soil water flow and reactive solute transport using the Lagrangian Soil Water and Solute Transport Model (LAST)

2020 ◽  
Author(s):  
Alexander Sternagel ◽  
Ralf Loritz ◽  
Wolfgang Wilcke ◽  
Erwin Zehe
Keyword(s):  
Solute Transport ◽  
Soil Water ◽  
Water Flow ◽  
Preferential Flow ◽  
Transport Model ◽  
Water Dynamics ◽  
Model Concept ◽  
Soil Matrix ◽  
Reactive Solute Transport ◽  
Water And Solute Transport

<p>Recently, we proposed an alternative model concept to represent rainfall-driven soil water dynamics and especially preferential water flow and solute transport in the vadose zone. Our LAST-Model is based on a Lagrangian perspective on the movement of water particles (Zehe and Jackisch, 2016) carrying solute masses through the subsurface which is separated into a soil matrix domain and a preferential flow domain (Sternagel et al., 2019). The preferential flow domain relies on observable field data like the average number of macropores of a given diameter, their hydraulic properties and their vertical length distribution. These data may either be derived from field observations or by inverse modelling using tracer data. Parameterization of the soil matrix domain requires soil hydraulic functions which determine the parameters of the water particle movement and particularly the distribution of flow velocities in different pores sizes. Infiltration into the matrix and the macropores depends on their respective moisture state and subsequently macropores are gradually filled. Macropores and matrix interact through diffusive mixing of water and solutes between the two flow domains which again depends on their water content and matric potential at the considered depths.</p><p>The LAST-Model was evaluated using tracer profiles and macropore data obtained at four different study sites in the Weiherbach catchment in south Germany and additionally compared against simulations using HYDRUS 1-D as benchmark model. The results generally corroborated the feasibility of the model concept and particularly the implemented representation of macropore flow and macropore-matrix exchange. We thus concluded that the LAST-Model approach provides a useful and alternative framework for simulating rainfall-driven soil water and solute dynamics and fingerprints of preferential flow.</p><p>This study presents an extension of the model allowing for the simulation of reactive solute transport. Transformation kinetics are considered by transferring mass from the parent to the child components in each water particle according to the corresponding reaction rates, which is limited by the compound solubility. A retardation coefficient is not helpful in the particle-based framework, as the solute mass is carried by the water particles and travels thus by default at the same velocity. Ad- and desorption are explicit represented through transfer of dissolved mass from the water particles at a given depth to surrounding adsorption sites of the soil solid phase and vice versa. This may either operate under rate-limited or non-limited conditions. Adsorbed solute masses will be considered to be degraded following first-order reaction kinetics. The retardation process delays the solute displacement and enables a suitable time scale for the degradation process, which must be smaller than the time scale for the re-mobilization of the solutes. The proposed extension will be benchmarked against observations of pesticide transport in soil profiles and at tile-drained field sites.</p><p> </p><p>Zehe, E., Jackisch, C.: A Lagrangian model for soil water dynamics during rainfall-driven conditions, Hydrol. Earth Syst. Sci., 20, 3511–3526, https://doi.org/10.5194/hess-20-3511-2016, 2016.</p><p> </p><p>Sternagel, A., Loritz, R., Wilcke, W., and Zehe, E.: Simulating preferential soil water flow and tracer transport using the Lagrangian Soil Water and Solute Transport Model, Hydrol. Earth Syst. Sci., 23, 4249–4267, https://doi.org/10.5194/hess-23-4249-2019, 2019.</p>

Historical development of soil-water physics and solute transport in porous media

2007 ◽  
Vol 7 (1) ◽  
pp. 59-66 ◽  
Author(s):  
D.E. Rolston
Keyword(s):  
Porous Media ◽  
Hydraulic Conductivity ◽  
Solute Transport ◽  
Soil Water ◽  
20Th Century ◽  
Water Flow ◽  
Unsaturated Soils ◽  
Transport Processes ◽  
Transport In Porous Media ◽  
Unsaturated Hydraulic Conductivity

The science of soil-water physics and contaminant transport in porous media began a little more than a century ago. The first equation to quantify the flow of water is attributed to Darcy. The next major development for unsaturated media was made by Buckingham in 1907. Buckingham quantified the energy state of soil water based on the thermodynamic potential energy. Buckingham then introduced the concept of unsaturated hydraulic conductivity, a function of water content. The water flux as the product of the unsaturated hydraulic conductivity and the total potential gradient has become the accepted Buckingham-Darcy law. Two decades later, Richards applied the continuity equation to Buckingham's equation and obtained a general partial differential equation describing water flow in unsaturated soils. For combined water and solute transport, it had been recognized since the latter half of the 19th century that salts and water do not move uniformly. It wasn't until the middle of the 20th century that scientists began to understand the complex processes of diffusion, dispersion, and convection and to develop mathematical formulations for solute transport. Knowledge on water flow and solute transport processes has expanded greatly since the early part of the 20th century to the present.

Characterization of hydraulic properties in the upper vadose zone for two aquifers in Slovenia

2021 ◽  
Author(s):  
Vesna Zupanc ◽  
Matjaž Glavan ◽  
Miha Curk ◽  
Urša Pečan ◽  
Michael Stockinger ◽  
...  
Keyword(s):  
Stable Isotopes ◽  
Soil Water ◽  
Water Flow ◽  
Vadose Zone ◽  
Hydraulic Properties ◽  
Isotope Composition ◽  
Transport Processes ◽  
Precipitation Event ◽  
Water Fluxes ◽  
Flow And Transport

<p>Environmental tracers, present in the environment and provided by nature, provide integrative information about both water flow and transport. For studying water flow and solute transport, the hydrogen and oxygen isotopes are of special interest, as their ratios provide a tracer signal with every precipitation event and are seasonally distributed. In order to follow the seasonal distribution of stable isotopes in the soil water and use this information for identifying hydrological processes and hydraulic properties, soil was sampled three times in three profiles, two on Krško polje aquifer in SE Slovenia and one on Ljubljansko polje in central Slovenia. Isotope composition of soil water was measured with the water-vapor-equilibration method. Based on the isotope composition of soil water integrative information about water flow and transport processes with time and depth below ground were assessed. Porewater isotopes were in similar range as precipitation for all three profiles.  Variable isotope ratios in the upper 60 cm for the different sampling times indicated dynamic water fluxes in this upper part of the vadose zone. Results also showed more evaporation at one sampling location, Brege. The information from stable isotopes will be of importance for further analyzing the water fluxes in the vadose zone of the study sties. <br>This research was financed by the ARRS BIAT 20-21-32 and IAEA CRP 1.50.18 Multiple isotope fingerprints to identify sources and transport of agro-contaminants.  </p>

Sign in / Sign up

Export Citation Format

Share Document