Soil Water Movement Research Articles

Influence of Annual Ryegrass (Lolium multiflorum) as Cover Crop on Soil Water Dynamics in Fragipan Soils of Southern Illinois, USA

Fragipans are dense subsurface soil layers that severely restrict root penetration and water movement. The presence of shallow fragipan horizons limits row crop production. We hypothesized that the roots of cover crop might improve soil physiochemical properties and biological activity, facilitating drainage and increasing effective soil depth for greater long-term soil water storage. To evaluate annual ryegrass as one component of a cover crop (CC) mix for promoting the characteristics and distribution of soil water, on-farm studies were conducted at Marion and Springerton in southern Illinois, USA. Soil samples were collected at 15 cm increments to 60 cm (Marion) and 90 cm (Springerton) depths during the fall of 2022. Both sites had low total soil carbon and nitrogen contents and acidic soil pH (≤6.4). A soil water retention curve was fitted using the van Genuchten equation. At Springerton, the CC treatment increased saturated (thetaS) and residual (thetaR) soil water contents above those of the no cover crop (NCC) at the 60–75 cm and 75–90 cm depths. Changes in volumetric soil water content were measured using a multi-depth soil water sensor for the Springerton site during late July to early August of the soybean growing phase of 2022; NCC had higher soil water than CC within the 0–15 cm depth, but CC had higher soil water than NCC at the 30–45 cm depth. These findings indicate that cover crop mix has the potential to improve soil water movement for soils with restrictive subsoil horizon, possibly through reducing the soil hydraulic gradient between the surface and restrictive subsurface soil layers.

A novel conceptual model coupling crop growth and soil water-heat-salt processes in arid area

In arid farmland, the soil water, heat, and salt movement tightly couple with crop growth. Accurately expressing the coupling relationships is essential to understanding the complex agricultural hydrological system and enhancing agricultural water productivity. A novel conceptual model coupling crop growth, soil water-heat-salt movement, and energy balance in arid area with shallow groundwater was developed which has potential for use in large heterogeneous irrigation districts. The model can calculate the crop-soil system through the effect of crop canopy and leaf area index (LAI) on potential evapotranspiration and the stress of soil moisture, temperature, and salinity status to root water uptake. The energy balance equation was used for calculating the upper boundary condition for soil temperature. The model was calibrated and validated using four years of field monitoring data for two crops located in a typical arid agricultural area in China. The model performed well in simulated field hydrology and crop growth processes. The model can capture soil water-heat-salt dynamics and model field evapotranspiration with RMSE below 1.14 mm/day, and crop transpiration with R2 of above 0.61. Furthermore, the model can accurately describe crop growth processes with R2 of 0.95 and RMSE of 0.49 for LAI, and R2 of 0.99, RMSE of 12.95 cm for crop height. With few parameters, our work supplies an alternative method for quantifying agricultural hydrological processes. In particular, the model has potential to simulate the regional crop growth and soil water-heat-salt processes for heterogeneous agricultural area. This method is helpful in determining the scientific irrigation management schedule.

Do Cover Crops Influence In Situ Soil Water Potential After Termination?

Soil water movement is energy-dependent and is influenced by various management practices. It can be understood by measuring the soil water potential (SWP); however, the influence of cover crops (CCs) on SWP is not currently well understood. The objective of this study was to assess the effects of CCs on SWP before and after termination in order to understand their effects on soil water availability for the subsequent cash crop. The experimental design was a completely randomized design with two levels of CCs (CCs vs. no cover crop [NC]) with three replicates. The SWP sensors were buried at 0–10, 10–20, and 20–30 cm depths before CCs were planted. Additionally, soil samples were collected at the aforementioned depths just before CC termination for soil organic carbon (SOC), bulk density (BD), and saturated hydraulic conductivity (Ksat) analysis. Results showed that CCs increased SOC and Ksat, and significantly lowered BD compared with NC management. Before termination, CC plots had significantly lower SWP values compared with NC management, suggesting that the transpirational needs of the CCs can lead to lower water content. After termination, CC management also resulted in lower SWP, suggesting that CCs can increase water-use efficiency by improving soil health parameters. However, effective planning is required for CC implementation, especially in semiarid and arid regions.

Simulating coefficient of soil moisture content uniformity of sprinkler irrigation systems using a COMSOL-3D model

Water distribution uniformity is important for assessing the hydraulic performance of sprinkler nozzles and designing sprinkler irrigation systems. However, current studies rarely consider the redistribution abilities of sprinkler water in the soil and judge whether a sprinkler system is qualified based only on the coefficient of uniformity on the ground (CUg), which leads to a serious underestimation of the actual irrigation effect of sprinkler irrigation projects. Hence, this study establishes and validates a COMSOL-3D model to simulate the movement of soil water during sprinkler irrigation. Moreover, the effects of the soil hydraulic and irrigation parameters on the coefficient of soil moisture content uniformity (CUs) were explored. The results showed that 1) the Nash-Sutcliffe efficiency coefficient (NSE) of the vertical wetting front (VWF) and soil moisture content (SMC) in the two sets of validation experiments was 0.965 and 0.862 and 0.867 and 0.900, respectively. 2) Owing to the differences in soil saturated hydraulic conductivity, after the end of irrigation, the VWF transport rates in vertical direction differed among the three textured soils, with loam exhibiting the highest rate, followed clay loam and silty clay. However, their CUs increased with increasing water transport time, indicating that the water distribution within the sprinkler-wetted zone became increasingly homogeneous as the sprinkler irrigation water was transported and diffused from the top to bottom in the soil. 3) Dimensional analyses showed that CUs were strongly influenced by the initial soil moisture content, irrigation time, and water transport time. 4) In Case 4, when CUg was reduced from 81.9 % to 60 %, CUs only decreased from 83.2 % to 81.4 % after 48 h of irrigation, indicating that the sprinkler irrigation water had good redistributive ability. This study provides a reference for rationally reducing the CUs of sprinkler irrigation projects and maximising project investments and planting interest.

Simulation algorithm of greenhouse soil water movement based on cellular automata

Soil moisture movement reveals the hydrological environment and irrigation characteristics of plants, which is very important and basic hydrological problem. A soil moisture movement algorithm based on cellular automata suitable for greenhouse was proposed to simulate the lateral flow and vertical penetration of soil water in different soil layers. Artemisia annua was used as an experimental plant, and two numerical simulations were set up to determine the minimum water replenishment and the degree of upward soil modification. To verify the performance of the model, corn was used as an experimental plant for application simulation. The numerical simulation results show that the effect of only water for plant cells is superior to soil cells, and the disturbance behavior of the underlying soil layer should be minimized during the cultivation process. The average MAE/RMSE of all soil layers is 0.86/1.06 and in deep layer (50 cm and 60 cm) is 0.53/0.64. It shows the model has a certain prediction and simulation ability, especially in deep soil layers. The proposed algorithm can simply calculate soil flow, set the minimum water replenishment, and evaluate the water replenishment efficiency, which can provide a theoretical reference for the water replenishment and soil replacement scheme.

Improved partitioning between matrix and macropore flow: Novel bimodal lognormal functions for water retention and hydraulic conductivity in pumice and non-pumice soils

Dual-porosity models have been shown to improve models of soil–water movement and enhance the water balance of structured soils. In this study we introduce novel, continuous, closed-form, bimodal, lognormal functions for soil–water retention, θ(ψ), and unsaturated hydraulic conductivity, K(ψ), enhancing traditional models and significantly improving predictions, particularly for pumice soils. These functions incorporate the thresholds for (a) water pressure, (b) soil water content, and (c) unsaturated conductivity, which accurately differentiates macropores from matrix pores. Validation using 313 observation points from laboratory data shows an increase in the Nash–Sutcliffe efficiency coefficient of θ(ψ) from 0.94 to 0.97, and for K(ψ) from 0.83 to 0.95. The model also requires six constant, semi-empirical tortuosity parameters and provides a physically constrained approach for the hydraulic parameters that reduces non-uniqueness risks. The derived functions yield improved predictions and enable the computation of macropore soil water content and flow contributions, with potential applications for the advancement of preferential flow modelling.

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

Interpretation of Soil Characteristics and Preferential Water Flow in Different Forest Covers of Karst Areas of China

Soil hydrology seriously affects the prevention of desertification in karst areas. However, water infiltration in the different soil layers of secondary forests and artificial forests in karst areas remains uncertain. This lack of clarity is also the factor that constrains local vegetation restoration. Therefore, monitoring and simulating the priority transport of soil moisture will help us understand the shallow soil moisture transport patterns after artificial vegetation restoration in the local area, providing a reference for more scientific restoration of the ecological environment and enhancement of carbon storage in karst areas. The integration of soil physical property assessments, computed tomography (CT) scanning, dye tracing studies, and HYDRUS-2D modeling was utilized to evaluate and contrast the attributes of soil macropores and the phenomenon of preferential flow across various forestland categories. This approach allowed for a comprehensive analysis of how the soil structure and water movement are influenced by different forest ecosystems and infiltration head simulations (5 mm, 15 mm, 35 mm, and 55 mm) to elucidate the dynamics of water movement across diverse soil types within karst regions, to identify the causes of water leakage due to preferential flow in secondary forests, and to understand the mechanisms of water conservation and reduction in artificial forests adopting a multifaceted approach. This study demonstrated that (1) the soil hydrological capacity of a plantation forest was 20% higher than a natural forest, which may be promoted by the clay content and distribution. (2) Afforestation-enhanced soils in karst regions demonstrate a significant capacity to mitigate the loss of clay particles during episodes of preferential flow and then improve the soil erosion resistance by about 5 times, which can effectively control desertification in karst area. (3) The uniform distribution of macropores in plantation forest soil was conducive to prevent water leakage more effectively than the secondary forest but was incapable of hindering the occurrence of preferential flow. The secondary forest had a very developed preferential flow phenomenon, and soil clay deposition occurred with an increase in depth. (4) Moreover, the results for preferential flow showed that the matrix flow depth did not increase with the increase in water quantity. Short-term and high-intensity heavy rainfall events facilitated the occurrence of preferential flow. Infiltration along the horizontal and vertical directions occurred simultaneously. These results could facilitate a further understanding of the contribution of the plantation to soil amelioration and the prevention of desertification in karst areas, and provide some suggestions for the sustainable development of forestry in karst areas where plantation restoration is an important ingredient.

Approximate analytical solution of a soil water movement equation under different ponding radii on the basis of numerical simulation

AbstractThis study addresses the problem of 2D soil water movement under ponding radii of 1, 2, and 3 cm. The soil water movement characteristics (shape parameters of the water content profile, ratio of horizontal wetting front to vertical wetting front, relationship between infiltration time and horizontal wetting front, and relationship between infiltration time and cumulative infiltration) under the above three kinds of water ponding radius were analyzed. On the basis of the assumption that the soil wetting body is a semi‐ellipse and the analytical solution of the 1D soil water movement equation at any angle, the approximate analytical solution of the 2D soil water movement equation under ponding conditions is optimized. The function relationships between infiltration time, wetting front, and cumulative infiltration are established. We applied the numerical data simulated by HYDRUS‐3D to validate the parameters in proposed analytical solutions and evaluated the relationships between the wetting front and hydraulic parameters. The results indicate that as the water ponding radius increases, the wetting body and 2D water content distribution becomes larger. When the water ponding radius was 2 cm, the numerical and analytical solution of 1D soil water distribution showed the best comparison results, and the model error was the smallest. The ratio of wetting fronts was linearly increased with the increase of air‐entry suction with R2 = 0.9969.

Differences in soil water movement between the dip and anti-dip slopes of a karst trough valley

Exploring temporal and spatial variations in soil water processes is crucial for identifying the ecological water supply capacity of soil and understanding hydrologic cycles within karst catchments. However, the existing body of knowledge related to soil water dynamics across dip and anti-dip slopes in karst landscapes remains limited. To address this gap, this study investigated soil water dynamics in a karst trough valley (KTV) in southwest China, using stable isotopes (δD and δ18O) from eight sampling sites and precipitation data from April 2017 to November 2018. Our findings suggest that seasonal variations in soil water isotopes are primarily influenced by climatic factors, with soil water lagging behind rainwater isotopes by one to four months due to differences in the residence time of precipitation in different soil layers. Spatially, our results reveal the coexistence of diffuse and preferential flows on both dip and anti-dip slopes. Specifically, on the dip slope, preferential flow dominates above 20 cm, while diffuse flow prevails below this depth. In contrast, on the anti-dip slope, preferential flow is prominent above 40 cm depth, while diffuse flow dominates beneath this threshold. These spatial patterns not only impact the ecological water and nutrient supply capacity of soil but they also contribute to distinct patterns in vegetative growth on both slopes. In turn, vegetation growth facilitates preferential flow pathways, establishing a positive feedback loop between soil water movement and vegetative cover. These findings may lead to a deeper understanding of soil water processes on the dip and anti-dip slopes of KTVs and provide crucial insights into the relationship between soil water and vegetative growth in similar karst landforms.

ESTIMATION OF SOIL WETTING PATTERN USING DIFFERENT IRRIGATION REGIMES IN CONJUCTION WITH SUPER ABSORBENT POLYMER

The design and management of irrigation systems require soil water movement and distribution knowledge. Finding the actual ability of drip irrigation to supply sufficient water uniformly and timely requires optimizing the operational management parameters. The objective of the study was to estimate the water distribution pattern with drip irrigation under different irrigation frequencies, durations and Belsap super absorbent polymer. A factorial experimental design layout of drip irrigation system was set up in Noiwet, Mogotio Sub County, Baringo County. The same quantity of water with three irrigation frequencies were applied once daily, after two days and three days and three durations which were full, half and a third the irrigation time. The total number of treatments were 18 with three replications making a total of 54 experiments. The wetted depths and widths of the soil were measured directly in the field using a ruler immediately after irrigation. ANOVA statistical analysis was carried out to determine the significance of the treatments on the soil wetting pattern of the clay soil. The results from the study revealed that the super absorbent polymer had no effect on the soil wetting pattern while the decrease in the irrigation frequency and increase in irrigation duration led to an enlargement in both the diameter and depth of the soil wetted zone. Information on depths and widths of wetted zone of soil under drip irrigation have a great significance in the design and management of drip irrigation system for delivering required quantity of water and chemicals to plants.

Characterization and prediction of soil hydraulic and solute transport parameters using a random forest model

AbstractTo effectively control nonpoint source pollution and predict its transport and load, understanding water and solute transport processes, patterns and mechanisms is essential. However, the measurement of water movement and solute transport parameters is usually a laborious and time‐consuming task. It is important to predict water movement and solute transport parameters from more readily available soil physical and chemical properties. In this study, a database of soil hydraulic and solute transport parameters containing information retrieved from 83 published studies on soil properties, land use, management measures, etc. was established to characterize and simulate soil water movement and solute transport processes. Our results showed that the soil particle composition was closely related to all soil hydraulic and solute transport parameters. As the soil texture changed from sand to clay, the soil residual water content (θr) obviously increased. The soil porosity was significantly positively correlated with θr, saturated water content (θs), van Genuchten parameter α, dispersity (λ) and retardation factor (R) and negatively correlated with van Genuchten parameter n (p < 0.01). The use of random forest models allowed the prediction of soil hydraulic and solute transport parameters by inputting common or even incomplete soil property parameters. The bulk density and particle composition jointly contributed 66% to the prediction of the soil hydraulic parameters and 44% to the prediction of the solute transport parameters. The pH exerted a notable impact on the solute transport parameters, especially dispersion coefficient (D), λ and R. The results may be useful in providing data support and facilitating the development of watershed nonpoint source pollution models.

Soil water movement may regulate soil water consumption and improve cotton yields under different cotton cropping systems

By quantifying the soil water movement (SWM) in crop planting systems, we can better understand the soil water consumption (SWC) and crop yield relationship; this finding is significant for determining the field water cycle and reducing agricultural water waste. In this paper, a case study was conducted on cotton production. Soil moisture sensors were set at depths of 10–110 cm under three cotton cropping systems (monoculture cotton (MC), wheat/delayed intercropped cotton (WIC), and wheat/direct-seeded cotton (WDC)) based on spatial grid methods; a geostatistical grid calculus was used to calculate SWM; and the crop and meteorological influence mechanisms on cotton lint yield were comprehensively analyzed. At the squaring stage, SWC and vertical SWM were significantly correlated with light, temperature and water conditions. At the flowering and boll development stage, SWC and vertical SWM were collectively affected by meteorological conditions and crops, and they were positively correlated with lint yield. The aboveground and belowground biomass accumulation at the flowering and boll development stage positively affected vertical SWM in and between cotton rows. Vertical SWM in cotton rows increased SWC in cotton rows. SWC in cotton rows and aboveground biomass positively impacted lint yield formation; SWC between rows negatively impacted lint yield. The SWC and vertical SWM between rows in the MC seedling stage exceeded those in cotton rows, and more precise irrigation at the seedling stage reduced water waste. The WIC horizontal SWC at the squaring and flowering and boll opening stages was relatively high, moving from the row midline to cotton row. A better SWC distribution in and between cotton rows promoted water utilization in the cotton rows; this method was feasible for improving cotton yield in diverse planting systems. The results could optimize precision irrigation management at different cotton growth stages and provide a theoretical reference for promoting sustainable agricultural production and climate adaptation.

High spatiotemporal asynchrony between water deficit and nitrate accumulation under apple orchards in deep loess deposits

Understanding the relationship between water and nitrate in soils is critical to agricultural production and aquifer vulnerability assessments, especially in regions with intensive agricultural activities and thick unsaturated zones. Deep soils usually record long-term variations in coupling water and nitrate processes, influenced by multiple environmental factors. As such, this study aims to investigate water and nitrate-nitrogen (NO3–−N) contents and their relationships in soil profiles of 13–23 m deep under farmland and neighboring apple orchards with different ages. Soil properties were considered to explore the individual effect of a single factor or combined effects of multiple factors by wavelet analysis and partial least squares path models. Compared with farmland, apple orchards had large water deficit and NO3–−N accumulation due to root water uptake and synthetic fertilizer inputs, with high spatiotemporal asynchrony in deep unsaturated zones. Specifically, water deficit and NO3–−N accumulation were respectively observed in relatively stable (5–15 m) and active layers (0–5 m, spatial asynchrony), with higher water deficit rates under young apple trees (e.g., 12 and 22 years old) but larger NO3–−N accumulation rates under old apple trees (e.g., 23–32 years old, temporal asynchrony), respectively. Under farmland and orchard, soil texture-organic carbon regulated soil water movement, while soil water-organic carbon affected nitrate transport and transformation processes, which further altered soil pH and electric conductivity conditions in deep vadose zones. The growth of apple trees leads to asynchronous water depletion and fertilizer accumulation in soils, and the > 20 years apple trees are limited by insufficient water despite the excessive fertilizer. The dried soil layers not only reduce potential groundwater recharge, but also prolong the residence time of nitrate bombs in thick unsaturated zones. This study provides essential information for sustainable management of vegetation and water resources, and benefits the complicated hydrological and biogeochemical cycles in thick vadose zones.

Effects of Burial Furrow Parameters on Soil Water Movement under Subsurface Stalk Composite Pipe Irrigation

Stalk composite pipe (SCP), a novel product used for subsurface irrigation, is primarily composed of crop stalks, soil, and water. In subsurface SCP irrigation (SSI), comprehending the effects of burial furrow parameters on soil water movement is critical. However, the precise effects are not known. Here, we aimed to assess the effects of varying burial furrow dimensions, including widths of 15, 20, and 25 cm, depths of 10, 15, and 20 cm, and backfilling-soil bulk densities of 1.1, 1.2, and 1.3 g cm−3, on wetted distance and soil water movement. We conducted an indoor soil bin experiment and numerical simulation using HYDRUS-2D in 2020–2021 in Zhengzhou City, Northern China. Wider burial furrows resulted in increased wetted distance and soil water content, ultimately leading to greater effect in a horizontal direction. In the horizontal profile, the average soil water content of the SSI treatments, which used burial furrow widths of 15, 20, and 25 cm, were 1.121, 1.230, and 1.280 times higher, respectively, than those of CK. The burial furrow depth had minimal effect on the wetted distance, whereas adjusting the burial furrow depth primarily affected the saturation depth. The speed of wetting-front migration was affected by the backfilling-soil bulk density. For a given soil and crop, it is advisable to determine first the burial furrow width and backfilling-soil bulk density. Subsequently, the burial furrow depth should be established based on the distribution of the crop’s roots. Our findings offer a scientific basis for using SSI.

Effects of soil heterogeneity and preferential flow on the water flow and isotope transport in an experimental hillslope

Soil water movement plays vital roles in hillslope runoff generation and groundwater and surface water interaction. However, there are still knowledge gaps about the impacts of soil heterogeneity and preferential flow on the internal water flow and transport process. In this study, the vertical soil heterogeneity focused on the variations in soil retention capacity, and the consideration of lateral preferential flow emphasized the higher hydraulic conductivity. We combined isotopic tracing and numerical modeling in an artificial hillslope, focusing on monitored processes of the artificial rainfall and isotopic tracing experiment. The results showed that the soil moisture quickly accumulated at the bottom of the hillslope during rainfall events, while the 2H enrichment occurred in the topsoil derived from enriched isotope injection in the second artificial rainfall. The evaporation process slowed down the mixing of new water in the topsoil and old water in the lower layer. We found that the vertical soil heterogeneity had significant influences on the internal water and isotope transport paths within the hillslope. The lateral preferential flow played an important role in the water flux and transport time to the seepage face. The coupling of isotopic tracing, which reflects the water transport and mixing with the hillslope, effectively improved the model simulation and mechanism analysis of hillslope water flow. Our findings provide new insights into the mechanisms governing soil water flow and transport dynamics in hillslopes, taking into account vertical soil heterogeneity and lateral preferential flow.

Effects of Coal Gasification Slag on the Migration of Cd2+ and Pb2+ in Soil

In recent years, coal gasification slag has been widely used as an adsorbent in many countries worldwide. To better understand remediation of soil heavy metal pollution by coal gasification slag, it is necessary to examine and analyze the influence mechanism of the adsorption characteristics of coal gasification slag on the migration of heavy metal ions in the soil. In this study, we used a soil water movement model and convection-dispersion equation and constructed a model of heavy metal migration after soil improvement by coal gasification slag; then, we evaluated and predicted the influence of heavy metal infiltration with liquid on groundwater after soil improvement. Combined with the soil column experiment, HYDRUS-1D software was used to predict the migration behavior of heavy metals in the soil. The results showed that the dispersion of the sample with coal gasification slag was smaller than that of the original soil sample, and the permeability coefficient of the soil was reduced by 59.34%; after adding coal gasification slag, the adsorption rates of Cd2+ and Pb2+ on the soil samples were higher than 95% and 99% respectively. The constructed heavy metal migration model was fitted well and predicted the migration behavior of heavy metals in the soil satisfactorily. The addition of coal gasification slag can effectively prevent the infiltration of Cd2+ and Pb2+ in the soil and can mitigate groundwater pollution. Therefore, using the adsorption characteristics of coal gasification slag to study the migration behavior of Cd2+ and Pb2+ in soil has important practical significance for soil adsorption and provides a new idea for the comprehensive utilization of coal gasification slag.

Matric Flux Potential in Time Independent Infiltration Problems from a Single Triangular and a Trapezoidal Irrigation Channel

In this paper, steady infiltration problems into a homogeneous soil from a single triangular and trapezoidal irrigation channel are considered. The governing equation is Richard's equation that represents the movement of water in unsaturated soil. It is a non-linear equation and can be solved by linearizing to become a modified Helmholtz equation. Dual Reciprocity Boundary Element Methods (DRBEM) are used in this study to numerically solve the modified Helmholtz equation. Therefore, by using a provided solution, the numerical Matric Flux Potential (MFP) is calculated. This method was applied to the homogeneous soil problem of stationer infiltration from triangular and trapezoidal single irrigation. Both numerical solutions were compared. The result show that the MFP value from the triangular irrigation is higher than the trapezoidal irrigation. This indicates that content water from the triangular irrigation channel is higher than the trapezoidal irrigation channel.

A novel vertex-centered finite volume method for solving Richards' equation and its adaptation to local mesh refinement

Accurate and efficient numerical simulations of soil water movement, as described by the highly nonlinear Richards' equation, often require local refinement near recharge or sink/source terms. In this paper, we present a novel numerical scheme for solving Richards' equation. Our approach is based on the vertex-centered finite volume method (VCFVM) and can be easily adapted to locally refined meshes. The proposed scheme offers some key features, including the definition of all unknowns over vertices of the primary mesh, expression of flux crossing dual edges as combinations of hydraulic heads at the vertices of the primary cell, and the capability to handle nonmatching meshes in the presence of local mesh refinement. For performance evaluation, soil water content and soil water potential simulated by the proposed scheme are benchmarked against results produced from HYDRUS (a widely used soil water numerical model) and the observed values in four test cases, including a convergence test case, a synthetic case, a laboratory experiment case and a field experiment case. The comparison results demonstrate the effectiveness and applicability of our scheme across a wide range of soil parameters and boundary conditions.

Методология влияния свободной порозности и водоотдачи на процессы, происходящие в почвогрунтах оросительных каналов с земляным руслом

Purpose of the research: to consider the influence of free porosity and water loss on the filtration of water from channels with an earthen bed. Explore existing calculation methods taking into account the parameters of free porosity and fluid loss. Collect calculation material that can form the basis of a program for calculating the filtration of channels with an earthen bed, taking into account the factors influencing the first and second stages of filtration in the channels. Materials and research methods: The work used theoretical methods for calculating filtration in channels with an earthen bed, presented the calculation of the loss of flow rate Qf for filtration per 1 km of length with a steady free surface for different types of channel cross-sections, considered the basic formulas proving that water supply changes with changes in transverse cross-section of the canal, the change in the stability of the canal slopes is shown, a test calculation of the canal for filtration strength according to the maximum pressure gradients is presented, the need to compare the excess of the critical filtration pressure gradient with the permissible pressure gradient is considered. Results and its discussion. The relationship between the total, capillary and minimum moisture capacity and water yield in the thickness of soils is presented, the dependence of free porosity on the groundwater level, and the change in water permeability of soils over time are considered. As a result of the research, material was collected on the influence of free porosity and water loss on the movement of water in soils and the influence on the increase in groundwater levels in the zone of influence of canals, a methodology was presented for the influence of free porosity and water loss on the processes occurring in the soils of irrigation canals in an earthen channel, which will allow its use for modeling these processes on a computer.