A Telemetric, Gravimetric Platform for Real-Time Physiological Phenotyping of Plant–Environment Interactions
, e61280. Publisher's VersionAbstract
Food security for the growing global population is a major concern. The data provided by genomic tools far exceeds the supply of phenotypic data, creating a knowledge gap. To meet the challenge of improving crops to feed the growing global population, this gap must be bridged.Physiological traits are considered key functional traits in the context of responsiveness or sensitivity to environmental conditions. Many recently introduced high-throughput (HTP) phenotyping techniques are based on remote sensing or imaging and are capable of directly measuring morphological traits, but measure physiological parameters mainly indirectly. This paper describes a method for direct physiological phenotyping that has several advantages for the functional phenotyping of plant–environment interactions. It helps users overcome the many challenges encountered in the use of load-cell gravimetric systems and pot experiments. The suggested techniques will enable users to distinguish between soil weight, plant weight and soil water content, providing a method for the continuous and simultaneous measurement of dynamic soil, plant and atmosphere conditions, alongside the measurement of key physiological traits. This method allows researchers to closely mimic field stress scenarios while taking into consideration the environment’s effects on the plants’ physiology. This method also minimizes pot effects, which are one of the major problems in pre-field phenotyping. It includes a feed-back fertigation system that enables a truly randomized experimental design at a field-like plant density. This system detects the soil-water-content limiting threshold (θ) and allows for the translation of data into knowledge through the use of a real-time analytic tool and an online statistical resource. This method for the rapid and direct measurement of the physiological responses of multiple plants to a dynamic environment has great potential for use in screening for beneficial traits associated with responses to abiotic stress, in the context of pre-field breeding and crop improvement.
The influence of surfactant-application method on the effectiveness of water-repellent soil remediation
. Geoderma 2020
. Publisher's VersionAbstract
Soil water repellency (SWR) has a substantial effect on soil–water hydrology: it hinders infiltration, leading to enhanced surface runoff and soil erosion, and causes preferential flow in the soil profile beyond that from the soil's natural heterogeneity. SWR is associated with soil organic matter content, the latter added to the soil by vegetation exudates, litter and residues, forest fires, and replacement of fresh water by treated wastewater for irrigation. Surfactants are surface-active substances composed of organic molecules with hydrophobic tails and hydrophilic heads that can reduce the surface tension (γ) of the aqueous solution, thereby reducing SWR, via adsorption to soil particles. Surfactants are commonly used to remediate water-repellent soils. We investigated the role of two surfactant-application methods on the efficacy of SWR remediation. Aqueous solutions of two commercial surfactants had a substantial effect on parameters used to characterize the persistence and severity of SWR. However, the efficacy of these surfactants in remediating sandy soils rendered water-repellent by irrigation with treated effluent was substantially affected by their application method. Whereas application of aqueous surfactant solution to the surface of water-repellent soil, the commonly used remediation method, formed finger-like plumes similar to those obtained for water application, bulbous-like plumes were formed when the soil was premixed with the aqueous surfactant solution prior to water application. These differences were attributed to the significant role of the rate-limited surfactant adsorption to the soil particles. © 2019 Elsevier B.V.
Dynamic physiological phenotyping of drought-stressed pepper plants treated with “productivity-enhancing” and “survivability-enhancing” biostimulants
. Frontiers in Plant Science 2019
. Publisher's VersionAbstract
The improvement of crop productivity under abiotic stress is one of the biggest challenges faced by the agricultural scientific community. Despite extensive research, the research-to-commercial transfer rate of abiotic stress-resistant crops remains very low. This is mainly due to the complexity of genotype × environment interactions and in particular, the ability to quantify the dynamic plant physiological response profile to a dynamic environment. Most existing phenotyping facilities collect information using robotics and automated image acquisition and analysis. However, their ability to directly measure the physiological properties of the whole plant is limited. We demonstrate a high-throughput functional phenotyping system (HFPS) that enables comparing plants’ dynamic responses to different ambient conditions in dynamic environments due to its direct and simultaneous measurement of yield-related physiological traits of plants under several treatments. The system is designed as one-to-one (1:1) plant–[sensors+controller] units, i.e., each individual plant has its own personalized sensor, controller and irrigation valves that enable (i) monitoring water-relation kinetics of each plant–environment response throughout the plant’s life cycle with high spatiotemporal resolution, (ii) a truly randomized experimental design due to multiple independent treatment scenarios for every plant, and (iii) reduction of artificial ambient perturbations due to the immobility of the plants or other objects. In addition, we propose two new resilience-quantifying-related traits that can also be phenotyped using the HFPS: transpiration recovery rate and night water reabsorption. We use the HFPS to screen the effects of two commercial biostimulants (a seaweed extract –ICL-SW, and a metabolite formula – ICL-NewFo1) on Capsicum annuum under different irrigation regimes. Biostimulants are considered an alternative approach to improving crop productivity. However, their complex mode of action necessitates cost-effective pre-field phenotyping. The combination of two types of treatment (biostimulants and drought) enabled us to evaluate the precision and resolution of the system in investigating the effect of biostimulants on drought tolerance. We analyze and discuss plant behavior at different stages, and assess the penalty and trade-off between productivity and resilience. In this test case, we suggest a protocol for the screening of biostimulants’ physiological mechanisms of action. © 2019 Dalal, Bourstein, Haish, Shenhar, Wallach and Moshelion.
Structure and hydraulic properties in soils under long-term irrigation with treated wastewater
. Geoderma 2019
, 90-98. Publisher's VersionAbstract
Secondary treated wastewater, a commonly used water resource in agriculture in (semi-)arid areas, often contains salts, sodium, and organic matter which may affect soil structure and hydraulic properties. The main objective of this study was to jointly analyse the effects of long-term irrigation with treated wastewater on physicochemical soil characteristics, soil structure, and soil water dynamics in undisturbed soils. X-ray microtomography was used to determine changes in macro-porosity (> 19 μm), pore size distribution, and pore connectivity of a sandy clay loam and a loamy sand. Differences in the pore network among soils irrigated with treated wastewater, fresh water that replaced treated wastewater, and non-irrigated control plots could be explained by changes in textural composition, soil physicochemical parameters, and hydraulic properties. In this study we showed that irrigation led to the development of a connected macro-pore network, independent of the studied water quality. The leaching of silt and clay particles in the sandy soil due to treated wastewater irrigation resulted in an increase of pores < 130 μm. While this change in texture reduced water retention, the unsaturated hydraulic conductivity was diminished by physicochemical alteration, i.e. induced water repellency and clay mineral swelling. Overall, the fine textured sandy clay loam was much more resistant to soil alteration by treated wastewater irrigation than the loamy sand. © 2018 Elsevier B.V.
Chapter 3 - Physical Characteristics of Soilless Media
. In Soilless Culture (Second Edition)
; Soilless Culture (Second Edition); Elsevier: Boston, 2019; pp. 33 - 112. Publisher's VersionAbstract
Similar to soils, soilless media are composed of three phases: solid, aqueous, and gaseous. The substantial differences between the former and the latter are related to the solid phase composition and associated texture and structure. While the solid phase of soils is mostly minerals, it is a mixture of minerals and organic matter for soilless culture, resulting in substantial differences in their physical and hydraulic properties. The physical and hydraulic properties of soil have been intensively studied, both experimentally and theoretically, the application and implementation of these knowledge and methods to soilless substrates is limited. This chapter aims to close the gaps in physics and hydraulic characterization and application between mineral soils and soilless media by applying terminology, methods, and approaches that are commonly used for the former to the latter. The current version of this Chapter extents the topic of particle wettability and its effect on the physical properties of soilless media. Although research on sales-media wettability and its effect on water retention and flow in the media has only recently initiated, its relevance is essential owing to the high organic matter content of these media and the relationship between the two.
Structure and hydraulic properties in soils under long-term irrigation with treated wastewater
. Geoderma 2019
, 90 - 98. Publisher's VersionAbstract
Secondary treated wastewater, a commonly used water resource in agriculture in (semi-)arid areas, often contains salts, sodium, and organic matter which may affect soil structure and hydraulic properties. The main objective of this study was to jointly analyse the effects of long-term irrigation with treated wastewater on physicochemical soil characteristics, soil structure, and soil water dynamics in undisturbed soils. X-ray microtomography was used to determine changes in macro-porosity (> 19 μm), pore size distribution, and pore connectivity of a sandy clay loam and a loamy sand. Differences in the pore network among soils irrigated with treated wastewater, fresh water that replaced treated wastewater, and non-irrigated control plots could be explained by changes in textural composition, soil physicochemical parameters, and hydraulic properties. In this study we showed that irrigation led to the development of a connected macro-pore network, independent of the studied water quality. The leaching of silt and clay particles in the sandy soil due to treated wastewater irrigation resulted in an increase of pores < 130 μm. While this change in texture reduced water retention, the unsaturated hydraulic conductivity was diminished by physicochemical alteration, i.e. induced water repellency and clay mineral swelling. Overall, the fine textured sandy clay loam was much more resistant to soil alteration by treated wastewater irrigation than the loamy sand.
Spatial variation of soil water repellency in a commercial orchard irrigated with treated wastewater
. Geoderma 2019
, 214 - 224. Publisher's VersionAbstract
The recognition of treated wastewater (TWW) as an alternative water resource is expanding in areas with a shortage of freshwater (FW). While many studies have been devoted to the effects of long-term irrigation with TWW on soil wettability and spatial flow variations in the soil profile, much less attention has been given to the spatial distribution of soil water repellency in the soil surface layer. This is the objective of the current study. Undisturbed soil samples (5 cm thick) were taken at 15-cm intervals parallel to a drip lateral in two adjacent plots of a commercial citrus orchard in central Israel. Each soil sample was sectioned into five consecutive 1-cm layers for which soil water repellency was determined by water drop penetration time method, and soil organic matter by loss-on-ignition method. Geostatistics and multivariate empirical mode decomposition were used to investigate the overall and scale-specific spatial variation of soil water repellency and its dependence on dripper intervals along the lateral. A high degree of soil water repellency with strong spatial variation was found in the surface soil after 4–6 years of TWW irrigation. Weak to moderate spatial dependence of soil water repellency with maximum autocorrelation distance of around 30 cm was discovered by geostatistical analysis. The spatial distribution of soil water repellency was considered to be greatly affected by the location of the drippers, being higher between adjacent drippers and lower underneath them. This soil water repellency distribution is presumed to result from ongoing lateral displacement of the amphiphilic substances in the TWW toward the outer edge of the wetted plume periphery. Multivariate empirical mode decomposition of the overall spatial variation of soil water repellency yielded three scale-specific variations with corresponding characteristic scales of 30 cm, 110 cm and 200 cm. Most of the soil water repellency variation was separated into the 30 cm and 110 cm spatial scales, which were correlated to processes related to the drippers and trees. Replacing TWW with FW for the reclamation of water-repellent soils partially alleviated the intensity of TWW irrigation-induced soil water repellency. Moreover, an inconsistency between the hot spots of water-repellency development between adjacent drippers and the areas that are effectively ameliorated by FW irrigation below the drippers could be developed and affect the spatial distribution of flow pattern in an a priori unpredictable way.
A High-Throughput Physiological Functional Phenotyping System for Time- and Cost-Effective Screening of Potential Biostimulants
. bioRxiv 2019
, 525592. Publisher's VersionAbstract
The improvement of crop productivity under abiotic stress is one of the biggest challenges faced by the agricultural scientific community. Despite extensive research, the research-to-commercial transfer rate of abiotic stress-resistant crops remains very low. This is mainly due to the complexity of genotype◻×◻environment interactions and in particular, the ability to quantify the dynamic plant physiological response profile to a dynamic environment.Most existing phenotyping facilities collect information using robotics and automated image acquisition and analysis. However, their ability to directly measure the physiological properties of the whole plant is limited. We demonstrate a high-throughput functional phenotyping system (HFPS) that enables comparing plants’ dynamic responses to different ambient conditions in dynamic environments due to its direct and simultaneous measurement of yield-related physiological traits of plants under several treatments. The system is designed as one-to-one (1:1) plant–[sensors+controller] units, i.e., each individual plant has its own personalized sensor, controller and irrigation valves that enable (i) monitoring water-relation kinetics of each plant–environment response throughout the plant’s life cycle with high spatiotemporal resolution, (ii) a truly randomized experimental design due to multiple independent treatment scenarios for every plant, and (iii) reduction of artificial ambient perturbations due to the immobility of the plants or other objects. In addition, we propose two new resilience-quantifying-related traits that can also be phenotyped using the HFPS: transpiration recovery rate and night water reabsorption.We use the HFPS to screen the effects of two commercial biostimulants (a seaweed extract—ICL-SW, and a metabolite formula—ICL-NewFo1) on Capsicum annuum under different irrigation regimes. Biostimulants are considered an alternative approach to improving crop productivity. However, their complex mode of action necessitates cost-effective pre-field phenotyping. The combination of two types of treatment (biostimulants and drought) enabled us to evaluate the precision and resolution of the system in investigating the effect of biostimulants on drought tolerance. We analyze and discuss plant behavior at different stages, and assess the penalty and trade-off between productivity and survivability. In this test case, we suggest a protocol for the screening of biostimulants’ physiological mechanisms of action.
ERT and salinity – A method to determine whether ERT-detected preferential pathways in brackish water-irrigated soils are water-induced or an artifact of salinity
. Journal of Hydrology 2019
, 35 - 45. Publisher's VersionAbstract
The identification of preferential flow generation and the tracking of their pathways in agricultural soils are challenging tasks, in particular while keeping the soil profile undisturbed. The noninvasive electrical resistivity tomography (ERT) method seems, a priori, to be an ideal tool for these purposes, but the simultaneous dependence of soil electrical resistivity (ER) on soil moisture content and salt concentration poses a problem when high-salinity water is involved. The dependence of these variables is expressed, for example, by Archie’s law. Irrigation with brackish or treated wastewater are two such complicated cases. The latter has been found to render soils water-repellent with inherent flow in preferential pathways. A method that resolves this complication is developed and applied for in-situ-measured ERT data. This method combines a model that simulates the preferential ER(t) as a series of interconnected well-mixed units (WMUs) with frequent ERT scans during soil wetting and subsequent drying. Following validation by comparison of its output to simulations using the Richards equation for flow and convection-dispersion equation for transport, the WMU model was used to simulate scenarios of simultaneous variations in moisture content, salinity, and ER(t) in a well-mixed soil volume (unit). The characteristic patterns acquired by these simulations were compared with in-situ measured ER(z,t) to determine whether the observed preferential ER pathways coincide with preferential flow or salinity pathways, or both. From the eight sections in the profile where this analysis was made for, five exhibit preferential flow, characterized by fast wetting front propagation that includes by-passing of a certain soil volume, while the wetting front velocity in the other three was moderate. We, therefore, concluded, based on the water content/salinity effect on ER, that the measured preferential ER pathways coincide with preferential flow pathways. Salinity has, for the current case, a minor, if any, effect on the ER spatial variation. Therefore, ERT is an advantageous mean of mapping of preferential flow pathways in the soil profile, and the method developed herein can be used for studying the water distribution in the soil profile. Particular irrigation design and management can be then used to ease the negative effects of preferential flow pathways on soil water availability to plant roots and the enhanced leaching of agrochemical below the active root zone.
Modeling gravity-driven fingering by using the moving boundary approach
. In EGU General Assembly Conference Abstracts
; EGU General Assembly Conference Abstracts; 2018; Vol. 20, pp. 4729.
Quantitative analysis of wetting front instabilities in soil caused by treated waste water irrigation
. Geoderma 2018
, 132 - 141. Publisher's VersionAbstract
Irrigation with treated waste water (TWW) is a common practice in agriculture, mainly in arid and semiarid areas as it provides a sustainable water resource available at all-season in general and at freshwater shortage in particular. However, TWW still contains abundant organic material which is known to decrease soil wettability, which in turn may promote flow instabilities that lead to the formation of preferential flow paths. We investigate the impact of long-term TWW irrigation on water wettability and infiltration into undisturbed soil cores from two commercially used orchards in Israel. Changes of water content during infiltration were quantitatively analysed by X-ray radiography. One orchard (sandy clay loam) had been irrigated with TWW for more than thirty years. In the other orchard (loamy sand) irrigation had been changed from freshwater to TWW in 2008 and switched back in some experimental plots to freshwater in 2012. Undisturbed soil cores were taken at the end of the dry and the rainy season to investigate the seasonal effect on water repellency and on infiltration dynamics in the laboratory. The irrigation experiments were done on field moist samples. A test series with different initial water contents was run to detect the influence on water movement at different wettabilities. In this study we show that the infiltration front stability is dependent on the history of waste water irrigation at the respective site and on the initial water content.
Risk-management strategies and transpiration rates of wild barley in uncertain environments
. Physiologia Plantarum 2018
, 412 - 428. Publisher's VersionAbstract
Regulation of the rate of transpiration is an important part of plants' adaptation to uncertain environments. Stomatal closure is the most common response to severe drought. By closing their stomata, plants reduce transpiration to better their odds of survival under dry conditions. Under mild to moderate drought conditions, there are several possible transpiration patterns that balance the risk of lost productivity with the risk of water loss. Here, we hypothesize that plant ecotypes that have evolved in environments characterized by unstable patterns of precipitation will display a wider range of patterns of transpiration regulation along with other quantitative physiological traits (QPTs), compared to ecotypes from less variable environments. We examined five accessions of wild barley (Hordeum vulgare ssp. spontaneum) from different locations in Israel (the B1K collection) with annual rainfall levels ranging from 100 to 900?mm, along with one domesticated line (cv. Morex). We measured several QPTs and morphological traits of these accessions under well-irrigated conditions, under drought stress and during recovery from drought. Our results revealed a correlation between precipitation-certainty conditions and QPT plasticity. Specifically, accessions from stable environments (very wet or very dry locations) were found to take greater risks in their water-balance regulation than accessions from areas in which rainfall is less predictable. Notably, less risk-taking genotypes recovered more quickly than more risk-taking ones once irrigation was resumed. We discuss the relationships between environment, polymorphism, physiological plasticity and fitness, and suggest a general risk-taking model in which transpiration-rate plasticity is negatively correlated with population polymorphism.
The moving-boundary approach for modeling gravity-driven stable and unstable flow in soils
. Water Resources Research 2017
, 344-360. Publisher's VersionAbstract
Abstract The Richards equation is unsuccessful at describing gravity-driven unstable flow with nonmonotonic water content distribution. This shortcoming is resolved in the current study by introducing the moving-boundary approach. Following this approach, the flow domain is divided into two subdomains with a sharp change in fluid saturation between them (moving boundary). The upper subdomain consists of water and air, whose relationship varies with space and time following the imposed boundary condition at the soil surface calculated by the Richards equation. The lower subdomain consists of an initially dry soil that remains constant. The location of the boundary between the two subdomains is part of the solution, rendering the problem nonlinear. The moving boundary solution was used after verification to demonstrate the effect of contact angle, soil characteristic curves and incoming flux on the dynamic water-entry pressure of the soil, which depends on the soil's wettability, incoming flux at the soil surface and the wetting front's propagation rate. Lower soil wettability hinders spontaneous invasion of the dry pores and, together with a higher input flux, induces water accumulation behind the wetting front (saturation overshoot). The wetting front starts to propagate once the pressure building up behind it exceeds the dynamic water-entry pressure. To conclude, the physically based novel moving-boundary approach for solving stable and gravity-driven unstable flow in soils was developed and verified. It supports the conjecture that saturation overshoot is a prerequisite for gravity-driven fingering.
High-throughput physiological phenotyping and screening system for the characterization of plant-environment interactions
. Plant J 2017
We present a simple and effective high-throughput experimental platform for simultaneous and continuous monitoring of water relations in the soil-plant-atmosphere continuum of numerous plants under dynamic environmental conditions. This system provides a simultaneously measured, detailed physiological response profile for each plant in the array, over time periods ranging from a few minutes to the entire growing season, under normal, stress and recovery conditions and at any phenological stage. Three probes for each pot in the array and a specially designed algorithm enable detailed water-relations characterization of whole-plant transpiration, biomass gain, stomatal conductance and root flux. They also enable quantitative calculation of the whole plant water-use efficiency and relative water content at high resolution under dynamic soil and atmospheric conditions. The system has no moving parts and can fit into many growing environments. A screening of 65 introgression lines of a wild tomato species (Solanum pennellii) crossed with cultivated tomato (S. lycopersicum), using our system and conventional gas-exchange tools, confirmed the accuracy of the system as well as its diagnostic capabilities. The use of this high-throughput diagnostic screening method is discussed in light of the gaps in our understanding of the genetic regulation of whole-plant performance, particularly under abiotic stress.
Induced heterogeneity of soil water content and chemical properties by treated wastewater irrigation and its reclamation by freshwater irrigation
. Water Resources Research 2017
, 4756-4774. Publisher's VersionAbstract
Abstract The recognition of treated wastewater (TWW) as an alternative water resource is expanding in areas with a shortage of freshwater (FW) resources. Today, most orchards in Israel are irrigated with TWW. While the benefits of using TWW for irrigation are apparent, evidence of its negative effects on soil, trees, and yield is accumulating. This study, performed in a commercial TWW-irrigated citrus orchard in central Israel, examined the effects of (1) soil-wettability decrease due to prolonged TWW irrigation on the spatial and temporal distribution of water content and associated chemical properties in the root zone; (2) the conversion of irrigation in half of the TWW-irrigated research plot to FW (2012) for soil reclamation. Electrical resistivity tomography surveys in the substantially water repellent soils revealed that water flow is occurring along preferential flow paths in both plots, leaving behind a considerably nonuniform water-content distribution. This was despite the gradual relief in soil water repellency measured in the FW plots. Four soil-sampling campaigns (spring and fall, 2014–2016), performed in 0–20 and 20–40 cm layers of the research plot, revealed bimodal gravimetrically measured water-content distribution. The preferential flow led to uneven chemical-property distribution, with substantially high concentrations in the dry spots, and lower concentrations in the wet spots along the preferential flow paths. The average salt and nutrient concentrations, which were initially high in both plots, gradually dispersed with time, as concentrations in the FW plots decreased. Nevertheless, the efficiency of reclaiming TWW soil by FW irrigation appears low.
A combination of stomata deregulation and a distinctive modulation of amino acid metabolism are associated with enhanced tolerance of wheat varieties to transient drought
. Metabolomics 2017
, 138. Publisher's VersionAbstract
Mediterranean winter crops are commonly and increasingly exposed to irregular rainfall and high temperatures, which lead to transient drought events of different degrees, adversely affecting growth and yield. Hence, exploring the diverse degrees of tolerance to drought existing in the crop and the molecular strategies behind it is pivotal for the development of ad hoc breeding programs.
Modeling the movement of water and solute through preferential flow paths and fractures
. In The Handbook of Groundwater Engineering: Third Edition
; The Handbook of Groundwater Engineering: Third Edition; 2016; pp. 509 - 530. Publisher's Version
Formation of soil–water repellency in olive orchards and its influence on infiltration pattern
, 1 - 11. Publisher's VersionAbstract
Water-repellent (hydrophobic) soils do not wet instantaneously, but only after some time (a few seconds to hours) of soil-particle contact with water. Some plant species can render soils hydrophobic but in this respect, olive trees have scarcely been examined. Measurements of water repellency in olive orchards of different ages in different locations in Israel using the water drop penetration time (WDPT) test have shown that soils tend to become hydrophobic, regardless of texture and structure. A comprehensive study was then performed for an irrigated young and mature olive grove and nearby uncultivated bare soil in the southern part of Israel. The study included intensive WDPT measurements, initial (repellency intensity) and rate of decrease (repellency persistence) for sessile drops placed on the soil surface, cumulative infiltration using tension disc infiltrometer, and monitoring flow in a transparent flow chamber packed with soils from the different plots. The soil from the mature olive plot was noticeably more water repellent than the young plot's soil, and both differed from the uncultivated soil that was fully wettable. The contact angle of a drop placed on the surface of a single layer of soil particles decreased exponentially with time, with a lower decay rate for the mature orchard soil. Cumulative infiltration had a convex pattern for wettable soils and a concave pattern for water-repellent ones. The difference in infiltration pattern was attributed to water/pressure buildup behind the wetting front as a result of the dynamic contact-angle-induced pore resistivity to wetting. The supplemental pressure, also known as dynamic water-entry pressure, increases the infiltration rate beyond that obtained by the capillary pressure per se. The significant correlation between soil sorptivity and the asymptotic infiltration rate, both calculated from the cumulative infiltration curves, and the WDPT substantiates the dependence of pressure overshoot and the rate at which the contact angle decreases prior to pore wetting. The considerable differences in plume shape, size, and internal saturation distribution between the wettable and water-repellent soils, indicating unstable flow in the latter, were also explained by the wettability-dependent water-entry pressure. The outcome of this study indirectly supports the findings that higher surface runoff and erosion are associated with no-till farming in olive orchards, due to the combination of no-till cropping and the near-surface accumulation of hydrophobic organic carbon compounds.
Kinetics of gravity-driven slug flow in partially wettable capillaries of varying cross section
. Water Resources Research 2016
, 8472-8486. Publisher's VersionAbstract
Abstract A mathematical model for slug (finite liquid volume) motion in not-fully-wettable capillary tubes with sinusoidally varying cross-sectional areas was developed. The model, based on the Navier-Stokes equation, accounts for the full viscous terms due to nonuniform geometry, the inertial term, the slug's front and rear meniscus hysteresis effect, and dependence of contact angle on flow velocity (dynamic contact angle). The model includes a velocity-dependent film that is left behind the advancing slug, reducing its mass. The model was successfully verified experimentally by recording slug movement in uniform and sinusoidal capillary tubes with a gray-scale high-speed camera. Simulation showed that tube nonuniformity has a substantial effect on slug flow pattern: in a uniform tube it is monotonic and depends mainly on the slug's momentary mass/length; an undulating tube radius results in nonmonotonic flow characteristics. The static nonzero contact angle varies locally in nonuniform tubes owing to the additional effect of wall slope. Moreover, the nonuniform cross-sectional area induces slug acceleration, deceleration, blockage, and metastable-equilibrium locations. Increasing contact angle further amplifies the geometry effect on slug propagation. The developed model provides a modified means of emulating slug flow in differently wettable porous media for intermittent inlet water supply (e.g., raindrops on the soil surface).