Publications

Capillarity and Evaporation Exacerbated Seepage Losses from Unlined Channels

Steady, essentially two-dimensional, tension-saturated seepage from a flat-bottomed soil channel in an arid environment is studied analytically. Physically, the action of capillary spreading, evaporative drive to the atmosphere both from the near-channel banks and capillary fringe boundary, Darcian resistance of the matrix and gravity are juxtaposed and result in trifurcation of infiltrated water into deep percolation and semi-infinite evaporation "wings" with two hinge points on the soil surface and two dividing streamlines. Mathematically, free boundary problem is solved by conformal mappings and the Polubarinova-Kochina boundary-value problem method. The dependence of deep percolation losses and evaporative return flow on the channel width, conductivity, static capillary rise height, and intensity of evaporation is found.

Dipolic Flows Relevant to Aquifer Storage and Recovery: Strack’s Sink Solution Revisited

Three-dimensional groundwater flow to a lake: an explicit analytical solution

New Analytical Solutions for Phreatic Darcian Flows Over Non-Planar Bedrocks

Optimal design of fibers subject to steady heat conduction

Evaporation-Induced Capillary Siphoning Through Hydraulically Connected Porous Domains: The Vedernikov–Bouwer Model Revisited

Seawater intrusion in the coastal aquifer of Wadi Ham, UAE

A transport model was developed to simulate the seawater intrusion in the aquifer system along the Kalbha and Fujairah coast of the United Arab Emirates. The model was used to simulate the salinity levels of the groundwater of Wadi Ham aquifer and its variation in time and space from January 1994 to March 2005. The area covering the coast of Gulf'of Oman in the study domain was taken as a constant concentration boundary with an average salinity (TDS) value of 35 000 mg/L. The effect of artificial recharge on seawater intrusion was evaluated. The results of the simulation indicated that the seawater intrusion is affected by the dry and wet conditions. During the dry years, the velocity vectors are directed from the Gulf of Oman to the aquifer causing severe intrusion problems. During the wet years when rainfall is relatively high and groundwater recharge is encountered from the ponding area of Wadi Ham dam, the velocity vectors are reversed.

The force of Forchheimer – the contributions of Philipp Forchheimer to groundwater hydrology

Minimizing Evaporation by Optimal Layering of Topsoil: Revisiting Ovsinsky's Smart Mulching‐Tillage Technology Via Gardner‐Warrick's Unsaturated Analytical Model and HYDRUS

Abstract Ovsinsky (1899, https://www.rulit.me/books/novaya‐sistema‐zemledeliya‐read‐193251‐1.html ) suggested and tested a water conserving soil no‐till technology for rain‐snow‐fed field crops in a semiarid environment in southern Russia. We model Ovsynsky's unsaturated flow fragment, in which 1‐D steady evaporation and evapotranspiration through a two‐layered soil from a horizontal static water table to a dry soil surface takes place. Gardner's exponential and algebraic functions are used for the unsaturated hydraulic conductivity‐suction head relations. The vertical evaporation flux depends on the dyads and triads (correspondingly) of the parameters of these functions, for example, the saturated hydraulic conductivity and the sorptive number of the two layers. The flux, as a function of the relative thickness of the upper stratum, is analytically found from the solution of one or two nonlinear equations. This relation can be nonmonotonic and exhibits either a minimum or maximum depending on whether this stratum is coarser or finer than the subjacent stratum fed from a horizontal isobar. HYDRUS‐1D simulations confirm these extrema. This explains the experimental results from the literature on mulching/tillage/soil crusting‐sealing, which can increase, decrease, or have no impact on evaporation from a shallow water table. Alterations of the soil's homogeneity to reduce evaporation losses can improve the hydrological balance of soil profiles.

Nonmonotonic moisture profile as a solution of Richards' equation for soils with conductivity hysteresis

The Saint‐Venant type isoperimetric inequalities for assessing saturated water storage in lacunary shallow perched aquifers

Abstract The Dupuit‐Forchheimer (DF) approximation for unconfined groundwater flows is reduced to 2D Poisson's equation (PE), the right hand side of which involves an intensive evapotranspiration (ET) rate. A shallow water table dips inward from a constant piezometric head boundary (closed curve) such that a “dry gap,” demarcated by another closed curve (unknown front), may emerge. Apriori estimates of the volume of the saturated zone and of the “dry gap” area are important for water resources management in drylands. We use the results from the theory of linear elasticity, torsion of elastic bars, for which PE is solved for the Prandtl function. Using the Poincaré metrics with the Gaussian constant curvature , isoperimetric inequalities are obtained for steady‐state DF flows. Conformal moments and confromal radii of the domains and 2‐D Hardy's type inequalities in domains, modeling the Saint Venant bar torsion problem, are involved in the obtained estimates. The studied boundary value problems (BVPs) are nonlinear for ET rates depending on the depth of the water table. In the DF model, the vadose zone (VZ) is considered as a “distributed sink” (similar to a standard “distributed source,” which models recharge to the water table in humid climates). The analytical VZ is collated with one obtained from a numerical solution to BVP for Richards’ equation in a 3‐D saturated‐unsaturated flow. BVPs for Richards’ equation in cylindrical domains, solved by HYDRUS software, give the pressure head, moisture content, Darcian velocity fields, and streamlines.