Soil Hydraulic Conductivity Research Articles

Prediction of absolute unsaturated hydraulic conductivity – comparison of four different capillary bundle models

Abstract. To model water, solute, and energy transport in porous media, it is essential to have accurate information about the soil hydraulic properties (SHPs), i.e., the water retention curve (WRC) and the soil hydraulic conductivity curve (HCC). It is important to have reliable data to parameterize these models, but equally critical is the selection of appropriate SHP models. While various expressions for the WRC are frequently compared, the capillary conductivity model proposed by Mualem (1976a) is widely used but rarely compared to alternatives. The objective of this study was to compare four different capillary bundle models in terms of their ability to accurately predict the HCC without scaling the conductivity function by a measured conductivity value. The four capillary bundle models include two simple models proposed by Burdine (1953) and Alexander and Skaggs (1986), which assume a bundle of parallel capillaries with tortuous flow paths, and two more sophisticated models based on statistical cut-and-random-rejoin approaches, namely those proposed by Childs and Collis-George (1950) and the aforementioned model of Mualem (1976a). To examine how the choice of the WRC parameterization affects the adequacy of different capillary bundle models, we utilized four different capillary saturation models in combination with each of the conductivity prediction models, resulting in 16 SHP model schemes. All schemes were calibrated using 12 carefully selected data sets that provided water retention and hydraulic conductivity data over a wide saturation range. Subsequently, the calibrated models were tested and rated by their ability to predict the hydraulic conductivity of 23 independent data sets of soils with varying textures. The statistical cut-and-random-rejoin models, particularly the Mualem (1976a) model, outperformed the simpler capillary bundle models in terms of predictive accuracy. This was independent of the specific WRC model used. Our findings suggest that the widespread use of the Mualem model is justified.

Effect of the Mid-Layer on the Diversion Length and Drainage Performance of a Three-Layer Cover with Capillary Barrier

The capillary barrier is a type of soil cover system commonly used in various geotechnical applications, such as limiting infiltration for slopes or landfills or providing cover for solid waste. It serves to prevent the movement of water through the soil layers by utilizing contrasting particle sizes. This paper focuses on investigating the effect of the granular layer on the performance of a three-layer cover with a capillary barrier, integrating the granular layer within clayey sand. The investigation involved one-dimensional infiltration tests utilizing four uniform granular soils with varying grain sizes. These tests were instrumental in calibrating soil water characteristic curves and hydraulic conductivity curves via back analysis. Subsequently, numerical analyses were conducted using a 15 m long model for each of the four distinct cover types. The results indicated that the fine gravel significantly improved the barrier performance beyond one-dimensional tests, owing to its high permeability and the influence of the slope. After the capillary barrier failure, the intermediate layers transitioned into efficient drainage layers, particularly in the gravel layer with the highest lateral drainage capacity. Clayey sand at the bottom delayed percolation, thereby supporting the conversion of the intermediate layer into an effective drainage component. Overall, the multi-layer system showed superior percolation performance compared to the clayey sand cover lacking a granular layer.

Anisotropy of soil hydraulic properties induced by repeated machinery traffic in a Typic Argiudoll of the Argentinean Pampas Region

The increase in heavy machinery traffic due to crop sequence intensification under no-tillage (NT) management can lead to soil compaction, modifying the soil pore configuration, with negative effects on water dynamics. The objective of this work is to evaluate the residual effects of repeated machinery traffic under different soil water contents on the anisotropy of soil hydraulic properties. The residual effect of repeated machinery traffic under different soil moisture conditions (partially dry, PD, and wet, W) on soil compaction was evaluated and compared with the control treatment (no traffic, T0). Field infiltration using tension disc infiltrometer was carried out, and in laboratory, in order to analyse the anisotropy of soil hydraulic properties infiltration experiments were carried out on vertical and horizontal undisturbed soil samples. No significant differences in field soil hydraulic conductivity K(h) were observed between treatments. The W treatment showed a decrease on the transmission and storage pores in the vertical direction, as compared to T0. Anisotropic behavior of soil hydraulic properties was observed for the compacted treatments. The W treatment showed an increment in the transmission porosity in the horizontal direction. In addition, the W treatment presented higher values of transmission and storage pores in the horizontal direction as compared with the vertical direction. From laboratory infiltration experiments, it was observed a strong reduction on the hydraulic conductivity measured at 0, − 3 and − 6 cm water pressure heads (K0, K3 and K6) in the vertical direction under W treatment. An anisotropic behavior for both compacted treatments was observed, with higher values of K0 in the vertical direction under the PD treatment and higher values of K0 and K3 in the horizontal direction under the W treatment. These results show that machinery traffic under wet conditions causes a disruption in soil porosity, not only with a reduction of the pore space, but also with a change in its configuration.

Characterisation and Mapping of Soils in Major Coffee Growing Regions of Uganda

There is limited information on the soil nutrient status, site specific fertilizer and agronomic recommendations for coffee in Uganda hence limiting its production and productivity. Therefore, this study was undertaken to characterize the soils under coffee farms and provide fertilizer and land management recommendations for coffee farming in selected districts in Uganda. 717 soil samples were collected from 45 coffee growing districts that were purposively selected to represent the major coffee growing regions of Eastern, Northern, Western, Mid-west and West Nile. 35 districts were sampled from Robusta coffee growing areas and 10 districts from Arabica coffee growing areas. Parameters considered in the laboratory analysis included; pH, total organic carbon, total nitrogen, available P, exchangeable Calcium (Ca), Sodium (Na), Magnesium (Mg), and Potassium (K) and micro elements: Iron (Fe), Zinc (Zn), Copper (Cu), Manganese (Mn) and Boron, soil texture, bulk density and hydraulic conductivity. Soil nutrient levels distribution maps for Robusta and Arabica coffee growing regions were generated in ArcGIS for the entire country. The soil chemical and physical properties were subjected to analysis of variance using Genstat 14th edition. The soil mapping results showed that, the overall average soil macro-nutrients concentrations were significantly different across regions (p < 0.05) with Eastern having the highest levels of macro elements (CEC of 19.28 meq/100 g, Base Saturation (BS) of 43.40%, pH of 5.78, N of 0.20%, K of 1.64 meq/100 g, P of 107.68 mg/kg and OM of 3.31%) followed by Western region. No significant difference (p > 0.05) was recorded for the micronutrients across the different regions, except Zn (p < 0.05). This study showed that Phosphorus and Potassium were generally moderate to high levels in most of the regions sampled while Nitrogen and organic matter were moderate to low. Results from this study provide a general picture of the nutrient status across all coffee growing regions in Uganda and highlight the required modifications for increased production and productivity.

Application of Compost and Manure in the Biopore Infiltration Hole to Improve Saturated Hydraulic Conductivity (Shc) of Soil in the Coffee Root Zone

Coffee is one of the primary sources of foreign exchange in Indonesia. The factors causing fluctuations in coffee production are climate change, high soil density, less organic material, which causes coffee roots to be unable to absorb water and nutrients optimally. The soil water movement measured by saturated hydraulic conductivity which represents the soil ability to transmit water. The application of compost and manure in the biopore infiltration hole (BIH) is an effective technology in reducing surface-runoff and organic wastes management, as well as soil characteristics improvement. This study was conducted using a Randomized Block Design with four treatments and four replications. The treatments were P1 as control (without BIH), P2 (BIH without compost), P3 (BIH+compost), and P4 (BIH+goat manure). Data analysis was carried out using the Ftest (ANOVA), and continued with Least Significant Difference (LSD) test at 5% level, when there was a significant difference. Observations of soil characteristics were carried out every two months in three depth of root zone, 0-20 cm, 20-40 cm, and 40-60 cm. While the observation of the Auger hole at a depth of 0-30 cm and 0-60 cm. Research variables include soil texture, soil bulk density, Saturated Hydraulic Conductivity of soil by Auger-Hole, Soil Organic Matter (SOC), pH, Cation Exchange Capasity (CEC), and coffee yield. Results showed that treatment of BIH+goat manure gave the best results, measured by improving soil hydraulic conductivity up to 40%. The highest coffee yield was found in the BIH+manure treatment to 3.29 t ha-1.

Experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering boundary effect

The boundary condition is a crucial factor affecting the permeability variation due to suffusion. An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted, where the hydraulic head difference (ΔH) varies, and the boundary includes non-loss and soil-loss conditions. Soil samples are filled into seven soil storerooms connected in turn. After evaluation, the variation in content of fine sand (ΔRf) and the hydraulic conductivity of soils in each storeroom (Ci) are analyzed. In the non-loss test, the soil sample filling area is divided into runoff, transited, and accumulated areas according to the negative or positive ΔRf values. ΔRf increases from negative to positive along the seepage path, and Ci decreases from runoff area to transited area and then rebounds in accumulated area. In the soil-loss test, all soil sample filling areas belong to the runoff area, where the gentle-loss, strengthened-loss, and alleviated-loss parts are further divided. ΔRf decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part, and Ci increases and then decreases along the seepage path. The relationship between ΔRf and Ci is different with the boundary condition. Ci exponentially decreases with ΔRf in the non-loss test and increases with ΔRf generally in the soil-loss test.

Combined influence of rotated transverse anisotropy in mechanical parameters and hydraulic conductivity of soils on reliability evaluation of an unsaturated slope under rainfalls

Natural slopes often exhibit tilted stratification with rotated transverse anisotropy in multiple soil properties, including the mechanical (e.g., friction angle φ and cohesion c) and hydraulic properties (e.g., saturated hydraulic conductivity ks). This phenomenon indicates that to achieve an accurate assessment of slope reliability under rainfall infiltration, the anisotropic spatial variation of the shear strength parameters and ks should be incorporated in a combined manner. Thus, this paper discusses the combined influence of rotated transverse anisotropy in the shear strength parameters and ks on the reliability and failure mechanism of an unsaturated slope under rainfalls. It is found that for a slope with tilted stratification, the failure mechanism and reliability estimator are dominantly influenced by the rotated transverse anisotropy in the shear strength parameters, while the influence induced by that in ks is slight. In particular, only incorporating the rotated transverse anisotropy in ks may lead to an incorrect reliability estimator of a slope with tilted stratification. Herein, the reliability index β would be estimated to be dramatically higher, and β of a slope with horizontal bedding would be higher than that of an anti-dip slope, which is inconsistent with the engineering experiences. Nonetheless, the rotated transverse anisotropy in ks should not be ignored in slope reliability assessment, because ignoring the rotated transverse anisotropy in ks would lead to overestimation of the reliability index β, which is adverse to the safety design of a slope.

Physical Properties of Retisol under Secondary Pulp and Paper Sludge Application

A positive effect of pulp and paper mill sludges as a rich source of organic substrates on soil properties was previously found for some types of sludge and soil. In this study, the effect of secondary pulp and paper sludge on water characteristics and thermal properties of Retisol, as well as the growth parameters of lettuce (Lactuca sativa L.), was tested on the basis of a pot experiment when watering plants with a 20 or 40% sludge solution. The sludge application enhanced plant growth with an increase in biomass accumulation of 21 and 53%, respectively, for 20 and 40% sludge treatments. When the sludge dose was increased from 0 to 40%, the number of leaves increased by 25%, and the leaf mass per area value increased by 42%. Due to the accumulation of more biomass in the shoots than in the roots, sludge causes a change in the allocation of plant biomass. A significant effect of the sludge application on soil particle and microaggregate sized compositions, as well as on the saturated soil hydraulic conductivity, was not found in this study. However, fitted soil water retention curves showed an increased soil water content in sludge-treated soil at all water content values exceeding field capacity. Secondary sludge application led to an increase in the saturated water content from 0.50 to 0.56 cm3 cm−3. The 40% sludge solution increased soil thermal conductivity from 0.92 to 0.98 W m−1 K−1 under saturated water content and from 0.83 to 0.92 W m−1 K−1 under field capacity. The thermal conductivity was higher in the sludge-treated than untreated soil due to a more pronounced positive effect of increased saturated water content than the negative effect of the increased organic matter content on heat transfer. The positive impact of secondary sludge application on both plant growth parameters and physical properties of Retisol, such as increased soil water-holding capacity and thermal conductivity coefficient confirms the possibility of using it to improve soil characteristics and plant productivity.

Spatiotemporal Heterogeneity of Water Conservation Function and Its Driving Factors in the Upper Yangtze River Basin

The water conservation function (WCF), as one of the most critical ecosystem services, has an important impact on the ecological sustainability of a region. Accurately characterizing the spatiotemporal heterogeneity of WCF and further exploring its driving factors are of great significance for river basin management. Here, the WCF of the upper Yangtze River basin (UYRB) from 1991 to 2020 was calculated using the water yield module in the Integrated Valuation of Ecosystem Service and Tradeoffs (InVEST) model. Also, we innovatively applied emerging hot spot analysis (EHSA), which could describe the location and pattern of historical changes more accurately, to investigate the spatiotemporal heterogeneity and evolution of WCF. Based on the Geographical Detector Model (GDM), the main driving factors of WCF and their interactions were revealed. The results showed the following: (1) the WCF in the UYRB experienced a temporal increase at a growth rate of 1.48 mm/a, while remarkable differences were observed across the change rates of sub-watersheds. (2) The spatial variation of the WCF showed a gradual increase from northwest to southeast. Interestingly, the Jinshajing River upstream (JSJU) source area with a low WCF showed an increasing trend (with diminishing cold spots). On the contrary, the downstream regions of the JSJU watershed (with intensifying cold spots) underwent a weakening WCF. (3) Among all driving factors, precipitation (q = 0.701) exhibited the most remarkable prominent impact on the spatial heterogeneity of the WCF. Additionally, the interaction of factors exhibited more explanatory power than each factor alone, such as precipitation and saturated soil hydraulic conductivity (q = 0.840). This research study is beneficial to water resource management and provides a theoretical basis for ecological restoration.

Drivers of field-saturated soil hydraulic conductivity: Implications for restoring degraded tropical landscapes

Water security represents a major challenge in East Africa, affecting the livelihoods of millions of people and hindering sustainable development. Predicted increases in rainfall intensity and variability are expected to exacerbate water insecurity and land degradation. Improving soil infiltrability is an effective strategy for addressing water insecurity and land degradation. Research on soil infiltrability is often highly localized; therefore, scientific understanding of the drivers of infiltrability on larger spatial scales is limited. The aim of this study was to understand the main drivers of infiltrability across five contrasting landscapes in Kenya. We measured field-saturated hydraulic conductivity (Kfs) in 257 plots and collected data for variables representing soil properties (sand content, soil organic carbon (SOC) and pH), land degradation (grazing pressure and presence of erosion), vegetation quantity (woody aboveground biomass), and vegetation quality (functional properties and diversity). We used generalized mixed-effects models to test for the effects of these variables on Kfs. Median Kfs for the five sites ranged between 23.8 and 101.8 mm h−1. We found that Kfs was positively associated with sand content (standardized effect 0.39), SOC content (0.15), and functional diversity of woody vegetation (0.09), while it had a negative relationship with the presence of erosion (−0.24) and grazing pressure (−0.09). Subsequently, we conclude that infiltrability can be enhanced through using land restoration strategies which specifically target parameters that affect Kfs. The results further support that Kfs is not solely dictated by inherent soil properties, and that management interventions which boost SOC, reduce erosion, and minimize unsustainable grazing can help address water scarcity by restoring soil hydrological function.

Ground Fire Legacy Effects on Water-Dynamics of Volcanic Tropical Soils

The forest floor is a critical component in maintaining the life cycles of forest ecosystems. It normally includes organic soil horizons, known as duff and litter, which are prone to be rapidly consumed after flaming and smoldering fires. This work aims to understand the legacy effects of surface and ground fires on the infiltration capacity of a volcanic forest soil. We studied five sites with fires recorded in the last 20 years. All of them are located in pine-oak forests of the volcanic mountain region in central Mexico with a temperate climate and Andic soil properties. Tension-infiltration tests were carried out to determine hydraulic conductivity and the number of active macropores. After each test, cores were taken to evaluate in a laboratory setting, where soil water repellency at different moisture concentrations and the integrative dynamic repellency index were determined. Field-saturated hydraulic conductivity was moderately high in all sites, with mean values of 13 and 42 mm·h−1 for burned and control plots, respectively. A non-linear relationship was found between recurrence and type of fires with the concentration of active pores and several dynamic water repellency parameters. This work confirmed the presence of latent combustion in these temperate neotropical forests. The changes in soil water repellency and hydraulic conductivity detected do not necessarily imply an exceeded soil infiltration capacity. However, many of the fires in this region are associated with increasing agricultural activities, so further studies are needed to determine if higher fire frequencies could exceed the resilience capacity of the soils triggering land degradation.

Uncertainty Quantification for the Noah‐MP Land Surface Model: A Case Study in a Grassland and Sandy Soil Region

AbstractLand surface models rely on a multitude of parameters to simulate land‐atmosphere interactions, but the parameter uncertainty can limit the reliability of model predictions. This study utilizes a parameter uncertainty quantification (UQ) framework to quantify and reduce the parameter uncertainty of the Noah‐MP land surface model in a grassland and sandy soil region in the Midwest of the USA. First, the sparse polynomial chaos expansion method which can quantify the interaction effect of parameters, is employed. A relatively small parameter sample size (i.e., 20 times of the number of parameters) was sufficient to identify the sensitive parameters; an additional sensitive parameter, the saturated soil hydraulic conductivity, was screened out compared to previous study. Then, based on the selected sensitive parameters, the weighted multi‐objective adaptive surrogate modeling‐based optimization algorithm is used as the parameter optimization method. The optimization results showed that the root mean square error of flux of latent heat (FLH) on about 82% of the total grids was reduced, and the number was about 57% for gross primary production (GPP) compared to the results using the original parameter settings, indicating that the Pareto parameter set by the UQ framework improved the Noah‐MP model in simulating FLH and GPP in a grassland and sandy soil region in the Midwest of the USA.

Influence of landscape position and climatic seasonality on soil water and gas conductivity properties in agricultural soils

Agricultural landscape management and climate seasonality can influence soil structure, hydraulic conductivity, and air permeability within the context of soil water and soil gas mobility. To investigate this, in situ and laboratory-based data were collected from three agricultural landscape positions within a watershed in eastern Ontario, Canada during a growing season. Macropore classification, water infiltration tests, and air permeability measurements were conducted in situ and standard soil characterizations were carried out on soil samples. Hydraulic conductivity of the soil matrix, based on grain size data, indicated that the highest values were consistently measured in the B horizon at each landscape setting. Macropores were found to be more abundant within uncultivated drainage ditch bank soils, compared to the adjacent cropped fields. Macropores in the ditch bank soils were exclusively consisted of circular biopores, while both circular and linear macropores were observed in the cultivated field soils. Air permeability, vertical hydraulic conductivity, and horizontal hydraulic conductivity were also greater in the uncultivated soils, relative to the cultivated soils. Field saturated hydraulic conductivity measurements offered evidence of anisotropy, likely due to the vertical nature of the macropore features. Macropore disposition and extent varied over the growing season, especially in the cultivated field soils where tillage and field trafficking are physically disruptive. Seasonality of macropore development will influence temporal changes in advection-based mass exchange of gas and water in the vadose zone. Modeling of mass exchange in agricultural soils should consider time variability in macroporosity to more realistically characterize infiltration and soil gas emissions.

A refined representation method of the particle size distribution curve based on grading entropy

The particle size distribution (PSD) curve is an important expression of the basic properties of the soil. The characteristic parameters such as median particle size D50, effective particle size D10, and some combinations of characteristic parameters can only represent some points of the PSD curve, failing to represent all information of the PSD curve. In this paper, a fraction characteristic parameter F that could reflect the fraction size variation of the PSD curve was introduced based on grading entropy. And a new presentation method of the A-B-F three parameters was proposed for a refined representation of the PSD curve. It was found that the newly constructed model not only better represents the differences in the width of PSD curves but also has a higher sensitivity than the A-B two-parameter model proposed by Lőrincz. For PSD curves with similar distributions, the new presentation method has a higher degree of discrimination than that of the D50-Cu-Cc three-parameter or the four-parameter combination Gc proposed by Arshad. Finally, the application of this new method in describing the spatial nonuniformity distribution of depositions and predicting the soil hydraulic conductivity was discussed. The research results can provide a reference for the refined representation of PSD curves.

Effect of leachate contamination in municipal solid waste on clay liner characteristics

Landfill leachate is generated as a consequence of water percolation through the solid waste, oxidation of the waste, and corrosion of the waste. Under-designed landfill sites allow the leachate to easily pass through the soil strata. To identify the effect of contaminants on the engineering properties of the clay liner, a pilot-scale lysimeter was developed. Experiments on two types of soils polluted with leachate were conducted which included silty clay and soil enhanced by chemical stabilization. Lime and silica fume are used as reactive barriers in landfill earthen liners. The leachate generated was monitored periodically and sampled for (pH, Electrical conductivity, Chemical oxygen demand, Chloride, Total dissolved solids, and total suspended solids). Lime (L), silica fume (SF), and lime-silica fume (L-SF) mix have been used for stabilizing soil. Three percentages were used for lime (2%, 4% and 6%) and three percentages were used for silica fume (3%, 5% and 7%) and the optimum percentage of silica fume (5%) was mixed with the proportions of lime. The optimum lime values are 4% from the research and 5% for silica fume. Soil samples were analyzed over a maturation period of 135 days for each run to identify the chemicals and their effect on the soil characteristics. The plasticity index of the stabilized soil continues to decline until it reaches 36% over the 135 days of the maturing period. Furthermore, an increase in silty clay soil hydraulic conductivity was observed from 3.3 × 10−6 cm/sec to 5 × 10−6 cm/sec at 135 days of maturation. Whereas there is a sharp decrease in hydraulic conductivity to achieve a maximum value in the stabilized soil (1.25 ×10−10 cm/sec). Based on laboratory test results, it can be concluded that soil stabilization can improve the geotechnical properties of landfill earthen liners due to pozzolanic reactions.

Infiltration Measurements during Dry Conditions in an Urban Park in Ljubljana, Slovenia

A thorough understanding of the hydrologic mechanisms that control the movement of water through the soil is essential for developing effective stormwater management strategies. Infiltration is critical for determining the amount of water entering the soil and controlling surface runoff. Spatial and temporal variations in soil properties strongly affect infiltration rates, which underscores the importance of evaluating field-specific values for hydraulic conductivity, which are also highly dependent on the chosen measurement and evaluation methods. The objective of this study is to determine and compare soil hydraulic conductivity under dry conditions using two field measurement techniques, namely the double-ring infiltrometer (DRI) and the mini-disk infiltrometer (MDI). The results demonstrate the importance of performing multiple replicates of infiltration tests, especially during the dry season, as the initial dry surface caused deviations in hydraulic conductivity estimates for both methods used (DRI and MDI). Significant spatial variability was observed within the radius of the test replicates over short distances (<1 m). In addition, experimental infiltration curves for a selected site were used to evaluate and compare soil hydraulic parameters through infiltration modeling. In general, the Philip, Green-Ampt, and Smith-Parlange theoretical models showed a better fit to the experimental DRI data than the semi-empirical Horton model.

Quantification of red soil macropores affected by slope erosion and sediment using computed tomography

AbstractSoil structure is an important factor interacting with soil erosion and sediment processes. However, few studies have focused on the relationship between soil macroporosity and soil erosion across different terrains. The aim of this study was to quantify and compare soil properties and macroporosity characteristics in collapsing gully areas and to explore their impact on the formation and development of collapsing gullies. Soil cores were excavated at different positions of a typical collapsing gully and then scanned to analyze soil macropores. Soil properties and saturated hydraulic conductivity were also investigated. The results showed that the contents of sand, silt, and clay, the mean weight diameter of aggregates, and the infiltrate rates varied at different positions. The valley had the greatest macroporosity (1.09% ± 0.33%), the number (5919 ± 703), volume (1468 ± 194 mm3), and surface area (10.4 ± 2.6 m2) of macropores, as well as the mean volume (16.8 ± 7.4 mm3) of macropores >1 mm3, whereas these indices were lowest at the slope (0.15% ± 0.14%, 1189 ± 747, 266 ± 188 mm3, 1.7 ± 1.4 m2, and 10.6 ± 2.9 mm3, respectively). The macroporosity and the number of macropore decreased with increasing depth but were also influenced by the erosion and sediment processes. The processes of sediment and the roots of vegetation also influenced the orientation of the macropores. Macropore characteristics at different sites of the collapsing gullies affected the soil water infiltration and hydraulic conductivity and further affected the processes of water erosion and mass erosion.

New procedure to estimate soil field capacity based on double ring infiltration data

Many models have been proposed to explain the water infiltration curve in soil, but according to the studies we have done so far, no model can simultaneously estimate the field capacity (FC) and soil hydraulic conductivity (K). The aim of this research is to present a new model that is defined based on the dynamic balance of water-soil-time and is able to estimate FC and soil hydraulic conductivity using cumulative field infiltration data. The studies were carried out in Varamin region with an area of 1370 square kilometers. For the penetration test, first, seven sides (A, B, C, D, E, F and G) were selected in 7 areas of agricultural land and with a double ring with an outer diameter (70 cm) and an inner diameter (30 cm) became. Also, to determine the soil texture, apparent specific gravity, PWP and FC, soil sampling was done at the depths of 0–10 and 20–30 cm. To determine FC moisture, pressure (−10 kPa) was used in a pressurized membrane device, and to determine volumetric moisture in PWP, sieved samples were placed on ceramic plates. Using a spray, the samples were wetted, so that they reached the saturation level, and the pressure was set to 1500 kPa. Examining the FC values showed that the maximum value of FC in site F (clay loam texture) is equal to 14.73 cm and the lowest value in site G (sandy clay loam texture) is equal to 6.8 cm by the model estimated by the results obtained in Laboratory confirmed. Meanwhile, the minimum saturated hydraulic conductivity of soil at site A was estimated to be 0.06 cm/min and the maximum at sites D and F was 0.28 cm/min. Time constant (Ksh) has shown completely different water infiltration rate in the first term of the model (Rfksh) compared to the second term of it (K), and no obvious relation has detected between them. But, their summation (k+Rfksh) was always equal to the maximum infiltration instantaneous rate (T→0). The research results showed that using the new model, the FC value varies between 6.8 and 14.73 cm. Also, the minimum and maximum instantaneous speed of water infiltration occurred in the site. An equal to 0.621 and in site D equal to 1.735 cm/min. The results of fitting the penetration model on cumulative penetration data showed that R2 varies between 0.994 and 0.999 and RMSE between 0.116 and 1.197. Therefore, the model is able to describe the penetration data well.

Comparison of Soil Hydraulic Properties Estimated by Steady- and Unsteady-Flow Methods in the Laboratory

In this study, soil hydraulic conductivity (K) and soil sorptivity (S) values estimated by applying various steady- and unsteady-flow methods using cumulative infiltration data of three disturbed soils (sandy loam, loam, clay) obtained from a disc infiltrometer in the laboratory at various negative pressure heads were compared. The steady-flow methods used were those of Ankeny et al. and Reynolds and Elrick as well as Logsdon and Jaynes, while the unsteady-flow methods were those of Haverkamp et al. (two-term (2T) and three-term (3T) infiltration equations) and Zhang. The method of White et al., which is a steady-flow method but also uses unsteady-flow infiltration data, was also examined. The results showed that the three steady-flow methods, as well as the Zhang equation, for values of the van Genuchten coefficient n > 1.35, tend to give similar values of K. The 2T infiltration equation with β = 0.6 provided hydraulic conductivity values greater than those estimated by the steady-state methods but gave negative K values in some cases. The values of the coefficients C1 and C2 of the 2T equation were affected by the infiltration time. The coefficient C1 increased while C2 decreased with increasing time when the cumulative linearization method (CL) was applied, but the change in C1 tended to be smaller than that in C2. The inverse solution of the 3T equation using the Excel Solver application for β = 0.75 and β = 1.6, when positive values of K were obtained, approached better the K values estimated by the steady-flow methods compared with those estimated using β = 0.6. Regarding the estimation of S from the unsteady-flow equations (2T, 3T, Zhang), comparable S values were obtained by all equations. The differences between the S values of the various methods are smaller compared to those of K, and S is less affected than K in terms of time. The problem of negative estimates of K might be attributed to the fact that the soils used in this study are classified as soils situated in the domain of lateral capillarity or are not completely homogeneous or soil compaction is observed at some depth. In the case where the soils are not completely homogeneous, the Sequential Infiltration Analysis (SIA) method with β = 0.75 corresponding to the soil types studied was proved to be effective in estimating K values.

New Method for Modelling Seasonal Variation in Resistance and Performance of Earthing Systems

The current techniques utilized for estimating seasonal fluctuations in earthing system resistance, including artificial neural networks (ANNs) and correlation/correction factors, rely on resistance records, soil resistivity measurements, and meteorological data collected across broad areas. However, they frequently fail to consider the impact of soil conditions and properties at the actual earthing location. As a solution, this research introduces a new method that models atmospheric conditions as soil suction and incorporates hydraulic soil properties (soil water retention characteristics and hydraulic conductivity) to estimate the seasonal changes in earthing resistance and performance. To illustrate this approach, this study constructs geometric models of vertical earthing rods for three homogeneous soil textures (clayey, silty, and sandy) utilizing COMSOL Multiphysics software. By coupling the differential equations governing electric current and water flow using Archie’s formula and solving numerically with the finite element method (FEM) for various soil suctions, this research reveals that soil water retention and resistivity variations are notably influenced by soil texture. Sandy soil displays higher variability, silt soil demonstrates moderate changes, while clayey soil exhibits lower fluctuation. By linking soil resistivity changes to soil suction and hydraulic properties, this innovative method predicts seasonal trends in earthing resistance and performance.