In Table 3, we list the input maps used in order to run Benfratello's method within a GIS framework, which are monthly precipitation, monthly average, maximum and minimum temperature, field capacity and land-cover, needed to modulate the maximum evapotranspiration demand. In order to obtain these maps, we proceeded as follows. Field capacity maps, for different land-cover scenarios, were obtained by multiplying the map of water content at field capacity by the map of the hydrologically active soil layer. The latter was obtained by assigning a rooted depth to each land-cover class, according to Foxx et al. (1984) and Walker (1989). In the nearby Gargano area, despite the deep rooted depth obtained by the referred literature, we preferred to assign a shallow hydrologically active soil layer according to Ancona et al. (2010) (Figure 5). Monthly precipitation maps were generated from measured precipitation data by means of Thiessen–Voronoi tessellation (see, e.g. Figure 6, top-right map). Monthly temperatures available at the micrometeorological stations were interpolated by a linear function of the elevation in order to provide the temperature maps (see, e.g. Figure 6, top-left map).