A Fast Regression-Based Approach to Map Water Status of Pomegranate Orchards with Sentinel 2 Data

: Midday stem water potential ( Ψ stem ) is an important parameter for monitoring the water status of pomegranate plants and for addressing irrigation management. However, Ψ stem ground surveys are time-consuming and difﬁcult to carry out periodically over vast areas. Remote sensing, speciﬁcally Copernicus Sentinel 2 data (S2), offers a promising alternative. S2 data are appropriate for Ψ stem monitoring due to their geometric, temporal and spectral resolutions. To test this hypothesis, two plots were selected within a pomegranate orchard in southern Italy. A pressure chamber was used to collect Ψ stem measurements on four days in summer 2021. Ground data were compared with the temporally closest S2 images with the aim of testing the effectiveness of remotely sensed imagery in estimating and mapping the Ψ stem of pomegranate plants. Regression models were applied with a limited number of ground observations. Despite limited ground observations, the results showed the promising capability of spectral indices (NDVI, NDRE and NDWI) and S2 bands in estimating (MAE ∼ = 0.10 MPa and NMAE < 10%) Ψ stem readings. To understand the dimensional relationship between S2 geometric resolution and the orchard pattern, predictive models were tested on both native S2 data and on denoised (unmixed) data, revealing that native data are more effective in predicting Ψ stem values.


Introduction
The pomegranate tree (Punica granatum L.) is a deciduous fruit tree native to central Asia that has been cultivated since ancient times; its diffusion in the Mediterranean area dates back to the third-second millennium BCE [1].Pomegranate cultivation is currently increasing in Italy, from barely 60 ha to more than 1400 ha in the last ten years.Apulia is one of the most important Italian regions for pomegranate cultivation, supplying about 34% of the national production [2].
Similarly to other fruits rich in polyphenols and antioxidant compounds, interest in pomegranate fruit and its derived products has been steadily increasing in recent years as a result of an emphasis on the effects of phytochemicals on human health [3][4][5][6].Specifically, it has been demonstrated that polyphenols are beneficial in terms of reducing stroke risk, diabetes incidence, blood pressure, systemic inflammation, cancer and cardiovascular diseases [7,8].Pomegranate fruits are a source of polyphenols and other active compounds with antioxidant and anti-inflammatory properties; has potential effects against metabolic disorders; and exerts cardiovascular protection, neuroprotective activity, hypoglycemic effects, tumor cell proliferation and angiogenesis interference [9][10][11][12].
The content of polyphenol compounds in fruit as well as other relevant traits with respect to fruit utilization are influenced by genotype [13][14][15], the environment and growing techniques [16][17][18][19][20].In warm, arid environments, such as those of southern Italy, the commercial and nutritional attributes of fruit are particularly conditioned by plant water status.This is often evaluated as the midday stem water potential (Ψ stem ) using a pressure bomb, the measurement method of which involves darkening the leaf and preventing its transpiration in order to balance the water status of the leaf and stem [21][22][23].Pomegranates are considered a rustic species capable of withstanding water deficits; hence, in order to save water, irrigation restrictions are desirable, especially during periods in which the fruit is less sensitive to water stress, such as the phase of linear fruit growth [19,[21][22][23].Consequently, monitoring the stem water potential in pomegranate orchards may be helpful for water management.However, Ψ stem surveys are time-consuming and difficult to periodically carry out over large plots.In this context, remote sensing could play an effective role in crop monitoring.The Sentinel-2 mission (S2) of the European Copernicus program provides optical imagery for land monitoring, the geometric and spectral features of which satisfy agricultural requirements [24].A broad range of vegetation and water indices (VIs and Wis, respectively) could be obtained by properly combining native S2 bands [25].Concerning the interpretation of spectral data for a prediction of agronomic parameters, according to the literature, VIs and WIs appear to be the most promising in retrieving crop-related information at the field level [26,27], assessing plant vigor [27], exploring climate change effects [28,29], monitoring and interpreting wildlife diseases [30], estimating crop nitrogen content [31] and production [32], supporting insurance activities for damage estimation [33,34], improving Common Agricultural Policy controls [35,36] and investigating crop water stress [37].Specifically, to remotely estimate stem water potential, several techniques can be applied, such as those involving thermal radiation, visible canopy, nearinfrared and shortwave infrared reflectance [38,39], and canopy radiometric temperature derived from an energy balance model based on airborne data [40].The abovementioned approaches have led to satisfactory results (RMSE = 0.36 MPa; RMSE = 0.11 MPa) but are all based on multiple-year time-series analyses.Time-series analyses are time-consuming and expensive; therefore, they are rarely used by farmers.
To the best of our knowledge, there are no fast, easy and cost-effective approaches to monitoring Ψ stem .Within this framework, the present study specifically intended to test the effectiveness of S2 free imagery in estimating and mapping Ψ stem in a pomegranate orchard located in Apulia, basing both regression model calibration and validation on a limited number of ground measures at varying water delivery levels.However, S2 imagery is provided at 10 and 20 m resolutions, depending on the spectral bands.In conventional orchards with 6 m spacing between rows, S2 resolution could be inadequate due to the inter-row detection, and therefore, the deductions could be inaccurate.Spectral unmixing is a procedure that could be used to extract the potential "pure" signal related to vegetation only, thus excluding the contributions of soil [41].To test the reliability of S2 data, models were developed both with native S2 imagery and "pure" S2 imagery derived from the spectral unmixing procedure.

Study Area and Pilot Plots
The present trial is part of a broader multiyear research project, which involves several aspects of pomegranate cultivation.The area of investigation (AOI) is located in the Apulia region (southeast Italy) and is centered on the organic commercial pomegranate orchard (cv.Acco) of the Durante farm (40 • 12 36 N; 18 • 04 07 E) in Nardò (Lecce province).Within this area, 2 adjacent pilot plots, both about 2.5 ha in size, were selected and labelled as plot N (north) and plot S (south), respectively, according to their relative position (Figure 1).The soil of both plots has a sandy-loam texture (USDA soil classification); the distance between trees is 1.75 m within rows and 6 m between rows; the tree canopy can be assumed to be a continuous strip about 3 m wide.The tree rows were established on raised beds and mulched by a white plastic foil; the trees were trained to transverse Y trellis system.During the trial, the soil between rows was managed by periodic plowings.
were established on raised beds and mulched by a white plastic foil; the trees were trained to transverse Y trellis system.During the trial, the soil between rows was managed by periodic plowings.Starting from late June, pilot plots N and S were drip-irrigated with a time frequency of 2-3 days.The irrigations were completed by 10:00 am.Water release was operated through emitters positioned at 50 cm intervals along a single-line pipe laid parallel to the row under the plastic mulch.In accordance with the broader research project, a specifi irrigation restriction was applied, in summer 2021, during the period of fruit linea growth, to assess its effect on the water status of pomegranate trees: plot S received abou 60% of the standard farm irrigation, about 20 m 3 /ha every 2-3 days, while plot N wa irrigated with the full amount of water on the same days.This period lasted from 28 June to 30 July.
From 18 to 20 July, unusual rainfalls occurred and provided about 62 mm o additional water.In the first half of August 2021, about 32 m 3 /ha of water was supplied every 2-3 days to both plots in accordance with standard farm management.The dates o drip irrigation and rainfall are provided in Table 1.All additional cultural practices wer uniformly applied to the two plots during the entire growing season, according to the farm's standard procedures.Starting from late June, pilot plots N and S were drip-irrigated with a time frequency of 2-3 days.The irrigations were completed by 10:00 am.Water release was operated through emitters positioned at 50 cm intervals along a single-line pipe laid parallel to the row under the plastic mulch.In accordance with the broader research project, a specific irrigation restriction was applied, in summer 2021, during the period of fruit linear growth, to assess its effect on the water status of pomegranate trees: plot S received about 60% of the standard farm irrigation, about 20 m 3 /ha every 2-3 days, while plot N was irrigated with the full amount of water on the same days.This period lasted from 28 June to 30 July.
From 18 to 20 July, unusual rainfalls occurred and provided about 62 mm of additional water.In the first half of August 2021, about 32 m 3 /ha of water was supplied every 2-3 days to both plots in accordance with standard farm management.The dates of drip irrigation and rainfall are provided in Table 1.All additional cultural practices were uniformly applied to the two plots during the entire growing season, according to the farm's standard procedures.

Ground Measurements
Ψ stem was measured using a pressure chamber (model 3005, Soilmoisture Equip.Corp, Santa Barbara, CA, USA).Measurements were taken for 9 trees per plot (each at least 10 m apart) as single replicates by sampling two small mature shoots located in the lower part of the canopy on two tertiary scaffolds on opposite sides of each tree.Before the measurements, the shoots were enclosed in two-layer bags (plastic inside and aluminum outside) for 2 h to reach equilibrium.The number of replicates was adjusted to take measurements within about 1 h (11:30 a.m. to 12:30 p.m.), aiming to keep environmental conditions as stable as possible to avoid changes that affect the readings of the physiological parameter Ψ stem .During the period of irrigation restriction, field surveys were carried out on three dates coinciding with the passage of the Sentinel 2 satellite over the study area: 28 July 2021, 7 August 2021 and 12 August 2021.Outside, that period, the same measurements were taken before the start of irrigation, that is, on 23 June.
The difference between the canopy cover of the trees of the two plots was investigated.Digital pictures of the tree shades were taken at solar noon on 23 June 2021, 28 July 2021, 7 August 2021 and 12 August 2021 and processed with the program GIMP (www.gimp.org,accessed on 24 March 2022), transforming the images into black (vegetation) and white (ground) pixels.The percentage of black pixels was indicative of the ground shaded area, which in turn, is known to be related to the tree leaf area.A GNSS receiver (Leica 1200, Leica Geosystem AG, Heerbrugg, Switzerland) was used to position the ground observations.The survey was conducted in the NRTK (Network Real Time Kinematic) VRS mode using the Puglia Region correction service (http://gps.sit.puglia.it/SpiderWeb/frmIndex.aspx, accessed on 24 March 2022).The average accuracy (3D) was 3 cm [42].
A total of 18 points were finally collected (see Figure 1 for the spatial distribution) for the 4 abovementioned dates.

Satellite Data
A single S2 tile, namely T33TYE, was used.Four S2 cloud-free images were obtained as close as possible to the dates of the field surveys (23 June 2021, 28 July 2021, 7 August 2021 and 12 August 2021).The images were downloaded from the Copernicus Scientific Open Data Hub (https://scihub.copernicus.eu/,accessed on 12 January 2022) as Level 2A products (at-the-ground reflectance calibrated and orthoprojected in WGS84/UTM).The Sentinel 2 acquisition time was set at around 10:30 a.m.local time.The S2 L2A data were provided as 100 km 2 tiles and contained a scene classification layer (SCL) alerting about shadowy, cloudy or faulty pixels.The S2 technical features are reported in Table 2.

Data Processing
The Ψ stem data from plots N and S were compared (i) for a single date and (ii) as aggregated dates to test the significances of both instantaneous differences, which refers to a single date, and differences over the trial period, respectively.A comparison of the group mean values was obtained using a t-test.The significance of differences was tested at p < 0.05.The data were processed by the PAST v. 3.24 software.A point vector file (shapefile format) was then generated from the available GNSS measurements to locate Ψ stem values in the plots at the right positions, making a direct comparison with S2 data possible.
During S2 data processing, the native bands with 20 m GSD were initially oversampled (bilinear resampling) to 10 m.This pixel size, if considered with respect to the abovementioned canopy and orchard structure, appears to introduce a high degree of noise in the pomegranate interpretation.In fact, the spectral signature of a single S2 pixel has to be considered a "mixed" one, resulting from the joint contribution of both pomegranate canopy and inter-row space (mostly bare soil in AOI); a scheme is given in Figure 2.
Consequently, a first issue to consider recovering the potential "pure" signal for the pomegranate canopy only.The hypothesis is that, with no pure signal recovered, all deductions would be inaccurate and, possibly, misleading.This task is somehow different form the ordinary spectral unmixing (SU) problem [41].In fact, ordinary SU relies on a priori knowledge of (i) the type of cover class that is expected to participate in the generation of the "recorded" mixed pixel and (ii) the pure spectrum of the cover classes involved.Ordinarily, SU gives an estimate of the percentage of contribution of each class to the recorded signal.
In this study, we needed to recover the pure signal of pomegranate (V) assuming that both the local percentage of area contribution from the two involved classes (pomegranate canopy and bare soil of inter-rows) and the spectrum of the pure local bare soil (S) at the pixel level are known.The spectral un-mixing problem for this study can, therefore, be described by Equation (1): where F i (x,y) is the local reflectance value of the i-th band from S2 data, S i is the reference reflectance value of the i-th band for the bare soil class, V i (x,y) is the desired estimate of the local pomegranate reflectance of the i-th band; and α(x,y) and β(x,y) are the numerical parameters representing the local partial contribution of S and V to the mixed signal F. They are assumed to vary locally but remain the same for all the bands.Consequently, a first issue to consider recovering the potential "pure" signal for the pomegranate canopy only.The hypothesis is that, with no pure signal recovered, all deductions would be inaccurate and, possibly, misleading.This task is somehow different form the ordinary spectral unmixing (SU) problem [41].In fact, ordinary SU relies on a priori knowledge of (i) the type of cover class that is expected to participate in the generation of the "recorded" mixed pixel and (ii) the pure spectrum of the cover classes involved.Ordinarily, SU gives an estimate of the percentage of contribution of each class to the recorded signal.
In this study, we needed to recover the pure signal of pomegranate (V) assuming that both the local percentage of area contribution from the two involved classes (pomegranate canopy and bare soil of inter-rows) and the spectrum of the pure local bare soil (S) at the pixel level are known.The spectral un-mixing problem for this study can, therefore, be described by Equation (1): where Fi(x,y) is the local reflectance value of the i-th band from S2 data, Si(x,y) is the reference reflectance value of the i-th band for the bare soil class, Vi(x,y) is the desired estimate of the local pomegranate reflectance of the i-th band; and α(x,y) and β(x,y) are the numerical parameters representing the local partial contribution of S and V to the mixed signal F. They are assumed to vary locally but remain the same for all the bands.
Since the unknown parameter in the problem is the pure V(x,y) spectrum, Equation (1) can be inverted according to Equation (2): Since the unknown parameter in the problem is the pure V(x,y) spectrum, Equation (1) can be inverted according to Equation (2): This makes it possible to obtain an estimate of the expected "pure" signal from the pomegranate canopy once the local contributions of V (β) and S (α) are known in terms of the percentage of area within the single S2 pixel.
Based on these premises, a reference spectrum (one for each of the dates considered) of "pure" bare soil was obtained by properly sampling the S2 images available at a location where a bare soil area (AOI2) with a size consistent with S2 geometric resolution was found based on a local inspection.AOI2 was found 1.8 km away from the pilot sites (Figure 3), showing a size of 0.8 ha and being characterized by the same sandy-loam texture (USDA soil classification) of the two pilot plots, which is very common in the area.
Parameters α(x,y) and β(x,y) were locally mapped once the orchard pattern (in terms of the size and direction of the rows) was known and a proper vector graticule was aligned to the S2 tile.The orchard rows were initially represented using a line vector layer interactively edited in QGIS.The rows were interpreted from a high-resolution 2019 aerial orthimage from AGEA (Italian Agency for Payments in Agriculture).The orchard owner confirmed that the orchard pattern did not change between 2019 and 2021.Lines representing the middle of the pomegranate rows were drawn and saved as vector shapefiles.Successively, a buffer of 1.5 m was generated around lines (B) to map the expected area covered by pomegranate canopy (about 3 m wide).
A 10 m squared graticule, perfectly fitting S2 pixels, was then generated (G) and clipped using (B).The resulting layer (C) contains only the parts of the graticule elements potentially filled with pomegranate canopy.For each, the corresponding area was computed using ordinary field calculation tools available in QGIS.The QGIS plugin MMQGIS v. 2021.9.10 was used to operate a spatial join aimed at transferring to (G) the value of the areal contribution of C elements falling within the same (G) feature.A further field computation involving G was finally carried out to derive the local value of β (the fraction of vegetated cover in a pixel) obtained by dividing the previous value by the area of a single S2 pixel (100 m 2 ).Consequently, α was obtained as (1-β).The entire process is graphically shown in Figure 4.
Horticulturae 2022, 8, x FOR PEER REVIEW 7 of 18 Based on these premises, a reference spectrum (one for each of the dates considered) of "pure" bare soil was obtained by properly sampling the S2 images available at a location where a bare soil area (AOI2) with a size consistent with S2 geometric resolution was found based on a local inspection.AOI2 was found 1.8 km away from the pilot sites (Figure 3), showing a size of 0.8 ha and being characterized by the same sandy-loam texture (USDA soil classification) of the two pilot plots, which is very common in the area.Parameters α(x,y) and β(x,y) were locally mapped once the orchard pattern (in terms of the size and direction of the rows) was known and a proper vector graticule was aligned to the S2 tile.The orchard rows were initially represented using a line vector layer interactively edited in QGIS.The rows were interpreted from a high-resolution 2019 aerial orthimage from AGEA (Italian Agency for Payments in Agriculture).The orchard owner confirmed that the orchard pattern did not change between 2019 and 2021.Lines representing the middle of the pomegranate rows were drawn and saved as vector shapefiles.Successively, a buffer of 1.5 m was generated around lines (B) to map the expected area covered by pomegranate canopy (about 3 m wide).
A 10 m squared graticule, perfectly fitting S2 pixels, was then generated (G) and clipped using (B).The resulting layer (C) contains only the parts of the graticule elements potentially filled with pomegranate canopy.For each, the corresponding area was computed using ordinary field calculation tools available in QGIS.The QGIS plugin MMQGIS v. 2021.9.10 was used to operate a spatial join aimed at transferring to (G) the value of the areal contribution of C elements falling within the same (G) feature.A further field computation involving G was finally carried out to derive the local value of β (the fraction of vegetated cover in a pixel) obtained by dividing the previous value by the area of a single S2 pixel (100 m 2 ).Consequently, α was obtained as (1-β).The entire process is graphically shown in Figure 4.The updated G layer was further completed with the local reflectance values from the S2 bands for all 4 dates considered, making it possible to couple spectral information (F) with the areal one (α and β).Given the previously mentioned reference spectra for bare soil, Equation (2) was applied for each G element and the local estimate of reflectance values of the "pure" pomegranate generated, including the ones at the position of the surveyed plants.
For the interpretation of spectra for predicting agronomic parameters (i.e., Ψ stem ), according to the literature, the adoption of appropriate vegetation indices (VIs) appeared to be the most promising one.VIs are in fact able to provide information on crop phenology, growth rate, physiological stresses, etc. [44].
For this study, 2 VIs were selected based on their renowned ability to investigate crop phenology, growth rate and chlorophyll content: (i) the normalized difference vegetation index (NDVI) and (ii) the normalized difference red-edge index (NDRE) [26,[45][46][47].
One additional water index, namely the normalized difference water index (NDWI) [47,48], was considered since water-content monitoring is the main issue of this research.
All of the selected indices, given their formulation (see Table 3), ranged between −1 and +1.In particular, higher NDVI and NDRE values corresponded to highly vegetated/active areas, and high NDWI values indicated high leaf water content or plant water status.
VIs were computed from both native and modified ("pure") spectra, thus generating two datasets of VIs (namely, VI M and VI U ) and NDWI (namely, NDWI M and NDWI U ).
Surveyed plant positions were used to extract the corresponding NDVI, NDRE and NDWI values to characterize the local spectral properties.The updated G layer was further completed with the local reflectance values from the S2 bands for all 4 dates considered, making it possible to couple spectral information (F) with the areal one (α and β).Given the previously mentioned reference spectra for bare soil, Equation (2) was applied for each G element and the local estimate of reflectance values of the "pure" pomegranate generated, including the ones at the position of the surveyed plants.
For the interpretation of spectra for predicting agronomic parameters (i.e., Ψ stem ) according to the literature, the adoption of appropriate vegetation indices (VIs) appeared to be the most promising one.VIs are in fact able to provide information on crop phenology, growth rate, physiological stresses, etc. [44].
For this study, 2 VIs were selected based on their renowned ability to investigate crop phenology, growth rate and chlorophyll content: (i) the normalized difference vegetation index (NDVI) and (ii) the normalized difference red-edge index (NDRE) [26,[45][46][47].
One additional water index, namely the normalized difference water index (NDWI [47,48], was considered since water-content monitoring is the main issue of this research. All of the selected indices, given their formulation (see Table 3), ranged between −1 and +1.In particular, higher NDVI and NDRE values corresponded to highly vegetated/active areas, and high NDWI values indicated high leaf water content or plan water status.

Spectral Index Name Formula
Normalized Difference Vegetation Index NDVI = Spectral Index Name Formula Normalized Difference Vegetation Index NDVI = Normalized Difference Water Index NDWI = To investigate the capability of obtaining spectral information from indices to distinguish ordinarily irrigated plants (plot N) from the ones receiving a reduced water volume (plot S), a t-test and F-test were performed with respect to the mean and variance values of the two groups for both native and pure datasets.A comparison was performed with reference to both the single dates (4 tests) and to the aggregated measures (i.e., 1 test including all observations for the 4 dates).

Relationships between Ground and Satellite Data through Regression Models
The Ψ stem values on surveyed plants were directly compared using regression models with both native and pure spectral indices and single bands.Additionally, using a similar approach proposed in [47], the following ratios were computed for both native and pure bands: The generic τ i ratio was tested against NDWI and against all of the available S2 bands for both the native and pure datasets.The 1st order polynomial of Equation ( 7) was used to model all of the relationships: where â and b are the "best" estimates of model coefficients, and x is the generic predictor (NDWI or single band).
An estimate of Ψ stem can be obtained from the ratios by inverting Equation (7) and by replacing τ i with τi (Equation ( 8)) where Ψstem is the estimated Ψ stem and VI is the generic vegetation index used to build the ratio.
To test the capability of the regression models of predicting Ψ stem , the mean absolute error (MAE, Equation ( 9)) and NMAE (Normalized MAE, Equation ( 10)) were used as quality measures to compare the estimates with the observed values of Ψ stem .
where |x obs − x est | is the absolute difference between the observed water potential (Ψ stem ) and the estimated one ( Ψstem ).
Given the low number of available observations, we decided to use them all to calibrate the predicting models and to approach the validation step using the Leave One Out strategy [49].In this way, all of the observations, alternatively, were used as members of the validation set to compute both MAE and NMAE.

Testing Ψ stem and Ground Shaded Area Differences between Plots S and N
The observed stem water potential values ranged from −0.78 MPa in plot N on 12 August to −1.66 MPa in plot S on 23 June.Ψ stem proved to always be lower in plot S than in plot N. A first analysis aimed to test if any significant difference affects the Ψ stem values from ground measures taken in plots S and N (Table 4).Comparisons at the single-date level proved that, on 23 June and 28 July 2021 (first and second dates), the Ψ stem differences between plots were not significant.Differences of about 0.06 and 0.04 MPa, respectively, were found, corresponding to percentage difference of about 4 and 3%.On 7 and 12 August 2022 (third and fourth dates); on the contrary, the Ψ stem difference between plots was found to be significant.Differences of 0.17 and 0.15 MPa, corresponding to percentage differences of about 14% and 12%, respectively, were found.Considering all of the dates together, the Ψ stem of plots N and S showed a significant difference of 0.11 MPa, corresponding to a percentage difference of about 8%.The portion of ground shaded by pomegranate canopies was different for the 2 plots.As an average of the four dates of measurements, the percentage of shaded soil in plot S exceeded that in plot N by about 17%.

Vegetation Indices Analyses
With reference to spectral measures from the satellite, the values for plots S and N were compared to test if any significant spectral difference could be found.Again, the tests were Similar results were obtained for the third and fourth surveys (7 August and 12 August) although NDWI M and NDWI U showed p-values very close to the significance limit of 0.05 (7 August: 0.07 and 0.06 for NDWI M and NDWI U , respectively; 12 August: 0.06 and 0.05 for NDWI M and NDWI U , respectively).In general, NDWIs showed lower values in plot N than in plot S.
Similarly, the t-test conducted for aggregated measures on NDVIs and NDREs showed no significant differences, while the NDWI M p-value was very close to the significance limit (p = 0.06) and NDWI U was significant (p = 0.03).Comparing these results with the ones in Table 4 confirms the ability of NDWI to detect differences in leaf water content.

Relating Ground to Satellite Data through Regression Models
Relating Ψ stem to spectral bands, VIs and NDWI through regression models led to poor results; in fact, no strong correlation was found for both the "pure" and "mixed" values.However, Ψstem was reasonably estimated when using τi .The results are reported in Tables 6 and 7.It is worth remembering that τi (estimate of τ i ) already summarizes the spectral content of both VIs and bands (see Table 3 and Equation ( 7)).For this reason, B8, B4, B5, NDVI and NDRE were not used as predictors (Equation (7)) of the corresponding τi .Table 6.Statistics concerning the τ 1 and τ 2 estimates by Equation (7).SE = standard error, r = Pearson correlation coefficient, R 2 = coefficient of determination.Predictors are native ("mixed") bands and spectral indices.B8, B4, B5, NDVI and NDRE were not considered (n/a).Table 6 shows the following good correlations: τ 1 vs. B2 M (R 2 = 0.70), τ 1 vs. B3 M (R 2 = 0.78) and τ 1 vs. NDWI M (R 2 = 0.74).Unsatisfying correlations were found between τ 1 and the other bands and spectral indices.
Table 7 reports the parameters and statistics of the model (Equation ( 7)) for τ i , with the "pure" bands and indices.The results demonstrate a general improvement in the estimates while working with de-noised spectral measures.The following R 2 values were, in fact, found for τ 3 , with B2 U , B3 U and NDWI U : 0.77, 0.82 and 0.75, respectively.
Whatever the approach (based on native or de-noised spectral measures), it can be noticed that the correlation of τ i with bands/indices is, generally, negative (r < 0).The only exception comes from NDRE, which once more, demonstrates the synthesis of features of vegetation different from NDVI.Negative correlation values suggest that, generally, plant water status increases when reflectance decreases (or spectral index), thus confirming the absorption capability of water.It is, therefore, reasonable to interpret the higher prediction power of the medium infrared bands as being correlated with their well-known sensitivity to the water content of surfaces.
With reference to the best performing models for τ i prediction (R 2 ≥ 0.7), the corre- sponding Ψ stem estimates were generated.The associated MAE and NMAE values were therefore computed (see Table 8) to test their prediction capability.
Table 8 shows that the Ψstem values obtained using τ 1 (B3 M ) and τ 2 (B11 M ) as predic- tors are the most accurate ones.Conversely, τ 3 and τ 4 , although showing better correlations with bands and indices (compare Tables 6 and 7), led to less accurate results.This means that native data appear to perform better when predicting local water potential conditions than de-noised ones.Although it remains to be confirmed with other more representative tests, this finding suggests that the joint effects of the soil-plant system and its related spectral response are probably more representative of the local water-related dynamics.
From an operational point of view, these findings can be translated in maps able to give an estimate of the local expected stem water potential.To make this evident for the case study, a map of Ψstem was generated for the four sample dates using the following two models: one based on τ 1 (B3 M ) and the other ones based on τ 2 (B11 M ).For each date, the maps of differences between the two approaches were also generated (Figure 5i-l).They show no negligible differences in reaching their maximum (about ±50%) during the growing season (up until 12 August).The value of differences cannot entirely be associated with the joint contribution of the estimated uncertainty of each of the maps compared (expected value is around 14% according to the Variance Propagation Law [50]).Consequently, the resulting differences highlight a low degree of robustness of the estimates, making them highly dependent on the predictor adopted.

Discussion
Midday stem water potential is well-known as a sensitive indicator of plant water status [23,51].In previous studies, the Ψstem of non-water-stressed pomegranate trees was found to vary from −0.73 to −0.98 MPa [52] or to decrease from −0.6 to −1.5 MPa midseason, when vapor pressure deficits increased [19].In the same studies, the Ψstem of stressed trees dropped to −1.98 MPa at the end of a water-withholding period [52], to −2.0 MPa when irrigation was managed with a sustained deficit or to −2.4 MPa when irrigation was managed with a regulated deficit [19].In the present study, Ψstem ranged from −0.78 MPa (plot N, 12 August) to −1.66 MPa (plot S, 23 June); these values, compared with the ones reported above, indicate a light level of water deficit for trees of both plots.This consideration is consistent with what was reported in a recent review on the effect of water deficits on pomegranate trees.According to this report, Ψstem values between −0.98 and −1.23 MPa indicate no/minimal water deficit, while values between −1.5 and −1.8 MPa indicate a mild/moderate level of water deficit, and values from −1.9 to −2.8 MPa account for a severe/strong water deficit.[20].Regarding the single-date data, the small difference between Ψstem of plots N and S found on 23 June and 28 July (i.e., before and at the end of the period in which the amount of irrigation in plot S was restricted) can be explained by  (a-d) and τ 2 (B11 M ) (e-h).Maps of percentage differences are also reported (i-l).Maps are given for the sampled dates: 23 June (a,e,i), 28 July (b,f,j), 7 August (c,g,k) and 12 August (d,h,l).Red dots represent the spatial distribution of surveyed pomegranate plants.The reference system is WGS 84 / UTM 33N, EPSG: 32633.

Discussion
Midday stem water potential is well-known as a sensitive indicator of plant water status [23,51].In previous studies, the Ψ stem of non-water-stressed pomegranate trees was found to vary from −0.73 to −0.98 MPa [52] or to decrease from −0.6 to −1.5 MPa mid-season, when vapor pressure deficits increased [19].In the same studies, the Ψ stem of stressed trees dropped to −1.98 MPa at the end of a water-withholding period [52], to −2.0 MPa when irrigation was managed with a sustained deficit or to −2.4 MPa when irrigation was managed with a regulated deficit [19].In the present study, Ψ stem ranged from −0.78 MPa (plot N, 12 August) to −1.66 MPa (plot S, 23 June); these values, compared with the ones reported above, indicate a light level of water deficit for trees of both plots.This consideration is consistent with what was reported in a recent review on the effect of water deficits on pomegranate trees.According to this report, Ψ stem values between −0.98 and −1.23 MPa indicate no/minimal water deficit, while values between −1.5 and −1.8 MPa indicate a mild/moderate level of water deficit, and values from −1.9 to −2.8 MPa account for a severe/strong water deficit [20].Regarding the single-date data, the small difference between Ψ stem of plots N and S found on 23 June and 28 July (i.e., before and at the end of the period in which the amount of irrigation in plot S was restricted) can be explained by the complete absence of irrigation for the first date and by the rainfall that occurred in the second half of July, flattening the tree water status for the second date.On the third (7 August) and fourth (12 August) dates of measurement, i.e., after the end of the differentiation of irrigation between the two plots, the lower Ψ stem values of plot S indicated a lower water status, which can be explained by the greater leaf surface attributable to these trees based on their canopy shadow.
As far as estimating the capability of these spectral measures is concerned, it was found that spectral signals from the green (B3) and the medium infrared (B11) bands are the best performing ones, making the generation of Ψ stem estimates with an MAE and an NMAE of about 0.11 MPa and 9%, respectively, possible.These results confirm the association of NDWIs and short-wave infrared (SWIR) bands with tree water status, as found in other studies [53,54].The approach followed in this test provides sufficiently reliable results despite having a limited amount of data (18 ground measurements for four dates), consistent with the objective of the present study: to develop an easy, fast and cost-effective methodology to map Ψ stem over the orchard.Other studies tried to estimate Ψ stem using a multi-seasonal approach, obtaining similar results using commercial high-resolution imagery (Planet and WorldView-2) [38,39], which are not compliant with the objective of the present study.
S2 imagery is provided at 10 (visible and infrared) and 20 m (red-edge and SWIR) resolution depending on the band.S2 geometric resolution could be inadequate when monitoring orchards with 6 m spacing between rows.To overcome this issue, spectral unmixing was performed in order to extract a "pure" signal from the canopy while excluding the bare soil contribution.Good model performances were achieved when native S2 data were used and slightly worse when de-noised measures were used.Although to be confirmed with other more representative tests, this finding suggests that the joint effects of the soil-plant system and its related spectral responses is probably more representative of the local water-related dynamics.If confirmed, this issue will strengthen, once more, the role of S2 L2A data even for describing the behavior of vegetated discontinued surface as orchards despite the expected noise from the inter-row areas.Both approaches (based on native and "pure" spectral bands) led to negative correlations between τ i and various bands/indices, suggesting that reflectance decreases (or spectral indices) when plant water status increases, thus confirming the absorption capability of water.Only NDRE led to positive correlations; therefore, this study confirms that NDRE expresses different vegetation features from NDVI or NDWI.
Regarding the degree of confidence of the Ψ stem estimates, it was proven to be unsatisfactory for operational purposes, since it is strictly dependent on the predictors selected: τ i (i.e., spectral index) and spectral bands used to model τ i (namely, B3 and B11).Nevertheless, some interesting indications come from our tests.The first one concerns the relationship between the spectral index involved in the τ i computation and the band used as a predictor.In fact, the results suggest that an NDRE-based τ can be better estimated from the mid-infrared bands; conversely, an NDVI-based τ can be better estimated form visible bands (specifically, the green one).
Conclusions: (i) Ψ stem can be satisfactorily estimated by the Ψstem equation, even with a limited amount of ground data; (ii) native data, despite noise, appear to perform better when predicting local Ψ stem values than de-noised ones; (iii) predictions based on band 3 (green @ 560) and band 11 (SWIR @ 1610 nm) are the most accurate ones; and (iv) Ψ stem estimates from different regression methods showed slightly differences in both value and spatial distribution.This alerts us that our deductions are probably not conclusive and, certainly, preliminary.Nevertheless, the approach appears to be promising, and further developments are desirable.

Figure 1 .
Figure 1.Apulia region (green lined), located in Southern Italy.Pomegranate plots N and S ar represented with blue and red lined polygons, respectively.Yellow dots represent the spatia distribution of surveyed pomegranate plants.The S2 T33TYE tile is represented as a blue square in the left frame (the reference system is WGS 84/UTM 33N, EPSG: 32,633).

Figure 1 .
Figure 1.Apulia region (green lined), located in Southern Italy.Pomegranate plots N and S are represented with blue and red lined polygons, respectively.Yellow dots represent the spatial distribution of surveyed pomegranate plants.The S2 T33TYE tile is represented as a blue square in the left frame (the reference system is WGS 84/UTM 33N, EPSG: 32,633).

Figure 2 .
Figure 2. (a) Sentinel 2 GSD (red lines) compared with an orchard pattern (within rows + between rows).Graphs (b-d) show the effects of the joint spectral contribution of the two classes possibly included in a single S2 pixel: (b) spectral signature of "pure" active vegetation; (c) spectral signature of "pure" bare soil (clay); (d) resulting mixed spectral signature.Spectral signatures (b-d) are simply demonstrative (from a ground spectro-radiometer) and used to show the expected spectral effect within a mixed pixel.

Figure 2 .
Figure 2. (a) Sentinel 2 GSD (red lines) compared with an orchard pattern (within rows + between rows).Graphs (b-d) show the effects of the joint spectral contribution of the two classes possibly included in a single S2 pixel: (b) spectral signature of "pure" active vegetation; (c) spectral signature of "pure" bare soil (clay); (d) resulting mixed spectral signature.Spectral signatures (b-d) are simply demonstrative (from a ground spectro-radiometer) and used to show the expected spectral effect within a mixed pixel.

Figure 3 .
Figure 3. Spatial disposition of AOI and AOI2.Pomegranate plots N and S are represented with blue-and red-lined polygons, respectively.The yellow-lined polygon refers to the non-vegetated field.(The reference system is WGS 84/UTM 33N, EPSG: 32,633).

Figure 3 .
Figure 3. Spatial disposition of AOI and AOI2.Pomegranate plots N and S are represented with blueand red-lined polygons, respectively.The yellow-lined polygon refers to the non-vegetated field.(The reference system is WGS 84/UTM 33N, EPSG: 32,633).

Table 1 .
Irrigation and rainfall during the trial period.

Table 1 .
Irrigation and rainfall during the trial period.

Table 8 .
Mean absolute errors (MAE) and normalized MAE values, as computed for the best performing models for the τ i estimate (n/a = result not compliant with model performance, i.e., R 2 of the model for τ i estimate < 0.7).