Peculiarities of Hydrodynamic and Electro-osmotic Flow of Water in Clay Soils

S
T
T. Yedoyan2,*
E
E. Neverova-Dziopak3
V
1Faculty of Construction, National University of Architecture and Construction of Armenia, 105 Teryan St., Yerevan 0009, Republic of Armenia.
2Faculty of Management and Technology, National University of Architecture and Construction of Armenia, 105 Teryan St., Yerevan 0009, Republic of Armenia.
3Faculty of Mining Surveying and Environmental Engineering, AGH University of Science and Technology, Kraków, Poland.
4Department of Mathematics, Construction Mechanics and Physics, National University of Architecture and Construction of Armenia, 105 Teryan St., Yerevan 0009, Republic of Armenia.

Background: Global climate change, accompanied by recurrent droughts in arid and semi-arid regions, necessitates the development of more effective approaches for regulating soil water regimes. Conventional hydrodynamic methods often perform poorly in soils with thick clay horizons and low permeability, where electrokinetic processes may provide alternative pathways for enhancing water movement.

Methods: This study evaluated the parameters of hydrodynamic and electro-osmotic flow in irrigated meadow clay soils of Armenia. Laboratory experiments were conducted under controlled hydrostatic pressure gradients and applied electrical field gradients to characterize filtration behaviour and electro-osmotic response.

Result: A strong curvilinear relationship was observed between filtration rate and pressure gradient (R = 0.983), indicating plastic-viscous behaviour of the soil solution. A similar curvilinear dependence was identified between the electro-osmotic coefficient and electrical field gradient (R = 0.884). Electro-osmotic effects substantially increased the water permeability when the electrical field strength exceeded 2.57 V cm-1. These findings demonstrated the potential of electro-osmotic phenomena to enhance water movement in heavy clay soils, offering promising applications for soil-water management under changing climatic conditions.

Efficient management of soil water regimes is critical for sustainable agriculture, irrigation planning and soil amelioration, particularly in heavy-textured, low-permeability soils (Kotorová and Mati, 2008; Lakshmi et al., 2016; Kumari et al., 2021; Tamiru et al., 2023). Conventional hydraulic methods, such as drainage and irrigation systems, often fail to ensure uniform water distribution in clay-rich horizons due to limited water movement (Frenkel et al., 1978; Kotorová et al., 2013; Machikowa et al., 2020; Balaji and Pandiarajan, 2022). This constraint poses serious challenges for crop productivity, irrigation efficiency and reclamation of saline or waterlogged soils, especially in arid and semi-arid regions (McNeal and Coleman, 1966; Kumar et al., 2015; Khoirunnisak et al., 2024). These issues are directly relevant to farmers, irrigation managers and policymakers concerned with agricultural water security. Electrokinetic processes, particularly electro-osmosis, offer a promising alternative by facilitating water movement through soil pores under an external electric field. In contrast to hydraulic flow, electro-osmotic transport operates independently of pressure gradients, enabling water redistribution in clay soils and improving permeability (Lockhart, 1986; Smollen and Kaffar, 1994; Miller et al., 1997; Zheng and Zhu, 2017; Butnan et al., 2024). Although electro-osmotic phenomena are well studied theoretically, their practical use in heavy soils remains limited, underscoring the need for experimental evaluation of water movement under combined hydraulic and electrical influences (Gucal and Khyamyalyaynen, 2020).
       
Traditionally, soil water movement is described by Darcy’s law (Darcy, 1856), which assumes a linear relationship between flow rate and pressure gradient. However, numerous studies have shown that this linearity often breaks down under certain conditions, such as high Reynolds numbers, low flow rates, or when fluids display non-ideal rheological behavior (Kotov and Nerpin, 1958; Nagy and Karadi, 1961; Nerpin and Chudnovsky, 1967). Experimental data indicated that flow frequently initiates only after a threshold pressure gradient is exceeded, underscoring the role of viscosity in flow initiation and behaviour (Gomboš, 2012). Depending on rheological properties, fluids are classified as Newtonian, non-Newtonian, or visco-plastic. These distinctions are particularly relevant for water transport in fine-textured soils, where capillary forces strongly influence hydraulic behaviour (Joshi, 2017; Qi et al., 2018).
       
Electro-osmosis provides an additional mechanism of water movement in which ions in the diffuse layer migrate toward electrodes under an external electric field, carrying water molecules and generating directed pore flow. In clay soils, electro-osmotic transport is not simply proportional to electric field strength, requiring modified models that incorporate soil structure and pore characteristics (Bondarenko, 1973). Comparisons of hydraulic and electro-osmotic mechanisms show that although both depend on pore diameter, electro-osmotic flow declines more slowly with decreasing pore size. Thus, applying direct electric current to clay-rich soils can significantly enhance permeability, offering an effective tool for water management and engineering reclamation practices (Stefanidis, 2021).
               
This study extends existing theoretical research by experimentally investigating water movement in heavy-textured soils under the combined influence of hydrostatic pressure and external electric fields. The objective is to identify the key hydrodynamic parameters governing these processes, thereby providing insights essential for optimizing soil water regulation and improving the efficiency of reclamation and irrigation systems.
Study area and soil sampling
 
The study focused on irrigated wet meadow-gray soils (classified as Anthrosols according to the WRB system) from the Masis region of Armenia. Soil samples were collected from 0-50 cm layer, air-dried, gently crushed and sieved through a 2 mm mesh.
       
Experiments were conducted during 2023-2024 at the Laboratory for Environmental Issues, Conservation of Water Resources and Their Efficient Use of the National University of Architecture and Construction of Armenia.
       
Physical, chemical and physicochemical properties of the soils were analyzed both with and without application of an electric current. Bulk density and total porosity in the soils of 0-50 cm layer were 1,300 kg m-3 and 53%, respectively (Table 1). Wilting point, field capacity and plant-available water were 14.2%, 29.3% and 15.1%, respectively.

Table 1: Physical and water-physical indicators of the (0-50 cm) soil.


       
The soil exhibited a light clay mechanical texture, with 48.6% physical clay, including 19.9% silt (Table 2). The top 0-50 cm soil contained 2.1% humus. Calcium carbonate content was high (21.5%), magnesium carbonate was 2.4% and exchangeable cations totaled 31.5 cmolc kg-1. The exchangeable sodium percentage (ESP) was 4.8%.

Table 2: Physical, chemical and physicochemical properties of the (0-50 cm) soil.


       
The indicators of water-soluble salts in the 0-50 cm soil layer are presented in Table 3. The soil is characterized by a total salt content of 0.196%, an electrical conductivity (EC) of 1.78 dS m-1, a pH of 7.5 and a water-soluble sodium concentration of 2.1 cmolc kg-1. Carbonate (CO32) was not detected.

Table 3: Chemical composition of soil water extract (0-50 cm layer).


 
Experimental setup
 
Laboratory experiments were carried out using a horizontally mounted soil column (1) with a length of 0.35 m (L) and diameter of 0.07 m (d), corresponding to a cross-sectional area of 0.00355 m2 (A). Sand-gravel filters (2) were installed at both ends to prevent mechanical suffusion. A feeding tank (3) at the inlet maintained a constant water level (H1), while filtrate was collected in a glass container (4) at the outlet (Fig 1).

Fig 1: Experimental setup.


       
The column was packed with dry meadow-brown soil compacted to the natural bulk density of 1,300 kg m-3. Two titanium electrodes, anode (5) and cathode (6), were embedded in the soil and connected to a direct current source to induce electroosmotic flow.
 
Calculation of hydraulic and electro-osmotic flow parameters
 
During the experiments, various pressure gradients (I) were applied through the feeding tank (1, Fig 1) and the corresponding filtrate volumes (W) were measured. Each experimental run was conducted in triplicate to ensure reproducibility. The pressure gradient (I), filtration rate (VI) and filtration coefficient (KI) were calculated according to Equations (1)-(3).





 
Where,
(H1 - H2)= The hydraulic pressure across the soil column of  length L.
t= The time at which the filtrate volume is collected.
A= The cross-sectional area of the column.
       
For electro-osmotic filtration, both hydraulic pressure gradients and electric field voltages were applied to the soil column. The electric field strength (E), the electro-osmotic velocity (VE) and the electro-osmotic coefficient (KE) were determined using the Equations (4) to (6):





   
Where,
U= The voltage applied to the electrodes.
QE= The volume of filtrate from electro-osmotic flow, which was further calculated using Equation (7):
 
        QE = QT - QI                  ...(7)
                                                                  
Where,
QT= The combined hydraulic and electro-osmotic flow.
QI = Hydraulic flow alone.
       
In soils with heavy mechanical texture, the relationship between filtration rate and pressure gradient often deviates from Darcy’s Law, exhibiting nonlinear behaviour. To account for this, both linear and nonlinear regression models were applied to the experimental data, expressed by Equations (8) and (9):
 
                  V = A + BI                           ...(8)
 
               V = aI3 + bI2 + cI + d                 ...(9)
 
Where,
A, B and a, b, c, d= Regression parameters obtained from statistical fitting of the experimental results.
               
All experimental data were processed using Origin 6.1 software to calculate mean values, determine regression parameters and verify the consistency and reliability of the results.  
Hydrodynamic movement of water in soil
 
The experimental results of hydraulic filtration in clay soils are summarized in Table 4, with replicate data included for reliability. These results indicate that, although a general proportionality exists between pressure gradient (I) and filtration rate (VI), the relationship is not strictly linear (Liu and Birkholzer, 2012). With increasing I, both V and the apparent filtration coefficient (KI) rise, however, the variation in KI demonstrates that Darcy’s law cannot fully describe water movement in these soils.

Table 4: Relationship between pressure gradient (I), filtration rate (VI) and filtration coefficient (KI).


       
Statistical evaluation of the relationship between V and I, using both linear and nonlinear regression models, is shown in Table 5 and illustrated in Fig 2. For linear regression, the correlation coefficient was R = 0.947, with a standard deviation (SD) of 1.16% and relative error (E) of 12.9%. In contrast, nonlinear regression produced a stronger correlation (R = 0.983), with lower SD (0.63%) and error (7.0%). These results clearly demonstrate that water flow in clay soils is better described by nonlinear rather than linear models, consistent with theoretical expectations for visco-plastic systems and similar to previously reported non-Darcian seepage behaviour in clay-rich porous media (Yin et al., 2024).

Table 5: Statistical indicators of linear and nonlinear regression analysis describing the relationship between filtration rate (VI) and pressure gradient (I).



Fig 2: Linear and nonlinear regression analysis of the relationship between filtration rate (V) and pressure gradient (I).


       
The corresponding regression dependencies are presented in Equations (8) and (9), with parameter values summarized in Table 6. In particular, the rectilinear regression yielded a mean hydraulic conductivity (KI) of 20.43 cm d{ ¹, while nonlinear regression produced a cubic parabolic form that accurately reflects the threshold behavior observed at low gradients. As shown in Equation (9), water movement begins only when I exceeds the limiting gradient I0, after which the flow rate increases rapidly.

Table 6: Numerical values of the parameters of equations (8) and (9).


       
Thus, regression analysis confirms that the filtration process in clay soils deviates from Darcy’s law, with water exhibiting visco-plastic properties. The cubic parabolic regression more accurately describes this nonlinear behavior and can therefore be applied in hydraulic and hydro-ameliorative calculations where precision is required.
 
Combined hydrodynamic and electroosmotic water movement
 
Analysis of the experimental data (Table 7) indicates that the combined action of hydraulic pressure and electric field markedly enhances the filtration rate. The increase is particularly pronounced at higher hydraulic gradients and elevated electric field intensities, highlighting the synergistic effect of electroosmosis in facilitating water movement through clay soils. For example, at I = 0.46 and E = 0 V cm-1, the filtration rate (VI) was 7.48 cm d-1, whereas under the same hydraulic conditions but with E = 2.57 V cm-1, the rate increased to 9.58 cm d-1. When the pressure gradient was raised to I = 0.80 and E = 3.43 V cm-1, the filtration rate reached 31.8 cm d-1.

Table 7: Filtration rates of clay soil under varying hydraulic pressure gradients and electric field intensities (cm d-¹).


       
The nonlinear dependence of total flow on electric field intensity is clearly illustrated in Fig 3. At lower electric field strengths (0-2.57 V cm-¹), the increase in filtration rate ranged from 2.1 to 3.3 cm d-1, corresponding to hydraulic gradients of 0.46-0.80. At higher electric fields (2.57-3.43 V cm-1), the rate of increase was much greater, ranging from 7.5 to 14.8 cm d-1. These results show that electroosmosis exerts disproportionately strong effect once the electric field exceeds a threshold value, a phenomenon also reported in controlled laboratory studies on clayey soils (Jayasekera, 2004).

Fig 3: Dependence of total filtration rate on pressure gradient and electric field voltage in clay soils.


       
To clarify the contribution of electro-osmosis, the ratio of total flow under combined hydraulic and electroosmotic forces (VI,E) to hydraulic flow alone (VI) was analyzed. The resulting relative values are summarized in Table 8 and illustrated in Fig 4. The data clearly demonstrate that the application of an electric field substantially enhances water movement through the soil, as reflected by the increasing VI, E/VI ratios with rising voltage intensity.

Table 8: Relative filtration rates (VI, E/VI) of clay soil under different hydraulic pressure gradients and electric field strengths.



Fig 4: Dependence of the relative flow ratio (VI,E/VI) on electric field strength (E).


       
Statistical analysis confirmed the robustness of this relationship, revealing a strong positive correlation between the relative increase in flow and the applied electric field (R = 0.93, SD = 0.129, N = 25, P<0.0001). This dependency is quantitatively described as:


The fitted regression curve is presented in Fig 4, highlighting the consistent and predictable influence of electroosmosis on soil permeability.
       
These results confirm the theoretical prediction, that electroosmotic effects are particularly significant in clay soils. Because hydraulic flow in capillary decreases proportionally to R4 with diameters, while electroosmotic flow decreases only to , the electric field becomes increasingly effective as particles becomes finer (Mitchell and Soga, 2005).
 
Electroosmotic coefficient (KE)
 
The dependence of KE on electric field intensity is shown in Fig 5, with regression statistics provided in Table 9.

Fig 5: Rectilinear and curvilinear dependence of the electro-osmosis coefficient (KE) on electric field strength (E).



Table 9: Statistical indicators of rectilinear and curvilinear regression analysis between the electro-osmosis coefficient and the gradient of the electric field.


       
Linear regression yielded a correlation coefficient of R = 0.846 and SD = 0.251, whereas nonlinear regression provided a better fit with R = 0.884 and SD = 0.228. The regression equations describing this relationship are:
 
                        KE = 10-5 E (1.374 - 0.590E + 0.096E2)             ...(11)      
 
                             KE = 10-5 (0.173 + 0.419 E)                    ...(12) 
 
Parameter values are summarized in Table 10. Within the experimental range of E = 0.57–3.43 Vcm-1, KE varied between 0.525 × 10-5 and 1.525 × 10-5 cm² s-1 V-1.

Table 10: Numerical values of the parameters of equations (11) and (12).


       
These results agree well with published ranges (4.91 × 10-6 to 1.57 × 10-5 cm2 s-1 V-1) reported by other researchers, confirming both the accuracy of the present experiments and the applicability of electroosmotic enhancement in heavy soils (Punia and Singh, 2018).
 
Implications for soil reclamation
 
The combined hydraulic and electroosmotic experiments demonstrate that the application of a direct electric current can significantly improve soil permeability in clay soils. This enhancement results from the fact that electroosmotic flow is less sensitive to pore size than hydraulic flow. Consequently, in clay soils where Darcy’s law fails and hydraulic movement is severely restricted, electroosmosis provides a powerful mechanism for accelerating water movement.
               
For reclamation of saline-sodic soils in Armenia, these findings are highly relevant. By reducing the time and water required for leaching, electro reclamation methods can improve the efficiency and sustainability of irrigation practices. The experimental confirmation of nonlinear relationships, threshold gradients and strong electroosmotic responses provides a robust scientific basis for implementing these technologies in practice. 
This study demonstrated that water movement in irrigated meadow–gray soils with clay mechanical composition does not fully confirm to Darcy’s law. Filtration rates within a pressure gradient range of 0.024 -0.76 exhibited both rectilinear (R = 0.947) and curvilinear (R = 0.983) dependencies on the applied gradient, with the presence of an initial threshold (I0) confirming the plastic-viscous behaviour of the soil solution.
       
The combined influence of hydrostatic pressure and electric fields revealed that electro-osmotic effects are negligible below 2.57 V cm-1 but become significant at higher field strengths, substantially enhancing soil permeability. Regression-based models established strong correlations between relative flow velocities and applied fields, as well as between the electro-osmosis coefficient (Ke) and electric field intensity (R = 0.84 - 0.88), highlighting the nonlinear character of electro-osmotic transport.
               
These findings emphasize the potential of electro-reclamation as an effective method for improving water movement in clay soils. Beyond providing a theoretical and experimental basis for engineering applications, the results hold practical implications for multiple stakeholders. For researchers, they contribute to refining models of soil-water interactions under coupled hydraulic-electrical fields. For engineers and irrigation planners, they suggest new design strategies for drainage and amelioration systems in clay soils. For policymakers and agricultural managers, they demonstrate the feasibility of incorporating electro-reclamation into water-saving and soil-rehabilitation programs, particularly in arid and semi-arid regions where efficient water use is a national priority.
The authors gratefully acknowledge the Laboratory for the Issues of Environment, Conservation of Water Resources and Their Efficient Use at the National University of Architecture and Construction of Armenia for providing the facilities and technical support required for conducting the laboratory experiments.
 
Disclaimers
 
The views and conclusions expressed in this article are solely those of the authors and do not necessarily represent the views of their affiliated institutions. The authors are responsible for the accuracy and completeness of the information provided, but do not accept any liability for any direct or indirect losses resulting from the use of this content.
The authors declare that there are no conflicts of interest regarding the publication of this article. No funding or sponsorship influenced the design of the study, data collection, analysis, decision to publish, or preparation of the manuscript.

  1. Balaji, K. and Pandiarajan, T. (2022). Performance assessment of soil moisture meter under sandy clay loam soil. Agricultural Reviews. 43(2): 243-248. doi: 10.18805/ag.R-2286.

  2. Bondarenko, N.F. (1973). Physics of Groundwater Movement. Gidrometeoizdat, Leningrad. pp. 215. (in Russian).

  3. Butnan, S., Sriraj, P. and Toomsan, B. (2024). Effects of electrolytes and soil-to-suspension ratios on pH in acidic coarse- textured soil. Indian Journal of Agricultural Research. 58(1): 120-125. doi: 10.18805/IJARe.AF-785.

  4. Darcy, H. (1856). Les Fontaines Publiques de la Ville de Dijon. Victor Dalmont, Paris. pp. 647.

  5. Frenkel, H., Goertzen, J. and Rhoades, J.D. (1978). Effects of clay type and content, exchangeable sodium percentage and electrolyte concentration on clay dispersion and soil hydraulic conductivity. Soil Science Society of America Journal. 42: 32-39.

  6. Gomboš, M. (2012). The Impact of Clay Minerals on Soil Hydrological Processes. In: Hydraulic Conductivity-Issues, Determination and Applications. IntechOpen. pp. 119-149. doi: 10.5772/ 47748.

  7. Gucal, M. and Khyamyalyaynen, V. (2020). Hydrodynamic parameters in various methods of soil electrokinetic treatment. E3S Web of Conferences. 174: 01062. doi: 10.1051/e3sconf/ 202017401062.

  8. Jayasekera, S. (2004). Electroosmotic consolidation and flow in clays. Australian Geomechanics. 39: 77-86.

  9. Joshi, R. (2017). Physical constraints of fine-textured heavy soils and their management-A review. Agricultural Reviews. 38(3): 216-222. doi: 10.18805/ag.v38i03.8981.

  10. Khoirunnisak, A., Prijono, S., Nopriani, S.L., Kemeta, T. and Faisal, M.A. (2024). Improvement of soil pore distribution and soil moisture content using organic matter addition technology. Indian Journal of Agricultural Research. 58(6): 1203-1209. doi: 10.18805/IJARe.AF-842.

  11. Kotorová, D. and Mati, R. (2008). Properties and moisture regime of heavy soils in relation to their cultivation. Cereal Research Communications. 36(5): 1751-1754.

  12. Kotorová, D., Jakubová, J. and Kováč, L. (2013). Soil structure and soil organic matter of selected soil types in different ecosystems. Agriculture (Poľnohospodárstvo). 59(1): 1-8. doi: 10.2478/agri-2013-0001.

  13. Kotov, A.I. and Nerpin, S.V. (1958). Water-resistant properties of clay soils and grounds and the nature of initial filtration gradients. Izvestiya AN SSR. 9: 106-109.

  14. Kumar, S.G., Sridhar, G. and Radhakrishnan, R. (2015). A case study of vacuum consolidation of soft clay deposit. Indian Geotechnical Journal. 45(1): 51-61. doi: 10.1007/ s40098-014-0107-5.

  15. Kumari, U., Swamy, K., Gupta, A., Karri, R.R. and Meikap, B.C. (2021). Global Water Challenge and Future Perspective. In: Green Technologies for the Defluoridation of Water. Elsevier. pp. 197-212. doi: 10.1016/B978-0-323-85768- 0.00002-6.

  16. Lakshmi, R.S., Sreelatha, T., Veerabhadrarao, K. and Venugopalarao, N. (2016). Effect of sugarcane monocropping on soil physical and chemical properties in texturally varied soils. Indian Journal of Agricultural Research. 155-159. doi: 10.18805/asd.v0iof.9622.

  17. Liu, H.H. and Birkholzer, J. (2012). On the relationship between water flux and hydraulic gradient for unsaturated and saturated clay. Journal of Hydrology. 475: 314-327. doi: 10.1016/j.jhydrol.2012.09.057.

  18. Lockhart, N.C. (1986). Electrode Watering of Fine Suspensions. In: Advances in Solid-Liquid Separation. Elsevier. pp. 241-274.

  19. Machikowa, T., Thong-ob, T. and Wonprasaid, S. (2020). Effect of soil moisture and determination of critical soil moisture contents of cassava. Indian Journal of Agricultural Research. 54(4): 483-488. doi: 10.18805/IJARe.A-491.

  20. McNeal, B.L. and Coleman, N.T. (1966). Effect of solution composition on soil hydraulic conductivity. Soil Science Society of America Proceedings. 30: 308-312. doi: 10.2136/sssaj 1966.03615995003000030007x.

  21. Miller, S., Sacchetta, C. and Veal, C. (1997). Electrode Watering of Waste Activated Sludge. In: AWWA 17th Federal Convention, Melbourne. pp. 302-308.

  22. Mitchell, J.K. and Soga, K. (2005). Fundamentals of Soil Behavior (3rd Edn.). John Wiley and Sons, New York. pp. 592.

  23. Nagy, I.V. and Karadi, G.U. (1961). Untersuchungen über den Gültigkeitsbereich des Gesetzes von Darcy. Österreichische Wasserwirtschaft. 13(12): 281-287.

  24. Nerpin, S.V. and Chudnovsky, A.F. (1967). Soil Physics. Nauka, Moscow. pp. 584. (in Russian).

  25. Punia, T. and Singh, A. (2018). Electroosmotic flow and electroosmosis in soil. International Journal of Fauna and Biological Studies. 5(2): 74-77.

  26. Qi, C., Qi, J. and Liu, J. (2018). Stability analysis method for rock slope with an irregular shear plane based on interface model. Advances in Civil Engineering. 2018: 8190908. doi: 10.1155/2018/8190908.

  27. Smollen, M. and Kaffar, A. (1994). Electroosmotically enhanced sludge-dewatering: pilot-plant study. Water Science and Technology. 30(8): 159-168. doi: 10.2166/wst.1994.0402.

  28. Stefanidis, K. (2021). Current Trends, Gaps and Future Prospects in E-Flow Science: Allocating Environmental Water Needs under a Changing World. In: Environmental Water Requirements in Mountainous Areas. [Dimitriou, E. and Papadaki, C. (Eds.)], Elsevier. pp. 201-234. doi: 10.1016/ B978-0-12-819342-6.00008-7.

  29. Tamiru, H., Bussa, B. and Geneti, T. (2023). Assessment of soil and water conservation practices and its implication on soil physicochemical properties at West Guji Zone, Gelana Woreda, Southern Ethiopia. Agricultural Science Digest. 43(3): 346-350. doi: 10.18805/ag.DF-456.

  30. Yin, Y., Cui, Z., Zhang, X., Song, J., Zhang, X., Chen, Y. and Dou, Z. (2024). Effects of clay content on non-linear seepage behaviors in sand-clay porous media based on low- field nuclear magnetic resonance. Water. 16: 883. doi: 10.3390/w16060883.

  31. Zheng, R. and Zhu, H.M. (2017). Experimental study on flocculation- accelerated deposition of dredger fill. Chinese Journal of Geotechnical Engineering. 39(2): 145-148. doi: 10.11779/CJGE2017S2036.

Peculiarities of Hydrodynamic and Electro-osmotic Flow of Water in Clay Soils

S
T
T. Yedoyan2,*
E
E. Neverova-Dziopak3
V
1Faculty of Construction, National University of Architecture and Construction of Armenia, 105 Teryan St., Yerevan 0009, Republic of Armenia.
2Faculty of Management and Technology, National University of Architecture and Construction of Armenia, 105 Teryan St., Yerevan 0009, Republic of Armenia.
3Faculty of Mining Surveying and Environmental Engineering, AGH University of Science and Technology, Kraków, Poland.
4Department of Mathematics, Construction Mechanics and Physics, National University of Architecture and Construction of Armenia, 105 Teryan St., Yerevan 0009, Republic of Armenia.

Background: Global climate change, accompanied by recurrent droughts in arid and semi-arid regions, necessitates the development of more effective approaches for regulating soil water regimes. Conventional hydrodynamic methods often perform poorly in soils with thick clay horizons and low permeability, where electrokinetic processes may provide alternative pathways for enhancing water movement.

Methods: This study evaluated the parameters of hydrodynamic and electro-osmotic flow in irrigated meadow clay soils of Armenia. Laboratory experiments were conducted under controlled hydrostatic pressure gradients and applied electrical field gradients to characterize filtration behaviour and electro-osmotic response.

Result: A strong curvilinear relationship was observed between filtration rate and pressure gradient (R = 0.983), indicating plastic-viscous behaviour of the soil solution. A similar curvilinear dependence was identified between the electro-osmotic coefficient and electrical field gradient (R = 0.884). Electro-osmotic effects substantially increased the water permeability when the electrical field strength exceeded 2.57 V cm-1. These findings demonstrated the potential of electro-osmotic phenomena to enhance water movement in heavy clay soils, offering promising applications for soil-water management under changing climatic conditions.

Efficient management of soil water regimes is critical for sustainable agriculture, irrigation planning and soil amelioration, particularly in heavy-textured, low-permeability soils (Kotorová and Mati, 2008; Lakshmi et al., 2016; Kumari et al., 2021; Tamiru et al., 2023). Conventional hydraulic methods, such as drainage and irrigation systems, often fail to ensure uniform water distribution in clay-rich horizons due to limited water movement (Frenkel et al., 1978; Kotorová et al., 2013; Machikowa et al., 2020; Balaji and Pandiarajan, 2022). This constraint poses serious challenges for crop productivity, irrigation efficiency and reclamation of saline or waterlogged soils, especially in arid and semi-arid regions (McNeal and Coleman, 1966; Kumar et al., 2015; Khoirunnisak et al., 2024). These issues are directly relevant to farmers, irrigation managers and policymakers concerned with agricultural water security. Electrokinetic processes, particularly electro-osmosis, offer a promising alternative by facilitating water movement through soil pores under an external electric field. In contrast to hydraulic flow, electro-osmotic transport operates independently of pressure gradients, enabling water redistribution in clay soils and improving permeability (Lockhart, 1986; Smollen and Kaffar, 1994; Miller et al., 1997; Zheng and Zhu, 2017; Butnan et al., 2024). Although electro-osmotic phenomena are well studied theoretically, their practical use in heavy soils remains limited, underscoring the need for experimental evaluation of water movement under combined hydraulic and electrical influences (Gucal and Khyamyalyaynen, 2020).
       
Traditionally, soil water movement is described by Darcy’s law (Darcy, 1856), which assumes a linear relationship between flow rate and pressure gradient. However, numerous studies have shown that this linearity often breaks down under certain conditions, such as high Reynolds numbers, low flow rates, or when fluids display non-ideal rheological behavior (Kotov and Nerpin, 1958; Nagy and Karadi, 1961; Nerpin and Chudnovsky, 1967). Experimental data indicated that flow frequently initiates only after a threshold pressure gradient is exceeded, underscoring the role of viscosity in flow initiation and behaviour (Gomboš, 2012). Depending on rheological properties, fluids are classified as Newtonian, non-Newtonian, or visco-plastic. These distinctions are particularly relevant for water transport in fine-textured soils, where capillary forces strongly influence hydraulic behaviour (Joshi, 2017; Qi et al., 2018).
       
Electro-osmosis provides an additional mechanism of water movement in which ions in the diffuse layer migrate toward electrodes under an external electric field, carrying water molecules and generating directed pore flow. In clay soils, electro-osmotic transport is not simply proportional to electric field strength, requiring modified models that incorporate soil structure and pore characteristics (Bondarenko, 1973). Comparisons of hydraulic and electro-osmotic mechanisms show that although both depend on pore diameter, electro-osmotic flow declines more slowly with decreasing pore size. Thus, applying direct electric current to clay-rich soils can significantly enhance permeability, offering an effective tool for water management and engineering reclamation practices (Stefanidis, 2021).
               
This study extends existing theoretical research by experimentally investigating water movement in heavy-textured soils under the combined influence of hydrostatic pressure and external electric fields. The objective is to identify the key hydrodynamic parameters governing these processes, thereby providing insights essential for optimizing soil water regulation and improving the efficiency of reclamation and irrigation systems.
Study area and soil sampling
 
The study focused on irrigated wet meadow-gray soils (classified as Anthrosols according to the WRB system) from the Masis region of Armenia. Soil samples were collected from 0-50 cm layer, air-dried, gently crushed and sieved through a 2 mm mesh.
       
Experiments were conducted during 2023-2024 at the Laboratory for Environmental Issues, Conservation of Water Resources and Their Efficient Use of the National University of Architecture and Construction of Armenia.
       
Physical, chemical and physicochemical properties of the soils were analyzed both with and without application of an electric current. Bulk density and total porosity in the soils of 0-50 cm layer were 1,300 kg m-3 and 53%, respectively (Table 1). Wilting point, field capacity and plant-available water were 14.2%, 29.3% and 15.1%, respectively.

Table 1: Physical and water-physical indicators of the (0-50 cm) soil.


       
The soil exhibited a light clay mechanical texture, with 48.6% physical clay, including 19.9% silt (Table 2). The top 0-50 cm soil contained 2.1% humus. Calcium carbonate content was high (21.5%), magnesium carbonate was 2.4% and exchangeable cations totaled 31.5 cmolc kg-1. The exchangeable sodium percentage (ESP) was 4.8%.

Table 2: Physical, chemical and physicochemical properties of the (0-50 cm) soil.


       
The indicators of water-soluble salts in the 0-50 cm soil layer are presented in Table 3. The soil is characterized by a total salt content of 0.196%, an electrical conductivity (EC) of 1.78 dS m-1, a pH of 7.5 and a water-soluble sodium concentration of 2.1 cmolc kg-1. Carbonate (CO32) was not detected.

Table 3: Chemical composition of soil water extract (0-50 cm layer).


 
Experimental setup
 
Laboratory experiments were carried out using a horizontally mounted soil column (1) with a length of 0.35 m (L) and diameter of 0.07 m (d), corresponding to a cross-sectional area of 0.00355 m2 (A). Sand-gravel filters (2) were installed at both ends to prevent mechanical suffusion. A feeding tank (3) at the inlet maintained a constant water level (H1), while filtrate was collected in a glass container (4) at the outlet (Fig 1).

Fig 1: Experimental setup.


       
The column was packed with dry meadow-brown soil compacted to the natural bulk density of 1,300 kg m-3. Two titanium electrodes, anode (5) and cathode (6), were embedded in the soil and connected to a direct current source to induce electroosmotic flow.
 
Calculation of hydraulic and electro-osmotic flow parameters
 
During the experiments, various pressure gradients (I) were applied through the feeding tank (1, Fig 1) and the corresponding filtrate volumes (W) were measured. Each experimental run was conducted in triplicate to ensure reproducibility. The pressure gradient (I), filtration rate (VI) and filtration coefficient (KI) were calculated according to Equations (1)-(3).





 
Where,
(H1 - H2)= The hydraulic pressure across the soil column of  length L.
t= The time at which the filtrate volume is collected.
A= The cross-sectional area of the column.
       
For electro-osmotic filtration, both hydraulic pressure gradients and electric field voltages were applied to the soil column. The electric field strength (E), the electro-osmotic velocity (VE) and the electro-osmotic coefficient (KE) were determined using the Equations (4) to (6):





   
Where,
U= The voltage applied to the electrodes.
QE= The volume of filtrate from electro-osmotic flow, which was further calculated using Equation (7):
 
        QE = QT - QI                  ...(7)
                                                                  
Where,
QT= The combined hydraulic and electro-osmotic flow.
QI = Hydraulic flow alone.
       
In soils with heavy mechanical texture, the relationship between filtration rate and pressure gradient often deviates from Darcy’s Law, exhibiting nonlinear behaviour. To account for this, both linear and nonlinear regression models were applied to the experimental data, expressed by Equations (8) and (9):
 
                  V = A + BI                           ...(8)
 
               V = aI3 + bI2 + cI + d                 ...(9)
 
Where,
A, B and a, b, c, d= Regression parameters obtained from statistical fitting of the experimental results.
               
All experimental data were processed using Origin 6.1 software to calculate mean values, determine regression parameters and verify the consistency and reliability of the results.  
Hydrodynamic movement of water in soil
 
The experimental results of hydraulic filtration in clay soils are summarized in Table 4, with replicate data included for reliability. These results indicate that, although a general proportionality exists between pressure gradient (I) and filtration rate (VI), the relationship is not strictly linear (Liu and Birkholzer, 2012). With increasing I, both V and the apparent filtration coefficient (KI) rise, however, the variation in KI demonstrates that Darcy’s law cannot fully describe water movement in these soils.

Table 4: Relationship between pressure gradient (I), filtration rate (VI) and filtration coefficient (KI).


       
Statistical evaluation of the relationship between V and I, using both linear and nonlinear regression models, is shown in Table 5 and illustrated in Fig 2. For linear regression, the correlation coefficient was R = 0.947, with a standard deviation (SD) of 1.16% and relative error (E) of 12.9%. In contrast, nonlinear regression produced a stronger correlation (R = 0.983), with lower SD (0.63%) and error (7.0%). These results clearly demonstrate that water flow in clay soils is better described by nonlinear rather than linear models, consistent with theoretical expectations for visco-plastic systems and similar to previously reported non-Darcian seepage behaviour in clay-rich porous media (Yin et al., 2024).

Table 5: Statistical indicators of linear and nonlinear regression analysis describing the relationship between filtration rate (VI) and pressure gradient (I).



Fig 2: Linear and nonlinear regression analysis of the relationship between filtration rate (V) and pressure gradient (I).


       
The corresponding regression dependencies are presented in Equations (8) and (9), with parameter values summarized in Table 6. In particular, the rectilinear regression yielded a mean hydraulic conductivity (KI) of 20.43 cm d{ ¹, while nonlinear regression produced a cubic parabolic form that accurately reflects the threshold behavior observed at low gradients. As shown in Equation (9), water movement begins only when I exceeds the limiting gradient I0, after which the flow rate increases rapidly.

Table 6: Numerical values of the parameters of equations (8) and (9).


       
Thus, regression analysis confirms that the filtration process in clay soils deviates from Darcy’s law, with water exhibiting visco-plastic properties. The cubic parabolic regression more accurately describes this nonlinear behavior and can therefore be applied in hydraulic and hydro-ameliorative calculations where precision is required.
 
Combined hydrodynamic and electroosmotic water movement
 
Analysis of the experimental data (Table 7) indicates that the combined action of hydraulic pressure and electric field markedly enhances the filtration rate. The increase is particularly pronounced at higher hydraulic gradients and elevated electric field intensities, highlighting the synergistic effect of electroosmosis in facilitating water movement through clay soils. For example, at I = 0.46 and E = 0 V cm-1, the filtration rate (VI) was 7.48 cm d-1, whereas under the same hydraulic conditions but with E = 2.57 V cm-1, the rate increased to 9.58 cm d-1. When the pressure gradient was raised to I = 0.80 and E = 3.43 V cm-1, the filtration rate reached 31.8 cm d-1.

Table 7: Filtration rates of clay soil under varying hydraulic pressure gradients and electric field intensities (cm d-¹).


       
The nonlinear dependence of total flow on electric field intensity is clearly illustrated in Fig 3. At lower electric field strengths (0-2.57 V cm-¹), the increase in filtration rate ranged from 2.1 to 3.3 cm d-1, corresponding to hydraulic gradients of 0.46-0.80. At higher electric fields (2.57-3.43 V cm-1), the rate of increase was much greater, ranging from 7.5 to 14.8 cm d-1. These results show that electroosmosis exerts disproportionately strong effect once the electric field exceeds a threshold value, a phenomenon also reported in controlled laboratory studies on clayey soils (Jayasekera, 2004).

Fig 3: Dependence of total filtration rate on pressure gradient and electric field voltage in clay soils.


       
To clarify the contribution of electro-osmosis, the ratio of total flow under combined hydraulic and electroosmotic forces (VI,E) to hydraulic flow alone (VI) was analyzed. The resulting relative values are summarized in Table 8 and illustrated in Fig 4. The data clearly demonstrate that the application of an electric field substantially enhances water movement through the soil, as reflected by the increasing VI, E/VI ratios with rising voltage intensity.

Table 8: Relative filtration rates (VI, E/VI) of clay soil under different hydraulic pressure gradients and electric field strengths.



Fig 4: Dependence of the relative flow ratio (VI,E/VI) on electric field strength (E).


       
Statistical analysis confirmed the robustness of this relationship, revealing a strong positive correlation between the relative increase in flow and the applied electric field (R = 0.93, SD = 0.129, N = 25, P<0.0001). This dependency is quantitatively described as:


The fitted regression curve is presented in Fig 4, highlighting the consistent and predictable influence of electroosmosis on soil permeability.
       
These results confirm the theoretical prediction, that electroosmotic effects are particularly significant in clay soils. Because hydraulic flow in capillary decreases proportionally to R4 with diameters, while electroosmotic flow decreases only to , the electric field becomes increasingly effective as particles becomes finer (Mitchell and Soga, 2005).
 
Electroosmotic coefficient (KE)
 
The dependence of KE on electric field intensity is shown in Fig 5, with regression statistics provided in Table 9.

Fig 5: Rectilinear and curvilinear dependence of the electro-osmosis coefficient (KE) on electric field strength (E).



Table 9: Statistical indicators of rectilinear and curvilinear regression analysis between the electro-osmosis coefficient and the gradient of the electric field.


       
Linear regression yielded a correlation coefficient of R = 0.846 and SD = 0.251, whereas nonlinear regression provided a better fit with R = 0.884 and SD = 0.228. The regression equations describing this relationship are:
 
                        KE = 10-5 E (1.374 - 0.590E + 0.096E2)             ...(11)      
 
                             KE = 10-5 (0.173 + 0.419 E)                    ...(12) 
 
Parameter values are summarized in Table 10. Within the experimental range of E = 0.57–3.43 Vcm-1, KE varied between 0.525 × 10-5 and 1.525 × 10-5 cm² s-1 V-1.

Table 10: Numerical values of the parameters of equations (11) and (12).


       
These results agree well with published ranges (4.91 × 10-6 to 1.57 × 10-5 cm2 s-1 V-1) reported by other researchers, confirming both the accuracy of the present experiments and the applicability of electroosmotic enhancement in heavy soils (Punia and Singh, 2018).
 
Implications for soil reclamation
 
The combined hydraulic and electroosmotic experiments demonstrate that the application of a direct electric current can significantly improve soil permeability in clay soils. This enhancement results from the fact that electroosmotic flow is less sensitive to pore size than hydraulic flow. Consequently, in clay soils where Darcy’s law fails and hydraulic movement is severely restricted, electroosmosis provides a powerful mechanism for accelerating water movement.
               
For reclamation of saline-sodic soils in Armenia, these findings are highly relevant. By reducing the time and water required for leaching, electro reclamation methods can improve the efficiency and sustainability of irrigation practices. The experimental confirmation of nonlinear relationships, threshold gradients and strong electroosmotic responses provides a robust scientific basis for implementing these technologies in practice. 
This study demonstrated that water movement in irrigated meadow–gray soils with clay mechanical composition does not fully confirm to Darcy’s law. Filtration rates within a pressure gradient range of 0.024 -0.76 exhibited both rectilinear (R = 0.947) and curvilinear (R = 0.983) dependencies on the applied gradient, with the presence of an initial threshold (I0) confirming the plastic-viscous behaviour of the soil solution.
       
The combined influence of hydrostatic pressure and electric fields revealed that electro-osmotic effects are negligible below 2.57 V cm-1 but become significant at higher field strengths, substantially enhancing soil permeability. Regression-based models established strong correlations between relative flow velocities and applied fields, as well as between the electro-osmosis coefficient (Ke) and electric field intensity (R = 0.84 - 0.88), highlighting the nonlinear character of electro-osmotic transport.
               
These findings emphasize the potential of electro-reclamation as an effective method for improving water movement in clay soils. Beyond providing a theoretical and experimental basis for engineering applications, the results hold practical implications for multiple stakeholders. For researchers, they contribute to refining models of soil-water interactions under coupled hydraulic-electrical fields. For engineers and irrigation planners, they suggest new design strategies for drainage and amelioration systems in clay soils. For policymakers and agricultural managers, they demonstrate the feasibility of incorporating electro-reclamation into water-saving and soil-rehabilitation programs, particularly in arid and semi-arid regions where efficient water use is a national priority.
The authors gratefully acknowledge the Laboratory for the Issues of Environment, Conservation of Water Resources and Their Efficient Use at the National University of Architecture and Construction of Armenia for providing the facilities and technical support required for conducting the laboratory experiments.
 
Disclaimers
 
The views and conclusions expressed in this article are solely those of the authors and do not necessarily represent the views of their affiliated institutions. The authors are responsible for the accuracy and completeness of the information provided, but do not accept any liability for any direct or indirect losses resulting from the use of this content.
The authors declare that there are no conflicts of interest regarding the publication of this article. No funding or sponsorship influenced the design of the study, data collection, analysis, decision to publish, or preparation of the manuscript.

  1. Balaji, K. and Pandiarajan, T. (2022). Performance assessment of soil moisture meter under sandy clay loam soil. Agricultural Reviews. 43(2): 243-248. doi: 10.18805/ag.R-2286.

  2. Bondarenko, N.F. (1973). Physics of Groundwater Movement. Gidrometeoizdat, Leningrad. pp. 215. (in Russian).

  3. Butnan, S., Sriraj, P. and Toomsan, B. (2024). Effects of electrolytes and soil-to-suspension ratios on pH in acidic coarse- textured soil. Indian Journal of Agricultural Research. 58(1): 120-125. doi: 10.18805/IJARe.AF-785.

  4. Darcy, H. (1856). Les Fontaines Publiques de la Ville de Dijon. Victor Dalmont, Paris. pp. 647.

  5. Frenkel, H., Goertzen, J. and Rhoades, J.D. (1978). Effects of clay type and content, exchangeable sodium percentage and electrolyte concentration on clay dispersion and soil hydraulic conductivity. Soil Science Society of America Journal. 42: 32-39.

  6. Gomboš, M. (2012). The Impact of Clay Minerals on Soil Hydrological Processes. In: Hydraulic Conductivity-Issues, Determination and Applications. IntechOpen. pp. 119-149. doi: 10.5772/ 47748.

  7. Gucal, M. and Khyamyalyaynen, V. (2020). Hydrodynamic parameters in various methods of soil electrokinetic treatment. E3S Web of Conferences. 174: 01062. doi: 10.1051/e3sconf/ 202017401062.

  8. Jayasekera, S. (2004). Electroosmotic consolidation and flow in clays. Australian Geomechanics. 39: 77-86.

  9. Joshi, R. (2017). Physical constraints of fine-textured heavy soils and their management-A review. Agricultural Reviews. 38(3): 216-222. doi: 10.18805/ag.v38i03.8981.

  10. Khoirunnisak, A., Prijono, S., Nopriani, S.L., Kemeta, T. and Faisal, M.A. (2024). Improvement of soil pore distribution and soil moisture content using organic matter addition technology. Indian Journal of Agricultural Research. 58(6): 1203-1209. doi: 10.18805/IJARe.AF-842.

  11. Kotorová, D. and Mati, R. (2008). Properties and moisture regime of heavy soils in relation to their cultivation. Cereal Research Communications. 36(5): 1751-1754.

  12. Kotorová, D., Jakubová, J. and Kováč, L. (2013). Soil structure and soil organic matter of selected soil types in different ecosystems. Agriculture (Poľnohospodárstvo). 59(1): 1-8. doi: 10.2478/agri-2013-0001.

  13. Kotov, A.I. and Nerpin, S.V. (1958). Water-resistant properties of clay soils and grounds and the nature of initial filtration gradients. Izvestiya AN SSR. 9: 106-109.

  14. Kumar, S.G., Sridhar, G. and Radhakrishnan, R. (2015). A case study of vacuum consolidation of soft clay deposit. Indian Geotechnical Journal. 45(1): 51-61. doi: 10.1007/ s40098-014-0107-5.

  15. Kumari, U., Swamy, K., Gupta, A., Karri, R.R. and Meikap, B.C. (2021). Global Water Challenge and Future Perspective. In: Green Technologies for the Defluoridation of Water. Elsevier. pp. 197-212. doi: 10.1016/B978-0-323-85768- 0.00002-6.

  16. Lakshmi, R.S., Sreelatha, T., Veerabhadrarao, K. and Venugopalarao, N. (2016). Effect of sugarcane monocropping on soil physical and chemical properties in texturally varied soils. Indian Journal of Agricultural Research. 155-159. doi: 10.18805/asd.v0iof.9622.

  17. Liu, H.H. and Birkholzer, J. (2012). On the relationship between water flux and hydraulic gradient for unsaturated and saturated clay. Journal of Hydrology. 475: 314-327. doi: 10.1016/j.jhydrol.2012.09.057.

  18. Lockhart, N.C. (1986). Electrode Watering of Fine Suspensions. In: Advances in Solid-Liquid Separation. Elsevier. pp. 241-274.

  19. Machikowa, T., Thong-ob, T. and Wonprasaid, S. (2020). Effect of soil moisture and determination of critical soil moisture contents of cassava. Indian Journal of Agricultural Research. 54(4): 483-488. doi: 10.18805/IJARe.A-491.

  20. McNeal, B.L. and Coleman, N.T. (1966). Effect of solution composition on soil hydraulic conductivity. Soil Science Society of America Proceedings. 30: 308-312. doi: 10.2136/sssaj 1966.03615995003000030007x.

  21. Miller, S., Sacchetta, C. and Veal, C. (1997). Electrode Watering of Waste Activated Sludge. In: AWWA 17th Federal Convention, Melbourne. pp. 302-308.

  22. Mitchell, J.K. and Soga, K. (2005). Fundamentals of Soil Behavior (3rd Edn.). John Wiley and Sons, New York. pp. 592.

  23. Nagy, I.V. and Karadi, G.U. (1961). Untersuchungen über den Gültigkeitsbereich des Gesetzes von Darcy. Österreichische Wasserwirtschaft. 13(12): 281-287.

  24. Nerpin, S.V. and Chudnovsky, A.F. (1967). Soil Physics. Nauka, Moscow. pp. 584. (in Russian).

  25. Punia, T. and Singh, A. (2018). Electroosmotic flow and electroosmosis in soil. International Journal of Fauna and Biological Studies. 5(2): 74-77.

  26. Qi, C., Qi, J. and Liu, J. (2018). Stability analysis method for rock slope with an irregular shear plane based on interface model. Advances in Civil Engineering. 2018: 8190908. doi: 10.1155/2018/8190908.

  27. Smollen, M. and Kaffar, A. (1994). Electroosmotically enhanced sludge-dewatering: pilot-plant study. Water Science and Technology. 30(8): 159-168. doi: 10.2166/wst.1994.0402.

  28. Stefanidis, K. (2021). Current Trends, Gaps and Future Prospects in E-Flow Science: Allocating Environmental Water Needs under a Changing World. In: Environmental Water Requirements in Mountainous Areas. [Dimitriou, E. and Papadaki, C. (Eds.)], Elsevier. pp. 201-234. doi: 10.1016/ B978-0-12-819342-6.00008-7.

  29. Tamiru, H., Bussa, B. and Geneti, T. (2023). Assessment of soil and water conservation practices and its implication on soil physicochemical properties at West Guji Zone, Gelana Woreda, Southern Ethiopia. Agricultural Science Digest. 43(3): 346-350. doi: 10.18805/ag.DF-456.

  30. Yin, Y., Cui, Z., Zhang, X., Song, J., Zhang, X., Chen, Y. and Dou, Z. (2024). Effects of clay content on non-linear seepage behaviors in sand-clay porous media based on low- field nuclear magnetic resonance. Water. 16: 883. doi: 10.3390/w16060883.

  31. Zheng, R. and Zhu, H.M. (2017). Experimental study on flocculation- accelerated deposition of dredger fill. Chinese Journal of Geotechnical Engineering. 39(2): 145-148. doi: 10.11779/CJGE2017S2036.
In this Article
Published In
Indian Journal of Agricultural Research

Editorial Board

View all (0)