Electrical resistivity imaging for monitoring soil water motion patterns under different drip irrigation scenarios

The use of hydrogeophysical methods provides insights for supporting optimal irrigation design and management. In the present study, the electrical resistivity imaging (ERI) was applied for monitoring the soil water motion patterns resulting from the adoption of water deficit scenarios in a micro-irrigated orange orchard (Eastern Sicily, Italy). The relationship of ERI with independent ancillary data of soil water content (SWC), plant transpiration ( T ) and in situ measurements of hydraulic conductivity at saturation ( K s , i.e., using the falling head method, FH) was evaluated. The soil water motion patterns and the maximum wet depths in the soil profile identified by ERI were quite dependent on SWC ( R 2 = 0.79 and 0.82, respectively). Moreover, ERI was able to detect T in the severe deficit irrigation treatment (electrical resistivity increases of about 20%), whereas this phenomenon was masked at higher SWC conditions. K s rates derived from ERI and FH approaches revealed different patterns and magnitudes among the irrigation treatments, as consequence of their different measurement scales and the methodological specificity. Finally, ERI has been proved suitable for identifying the soil wetting/drying patterns and the geometrical characteristics of wet bulbs, which represent some of the most influential variables for the optimal design and management of micro-irrigation systems.


Introduction
In semi-arid and arid regions, irrigation management practices depend on the accurate characterization of temporal and spatial soil water content (SWC) dynamics (Vereecken et al. 2008). In general, the soil texture and the soil hydraulic characteristics represent the two main drivers of SWC changes and soil water infiltration (Campbell and Norman 1998). Nowadays, the most common methods adopted to measure SWC distribution at the root-zone level (e.g., time domain and/or frequency domain sensors, neutron probes) present several limitations (Robinson et al. 2008). The main drawback of these SWC methods concerns their pointbased nature and specifically the non-representativeness and restricted sampling volume (i.e., 10-100 cm 3 ). In addition, these measurements strongly depend on sensors location in the soil profile (e.g., Bogena et al. 2015;Robinson et al. 2008). Furthermore, their installation requires drilling that may cause disturbance of soil structure and water flow regime. However, even if the SWC variation is recognized as a proxy of soil water infiltration (Brindt et al. 2019), identifying water flow pathways due to irrigation in the unsaturated soil profile where plant roots are mainly distributed is a complex challenge Zhang et al. 2017). In fact, very large variations in SWC exist throughout the root zone especially under localized irrigation (García-Tejera et al. 2017), mainly due to the preferential root growth in the wetted bulbs (Klepper 1991). Several methods have been developed and applied for the infiltration rates characterization. Specifically, the soil hydraulic conductivity at saturation (K s ) can be obtained by field and laboratory methods. In situ K s measurements, based on constant head (CH) or falling head (FH) methods, are generally preferred instead of laboratory determinations. In fact, the small-sized soil samples handled in laboratory may not be representative of the field conditions (Ibrahim and Aliyu 2016). The high spatial and temporal variability of K s requires, however, a huge number of field measurements for a comprehensive understanding of the infiltration process. Moreover, it is known that the different in situ K s methods provide distinct infiltration rates even under the same soil conditions .
The use of hydrogeophysical methods may contribute to solve some of the above-mentioned limitations for capturing lateral SWC changes and potentially for determining the relative soil water motion. It has been demonstrated how especially electrical resistivity (ER) methods, such as electrical resistivity tomography (ERT) or electrical resistivity imaging (ERI), can successfully image the SWC dynamics, as a derived soil property (e.g., Bertermann and Schwarz 2018;Binley et al. 2015;Robinson et al. 2012;Schwartz et al. 2008). Moreover, ER surveys involve the use of several electrodes in galvanic contact with the soil surface and/or subsurface, which allows significant flexibility for the configuration of ER acquisitions at the field scale . Even though, the most suitable ER configuration needs to be tailored as function of the target problem, and in terms of sensitivity (e.g., horizontal and vertical scales) that can be practically achieved (Ronczka et al. 2015). In this sense, recent ER applications have been conducted in micro-irrigated field conditions for evaluating the wetted fraction area and the irrigation fronts (e.g., Cassiani et al. 2015;Hardie et al. 2018;Moreno et al. 2015;Vanella et al. , 2019. Moreover, the rapid development of timelapse ER measurement systems permits to explore the time domain of eco-hydrological dynamics processes with high accuracy (e.g., Jayawickreme et al. 2014;Singha et al. 2015;Williams et al. 2017). However, the use of these approaches for determining the K s have not yet been exploited, despite it may contribute for better defining the geometry of the wet bulbs, which represents one of the most influential variables for the design and management of drip irrigation systems (Arbat et al. 2013).
In this study, time-lapse ERI surveys were performed to identify the drying/wetting patters of the unsaturated soil profile following the application of deficit irrigation (DI) strategies integrated with drip irrigation systems. ERI surveys have provided important insights to describe the characteristics of the subsoil in terms of SWC, exploring the ER sensitivity at the different levels of SWC. Therefore, the aims of the study were (1) to identify and characterize the water flow paths occurring during micro-irrigation events, (2) to individuate the mechanisms of mass exchange in the subsoil (i.e., wetting and drying patterns) due to irrigation and plant transpiration (T), and (3) to assess the potentiality of ERI for the soil water motion identification. The general objectives have been achieved by: (a) identifying ER distribution as function of the different SWC conditions during the irrigation season; (b) tracking and assessing the soil drying and wetting patterns through the ER changes during irrigation events supplied at different crop evapotranspiration (ET c ) rates (100, 75 and 50% of ET c ); and (c) evaluating the potential use of ERI in determining K s by assessing the temporal evolution of the wet bulbs under different irrigation scenarios (full and DI).

Irrigation setup and ancillary data
The experimental trial was conducted in an orange orchard planted with 12-year-old trees (Citrus sinensis (L.) Osbeck) cv 'Tarocco Sciara' C1882 grafted on Carrizo citrange rootstock [Poncirus trifoliata (L.) Raf. × C. sinensis (L.) Osbeck], located in Southern Italy (Sicily) and managed by CREA-OFA (Centro di Ricerca Olivicoltura, Frutticoltura e Agrumicoltura of the Italian Council for Agricultural Research and Agricultural Economics Analyses). This area presents semi-arid climate conditions, with mean air temperature and relative humidity values of 18.2 °C and 77%, respectively; and annual rainfall and reference evapotranspiration (ET 0 ) values of 1300.7 and 656.6 mm, respectively, for the year 2019. At the orchard under study, trees were trained with a rounded shape with plant density of 380 trees per hectare for 4 × 6 m spacing. The effective rooting depth ranged between 30 and 40 cm (Cassiani et al. 2015;. The effects of DI strategies on the soil drying/wetting patterns of the subsoil were compared with a control treatment (T1), which receives 100% of ET c by two surface drip lines, located directly close to the trunk, each characterized by a flow rate of 4 L h −1 per single-emitter (spaced 0.6 m on each drip line) with a total number of 12 emitters per tree (48 L h −1 ). The FAO-56 approach (Allen et al. 1998) was implemented for estimating the daily ET c . Reference ET 0 , estimates using the Penman Monteith equation, was multiplied by the FAO-56 K c for citrus adjusted by a localized coefficient (i.e., 0.7) calculated as ratio between the canopy cover size (m 2 ) and the area pertaining to each tree (24 m 2 ) (Consoli et al. 2006(Consoli et al. , 2017Fereres et al. 1981;Longo-Minnolo et al. 2020;Saitta et al. 2020). ET 0 values were obtained using hourly and daily meteorological data provided by a weather station located 2 km far from the field site and managed by the Servizio Informativo Agrometeorologico Siciliano (SIAS).
The experimental design consists in a randomized block scheme with 4 irrigation treatments replicate in 3 times; each replica contains 24 trees (e.g., details are given in Consoli et al. 2014Consoli et al. , 2017Puglisi et al. 2019). The adopted DI strategies have the following characteristics: -Sustained deficit irrigation (T2), receiving 75% of ET c by two dripper lines buried at 0.35 m from the soil surface, characterized by a flow rate of the single-emitter of 4 and 2 L h −1 (spaced 0.6 m on each drip line), emitting a total of 36 L h −1 . This strategy allowed to neglect the evaporation losses (e.g., about 25%, according to Consoli and Papa 2013), being the irrigation doses similar to those provided in T1; -Regulated deficit irrigation (T3), receiving an irrigation amount of 100% ET c until DOY (day-of-the-year) 217 in 2019, and 50% ET c . from DOY 218 till the end of the of the irrigation season, emitting a total of 36 or 24 L h −1 per tree. The irrigation systems specifics are the same as in T1; -Partial root-zone drying (T4), receiving 50% of ET c by two surface drip lines located 0.35 m from the trunk line, each characterized by a flow rate of the single-emitter of 4 L h −1 (spaced 0.6 m on each drip line); the drip lines are activated alternatively every week, thus wetting only half of the orange tree root system, emitting a total of 24 L h −1 per tree.
The soil texture at the field site is sandy loam, with 65, 12 and 23% of sand, clay and silt, respectively, and average bulk density of 1.25 g cm −3 (Aiello et al. 2014;D'Emilio et al 2018). SWC for field capacity (FC, log of the pressure in hPa, pF = 2) and wilting point (WP, pF = 4.2) were determined using a sandbox and a pressure plate apparatus as described in Consoli et al. (2017). Ancillary data of SWC (Decagon, Inc., Pullman, WA, USA) and T (Tranzflo NZ Ltd., Palmerston North, NZ) were used for monitoring at hourly scale the soil-water-plant exchanges processes occurring at each irrigation treatment (i.e., full and DI). In particular, 5 SWC probes were installed at the field at 0.3 m below the soil surface, i.e., one for T1-T3 and 2 at both sides of T4 (West and East). Additionally, two trees per treatment (8 in total) were instrumented with a sap flow sensor, located at a height of 0.4 m on the tree trunks, adopting the heat pulse (HP) method and an hoc corrections for wounding effects were applied using 0.48 and 0.33 as fractions of wood and water in the sapwood, respectively (Saitta et al. 2020).
Details on sensors installed at the field site are reported in Mary et al. (2019), , and Vanella et al. ( , 2019.

ERI data acquisition
ERI surveys were carried out during the 2019 irrigation season (DOY, 168-278) using 2 ERI arrays, shown in Fig. 1a, b. The ERI arrays (with length of 10.65 m) covered simultaneously two irrigation treatments (i.e., T1-T2 in Fig. 1a and T3-T4 in Fig. 1b). The ERI arrays consisted of 72 electrodes (stainless steel rods of about 0.15 m, with diameter of 0.03 m) buried for 2/3 of their length into the soil surface with a spacing of 0.15 m. ERI dataset were acquired by a ten-channel Syscal Pro resistivity meter (IRIS Instruments, Orleans, France). The electrode acquisition scheme was a full dipole-dipole skip-2 (with 5,000 direct and reciprocal measurements), because of its inherent strength in solving ER lateral changes (Samoüelian et al. 2005). The high spatial coverage of the adopted ERI configuration permitted to reach depths of investigation of about 1 m (Fig. 1).
The average time for each ERI dataset acquisition was about 25 min (Table 1). A pulse duration of 250 ms for each measurement cycle and a target of 50 mV for potential readings were set as criteria for current injection. The Fig. 1 Electrical resistivity imaging (ERI) arrays at the field site: a refers to acquisitions performed at T1 (full irrigation) and T2 (sustained deficit irrigation); and b at T3 (regulated deficit irrigation) and T4 (partial root-zone drying) ERI monitoring was performed on different temporal scale: seasonal term and short term. The seasonal-term ERI monitoring allowed the evaluation of the background ER (i.e., at t 0 , initial condition, no irrigation) during three periods of the irrigation season (June, July and September; Table 1). The short-term ERI monitoring was performed for assessing the temporal evolution of wet bulbs through the acquisition of ERI dataset, with high temporal resolution (time-lapse mode), specifically during (t 1 -t 5 ) and after (t 6 -t 7 ) the irrigation event. Details on ERI acquisitions duration and irrigation timing are reported in Table 1.

ERI data processing
ERI background and time-lapse dataset were processed with the freeware R 2 code (v3.1) (Binley 2016) to obtain forward/ inverse solution for two-dimensional (2D) current flow in a finite element mesh. As defined by Binley and Kemna (2005) and Binley (2015), the inverse solution is based on a regularized objective function combined with weighted least squares (an 'Occams' type solution). A 2D triangular mesh generated in Gmsh software (Geuzaine and Remacle 2009), consisting of 5,085 elements and 2,621 nodes, was adopted for the ERI background and time-lapse inversions. They were performed at 10 and 5% error level, respectively . The reconstruction of the 2D ERI imagery was performed using ParaView software (v3.8.1).
ER changes (in percentage terms, %) were assessed by running the inversion of the ratio between the ERI dataset referred to selected time periods (e.g., t 1 , …, t 7 , during and after irrigation, Table 1) and the background ERI dataset (t 0 , Table 1), as follows: Table 1 Time schedule of the seasonal and short-term electrical resistivity imaging (ERI) acquisitions carried out at the field site during irrigation timing (t) for the different day-of-the-year (DOY): t 0 , denotes the initial condition with no irrigation; t 1 -t 5 , refer to the irrigation phase and; t 6 -t 7 refer to after the end of the irrigation Starting and ending times (hh:mm) are given in local time. T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial root-zone drying strategies, respectively where d r is the resistance ratio ( Ω ), d t e d 0 ( Ω ) are the resistance dataset of selected time periods (t 1 -t 7 ) and of the initial condition (t 0 ), and F( ohm ) is the resistance value ( Ω ) obtained by running the forward model for an arbitrarily chosen ER (i.e., 100 Ωm). This procedure allowed to identifying the ER changes (%) compared to the initial condition (t 0 , no irrigation) and, thus, to evidence wetting or drying soil patterns (e.g., the threshold corresponding to a decline/increase in ER was set equal to or greater than 10%). ER changes (%) were mainly related to variations in SWC occurring in T1-T4, assuming that further variables including soil temperature, salinity, and composition and arrangement of soil particles, vary minimally during the ERI short-term monitoring (Samoüelian et al. 2005).

ERI-based K s rates
The soil water motion rate derived by ERI imagery (K s, ERI ) was calculated as the ratio between the maximum ERI-based wetting depth (d ERI ) and the time between two consecutive instants within an irrigation event (t 1 -t 5 , Table 1), as follows: where, Δd ERI (cm) is the difference between the maximum depths reached by the wet bulb at time t i and t i-1 , and Δt (s) is the difference between t i and t i−1 during the irrigation event (t 1 -t 5 ).

Hydraulic conductivity at saturation measurements
The falling head method (FH) was used to measure the soil saturated hydraulic conductivity (K s,FH ) at the irrigation treatments (NAVFAC 1986). The K s,FH values were retrieved according to the procedure described by Caselles-Osorio and García (2006) and Pedescoll et al. (2009), consisting in the measurement of the travel time of a water column that moved vertically along an impervious permeameter driven into the soil.
The K s,FH apparatus was a steel tube with height equals to 0.65 m and an internal diameter of 0.10 m. To supply a water pulse mode through the tube, a ball valve was connected with another 0.65 m polyethylene terephthalate tube with a capacity of 6.6 L (Fig. 2). At each treatment, the tube was placed at a distance of 0.5 m from the tree trunk, in a hole pre-drilled into the soil surface (0.3 m depth); and then it was filled with water. A pressure probe (STS-Sensor Technik Sirnach, AG) was inserted inside the tube to measure the pressure (or water heights) variation in time during the soil water motion. The pressure probe worked through a data logger system (CR 200-R, Campbell Scientific), connected to a laptop, that recorded pressure data every minute up to the entire duration of each K s,FH measurement (set at 60 min). A total of 5 repetitions per treatment unit (T1-T4) were collected during the 2019 irrigation season. Measurements of K s,FH were performed in the tree row adjacent to the one where ERI surveys were carried out. The relationship between water level into the tube and time is represented by a negative exponential curve, and its slope is related to the K s,FH (m s −1 ), as follows: where, d is the diameter of the tube (m); L is the buried length of the tube (m); t is time (s); h 1 and h 2 are the heights of the water level (m) inside the tube at time 1 and time 2 (s), respectively.
To derive the K s,FH values, the best fit between the observed (h obs ) and modeled (h mod ) heights of the water level was solved using the ordinary last square method, through an iterative non-linear procedure implemented in Excel solver (Frontline Systems, Incline Village, NV), following Eq. 4: where h obs (m) is the height of the water level observed inside the tube at time t during the experiment; h mod (m) is the corresponding modeled water level calculated by inverting Eq. (3).

Statistical analyses
The goodness of the relationships between the mean SWC and ER decreasing (%) observed at T1-T4, as well as those retrieved between the mean SWC and the ERI-derived-wet bulb depths, was identified on the basis of the coefficient of determination (R 2 ). Discrepancies in K s values among T1-T4 were assessed by performing one-way analyses of variance (ANOVA) (for both K s,ERI and K s,FH ) and the treatment mean values were compared each other adopting the Fisher's least significant difference test (LSD) at 0.05 significance level (p ≤ 0.05).

General weather patterns, transpiration, irrigation and SWC
During the 2019 irrigation season, the cumulative ET 0 and rainfall values were 643 and 91 mm, respectively (Fig. 3a,  c); whereas, the cumulative irrigation amounts were 317, 238, 237 and 159 mm for T1, T2, T3 and T4, respectively. The daily ET 0 and ET c values observed during the seasonal ERI campaigns in June, July and September, 2019 (Table 1) were 6.79, 6.51, and 3.83 mm day −1 and 3.32, 3.18 and 1.87 mm day −1 , respectively (Fig. 3a). As for daily ET 0 and ET c , a similar decreasing temporal trend throughout the 2019 irrigation season was observed in terms of daily T rates (Fig. 3b). At the irrigation treatment level, daily T rates shown average values of 1.55 (± 0.19), 1.10 (± 0.18), 0.86 (± 0.12) and 1.18 (± 0.18) mm day −1 in T1, T2, T3 and T4, respectively. In general, the SWC conditions during the irrigation period ranged between the FC (0.28 m 3 m −3 ) and the WP (0.14 m 3 m −3 ) values for the soil under study (showing a water-holding capacity of 0.14 m 3 m −3 ) (Fig. 3c). During the seasonal ERI campaigns, hourly SWC values at the initial condition (t 0 ), ranged from 0.17 to 0.25 m 3 m −3 in June, from 0.21 to 0.25 m 3 m −3 in July, and from 0.20 to 0.25 m 3 m −3 in September, showing SWC higher in T1 and T2 (0.24 ± 0.01 m 3 m −3 ) than that in T3 and T4 (0.20 ± 0.02 m 3 m −3 ). Figure 4 shows the background ER tomograms at the initial condition (t 0 , i.e., no irrigation) in T1-T4, within the seasonal-term ERI monitoring (Table 1).

Background ERI images
Mean ER values (Ω m) showed a decreasing trend of about 16% at all the treatments from the beginning (June) to the end (September) of the 2019 irrigation season ( Fig. 5a-d). This ER decreasing pattern was higher in T1 (Fig. 5a) and T2 (Fig. 5b) (19% and 20%, respectively) and lower in T3 (Fig. 5c) and T4 (Fig. 5d) (10% and 14%, respectively). Figures 6, 7 and 8 show the time-lapse inversions of the short-term ERI dataset acquired during (t 1 -t 5 ) and after (t 6 -t 7 ) the irrigation events, compared to the initial condition (t 0 , Fig. 4). The mean ER changes (%) observed in T1-T4 during the short-term ERI campaigns are reported in Fig. 9.   Fig. 3 Daily temporal patterns of a reference (ET 0 ) and crop (ET c ) evapotranspiration rates (mm day −1 ); b transpiration (T) (mm day −1 ), and c soil water content (SWC) conditions (m 3 m −3 ), irrigation and rainfall (mm) at the field site from day-of-the-year (DOY) 162-278 (2019). The arrows indicates the periods of the seasonal electrical resistivity imaging (ERI) monitoring. T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial root-zone drying strategies, respectively; FC and WP stand for field capacity and wilting point, respectively Evident contrasts in ER changes (%) became readily apparent for the different irrigation events during the 2019 irrigation season, representing ER increasing (yellow areas) and decreasing (blue areas) patterns. The ER increasing pattern was detected mainly in the deeper soil layers of T4 in June (Figs. 6, 9d), reaching ER increases of about 20% and involving in average 28% of the entire ERI transect. However, the ER decreasing trends represent the most predominant phenomenon, acting especially in the shallow soil layers due to irrigation (Figs. 6,7,8,9). Generally, irrigation resulted in the formation of individual rounded wet bulbs that became elliptical with time and formed a continuous "wet band" by their overlapping. However, the shape of the wet bulbs identified by ERI (Figs. 6, 7, 8) was dependent on the background SWC (t 0 ) in T1-T4, and on the different flow rates emitted (i.e., 8 L h −1 in T1; 6 L h −1 in T2; 8 L h −1 in T3 until DOY 217 and then 4 L h −1 , 4 L h −1 for T4). Figure 10a, b shows that wet bulbs were identified better when the initial SWC was lower, i.e., in T3 and T4; whereas, their identification was less clear during wetter initial conditions (as in T1 and T2). This behavior was also observed when comparing the decreasing patterns during the seasonal ERI monitoring, being the degree of definition of these decreases less detectable from June to September, 2019. Specifically, a good relationship was Mean values (and standard error) of electrical resistivity (ER) tomograms (Ω m) at background (t 0 ) for T1 (full irrigation), T2 (sustained deficit irrigation), T3 (regulated deficit irrigation) and T4 (partial root-zone drying); n corresponds to the number of cells of each tomogram observed between the mean hourly SWC measured during the ERI acquisitions and the observed average ER % decrease, with R 2 of 0.79 (Fig. 10a). A good relationship was also observed between the mean hourly SWC and the ERI-derived maximum depth of the wet bulbs with R 2 of 0.82 (Fig. 10b).

Emitting rates-wet bulbs relationship
The increasing of the emitter flow rate raised the horizontal radius of the wet bulbs (T1-T3 versus T4). In fact, when the same flow rate was supplied in T1 and T3 (i.e., 8 L h −1 ), the horizontal radius was similar in June and July, reaching nearly half of the distance between the emitters (i.e., 0.6 m) at time t 2 (Figs. 6, 7). Such treatments were, thus, characterized by a continuous horizontal band of SWC along the irrigation lines. On September, when T3 was supplied as T4 (4 L h −1 ), the wetting fronts of T3-T4 showed smaller horizontal radius (Fig. 8).
On the other hand, increasing the emitter flow rate did not raise the vertical radius of the wet bulbs. In fact, it was quite similar and oscillating around 0.3 m depth from the soil surface in T1, T3 and T4. After the time t 3 , a prevalent lateral water movement was observed in the irrigation treatments. A stationary pattern was observed at the end of irrigation events (t 6 -t 7 ), with a decline of magnitude of ER decreasing in average 5% less negative than in t 5 (Fig. 9). No vertical preferential flow patterns were detected. Figure 11 shows the mean (± standard error) K s values (μm s −1 ) obtained from ERI (K s,ERI ) and FH (K s,FH ) approaches. In Fig. 11a, it is inferred that T2 had the lowest K s,ERI value (6.43 μm s −1 ) showing significant differences with T1 (25.03 μm s −1 ). The highest K s,ERI values were found in the surface DI treatments (T3-T4) with mean K s,ERI

Fig. 6
Electrical resistivity (ER) change (%) observed during the irrigation phases and after the irrigation event (t 1 -t 7 ) compared to the initial condition (t 0 ) in June, 2019. T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial rootzone drying strategies, respectively

Fig. 7
Electrical resistivity (ER) change (%) observed during the irrigation phases and after the irrigation event (t 1 -t 7 ) compared to the initial condition (t 0 ) in July, 2019. T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial rootzone drying strategies, respectively values of 36.33 μm s −1 , not showing significant differences among them.
Nevertheless, the magnitudes of K s,ERI and K s,FH were quite different, presenting the FH measurements one order of magnitude higher than the ERI approach.

Discussion
In this study, the ERI technique was applied to identify the drying/wetting patters of the unsaturated soil profile following the application of DI strategies integrated with drip irrigation systems. The 2D ERI surveys provided important insights to describe the characteristics of the subsoil in terms of SWC, exploring the ER sensitivity at the different levels of SWC; this aspect has been scarcely investigated in Fig. 8 Electrical resistivity (ER) change (%) observed during the irrigation phases and after the irrigation event (t 1 -t 7 ) compared to the initial condition (t 0 ) in September, 2019. T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial root-zone drying strategies, respectively

Fig. 9
Mean electrical resistivity (ER) change (%) observed during the irrigation phases (t 1 -t 5 ) and after the irrigation event (shaded area; t 6 -t 7 ) in comparison to the initial condition (t 0 , no irrigation period) for the ERI monitoring period (June, July and September, 2019). T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial root-zone drying strategies, respectively previous studies (Mary et al. 2019). In fact, following the homogeneous structure of the soil and the climatic conditions at the field site (Fig. 3), it was shown that the driving factor explaining the decreasing trend in ER is closely related to SWC changes induced by drip irrigation. This consideration is confirmed by independent SWC measurements that showed a transition pattern from drier soil conditions (higher ER values) to wetter conditions (lower ER values), going from June to September during the 2019 irrigation season. It was also confirmed, comparing the seasonal ERI images, that those treatments that received more water (T1 and T2) showed lower ER values than those subjected to more severe DI strategies (T3 and T4) (Figs. 3, 4). This observation is in line with the results of Cassiani et al. (2015) and Bertermann and Schwarz (2018), which have achieved robust relationships between the ER and the different levels of SWC at the laboratory scale, recognizing the SWC as the most influential soil property. This points out the need of using ancillary data for characterizing the field condition for a realistic interpretation of ERI data .
The SWC deficit, detected by the increase in ER that occurs in the deeper soil layers of T4 especially in June, was attributed to the combination of higher T rates (Fig. 3b) and lower SWC (Fig. 3c) recorded in this treatment in comparison with the other DI treatments. This ER increasing pattern (up to 20%) was not observed in the other DI treatments, since they received more irrigation and, therefore, root-water uptake (RWU) process was partially masked by the higher initial SWC conditions. Similar observations were obtained by other authors (Mares et al. 2016;Mary et al. 2019;, who detected a strictly alignment between the high T rates derived by sap flow measurements and the decline of SWC retrieved by ERI due to the RWU. However, uncertainties in the estimation of T rates cannot be excluded and need to be addressed (Flo et al. 2019;Motisi et al. 2012). In this regard, a meta-analysis carried out by Flo et al. (2019), using data of 16 studies and 21 species, evidenced that the average accuracy deviation for sap flow Fig. 10 Relationships between the mean soil water content (SWC; m 3 m −3 ) and a electrical resistivity (ER) decreasing (%); and b the electrical resistivity imaging (ERI)-derived wet bulbs depths (m). T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial root-zone drying strategies, respectively

Fig. 11
Mean (± standard error) hydraulic conductivity at saturation (K s ) values (μm s −1 ) derived by a electrical resistivity imaging (ERI) (K s,ERI ); and b falling head (FH) method (K s,FH ). T1-T4 refer to full irrigation, sustained deficit irrigation, regulated deficit irrigation and partial root-zone drying strategies, respectively. Different letters indicate K s significant differences among treatments according to Fisher's least significant difference test (LSD) (p ≤ 0.05) HP method was about 14%. Probes installation, wounding, scaling sap flow variability and variation in wood parameters were highlighted as the main uncertainties sources influencing sap flow HP method accuracies (Flo et al. 2019;Intrigliolo and Castel 2007;Intrigliolo et al. 2009;Motisi et al. 2012). On the other hand, ER decreasing patterns obtained by the short-term ERI analysis allowed delineating the wet bulbs in different irrigation scenarios (full irrigation versus DI). Small differences were observed between the vertical extent of the wet bulbs formed as a result of the different flow rates (8-4 L h −1 ), reaching average depths ranging from 0.16 to 0.46 m from the soil surface. This confirmed that the vertical flow component depends mainly on the soil texture (i.e., sandy loam). On the contrary, increasing the flow rate from 4 to 8 L h −1 had a more evident influence on the horizontal component, raising the horizontal radius of the wet bulbs. These findings were in agreement with those reached by Elaiuy et al. (2015), who suggested the use of a wider spacing among drippers when the emitting rates increased from 1.0 to 1.6 L h −1 . Moreover, the absence of preferential vertical flow detected by ERI demonstrates the efficiency of the microrrigation system and the appropriate irrigation scheduling operating at the field (i.e., no percolation flow means that the whole SWC is available for the RWU). In this scenario, ERI may be considered as a valid tool for evaluating the irrigation system emission uniformity, aiming at monitoring the discharged flow rates of several drippers at the same time (Rossi et al. 2013). Furthermore, the study demonstrates the ability of ERI to detect wet bulbs for SWC between 0.18 and 0.27 m 3 m −3 , providing useful information on the DI regime. For example, the SWC observed in T2 did not allow identifying well-defined wet bulbs due to the high initial SWC.
The approaches applied to determine K s in the different irrigation scenarios (K s,ERI and K s,FH ) have revealed different patterns, which also vary in their magnitude, mainly due to the different measurement scales and methodologies adopted. More specifically, K s,ERI presented significantly lower values for those treatments with higher SWC (T1 and T2 vs T3 and T4; Fig. 11), whereas no significant differences were observed in K s,FH . This could be due to the FH method that started measuring from 30 cm above the soil surface, thus neglecting the upper layer where, as identified by the ERI technique, most of the SWC is present, unifying the effect of the irrigation strategies on K s,FH . Thus, the soil sampled with FH may not satisfied the saturated conditions at the base of the HF method, which could explain the higher order of magnitude of K s,FH compared to K s,ERI . Additionally, K s,ERI (Eq. 2) performance mainly depends on the ability of ERI in detecting the wet bulbs which in this study seems to be strictly influenced by the degree of soil saturation. Otherwise, FH approach determines K s,FH from the difference of water heights in the permeameter with time, which can be the result not only of a prevalent vertical flow but also of a horizontal one. Even if no preferential vertical flow were individuated by ERI, it is not possible to exclude the presence of preferential horizontally flow due to the roots distribution in T1-T4, that could explain the higher order of magnitude of K s,FH values compared to K s,ERI . Several authors have related the higher infiltration capacity of irrigated soils to the formation of macro-pores due to roots activity, which can contribute almost 85% of the total infiltration variation (Bronick and Lal 2005;Cameira et al. 2003;Wu et al. 2017). However, some authors have pointed out some limitations in the in situ K s determination, as the necessity of a high number of measurements in order to accurately characterize and represent the wide K s variability in space and time ). In particular, K s,FH estimates depend on the initial SWC, on the height of the water released in the soil and on the duration of the soil water motion process (Alagna et al. 2016). Furthermore, the FH methodology provides point-based K s estimations that may not be representative of the entire system, which is partially solved with the use of the ERI technique. However, despite the limitations above-mentioned, the results of this study indicate the potential use of the FH methodology to evaluate the K s on the least disturbed layers of the soil (depth > 0.3 m) while the ERI surveys proved to be more useful to determine the effect of irrigation treatments on K s estimates in unsaturated soils.

Conclusion
The study herein presented highlighted that hydrogeophysical techniques can play an important role in supporting irrigation management strategies for giving insights on the efficiency of the irrigation systems as well as for characterizing the soil water dynamics.
On the one hand, ERI has permitted to identify the water flow paths under different drip irrigation scenarios by delineating the wet soil bulbs. The resulting ERI-derived wet soil bulbs have reached similar depths in all irrigation treatments, regardless of flow rates, suggesting that the vertical component of the flow is prevalent and depends mainly on the soil structure. In addition, ERI technique has allowed to evaluate the efficiency of the microirrigation through the simultaneous monitoring of the flow supplied by the different drippers, enabling the identification of potential failures in the irrigation system. Furthermore, the absence of preferential vertical flows detected by ERI indicates that the irrigation schedule is appropriate. On the other hand, the identification of subsoil ER patterns has permitted to capture the soil drying effect related to the RWU process, especially in the irrigation strategy characterized by the most severe deficit. Finally, ERI resulted more useful for determining the K s in unsaturated soils than FH methodology, which provided K s information mainly referable to the less disturbed soil layers.