Eco-hydrological modeling of soil wetting pattern dimensions under drip irrigation systems

Reliable information on the horizontal and vertical dimensions of the wetted soil beneath a point source is critical for designing accurate, cost-effective, and efficient surface and subsurface drip irrigation systems. Several factors, including soil properties, initial soil conditions, dripper flow rate, number of drippers, spacing between drippers, irrigation management, plant root characteristics, and evapotranspiration, influence the dimensions and shape of wetting patterns. The objective of this study was to briefly review previous studies, collect the analytical, numerical, and empirical models developed, and evaluate the effectiveness of the most common empirical method for predicting the dimensions of soil wetted around drippers using measured data from field surveys. With this review study, we aim to promote a better understanding of soil water dynamics under point-source drip irrigation systems, help improve soil water dynamics under point-source drip irrigation systems, and identify issues that should be better addressed in future modeling efforts. A drip irrigation system was configured with three different emitters with different capacities (2, 4, and 8 l h-1) in the point source to determine the soil wetting front under the point source. The five most selected empirical equations (Al-Ogaidi, Malek and Peters, Amin and Ekhmaj, Li and Schwartzman and Zur) were statistically analyzed to test the efficiency in sandy loam soil. According to the results of the field investigation, statistical comparisons of the empirical models with the field investigation data were performed using the mean absolute error (MAE), root mean square error (RMSE), Nash–Sutcliffe model efficiency (CE), and coefficients of determination (R2). The advanced simulation of the wetting front was used based on the best accuracy of the selected empirical model. In general, the Li model (MAE, RMSE, EF, and R2 were 0.698 cm, 0.894 cm, 0.970 cm2 cm-2, and 0.970, respectively, for the wetted soil width and 1.800 cm, 1.974 cm, 0.927 cm2 cm-2, and 0.986, for the vertical advance) proved to be the best after statistical analysis with field data.


Introduction
In India, agriculture is the sector that consumes most of the water resources (almost 80%) for environmental reasons [1]. As the population density is increasing daily and water consumption is increasing, the irrigation system needs to be modified to produce more with less water and more food, while the agricultural productivity of soil and water resources is suffering. Microirrigation technology is widely used and applied worldwide [2]. In areas where water resources are becoming scarce and problematic, drip irrigation technology leads to "more yield per drop of water" for more sustainable agriculture [3]. This approach is one of the most effective methods to deliver water and nutrients to the root zone of plants [4,5]. The water savings recorded compared to other surface irrigation methods could be as high as 50% for a drip irrigation system [6]. Table 1 Calculated wetted volume dimensions (Wetted depth and wetted radius, units are m) in the 18 soils from Group 2 after application of 1.65 and 6.6 L of water (at 1.65 l h − 1 ) [19]. Statistical 19 48 Where t is time, z + is vertical distance downward, z-is upward distance, and r is radial distance.
D.K. Vishwakarma et al. subsurface drip irrigation. At a depth of 30 cm in sandy soil, 94% of the water was applied below the emitter. In silty soil, the wetting pattern was nearly spherical. The horizontal and vertical dimensions of the wetting front were measured at 28 and 20 cm, respectively, for sandy loam soil in sandy irrigated palms at irrigation rates of 1, 1.8, and 3.4 l h-1 [21]. The soil wetting patterns studied by Ainechee et al. [22] under irrigation sprinkling showed that the amount of water increases the wetted width and depth. The maximum wetted width and depth were found in sandy soils for all irrigation water volumes, followed by silty loam and clay soils. Siyal and Skaggs [23] found that soil texture has a greater influence on wetting geometry due to its relationship with soil hydraulic conductivity and water retention. Greater horizontal expansion in the wetted dimensions usually occurs in the fine-textured soils or the fine-textured layers of the layered soils. Nafchi et al. [24] studied the geometry and volume of the wetted soil during drip irrigation (Table 2). For soils with fine texture, the average diameter of the wetted soil volume was larger than for soils with coarse texture. The average diameter of wetted soil volume for fine textured soils was greater than for coarse-textured soils.
Yasin et al. [25] have shown that a slight increase in horizontal feed and an almost nonexistent increase in vertical feed enhance the changes in bulk density and soil depth while diagonally. They also showed that the increase in vertical advance is more significant than that in horizontal advance, that it can be continued diagonally by increasing the application rate, and that all horizontal advances, vertical advances, and diagonal advances of the wetting front increase when the same amount of water is applied with the decrease in application rate.
Fan et al. [26] investigated and showed that the geometry of the soil-wetting pattern exhibits slight but significant size changes. Moisture is generally "ellipsoidal" in its wetting pattern and shape and in the volumetric soil water content contour. Soil texture has a significant effect on the characteristics of the wetting patterns. The vertical and horizontal frontal wetting distance and soil moisture volume decrease as the clay concentration of the soil increases.

Effect of the flow rate of dripper on soil moisture distribution
The radius of the wetted surface under spot drip irrigation on Malaysian sandy soil was studied by Ekhamaj et al. [27]. The study showed that soil moisture decreased with depth downward and horizontally from the point source of application because moisture was mainly maintained within 0.2 m of the point source. The wetted soil depth of all application rates was reported to be greater than the maximum wetted soil radius at the surface.
Thabet and Zayani [28] studied the wetting patterns in loamy, sandy soils under point source drip irrigation with two emitter flow rates, i.e., 1.5 and 4 l h-1. The study showed a higher vertical movement of soil moisture at a higher flow rate of 40 cm at 1.5 l h-1 after 6 h of water application. In contrast, it was 52.5 cm after 6 h of water application at 4 l h-1. However, the maximum wetting radius (30 cm) was recorded at the lower flow rate (1.5 l h-1) and the minimum (22 cm) at the higher delivery rate (3 h of application) (4 l h-1). Amer et al. [29] studied the geometric wetting of soil at different flow rates of 2, 4, 8, 16, and 24 l h-1 and different operating hours of irrigation from a point source of drip irrigation for cucurbits and grapes in an Egyptian district. The research study showed that the wetted soil depth almost doubled only after irrigation. However, the width of the wetted soil did not increase significantly. After water redistribution, the soil was almost twice as wide. Higher flow rates increased the width of soil wetting but decreased the wetted soil depth. Kandelous et al. [10] studied the distribution of soil water content between two SDI emitters. The study showed that the higher the flow rate of the emitters, the faster the horizontal wetting front progressed. However, the speed of the vertical wetting front is relatively unremarkable. Molavi et al. [30] recorded the coordinates of the soil-wetted bulb as the flow function of the emitters, time of water application, volumetric volume variation, and saturation of soil hydraulic conductivity. Saxena et al. [31] investigated that the maximum radius of wetted soil at and below the soil surface increased with the flow rate of the emitter. However, the wetted soil depth showed the same trend for Vertisols.

Effect of the flow rate of dripper on soil moisture distribution
The wetted soil's width and depth depend on the irrigation duration [32], and both increase with time. In clayey mixed sand soils, the wetted soil depth was greater than the wetted width. Kyada and Munjarappa [33] reported that wetted width increased at the beginning of irrigation in clayey soils. However, the increase was relatively slow. The change in moisture distribution patterns for soil (wetted width and wetted soil depth) is summarized in Table 3 as a function of irrigation duration in different agroecological situations.

Table 2
Variation of wetted depth and width with soil texture for drip irrigation (Nafchi et al. [24]

The main research aims and constraints of this paper
The purpose of this study is to evaluate and summarize the various models that have been developed to determine the wetting process resulting from drip irrigation in the hope of achieving the specific objectives: (a) delineate between analytical, semi-analytical, numerical, and empirical models for determining the wetting pattern resulting from drip irrigation; (b) improve understanding of the soil moisture distribution pattern as a function of soil hydraulic properties, drip runoff rates, spacing, irrigation amounts and frequency, plant water uptake rates, and root distribution patterns, and provide planners with current information on the subject; (c) a brief compilation of existing models for estimating soil wetting geometry; (d) a compilation of different wetting fronts under different irrigation durations in drip irrigation; (e) an estimate of the wetting pattern in clay soils under a spot drip irrigation; and (f) an evaluation of the accuracy of five commonly used empirical models, namely Al-Ogaidi, Malek and Peters, Amin and Ekhmaj, Li and Schwartzman, and Zur with observed field studies.

Description of various models
Researchers have developed several methods that define infiltration from a point/line source of drip irrigation and use them to predict wetting dimension patterns of the moistened soil (above or below the soil surface) for efficient, effective, economical, and accurate design, operation, and management of drip irrigation systems [35,36]. There are three main classes of models for evaluating soil wetting dimensions, namely (i) analytical, (ii) numerical, and (iii) empirical. Analytical models are developed in the mathematical models with a very closed form of a solution; that is, a mathematical, analytical function can explain how the equation of a system can change the solution. A computer model was created to simulate and replicate the mechanisms of a particular system. Empirical modeling is a collective term for practices that involve observations and experiments to create models. An empirical simulation is a form of simulation in which prototypes are developed using specific concepts. Analytical, empirical, and numerical models can help estimate soil moisture movement in the point source of drip irrigation for a given soil type, environmental setting, and various design parameters. Therefore, they can save time and financial capital compared to field tests. This study mainly aimed to collect the developed empirical analytical and numerical models and computer-coded programs to predict the coordinate of the wetting front. Several empirical equations were collected in the study.

Numerical models
Several models use the basic Richards equations as the governing equations of the water flow equation through porous media in the two-dimensional (2D) form to predict soil water pressure or moisture content distribution in unsaturated soil [37]. The hydraulic conductivity of soil under unsaturated movement conditions is nonlinear [38]. It depends on soil water pressure, converted to soil water content using water retention functions. Sometimes numerical flow models are considered inefficient and impractical for planning purposes. The following equation (1), an adapted form of Richard's equations, can be used to explain axisymmetric flow in a variably saturated, rigid, isotropic porous medium: where, θ = volumetric soil water content (L 3 /L 3 ); h = soil pressure head (L); t = time (T); s = radial coordinate (L); z = vertical coordinate (L); K = unsaturated soil hydraulic conductivity (L/T). Many researchers have studied and solved the above formula both numerically and analytically. In solving this equation, choosing the limits and initial conditions is crucial. Philip [39] and Gardner and Mayhugh [40] presented a new equation for Richard's equation, which they recommended use by solving the hydraulic conductivity (Eq. (2)).
With the help of "WetUp," a user-friendly computer program, the results presented by Philip [39] were used to compute the sizes of soil-wetting patterns [42].

Analytical models
These models are used to predict dimensions of soil wetting under the point source of drip irrigation, usually by solving the governing equation of the water flow equation under certain conditions. Cook et al. [42] developed the Microsoft Windows OS-based software program "WetUp" (v1.5) to perform the visualization of wetting patterns. It is user-friendly, environmentally friendly, and freely available to users. The coded program estimates the dimensions of soil wetting patterns in different soil textures with different physical and hydraulic soil properties for surface or subsurface (buried) point sources (emitters). The software includes a database of preconfigured soil types, dripper flow rates (from 0.503 to 2.7 l h-1), irrigation durations (from 1 to 24 h), initial soil moisture conditions (preset by 3, 6, and 10 m suction heights), and dripper location (above or below ground). Philip's [39] approach simulates runoff from an above-and below-ground point source in WetUp. The approach defines the water travel time based on a quasi-linear analysis of a constant three-dimensional unsaturated water flow. The distance between the buried source and the wetted perimeter at a given time can be calculated with equation (3) using Philip's equation (30) [39].
where r and ∅ denote the polar coordinates of a sphere (s = r sin∅, z = r cos∅), R = αr/2 in which α denotes the reciprocal of the macroscopic capillary length scale proposed by White and Sully [41], and L(x) denotes the dilogarithm defined by White and Sully [41] in Eq. (4): where x is a dummy variable, and the dimensionless time T is given by Eq (5): Where q is the emitter flow rate, t is irrigation time, Δθ is the average change in volumetric soil water content in the wetting zone of soil, and α is the reciprocal of the macroscopic capillary length scale. For a surface source, T is related to the wetted perimeter by equation (44) of Philip [39]: In emerging the theory, we deliberate a source of strength q (L 3 T − 1 ) positioned at (s, z) = (0, 0). The radiated distance in the source plane (z = 0) and the maximum vertical distance (s = 0). Eqs. (3) and (6) successfully solved to get estimates of R for a given set of T values using Maple 7 algorithms that were created in the MATLAB programming language [43]. A binary splitting process was used to reduce the relative difference [(LHS-RHS)/LHS] between the left and right halves of equations (3) and (6) to less than 1 × 10 − 6 , as described in the description of this method.
WetUp has been compared to other empirical and numerical solutions for approximating soil wetting dimensions [44]. Its predictions of the geometry of the soil-wetting pattern were less accurate than the empirical equation of Amin and Ekhmaj [45]. The results of the numerical model developed using Hydrus-2D by Simunek et al. [46]. Cook et al. [47] also reported that WetUp underestimated lateral soil wetting at large applied water volumes, especially for soils with coarse textures.
Emitter flow rate, the total volume of water applied, and Saturated hydraulic conductivity.

Warrick [55]
Total added water volume and wetted soil volume.

4.
Zur Wetted soil width, depth, and depth between maximal lateral front advance and soil surface level.
The total volume of applied water.

6.
Li et al. [57] a. Sandy Soil The total volume of applied water.

7.
Ekhmaj et al. [27] a. For Flow rate and total volume of applied water.

Amin and Ekhmaj [58]
The average change in volumetric water content within the wetted zone, the total volume of applied water, emitter flow rate, and saturated soil hydraulic conductivity. 9.
Acar et al. [59] x 2 Wetted soil volume, maximal lateral wetting front advance, depth between maximal lateral wetting front advance and vertical wetting front advance, and depth between maximal lateral front advance and soil surface level.

10.
Ainechee et al. [22] a. Sandy soil Emitter flow rate, the total volume of water applied, and Saturated hydraulic conductivity.

11.
Malek and Peters Emitter flow rate, saturated hydraulic conductivity, time duration of irrigation, the average change in volumetric water content within the wetted zone, and bulk density of soil. 12.
Molavi et al. [30] a. Physical properties of soil, emitter flow rate, and duration of irrigation. 14.

Empirical model
Several empirical models have been developed based on assumptions and drip conditions [48]. Interpretations from field observations have been combined with linear or nonlinear regression and dimensional analysis to create empirical models that can be tested for accuracy [49,50]. The rationale of the empirical models was expressed to estimate one of the parameters such as soil moisture content, the radius of wetted soil bulb, volume of wetted soil bulb, depth of soil wetting, boundaries and shape of wetted soil volume, which is a function of the total amount of water applied, irrigation measure, soil porosity, and soil physical and hydraulic properties (soil bulk density, percent sand, silt, and clay, average change in soil water content, initial soil moisture content, saturated soil hydraulic conductivity, or steady-state infiltration rate) using regression analysis (linear and nonlinear), dimensional analysis, and artificial neural networks (ANNs) from field investigations [2,51].
These methods are used when there is no knowledge of the hydraulic properties of the soil, and mathematical simplicity is essential. Some simple models have been developed for sizing soil-wetting patterns based on volume balance and flow geometry [52,53]. The empirical models developed in the past for surface and subsurface drip irrigation are listed in Tables 4 and 5 with their input parameters for estimating soil wetting.

Study site
The experiment was conducted on an experimental farm of a high-density apple orchard of Sher-e-Kashmir University of Agricultural Sciences & Technology of Kashmir, Shalimar Campus (74.87o E longitude, 34.14o N latitude and 1606 m above mean sea level) for observation and measurement of wetting front (Fig. 1a).
The farm was laid out with a gravity drip irrigation system with a 2, 4 and 8 l h-1 flow rate. The most important and durable soil property is the texture of the wetting front. In order to compare the results of the empirical models for the simulated dimensions of the wetting patterns with the results of other empirical models, an observation data set under drip irrigation is required. Soil samples were collected from the experimental field at 0-15, 15-30, 30-45, and 45-60 cm soil depths. The particle size distribution was analyzed using a sieve set and a calibrated Bouyoucos hydrometer according to the methodology proposed by Bouyoucos [78] and sieve analysis (Gee and Or, [79]) and is shown in Table 6.
The system has three side channels, one for 2 l h-1, one for 4 l h-1, and one for 8 l h-1 emitter discharge rates; for each side channel, three emitters are used for three replicates, and the average value is taken for each. In addition, each emitter and the side ducts were placed far away from the other emitters to eliminate the possibility of other emitters interfering with the wetting front of each emitter. A cross-section of soil approximately 0.50 m wide and 0.50 m deep was sharply delineated and filed after 24 h of each irrigation and t. The width and depth of the wetted soil front in each section were measured (Fig. 1b).

New observations and models accuracy test
Knowledge of crop water requirements, wetting pattern, wetted bale size, and multiple spray devices is necessary for optimal irrigation practices. The performance of the five most popular empirical models (Al-Ogaidi [66], Malek and Peters [60], Amin and Emitter flow rate, the total volume of water applied, and Saturated hydraulic conductivity.

18.
Al-Ogaidi et al. Emitter flow rate, duration of irrigation, soil bulk density, initial soil moisture content, saturated soil hydraulic conductivity, percentage of sand silt and clay. 20.
Emitter flow rate, duration of irrigation, soil bulk density, initial soil moisture content, saturated soil hydraulic conductivity, percentage of sand silt and clay.
Emitter flow rate, the total volume of water applied, and hydraulic conductivity of soil and emitter installation depth.

5.
Aldhfees et al. [ Emitter flow rate, the total volume of water applied, Saturated hydraulic conductivity, and emitter installation depth. Ekhmaj [58], Li [57], and the model of Schwartzman and Zur [52]) is evaluated using the measured wetting front during field investigation in sandy loam soils. The study's results will be useful in drip irrigation design, especially in selecting lateral and manifold spacing, runoff rate, and application time.

Phull and Babar
An experimental test was conducted to measure the soil wetting front. Based on the data collected during the field investigation, the models based on the calculation of the dimensions of the wetting pattern using selected empirical models were tested. Some statistical criteria were used to evaluate and illustrate the performance of the existing empirical models. Some statistical criteria were used to Where N denotes the total number of observations; Wetting Front Foreasted,i denotes ith simulated/forecast wetting front values; Wetting Front Foreasted,i denotes the mean value of simulated/forecast wetting front data; Wetting Front Measured,i = ith measured wetting front values and Wetting Front Measured,i = denotes the mean value of measured wetting front data.

Model accuracy
Comparing all empirical models, the five most popular empirical models, namely Al-Ogaidi, Malek and Peters, Amin and Ekhmaj, Li and Schwartzman, and Zur were selected and compared with the observed field investigation (model equations in Table 1, S. No. 2, 6, 8,11, and 18). This model was found to be the best in this investigation. The model was considered an ideal estimate if it had a lower RMSE, MAE, and CE value and a higher R2 value [48]. A summary of model performance, including statistical parameters such as MAE, RMSE, CE, and R2, is provided in Tables 7 and 8.  Tables 7 and 8 lists the goodness-of-fit test indices for all treatments compared in the study for each of the selected empirical models. In contrast to the horizontal wetting front, good agreement was found between the simulated and measured horizontal wetting fronts.   As a result of the statistical evaluation criteria, the Li model outperformed the selected empirical models with the lowest MAE, RMSE, and the highest NSE and R2 among the selected empirical models. The horizontal wetting front model, followed by the models of Amin and Ekhmaj, Al-Ogaidi, Malek and Peters, and Schwartzman and Zur, performed the best ( Table 7). The Li model performed better than the other empirical models among the selected models in terms of its statistical criteria, i.e., lowest MAE and RMSE, and has the highest CE and R2 among the selected models, followed by the model of Schwartzman and Zur, the model of Malek and Peters, the model of Amin and Ekhmaj, and the model of Al-Ogaidi for the vertical wetting front ( Table 8).
The efficiency of Li's empirical model was very high compared with all other empirical models (Tables 7 and 8). The empirical model was developed specifically for the clay soil. It contributes to a more accurate prediction of the wetted soil dimensions because the model considers the conditions in the clay soil and adjusts the relationship accordingly.
The main explanation for these results could be the consideration of the physical and hydraulic properties of the soil, such as the bulk and particle density, the percentage of clay, silt, and sand, the saturated hydraulic conductivity, and other optimistic field values such as the initial moisture content of the soil. Fig. 2a-c shows the measured and simulated propagation of the wetting front (horizontal and vertical) under the different discharge rates of the emitter (2, 4, and 8 l h-1). In Fig. 2, the difference between the value predicted by the model and the observed value using the model of Schwartzman and Zur is slightly larger for the initial time and slightly smaller for the wetted soil width. Compared to other empirical models, Schwartzman and Zur predicted the overestimation to a greater extent than other empirical models. For clayey soils, an inverse trend was observed for the depth of vertically wetted soil, and this trend is very similar to the all-emitter discharge rates with the measured vertical advance of a wetting front.
The radius and depth of soil surface wetting are slightly underestimated by the Li model, as shown in Fig. 2. Only minor discrepancies were found between the simulated and observed values, except for the model of Schwartzman and Zur, which underestimates the radius of surface wetting [60]. Considering the averaged distribution of plant roots, a designer could use these empirical relationships to find an appropriate flow rate and spacing of emitters. Then, the horizontal extent of the soil-wetting front could be easily determined.

Relationship between wetted soil width and depth
The width and depth of the wetted zones of the soil play an important role in the design of drip irrigation systems. Thus, for different irrigation durations, the wetted soil's width and the wetted soil's depth were determined for drippers with a capacity of 2, 4, and 8 l h-1. The following wetting curves are obtained from the field tests and the results obtained (Fig. 3a-c). The effect of different duration of irrigation and different flow rates of the emitters on the wetted dimensions using different flow rates of the emitters (i.e., 2, 4, and 8 l h-1) are shown in Fig. 3. The radius and depth of the wetted zone increased in all cases when the flow rate of the emitter increased in all cases.
Tables 9 and 10 summarize the effect of the different types of emitters on the measured wetted soil width or depth. As shown in Table 9, the width of the wetted soil was measured and was 14.65, 16.65, and 20.62 cm at an emitter discharge rate of 2 l h-1; rate: 16.20, 20.25, 23.90 cm at an emitter application rate of 4 l h-1; and 20.00, 24.50, 28.81 cm at an emitter application rate of 8 l h-1 with 30-, 60-, and 120-min application times. Similarly, in Table 10, the wetted soil depth was measured and was 13.10, 16.20, and 20.44 cm at an application rate of 2 l h-1; 15.10, 21.50, and 26.00 cm at an application rate of 4 l h-1; and 19.40, 25.00, and 31.00 cm at an application rate of 8 l h-1 with 30-, 60-, and 120-min application times, respectively.
The width and depth of the wetted soil were predicted using Li models for different irrigation durations (2 l h-1, 4 l h-1, and 8 l h-1) (see Fig. 4). The similar increase in the radius of the wetted soil was initially rapid and then slowed over time with the measured value, as shown in Fig. 4. A similar trend was also observed for 4 l h-1 and 8 l h-1 emitter flow rates.
The predicted moist soil width 10 min after the start of irrigation was found to be 10.41 cm. The highest value after 3 h was 21.82 cm for 2 l h-1 emitter (Fig. 4a). Similarly, the wetted soil depth 10 min after the start of irrigation period was 7.95 cm, the maximum after 3 h, i.e. 24.56 cm, as shown in Fig. 4b.  In contrast, Fig. 4b shows that at a flow rate of 4 l h-1 of the emitter, the width of the wetted soil was 12.56 cm 10 min after the start of irrigation, and 3 h later, the maximum width of the wetted soil was 26.34 cm. Similarly, the depth of the wetted soil was 9.95 cm and 30.75 cm, respectively. From Fig. 4c, it can be seen that at a flow rate of 8 l h-1, after 10 min of irrigation, the width of the wetted soil was found to be 15.16 cm, and 3 h later, the maximum width of the wetted soil was 30.86 cm. Similarly, the depth of the wetted soil 10 min after irrigation was 12.46 cm, and maximum 36.76 cm after 3 h.
The higher the water applied, the more vertical wetting front was found. Similar results were found by Mostaghimi et al. [82]. When the irrigation duration was increased from 1 h to 5 h, and the flow rate was increased from 0.6 l h-1 to 3.0 l h-1, more vertical wetting of the soil was found than horizontal wetting front [83]. However, there was little progress in the vertical front when the irrigation duration and sprinkler flow rate increased. Therefore, irrigation duration needs to be improved to prevent water loss due to   deep percolation, mainly with a long time scale for the higher flow rate and, for drip irrigation, mainly with a short time scale for the slightly lower flow rate. The relationship between the wetted soil front and the elapsed time at different flow rates of the emitters is shown in Fig. 4a-c.
In addition, the result of the models considered is illustrated by plotting the observed versus predicted wet dimensions of the soil for each model with a 1:1 axis (Fig. 5a). To illustrate the model results, a linear line was inserted in Fig. 5a, showing that the performance of the Li model, followed by the Al-Ogaidi, Malek, and Peters models, is reasonable because all points have a uniform distribution on a best-fit line (1:1) of wetted soil width. To illustrate the results of the model, a linear fit line was included in Fig. 5a, showing that the performance of the Li model, followed by the Amin and Ekhmaj and Al-Ogaidi models, is adequate since all points have a uniform distribution on a best-fit line (1:1) in the wetted soil depth. Malek and Peter's model was very influential in determining the advance of the wet front, especially at the wet depth, as shown in Fig. 5b. An empirical model by Schwartzman and Zur showed lower performance in determining wetted width and depth.

Conclusion
The distribution of soil water beneath a point source of drip is a three-dimensional problem. Accurate estimation of the wetted dimensions of the soil under drip irrigation can help designers determine the most efficient flow rates and spacing of emitters. The geometry of wetted soil dimensions depends primarily on the physical properties of the soil, drip irrigation parameters such as spacing of drip emitters, flow rate, number of drip emitters per plant or tree, amount of water applied, irrigation frequency, and plant evapotranspiration, and other factors. All of these factors must be considered when designing drip irrigation systems. Several analytical, numerical, and empirical models have been developed to predict the dimensions of wetted soil under a spot surface or subsurface irrigation system. These models are based on experimental results for different soil types and irrigation flow rates. They have yet to be calibrated and should be validated on other site-specific conditions. Wetting front patterns are generated using either an analytical or a numerical model. However, this method takes more time and is tedious and challenging.
Using the above different models, predicting the geometry of the wetted soil dimensions under certain conditions is beneficial when designing the drip irrigation system. The author selected five empirical models based on the most popular ones: i.e., Al-Ogaidi, Malek and Peters, Amin and Ekhmaj, Li & Schwartzman, and Zur models. It is very environmentally friendly and simple to determine the moisture geometry of the soil under spot drip irrigation.
Although no studies have yet been conducted under semiarid conditions to directly quantify the effects of these irrigation systems on tree water status, tree orchard water use, and fruit yield and quality over many years, future research can be conducted on these aspects. In addition to this study, further improvements can be made through practical integration with field work and additional modeling and irrigation scenario analysis to expand this work further. For better prediction accuracy, a wide range of emitter discharge rates should be analyzed better to predict emitter spacing under different irrigation system operating conditions. It is recommended that field trials be conducted to investigate if and how nitrogen and potassium move with the waterfront in the soil based on the timing of application and successive applications.

Author contribution statement
Dinesh Kumar Vishwakarma and Rohitashw Kumar: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data; Wrote the paper. </p> Arvind Singh Tomar, and Alban Kuriqi: Analyzed and interpreted the data; Contributed reagents, materials, analysis tools, or data; Wrote the paper. </p>

Data availability statement
Data will be made available on request.

Additional information
No additional information is available for this paper.

Funding statement
The financial support was received through ICAR-AICRP on Plasticulture Engineering and Technology-sponsored projects during the experiment.

Declaration of competing interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgment
This study would not be possible without the Lab facility and resources of ICAR AICRP College of Agricultural Engineering and Technology's SKUAST-K Jammu and Kashmir. The authors thank the anonymous reviewers for their insightful comments and the editor, who considerably enhanced the manuscript with their thoughtful critiques and improvement. Alban Kuriqi is grateful for the Foundation for Science and Technology's support through funding UIDB/04625/2020 from the research unit CERIS.

WA
Soil surface wetted area (L 2 ) Q Measured emitter flow rate due to pressure build-up in the soil surrounding emitters (L 3 T − 1 ) I Infiltration rate at the time soil surface wetted area became constant (LT − 1 ) V oz Wetted soil volume (L 3 ) a ½ of the maximal lateral wetting front advance (L) b Depth between maximal lateral wetting front advance and vertical wetting front advance (L) h Depth between maximal lateral front advance and soil surface level (L) W Wetted soil width (L) R Wetted soil radius (L) Z Wetted soil depth (L) K s Saturated hydraulic conductivity (L 1 T − 1 ) V Total volume of applied water (L 3 ) q Emitter flow rate (L 3 T − 1 ) a 1 ,..a 8 Empirical coefficients (dimensionless) WP Any wetted parameter (L) b 1 ,..b 8 Empirical coefficients ( Wetting front for intersection zone in horizontal (L) Z if wetting front for intersection zone in a vertical direction (L) t 1 Intersection time of wetting front in a double source of an emitter (T) t 2 Final time of wetting front in the double source of the emitter (T) R 0 cylindrical saturated zone with radius (L) γ Euler's constant (0.577216) R e Equivalent radius of the cavity surrounding a single emitter as a function of any soil water pressure head at the source (P) P Soil water pressure head at the source Φ 0 Soil matric potential (LMT − 2 ) L Length of ceramic pipe (L) K se Saturated hydraulic conductivity of the ceramic pipe (LT − 1 ) h Pressure head (L) θ Volumetric soil water content, (L 3 L − 3 )