Next Article in Journal
Removal of Phosphate from Wastewater with a Recyclable La-Based Particulate Adsorbent in a Small-Scale Reactor
Next Article in Special Issue
CSM-CERES-Wheat Sensitivity to Evapotranspiration Modeling Frameworks under a Range of Wind Speeds
Previous Article in Journal
The Impacts of Precipitation on Fluorescent Dissolved Organic Matter (FDOM) in an Urban River System
Previous Article in Special Issue
Searching for Sustainable-Irrigation Issues of Clementine Orchards in the Syrian Akkar Plain: Effects of Irrigation Method and Canopy Size on Crop Coefficients, Transpiration, and Water Use with SIMDualKc Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Spatio-Temporal Dynamics of Deep Percolation Using Crop Evapotranspiration Derived from Earth Observations through Google Earth Engine

1
Department of Biosystems Engineering, Instituto Superior de Agronomia, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisbon, Portugal
2
LEAF-Linking Landscape, Environment, Agriculture and Food-Research Center, Associated Laboratory TERRA, Instituto Superior de Agronomia, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisbon, Portugal
*
Author to whom correspondence should be addressed.
Water 2022, 14(15), 2324; https://doi.org/10.3390/w14152324
Submission received: 21 June 2022 / Revised: 15 July 2022 / Accepted: 25 July 2022 / Published: 27 July 2022

Abstract

:
Excess irrigation may result in deep percolation and nitrate transport to groundwater. Furthermore, under Mediterranean climate conditions, heavy winter rains often result in high deep percolation, requiring the separate identification of the two sources of deep percolated water. An integrated methodology was developed to estimate the spatio-temporal dynamics of deep percolation, with the actual crop evapotranspiration (ETc act) being derived from satellite images data and processed on the Google Earth Engine (GEE) platform. GEE allowed to extract time series of vegetation indices derived from Sentinel-2 enabling to define the actual crop coefficient (Kc act) curves based on the observed lengths of crop growth stages. The crop growth stage lengths were then used to feed the soil water balance model ISAREG, and the standard Kc values were derived from the literature; thus, allowing the estimation of irrigation water requirements and deep drainage for independent Homogeneous Units of Analysis (HUA) at the Irrigation Scheme. The HUA are defined according to crop, soil type, and irrigation system. The ISAREG model was previously validated for diverse crops at plot level showing a good accuracy using soil water measurements and farmers’ irrigation calendars. Results show that during the crop season, irrigation caused 11 ± 3% of the total deep percolation. When the hotspots associated with the irrigation events corresponded to soils with low suitability for irrigation, the cultivated crop had no influence. However, maize and spring vegetables stood out when the hotspots corresponded to soils with high suitability for irrigation. On average, during the off-season period, deep percolation averaged 54 ± 6% of the annual precipitation. The spatial aggregation into the Irrigation Scheme scale provided a method for earth-observation-based accounting of the irrigation water requirements, with interest for the water user’s association manager, and at the same time for the detection of water losses by deep percolation and of hotspots within the irrigation scheme.

1. Introduction

Given the climatic conditions of the Mediterranean regions, groundwater reserves play an important role in interannual regularization, minimizing the impact of climate variability, particularly droughts. However, the conservation of groundwater, in terms of quantity and quality, is a major environmental challenge [1]. The conservation of quantity requires the maintenance of its recharge, while the conservation of quality requires the existence of low concentrations of nitrates and other elements in the flow of water that reaches the aquifer. The challenge increases if we consider the additional effect of climate change, which could lead to a decrease in aquifer recharge, increased crop irrigation needs, and rising sea levels [2,3]. Due to the characteristics of the Mediterranean climate, where the hot season is also the driest, irrigation is mandatory for achieving high productivity and competitiveness in the agricultural sector. However, intensive land use in terms of irrigation and fertilizer consumption in areas overlying aquifers can lead to risks of contamination, as water, when in excess, or poorly distributed, can transport nitrates into the deepest layers [4]. This situation is particularly important in the Nitrate Vulnerable Areas (ZVN), where nitrate concentration in groundwater has reached/exceeded the limit imposed by the European community (50 mg NO3 L−1). It is the deep percolation of water, resulting from precipitation and irrigation, that serves as a transport of nitrates to groundwater. Thus, the mismanagement of water and fertilizers plays an important role in modifying the natural recharge of aquifers and in the groundwater contamination [5].
Some studies present estimates of aquifer recharge as a percentage of the annual rainfall [6]. However, the deep percolation and recharge associated with the practice of irrigation has been subject of reduced attention. In irrigated areas, in the absence of precipitation, recharging during the summer months can be significant depending on the water balance in the crop rootzone, creating hotspots for potential aquifer contamination with nitrates from fertilizers, organic matter mineralization, etc.
The water balance in the crop rootzone, with the objective of calculating deep percolation, requires the determination of the crop water requirements which in turn relies on the correct estimation of their biophysical parameters, such as crop coefficients and the durations of phenological phases. Although there is information in the literature, as in the recent review by Pereira et al. [7,8], these authors state that these parameters may vary with several factors so regional characterization is advised.
Satellite remote sensing (RS) data have been widely used for phenological monitoring of vegetation in large geographic domains, using various observation methods [9,10,11,12]. Emerging techniques are based on time series analysis of remote detection multispectral imaging which can be used to estimate the crops’ phenological phases (e.g., beginning of the cycle, peak vegetative development, and the duration of phenological phases) [13,14,15]. As recently discussed by Mahlayeye et al. [16], there is a need to further study the cropping patterns at larger scales such as irrigation schemes and watersheds.
Remote sensing temporal data with high spatial and temporal resolution can be acquired, for example, from the Sentinel-2 satellite, which is an Earth Observation (EO) mission of the European Union, the Copernicus program, that systematically acquires high-resolution optical images. Together, Sentinel-2 satellites have high temporal resolution (5 days) and spatial resolutions of 10, 20, and 60 m. These images are made available free of charge by the European Space Agency (ESA) and can be downloaded on various platforms, such as the Google Earth Engine (GEE), which makes these images available in its library as a collection. It has two Sentinel-2 collections, the “COPERNICUS/S2” which corresponds to the reflectance at the top of the atmosphere and the “COPERNICUS/S2_SR” for the surface reflectance [17,18].
Vegetation indices (VI) correspond to an image treatment technique that is applied to extract information such as crop planting density and vegetation vigor (e.g., [19]). Based on these, the values of the actual crop coefficients (Kc act), planting dates, and the duration of the phenological phases can be estimated. They can also be used to identify the soil cropping patterns in a given period [20,21,22]. The Normalized Difference Vegetation Index (NDVI) is a widely used index to characterize crop phenology and it is calculated using the red and near infrared bands present in most sensors [23,24]. Several authors have proposed correlations between NDVI and Kc act based upon field data, establishing a linear relationship between the two variables [10,20,25,26,27], although this relation is affected by other factors including crop architecture. The approach proposed by Calera-Belmonte et al. [20] has been widely used both at the plot scale and at the regional level for obtaining the actual Kc, the length of the crop phenological phases, and the identification of crop types, with the aim of determining the irrigation water requirements [10,13]. The Kc act values are then used to calculate the actual crop evapotranspiration and subsequently perform the soil water balance [26,28,29,30].
The availability of free satellite images with high spatial and temporal resolution makes the use of this technology very interesting for obtaining crop data to improve the estimation of actual ETc. Furthermore, the Google Earth Engine (GEE) [31] which is a “cloud-based” platform for planetary-scale geospatial analysis that combines computing capabilities to the storage and processing of satellite imagery, enables the access to publicly available RS data, such as Landsat 4–8 or Sentinel 1–2, and its data processing and visualization. Thus, GEE allows performance of this type of study without the acquisition of high-performance computers, enabling a wider adoption of EO data on the estimation of ETc act. This service contains several libraries with different types of data and allows the processing of this data on their own servers. Datasets can be acquired and manipulated with a JavaScript code editor platform [31]. GEE allows for rapid data processing, namely the automation of the extraction of temporal series of vegetation indices (e.g., NDVI), without the need for a large data storage capacity, or high-performance computers [31,32,33]. These characteristics of the platform made it possible to overcome important limitations to a greater adoption of high-resolution satellite images, as is the case of Sentinel-2, for crop data retrieval [31]. Vegetation indices such as NDVI can be extracted from satellite images on the platform, and they can be viewed in time series charts or transferred to external processing. GEE was successfully used to estimate crop evapotranspiration [34,35,36], crop productivity [37,38], to cartograph crop areas [39,40,41], for supporting irrigation scheduling [42,43], and water management [44].
Aiming at performing the spatial analysis of the deep percolation and detecting hotspots at an irrigation scheme level, it is necessary to integrate soil water balance modeling with spatially distributed information in a GIS environment. To perform this integration the following steps are required: (a) development of the GIS database, climatic database, soil database, and other databases (administrative limits, irrigation, topographic and land use information); (b) integration of the database and the customization of the GIS; (c) integration of the irrigation simulation model. Soil water balance modeling can be completed in a GIS environment, using a vector approach through the definition of homogeneous analysis units (soil × crop × irrigation system), or using data grids (raster) for climate, soil, and crops [45].
The question to be answered by the present study is: are there hotspot areas in the Caia Irrigation Scheme with high deep percolation due to irrigation practices? To answer this question, the following objectives were established: (1) estimate crop evapotranspiration using remote sensing images to improve the definition of the crop stages; (2) define homogeneous analysis units (soil × crop × irrigation system) in a GIS environment, for which the soil water balance is modeled; (3) estimate the deep percolation associated with the irrigation and with the precipitation along the crop season, and also off season, after crop harvesting; and (4) perform the spatial analysis of deep percolation and identify its potential causes.

2. Materials and Methods

2.1. Caia Irrigation Scheme Case Study

2.1.1. Localization and General Characteristics

The Caia Irrigation Scheme (CIS) is located in the Alentejo region, southern Portugal (Figure 1). The total irrigated area of the CIS, which is around 7000 ha, is included in the Elvas Vulnerable Zone to Nitrates. The CIS serves 887 farmers with water from the Caia Reservoir, located in the Caia River, which also provides water for domestic consumption.
Water is distributed to farms in open channel, with upstream control, and a fixed rotation scheme where each farmer receives a maximum amount of water per hectare.

2.1.2. Climate and Soil Characterization

According to the Köppen classification [46], the climate in the region of study is a Csa—temperate Mediterranean with hot and dry summers. Weather data used in the current study was collected from a nearby weather station (38°54′56″ N, 7°05′56″ W), located under reference site conditions. The long term (2003–2021) average annual precipitation is 519 mm, mostly occurring in winter, while the average annual reference crop evapotranspiration (FAO-PM, [47]), which represents the climatic demand, is 1404 mm. Figure 2 shows the monthly average values of the reference crop evapotranspiration (ETo) and precipitation for the four studied years. It clearly shows a water shortage for the crops during the spring–summer seasons. When comparing the conditions of the studied years with the long-term data, it is shown that 2018 was a humid year with 835 mm of precipitation while 2019 was the driest one with 317 mm. The annual ETo ranged from 1290 to 1375 mm during the studied years.
According to the World Reference Base of Soil Resource [48], the predominant soils in the irrigated perimeter are Fluvissols (44.9%), Luvissols (29.6%), and Calcissols (18.7%).

2.1.3. Main Crops and Irrigation Systems

In 2002, the main crops in CIS were maize (Zea mays L.) occupying 49% of the area, wheat (Triticum aestivum L.) with 17%, sunflower (Helianthus annuus L.) with 7%, tomato (Solanum lycopersicum L.) with 6%, and olive orchards (Olea europaea L.) (super-intensive and hedgerow) with 4% of the total area [49]. However, in the last few years (Figure 3) there has been a significant conversion from areas with annual crops to permanent crops such as olive orchards and nut-tree orchards, including walnut and almond.
These areas have been stabilizing in the last years (Figure 3). According to the Caia water user’s association, in 2020 the irrigated area was 8559.4 ha, because it includes areas that are outside the irrigation scheme but that receive water from the CIS. The main irrigation system is drip with 82% of the irrigated area, followed by sprinkler irrigation with 17.1%, mainly center pivots, while surface irrigation has a very small share of the irrigated area.

2.2. Data Collection and Georeference Database Building

A modeling tool (ISAREG model [50]) was used to support the estimation of the crops water consumption (ETc) and of the deep percolation associated with both irrigation and precipitation along the studied years (off-season was included in the study). To feed the modeling tool, a georeferenced database with diverse information was built. Due to the diversity of data sources, they were edited and converted to the same spatial and temporal resolution.
Data used for the estimation of the deep percolation are summarized in Table 1 and the data sources used to collect the various types of information are also included.

2.2.1. Crop Characteristics Database

The data relative to crop characteristics included the identification of the crop type, spatial distribution of the diverse cropping systems, crop growth stages dates, and crop coefficients for the diverse crop stages. These data were collected from diverse sources as detailed in Table 1.
As previously stated, due to its capabilities the Google Earth Engine (GEE) platform was used in the current study.
  • Satellite imagery processing and crop data collation
Aiming at decreasing the uncertainties in the characterization of the cropping systems and crop growing seasons, Sentinel-2 satellite imagery with high spatial (10 × 10 m) and temporal (5 days) resolution was used. All Sentinel-2 images available on GEE platform for the CIS region between October 2016 and October 2020 were used, with a cloud cover with less than 35%, which corresponds to a compromise between the number of images available to produce the VI time series and its quality. The Sentinel-2 information was then processed in the platform Google Earth Engine using the Java script code editor.
The flowchart presented in Figure 4 shows the global process of selection, validation, and pre-processing of information in order to build the database for crops.
2.
Crop identification and spatialization
Firstly, a layer containing the spatial distribution of the different types of crops was created in QGIS; it included the vector information of the cultivated fields. This layer was based on field information provided by the Caia water user’s association (WUA). In the following and based upon the land use map (COS) (Table 1) the polygons of permanent crops (olive orchards, orchards, and vineyards) were extracted, based on the assumption that these crops vary little from year to year within the studied years (Figure 3). For producing the NDVI [53] for each crop and each field the series of the spectral bands of the Sentinel-2 sensors B4 (red) and B8 (near infrared) and QA60 (cloud mask) were used from the image collection “COPERNICUS/S2_SR”, which corresponds to Level-2 data (orthorectified Bottom-Of-Atmosphere (BOA) reflectance) [18]. The GEE allowed the extraction of the average NDVI per pixel for the period March–May which allowed for the identification of the winter crops (Figure 5a) while the average NDVI for the period May–August allowed identifying the spring–summer crops (Figure 5b). During summer, the NDVI map allowed the precise identification of the irrigated areas. Figure 5c shows an example of the identification of the actual irrigated areas using the Sentinel-2 images. Images were previously processed to discard the pixels that correspond to small roads separating orchard blocks, roads, and edges as the soil conditions in these pixels are not the same as agricultural zones without vegetation; thus, avoiding its influence on the vegetation indices calculation. In addition, the crop NDVI temporal profile analysis allowed identifying the plot fields that did not match the typical crop NDVI pattern, which were removed from its respective crop area [14].
The cropping pattern in CIS includes more than 30 types of crops, and therefore crops were grouped according to their relative importance in terms of cropping area, crop cycle, and irrigation water requirements.
3.
Locally adjusted crop growing stages
Subsequently, the processed satellite images were used to obtain the observed planting dates and the diverse crop growth stage lengths, since the poor characterization of the crop cycle is one of the main sources of uncertainty in the estimation of crop water requirements (e.g., [54]).
The polygons corresponding to the fields with the different crop groups were imported into the GEE platform and the average values of NDVI for each polygon were extracted from the image collections “COPERNICUS/S2” and “COPERNICUS/S2_SR”.
In the following, the NDVI temporal profile for each crop was obtained from the spectral signature of the fields identified with a certain crop; the actual Kc curves (affected by diverse stresses, e.g., water and nutrients) were derived from those profiles using empirical equations [20,25,27]. Atypical NDVI profiles were identified and discarded.
The NDVI temporal profile for each crop allowed estimating the average sowing/planting dates and the lengths of the growing stages. Thus, allowing defining the crop growing stages required for constructing the Kc segmented curve [47] for each crop. These were assumed to be representative for the entire CIS.
4.
Crop coefficients along the growing seasons
The standard potential Kc values for the diverse crop growing stages (initial, mid-season, and end-season) were collected from the tabulated values recently updated by Pereira et al. [7,8] and Rallo et al. [51]. For each year, the Kc values for the initial stage were those provided by Allen et al. [47] but adjusted for soil wetted fraction (irrigation or precipitation). While the crop coefficients for the mid-season and end-season values were adjusted to the prevailing local climate conditions (RHmin and u2) in each studied year using an algorithm based upon the methodologies described in Allen et al. [47].
The crop data used in the computations is presented in Table A1 and Table A2 in Appendix A.

2.2.2. Soil Characteristics Database

Due to the large spatial variability of CIS soils, they were grouped into classes of land suitability for irrigation. These classes result from the integration of three characteristics: soil texture, slope, and permeability. Class I includes soils with excellent suitability for irrigation, with sandy loam or clay texture, good structure and permeability, and a maximum slope of 2%. Soils in class II present good suitability for irrigation, are generally deep soils, and moderately drained. Soils within class III present reasonable suitability for irrigation and are shallow soils with heavy texture, with slopes varying between 4% and 6%. The soil textures were extracted form Portuguese Soil Map (CSP). Slopes were obtained from the digital terrain model (MDT), built using information on contour lines, elevation points, and the hydrographic network, with the support of QGIS. The permeability classes as defined by ARS/USDA [55] were used.

2.2.3. Homogenous Unit of Analysis

Using the information from the previously referred databases, diverse homogenous units of analysis (HUA) were defined at the CIS and the water balance was independently modeled for these units. Each HUA corresponds to an area where the following three features influencing the water balance terms do not spatially change: land suitability for irrigation, crop, and irrigation system. The HUA were obtained in QGis by intersecting the layers containing the three features above.

2.3. Soil Water Balance Calculation

2.3.1. Modeling Strategies and Field Validation

The soil water balance (SWB) was modeled for each HUA using the ISAREG software. The ISAREG model calculates crop evapotranspiration according to the FAO methodology and the SWB follows a reservoir approach. It was used because it has very good accuracy when adequately calibrated and validated against observations of the soil water (e.g., [56]) and/or of crop ET and because it has been applied to diverse climatic and soil conditions. Examples of application include those to field (e.g., [57,58,59]), vegetable crops [60,61], pastures and grasslands [62,63], and orchards [64,65]. Furthermore, the model has been used for the assessment of climate change on irrigation water requirements (e.g., [66,67]) and for irrigation management at irrigation scheme [68] and basin level [69].
The flowchart presented in Figure 6 shows how the water balance is calculated by ISAREG. The SWB was separated in two components, one for the irrigation season, to account for the percolation created by irrigation, and the other for the off-season period, to account for the percolation originated by autumn–winter precipitation. The general water balance equation (Equation (1)) is applied to the crop rootzone:
DP = P + I − ETc act − ΔS
where DP is deep percolation, P is precipitation, I is irrigation (0 for the off-season period), ETc act is actual crop evapotranspiration, and ΔS is the storage variation. Runoff was considered = 0 due to the large scale of application. All terms are in mm.
The information relative to the crop characteristics, derived from EO, along with the farmers actual irrigation calendars were introduced in the modeling tool. In a few cases, the information relative to the irrigation schedules was not available; thus, to overcome the issue the WUA provided for information regarding the total water consumption per field. This information enabled to estimate the irrigation schedules per crop based upon the assumption that the farmers behave similarly in terms of irrigation depths and timing.
Several fields with the most representative crops within the irrigation district (Table 2) were used for validating the soil water balance simulations by comparing the soil water storage for the entire rootzone (ASW) predicted by ISAREG with soil water storage obtained from soil moisture data obtained from continuous monitoring with capacitive probes provided by the Portuguese Association of Precision Agriculture (APAP).
A set of goodness-of-fit indicators was used to assess the performance of the ISAREG model in estimating the soil water storage (ASW, mm). These indicators are: (i) the regression coefficient (b0) describing a linear regression forced through the origin [70] between measured (ASWOBS) and model-simulated ASW (ASWSIM) values; (ii) the percent bias (PBIAS, %), which measures the average tendency of the simulated data to be larger or smaller than their corresponding observations (Gupta et al., 1999); (iii) the coefficient of determination (R2) of the ordinary least squares regression; (iv) the root mean square error (RMSE, mm); (v) the normalized RMSE (NRMSE, %), defined as the ratio between RMSE and the mean of observations O ¯ ; and (vi) the modeling efficiency (EF, non-dimensional), which is an indicator proposed by Nash and Sutcliffe [71] used to assess the magnitude of the mean square error (MSE = RMSE2) relative to the observed data variance [72].

2.3.2. Irrigation Water Requirements

Irrigation water requirements (IWR) were calculated by dividing the net irrigation water requirements calculated with ISAREG, by the average irrigation system efficiency (drip = 0.9; center pivot = 0.8; traveling gun = 0.75; paddy rice basins = 0.7). In the case of permanent crops irrigated by drip systems, the value was multiplied by the reduction factor 0.8 to account for the reduction in the irrigated area at the soil surface [73].

3. Results

3.1. NDVI Temporal Profiles and Crop Coefficient Curves

3.1.1. Permanent Crops

The NDVI profiles for the permanent crops, e.g., vineyards, orchards, and olive groves are shown in Figure 7. The NDVI profiles allowed distinguishing plots with young orchards from those in full production. For the orchards, the end of dormancy was identified to be by March (DOY 90), while leaf senescence occurred by October (DOY 290). The NDVI increase during spring may be associated with the interference of wet soil upon the reflectance [74]. As shown in Figure 7, the presence of clouds along the crops’ seasons affected the precise identification of the various crop growth stages. For the case of the olive groves, the NDVI profiles allowed to distinguish between hedgerow and intensive groves. For the hedgerow, the NDVI is high, generally ranging from 0.3 to 0.6, and is minimally influenced by the presence of weeds in the inter-row spacing, while for the intensive groves, due to its management and higher space between trees the weeds develop particularly during spring and autumn; thus, presenting higher NDVI values. Overall, for all permanent crops, the high NDVI values at the beginning of the year may be associated with the presence of weeds that grow in the inter-row spacing due to winter and spring precipitation events.

3.1.2. Winter Crops

As per the permanent crops, the lack of images for some dates due to the presence of clouds prevented the adequate characterization of the crop growth stages for the winter crops (garlic, cereals, pastures) (Figure 8); thus, data was complemented with that available in the literature (Table A1 in Appendix A).

3.1.3. Annual Spring Crops

For the spring crops, the results of the NDVI temporal profiles allowed for an accurate estimation of the planting/sowing dates as well as for the duration of the crop growth stages; thus, allowing improvement of the definition of the crop coefficient curves, and therefore increasing the accuracy of the water requirement computations. Figure 9 presents examples of the NDVI time series obtained for several spring crops along the 2016/2017 season. The first attempt for the definition of the FAO Kc curves was to adjust it to the mean values of the Kc NDVI after applying the Calera Belmonte et al. [20] relationship to the NDVI profiles; however, this procedure created, in some cases, a flattening of the Kc curve. The use of the average Kc curves showed to be not appropriate for the definition of the crop growth stages because farmers used different varieties and sowing/planting dates. Therefore, a different approach was used in the current study, consisting of adjusting the Kc FAO curves, particularly the duration of the crop phases, to the Kc-NDVI clouds of points. Figure 9 shows the results of the Kc curves obtained for various crops while Table A2, in Appendix A, presents the values for the initial, mid- and end-season actual Kc for all crops in all the considered years.
Abnormal behavior of NDVI occurred mainly during the initial stage resulting in a less accurate determination of the duration of this stage. This behavior occurs mainly for the lower NDVI values, due to instability associated with soil background influence during the initial stage. Similar findings were reported by Cao et al. [74] and in the recent review by Pôças et al. [10]. Other sources of uncertainty in the definition of the Kc-NDVI curves is the use of empirical formulas calibrated and validated for other regions [10,21,27]. Similar findings were reported by Vilar et al. [75] using Landsat-8 images and by Rolim et al. [14] using images from SPOT-5 TAKE-5. Table 3 presents the length of the crop growth stages and the seeding/planting dates for the spring crops, estimated from the NDVI temporal profiles. Results show the adequacy of using the NDVI temporal profiles for improving the characterization of the crop cycles.
For tomato crop, the use of plastic mulch was assumed. In the case of pepper and melon crops, plastic mulch was not considered since most farmers in the CIS do not use it. For the winter vegetables, the Kc values that were ingested in the modeling tool correspond to the most representative crop, i.e., the broccoli because they are the main crop in terms of cropped area. The same reasoning was applied to pepper, which were considered to represent the spring–summer vegetables, while almonds were selected for representing the group of dry fruits. For the pastures and the autumn–winter forages, the winter cereal data were used because they presented similar crop cycles. For the irrigated pastures, an average Kc value of 0.75 was considered for the entire irrigation season.
The cropping patterns in the CIS, for the four studied seasons, after EO data validation against the data provided by the WUA are presented in Figure 10. The maps produced in the GEE platform with the average NDVI per pixel, for each growing period, allowed to identify the areas with winter and spring crops. During summer, the NDVI map allowed the precise identification of the irrigated areas. As previously stated, during the entire study period, olive groves were the predominant crop in CIS spreading across the irrigated perimeter. Winter cereals also occupy an important area, but they are preferentially located on the soil classes of low suitability for irrigation. Maize and fodder crops are the most representative annual spring–summer crops.

3.2. Irrigation Systems

The spatial distribution of the irrigation systems in the CIS for the studied years of 2017–2020 is presented in Figure 11. This information was provided for 2020 by the water users association, while, for 2017–2019, it was estimated from the high-resolution satellite images. As depicted in Figure 11, drip irrigation systems were predominant, followed by the sprinkler irrigation system with center pivots while travelling guns and surface irrigation have little significance.

3.3. Land Suitability for Irrigation

As previously explained, the classes of land suitability for irrigation, which is one of the layers needed for the definition of the HUA, were obtained by overlaying information relative to soil use capacity, land slope, and soil permeability. Figure 12 shows the spatial distribution of the land suitability classes for irrigation within the CIP. Results show that most of the CIS area consists of class III soils, presenting moderate limitations for irrigation, while classes II and I, which present lower limitations for irrigation, are mainly located along the banks of the Caia River. Results also show that the class with the lowest suitability for irrigation (exc) occupies 15% of the total area. Table 4 presents the average soil properties used to characterize the different land suitability classes that were inputted in the modeling tool. These properties were determined considering the representativeness of each soil type within each class.

3.4. Homogeneous Units of Analysis

The produced homogeneous units of analysis (HUA), obtained by intersecting in QGIS the layers containing information of the spatial distribution of crop types, the land suitability for irrigation classes and the irrigation systems, are depicted in Figure 13. The total numbers of HUA were 104, 88, 112, and 109, respectively, for 2016/2017, 2017/2018, 2018/2019, and 2019/2020 (Figure 13). In a certain area, the changes in the HUA from one year to the other are mainly due to changes in the cultivated crops. An example of the detailed description of the HUA can be found in Table A3 in the Appendix A. As previously pointed out, the soil water balance model was run independently for each HUA and the results were spatialized using the GIS platform.

3.5. Soil Water Balance Modeling Results

3.5.1. Model Validation for Soil Water Storage Predictions

The soil water balance simulations were validated using soil moisture time series obtained, for the study years, in three farms planted with maize, processing tomato, olive groves. Based upon the soil moisture measured within the rootzone with capacitance probes presenting sensors each 10 cm, the soil water storage (ASW) was calculated and compared with the simulated values, as shown in Figure 14. It can be observed that the model was able to adequately predict the behavior of the ASW along the seasons for the different crops and soils. In the example, maize was cropped in a loamy sand, tomato in a silty loam, and olives in a clay loam. Furthermore, the model was able to simulate well the impacts on ASW of the diverse irrigation schedules applied by the farmers along the different seasons.
The goodness-of-fit indicators relative to the ISAREG model test using farmers’ field observations of ASW are presented in Table 5. Results confirmed the good accuracy of the model in predicting the observed ASW after proper model parameterization. For maize crops, the model shows no tendency for under- or over-estimation of ASW with b0 close to 1.00 and quite low PBIAS (0 < PBIAS < 1.5%). Results also show that the estimated ASW values were statistically close to the measured ones. The R2 is high, indicating that the model could explain most of the observed variance. The errors of estimates are quite small: RMSE ranging from approximately 3 to 7 mm and the NRMSEs are lower than 6%. Finally, the EF values are higher than 0.76; thus, indicating that the variance of residuals is much lower than the variance of the measured ASW data. Similar results were obtained for maize with other applications of SWB models in Portugal (e.g., [76,77]).
Very good modeling results were obtained using the two data sets of ASWOBS in tomato (Table 5) with b0 close to 1.0 and low PBIAS values; thus, the model shows no trend for under- or over-estimation. R2 values range from 0.87 to 0.90; thus, higher than those obtained for maize. Estimation errors were low, with RMSE close to 6 mm and NRMSE lower than 4%. EF results are high (EF > 0.86) confirming that the variance of residuals is much lower than the variance of measured data. These results are comparable with those in the literature for the Mediterranean (e.g., [78]).
For the olive grove, the model was run for the growing and non-growing seasons showing very good accuracy of the model along the years (Figure 14b). Results in Table 5 show just a slight tendency for the model to over-estimate the ASWOBS by 2018 (b0 = 1.03 and PBIAS = 4%) while in the following years the model showed no tendency to under- or over-estimate. Errors of estimates were also relatively low with RMSE < 19 mm and representing less than 8% of the mean of the ASWOBS. The EF values are generally high. Few studies are available in literature relative to the use of SWB models for predicting ASW in olive orchards. The results of the current study are better than those reported using the WABOL model for an olive orchard in Spain [79].
Overall, results show that the parameterization of the ISAREG model was adequate and that the model was able to simulate the soil water storage along the seasons for the three studied crops. Thus, it can be assumed that the model may be further explored for estimating the soil water balance for the other crops present in the CIS.

3.5.2. Irrigation Water Requirements

The spatio-temporal distribution of crop irrigation water requirements (IWR), computed for each HUA (Figure 13), is presented in Figure 15, while Figure 16a,b shows the IWR per crop and per year, respectively.
The crop with the highest IWR was paddy rice due to the use of flooding irrigation, with lower application efficiencies than the other systems resulting from the large evaporation losses. Regarding the crops irrigated by pressurized systems, maize presents the highest IWR, with values above 750 mm during the 2017 and 2019 seasons, while lower values of IWR were estimated for 2018 and 2020. The higher IWR by 2017 and 2019 relates to the occurrence of less precipitation during the spring combined with higher ETo (Figure 16). The relatively large amplitudes of the boxes associated with each crop type (Figure 16a) are related to crop management practices (e.g., plant density, varieties) as well as the cultivation in different soil classes, which influences the water balance results. The second highest IWR are those for fodder crops, followed by tomato for industry and vegetable crops. The high IWR of fodder crops are, in part, associated with the lower efficiency of the traveling gun irrigation system, as compared with the center pivot and drip irrigation (see Section 2.3). For the lowest IWR, there are the olive groves, averaging 275 mm, and legume crops averaging 180 mm. The IWR are also influenced by the soil depth that constraints the soil available water for the crops. In the case of spring–summer vegetable crops, IWR showed a reduced spatial variability, which may be due to the shallow rooting system, so the depth of the soil does not impose limitations on the root system and therefore does not affect the IWR. Differently, in the case of maize and olive groves, there was a greater spatial variation in the IWR. In the olive grove, this spatial variation was most noticeable in 2017, because it was a precipitation scarce year. Thus, in the cases that orchards were installed in shallower soils, the IWR were higher, and irrigation started by early March.
Figure 16b shows the IWR for each irrigation campaign, considering all crops together. The higher IWR in 2017 relates to both the available water from precipitation as well as the climatic demand (ETo). On one hand, during the spring of 2017, precipitation was lower than the historical average for the months of April and May, when the spring crops are planted (Figure 2). On the other hand, the climatic demand (ETo) was higher than the average for April, June, and September, increasing crop evapotranspiration. During the winter of 2019, precipitation was very low, leading to a small amount of water storage in the soil profile at the time of planting in the spring, therefore, the irrigation season began earlier, as in 2017, leading to higher irrigation requirements when compared with 2018 and 2020.

3.5.3. Deep Percolation

The results relative to the deep percolation associated with the irrigation events for the set of selected crops is shown in Figure 16c. The null values correspond to the dry years of 2017 and 2019. Results show that the median deep percolation value is lower than the average value (right skewed) which results from a longer tail in the high deep percolation values relative to the low ones. The relatively large amplitude of the boxes is due to the different classes of soil suitability for irrigation and, in the case of legume crops, due to the different irrigation systems used (center pivot, travelling gun, and drip irrigation) and respective application efficiencies. For the same crop, soil class, and year, the HUA with a traveling gun system presented the highest deep percolation values, 32% higher than for those using drip irrigation systems.
Considering all the years, crops and soil classes, deep percolation associated with irrigation averages 11 ± 3% of the total deep percolation that occurred during the year.
Crops associated with higher values of deep percolation resulting from irrigation are maize and vegetable crops, and at the same time, these crops present the largest variability of the results when all years and soil classes are considered. Farmers tend to apply more water to these crops since they are the ones presenting the highest economic returns. The opposite is found for sunflower and fodder crops. The irrigation schemes for olive groves are characterized by a period of low or null application of water, during summer dormancy, but prior to this period, farmers tend to irrigate more and even beyond soil field capacity and therefore some deep percolation occurs [67,80].
Regarding land suitability for irrigation, the classes with higher amounts of deep percolation associated with the irrigation events are III and excluded, due to the shallower depth of the soil profile and lower water retention capacity (Figure 16d).
The spatial variability of deep percolation that occurs during the periods when the crop is in the field is shown in Figure 17a. Overall, during the crop season, the sum of irrigation and precipitation are responsible for 29 ± 13% of the deep percolation, while the precipitation that occurred during the spring and autumn (off-season), typical of the Mediterranean climate was responsible for the higher share of total deep percolation. The spatialization of the deep percolation calculated for each HAU shows that there are some hotspots for deep percolation during the studied years. In this study, it was considered that deep percolation hotspots were associated with HUA with deep percolation higher than 200 mm during the crop season. The 2016/2017 season presented a low hotspot area corresponding to 10% of the total irrigated area, although 26% of the irrigated area presents deep percolation higher than 150 mm during the crop-growing season. Differently during the 2017/2018 season, the deep percolation hotspot area increased to 21%. Results show that 55% of these hotspot areas corresponding HUAs include soils with land suitability to irrigation from classes exc and III, mainly due to shallow depth and lower water retention capacity. In these HUA the crop is not a differentiating factor. However, 26 and 18% of the HUAs include classes II and I, classified as good and excellent, respectively.
In this case, crops with higher irrigation water requirements and higher economic return as maize and spring vegetables are predominant. The 2018/2019 season presents a low deep percolation hotspot area of 6%, like that for 2016/2017. During 2019/2020, 10% of the area is a deep percolation hotspot, from which 56% present suitability classes exc and III and all types of crops, while for the 44% of hotspots on soil classes I and II the crops are maize and pastures. On average, irrigation was responsible for 26 ± 8% of deep percolation during the crop seasons.
Figure 17b presents the deep percolation that occurred during the off-season, after crop harvesting; thus, resulting from precipitation only. The temporal pattern shows the inter-annual variation in precipitation, which in the case of the humid year of 2018 resulted in almost the entire CIS area being a hotspot. Spatially, the areas with lower deep percolation correspond to land suitability classes very good and excellent, i.e., classes I and II. On average, during 2018, deep percolation represented 54 ± 6% of the average annual value. High deep percolation during autumn–winter is favorable since it recharges the Caia aquifer. However, the Caia Irrigation Scheme is within a designated vulnerable zone, overlaying a nitrate-contaminated aquifer. Therefore, if a large amount of N remains in the soil profile after the harvest of the spring–summer crops, it will be leached to the aquifer with the first autumn precipitation. The soil water balance model predicted an amount of deep percolation associated to precipitation of about 23 ± 6% of the total precipitation, which agrees with other studies performed in the area [6].
Deep percolation in the referred hotspots leads to low efficiency of irrigation. Furthermore, it affects the aquifer quality, since drainage water is responsible for the transport of nitrates resulting from fertilization and crop residue mineralization.

4. Conclusions

An integrated methodology was developed to estimate the spatio-temporal dynamics of deep percolation with the farm information and crop coefficient curves being derived from high-resolution satellite images processed on the Google Engine platform.
The identification and characterization of crop cycles, fundamental for the quantification of crop irrigation water requirements, was based on the treatment of several high-resolution Sentinel-2 temporal images on the Google Earth Engine platform. In a first analysis of the NDVI temporal profiles extracted for each plot, in the GEE platform, the atypical profiles for each crop were identified and discarded. This method also allowed improvement of the identification of the crop growth stages and the planting dates at the different plots. This information was used to feed the soil water balance model, ISAREG, used to compute deep percolation below the root zone, for each homogenous unit of analysis.
The use of the proposed integrated methodology showed a good accuracy when it was validated using observed soil water storage data.
Results from the spatio-temporal analysis of the computed deep percolation along the years allowed to answer the question raised at the beginning of the present study which was: are there hotspot areas in the Caia Irrigation Scheme with high deep percolation due to irrigation practices? The results show that deep percolation during the crop season was quite variable among years, with average values ranging between 77.3 and 119.2 mm. However, these values present a very high spatial variability, ranging from 0 to more than 200 mm. Deep percolation was higher during the 2017/2018 campaign mainly resulting from precipitation due to the occurrence of a high amount of rainfall in a low number of events, which in parallel led to reduced irrigation water requirements. Deep percolation associated with irrigation averaged 11 ± 3% of total deep percolation. On one hand, the deep percolation hotspots associated with the irrigation events correspond to soils with low land suitability for irrigation, independently of the crop, probably due to the shallower depth of the soil profile and median low water retention capacity. On the other hand, maize and spring vegetables stand out when the hotspots correspond to land classes I and II.
Overall, during the crop season, irrigation plus precipitation was responsible for 29 ± 13% of the deep percolation, while the precipitation that occurred during the spring and autumn, typical of the Mediterranean climate, was responsible for the higher share. On average, during the off-season period, deep percolation corresponded to 54 ± 6% of the average annual precipitation. Thus, aiming at controlling the nitrates leaching in this vulnerable area, farmers need to have appropriate support to manage the irrigation scheduling according to the soil, crop, and weather conditions as well as to manage fertilization.

Author Contributions

Conceptualization, M.d.R.C. and J.R.; methodology, A.F., M.d.R.C., J.R. and P.P.; data curation, A.F.; writing—original draft preparation, A.F. and M.d.R.C.; writing—review and editing, P.P. and J.R.; project administration and funding acquisition, M.d.R.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by HubIS Project (PRIMA/0006/2019). The support of the Portuguese Foundation for Science and Technology (FCT) through the grants attributed to the research Unit UID/AGR/04129/2020 (LEAF) LEAF (UID/AGR/04129/2019), to J. Rolim (DL 57/2016/CP1382/CT0021), and to P. Paredes (DL 57/2016/CP1382/CT0022) is acknowledged.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

The Centro Operativo de Tecnologias do Regadio (COTR) is acknowledged for providing the meteorological data. The Portuguese Association of Precision Agriculture (APAP) is acknowledged for providing the soil moisture data measured in the farmers’ fields.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Crop stage duration and seeding/planting dates obtained from the literature.
Table A1. Crop stage duration and seeding/planting dates obtained from the literature.
CropCrop Stage Duration (Days)Planting DateLiterature
InitialDevelopmentMidLate
Garlic7030472701/Dec[6,81]
Broccoli3545401501/Nov.[47]
Chickpea2030351515/Mar
Winter cereals30140402015/Oct
Olive groves3090609001/Mar
Orchards30501303001/Apr
Vineyards2050756001/Mar
Almonds32651404123/Jan[82]
Table A2. Actual crop coefficients.
Table A2. Actual crop coefficients.
Kc act iniKc act midKc act end
201720182019202020172018201920202017201820192020
Garlic0.70---1.00---0.65- -
Tomato0.50-0.500.501.10-1.041.050.90-0.840.85
Melon0.40---0.99---0.79- -
Sweet peppers0.500.50.50.51.091.031.051.051.040.981.001.00
Broccoli0.400.40.400.401.040.940.961.001.040.94 1.00
Sunflower0.350.350.350.351.191.131.151.150.340.30.30.30
Rapeseed0.350.350.350.351.031.021.051.000.280.350.30.25
Chickpeas-0.40.400.40-0.980.991.00-0.350.340.35
Rice1.051.051.051.051.201.21.201.201.091.051.051.05
Maize grain0.300.30.300.301.201.21.191.200.360.30,390.30
Maize silage and sorghum0.300.30.300.301.201.061.141.151.010.860.940.95
Winter cereals0.300.30.300.301.091.081.111.060.340.330.360.31
Fodder crop multiple cuts0.400.40.400.400.990.890.910.950.940.840.860.90
Pastures0.750.750.750.750.750.750.750.750.750.750.750.75
Almond0.400.40.400.400.920.820.840.850.670.570.590.60
Orchards0.450.450.450.451.060.970.991.000.760.670.690.70
Vineyards0.300.30.300.300.750.670.690.700.600.520.540.55
Table A3. Description of the homogeneous units of analysis for the 2019/2020 season.
Table A3. Description of the homogeneous units of analysis for the 2019/2020 season.
IDSoil TypeIrrigation SystemCropIDSoil TypeIrrigation SystemCrop
0Exc.n.a.n.a.55IIDripChickpeas and peas
1Exc.FallowFallow56IIDripWinter vegetables
2Exc.DripAlmond, walnut, and pistachio57IIDripSummer vegetables
3Exc.DripChickpeas and peas58IIDripOlive groves
4Exc.DripWinter vegetables59IIDripOther crops
5Exc.DripSummer vegetables60IIDripOrchards
6Exc.DripOlive groves61IIDripTomato
7Exc.DripOther crops62IIDripVineyards
8Exc.DripPastures63IITravelling gunChickpeas and peas
9Exc.DripOrchards64IICenter pivotWinter cereals
10Exc.DripTomato65IICenter pivotRapeseed
11Exc.DripVineyards66IICenter pivotFodder crops
12Exc.Travelling gunChickpeas and peas67IICenter pivotSunflower
13Exc.Center pivotWinter cereals68IICenter pivotChickpeas and peas
14Exc.Center pivotRapeseed69IICenter pivotWinter vegetables
15Exc.Center pivotFodder crops with multiple cuts70IICenter pivotSummer vegetables
16Exc.Center pivotFodder crops71IICenter pivotMaize grain
17Exc.Center pivotSunflower72IICenter pivotMaize silage and sorghum
18Exc.Center pivotChickpeas and peas73IICenter pivotPastures
19Exc.Center pivotWinter vegetables74IIFlooded paddiesRice
20Exc.Center pivotMaize grain75IIRainfed cropsFodder crops
21Exc.Center pivotMaize silage and sorghum76IIRainfed cropsMaize grain
22Exc.Center pivotPastures77IIRainfed cropsOlive groves
23Exc.Flooded paddiesRice78IIIn.a.n.a.
24Exc.Rainfed cropsFodder crops79IIIFallowFallow
25Exc.Rainfed cropsAlmond, walnut, and pistachio80IIIDripAlmond, walnut, and pistachio
26Exc.Rainfed cropsMaize grain81IIIDripChickpeas and peas
27Exc.Rainfed cropsOlive groves82IIIDripWinter vegetables
28Exc.Rainfed cropsOther crops83IIIDripSummer vegetables
29Exc.Rainfed cropsVineyards84IIIDripMaize grain
30In.a.n.a.85IIIDripOlive groves
31IFallowFallow86IIIDripOther crops
32IDripAlmond, walnut, and pistachio87IIICenter pivotPastures
33IDripChickpeas and peas88IIIDripOrchards
34IDripSummer vegetables89IIIDripTomato
35IDripOlive groves90IIIDripVineyards
36IDripOther crops91IIITravelling gunChickpeas and peas
37IDripOrchards92IIICenter pivotWinter cereals
38IDripTomato93IIICenter pivotRapeseed
39IDripVineyards94IIICenter pivotFodder crops with multiple cuts
40ITravelling gunChickpeas and peas95IIICenter pivotFodder crops
41ICenter pivotWinter cereals96IIICenter pivotSunflower
42ICenter pivotFodder crops with multiple cuts97IIICenter pivotChickpeas and peas
43ICenter pivotFodder crops98IIICenter pivotWinter vegetables
44ICenter pivotSunflower99IIICenter pivotMaize grain
45ICenter pivotSummer vegetables100IIICenter pivotMaize silage and sorghum
46ICenter pivotMaize grain101IIICenter pivotPastures
47ICenter pivotMaize silage and sorghum102IIIFlooded paddiesRice
48ICenter pivotPastures103IIIRainfed cropsFodder crops
49IFlooded paddiesRice104IIIRainfed cropsAlmond, walnut, and pistachio
50IRainfed cropsFodder crops105IIIRainfed cropsMaize grain
51IRainfed cropsOlive groves106IIIRainfed cropsOlive groves
52IIn.a.n.a.107IIIRainfed cropsOther crops
53IIFallowFallow108IIIRainfed cropsVineyards
54IIDripAlmond, walnut, and pistachio
Note: n.a.: not applicable since it concerns non-agricultural areas.

References

  1. Arauzo, M.; Martínez-Bastida, J.J. Environmental factors affecting diffuse nitrate pollution in the major aquifers of central Spain: Groundwater vulnerability vs. groundwater pollution. Environ. Earth Sci. 2015, 73, 8271–8286. [Google Scholar] [CrossRef] [Green Version]
  2. Ferreira, J.P.L.; Chachadi, A.G.; Diamantino, C.; Henriques, M.J. Assessing aquifer vulnerability to seawater intrusion using GALDIT method: Part 1-Application to the Portuguese Aquifer of Monte Gordo. In Water in Celtic Countries: Quantity, Quality and Climate Variability, Proceedings of the Fourth Inter Colloquium on Hydrology and Management of Water Resources, Guimares, Portugal, 11–13 July 2005; Lobo Ferreira, J.P., Viera, J.M.P., Eds.; IAHS Press: Wallingford, UK, 2007; pp. 161–171. ISBN 978-1-901502-88-6. [Google Scholar]
  3. Stigter, T.Y.; Nunes, J.P.; Pisani, B.; Fakir, Y.; Hugman, R.; Li, Y.; Tomé, S.; Ribeiro, L.; Samper, J.; Oliveira, R.; et al. Comparative assessment of climate change and its impacts on three coastal aquifers in the Mediterranean. Reg. Environ. Chang. 2014, 14, 41–56. [Google Scholar] [CrossRef]
  4. Carneiro, J.; Coutinho, J.; Trindade, H. Nitrate leaching from a maize×oats double-cropping forage system fertilized with organic residues under Mediterranean conditions. Agric. Ecosyst. Environ. 2012, 160, 29–39. [Google Scholar] [CrossRef]
  5. Poch-Massegú, R.; Jiménez-Martínez, J.; Wallis, K.; de Cartagena, F.R.; Candela, L. Irrigation return flow and nitrate leaching under different crops and irrigation methods in Western Mediterranean weather conditions. Agric. Water Manag. 2014, 134, 1–13. [Google Scholar] [CrossRef]
  6. Almeida, C.; Mendonça, J.J.L.; Jesus, M.R.; Gomes, A.J. Sistemas Aquíferos de Portugal Continental; Centro de Geologia & Instituto da Água: Lisboa, Portugal, 2000; Volume 3, pp. 432–661. [Google Scholar]
  7. Pereira, L.S.; Paredes, P.; López-Urrea, R.; Hunsaker, D.J.; Mota, M.; Shad, Z.M. Standard single and basal crop coefficients for vegetable crops, an update of FAO56 crop water requirements approach. Agric. Water Manag. 2020, 243, 106196. [Google Scholar] [CrossRef]
  8. Pereira, L.S.; Paredes, P.; Hunsaker, D.J.; López-Urrea, R.; Shad, Z.M. Standard single and basal crop coefficients for field crops. Updates and advances to the FAO56 crop water requirements method. Agric. Water Manag. 2021, 243, 106466. [Google Scholar] [CrossRef]
  9. Xiao, Y.; Zhan, Q. A Review of Remote Sensing Applications in Urban Planning and Management in China; Joint Urban Remote Sensing Event: Shanghai, China, 2009; pp. 1–5. [Google Scholar] [CrossRef]
  10. Pôças, I.; Calera, A.; Campos, I.; Cunha, M. Remote sensing for estimating and mapping single and basal crop coefficients: A review on spectral vegetation indices approaches. Agric. Water Manag. 2020, 233, 106081. [Google Scholar] [CrossRef]
  11. Diao, C. Remote sensing phenological monitoring framework to characterize corn and soybean physiological growing stages. Remote Sens. Environ. 2020, 248, 111960. [Google Scholar] [CrossRef]
  12. Gao, F.; Zhang, X. Mapping crop phenology in near real-time using satellite remote sensing: Challenges and opportunities. J. Remote Sens. 2021, 2021, 8379391. [Google Scholar] [CrossRef]
  13. Navarro, A.; Rolim, J.; Miguel, I.; Catalão, J.; Silva, J.; Painho, M.; Vekerdy, Z. Crop monitoring based on SPOT-5 Take-5 and sentinel-1A data for the estimation of crop water requirements. Remote Sens. 2016, 8, 525. [Google Scholar] [CrossRef] [Green Version]
  14. Rolim, J.; Navarro, A.; Vilar, P.; Saraiva, C.; Catalao, J. Crop data retrieval using earth observation data to support agricultural water management. Eng. Agric. 2019, 39, 381–390. [Google Scholar] [CrossRef]
  15. Ballesteros, R.; Moreno, M.A.; Barroso, F.; González-Gómez, L.; Ortega, J.F. Assessment of Maize Growth and Development with High-and Medium-Resolution Remote Sensing Products. Agronomy 2021, 11, 940. [Google Scholar] [CrossRef]
  16. Mahlayeye, M.; Darvishzadeh, R.; Nelson, A. Cropping Patterns of Annual Crops: A Remote Sensing Review. Remote Sens. 2022, 14, 2404. [Google Scholar] [CrossRef]
  17. European Space Agency Observing the Earth. Available online: http://www.esa.int/Applications/Observing_the_Earth/ESA_for_Earth (accessed on 29 May 2021).
  18. Plataforma Google Earth Engine. Available online: https://code.earthengine.google.com/ (accessed on 15 July 2021).
  19. Mulla, D.J. Twenty five years of remote sensing in precision agriculture: Key advances and remaining knowledge gaps. Biosyst. Eng. 2013, 11, 358–371. [Google Scholar] [CrossRef]
  20. Belmonte, A.C.; Jochum, A.M.; García, A.C.; Rodríguez, A.M.; Fuster, P.L. Irrigation management from space: Towards user-friendly products. Irrig. Drain. Syst. 2005, 19, 337–353. [Google Scholar] [CrossRef]
  21. D’Urso, G.; Belmonte, A.C. Operative approaches to determine crop water requirements from earth observation data: Methodologies and applications. AIP Conf. Proc. 2006, 852, 14–25. [Google Scholar] [CrossRef]
  22. Bégué, A.; Arvor, D.; Bellon, B.; Betbeder, J.; de Abelleyra, D.; Ferraz, R.P.D.; Lebourgeois, V.; Lelong, C.; Simões, M.; Verón, S.R. Remote sensing and cropping practices: A review. Remote Sens. 2018, 10, 99. [Google Scholar] [CrossRef] [Green Version]
  23. Rouse, W.; Haas, R.; Scheel, J.; Deering, W. Monitoring Vegetation Systems in Great Plains with ERST, Proceedings of the Third ERTS Symposium, NASA SP-351; US Government Printing Office: Washington, DC, USA, 1973; pp. 309–317.
  24. Zeng, L.; Wardlow, B.D.; Xiang, D.; Hu, S.; Li, D. A review of vegetation phenological metrics extraction using time-series, multispectral satellite data. Remote Sens. Environ. 2020, 237, 111511. [Google Scholar] [CrossRef]
  25. González-Piqueras, J. Evapotranspiration de la Cubierta Vegetal mediante la Determinación del Coeficiente de Cultivo por Teledetección. Extensión a Escala Regional: Aquífero 08.29 Mancha Oriental. Ph.D. Thesis, Universitat de Valencia, Valencia, Spain, 2006. Available online: http://hdl.handle.net/10803/10340 (accessed on 29 May 2021).
  26. D’Urso, G.; Richter, K.; Calera, A.; Osann, M.A.; Escadafal, R.; Garatuza-Pajan, J.; Vuolo, F. Earth Observation products for operational irrigation management in the context of the PLEIADeS project. Agric. Water Manag. 2010, 98, 271–282. [Google Scholar] [CrossRef]
  27. Toureiro, C.; Serralheiro, R.; Shahidian, S.; Sousa, A. Irrigation management with remote sensing: Evaluating irrigation requirement for maize under Mediterranean climate condition. Agric. Water Manag. 2017, 184, 211–220. [Google Scholar] [CrossRef]
  28. Consoli, S.; Vanella, D. Mapping crop evapotranspiration by integrating vegetation indices into a soil water balance model. Agric. Water Manag. 2014, 143, 71–81. [Google Scholar] [CrossRef]
  29. Odi-Lara, M.; Campos, I.; Neale, C.M.U.; Ortega-Farías, S.; Poblete-Echeverría, C.; Balbontín, C.; Calera, A. Estimating evapotranspiration of an apple orchard using a remote sensing-based soil water balance. Remote Sens. 2016, 8, 253. [Google Scholar] [CrossRef] [Green Version]
  30. Awada, H.; Di Prima, S.; Sirca, C.; Giadrossich, F.; Marras, S.; Spano, D.; Pirastru, M. A remote sensing and mod-eling integrated approach for constructing continuous time series of daily actual evapotranspiration. Agric. Water Manag. 2022, 260, 107320. [Google Scholar] [CrossRef]
  31. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  32. Mutanga, O.; Kumar, L. Google Earth Engine Applications. Remote Sens. 2019, 11, 591. [Google Scholar] [CrossRef] [Green Version]
  33. Amani, M.; Ghorbanian, A.; Ahmadi, S.A.; Kakooei, M.; Moghimi, A.; Mirmazloumi, S.M.; Moghaddam, S.H.A.; Mahdavi, S.; Ghahremanloo, M.; Parsian, S.; et al. Google earth engine cloud computing platform for remote sensing big data applications: A comprehensive review. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 5326–5350. [Google Scholar] [CrossRef]
  34. Mhawej, M.; Faour, G. Open-source Google Earth Engine 30-m evapotranspiration rates retrieval: The SEBALIGEE system. Environ. Model. Softw. 2020, 133, 104845. [Google Scholar] [CrossRef]
  35. Kilic, A.; Allen, R.G.; Blankenau, P.; Revelle, P.; Ozturk, D.; Huntington, J. Global production and free access to Landsat-scale Evapotranspiration with EEFlux and eeMETRIC. In Proceedings of the 6th Decennial National Irrigation Symposium, American Society of Agricultural and Biological Engineers, San Diego, CA, USA, 6–8 December 2021; pp. 2020–2038. [Google Scholar] [CrossRef]
  36. Wu, F.; Wu, B.; Zhu, W.; Yan, N.; Ma, Z.; Wang, L.; Lu, Y.; Xu, J. ETWatch Cloud: APIs for regional actual evapotranspiration data. Environ. Model. Softw. 2021, 145, 105174. [Google Scholar] [CrossRef]
  37. He, M.; Kimball, J.S.; Maneta, M.P.; Maxwell, B.D.; Moreno, A.; Beguería, S.; Wu, X. Regional crop gross primary productivity and yield estimation using fused landsat-MODIS data. Remote Sens. 2018, 10, 372. [Google Scholar] [CrossRef] [Green Version]
  38. Mandal, D.; Kumar, V.; Bhattacharya, A.; Rao, Y.S.; Siqueira, P.; Bera, S. Sen4Rice: A processing chain for differentiating early and late transplanted rice using time-series sentinel-1 SAR data with google earth engine. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1947–1951. [Google Scholar] [CrossRef]
  39. Xiong, J.; Thenkabail, P.S.; Tilton, J.C.; Gumma, M.K.; Teluguntla, P.; Oliphant, A.; Congalton, R.G.; Yadav, K.; Gorelick, N. Nominal 30-m cropland extent map of continental Africa by integrating pixel-based and object-based algorithms using Sentinel-2 and Landsat-8 data on google earth engine. Remote Sens. 2017, 9, 1065. [Google Scholar] [CrossRef] [Green Version]
  40. Aguilar, R.; Zurita-Milla, R.; Izquierdo-Verdiguier, E.; de By, R.A. A cloud-based multi-temporal ensemble classifier to map smallholder farming systems. Remote Sens. 2018, 10, 729. [Google Scholar] [CrossRef] [Green Version]
  41. Liu, L.; Xiao, X.; Qin, Y.; Wang, J.; Xu, X.; Hu, Y.; Qiao, Z. Mapping cropping intensity in China using time series Landsat and Sentinel-2 images and Google Earth Engine. Remote Sens. Environ. 2020, 239, 111624. [Google Scholar] [CrossRef]
  42. Calera, A.; Campos, I.; Osann, A.; D’Urso, G.; Menenti, M.; Calera, A.; Campos, I.; Osann, A.; D’Urso, G.; Menenti, M. Remote Sensing for Crop Water Management: From ET Modelling to Services for the End Users. Sensors 2017, 17, 1104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Venancio, L.P.; Eugenio, F.C.; Filgueiras, R.; Da Cunha, F.F.; Dos Santos, R.A.; Ribeiro, W.R.; Mantovani, E.C. Mapping within-field variability of soybean evapotranspiration and crop coefficient using the Earth Engine Evaporation Flux (EEFlux) application. PLoS ONE 2020, 15, e0235620. [Google Scholar] [CrossRef] [PubMed]
  44. Melton, F.S.; Huntington, J.L.; Grimm, R.; Herring, J.; Hall, M.; Rollison, D.; Erickson, T.; Allen, R.G.; Anderson, M.; Fisher, J.B.; et al. 2021: OpenET: Filling a critical data gap in water management for the western United States. J. Am. Water Resour. Assoc. 2021, 1–24. [Google Scholar] [CrossRef]
  45. Todorovic, M.; Steduto, P. A GIS for irrigation management. Phys. Chem. Earth Parts A/B/C 2003, 28, 163–174. [Google Scholar] [CrossRef]
  46. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World Map of the Köppen-Geiger climate classification updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef]
  47. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements. Irrig. Drain. 1998, 56, 300. [Google Scholar]
  48. FAO. International Soil Classification System for Naming Soils and Creating Legends for Soil Maps. In World Reference Base for Soil Resources; FAO: Roma, Italy, 2014. [Google Scholar]
  49. Associação Beneficiários do Caia. Report 2013. Available online: http://www.abcaia.pt/index.php/associacao/relatorio-abcaia/relatorio-2013 (accessed on 16 April 2021).
  50. Teixeira, J.L.; Pereira, L.S. ISAREG, an irrigation scheduling model. ICID Bull. 1992, 41, 29–48. [Google Scholar]
  51. Rallo, G.; Paço, T.A.; Paredes, P.; Puig-Sirera, À.; Massai, R.; Provenzano, G.; Pereira, L.S. Updated single and dual crop coefficients for tree and vine fruit crops. Agric. Water Manag. 2021, 250, 106645. [Google Scholar] [CrossRef]
  52. Pereira, L.S. Necessidades de Água e Métodos de Rega; Publicaçöes Europa-América: Lisboa, Portugal, 2004; p. 313. ISBN 5601072370609. [Google Scholar]
  53. Keller, J.; Bliesner, R.D. Sprinkle and Trickle Irrigation; The Blackburn Press: Caldwell, NJ, USA, 2000; ISBN 1-930665-19-9. [Google Scholar]
  54. Allen, R.G.; Wright, J.L.; Pruitt, W.O.; Pereira, L.S.; Jensen, M.E. Water requirements. In Design and Operation of Farm Irrigation Systems, 2nd ed.; Hoffman, G.J., Evans, R.G., Jensen, M.E., Martin, D.L., Elliot, R.L., Eds.; ASABE: St. Joseph, MI, USA, 2007; pp. 208–288. [Google Scholar]
  55. Renard, K.G. Predicting Soil Erosion by Water: A Guide to conservation planning with the Revised Universal Soil Loss Equation (RUSLE); Agriculture Research Service, United States Government Printing: Washington, DC, USA, 1997; 384p.
  56. Cholpankulov, E.D.; Inchenkova, O.P.; Paredes, P.; Pereira, L.S. Cotton irrigation scheduling in central Asia: Model calibration and validation with consideration of groundwater contribution. Irrig. Drain. 2008, 57, 516–532. [Google Scholar] [CrossRef] [Green Version]
  57. Stulina, G.; Cameira, M.R.; Pereira, L.S. Using RZWQM to search improved practices for irrigated maize in Fergana, Uzbekistan. Agric. Water Manag. 2005, 77, 263–281. [Google Scholar] [CrossRef]
  58. Popova, Z.; Pereira, L.S. Modelling for maize irrigation scheduling using long term experimental data from Plovdiv region, Bulgaria. Agric. Water Manag. 2011, 98, 675–683. [Google Scholar] [CrossRef]
  59. Wu, Y.; Liu, T.; Paredes, P.; Duan, L.; Pereira, L.S. Water use by a groundwater dependent maize in a semi-arid region of Inner Mongolia: Evapotranspiration partitioning and capillary rise. Agric. Water Manag. 2015, 152, 222–232. [Google Scholar] [CrossRef] [Green Version]
  60. Sousa, V.; Pereira, L.S. Regional analysis of irrigation water requirements using kriging: Application to potato crop (Solanum tuberosum L.) at Trás-os-Montes. Agric. Water Manag. 1999, 40, 221–233. [Google Scholar] [CrossRef]
  61. Chaterlán, Y.; León, M.; Duarte, C.; López, T.; Paredes, P.; Pereira, L.S. Determination of crop coefficients for horticultural crops in Cuba through field experiments and water balance simulation. Acta Hortic. 2011, 889, 475–482. [Google Scholar] [CrossRef]
  62. Cancela, J.J.; Cuesta, T.S.; Neira, X.X.; Pereira, L.S. Modelling for improved irrigation water management in a temperate region of Northern Spain. Biosyst. Eng. 2008, 99, 587–597. [Google Scholar] [CrossRef]
  63. Wu, Y.; Liu, T.; Paredes, P.; Duan, L.; Wang, H.; Wang, T.; Pereira, L.S. Ecohydrology of groundwater-dependent grasslands of the semi-arid Horqin sandy land of Inner Mongolia focusing on evapotranspiration partition. Ecohydrology 2016, 9, 1052–1067. [Google Scholar] [CrossRef]
  64. Alba, I.; Rodrigues, P.N.; Pereira, L.S. Irrigation scheduling simulation for citrus in Sicily to cope with water scarcity. In Tools for Drought Mitigation in Mediterranean Regions; Rossi, G., Cancelliere, A., Pereira, L.S., Oweis, T., Shatanawi, M., Zairi, A., Eds.; Kluwer: Dordrecht, The Netherlands, 2003; pp. 223–242. [Google Scholar] [CrossRef]
  65. Chaterlán, Y.; Hernández, G.; López, T.; Martínez, R.; Puig, O.; Paredes, P.; Pereira, L.S. Estimation of the papaya crop coefficients for improving irrigation water management in south of Havana. Acta Hortic. 2012, 928, 179–186. [Google Scholar] [CrossRef]
  66. Valverde, P.; Serralheiro, R.; de Carvalho, M.; Maia, R.I.; Oliveira, B.; Ramos, V. Climate change impacts on irrigated agriculture in the Guadiana River basin (Portugal). Agric. Water Manag. 2015, 152, 17–30. [Google Scholar] [CrossRef] [Green Version]
  67. Branquinho, S.; Rolim, J.; Teixeira, J.L. Climate Change Adaptation Measures in the Irrigation of a Super-Intensive Olive Orchard in the South of Portugal. Agronomy 2021, 11, 1658. [Google Scholar] [CrossRef]
  68. Zaccaria, D.; Oueslati, I.; Neale, C.M.U.; Lamaddalena, N.; Vurro, M.; Pereira, L.S. Flexible delivery schedules to improve farm irrigation and reduce pressure on groundwater: A case study in southern Italy. Irrig. Sci. 2009, 28, 257–270. [Google Scholar] [CrossRef]
  69. Victoria, F.B.; Viegas Filho, J.S.; Pereira, L.S.; Teixeira, J.L.; Lanna, A.E. Multi-scale modeling for water resources planning and management in rural basins. Agric. Water Manag. 2005, 77, 4–20. [Google Scholar] [CrossRef]
  70. Eisenhauer, J.G. Regression through the origin. Teach. Stat. 2003, 25, 76–80. [Google Scholar] [CrossRef]
  71. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models: Part 1 A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  72. Legates, D.R.; McCabe, G.J. Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 1999, 35, 233–241. [Google Scholar] [CrossRef]
  73. Raposo, J.R. A Rega-Dos Primitivos Regadios às Modernas Técnicas de Rega; Fundação Calouste Gulbekian: Lisbon, Portugal, 1996. [Google Scholar]
  74. Cao, R.; Chen, Y.; Shen, M.; Chen, J.; Zhou, J.; Wang, C.; Yang, W. A simple method to improve the quality of NDVI time-series data by integrating spatiotemporal information with the Savitzky-Golay filter. Remote Sens. Environ. 2018, 217, 244–257. [Google Scholar] [CrossRef]
  75. Vilar, P.; Navarro, A.; Rolim, J. Utilização de Imagens de Deteção Remota para Monitorização das Culturas e Estimação das Necessidades de Rega. In Proceedings of the VIII Conferência Nacional de Cartografia e Geodesia, Amadora, Portugal, 29–30 November 2015; pp. 1–8. [Google Scholar]
  76. Rosa, R.D.; Paredes, P.; Rodrigues, G.C.; Fernando, R.M.; Alves, I.; Allen, R.G.; Pereira, L.S. Implementing the dual crop coefficient approach in interactive software: 2. Model testing. Agric. Water Manag. 2012, 103, 62–77. [Google Scholar] [CrossRef]
  77. Paredes, P.; Rodrigues, G.C.; Alves, I.; Pereira, L.S. Partitioning evapotranspiration, yield prediction and economic re-turns of maize under various irrigation management strategies. Agric. Water Manag. 2014, 135, 27–39. [Google Scholar] [CrossRef]
  78. Zairi, A.; El Amami, H.; Slatni, A.; Pereira, L.S.; Rodrigues, P.N.; Machado, T. Coping with drought: Deficit irrigation strategies for cereals and field vegetable crops in Central Tunisia. In Tools for Drought Mitigation in Mediterranean Regions; Rossi, G., Cancelliere, A., Pereira, L.S., Oweis, T., Shat-anawi, M., Zairi, A., Eds.; Kluwer: Dordrecht, The Netherlands, 2003; pp. 181–201. [Google Scholar]
  79. Abazi, U.; Lorite, I.J.; Cárceles, B.; Raya, A.M.; Durán, V.H.; Francia, J.R.; Gómez, J.A. WABOL: A conceptual water balance model for analyzing rainfall water use in olive orchards under different soil and cover crop management strategies. Comp. Electron. Agric. 2013, 91, 35–48. [Google Scholar] [CrossRef] [Green Version]
  80. Fernández, S.C.; Gallardo, J.R.; Mayorga, A.A.V. Fundamentos de Teledetección Espacial; Ediciones Rialp: Ciudad Real, Spain, 1996. [Google Scholar]
  81. Villalobos, F.J.; Testi, L.; Rizzalli, R.; Orgaz, F. Evapotranspiration and crop coefficients of irrigated garlic (Allium sativum L.) in a semi-arid climate. Agric. Water Manag. 2004, 64, 233–249. [Google Scholar] [CrossRef]
  82. Moita, R.A.D. Avaliação das Necessidades de rega de um Amendoal na Área de Influência do Alqueva. Master’s Thesis, Instituto Superior de Agronomia, Lisbon, Portugal, 2021. [Google Scholar]
Figure 1. Localization of the Caia Irrigation Scheme in Portugal and in Europe.
Figure 1. Localization of the Caia Irrigation Scheme in Portugal and in Europe.
Water 14 02324 g001
Figure 2. Monthly reference evapotranspiration (ETo) (markers) and precipitation (columns) for the four studied years.
Figure 2. Monthly reference evapotranspiration (ETo) (markers) and precipitation (columns) for the four studied years.
Water 14 02324 g002
Figure 3. Areas occupied by the different crops for the period 2017–2020 as a percentage of total cultivated area (OG—olive groves; WC—winter cereals; M—maize; O—orchards; PT—processing tomato; V—vegetables; FC—fodder crops; Vy—vineyards; PR—paddy rice; P—pastures; R—rapeseed).
Figure 3. Areas occupied by the different crops for the period 2017–2020 as a percentage of total cultivated area (OG—olive groves; WC—winter cereals; M—maize; O—orchards; PT—processing tomato; V—vegetables; FC—fodder crops; Vy—vineyards; PR—paddy rice; P—pastures; R—rapeseed).
Water 14 02324 g003
Figure 4. Flowchart showing the global methodology used for crop data processing.
Figure 4. Flowchart showing the global methodology used for crop data processing.
Water 14 02324 g004
Figure 5. Average NDVI for: (a) the period March–May; (b) the period May–August; (c) identification of the irrigated areas within the farmers’ fields.
Figure 5. Average NDVI for: (a) the period March–May; (b) the period May–August; (c) identification of the irrigated areas within the farmers’ fields.
Water 14 02324 g005
Figure 6. Flowchart with the modeling strategy used for the soil water balance and the deep percolation predictions (Kc is the crop coefficient, Zr is the rooting depth, p is the soil water depletion fraction for no stress, θFC is the soil water content at field capacity, θWP is the soil water content at permanent wilting point, Z is the soil depth, Bd is the soil bulk density, TAW is the total available soil water, DP is the deep percolation).
Figure 6. Flowchart with the modeling strategy used for the soil water balance and the deep percolation predictions (Kc is the crop coefficient, Zr is the rooting depth, p is the soil water depletion fraction for no stress, θFC is the soil water content at field capacity, θWP is the soil water content at permanent wilting point, Z is the soil depth, Bd is the soil bulk density, TAW is the total available soil water, DP is the deep percolation).
Water 14 02324 g006
Figure 7. Example of the NDVI temporal profiles for the main permanent crops in the Caia Irrigation Scheme for 2017.
Figure 7. Example of the NDVI temporal profiles for the main permanent crops in the Caia Irrigation Scheme for 2017.
Water 14 02324 g007
Figure 8. Example of the NDVI temporal profiles for annual winter crops and pastures for 2017.
Figure 8. Example of the NDVI temporal profiles for annual winter crops and pastures for 2017.
Water 14 02324 g008
Figure 9. NDVI time series profiles and respective actual crop coefficients (Kc act) for the main annual spring crops in the Caia Irrigation Scheme for 2017 (lines—FAO56 [47] Kc segmented curves; markers—clouds of Kc-NDVI values determined according to Calera-Belmonte et al. [20].
Figure 9. NDVI time series profiles and respective actual crop coefficients (Kc act) for the main annual spring crops in the Caia Irrigation Scheme for 2017 (lines—FAO56 [47] Kc segmented curves; markers—clouds of Kc-NDVI values determined according to Calera-Belmonte et al. [20].
Water 14 02324 g009
Figure 10. Spatio-temporal distribution of crops within the Caia Irrigation Scheme for the studied years of 2017–2020.
Figure 10. Spatio-temporal distribution of crops within the Caia Irrigation Scheme for the studied years of 2017–2020.
Water 14 02324 g010
Figure 11. Spatio-temporal distribution of the irrigation systems within the Caia Irrigation Scheme for 2017–2020.
Figure 11. Spatio-temporal distribution of the irrigation systems within the Caia Irrigation Scheme for 2017–2020.
Water 14 02324 g011
Figure 12. Classes of land suitability for irrigation within the Caia Irrigation Scheme (I—excellent; II—good; III—medium; exc—not suitable).
Figure 12. Classes of land suitability for irrigation within the Caia Irrigation Scheme (I—excellent; II—good; III—medium; exc—not suitable).
Water 14 02324 g012
Figure 13. Homogeneous units of analysis within the Caia Irrigation Scheme along the studied years.
Figure 13. Homogeneous units of analysis within the Caia Irrigation Scheme along the studied years.
Water 14 02324 g013
Figure 14. Comparison between soil water storage simulated and measured in the farmers’ fields along the studied seasons: (a) maize; (b) olive groves; and (c) processing tomato.
Figure 14. Comparison between soil water storage simulated and measured in the farmers’ fields along the studied seasons: (a) maize; (b) olive groves; and (c) processing tomato.
Water 14 02324 g014
Figure 15. Irrigation water requirements within the diverse homogeneous units of analysis (HUA).
Figure 15. Irrigation water requirements within the diverse homogeneous units of analysis (HUA).
Water 14 02324 g015
Figure 16. Box and wishers plots showing the lower quartile (Q1), the median (Q2), the upper quartile (Q3), the mean (×), and the minimum and maximum values (◦ outlier values are also depicted) relative to the 2017–2020 period for: (a) irrigation water requirements for the crops in the Caia Irrigation Scheme; (b) irrigation requirements for the different studied years; (c) deep percolation from irrigation per crop; (d) deep percolation from irrigation per class of land suitability for irrigation (Mz—maize; OG—olive groves; Sf—sunflower; Hsc—spring vegetable crops; To—tomato; Lc—legume crops; Fc—fodder crops; I, II, III, and EXC are classes of decreasing land suitability for irrigation).
Figure 16. Box and wishers plots showing the lower quartile (Q1), the median (Q2), the upper quartile (Q3), the mean (×), and the minimum and maximum values (◦ outlier values are also depicted) relative to the 2017–2020 period for: (a) irrigation water requirements for the crops in the Caia Irrigation Scheme; (b) irrigation requirements for the different studied years; (c) deep percolation from irrigation per crop; (d) deep percolation from irrigation per class of land suitability for irrigation (Mz—maize; OG—olive groves; Sf—sunflower; Hsc—spring vegetable crops; To—tomato; Lc—legume crops; Fc—fodder crops; I, II, III, and EXC are classes of decreasing land suitability for irrigation).
Water 14 02324 g016
Figure 17. Spatio-temporal distribution of deep percolation: (a) during the crop season; and (b) during off-crop season. The spatial aggregation into the irrigation scheme scale provided a method for earth-observation-based accounting of the irrigation water requirements, with interest for the water user’s association manager, and at the same time for the detection of water losses by deep percolation and of deep percolation hotspots (deep percolation higher than 200 mm) within the Caia Irrigation Scheme.
Figure 17. Spatio-temporal distribution of deep percolation: (a) during the crop season; and (b) during off-crop season. The spatial aggregation into the irrigation scheme scale provided a method for earth-observation-based accounting of the irrigation water requirements, with interest for the water user’s association manager, and at the same time for the detection of water losses by deep percolation and of deep percolation hotspots (deep percolation higher than 200 mm) within the Caia Irrigation Scheme.
Water 14 02324 g017aWater 14 02324 g017b
Table 1. Data sources used in the present study.
Table 1. Data sources used in the present study.
Data SetObservationsSource
SoilsSoil water holding capacity and textural characteristics; land use capacity Map 1:25,000Portuguese Soil Map (CSP) and Land Use Capacity (DGADR—Ministry of Agriculture)
WeatherDaily weather data (2002–2020): maximum and minimum temperatures (°C), maximum and minimum relative humidity (%), solar radiation (kJ·m−2·dia−1), wind speed (m·s−1), and precipitation (mm)Meteorological station of Elvas (38°54′56″ N, 7°05′56″ W, 202 m a.s.l) (COTR)
TopographyMilitary Map of Portugal (1:25,000)Instituto Geográfico do Exército
Hydrography
and altimetry
DTM (Digital Terrain Model) for the slope
Contour lines, elevation points, and water flow lines
Instituto Geográfico do Exército
Administrative limitsCaia Irrigation Scheme limitsOfficial Administrative Map of Portugal and Irrigation Scheme Maps of Portugal—Ministry of Environment
Nitrate Vulnerable Zone limitsNitrates Vulnerable Zones limits Map (1:25,000)Ministry of Agriculture
CropsCropping patterns, crop phenological stagesSentinel-2 images (Level-2A), Google Earth Satellite Images (QMS)
Crop coefficients[7,8,51]
Land cover classesLand use map (COS)
Land cover classes, Corine Land Cover (CLC) 2018 (spatial resolution 20 m)
Crop plots map
Ministry of Environment
Copernicus Program (EEA, JRC)Caia Water Users Association
Farm/plot dataFarm/plot identification (WMS/WFS format)Ministry of Agriculture
GroundwaterWater table depthWater Resources National Information System—Ministry of Environment
IrrigationCrop systems, irrigation calendars, irrigation systems, and soil moisture dataCaia Water Users Association
Irrigation systemsIrrigation systems efficiency[52]
Table 2. Farmers’ field locations and characteristics used to feed and to validate the water balance model simulations, Caia Irrigation Scheme.
Table 2. Farmers’ field locations and characteristics used to feed and to validate the water balance model simulations, Caia Irrigation Scheme.
FieldCropArea (ha)Irrigation SystemSoil TextureType of Field DataYears with DataData Provider Water 14 02324 i001
AMaize FAO60058.8PivotLoamy sandIrrigation amounts and frequency
Soil moisture
2018, 2019, 2020APAP
BMaize FAO20037.4PivotLoamIrrigation amounts and frequency2017ABCaia
COlive grove_126.8DripClay loamIrrigation amounts and frequency
Soil moisture
2018, 2019, 2020APAP
DOlive grove_2396.3DripClay loamIrrigation amounts and frequency2017ABCaia
EProcessing tomato_120.2DripSilty loamIrrigation amounts and frequency
Soil moisture
2019, 2020APAP
FProcessing tomato_228.7DripLoamIrrigation amounts and frequency2017ABCaia
Note: APAP—Portuguese association of precision agriculture; ABCaia—Caia water users association.
Table 3. Example of the average lengths of the spring crop growth stages and of the seeding/planting dates estimated from NDVI temporal profiles for 2017.
Table 3. Example of the average lengths of the spring crop growth stages and of the seeding/planting dates estimated from NDVI temporal profiles for 2017.
CropLengths of Crop Growth Stages (Days)Seeding/Planting DateObs.
InitialDevelopmentMidLate
Paddy rice3027501325/Apr.
Tomato3037501325/Apr.
Melon3527332730/Apr.
Bell pepper3060301302/May
Rapeseed3030402502/Apr.
Sunflower2532403020/May
Grain maize2027533725/Apr.
Silage maize2230351213/May
Fodder crop multi-cuts3011720---16/Nov.1st cut
63120---2nd cut
31613---3rd cut
Table 4. Average soil properties for each class of land suitability for irrigation.
Table 4. Average soil properties for each class of land suitability for irrigation.
Land Suitability for IrrigationθFC
(g·g−1)
θWP
(g·g−1)
Z
(cm)
Bd
(g·cm−3)
I34.7217.021201.22
II30.8413.46801.23
III29.1212.52651.25
exc30.8314.18501.25
Note: θFC is the soil moisture content at field capacity, θWP is the soil moisture content at the permanent wilting point, Z is the soil depth, and Bd is the bulk density.
Table 5. “Goodness-of-fit” indicators obtained by comparing measured and simulated available soil water for the diverse crops during the 2018–2020 season parameters.
Table 5. “Goodness-of-fit” indicators obtained by comparing measured and simulated available soil water for the diverse crops during the 2018–2020 season parameters.
CropYearObservationsb0R2PBIAS
(%)
RMSE
(mm)
NRMSE
(%)
EF
Maize2018571.000.77−0.13.12.50.76
2019711.000.820.66.85.90.80
2020781.010.841.25.85.00.83
Olive grove20183281.030.854.018.67.30.80
20193641.010.931.113.65.60.92
20203481.000.80−0.217.56.20.77
Tomato for20191210.990.90−1.25.53.40.88
processing2020990.990.87−0.76.13.90.86
Note: b0—coefficient of regression forced to the origin; R2—coefficient of determination; PBIAS—percent bias; RMSE—root mean square error; NRMSE—normalized root mean square error; EF—modeling efficiency.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ferreira, A.; Rolim, J.; Paredes, P.; Cameira, M.d.R. Assessing Spatio-Temporal Dynamics of Deep Percolation Using Crop Evapotranspiration Derived from Earth Observations through Google Earth Engine. Water 2022, 14, 2324. https://doi.org/10.3390/w14152324

AMA Style

Ferreira A, Rolim J, Paredes P, Cameira MdR. Assessing Spatio-Temporal Dynamics of Deep Percolation Using Crop Evapotranspiration Derived from Earth Observations through Google Earth Engine. Water. 2022; 14(15):2324. https://doi.org/10.3390/w14152324

Chicago/Turabian Style

Ferreira, Antónia, João Rolim, Paula Paredes, and Maria do Rosário Cameira. 2022. "Assessing Spatio-Temporal Dynamics of Deep Percolation Using Crop Evapotranspiration Derived from Earth Observations through Google Earth Engine" Water 14, no. 15: 2324. https://doi.org/10.3390/w14152324

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop