Comparative Non-Darcian Modeling of Subsurface Preferential Flow Experimental Observations in a Riparian Buffer

2021 ◽  
Vol 64 (6) ◽  
pp. 1867-1881
Author(s):  
Enrique Orozco-López ◽  
Rafael Muñoz-Carpena
Keyword(s):  
Vadose Zone ◽  
Preferential Flow ◽  
Riparian Buffers ◽  
Riparian Buffer ◽  
Model Parameters ◽  
List Type ◽  
Dispersive Wave ◽  
Macropore Flow ◽  
Modeling Framework ◽  
Soil Matrix

HighlightsHigh ecohydrological activity drives macropore prevalence in riparian buffers.An abundance of macropore flow (MF) was confirmed in a field riparian buffer in Kenya.Source-response (SR) and multilayer kinematic diffusive wave (MKDW) MF models are compared.A novel MKDW modeling framework efficiently identifies and predicts preferential flow in riparian buffers.Abstract. The significant ecohydrological activity typical of riparian buffers makes them potential hotspots of macropores, i.e., structured preferential flow pathways, through the soil vadose zone. The prevalence of these preferential pathways can allow transported contaminants to bypass the soil matrix and quickly reach a seasonal shallow water table and the adjacent surface waterbody. This quick transport can ultimately limit the role of riparian buffers for runoff pollution control. Currently, there are no management tools that incorporate macropore flow (MF) when assessing riparian buffer performance. The objective of this study was to experimentally quantify and mathematically simulate macropore flow and arrival time in a riparian buffer under field conditions. Three infiltration experiments were conducted with a grid of 20 time-domain transmission (TDT) dielectric soil moisture sensors along a field riparian buffer transect in Kenya to quantify the presence of macropore flow and to test two non-Darcian soil MF models, including the source-responsive (SR) model and the modified kinematic-dispersive wave (MKDW) model developed in this study, by adding a user-defined multilayer convection scheme and a new hysteresis function between water flux and content. The abundance of MF in the riparian buffer was corroborated experimentally. Modeling results showed that the MKDW model was an efficient (average NSE of 0.937 and 0.721 for calibration and testing, respectively), flexible, and robust method to identify and represent non-linear and non-sequential MF signals at any soil depth and antecedent conditions. The SR model was computationally inexpensive and provided good calibration results (NSE = 0.867) but required piecemeal recalibration of the travel time and maximum water content at each layer and yielded lower performance in testing. The Akaike (AIC) and Bayesian (BIC) information criteria showed that MKDW outperformed SR when accounting for the trade-off between model complexity and efficiency. The results support further research focused on independent characterization of model parameters at the field scale, and the inclusion of MKDW in holistic riparian buffer management and decision-support tools such as VFSmod. Keywords: Kinematic-dispersive wave, Macropore flow, Numerical modeling, Preferential flow, Riparian vadose zone.

scholarly journals Quantifying the Importance of Preferential Flow in a Riparian Buffer

2021 ◽  
Vol 64 (3) ◽  
pp. 937-947
Author(s):  
Lucie Guertault ◽  
Garey A. Fox ◽  
Todd Halihan ◽  
Rafael Muñoz-Carpena
Keyword(s):  
Preferential Flow ◽  
Riparian Buffers ◽  
Riparian Buffer ◽  
Design Tools ◽  
Soil Matrix ◽  
Innovative Method ◽  
Matrix Flow ◽  
Buffer Design ◽  
Flow Pathways ◽  
Preferential Flow Pathways

HighlightsRiparian buffers and vegetative filter strips are uniquely susceptible to preferential flow.An innovative method is proposed to partition infiltration into matrix and macropore domains.Riparian buffer matrix and plot-scale infiltration experiments were simulated with HYDRUS-1D and VFSMOD.Preferential flow accounted for 32% to 47% of infiltration depending on hydrologic conditions.Preferential flow mechanisms should be incorporated into riparian buffer design tools and models.Abstract. Riparian buffers are uniquely susceptible to preferential flow due to the abundance of root channels, biological activity, and frequent wetting and drying cycles. Previous research has indicated such susceptibility and even measured the connectivity of preferential flow pathways with adjacent streams and rivers. However, limited research has attempted to partition the riparian buffer infiltration between matrix and preferential flow domains. The objectives of this research were to develop an innovative method to quantify soil matrix infiltration at the plot scale, develop a method to partition infiltration into matrix and macropore infiltration at the plot scale, and then use these methods to quantify the significance of macropore infiltration at a riparian buffer site. This research further demonstrated the importance of considering preferential flow processes in design tools and models to evaluate riparian buffer effectiveness. Sprinkler and runon field experiments were conducted at an established riparian buffer site with sandy loam soil. Trenches were installed and instrumented with soil moisture sensors along the width of the riparian buffer (i.e., along the flow path toward the stream) for detecting non-uniform flow patterns due to preferential flow. Riparian buffer parameters, including soil hydraulic parameters, were estimated using HYDRUS-1D for the sprinkler experiments and VFSMOD for the runon experiments. This research partitioned the infiltration into matrix and preferential flow domains by assuming negligible exchange of water between the soil matrix and preferential flow pathways in comparison to the magnitude of soil matrix flow. For these experimental conditions with 0.20 to 0.48 L s-1 of runon and initial soil water contents of 0.29 to 0.32 cm3 cm-3, preferential flow accounted for at least 27% to 32% of the total runon water entering the riparian buffer. This corresponded to approximately 32% to 47% of the total infiltration. While increasing the riparian buffer plot soil hydraulic conductivity in single-porosity models can adequately predict the total infiltration and therefore the surface outflow from the buffer, design tools and models should specifically consider preferential flow processes to improve predictive power regarding the actual infiltration processes and correspondingly the non-equilibrium flow and solute transport mechanisms. Keywords: Flow partitioning, HYDRUS, Matrix flow, Preferential flow, Riparian buffer, VFSMOD.

scholarly journals Perspective: Preferential Flow in Riparian Buffers: Current Research and Future Needs

2021 ◽  
Vol 64 (6) ◽  
pp. 1907-1911
Author(s):  
Derek M. Heeren ◽  
Lucie Guertault ◽  
Kyle Mankin
Keyword(s):  
Field Experiments ◽  
Preferential Flow ◽  
Light Transmission ◽  
Management Practice ◽  
Fate And Transport ◽  
Riparian Buffers ◽  
Riparian Buffer ◽  
List Type ◽  
Buffer Strip ◽  
Future Needs

HighlightsPreferential flow (PF) can critically reduce riparian buffer contaminant removal efficiency.This collection presents research on PF measurement, visualization, modeling, and contaminant transport impacts.Future needs include tools to identify landscape-scale PF areas and conservation practices.Future models for research and practice should account for PF in riparian buffers.Abstract. Preferential flow in riparian buffers can substantially compromise their effectiveness in reducing contaminants from overland runoff. The objective of this article is to introduce a collection of five articles on current research into subsurface preferential flow measurement, visualization, modeling, and impacts on contaminant fate and transport at scales ranging from the subsurface pore scale to the plot scale to the watershed scale. This collection presents selected works from a broader invited session on “Preferential flow and piping in riparian buffers” at the 2020 ASABE Annual International Meeting. Major findings include: new methodologies, such as light transmission and geophysics, to characterize subsurface preferential flow; an infiltration partitioning approach to quantify preferential flow from field experiments; a kinematic dispersive wave model to effectively simulate subsurface preferential flow; and the significant impact of surface concentrated flow pathways on pesticide fate and transport both upstream and within a riparian buffer. Future work is needed to develop methods and tools to identify PF areas and management solutions within a landscape, and to update both research and design models to better quantify and account for PF processes. Keywords: Best management practice, Buffer strip, Agricultural conservation practice, Filter strip, Macropore, Nonpoint-source pollution.

scholarly journals Detecting Macropore Fingering Using Temporal Electrical Resistivity Imaging

2021 ◽  
Vol 37 (5) ◽  
pp. 861-870
Author(s):  
Todd Halihan ◽  
John Paul Hager ◽  
Lucie Guertault ◽  
Garey A. Fox
Keyword(s):  
Electrical Resistivity ◽  
Preferential Flow ◽  
Field Evaluation ◽  
Low Noise ◽  
Riparian Buffers ◽  
Electrical Resistivity Imaging ◽  
List Type ◽  
Soil Matrix ◽  
Resistivity Imaging ◽  
Field Setting

HighlightsSingle macropores can be detected using temporal electrical resistivity imaging under controlled conditions.Macropore flow can be detected based on preferentially wetted fingers of increased conductance.Macropore activation does not appear to require saturated surface conditions to induce preferential flow.Abstract. Riparian soils are uniquely susceptible to the formation of macropores, which are hypothesized to promote fast transport of water and contaminants through upper soil layers. Electrical Resistivity Imaging (ERI) can locate spatial heterogeneities in soil wetting patterns and evaluate differences due to vegetation, thus optimizing the design of riparian buffers. Temporal ERI (TERI) imaging was conducted in a fine and coarse field setting with artificial macropores to evaluate flow under unsaturated simulated rainfall conditions and saturated infiltrometer conditions. While single macropores are detectable using TERI datasets, the results in an average field setting would detect the wetted area surrounding a macropore, not the macropore itself. The results were similar for both the primary fine grain soil site in Oklahoma as well as the coarse grain site in North Carolina. TERI data indicated that without artificial conditions with low noise conditions, a single macropore would not be detected, a wetted zone would be the best detection. In ordinary field evaluation of natural macropores, the TERI technique would detect the wetted zone around a macropore similar to a high hydraulic conductivity zone in a heterogeneous soil matrix. Finally, the results confirmed that macropore activation does not require saturated conditions to generate preferential flow. Keywords: Hydrogeophysics, Preferential flowpaths, Riparian buffers, Temporal electrical resistivity imaging.

DRAINMOD-P: A Model for Simulating Phosphorus Dynamics and Transport in Drained Agricultural Lands: II. Model Testing

2021 ◽  
Vol 64 (6) ◽  
pp. 1849-1866
Author(s):  
Manal H. Askar ◽  
Mohamed A. Youssef ◽  
Dean L. Hesterberg ◽  
Kevin W. King ◽  
Aziz Amoozegar ◽  
...  
Keyword(s):  
Surface Runoff ◽  
Preferential Flow ◽  
Model Performance ◽  
Water Quality Modeling ◽  
List Type ◽  
Macropore Flow ◽  
Major Pathway ◽  
Phosphorus Loss ◽  
Particulate Form ◽  
Subsurface Drains

HighlightsDRAINMOD-P was tested using a dataset from a drained field with desiccation cracks.Surface and subsurface phosphorus losses were mainly in the particulate form.Surface runoff was a major pathway for phosphorus loss in this field.The model performance in predicting edge-of-field phosphorus loss is promising.Abstract. The recently developed phosphorus (P) model DRAINMOD-P was tested using a four-year dataset from a subsurface-drained field in northwest Ohio with significant potential for desiccation cracking or preferential flow. The model satisfactorily predicted subsurface drainage discharge, with a monthly Nash-Sutcliffe efficiency (NSE) of 0.59 and index of agreement (IOA) of 0.89. Lack of annual water budget closure was reported and was likely caused by uncertainty in measured surface runoff and/or modeling approaches representing macropore flow. More than 80% of predicted surface and subsurface P losses were in the particulate form. Surface runoff was the major pathway for P loss, contributing 78% of predicted total P (TP) load. On average, predicted macropore flow represented about 15% of drainage discharge and contributed 21% of DRP loss via subsurface drains. The performance of DRAINMOD-P in predicting monthly dissolved reactive P and TP losses through subsurface drains can be rated as poor (NSE = 0.33 and IOA = 0.60) and very good (NSE = 0.81 and IOA = 0.95), respectively. DRAINMOD-P demonstrated potential for simulating P fate and transport in drained cropland. More testing is needed to further examine newly incorporated hydrological and biogeochemical components of the model. Keywords: Agricultural drainage, Edge-of-field phosphorus load, Macropore flow, Phosphorus model, Sediment yield, Water quality modeling.

scholarly journals Effects of preferential flow on snowmelt partitioning and groundwater recharge in frozen soils

2019 ◽  
Vol 23 (12) ◽  
pp. 5017-5031 ◽  
Author(s):  
Aaron A. Mohammed ◽  
Igor Pavlovskii ◽  
Edwin E. Cey ◽  
Masaki Hayashi
Keyword(s):  
Soil Moisture ◽  
Groundwater Recharge ◽  
Preferential Flow ◽  
Frozen Soil ◽  
Time Lags ◽  
Runoff Generation ◽  
Frozen Ground ◽  
Soil Matrix ◽  
Frozen Soils ◽  
Flow Through

Abstract. Snowmelt is a major source of groundwater recharge in cold regions. Throughout many landscapes snowmelt occurs when the ground is still frozen; thus frozen soil processes play an important role in snowmelt routing, and, by extension, the timing and magnitude of recharge. This study investigated the vadose zone dynamics governing snowmelt infiltration and groundwater recharge at three grassland sites in the Canadian Prairies over the winter and spring of 2017. The region is characterized by numerous topographic depressions where the ponding of snowmelt runoff results in focused infiltration and recharge. Water balance estimates showed infiltration was the dominant sink (35 %–85 %) of snowmelt under uplands (i.e. areas outside of depressions), even when the ground was frozen, with soil moisture responses indicating flow through the frozen layer. The refreezing of infiltrated meltwater during winter melt events enhanced runoff generation in subsequent melt events. At one site, time lags of up to 3 d between snow cover depletion on uplands and ponding in depressions demonstrated the role of a shallow subsurface transmission pathway or interflow through frozen soil in routing snowmelt from uplands to depressions. At all sites, depression-focused infiltration and recharge began before complete ground thaw and a significant portion (45 %–100 %) occurred while the ground was partially frozen. Relatively rapid infiltration rates and non-sequential soil moisture and groundwater responses, observed prior to ground thaw, indicated preferential flow through frozen soils. The preferential flow dynamics are attributed to macropore networks within the grassland soils, which allow infiltrated meltwater to bypass portions of the frozen soil matrix and facilitate both the lateral transport of meltwater between topographic positions and groundwater recharge through frozen ground. Both of these flow paths may facilitate preferential mass transport to groundwater.

scholarly journals Prediction of leaching and groundwater contamination by pesticides

2006 ◽  
Vol 78 (5) ◽  
pp. 1081-1090 ◽  
Author(s):  
Werner Kördel ◽  
Michael Klein
Keyword(s):  
Maximum Concentration ◽  
Management Practices ◽  
Preferential Flow ◽  
Soil Depth ◽  
Sandy Soils ◽  
Field Studies ◽  
Macropore Flow ◽  
Eu Directive ◽  
Present Information ◽  
Silty Soils

Herein, we describe how pesticide leaching is assessed in Europe in order to fulfill EU Directive 91/414. The assessment schemes were developed to protect groundwater from unacceptable effects caused by pesticide use. They presently focus on chromatographic flow processes, which are dominant in sandy soils. Nevertheless, important regions in Europe are characterized by structured soils where transport through macropores is most relevant.Comparison of parallel field studies with isoproturon performed in sandy and silty soils showed that maximum concentration in the structured soil at a soil depth of 1 m may exceed respective concentrations in sandy soils by a factor of 60. Similar results were obtained by lysimeter studies using silty soil cores with maximum concentration of 40 μg/l at the soil bottom. These results demonstrate that preferential flow is more the rule than the exception in well-structured fine-textured soils, and pesticide losses via macropore flow may exceed losses via matrix transport considerably. All present information available for macropore flow suggest the need for greater regional assessments. Other recommendations include analysis of the influence of different soil management practices on the formation of macropores.

Dynamic Modeling of Channel Formation During Fluid Injection Into Unconsolidated Formations

SPE Journal ◽  
2015 ◽  
Vol 20 (04) ◽  
pp. 689-700 ◽  
Author(s):  
S.. Ameen ◽  
A. Dahi Taleghani
Keyword(s):  
Preferential Flow ◽  
Volume Fraction ◽  
Injection Pressure ◽  
Internal Erosion ◽  
Fluid Injection ◽  
Critical Threshold ◽  
Model Parameters ◽  
Flow Paths ◽  
Preferential Flow Paths ◽  
Stress Changes

Summary Injectivity loss is a common problem in unconsolidated-sand formations. Injection of water into a poorly cemented granular medium may lead to internal erosion, and consequently formation of preferential flow paths within the medium because of channelization. Channelization in the porous medium might occur when fluid-induced stresses become locally larger than a critical threshold and small grains are dislodged and carried away; hence, porosity and permeability of the medium will evolve along the induced flow paths. Vice versa, flowback during shut-in might carry particles back to the well and cause sand accumulation inside the well, and subsequently loss of injectivity. In most cases, to maintain the injection rate, operators will increase injection pressure and pumping power. The increased injection pressure results in stress changes and possibly further changes in channel patterns around the wellbore. Experimental laboratory studies have confirmed the presence of the transition from uniform Darcy flow to a fingered-pattern flow. To predict these phenomena, a model is needed to fill this gap by predicting the formation of preferential flow paths and their evolution. A model based on the multiphase-volume-fraction concept is used to decompose porosity into mobile and immobile porosities where phases may change spatially, evolve over time, and lead to development of erosional channels depending on injection rates, viscosity, and rock properties. This model will account for both particle release and suspension deposition. By use of this model, a methodology is proposed to derive model parameters from routine injection tests by inverse analysis. The proposed model presents the characteristic behavior of unconsolidated formation during fluid injection and the possible effect of injection parameters on downhole-permeability evolution.

scholarly journals Dye staining and excavation of a lateral preferential flow network

2009 ◽  
Vol 13 (6) ◽  
pp. 935-944 ◽  
Author(s):  
A. E. Anderson ◽  
M. Weiler ◽  
Y. Alila ◽  
R. O. Hudson
Keyword(s):  
Slope Stability ◽  
Solute Transport ◽  
Preferential Flow ◽  
Runoff Generation ◽  
Flow Paths ◽  
Soil Matrix ◽  
Preferential Flow Paths ◽  
Dye Staining ◽  
Hillslope Scale ◽  
Surrounding Soil

Abstract. Preferential flow paths have been found to be important for runoff generation, solute transport, and slope stability in many areas around the world. Although many studies have identified the particular characteristics of individual features and measured the runoff generation and solute transport within hillslopes, very few studies have determined how individual features are hydraulically connected at a hillslope scale. In this study, we used dye staining and excavation to determine the morphology and spatial pattern of a preferential flow network over a large scale (30 m). We explore the feasibility of extending small-scale dye staining techniques to the hillslope scale. We determine the lateral preferential flow paths that are active during the steady-state flow conditions and their interaction with the surrounding soil matrix. We also calculate the velocities of the flow through each cross-section of the hillslope and compare them to hillslope scale applied tracer measurements. Finally, we investigate the relationship between the contributing area and the characteristics of the preferential flow paths. The experiment revealed that larger contributing areas coincided with highly developed and hydraulically connected preferential flow paths that had flow with little interaction with the surrounding soil matrix. We found evidence of subsurface erosion and deposition of soil and organic material laterally and vertically within the soil. These results are important because they add to the understanding of the runoff generation, solute transport, and slope stability of preferential flow-dominated hillslopes.

Sign in / Sign up

Export Citation Format

Share Document