Analyzing the role of soil and vegetation spatial variability in modelling hydrological processes for irrigation optimization at large scale

The main purpose of this paper was to study the effect of spatial variability of soil hydraulic properties (HP) and vegetation parameters (VP) (e.g., leaf-area index, LAI, and crop coefficient, Kc) on modelling agro-hydrological processes and optimising irrigation volumes at large scale. Based on this analysis, the effect of partly overlooking the spatial variability of soil HP and/or VP inputs was verified on a 140 ha irrigation sector in “Sinistra Ofanto” irrigation system in Apulia Region, Southern Italy. Five soil profiles were excavated and the HP were measured in all the soil horizons. Additionally, measurements of soil HP were taken in the surface soil layer in ninety sites distributed over the whole irrigation sector. All the HP measurements were carried out using tension infiltrometer. Remote sensing applications were used to obtain LAI and Kc using European Space Agency’s (ESA) Sentinel-2 images with 10 m resolutions. First, distributed (on ninety polygons with an average area of about 1.5 ha) optimal irrigation volumes and related deep percolation volumes at a depth of 80 cm, were computed using an agro-hydrological model and accounting for the actual observed variability of soil HP and VP inputs. The sector scale irrigation and deep percolation volumes were obtained by aggregating the distributed irrigation volumes. This was considered as the reference scenario (hereafter DVS—Detailed Variability Scenario). Then, reduced variability scenarios (hereafter RVS—Reduced Variability Scenario) were considered, where the information on the actual spatial variability of the soil HP and VP was gradually overlooked to find the minimum data set needed to still have sector scale irrigation volumes and related deep percolation volumes comparable to those obtained under the DVS. Results showed that overlooking VP (RVS-VP) variability did not significantly change the optimal irrigation volumes and the deep percolation fluxes. By contrast, neglecting the HP variability (RVS-HP) showed significant effects on both the irrigation and percolation volumes compared to the DVS. The main practical finding was that, at least for the area investigated in this study, hydraulic characterization of one soil profile in an area of approximately 30 ha provides sector scale irrigation volumes and percolation fluxes comparable to those obtained under the DVS, thus by accounting for all the observed local variability.