Investigating a microbial approach to water conservation: Effects of Bacillus subtilis and Surfactin on evaporation dynamics in loam and sandy loam soils

Semi-arid regions faced with increasingly scarce freshwater resources must manage competing demands in the food-energy-water nexus. A possible solution modifies soil hydrologic properties using biosurfactants to reduce evaporation and improve water retention. In this study, two different soil textures representative of agricultural soils in Kansas were treated with a direct application of the biosurfactant, Surfactin, and an indirect application via inoculation of Bacillus subtilis. Evaporation rates of the wetted soils were measured when exposed to artificial sunlight (1000 W/m2) and compared to non-treated control soils. Experimental results indicate that both treatments alter soil moisture dynamics by increasing evaporation rates by when soil moisture is plentiful (i.e., constant rate period) and decreasing evaporation rates by when moisture is scarce (i.e., slower rate period). Furthermore, both treatments significantly reduced the soil moisture content at which the soil transitioned from constant rate to slower rate evaporation. Out of the two treatments, inoculation with B. subtilis generally produced greater changes in evaporation dynamics; for example, the treatment with B. subtilis in sandy loam soils increased constant rate periods of evaporation by 43% and decreased slower rate evaporation by 49%. In comparing the two soil textures, the sandy loam soil exhibited a larger treatment effect than the loam soil. To evaluate the potential significance of the treatment effects, a System Dynamics Model operationalized the evaporation rate results and simulated soil moisture dynamics under typical daily precipitation conditions. The results from this model indicate both treatment methods significantly altered soil moisture dynamics in the sandy loam soils and increased the probability of the soil exhibiting constant rate evaporation relative to the control soils. Overall, these findings suggest that the decrease in soil moisture threshold observed in the experimental setting could increase soil moisture availability by prolonging the constant rate stage of evaporation. As inoculation with B. subtilis in the sandy loam soil had the most pronounced effects in both the experimental and simulated contexts, future work should focus on testing this treatment in field trials with similar soil textures.

Salt Accumulation during Cropping Season in an Arid Irrigation Area with Shallow Water Table Depth: A 10-Year Regional Monitoring

Nowadays, irrigation takes up about 35% of agricultural water consumption worldwide, and irrigation induced secondary soil salinizationsalinization affects the crop production and sustainable development of arid irrigation areas globally. However, the regular pattern of salt accumulation in the root zone during the cropping season and the contributions of its attribute factors are still unclear. Therefore, a 10-year monitoring was conducted in the Hetao Irrigation District to reveal the soil salt accumulation during the cropping season and to relate it to influential factors, including potential crop evapotranspiration, water input (field irrigation + precipitation) and water table depth. It was found that under the climate conditions and water-saving irrigation measures of the investigated 10-year period, (1) the salt accumulated during the cropping season could be effectively leached by autumn irrigation and the root zone soil could remain suitable for crop germination, (2) the cropping season water deficit (potential crop evapotranspiration − field irrigation − precipitation) showed strong correlation with the cropping season salt accumulation, and (3) maintaining the cropping season average water table depth larger than a critical depth (roughly 3 m) might be the most economical way to alleviate salt accumulation. Therefore, it is recommended to balance the salt leaching and the water table depth controlling in the future water-saving irrigation management practices.

Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis

Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11242-021-01644-9.

Calculation of Steady-State Evaporation for an Arbitrary Matric Potential at Bare Ground Surface

Evaporation from soil columns in the presence of a water table is a long lasting subject that has received great attention for many decades. Available analytical studies on the subject often involve an assumption that the potential evaporation rate is much less than the saturated hydraulic conductivity of the soil. In this study, we develop a new semi-analytical method to estimate the evaporation rate for an arbitrary matric potential head at bare soil surface without assuming that the potential evaporation rate is much less than the saturated hydraulic conductivity of the soil. The results show that the evaporation rates calculated by the new solutions fit well with the HYDRUS-1D simulation. The new solutions also can reproduce the results of potential evaporation rate calculated from previous equations under the special condition of an infinite matric potential head at bare soil surface. The developed new solutions expand our present knowledge of evaporation estimation at bare ground surface to more general field conditions.

Journey of water in pine cones

Pine cones fold their scales when it rains to prevent seeds from short-distance dispersal. Given that the scales of pine cones consist of nothing but dead cells, this folding motion is evidently related to structural changes. In this study, the structural characteristics of pine cones are studied on micro-/macro-scale using various imaging instruments. Raindrops fall along the outer scales to the three layers (bract scales, fibers and innermost lignified structure) of inner pine cones. However, not all the layers but only the bract scales get wet and then, most raindrops move to the inner scales. These systems reduce the amount of water used and minimize the time spent on structural changes. The result shows that the pine cones have structural advantages that could influence the efficient motion of pine cones. This study provides new insights to understand the motion of pine cones and would be used to design a novel water transport system.

Impact of heterogeneity on oxygen transfer in a fluctuating capillary fringe

We performed quasi-two-dimensional flow through laboratory experiments to study the effect of a coarse-material inclusion, located in the proximity of the water table, on flow and oxygen transfer in the capillary fringe. The experiments investigate different phases of mass transfer from the unsaturated zone to anoxic groundwater under both steady-state and transient flow conditions, the latter obtained by fluctuating the water table. Monitoring of flow and transport in the different experimental phases was performed by visual inspection of the complex flow field using a dye tracer solution, measurement of oxygen profiles across the capillary fringe, and determination of oxygen fluxes in the effluent of the flow-through chamber. Our results show significant effects of the coarse-material inclusion on oxygen transfer during the different phases of the experiments. At steady state, the oxygen flux across the unsaturated/saturated interface was considerably enhanced due to flow focusing in the fully water-saturated coarse-material inclusion. During drainage, a zone of higher water saturation formed in the fine material overlying the coarse lens. The entrapped oxygen-rich aqueous phase contributed to the total amount of oxygen supplied to the system when the water table was raised back to its initial level. In case of imbibition, pronounced air entrapment occurred in the coarse lens, causing oxygen to partition between the aqueous and gaseous phases. The oxygen mass supplied to the anoxic groundwater following the imbibition event was found to be remarkably higher (approximately seven times) in the heterogeneous system compared with a similar experiment performed in a homogeneous porous medium.

Structure of drying fronts in three-dimensional porous media

Evaporation in a three-dimensional (3D) porous medium, a sand column saturated by water, was studied using synchrotron x-ray tomography. Three-dimensional images of the medium with a resolution of 7 μm were obtained during the evaporation. The entire column was scanned seven times, resulting in nearly 10(4) 2D cross sections and illustrating the spatial distribution of air, liquid, and solid phases at the pore scale. The results were analyzed in order to gain new insights and better understanding of the characteristics of the drying front that was formed when the liquid-filled pores were invaded by air, as well as the structure of the liquid phase as it was dried. The analysis indicates that the liquid phase has a self-similar fractal structure, with its fractal dimension D(f) in all the cross sections being a function of the water content or saturation. In addition, D(f) for the 3D liquid structure, as well as its density correlation function, were computed using the 3D images. A crossover length scale ξ was identified that separates the fractal regime from the compact geometry. For length scales r>ξ, the density correlation function approaches asymptotically the water content of the porous medium. The drying front is shown to be rough and multi-affine, rather than self-affine. Its properties were also computed using the 3D images. The roughness characteristics agree with those for imbibition in porous media, but not with those of fracture surfaces and crack lines.

Spreading profile of evaporative liquid drops in thin porous layer

Spreading of evaporative liquid drops in a thin porous layer has been studied. The entire spreading process can be divided into three distinct phases according to the change of the wetted porous region size. The first phase is characterized by the expansion of the wetted porous region and shrinking of the liquid drop. Contact line pinning is observed in the wetted porous region in the second phase even with the liquid drop totally absorbed into the porous layer. The third phase sees the shrinkage of the wetted porous region until it is not observable. Based on these observations, a model is devised to simulate the spreading of a liquid drop under the studied conditions. Partial differential equations are used to describe the relation between liquid drop volume and other important parameters of a fluid flow, including maximum wetted region diameter achieved, time taken to reach each spreading process phase, and evaporation rate. Calculated results are in good agreement with the experimental data.