Evaluation of Normalized Difference Water Index as a Tool for Monitoring Pasture Seasonal and Inter-Annual Variability in a Mediterranean Agro-Silvo-Pastoral System

: Extensive animal production in Iberian Peninsula is based on pastures, integrated within the important agro-silvo-pastoral system, named “montado” in Portugal and “dehesa” in Spain. Temperature and precipitation are the main driving climatic factors affecting agricultural productivity and, in dryland pastures, the hydrological cycle of soil, identiﬁed by soil moisture content (SMC), is the main engine of the vegetation development. The objective of this work was to evaluate the normalized difference water index (NDWI) based on Sentinel-2 imagery as a tool for monitoring pasture seasonal dynamics and inter-annual variability in a Mediterranean agro-silvo-pastoral system. Forty-one valid NDWI records were used between January and June 2016 and between January 2017 and June 2018. The 2.3 ha experimental ﬁeld is located within the “Mitra” farm, in the South of Portugal. Soil moisture content, pasture moisture content (PMC), pasture surface temperature (T ir ), pasture biomass productivity and pasture quality degradation index (PQDI) were evaluated in 12 satellite pixels (10 m × 10 m). The results show signiﬁcant correlations ( p < 0.01) between NDWI and: (i) SMC (R 2 = 0.7548); (ii) PMC (R 2 = 0.8938); (iii) T ir (R 2 = 0.5428); (iv) biomass (R 2 = 0.7556); and (v) PQDI (R 2 = 0.7333). These ﬁndings suggest that satellite-derived NDWI can be used in site-speciﬁc management of “montado” ecosystem to support farmers’ decision making.


Introduction
Extensive animal production in Portugal, especially of cows and sheep, is based on pastures integrated within the context of "montado", an agro-silvo-pastoral system, mainly under a sparse canopy of holm oak (Quercus ilex ssp. rotundifolia Lam.) [1]. The tree cover is mostly a result of natural regeneration and the evergreen canopy can shade 25% to 100% of the land. The trees can be valuable as fodder crop, providing acorns and leafy branches in autumn and winter, when the herbage production is low. The understory consists of evergreen shrubs and a ground layer covered by annual herbaceous species that dry and die off at the end of spring [2]. Permanent pastures constitute the herbaceous layer, a major source of feed for livestock, and are of major importance from an ecological point of view and as a carbon sink [3]. These are pastures with low productivity growing on poor and marginal soils, involving little use of fertilizers, and thus only capable of sustaining low animal allotments, resulting in very small economic returns for the farmers [4]. However, there is a growing pressure to feed the On the other hand, measuring productivity and pasture quality requires cutting and laboratory processing of representative samples, which cannot be carried out by farmers for large and heterogeneous plots. Various published works have focused on the "Grassmaster II" probe [9], which provides very acceptable estimates of pasture productivity, although has the main limitation of requiring an on site survey and measurement of discreet points. The active optical sensors (OAS), such as "OptRx", can be mounted on mobile platforms and associated with a GNSS (global navigation satellite systems), antenna [10], and used to measure NDVI (normalized difference vegetation index; Equation (1), which is the normalized difference between the near infrared (NIR) and visible red reflectance, corresponding respectively to the strong chlorophyll absorption region and high reflectance plateau of vegetation canopies [6]. This index is responsive to changes in both the chlorophyll content and the intracellular spaces in the spongy mesophyll of plant leaves. The anatomy and physiology of the leaf absorbs radiation from wavelengths around the red area of the spectrum and reflects the near infrared [3]. Higher NDVI values reflect greater vigor and photosynthetic capacity (or greenness) of vegetation canopy, whereas lower NDVI values are reflective of vegetative stress resulting in chlorophyll reductions and changes in the leaves' internal structure due to wilting [3][4][5][6][7][8]. The two aforementioned channels used in NDVI can sense through different depths of vegetation canopies: near infrared channel can "see" through roughly eight leaf layers, while the red channel sees only one leaf layer or less because of the strong chlorophyll absorption [11].
where NIR-Near infrared; RED-red. These OAS provide a very interesting response in terms of crude protein (CP), neutral detergent fiber (NDF) and pasture quality degradation index (PQDI) estimates [10], although they require the presence of a technician and running a mobile platform in repeated swats to adequately cover the field in order to obtain representative maps of the spatial variability of the parameters being measured.
As an alternative to this set of proximal sensors, satellite imagery (such as Sentinel-2), provided with very interesting temporal and spatial resolutions, is currently being studied by many research teams for very promising applications related to agriculture. Vegetation indexes, resulting from empirical combinations of reflectance in various wavelengths, have proved very useful for detecting changes in biomass status as measured by changes in chlorophyll activity (biomass vegetation indexes, BVI) and vegetation water content (moisture vegetation indexes, MVI) [2]. Indexes obtained by remote sensing have, in comparison to proximal sensors, important advantages in terms of the scale of the response, speed of processing, and low cost [12]. The major disadvantage of optical satellite imagery is the dependence on a clear, cloud-free view of the object of interest, which is especially challenging in temperate and rainy regions [13], leading to studies on complementing remote sensing by satellite, with proximal sensing, specially in the Mediterranean-type "montado" ecosystem, where the tree canopy obstructs the access to the pasture under the trees.
As already mentioned, NDVI is an index primarily related to canopy chlorophyll content [14] and one of the most widely used indexes for the remote sensing of vegetation [3,11]. Nonetheless, the main engine for the development of the dryland pastures is the hydrological cycle of the soil, and thus the temporal tracking of vegetation mass with spectral measures has been widely investigated, especially in water-limited regions [2], which is the case of the Mediterranean region. The accurate assessment of the seasonal dynamic of drought in these regions by the use of remote-sensing and field observations is essential to determination of the major constraints of such ecosystems [2].
The normalized difference water index (NDWI, Equation (2); [11]) is an index which can be obtained from band 8A (B8A; 865 nm) and band 11 (B11; 1610 nm) of Sentinel-2 spectral data and which can be of great interest as a support to decision making in terms of pasture and grazing management. It was initially proposed for describing the status of the vegetation water content over large areas from space: the NDWI values over areas with dead grass are negative, while those over green vegetation are positive [11]. The NDWI is a measure of liquid water molecules in vegetation canopies that interact with the incoming solar radiation [11], specially conceived for the estimate of soil moisture and canopy water content [5,14]. The NDWI is often a function of local climate and soil properties controlling water availability [5], and sensitive to changes in liquid water because it incorporates a short-wave infrared (SWIR) band (B8A; Equation (2)). This green vegetation spectra region is dominated by water absorption effects that capture important information on seasonally variable water status [5,15]. The NDWI responds to changes in both the water content (absorption of B8A radiation) and spongy mesophyll (reflectance of B11 radiation) in vegetation canopies, respectively [11].
The objective of this work was to evaluate NDWI based on Sentinel-2 imagery as a tool for monitoring pasture seasonal dynamics and inter-annual variability in a Mediterranean agro-silvo-pastoral system.  canopies that interact with the incoming solar radiation [11], specially conceived for the estimate of soil moisture and canopy water content [5,14]. The NDWI is often a function of local climate and soil properties controlling water availability [5], and sensitive to changes in liquid water because it incorporates a short-wave infrared (SWIR) band (B8A; Equation 2). This green vegetation spectra region is dominated by water absorption effects that capture important information on seasonally variable water status [5,15]. The NDWI responds to changes in both the water content (absorption of B8A radiation) and spongy mesophyll (reflectance of B11 radiation) in vegetation canopies, respectively [11].

Materials and Methods
The objective of this work was to evaluate NDWI based on Sentinel-2 imagery as a tool for monitoring pasture seasonal dynamics and inter-annual variability in a Mediterranean agro-silvo-pastoral system. Chronological diagram of data acquisition of daily mean temperature, daily accumulated rainfall, soil moisture content (SMC), pasture surface temperature (Tir), pasture productivity and quality, normalized difference vegetation index (NDVI) and normalized difference water index (NDWI) Sentinel-2 satellite-based.
The relief of the field is slightly wavy, with altimetry between 215 and 219 m, which is representative of the typical pattern in the region. Chronological diagram of data acquisition of daily mean temperature, daily accumulated rainfall, soil moisture content (SMC), pasture surface temperature (T ir ), pasture productivity and quality, normalized difference vegetation index (NDVI) and normalized difference water index (NDWI) Sentinel-2 satellite-based.
The relief of the field is slightly wavy, with altimetry between 215 and 219 m, which is representative of the typical pattern in the region.

Climate
The Mediterranean climate is characterized by dry summers, with less than 10 mm average monthly rainfall, and mild or moderately cold winters, with an unpredictable rainfall regime. The monthly average temperature range is 8 °C to 26 °C and the minimum temperatures are close to 0°C in the winter time and maximum temperatures are close to 40 °C in the summer. The annual rainfall range is 400 to 600 mm and it mainly takes place between October and March.
The experiment was carried out between September 2015 and June 2018. Figure 3 shows the corresponding daily thermo-pluviometric diagram of the "Mitra" meteorological station. This period covers the complete seasonal cycle, and therefore, a wide variety of climate and soil water content situations taking into account the marked irregularity of inter-annual distribution of the rainfall in this region.
The accumulated rainfall between September and June was 544.1 mm in 2015-2016, 412.3 mm in 2016-2017 and 610 mm in 2017-2018. Particularly in the months between February and June, the accumulated rainfall was 327 mm in 2016, 167 mm in 2017 and 461 mm in 2018, which allow us to classify these three springs as average, dry and wet, respectively. From the end of June onwards, a significant drop in the rainfall and the higher temperatures naturally lead to the end of the pasture's productive cycle, which only resumes the following autumn.

Climate
The Mediterranean climate is characterized by dry summers, with less than 10 mm average monthly rainfall, and mild or moderately cold winters, with an unpredictable rainfall regime. The monthly average temperature range is 8 • C to 26 • C and the minimum temperatures are close to 0 • C in the winter time and maximum temperatures are close to 40 • C in the summer. The annual rainfall range is 400 to 600 mm and it mainly takes place between October and March.
The experiment was carried out between September 2015 and June 2018. Figure 3 shows the corresponding daily thermo-pluviometric diagram of the "Mitra" meteorological station. This period covers the complete seasonal cycle, and therefore, a wide variety of climate and soil water content situations taking into account the marked irregularity of inter-annual distribution of the rainfall in this region.
The accumulated rainfall between September and June was 544.1 mm in 2015-2016, 412.3 mm in 2016-2017 and 610 mm in 2017-2018. Particularly in the months between February and June, the accumulated rainfall was 327 mm in 2016, 167 mm in 2017 and 461 mm in 2018, which allow us to classify these three springs as average, dry and wet, respectively. From the end of June onwards, a significant drop in the rainfall and the higher temperatures naturally lead to the end of the pasture's productive cycle, which only resumes the following autumn. In December 2016, 12 PVC tubes were installed in the ground to a depth of 0.20 m. Soil moisture content measured at 0.20 m depth is regarded as a proxy of the water availability for the plant photosynthetic activity due to their shallow root depth [5].   Pasture surface temperature measurements were carried out in each sampling point monthly, always at same time (9 a.m.), between January 2017 and June 2018. An infrared camera (ThermaCAM TM FLIR Systems, Wilsonville, OR, USA) was used. Thermography images were analyzed using ResearchIR ® 3.0 (FLIR Systems, Wilsonville, OR, USA) software. Processed data from the infrared images were exported to a spreadsheet for calculation of the mean and standard deviation of pasture surface temperature. At each sampling point, pasture was cut at 1 to 2 cm above ground level and collected. It was then (i) weighted immediately to obtain fresh mass; (ii) dehydrated (72 h at 65 • C); and (iii) weighted again, in order to obtain pasture biomass production (kg·ha −1 ) and pasture moisture content (PMC, %). The dehydrated samples were analysed (both in % of dry matter) in order to determine the content of CP and NDF [18]. Based on these values, the PQDI, ratio between NDF and CP, was calculated.

Vegetation Multispectral Measurements by Remote Sensing
Sentinel-2 optical images (freely available from the European Space Agency, ESA), were used. For this work, Sentinel-2 band 4 (B4; 10 m spatial resolution; 665 nm), band 8 (B8; 10 m spatial resolution; 842 nm), band 8A (B8A; 20 m spatial resolution; 865 nm) and band 11 (B11; 20 m spatial resolution; 1610 nm), atmospherically corrected imagery, were extracted from Copernicus data hub and used to calculate NDVI (Equation (1)) and NDWI (Equation (2)). This data was obtained for the 12 geo-referenced pixels where the measurements of SMC and T ir were carried out and pasture samples were collected. In order to make the reconstruction of historical NDVI and NDWI trends, time series NDVI and NDWI records were retrieved between January and June 2016 and between January 2017 and June 2018. A preliminary processing was carried out on these records to remove outliers due to the presence of clouds. Only the images without presence of clouds were used in the analysis.

Statistical Analysis of the Data
The mean and standard deviation were calculated for all parameters (biomass, PMC, CP, NDF, PQDI, SMC, T ir , NDVI and NDWI) data.
Regression analysis (p < 0.05) was performed between (i) NDWI and NDVI (for all valid records between January and June 2016 and between January 2017 and June 2018; Figure 1); (ii) NDWI, SMC and T ir (for all monthly sampling date, between January 2017 and June 2018; Figure 1); (iii) NDWI and pasture parameters (biomass, PMC and PQDI; for all monthly sampling date between February or March and June 2016-2018; Figure 1), using the proportion of explained variance as estimated by the coefficient of determination (R 2 ), using 'MSTAT-C' software.

Results and Discussion
3.1. Evolution of Soil Moisture Content (SMC) and Pasture Surface Temperature (T ir ) Table 1 shows the evolution of the mean SMC and T ir of the "Mitra" experimental field between January 2017 and June 2018. Between January and June 2017, as the seasons succeeded one another (winter: January-March; and spring: April-June), there was a tendency for a decrease in SMC and an increase in T ir , which reflects the natural effect of climate conditions, specially the air temperature/ rainfall pair (see Figure 3). Soil moisture is a critical component in land surface-atmospheric processes, and prolonged soil moisture deficits often lead to drought induced vegetation stress [8]. Pasture surface temperature retrieved from thermal infrared (T ir ) provides information about temperature extremes encountered during crop development in this region: extremely low in winter and extremely high in summer, and the consequences on stress and stomatal opening/closure process [19]. In the summer period, especially in the months of July and August, T ir reaches maximum values, while soil moisture has an opposite behaviour. After the arrival of autumn and during winter there is a gradual decrease in T ir and an increase in SMC. Nonetheless, in spring of 2018, the concentration of rainfall during the months of March, April and May (see Figure 3) delayed the normal decrease in SMC and even induced an irregular behaviour of T ir . The combination of these conditions (namely the persistence of relatively high SMC values and the less drastic raise in T ir ) impacts the evolution of the growth cycle of dryland pastures, both in terms of productivity, and quantity [4]. According to Wessels et al. [20], SMC is the main factor determining the composition, distribution and growth status of vegetation in arid and semi-arid regions. The temporal variability of SMC is a relevant factor, namely in the management of irrigation in intensive production systems (irrigated pastures and forages) and in the choice of the best opportunity for fertilizer topdressing, since it requires soil moisture in order to have a real impact on the pasture. Figure 4 shows a negative significant correlation between SMC and T ir (R 2 = 0.6656; p < 0.01) over sampling times (between January 2017 and June 2018).
The combined effect of air temperature and rainfall distribution has an inverse effect on SMC and T ir ; the two factors that control the vegetative development: a gradual decrease in rainfall and a gradual increase in air temperature tend to decrease SMC and increase T ir , and vice versa. For example, since photosynthesis of plants from temperate regions is optimized at 20-25 • C [21], it is possible to identify three physiological periods in the Mediterranean region: (i) winter period, in which transpiration is minimum and the stomata open only to allow gas exchange; (ii) the spring and autumn period, in which the soil moisture is sufficient for plants to control the leaf temperature through The combined effect of air temperature and rainfall distribution has an inverse effect on SMC and Tir; the two factors that control the vegetative development: a gradual decrease in rainfall and a gradual increase in air temperature tend to decrease SMC and increase Tir, and vice versa. For example, since photosynthesis of plants from temperate regions is optimized at 20-25 °C [21], it is possible to identify three physiological periods in the Mediterranean region: (i) winter period, in which transpiration is minimum and the stomata open only to allow gas exchange; (ii) the spring and autumn period, in which the soil moisture is sufficient for plants to control the leaf temperature through transpiration; and (iii) the summer period, in which soil moisture is insufficient for control of the leaf temperature by the plants, leading to senescence and crop maturity.

Evolution of Pasture Parameter Patterns Throughout the Vegetative Cycle
According to Zhao et al. [22], grassland biomass and quality are the main direct indicators of grassland's ecological status. Table 2 shows the evolution of the mean pasture parameters (biomass, PMC, CP, NDF and PQDI) of the "Mitra" experimental field in three consecutive springs (2016-2018).

Evolution of Pasture Parameter Patterns Throughout the Vegetative Cycle
According to Zhao et al. [22], grassland biomass and quality are the main direct indicators of grassland's ecological status. Table 2 shows the evolution of the mean pasture parameters (biomass, PMC, CP, NDF and PQDI) of the "Mitra" experimental field in three consecutive springs (2016-2018). Taking into account the marked climatic seasonality in the Mediterranean area, a three-year period was selected to cover diverse seasonal cycles, and therefore, study a wide variety of climate and, consequently, of SMC situations. Figure 5a shows the evolution of pasture productivity (biomass) in spring of 2016, 2017 and 2018. As expected in dryland pastures, this evolution reflects the effect of rainfall distribution (Figure 6a,b). It is evident that the most productive spring (2016) was one in which rainfall was evenly distributed in autumn, winter and spring (2015-2016).  The agricultural year September 2017-August 2018 can be considered "a relatively unbalanced year" since despite an accumulated rainfall of 612 mm, the autumn was relatively dry, with only 52 mm of accumulated rainfall between September and November ( Figure 6 (b)). This lack of rainfall at the start of pasture growth compromised its autumn development and, consequently, the spring development. Thus, the high concentration of rainfall in the spring (with 361 mm of accumulated rainfall between March and May, Figure 6 (b)), proved to be too late, and did not lead to an increase in productivity due to the advanced stage of plant development. This accumulated rainfall and the consequent increase in the soil moisture, nonetheless prolonged the growth stage of the pasture, causing the definite decrease in productivity to happen only after DOY 150 (in June; Figure 5 (a)).
These results indicate that climatic factors such as precipitation and temperature affect the spatial distribution and inter-annual fluctuation of biomass [23]. According to Gao et al. [24], the year-to-year biomass production pattern is generally consistent with the inter-annual variation in The productivity in February-March (at the end of winter) is the baseline resulting from the plant growth provided by the meteorological conditions of autumn, since the low temperatures of winter (between December and February in Mediterranean region), with minimum temperature below 10 • C, lead to a period of plant dormancy which affects grassland biomass [16].
The agricultural year of September 2015-August 2016, with an accumulated rainfall of 547 mm, can be considered "a relatively balanced year". The peak rainfall in October allowed the germination and growth of the pasture. The regular rainfall during winter and spring (with a precipitation peak in May), ensured sufficient SMC which, in association with favourable air temperature, lengthened the growth cycle of the pasture. After the production peak of spring, which in 2016 occurred at the beginning of May (Table 2), an important decrease of PMC was observed, which reflects the corresponding loss of pasture quality. The biomass production curve in 2016 is typical of years with a balanced rainfall distribution under the Mediterranean conditions of Alentejo [4].
The agricultural year September 2016-August 2017, with an accumulated rainfall of only 421 mm can be considered "a relatively dry year". Despite the rainy autumn (with 171 mm of accumulated rainfall between September and November; Figure 6b), which provided the necessary conditions for initial development of pasture, the winter and spring were relatively dry (with only 236 mm of rainfall in the six-month period of December to May), which compromised the peak vegetative growth and pasture productivity in the spring. The decrease in productivity started in the month of April, from DOY 100 (Figure 5a).
The agricultural year September 2017-August 2018 can be considered "a relatively unbalanced year" since despite an accumulated rainfall of 612 mm, the autumn was relatively dry, with only 52 mm of accumulated rainfall between September and November (Figure 6b). This lack of rainfall at the start of pasture growth compromised its autumn development and, consequently, the spring development. Thus, the high concentration of rainfall in the spring (with 361 mm of accumulated rainfall between March and May, Figure 6b), proved to be too late, and did not lead to an increase in productivity due to the advanced stage of plant development. This accumulated rainfall and the consequent increase in the soil moisture, nonetheless prolonged the growth stage of the pasture, causing the definite decrease in productivity to happen only after DOY 150 (in June; Figure 5a).
These results indicate that climatic factors such as precipitation and temperature affect the spatial distribution and inter-annual fluctuation of biomass [23]. According to Gao et al. [24], the year-to-year biomass production pattern is generally consistent with the inter-annual variation in the growing season precipitation, but an opposite coupled pattern is observed with the growing season temperature.
This inter-annual variability introduces an uncertainty factor for the farm manager, and can affect the dynamic management of grazing and fertilization.
The effect of the agricultural year (in terms of rainfall distribution) is also visible in the quality of the pasture (CP and NDF, Table 2). According to Sala and Parton [25], precipitation is an essential factor in controlling biodiversity and ecosystem stability of terrestrial biomes, especially in the case of arid and semiarid ecosystems. Figure 5b shows the evolution of PQDI in the spring of 2016, 2017 and 2018. Since PQDI is the result of the ratio between NDF and CP, its evolution, although variable from year to year, shows the general tendency for a decrease in the quality of pasture, as the growth cycle of Mediterranean pastures evolves [4]. The degradation of the pasture quality in June 2017 (PQDI = 14.9) at the end of the spring of a dry year is particularly significant, as it emphasizes the importance of rainfall for prolonging the nutritive value of pasture. According to Gu et al. [8] soil moisture is a critical component in land surface atmospheric processes and prolonged soil moisture deficits often lead to drought induced vegetation stress. This stage of greater degradation of pasture quality at the end of spring, also observed by Chakround et al. [2], begins, on average, from the beginning of the month of May (DOY 120; Figure 5b), but can begin at different times depending on the weather. It is a key moment for the farm manager in terms of anticipation of the need to provide feed supplements [10] and a determining factor in terms of grazing dynamics. Table 3 shows the evolution of the mean NDWI and NDVI at the first semester of 2016, all the year of 2017, and in the first semester of 2018. A total of 41 valid records are presented (a preliminary screening was carried out to remove outliers due to existence of clouds). The highest temporal frequency (5 days) of records observed in 2018 (relative to the temporal resolution of 10 days in 2016 and 2017) provided a larger number of valid records (without cloud interference) in this year, increasing the interest and the potential of these indexes in monitoring crop development. Despite inter-annual differences, the general trend consists of a rise in NDWI and NDVI values at the beginning of spring, reaching a maximum value during April (DOY 90-120), followed by a clear decrease. This pattern, typical of dryland ecosystems, reflects the conjugation of the effects of air temperature and precipitation and their influence on SMC, PMC and vegetative vigour. Vegetation index, NDVI, is responsive to changes in both the chlorophyll content and the intracellular spaces in spongy mesophyll of plant leaves: higher NDVI values reflect greater vigour and photosynthetic capacity of vegetation, whereas lower NDVI values are reflective of vegetative stress resulting in chlorophyll reductions and changes in the leaves' internal structure due to wilting [8]. On the other hand, NDWI responds to changes in both the water content and spongy mesophyll in vegetation that interacts with the incoming solar radiation [11]. Figure 7 shows the strong and significant correlation between NDWI and NDVI (R 2 = 0.8621; p < 0.01) in the set of 41 sampling moments considered in Table 3 (between January and June 2016 and between January 2017 and June 2018), as was also observed in other works [8,11]. Serrano et al. [10] obtained strong and significant correlations between NDVI, which is mainly related to canopy chlorophyll content [14] and pasture quality (NDF, CP and PQDI), showing that this index can be of great interest in creating a system of alarms for the farm manager indicating the need to provide supplementary animal feed. In view of the great influence of SMC and of T ir on the physiological processes that are responsible for biomass production, a more specific index was evaluated in this work. Water index (NDWI) is sensitive to changes in liquid water content of vegetation and, thus, must be viewed not as a replacement of NDVI, but as a way to complement it, as was also proposed by [11]. Continued evaluation of satellite-derived NDWI for vegetation drought monitoring using ground observations (e.g., SMC) is required to better understand how this index responds to soil moisture fluctuations, which is ultimately tied to drought stress in plants [8]. According to Sanchez-Ruiz et al. [5], SMC measured at 10-20 cm depth is regarded as a proxy of the water availability for pasture photosynthetic activity, due to its shallow root depth. Figure 8 shows the temporal tracking of SMC, Tir and spectral NDWI measurements between January 2017 and June 2018. A great similarity can be observed between the NDWI and SMC profiles, while Tir has an inverse tendency. If we take as reference the minimum Tir of 10 °C (shown as a red dotted line in Figure 8), below which plants enter dormancy [4], we can determine the most favourable periods (Tir > 10 °C and SMC > 10%) for the plants to access the nutrients needed for their growth. Figure 8 shows that these occasions take place in March (spring of 2017) or April (spring 2018), which are both associated with peaks of NDWI. Thus this index seems to serve as an interesting indicator for supporting decision making by the farm manager in terms of the timing of application of topdressing fertilizer, which is a very common practice in the region for improving the pastures. Continued evaluation of satellite-derived NDWI for vegetation drought monitoring using ground observations (e.g., SMC) is required to better understand how this index responds to soil moisture fluctuations, which is ultimately tied to drought stress in plants [8]. According to Sanchez-Ruiz et al. [5], SMC measured at 10-20 cm depth is regarded as a proxy of the water availability for pasture photosynthetic activity, due to its shallow root depth. Figure 8 shows the temporal tracking of SMC, T ir and spectral NDWI measurements between January 2017 and June 2018. A great similarity can be observed between the NDWI and SMC profiles, while T ir has an inverse tendency. If we take as reference the minimum T ir of 10 • C (shown as a red dotted line in Figure 8), below which plants enter dormancy [4], we can determine the most favourable periods (T ir > 10 • C and SMC > 10%) for the plants to access the nutrients needed for their growth. Figure 8 shows that these occasions take place in March (spring of 2017) or April (spring 2018), which are both associated with peaks of NDWI. Thus this index seems to serve as an interesting indicator for supporting decision making by the farm manager in terms of the timing of application of topdressing fertilizer, which is a very common practice in the region for improving the pastures.  Figure 9 shows the relationship between SMC and NDWI and between Tir and NDWI along the sampling times (between January 2017 and June 2018). The significant correlation between NDWI and SMC (R 2 = 0.7548; p < 0.01) shows that satellite-derived NDWI is valid for monitoring SMC, which confirms previous works [2,8,11]. These results corroborate the findings of Wang et al. [15] which suggests that shortwave infrared (SWIR) band of the "water index" (NDWI) captures important information on seasonally variable water status throughout the vegetative cycle of pastures. This relation is particularly interesting in water-limited regions, which is the case of the Mediterranean region, where it is essential to determine the main constraints of the ecosystem [2].

Correlation between NDWI and NDVI, SMC, T ir , Pasture Moisture Content (PMC), Biomass and Pasture Quality Degradation Index (PQDI)
Having verified the correlation of NDWI with SMC, and evaluated NDWI as a "water index" [8,11], associated to canopy water content and soil moisture [14], sensitive to the total amount of liquid water in the leaves [11], the next step is to evaluate the correlation between this index and PMC ( Figure 10). The great similarity of patterns of the temporal variability of PMC and spectral NDWI (average monthly values measured between February or March and June of three successive years, 2016-2018) indicates that NDWI describes the status of the vegetation moisture content, which is often a function of local climate and soil properties controlling water availability [5]. These results are in agreement with the considerations of Sanchez-Ruiz et al. [5] that the spectral signature of vegetation in NIR and SWIR bands can be associated with the plant water status. The NDWI is sensitive to changes in leaf water content because green vegetation spectra in the SWIR region is dominated by water absorption effects [5], whereby, the quantity of leaf material can be known in order to estimate vegetation dryness [26].  Figure 9 shows the relationship between SMC and NDWI and between T ir and NDWI along the sampling times (between January 2017 and June 2018). The significant correlation between NDWI and SMC (R 2 = 0.7548; p < 0.01) shows that satellite-derived NDWI is valid for monitoring SMC, which confirms previous works [2,8,11]. These results corroborate the findings of Wang et al. [15] which suggests that shortwave infrared (SWIR) band of the "water index" (NDWI) captures important information on seasonally variable water status throughout the vegetative cycle of pastures. This relation is particularly interesting in water-limited regions, which is the case of the Mediterranean region, where it is essential to determine the main constraints of the ecosystem [2].
Having verified the correlation of NDWI with SMC, and evaluated NDWI as a "water index" [8,11], associated to canopy water content and soil moisture [14], sensitive to the total amount of liquid water in the leaves [11], the next step is to evaluate the correlation between this index and PMC ( Figure 10). The great similarity of patterns of the temporal variability of PMC and spectral NDWI (average monthly values measured between February or March and June of three successive years, 2016-2018) indicates that NDWI describes the status of the vegetation moisture content, which is often a function of local climate and soil properties controlling water availability [5]. These results are in agreement with the considerations of Sanchez-Ruiz et al. [5] that the spectral signature of vegetation in NIR and SWIR bands can be associated with the plant water status. The NDWI is sensitive to changes in leaf water content because green vegetation spectra in the SWIR region is dominated by water absorption effects [5], whereby, the quantity of leaf material can be known in order to estimate vegetation dryness [26].      Figure 11 shows the relationship between NDWI and PMC and between NDWI and PQDI. Pasture moisture content is a parameter indicative of plant vegetative vigour, with a tendency to decease at the end of the growth cycle. In the case of "relatively dry year" (2017), PMC in June was extremely low (14.7 ± 4.8%), corresponding, as mentioned above, to the highest value of PQDI (14.9 ± 2.5). These two parameters (PMC and PQDI), have, thus an inverse behaviour in the relation with NDWI (Figure 11), although both with significant determination coefficients (R 2 = 0.8938; p < 0.01 to PMC; R 2 = 0.7333; p < 0.01 to PQDI). Figure 11 shows the relationship between NDWI and PMC and between NDWI and PQDI. Pasture moisture content is a parameter indicative of plant vegetative vigour, with a tendency to decease at the end of the growth cycle. In the case of "relatively dry year" (2017), PMC in June was extremely low (14.7 ± 4.8%), corresponding, as mentioned above, to the highest value of PQDI (14.9 ± 2.5). These two parameters (PMC and PQDI), have, thus an inverse behaviour in the relation with NDWI (Figure 11), although both with significant determination coefficients (R 2 = 0.8938; p < 0.01 to PMC; R 2 = 0.7333; p < 0.01 to PQDI). PQDI Figure 12 shows the relationship between NDWI and biomass. The resulting curve has a tendency for over-saturation, and consequently a decline in the accuracy/sensitivity of this index at high values of biomass (above 20000 kg·ha −1 ), an aspect documented in various works related to NDVI ( [10,22,[27][28][29][30]. These results show that NDWI based on remote sensing can be an interesting tool to estimate pasture biomass (R 2 = 0.7556; p < 0.01). Pasture biomass is a vital input for the management of dairy systems, allowing the calculation of feed budget, based on which decisions are made for farm practices such as rotational lengths and supplementary feeding, leading to profitability and a sustainable use of pasture resources [31].
In comparison with the use of NDVI in pastures [10], these results suggest that NDWI has a more broad application potential, and is appropriate for monitoring drought stress in vegetation, through plant available water (SMC), light-harvesting capacity (biomass) or crop progress and vegetation health (PMC and PQDI, respectively) [19]. The response to the level of SMC and Tir is particularly important, since these parameters are drivers of dryland pasture growth [28].  Figure 12 shows the relationship between NDWI and biomass. The resulting curve has a tendency for over-saturation, and consequently a decline in the accuracy/sensitivity of this index at high values of biomass (above 20,000 kg·ha −1 ), an aspect documented in various works related to NDVI ( [10,22,[27][28][29][30]. These results show that NDWI based on remote sensing can be an interesting tool to estimate pasture biomass (R 2 = 0.7556; p < 0.01). Pasture biomass is a vital input for the management of dairy systems, allowing the calculation of feed budget, based on which decisions are made for farm practices such as rotational lengths and supplementary feeding, leading to profitability and a sustainable use of pasture resources [31].
In comparison with the use of NDVI in pastures [10], these results suggest that NDWI has a more broad application potential, and is appropriate for monitoring drought stress in vegetation, through plant available water (SMC), light-harvesting capacity (biomass) or crop progress and vegetation health (PMC and PQDI, respectively) [19]. The response to the level of SMC and T ir is particularly important, since these parameters are drivers of dryland pasture growth [28].
In summary, the practical interest of this study is related with the integration of satellite-derived NDWI and NDVI in site-specific management of the "montado" ecosystem This information can support the decision making process when the farmer intends to intensify the production (irrigation or fertilization management, animal feed supplementation or dynamic grazing management), based on soil moisture content (SMC) and pasture productivity (biomass) and quality (CP, NDF, PQDI and PMC) monitoring. In summary, the practical interest of this study is related with the integration of satellite-derived NDWI and NDVI in site-specific management of the "montado" ecosystem This information can support the decision making process when the farmer intends to intensify the production (irrigation or fertilization management, animal feed supplementation or dynamic grazing management), based on soil moisture content (SMC) and pasture productivity (biomass) and quality (CP, NDF, PQDI and PMC) monitoring.

Conclusions
In the Iberian Peninsula, dryland natural pastures are the basis of animal feeding under extensive production regimes. The annual climatic seasonality, associated to inter-annual irregularity of precipitation in this semi-arid region, affects soil moisture content and surface pasture temperature, driver factors of vegetation development. This work demonstrates that remote sensing can provide valuable large-scale spatial and temporal information about yield-limiting factors and crop response under current and variable climate conditions. The normalized difference water index (NDWI) shows strong and significant correlations with soil moisture content (R 2 = 0.7548), pasture moisture content (R 2 = 0.8938), surface pasture temperature (R 2 = 0.5428), biomass productivity (R 2 = 0.7556) and pasture quality degradation index (R 2 = 0.7333). These findings suggest that satellite-derived NDWI can be used in site-specific management of the "montado" ecosystem, as an indicator of critical importance for monitoring and predicting deterioration/restoration of pasture and to support decision making related to agronomic practices for the sustainable use of grassland resources. These decisions include differential soil amendment or fertilization, pasture irrigation, or dynamic grazing. This can be an interesting research topic allowing the assessment of dryness or drought stress on vegetation in Mediterranean-type ecosystems through plant available water, light-harvesting capacity or vegetation health, with applications in other important crops in the region, such as the vineyard or the olive grove. However, given the complexity and number of