A framework for irrigation performance assessment using WaPOR data: the case of a sugarcane estate in Mozambique
- 1Department of Land and Water Management, IHE Delft Institute for Water Education, 2611 AX Delft, the Netherlands
- 2Water Management Group, Delft University of Technology, 2600 AA Delft, the Netherlands
- 3Water Resources Management Group, Wageningen University & Research, 6700 AA Wageningen, the Netherlands
- 4Xinavane estate, Xinavane, Mozambique
- 5Inclusive Agri-Food Value Chain Development Programme (PROCAVA) – FDA, Maputo, Mozambique
- 6VanderSat B. V., Agri, Food and Commodity unit, Wilhelminastraat 43a, 2011 VK Haarlem, the Netherlands
Correspondence: Abebe D. Chukalla (email@example.com)
The growing competition for finite land and water resources and the need to feed an ever-growing population require new techniques to monitor the performance of irrigation schemes and improve land and water productivity. Datasets from FAO's portal to monitor Water Productivity through Open access Remotely sensed derived data (WaPOR) are increasingly applied as a cost-effective means to support irrigation performance assessment and identify possible pathways for improvement. This study presents a framework that applies WaPOR data to assess irrigation performance indicators, including uniformity, equity, adequacy, and land and water productivity differentiated by irrigation method (furrow, sprinkler, and centre pivot) at the Xinavane sugarcane estate, Mozambique. The WaPOR data on water, land, and climate are in near-real time and spatially distributed, with the finest spatial resolution in the area of 100 m. The WaPOR data were first validated agronomically by examining the biomass response to water, and then the data were used to systematically analyse seasonal indicators for the period 2015 to 2018 on ∼8000 ha. The WaPOR-based yield estimates were found to be comparable to the estate-measured yields with ±20 % difference, a root mean square error of 19±2.5 t ha−1 and a mean absolute error of 15±1.6 t ha−1. A climate normalization factor that enables the spatial and temporal comparison of performance indicators are applied. The assessment highlights that in Xinavane no single irrigation method performs the best across all performance indicators. Centre pivot compared to sprinkler and furrow irrigation shows higher adequacy, equity, and land productivity but lower water productivity. The three irrigation methods have excellent uniformity (∼94 %) in the four seasons and acceptable adequacy for most periods of the season except in 2016, when a drought was observed. While this study is done for sugarcane in one irrigation scheme, the approach can be broadened to compare other crops across fields or irrigation schemes across Africa with diverse management units in the different agroclimatic zones within FAO WaPOR coverage. We conclude that the framework is useful for assessing irrigation performance using the WaPOR dataset.
Increasing agricultural production to feed the growing global population can be achieved through either expanding agricultural land or by increasing productivity of the existing agricultural areas. With growing competition and scarcity of finite water and land resources, as well as the environmental and social costs of expanding agricultural land (Hess et al., 2016), improving irrigation performance indicators including land and water productivity has a clear preference.
The increasing global demand for sugar is also reflected in the steady increase in sugarcane production in Mozambique at an average annual rate of 10 % (FAO, 2019). The majority of this increase comes from expanding agricultural land (Hess et al., 2016). Whilst Moraes et al. (2018) estimate there is a vast potential for expanding sugarcane production in Mozambique (∼15 % of the land area is suitable for sugarcane production), water and land resources in the country are under increasing strain due to land degradation (Sutton et al., 2016), sectoral competition, and climate effects (e.g. drought and flood) (Van der Zaag and Carmo Vaz, 2003; Arndt et al., 2011). With the land productivity well below the global average (Binswanger-Mkhize and Savastano, 2017; Nkamleu, 2011), and amongst the lowest in the southern African region (Johnson et al., 2014), there is an opportunity to meet the demand without expanding the agricultural land. Thus, raising sugarcane productivity per unit of land and water on existing croplands needs to be explored by conducting irrigation performance assessment.
Monitoring irrigation performance indicators is key in checking general health, comparing the spatial and temporal performances of the scheme, and looking for causes and providing corrective action that aims at improving overall service provision and productivity (Molden et al., 1998; Bos et al., 2005). Traditional irrigation performance assessment considers indicators that can be categorized as (i) water balance, water service, and maintenance, (ii) environment, and (iii) economic indicators. The water balance, water service, and maintenance indicators are water-flux- and production-based indicators. The water-delivery- and production-based indicators include uniformity (evenness of water distribution within fields), equity (uniformity of water distribution between fields), adequacy (sufficiency of irrigation delivery compared to the requirement), land productivity (production per unit area), water productivity (production per unit water use), and efficiency (the fraction of productive water use) (Molden and Gates, 1990; Bos, 1997; Molden et al., 1998). These irrigation performance indicators were assessed using field data such as flow (discharge), crop yield, and plot level water consumption estimates using lysimeters or crop models (Araya et al., 2011; Dejen, 2015; Edreira et al., 2018).
Recent developments and improvements of remote-sensing (RS) products offer a viable alternative (Bastiaanssen et al., 1996; Karimi et al., 2011). RS-derived data have been increasingly applied as a cost-effective means for irrigation performance assessment. RS-derived irrigation performance assessment is based on production and actual water consumption, the latter of which is fairly considered the net outcome and result of effective rainfall and irrigation, allowing for a hydrological assessment and quantification of the net water abstracted by irrigated crops. In addition, it provides spatially distributed data, covers long periods and wide areas, and can be done retrospectively (Bastiaanssen et al., 1996; Karimi et al., 2011). Field data, in contrast, do not represent the spatial variation across an irrigation system well and are costly to obtain (Bastiaanssen et al., 2000). Traditional and RS-based performance assessments are complementary as the former has strength in observing the horizontal water fluxes such as discharges, while the latter has strength in observing high-resolution vertical water fluxes and biomass production.
Earlier studies provide insight into the application of RS-derived data to assess irrigation performance indicators. In this research, the earlier RS-based irrigation performance assessment studies are strengthened by considering a simple consistency check to validate the RS-derived data for established biomass response to water consumption (Steduto and Albrizio, 2005) and by introducing a comprehensive framework that guides the step-by-step translation of RS-derived datasets into irrigated agricultural performance indicators. In addition, the current study introduces a climate normalization factor that enables the spatial and seasonal comparison of irrigation performance indicators. The climate normalization factor is applied to distinguish climatic factors from agricultural management factors in their effect on irrigation performance.
This study first evaluates the FAO's portal to monitor Water Productivity through Open access Remotely sensed derived data (WaPOR) for consistency based on the established agronomic principle (biomass response to water consumption). It is then used to develop a framework to assess irrigation performance indicators, including adequacy, uniformity, equity, and land and water productivity. This framework is then used to assess the irrigation performance at Xinavane sugarcane estate differentiated by irrigation method.
2.1 Study area
The study focusses on one of the largest sugarcane estates in Maputo province in Mozambique, the Xinavane estate. The estate is located on the banks of the Incomati River, approximately 136 km northwest of Maputo. This region is characterized by optimal conditions for sugarcane production in terms of climate, soils, and water availability. With a seasonal long-term average precipitation of 721 mm yr−1 (den Besten et al., 2020), sugarcane production requires irrigation water, especially during the dry season, supplied by the Incomati River.
The most important water infrastructure in the Incomati Basin in Mozambique is the Corumana Dam, which was built for improving flood control and regulating downstream irrigation abstractions (including Xinavane) and hydropower production (de Boer and Droogers, 2016). Xinavane sugarcane estate, despite receiving allocations from the dam, remains largely vulnerable to climate variability. During a recent drought in 2016, reservoir levels in the Corumana Dam dropped drastically, and little water was available for irrigation in the Xinavane sugarcane estate. This resulted in a significant reduction in sugarcane production in 2016 compared to previous years (Tongaat Hullet, 2018). Such events are expected to continue to occur. To partially address this, Mozambique put drought mitigation measures in place for the Xinavane area, including the construction of the new Moamba Major Dam (760×106 m3) and the heightening of the Corumana Dam wall, which will result in a capacity increase from 879×106 m3 to 1260×106 m3 (Tongaat Hullet, 2018).
The widely used irrigation methods at the Xinavane sugarcane estate are furrow, overhead sprinkler (hereinafter referred to as sprinkler), and centre pivot irrigation (Fig. 1). A total of 8027 ha categorized into 387 georeferenced fields and three irrigation application methods are considered in our analysis. Furrow, sprinkler, and centre pivot irrigation covers 3343, 3629, and 1055 ha, respectively. The average field size under furrow, sprinkler, and centre pivot irrigation methods is 17, 18.3, and 55.8 ha, respectively. All fields in the sample are operated and managed by the estate; fields operated by outgrowers were excluded from the analyses.
2.2 WaPOR datasets
Datasets from FAO's portal to monitor Water Productivity through Open access Remotely sensed derived data (WaPOR; URL: https://wapor.apps.fao.org/home/WAPOR_2/1, last access: 15 December 2021) are used for the analyses as it provides the required layers to estimate both land and water productivity. The database covers Africa and the Near East regions in near-real time for the period from 2009 to date (2021) (FAO, 2020c). WaPOR datasets are available at the continental scale (Level 1 at 250 m), country scale (Level 2 at 100 m), and project level (Level 3 at 30 m). The latest WaPOR version (WaPOR v2.1) is an improvement from WaPOR v1.0, following the quality assessments by IHE Delft and ITC (Mul and Bastiaanssen, 2019; FAO, 2020a). The methodology used for compiling the actual evapotranspiration of WaPOR is based on the ETLook method (Bastiaanssen et al., 2012) and further developed by the FRAME consortium (the full description of the methodology is provided in FAO, 2020b). WaPOR v2.1 was found suitable for inter-plot comparison of irrigation performance indicators for plots larger than 2 ha (Blatchford et al., 2020).
At Xinavane, the finest resolution of the WaPOR data is 100 m (Level 2). The WaPOR Level 2 datasets used in this study include layers for actual evaporation (E), transpiration (T), and net primary production (NPP) at a dekadal (10 d) timescale. In addition, dekadal precipitation at 5 km resolution, dekadal reference evapotranspiration at 20 km resolution, and annual land cover classification (LCC) at 100 m resolution were used. The precipitation (P) and reference evapotranspiration (RET) datasets were resampled to 100 m resolution using the nearest-neighbour resampling techniques (GDAL, 2021). An overview of the WaPOR data used in the analyses is presented in Table 1.
Although there is a continuous WaPOR L2 dataset (100 m) available from 2009 to date (2021), only the data from 2014 are derived from the PROBA-V satellite. The data prior to 2014 are derived from resampled Level 1 (250 m) data, which are obtained from the MODIS satellite. Since this creates a discontinuity in the data as observed by Chukalla et al. (2020b), the pre-2014 data have been discarded in this analysis, and only data starting from the 2014–2015 growing season onwards have been selected.
2.3 A framework for assessing irrigation performance using WaPOR data
Figure 2 shows the flowchart describing the approach to assessing WaPOR-based irrigation performance indicators at the Xinavane sugarcane estate. Irrigation performance indicators are derived from WaPOR and field data in three main steps. First, actual evapotranspiration (ET), reference evapotranspiration (RET), and net primary production (NPP) layers of FAO WaPOR are preprocessed to match the spatial resolution, remove non-crop pixels using crop map or land cover classification (LCC), and undergo a quality check. Second, the seasonal ETa (ETa,s), seasonal potential evapotranspiration (ETp,s), and seasonal NPP (NPPs) are calculated from their respective WaPOR layers between the start of the season (SOS) and end of the season (EOS) for each plot. ETp,s is derived from RET and crop coefficient (Kc). Finally, the irrigation performance indicators are analysed. At this stage, NPPs is translated to above-ground biomass (hereafter referred to as biomass – B) using crop-specific information – above over total biomass (AOT) for non-root corps or below over total for root and tuber crops, light use efficiency correction factor (fc), and moisture content of fresh biomass (mc). The biomass is multiplied by harvest index (HI) to derive the crop yield. The remainder of this section describes the input data and equations used in each step in more detail.
2.3.1 Seasonal water consumption and crop yield
The sugarcane estate operates on a ratooning system. Thus, the start of the growing season (1 d after harvesting) and end of season (harvesting date) vary per field. The actual growing period of each field was used to calculate the production per unit of land and per unit of water consumed. The average length of the growing season is 347±32 d. This study covers four growing seasons: season 1 (2014/15), season 2 (2015/16), season 3 (2016/17), and season 4 (2017/18), reported as 2015, 2016, 2017, and 2018, respectively, i.e. the year the fields are harvested (Fig. 3).
Seasonal water consumption
Actual water consumption refers to the amount of water that is depleted from the root zone through the process of transpiration by a crop and direct evaporation from the soil represented by WaPOR E+T(ETa). The seasonal ETa is the total actual water consumption during the cropping season.
The seasonal NPP layer from WaPOR, accumulated over the crop-growing period (Fig. 3), is converted to above-ground biomass (B) in kilograms per hectare (kg ha−1) and crop yield (Y) in kilograms per hectare using Eqs. (1) and (2) (Mul and Bastiaanssen, 2019):
where mc [–] is the moisture content of the fresh biomass, fc [–] is the light use efficiency (LUE) correction factor calculated by dividing the LUE of the crop (in this case sugarcane) by the LUE of a generic crop type that WaPOR NPP layer uses (2.7 MJ g−1 biomass; FAO, 2018, 2020b), and AOT [–] is the ratio of above-ground over total biomass. B and Y can be expressed in tonnes per hectare (t ha−1), by dividing the amount in kilograms per hectare by 1000. Crop yield is calculated by multiplying the biomass by the harvest index (HI [–]):
In absence of field data, literature was consulted to estimate these crop parameters. Table 2 presents the values and the source of the parameters.
The WaPOR-based sugar cane yield was validated with sugarcane yields as measured by the Xinavane estate for four seasons on 387 fields. In addition, the WaPOR-based biomass and water consumption were checked for consistency with agronomic principles. An increasingly strong linear relationship is expected between biomass and evapotranspiration (Steduto and Albrizio, 2005), between biomass and transpiration (De Wit, 1958), and between biomass and normalized transpiration (Steduto and Albrizio, 2005), whereby the normalized transpiration is the sum of the daily ratio of transpiration over reference evapotranspiration over the crop season (Steduto et al., 2007).
2.3.2 Performance assessment indicators
The irrigation performance indicators selected for this study are uniformity, equity, adequacy, and productivity; these were selected as these could be assessed (sometimes with a slight modification) using the WaPOR data. These performance indicators are further explained below, and the set of equations for water-consumption- and production-based performance indicators are presented in Table A1.
Uniformity measures the evenness of water consumption within an irrigated field. It is calculated by assessing the coefficients of variation (CV) of seasonal ETa within a field. Thus, uniformity is 1 minus the CV (Ascough and Kiker, 2002). It serves as a measure for the heterogeneity of soil water storage capacity and thus water storage efficiency in a field. It can serve as a proxy for irrigation distribution uniformity (Burt et al., 1997) in farms where the management is central and the same level of inputs is consistently applied (e.g. variable rate input application is not practised). Other factors like soil type, fertility, pest, and crop variety can also affect actual water consumption and thus uniformity. Thus, the CV of seasonal ETa indicates the combined effect of all factors (water, fertility, pests, diseases, and salinity).
According to Pitts et al. (1996), the acceptable standard uniformity of irrigation application distribution for centre pivot, sprinkler, drip, and furrow irrigation methods is 75 %, 75 %, 85 %, and 65 %, respectively. The distribution uniformity exceeding the standard threshold is considered excellent.
Equity measures the evenness of water consumption between fields within an irrigation scheme with a homogenous crop, which could be a proxy for an even distribution of water to the different irrigated fields. It is calculated as the CV of the average ET of each field, which is an indication of equity in the scheme. A CV of 0 % to 10 % is defined as good equity, CV of 10 % to 25 % as fair equity, and CV > 25 % as poor equity (Bastiaanssen et al., 1996; Karimi et al., 2019).
Adequacy (A) is the measure of the degree of agreement between the actual water use and crop water requirement (Bastiaanssen and Bos, 1999; Clemmens and Molden, 2007). Adequacy is estimated as the ratio of seasonal actual evapotranspiration (ETa,s) over seasonal potential evapotranspiration (ETp,s) (Kharrou et al., 2013; Karimi et al., 2019). Potential evapotranspiration is the maximum crop evapotranspiration from disease-free and well-fertilized cropped fields under optimum soil water conditions; it is calculated by multiplying reference evapotranspiration by the crop coefficient (Allen et al., 1998). The seasonal ETp,s of sugarcane is aggregated from the monthly value of crop coefficient multiplied by the reference evapotranspiration (Table A2). Good adequacy performance is defined for the range of , acceptable range , and poor performance A≤0.68 (Karimi et al., 2019).
Productivity is a measure of benefit generated per unit of resource used. The benefit could be biophysical, economic, and/or social; the resource base could be consumed or supplied water or land covered by the crop (Zwart and Bastiaanssen, 2004; Hellegers et al., 2009; Karimi et al., 2011). This study focussed on biophysical production per unit of land or evapotranspiration, also known as land and water productivity.
Land productivity is defined as biomass production or crop yield per unit of land. For water, we similarly distinguish biomass water productivity (WPb) and crop yield water productivity (WP). WPb is defined as the ratio of biomass over seasonal ETa,s, whereas WP is defined as the yield over ETa,s. Since for sugarcane we use a harvest index of 1, WPb is equal to WP here.
Spatial–temporal variations can be caused by both management practices and climate. Figure B1 shows a correlation between water productivity and reference evapotranspiration (r2 of 0.5, 0.7, and 0.8 for furrow-, sprinkler-, and centre-pivot-irrigated fields, respectively). The correlation between actual evapotranspiration and reference evapotranspiration (Fig. B2) is even stronger (r2>0.8). Thus, to exclude the climate-related factor, we normalized the water productivity and evapotranspiration using a climate normalization factor. This is defined as the ratio of the weighted average reference evapotranspiration (weighted based on the field size and growing length of the fields) to the reference evapotranspiration at the field (Eq. 3).
where fnorm [–] is the normalizing factor for the selected indicator, is weighted average reference evapotranspiration, and RETi is reference evapotranspiration at a field in millimetres per season.
2.4 Consistency check of WaPOR data
Figure 4 shows the relationship between biomass (B; WaPOR-derived and WaPOR-observed) and water consumption of irrigated fields categorized by irrigation method for the year 2018 (with the Supplement, Fig. S1, showing the other 3 years from 2015 to 2017). In furrow- and sprinkler-irrigated fields, the WaPOR-derived biomass and actual evapotranspiration show a high correlation (a minimum r2 of ∼0.83 (n≈150) in 2015, 2017, and 2018 and r2≈0.63 in the relatively dry year of 2016), indicating consistency between the two independently generated datasets. For the centre-pivot-irrigated fields, r2 is much lower, with a value of ≈0.6 in 2015, 2016, and 2017 and the lowest r2 of 0.2 (n≈19) in 2018. The low number of fields irrigated by centre pivots may have contributed to the low correlation. Moreover, the estate-observed yield at Xinavane sugar estate versus ETa shows a high spread and thus a low correlation (r2≈0.13).
Table S1 in the Supplement provide the analyses of the relationship between biomass and transpiration and biomass and normalized transpiration for the entire period of analyses (2015–2018). In contrast to expectations based on agronomic principles, the correlation decreases when considering biomass and transpiration (∼0.80) and biomass and normalized transpiration () (∼0.71) (see the Supplement). The accuracy of the evaporation and transpiration split in WaPOR is therefore questioned; this was also observed by Mul and Bastiaanssen (2019). Further analyses will therefore only focus on indicators that use evapotranspiration, not evaporation and transpiration, as input. For instance, the beneficial fraction (i.e. the ratio of transpiration over evapotranspiration) is not included in the analysis. Yet, two tests based on WaPOR-derived biomass and total actual evapotranspiration (ETa) have confirmed the agronomic expectations (Table S2). The first is that the correlation coefficient of the linear regression line passing through the origin for the biomass vs. normalized actual water consumption is higher than that of the correlation coefficient for the biomass vs. actual water consumption. Second, the crop water productivity normalized by reference evapotranspiration (WP*) is confirmed to be conservative and within the range of values for C4 crops (30–35 g m2), including sugarcane (Steduto et al., 2007, 2009).
3.1 Seasonal water consumption
Figure 5 shows the seasonal actual and potential evapotranspiration and seasonal precipitation at Xinavane sugarcane estate, distinguished by the three irrigation application methods. The four-season (2015 to 2018) average precipitation is 640 mm per season and ranges from the minimum of 500 mm per season in 2016 to the maximum precipitation of 875 mm per season in 2017. The four-season average ETa at Xinavane is 1350 mm per season, and its average seasonal values range between 1255 mm per season in 2018 at furrow-irrigated fields to 1533 mm per season in 2016 at fields irrigated using centre pivots. In the four seasons the ETa is significantly the highest (P value < 0.05) at fields irrigated using centre pivots, followed by sprinklers and furrows (Table A4).
The high average ETa over Xinavane irrigation scheme in 2016 coincides with the reported drought year. This mainly manifested itself with high ETpot as the annual precipitation that falls within the command area was not much lower than in 2015 and 2018. After normalization for climate variation, the normalized ETa is actually lowest for 2016, indicating higher water deficit (lowest actual per unit of potential evapotranspiration), with the drought having more impact on sprinkler and furrow irrigation than on centre pivot irrigation. Despite the ETa being the highest in 2016, when normalized by climate, the results show that 2016 experiences the highest water deficit. The four-season average actual water consumption of centre pivots remains the highest followed by sprinklers and furrows, except for 2016, when the sprinkler-normalized ETa is at the same level as furrow ETa (Fig. 6). This indicates that the sprinkler system was more affected by the drought conditions in 2016 compared to the other systems.
3.2 Performance of irrigation delivery
The uniformity of water consumption within the fields is ∼94 % for all three irrigation methods (Fig. 7). The calculated uniformity is above the standard values per irrigation method and is therefore considered excellent. Centre pivots show an even higher uniformity than the other irrigation methods.
The average seasonal coefficient of variation (CV) of ETa,s among fields irrigated by the same irrigation method is 15 % (Fig. 8). Fields irrigated using furrows, with a CV of 18 %, have the highest heterogeneity in water consumption compared to areas irrigated using the sprinkler (CV = 14 %) and centre pivot irrigation method (CV = 13 %). The coefficient of variation of water consumption between fields irrigated by a particular irrigation method and thus equity of water use among the fields is considered fair.
The four-season average adequacy varies spatially across the Xinavane irrigation scheme, with visible differences between fields irrigated using centre pivots compared to fields irrigated using furrows and sprinklers for the period analysed. Figure 9 shows the highest adequacy for fields irrigated using centre pivots (0.75) followed by fields irrigated using sprinklers and furrows (∼0.69). In the study period, the adequacy performance at fields under centre pivots falls in the acceptable range (from 0.68 and 0.8) for sugarcane (Karimi et al., 2019). The adequacy in fields under sprinkler and furrow irrigation also is acceptable except in the year 2016, which is recognized as a drought year, when adequacy was poor.
3.2.4 Land productivity
The 4-year seasonal average WaPOR-based yield is 89 t ha−1 (86 t ha−1 for fields irrigated using furrows, 88 t ha−1 for areas irrigated using sprinklers, and 93 t ha−1 for fields irrigated using centre pivots). For all years (except 2017), the highest sugarcane yield (land productivity) at Xinavane is found in fields irrigated by centre pivots, followed by fields irrigated by sprinkler and furrow irrigation methods (Fig. 10).
The 4-year seasonal WaPOR yield is in the same order of magnitude compared to the estate-measured sugarcane yield: 86 t ha−1 vs. 81.4 t ha−1, 88 t ha−1 vs. 93 t ha−1, and 93 t ha−1 vs. 99 t ha−1 for fields irrigated using furrow, sprinkler, and centre pivot irrigation methods, respectively. Part of the minor discrepancy between the WaPOR- and estate-measured yield could be due to the selection of crop parameters such as harvest index and moisture content. Yet, the comparison between both yields shows acceptable statistics (Table A3 in the Appendix), with a root mean square error of 19±2.5 t ha−1 and mean absolute error of 15±1.6 t ha−1.
Whilst the average values for WaPOR-based yields are of the same magnitude as the estate-observed data (65 % of yield differences at the fields are within ±20 %), WaPOR overestimates relatively low yields (marks on scatter plot above 1:1 line) and underestimates relatively high yields (marks on scatter plot below 1:1 line) (Fig. 11). WaPOR yields thus show a marked less variation in yields than reported by the estate.
3.2.5 Water productivity
The seasonal and four-season average water productivity at Xinavane is shown in Fig. 12. The four-season average water productivity is the highest for furrow-irrigated fields (6.9 kg m−3), compared to the values for fields irrigated with sprinklers (6.7 kg m−3) and centre pivots (6.6 kg m−3). One of the reasons for such differences is the fraction of ETa being utilized for productive purposes (transpiration) compared to non-productive evaporation. Raes et al. (2013) report that centre pivot and sprinkler irrigation wets 100 % of the field compared to furrow irrigation that wets ∼80 % of the field and thus results in higher evaporation rates, which is in line with our observations.
The large variation of WP over the years (Fig. 12) is also apparent after normalization for climate variation (Fig. 13). The normalized WP is highest in a relatively dry year (2016) compared to the other 3 years; this is opposite to WP, where 2016 has the lowest WP. It indicates that climate-related parameters expressed through potential evapotranspiration have a large impact on the WP. The normalized WP shows the variations which are related to management practices; during the drought of 2016, the Xinavane estate practised deficit irrigation, which is reflected in the high normalized WP values.
4.1 The framework
The presented framework was used to conduct an irrigation performance assessment using WaPOR data. Our analysis shows that fields irrigated using centre pivots have the highest equity, adequacy, and land productivity, followed by fields irrigated using sprinklers and furrows. This outcome agrees with the conclusion by Karimi et al. (2019), who assessed performance of irrigated sugarcane in Eswatini (Swaziland) by differentiating areas according to management regimes including irrigation methods. The adequacy performance under the three irrigation methods was generally acceptable, except in 2016, when performance of all three irrigation methods was poor. Fields under centre pivots do, however, have the lowest water productivity, followed by sprinkler and furrow irrigation, which is contrary to the finding by Karimi et al. (2019), who reported the WP of centre pivots to exceed that of furrow irrigation. In fact, it is claimed that pressurized irrigation (sprinklers and centre pivots) improves uniform distribution and application efficiency of irrigation water and increases crop yield (Playán and Mateos, 2006). Yet, these irrigation methods increase seasonal evaporation (Playán and Mateos, 2006), which could be due to differences in percentage of land wetted. Our findings show that the uniformity of water consumption on the fields under the three irrigation methods is reasonably comparable and high (∼94 %), which can be regarded as excellent according to the standard set by Pitts et al. (1996). The high uniformity of water consumption in furrow-irrigated fields is in the same range as that of centre pivots and sprinklers, which is unlike what was found in South Africa (Griffiths and Lecler, 2001).
The results of normalization for climate differences of the water consumption and water productivity allows for comparison of the results under different climate conditions (different years). While the ranking for the different irrigation technologies according to the indicators remains the same, it clearly shows the impact of the climate. In particular, during the drought year of 2016 when the potential evapotranspiration was relatively high, the normalized water consumption was low, indicating higher water deficit compared to the other years. The impact on sprinkler-irrigated fields was the highest. On the other hand, the normalized WP during 2016 was the highest of all the years, even though the WP was lowest for the same biomass in 2016, indicating the climate having a large impact on non-beneficial evaporation.
This study shows that the presented framework offers a systematic approach to assess irrigation performance indicators using WaPOR and field data. Five WaPOR-derived irrigation performance indicators, namely uniformity, equity, adequacy, and land and water productivity, are used to monitor the quality of the irrigation and agronomic services. Our framework builds on earlier studies that assess irrigation performance indicators based on RS (Karimi et al., 2019; Blatchford et al., 2020) and provides a comprehensive and simple step-by-step framework to conduct an agronomic validation of RS-derived products using WaPOR data. The approaches in the framework are scripted with Python in Jupyter Notebooks, which can be run on local machines and Google Colaboratory (Colab). The Python script is published together with observed yield data in GitHub (Chukalla et al., 2020a). It shows that with limited field information (crop type and cropping season) and some parameters obtained from the literature, the analyses can be implemented.
4.1.1 Limitations of the WaPOR database
The linear relationship between the independently derived WaPOR biomass and water consumption agrees with the expected agronomic principles (De Wit, 1958; Steduto and Albrizio, 2005). However, the correlation coefficient of the biomass versus actual evapotranspiration is higher than the correlation coefficient of the biomass versus transpiration and biomass versus normalized transpiration. This implies an inaccurate estimation of transpiration (T) and evaporation (E) in WaPOR. WaPOR separates the available energy into T and E using a factor α⋅LAI, where α is the light extinction factor (FAO, 2018; Mul and Bastiaanssen, 2019). A review on values for α shows large differences between different land use classes and within land use classes (Zhang et al., 2016). Thus, WaPOR applying only one fixed value for α could have serious implications for the use of the T and E layers of WaPOR, such as in quantifying beneficial fraction (the ratio of transpiration over evapotranspiration).
Even though the analyses seem to be consistent with the understanding of how the different irrigation technologies perform, there are some known limitations of RS and WaPOR data in particular, which need to be mentioned here. These may stem from (i) the Land Surface Temperature (LST) data used by WaPOR (which are taken from MODIS and have a resolution of 1 km, used to derive moisture stress and thus to calculate the actual evapotranspiration and net primary production; this could be the cause for the reduced variation of WaPOR biomass data and may affect the spatial variation of evapotranspiration as well); (ii) land cover noise of non-sugarcane land use such as farm roads and irrigation and drainage infrastructures within a pixel; (iii) the number of cloud-free RS images on which the analysis and numerical interpolation are based (the fewer the cloud-free images, the poorer the data quality and the higher the uncertainty in the indicators one can expect); (iv) the time of day when the images are taken (determinant for which part of the daily ET curve is monitored and the time of day the water stress is more or less severe); and (v) the angle of image capture and its correction function.
The methods used in WaPOR for data production and statistical methods for the reconstruction of missing values are, however, on a par with those used in other RS-based products for monitoring agro-hydrological parameters developed by the scientific community. As such, some of these limitations are inherent to the use of remote sensing in general. Yet, our analysis shows consistency between the different datasets.
4.1.2 Limitation of the crop-related information
Crop-specific parameters such as harvest index, the moisture content of the fresh yield, and the ratio of above-ground over total biomass were fixed values and determined using literature and fieldwork in Ethiopia. However, it is known that these crop parameters can vary significantly based on climatic or field management conditions. Other variations may stem from differential exposure to pests and diseases and soil and rooting conditions caused by waterlogging (den Besten et al., 2021) and soil salinity, which are not catered for. We were unable to determine how much these assumptions affect the results. All these factors are potential sources of (slight) deviations in the numerical output of WaPOR that may lead to over- and underestimations of crop yield and WP.
Having noted this, we did perform a validation of the WaPOR biomass data using observed harvested cane data of more than 300 fields over four seasons. WaPOR biomass data for ∼65 % of the field level comparison differed within a ±20 % range. The comparison between the estate-measured yield and WaPOR biomass showed acceptable statistics (Table A3).
4.2 The way forward
Investments in high-quality public domain global and regional remote-sensing data products for water and lands uses, e.g. WaPOR datasets, have made it possible to conduct spatio-temporal analysis of irrigation performance at multiple scales from an irrigation scheme to district scale, basin scale, and the whole country. This provides a great advantage, especially in areas where both water and land resources are scarce and in situ data are scant. This study presents a RS-based assessment framework and showcases the power of using the WaPOR dataset in providing spatial and temporal irrigation performance indicators. Such information cannot be generated with the data collected traditionally (point data) or would come at a significant cost.
Yet, accurate interpretation of the results, diagnosis of the causes of the performance variation, and formulation of practical solutions cannot be done unless the WaPOR analyses and results are complemented with observed data of field conditions (e.g. the level of water and nutrient inputs, waterlogging, and salinity levels) that can help explore the constraints. Though this limitation puts a disclaimer on our findings, the procedures in this study can provide a useful reference for similar future studies.
Subsequent studies could additionally consider socio-economic performance indicators, such as social water productivity (e.g. employment per unit water or land use) and economic water productivity (economic return per unit water or land use), which could help to implement comprehensive performance assessment of irrigation schemes.
Remote-sensing datasets are increasingly applied as an innovative tool for monitoring the performance of irrigation schemes in order to improve land and water productivity amid the growing competition for finite and even dwindling resources (land and water). In this study, first, the remotely sensed FAO WaPOR dataset was successfully validated by comparing WaPOR-derived sugarcane yield with field observations, as well as agronomically. The yield response to water confirms agronomic expectations: (i) the correlation between biomass and actual water consumption normalized for climate is stronger than the correlation between biomass and actual water consumption, and (ii) the water productivity of sugarcane normalized by reference evapotranspiration falls within the conservative values reported for C4 crops. Second, the WaPOR-derived datasets were applied to assess irrigation performance indicators, including uniformity, equity, adequacy, and land and water productivity at Xinavane sugarcane estate, segmented by irrigation method. We conclude that the systematic approach demonstrated in the current study can serve as a framework to operationalize the use of WaPOR-derived data and other increasingly available RS-derived products for irrigation performance monitoring and assessment.
The comprehensive WaPOR-based irrigation performance assessment in this sugarcane state finds that fields irrigated by centre pivots have the highest adequacy, land productivity, and equity, followed by sprinkler- and furrow-irrigated fields, but the lowest water productivity.
We identified that part of the spatial and seasonal variation of indicators, water productivity, and seasonal water consumption in particular is explained by non-climatic factors that can be influenced by management interventions. Investigating the root causes of the land productivity variation and whether proper management of inputs and controlling of salinity and drainage could improve productivity and the overall performance require further study, including field-based observations.
* where SOS and EOS are the start of season and end of season, ETa,s is seasonal actual evapotranspiration, ETp,s and ETp,m are seasonal and monthly potential evapotranspiration, RETm is monthly reference evapotranspiration, kc,m is the crop coefficient, and NPPs is seasonal net primary production.
* Sum is the product of count [–] and average [mm per season].
The Python code in the Jupyter Notebook and Google Colaboratory (Colab) environment that is used for geospatial analysis of the current study can be freely obtained from the GitHub repository of the Water Accounting and Water Productivity group at https://doi.org/10.5281/zenodo.4641360 (Chukalla et al., 2020a).
All the remote-sensing-derived WaPOR datasets including evaporation, transpiration, net primary production, precipitation, reference evapotranspiration, and land cover maps can be freely obtained from the WaPOR portal (https://wapor.apps.fao.org/home/WAPOR_2/1; FAO, 2020c). The boundary (shapefile) of the case study and the estate-observed yield data can be acquired from the GitHub repository of the Water Accounting and Water Productivity group at https://doi.org/10.5281/zenodo.4641360 (Chukalla et al., 2020a).
The supplement related to this article is available online at: https://doi.org/10.5194/hess-26-2759-2022-supplement.
ADC, MLM, and PK conceptualized the paper and designed the methodology; ADC, with input from MLM and PK, scripted the Python code; PvdZ, EM, EM, and NdB supported the curation and interpretation of the field data; ADC performed the formal analysis and visualization; ADC and MLM wrote the original draft of the paper; MLM, PvdZ, GvH, PK, and ADC reviewed and edited the paper; and ADC, GvH, MLM, and PvdZ responded to the reviewers' comments.
At least one of the (co-)authors is a member of the editorial board of Hydrology and Earth System Sciences. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors would like to express their gratitude to Tongaat Hulett for access to useful data on Xinavane estate farm.
This research has been supported by the Water Productivity Improvement in Practice (Water-PIP) project, which is funded by the IHE Delft Water and Development Partnership Programme (DUPC2) under the programmatic cooperation between the Directorate General for International Cooperation (DGIS) of the Ministry of Foreign Affairs of the Netherlands and IHE Delft (ID: DGIS Activity DME0121369).
This paper was edited by Dominic Mazvimavi and reviewed by two anonymous referees.
Allen, R. G., Pereira, L. S., Raes, D., and Smith, M.: Crop evapotranspiration – Guidelines for computing crop water requirements, FAO Irrigation and drainage paper 56, FAO, Rome, http://www.fao.org/docrep/x0490e/x0490e00.htm (last access: 15 April 2020), 1998.
Araya, A., Stroosnijder, L., Girmay, G., and Keesstra, S.: Crop coefficient, yield response to water stress and water productivity of teff (Eragrostis tef (Zucc.), Agr. Water Manage., 98, 775–783, 2011.
Arndt, C., Strzepeck, K., Tarp, F., Thurlow, J., Fant, C., and Wright, L.: Adapting to climate change: an integrated biophysical and economic assessment for Mozambique, Sustainabil. Sci., 6, 7–20, https://doi.org/10.1007/s11625-010-0118-9, 2011.
Ascough, G. and Kiker, G.: The effect of irrigation uniformity on irrigation water requirements, Water SA, 28, 235–242, 2002.
Bastiaanssen, W., Cheema, M., Immerzeel, W., Miltenburg, I., and Pelgrum, H.: Surface energy balance and actual evapotranspiration of the transboundary Indus Basin estimated from satellite measurements and the ETLook model, Water Resour. Res., 48, W11512, https://doi.org/10.1029/2011WR010482, 2012.
Bastiaanssen, W. G. and Bos, M.: Irrigation performance indicators based on remotely sensed data: a review of literature, Irrig. Drain. Syst., 13, 291–311, 1999.
Bastiaanssen, W. G., Van der Wal, T., and Visser, T.: Diagnosis of regional evaporation by remote sensing to support irrigation performance assessment, Irrig. Drain. Syst., 10, 1–23, 1996.
Bastiaanssen, W. G., Molden, D. J., and Makin, I. W.: Remote sensing for irrigated agriculture: examples from research and possible applications, Agr. Water Manage., 46, 137–155, 2000.
Binswanger-Mkhize, H. P. and Savastano, S.: Agricultural intensification: The status in six African countries, Food Policy, 67, 26–40, https://doi.org/10.1016/j.foodpol.2016.09.021, 2017.
Blatchford, M., M Mannaerts, C., Zeng, Y., Nouri, H., and Karimi, P.: Influence of Spatial Resolution on Remote Sensing-Based Irrigation Performance Assessment Using WaPOR Data, Remote Sens., 12, 2949, https://doi.org/10.3390/rs12182949, 2020.
Bos, M. G.: Performance indicators for irrigation and drainage, Irrig. Drain. Syst., 11, 119–137, 1997.
Bos, M. G., Burton, M. A., and Molden, D. J.: Irrigation and drainage performance assessment: practical guidelines, CABI publishing, https://doi.org/10.1079/9780851999678.0000, 2005.
Burt, C. M., Clemmens, A. J., Strelkoff, T. S., Solomon, K. H., Bliesner, R. D., Hardy, L. A., Howell, T. A., and Eisenhauer, D. E.: Irrigation performance measures: efficiency and uniformity, J. Irrig. Drain. Eng., 123, 423–442, 1997.
Chukalla, A. D., Mul, M., Tran, B., and Karimi, P.: Standardized protocol for land and water productivity analyses using WaPOR (v1.1), Zenodo [code and data set], https://doi.org/10.5281/zenodo.4641360, 2020a.
Chukalla, A. D., Mul, M., van Halsema, G., van der Zaag, P., Uyttendaele, T., and Karimi, P.: Water Productivity Analyses using the WaPOR Database. A Case Study in Xinavane, Mozambique, IHE Delft Institute for Water Education, Delft, the Netherlands, https://waterpip.un-ihe.org/sites/waterpip.un-ihe.org/files/ihe_waterpip_wp1_xinavane_reduced.pdf (last access: 15 July 2021), 2020b.
Clemmens, A. and Molden, D.: Water uses and productivity of irrigation systems, Irrig. Sci., 25, 247–261, 2007.
de Boer, F. and Droogers, P.: Water allocation models for the Incomati River Basin, Mozambique, FutureWater, Wageningen, https://www.futurewater.nl/wp-content/uploads/2016/06/WAM_Incomati_ARA-Sul_FINAL.pdf (last access: 15 July 2021), 2016.
Dejen, Z. A.: Hydraulic and operational performance of irrigation schemes in view of water saving and sustainability: sugar estates and community managed schemes In Ethiopia, CRC Press/Balkema, ISBN 9781138027671, 2015.
den Besten, N., Kassing, R., Muchanga, E., Earnshaw, C., de Jeu, R., Karimi, P., and van der Zaag, P.: A novel approach to the use of earth observation to estimate daily evaporation in a sugarcane plantation in Xinavane, Mozambique, Phys. Chem. Earth Pt. A/B/C, 124, 102940, https://doi.org/10.1016/j.pce.2020.102940, 2020.
den Besten, N., Steele-Dunne, S., de Jeu, R., and van der Zaag, P.: Towards Monitoring Waterlogging with Remote Sensing for Sustainable Irrigated Agriculture, Remote Sens., 13, 2929, https://doi.org/10.3390/rs13152929, 2021.
De Wit, C. T.: Transpiration and crop yields, 0372-6223, Wageningen University, https://edepot.wur.nl/186445 (last access: July 2021), 1958.
Edreira, J. I. R., Guilpart, N., Sadras, V., Cassman, K. G., Van Ittersum, M. K., Schils, R. L., and Grassini, P.: Water productivity of rainfed maize and wheat: A local to global perspective, Agr. Forest Meteorol., 259, 364–373, 2018.
FAO: WaPOR Database Methodology: Level 2, Remote Sensing for Water Productivity Technical Report: Methodology Series, Licence: CC BY-NC-SA 3.0 IGO, FAO, Rome, 63 pp., http://www.fao.org/3/i8225en/I8225EN.pdf (last access: 15 July 2021), 2018.
FAO: Food Outlook – Biannual Report on Global Food Markets, Licence: CC BY-NC-SA 3.0 IGO, FAO, Rome, https://reliefweb.int/attachments/b31d1e40-f571-362e-97d4-faca44f2b0af/CA4526EN.pdf (last access: 22 June 2021), 2019.
FAO: WaPOR V2 quality assessment – Technical Report on the Data Quality of the WaPOR FAO Database version 2, Rome, 89 pp., http://www.fao.org/fileadmin/templates/est/COMM_MARKETS_MONITORING/Sugar/Documents/sugar_assessment_food_outlook_may_2019.pdf (last access: 22 June 2021), 2020a.
FAO: WaPOR database methodology: Version 2 release, April 2020, Rome, 78 pp., https://doi.org/10.4060/ca9894en, 2020b.
FAO: WaPOR: The FAO portal to monitor Water Productivity through Open access of Remotely sensed derived data, FAO, Rome, Italy, https://wapor.apps.fao.org/home/WAPOR_2/1 (last access: 23 October 2020), 2020c.
GDAL: GDAL/OGR contributors Geospatial Data Abstraction software Library, Open Source Geospatial Foundation, https://gdal.org, last access: June 2021.
Griffiths, B. and Lecler, N.: Irrigation system evaluation, Proc. S. Afr. Sug. Technol. Ass., 75, 58–67, 2001.
Hellegers, P., Soppe, R., Perry, C., and Bastiaanssen, W.: Combining remote sensing and economic analysis to support decisions that affect water productivity, Irrig. Sci., 27, 243–251, 2009.
Hess, T. M., Sumberg, J., Biggs, T., Georgescu, M., Haro-Monteagudo, D., Jewitt, G., Ozdogan, M., Marshall, M., Thenkabail, P., and Daccache, A.: A sweet deal? Sugarcane, water and agricultural transformation in Sub-Saharan Africa, Global Environ. Change, 39, 181–194, 2016.
Johnson, M. E., Benin, S., You, L., Diao, X., Chilonda, P., and Kennedy, A.: Exploring strategic priorities for regional agricultural research and development investments in Southern Africa, Intl. Food Policy Res. Inst., http://ebrary.ifpri.org/cdm/ref/collection/p15738coll2/id/127994 (last access: 15 May 2020), 2014.
Karimi, P., David, S., and Bastiaanssen, W. G.: Mapping crop water productivity in the Nile basin through combined use of remote sensing and census data, in: ICID 21st Congress on Irrigation and Drainage: Water Productivity towards Food Security, 15–23 October 2011, Tehran, Iran, New Delhi, India, 137–148, 30-A, 2011.
Karimi, P., Bongani, B., Blatchford, M., and de Fraiture, C.: Global satellite-based ET products for the local level irrigation management: An application of irrigation performance assessment in the sugarbelt of Swaziland, Remote Sens., 11, 705, https://doi.org/10.3390/rs11060705, 2019.
Kharrou, M. H., Le Page, M., Chehbouni, A., Simonneaux, V., Er-Raki, S., Jarlan, L., Ouzine, L., Khabba, S., and Chehbouni, G.: Assessment of equity and adequacy of water delivery in irrigation systems using remote sensing-based indicators in semi-arid region, Morocco, Water Resour. Manage., 27, 4697–4714, 2013.
Mehareb, E. and Galal, M. O.: Ratooning ability of some promising Egyptian sugarcane varieties, J. Sugarcane Tech., 32, 50–59, 2017.
Molden, D. J. and Gates, T. K.: Performance measures for evaluation of irrigation-water-delivery systems, J. Irrig. Drain. Eng., 116, 804–823, 1990.
Molden, D. J., Sakthivadivel, R., Perry, C. J., and De Fraiture, C.: Indicators for comparing performance of irrigated agricultural systems, https://hdl.handle.net/10535/4327 (last access: 15 June 2021), 1998.
Moraes, R. A., Beauclair, E. G. F. D., Gomes, F. H., Moreira, M. M. R., and Maia, R. D.: Identification and estimation of the sugarcane production potential of Mozambique, Afr. J. Agricult. Res., 13, 1924–1931, 2018.
Mul, M, and Bastiaanssen, W.: WaPOR quality assessment: Technical report on the data quality of the WaPOR FAO database version 1.0, Rome, 134 pp., http://www.fao.org/3/ca4895en/CA4895EN.pdf (last access: 1 February 2021), 2019.
Nkamleu, G.: Extensification versus intensification: revisiting the role of land in African agricultural growth, African Economic Conference, Africa Development Bank, Addis Ababa, Ethiopia, https://archive.uneca.org/sites/default/files/uploaded-documents/AEC/2011/nkamleu-extensification_versus_intensification_1.pdf (last access: 21 May 2022), 2011.
Pitts, D., Peterson, K., Gilbert, G., and Fastenau, R.: Field assessment of irrigation system performance, Appl. Eng. Agricult., 12, 307–313, 1996.
Playán, E. and Mateos, L.: Modernization and optimization of irrigation systems to increase water productivity, Agr. Water Manage., 80, 100–116, 2006.
Raes, D., Steduto, P., and Hsiao, T. C.: Reference manual, Chapter 2, AquaCrop model, Version 4.0, Food and Agriculture Organization of the United Nations, Rome, Italy, https://www.fao.org/3/bc122e/bc122e.pdf (last access: 20 June 2020), 2013.
Smith, M.: CROPWAT: A computer program for irrigation planning and management, 46, Food and Agriculture Organization, Rome, Italy, https://books.google.nl/books?id=p9tB2ht47NAC&printsec=frontcover&source=gbs_ge_summary_r&cad=0 (last acess: 22 June 2021), 1992.
Steduto, P. and Albrizio, R.: Resource use efficiency of field-grown sunflower, sorghum, wheat and chickpea: II. Water use efficiency and comparison with radiation use efficiency, Agr. Forest Meteorol., 130, 269–281, 2005.
Steduto, P., Hsiao, T. C., and Fereres, E.: On the conservative behavior of biomass water productivity, Irrig. Sci., 25, 189–207, 2007.
Steduto, P., Hsiao, T. C., Raes, D., and Fereres, E.: AquaCrop – The FAO crop model to simulate yield response to water: I. Concepts and underlying principles, Agron. J., 101, 426–437, 2009.
Sutton, P. C., Anderson, S. J., Costanza, R., and Kubiszewski, I.: The ecological economics of land degradation: Impacts on ecosystem service values, Ecol. Econ., 129, 182–192, https://doi.org/10.1016/j.ecolecon.2016.06.016, 2016.
Tongaat Hullet: 2018 Integrated annual report: Value creation for all stakeholders through an all-inclusive approach to growth and development, Tongaat Hulett, South Africa, https://www.tongaat.com/annual_reports/ar_2018/downloads/AR-2018-FULL.pdf (last acess: 10 May 2020), 2018.
Van der Zaag, P. and Carmo Vaz, Á.: Sharing the Incomati waters: cooperation and competition in the balance, Water Policy, 5, 349–368, 2003.
Villalobos, F. J. and Fereres, E.: Principles of agronomy for sustainable agriculture, Springer, New York, USA, ISBN 978-3-319-46116-8, 2016.
Yilma, W. A.: Computation and spatial observation of water productivity in Awash River Basin, UNESCO-IHE, Delft, https://ihedelftrepository.contentdm.oclc.org/digital/collection/masters1/id/50110/ (last access: 15 July 2020), 2017.
Zhang, Y., Peña-Arancibia, J. L., McVicar, T. R., Chiew, F. H., Vaze, J., Liu, C., Lu, X., Zheng, H., Wang, Y., and Liu, Y. Y.: Multi-decadal trends in global terrestrial evapotranspiration and its components, Scient. Rep., 6, 1–12, 2016.
Zwart, S. J. and Bastiaanssen, W. G.: Review of measured crop water productivity values for irrigated wheat, rice, cotton and maize, Agr. Water Manage., 69, 115–133, 2004.