
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)12381-X
10.1016/j.heliyon.2024.e36350
e36350
Research Article
Comprehensive framework for interpretation of WaPOR water productivity
Veysi Shadman sh.veysi@areeo.ac.ir
a⁎
Galehban Eslam b
Nouri Milad a
Mallah Sina a
Nouri Hamideh c
a Soil and Water Research Institute, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran
b Department of Remote Sensing and GIS, Faculty of Geography, University of Tehran, Tehran, Iran
c Water Science and Policy Coordination, DCCEEW, Australian Government, Australia
⁎ Corresponding author. sh.veysi@areeo.ac.ir
15 8 2024
30 8 2024
15 8 2024
10 16 e363504 3 2024
14 8 2024
14 8 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
This study presents a comprehensive framework for analyzing water productivity products provided by the FAO Water Productivity Open-access portal (WaPOR), focusing on various crops cultivated in both rainfed and irrigated areas within a semi-arid basin in Iran. Two indices, namely Gross Water Productivity (GWP) and Net Water Productivity (NWP), were introduced to quantify water productivity across crop fields. However, these indices may mislead decision-makers, because they aggregate water productivity for all crops and exacerbate the challenges posed by water scarcity. Therefore, mapping crop types seems necessary to enhance the interpretation of these indices and develop a dimensionless index for comparing different crops. The results demonstrated a fundamental change when comparing dimensionless water productivity with GWP and NWP products. Surprisingly, some pixels initially exhibiting high water productivity ranked as low water-productive pixels based on the derived dimensionless index, and vice versa. Based on dimensionless indicators, rainfed crops, particularly rainfed cereals, ranked as the most water-productive crops. The areas with dimensionless values below 0.5 warrant heightened attention to curtail non-beneficial water consumption and elevate water productivity. This research emphasizes the significance of mapping cultivation types as supplementary layers to facilitate precise, data-driven decision-making and enable comparisons of crops based on dimensionless water productivity indices.

Keywords

Crop type
Dimensionless indices
Machine learning
Support layers
Water scarcity
==== Body
pmc1 Introduction

Water scarcity is commonly defined as a situation where demand for water surpasses available supply [1]. In a broader context, economic water scarcity arises from water mismanagement, weak governance, insufficient investment in water infrastructure, or a lack of human capacity, all of which constrain access to water [2,3]. Agricultural economic water scarcity was also defined as the lack of irrigation not due to hydrological limitations but rather because of restricted institutional and economic capacity [4].

Water scarcity is a critical limit that endangers food security in drylands and is expected to worsen due to climate change and population growth by 2050 [[5], [6], [7]]. Meeting the challenge of producing more food for more people under harsher climatic conditions will require optimizing the water-food-energy nexus, improving water productivity (WP), and expanding cropland. However, land expansion is not a practical solution in regions already experiencing water stress. Therefore, increasing WP has been recognized as a practical strategy to achieve food security [8] and aligns with the United Nations Sustainable Development Goals (SDG 6.4) [9,10].

The physical WP is simply equals the ratio of crop yield to water consumed. Multiple definitions of water productivity (for WP) are available, with another important one being economic water productivity. It can be defined as the ratio of the value of agricultural benefit to the amount of water consumed or used in the production process. This metric assesses the economic efficiency of water use in agriculture by considering the monetary return obtained from the water input [11]. Nutritional WP also offers a valuable metric for assessing WP by quantifying the nutrient output generated per unit of water used [12]. Nevertheless, physical WP, commonly referred to as “crop per drop”, remains the most widely recognized definition of water productivity.

Despite its simplicity, several challenges are associated with accurate interpretation and calculation of physical WP. One of the challenges is the misinterpretation of WP, which assumes that increasing its value leads to water savings in water-scarce regions. However, improving WP does not necessarily save water but instead can lead to an increase in water consumption, known as the Jevons’ paradox [13,14].

The scale dependency of WP is one of the main reasons causing miscalculation and misinterpretation. Since actual evapotranspiration (ET) is hardly measured at the field scale, WP is often estimated using inflow (irrigation + effective rainfall). WP is, therefore, obtained based on the water inflow in fields. In this case, the summation of water inflow volumes at fields situated in a given region may not be equal to total inflows at the catchment scale. This can be ascribed to the fact that a portion of return flows (non-consumptive recoverable deep percolation + runoff) in a field may become inflow in an adjacent field [15]. This explains why water-saving at the field scale cannot necessarily translate into water-saving at the regional scale, known as real water-saving, particularly in closed catchments. In other words, by enhancing the WP values at the field scale, one may expect WP increment but not real water-saving at the basin-scale. Therefore, this misinterpretation may result in misguided and flawed policies in water-stressed areas wherein water-saving should be at the top of the agenda.

Comparing WP of different crops cultivated in a region is another pitfall [16]. In some cases, water productivity is represented for a large area (e.g., a basin scale), which is common for policymakers to calculate regional WP by dividing the total crop production to total water use in a basin [17]. However, the numerator encompasses a broad spectrum of crops, and aggregating all crop productions to compute WP at a regional scale, may be misleading [18]. Neglecting other factors than water affecting the WP can also cause misinterpretation of WP values [19]. For instance, WP of a given crop in a specific region is affected by crop type, cultivar, climatic factors, management, soil condition, etc. [20,21]. When defining the WP gap, it is important to note that these gap values should be treated regionally, and caution should be exercised when recommending or comparing WP values obtained under different regions. WP obtained for a high-yielding humid part may not be suitable as a norm for a low-yielding hyper-arid area.

Alternative gridded datasets have been widely applied in water management to provide invaluable data in data-poor areas [22,23]. These datasets can also be helpful in determining and interpreting WP, as they provide necessary data (e.g., evapotranspiration components and biomass) and information on other factors affecting WP, such as soil, climate, topography, and management [[24], [25], [26], [27], [28], [29]]. However, uncertainties in gridded products can affect the accuracy of WP results, which can impact decision-making processes. Therefore, caution must be exercised when using gridded products to estimate WP in different basins with varying crop cultivation under rainfed and irrigated conditions.

There are limited studies on the application of gridded datasets to analyze WP. One such study, referenced as [30], suggested calculating dimensionless WP for wheat, rice, and maize on national and global scales. Initially, they merged crop dominance maps [31] with Monthly Irrigated and Rainfed Crop Areas from MIRCA-2000 [32] at a 10 km spatial resolution. Subsequently, employing satellite-derived Normalized Difference Vegetation Index (NDVI) data, they refined and disaggregated global maps of cereals to a 1 km spatial resolution. Then, they developed a formula to estimate WP based on remotely-sensed parameters such as the fraction of Absorbed Photosynthetically Active Radiation and surface albedo. Notably, their methodology relied on broad spatial and temporal scales for crop mapping, utilizing NDVI data at 10-day intervals spanning from 1998 to 2008. Consequently, their findings reflect an average condition over a decade, overlooking annual variations due to weather anomalies or anthropogenic impacts. This approach is best suited for regions with contiguous crop coverage, typified by large fields aligning with the 100 ha spatial resolution (1 km × 1 km) of crop type maps. However, it may prove inadequate for finer-scale assessments in regions with fragmented crop coverage or small fields indicative of subsistence agriculture. Furthermore, the use of satellite-derived crop maps employing unsupervised mapping techniques may introduce significant uncertainty, particularly in regions with small-scale fields and notable spatial variations in crop types.

Given these limitations, there is a pressing need for WP studies conducted at finer spatial scales to accurately capture both the spatial and annual variations in WP. This need is particularly relevant in regions like Iran, where small-scale farming predominates, and neighboring fields often exhibit diverse crop types. To address this research gap, our study aims to estimate WP within a catchment characterized by small-scale fields and annual variations, utilizing gridded products. Additionally, we introduce a dimensionless WP metric to enhance the interpretation of our results for water management purposes.

2 Material and methods

2.1 Study area

Honam is a catchment located in the Lorestan province of western Iran, covering 14,000 ha. The common crops are beans, chickpeas, sugar beet, and winter grain (i.e., wheat, barely). The average elevation of this catchment is approximately 1600 m above sea level, and the region receives an average annual precipitation of 500 mm. Fig. 1.Fig (1) Location of the Honam catchment.

Fig (1)

The crop calendar in the Honam region outlines the schedules for planting and harvesting the main crops. Among the various crops grown in the region wheat and barley are categorized as fall crops. Sugar beet is classified as summer crops. Alfalfa and clover are perennial crops that continue to grow for several years after planting. Orchards are considered permanent crops. It is important to note that in the Honam catchment area, there is no double-cropping on each farmland. Fig. 2 shows the crop calendar of Honam catchment for main crops.Fig (2) Crop calendar of the main crops in Honam catchment.

Fig (2)

2.2 Dataset and preprocessing

Land Use and Land Cover (LULC) changes using the European Space Agency's (ESA) high-resolution satellite imagery with a spatial resolution of 10 (m) were surveyed in this research (https://www.esa-landcover-cci.org/) for detection of cropland areas. Preprocessing and extraction of cropland were performed within Google Earth Engine (GEE) Fig. 3. Therefore, crop land areas delineated from other land use categories.Fig (3) Maps of Honam 10 m resolution, Land Use Land Cover map.

Fig (3)

Meanwhile, cropland was reprojected to a 250 (m) as target resolution. This is necessary because the Gross Water Productivity (GWP) and Net Water Productivity (NWP) were obtained from the open-access WaPOR dataset, which have a pixel size of 250 × 250 (m). The WaPOR dataset (https://wapor.apps.fao.org/) provides up-to-date information on the Near East and Africa region, including the Honam catchment. The WaPOR dataset contains several components related to water productivity, including the ratio of biomass production to total evapotranspiration, evaporation from soil, canopy transpiration, and interception. GWP and NWP take into account the water consumed solely for evapotranspiration plus interception and transpiration respectively.(1) GWP=TBPE+T+I

(2) NWP=TBPT

where TBP is total biomass production in kgDM/ha. E is evaporation, T is transpiration and I is interception, all in mm. ETIa-WaPOR is derived from the ETLook model. ETLook uses a simplified Penman-Monteith model to estimate ETIa, which is based on input data such as vegetation indices, land surface temperature, and meteorological data [33]. This model separates evaporation, transpiration, and interception using equations (3), (4), (5).(3) E=Δ(Rnsoil−G)+ρairCp(es−ea)rasoilΔ+γ(1+rssoilrasoil)

(4) T=Δ(Rncanopy−G)+ρairCp(es−ea)racanopyΔ+γ(1+rscanopyracanopy)

(5) I=0.2ILAI(1−11+CvegPCP0.2ILAI)

where E and T (mm.day−1) are the evaporation and transpiration, respectively, Rn (MJ m−2 day−1) of the soil (Rn,soil) and canopy (Rn,canopy) is the net radiation, G (MJ m−2 day−1) is the ground heat flux, ρa (kg.m−3) denotes the density of air, CPP (MJ kg−1°C−1) stands for the specific heat of air, (es-ea) (kPa) represents the vapor pressure deficit (VPD), ra and rs (s/m) stand for the aerodynamic soil or canopy resistance, respectively, Δ (kPa°C−1) is the slope of the curve relating saturated water vapor pressure to air temperature, and γ denotes the psychrometric constant (kPa°C−1). Interception (mm.day−1) is a function of vegetation cover, LAI and precipitation [34]. The summary characteristics of the dataset used in the research are shown in Table .1.Table 1 Data employed in this research.

Table 1Data type	symbol	Source	Spatial resolution	Temporal resolution	
Digital elevation model	DEM	SRTM	90 (m)	*	
Land use Land Cover	LULC	ESA	10 (m)	Annual	
Gross biomass water productivity	GWP	WaPOR	250 (m)	Annual	
Net biomass water productivity	NWP	WaPOR	250 (m)	Annual	
Precipitation	*	CHIRPS	5500 (m)	Daily	

2.3 Crop type mapping

Various remote sensing data classification methods are used to identify crop types, including Support Vector Machines [35], decision trees [36], and Random Forest [37]. Random Forest is a popular classifier that is well-suited for data distribution [38]. Therefore, the Random Forest method was used to distinct irrigated and rainfed areas and detect different crop types using vegetation indices derived Landsat-8 data coded in Google Earth Engine (GEE). Accordingly, forty-two Landsat-8 images, each with a spatial resolution of 30 m, captured in 2020 and exhibiting less than 5 % cloud coverage, were chosen from the GEE online platform to comprehensively cover the entire study area. The map from a 30-m resolution obtained from Landsat was resampled to 250 m using nearest neighbor [22,23]. It is worth noting that given the scale of field areas in our study region, upscaling the resolution from 30 m to 250 m appears reasonable, as supported by field observations. Each pixel contains either spring or winter crops, indicating that each grid is cultivated with only one crop per year. The maps for winter and spring crops were meticulously separated and integrated to derive the annual crop map. Fallow pixels were excluded to maintain accuracy. Eight vegetation indices have been used for supervising different crops, as presented in Table .2.Table (2) Vegetation indices used from Landsat 8 for crop-mapping.

Table (2)Index	Equation	
Normalized Difference Vegetation Index (NDVI)	(6) (RNIR−RRed)/(RNIR+RRed)

	
Green Normalized Difference Vegetation Index (GNDVI)	(7) (RNIR−RGreen)/(RNIR+RGreen)

	
Green Red Vegetation Index (GRVI)	(8) (RGreen−RRed)/(RGreen+RRed)

	
Green-blue vegetation index (GBVI)	(9) (RGreen−RBlue)/(RGreen+RBlue)

	
Normalized Difference Moisture Index (NDMI)	(10) (RNIR−RSwir1)/(RNIR+RSwir1)

	
Normalized Burn Ratio 2 (NBR2)	(11) (RSwir1−RSwir2)/(RSwir1+RSwir2)

	
Enhanced vegetation index (EVI)	(12) 2.5(RNIR−RRed)%((RNIR+6×RRed‐7.5×RBlue)+1

	
Green Chlorophyll Vegetation Index (GCVI)	(13) (RNIR−RGreen)−1

	

Bands 5, 4, 6 and 7 of Landsat-8 were represented as RNIR, RRed, RSwir1, and RSwir2, respectively, while bands 2 and 3 were represented as RBlue and RGreen, respectively.

Field sample collection is a pivotal component of the crop mapping procedure, demanding timely and meticulous execution. To ensure accuracy, we first identified hot-spot agricultural zones where various major crops are cultivated. Subsequently, 286 field sample collections were concentrated in these identified areas. Field sample collectors were equipped with appropriate instruments (e.g., GPS devices) organized for ground sampling. Samples were divided into two groups: 132 points in irrigated areas and 154 in rain-fed cultivation. 70 % of the sample data was randomly selected for training the RF classification algorithms, while the remaining 30 % was used for testing and validation. The classification performance was independently validated for each crop type using a confusion matrix. The accuracy of the classification was assessed using Producer's Accuracy (PA), User Accuracy (UA), and Kappa coefficient (K) [39]:(14) UA=100%−CommissionError

(15) PA=100%−OmissionError

(16) Kappa=N∑i=1nXii−N∑i=1nX+iXi+N2−∑i=1nX+iXi+

Where Omission Errors occur if a pixel is left out of the class type being assessed, Commission Errors happen if a pixel is incorrectly involved in the class type being assessed, X + i: marginal sum of the columns, Xi+: marginal sum of the rows, N: total number of pixels in the confusion matrix, n: number of classes, and Xii: correctly classified pixels in each class.

2.4 Dimensionless water productivity

Water productivity in agriculture can be challenging due to the differences in the dimensions of the numerator and denominator [30]. To ensure a valid comparison between crops, it is essential to compare each crop within its bright and hot spots [40]. Bright and hot spots refer to crops performing better or worse than their neighbors under similar conditions. The threshold for defining better performance is arbitrary, and different studies have used different percentile values. For instance Ref. [41], used the 95th percentile by selecting the upper 5 %, while [42] defined productivity as split into four different quadrants. The process involved the detection of bright and hot pixels for each crop based on the crop type of the region, which was then used to calculate a dimensionless water productivity index (DWPI). This approach scales down features between 0 and 1. To make water productivity a useful performance indicator within this research, dimensionless GWP and NWP in Eqs. (17), (18) have been defined.(17) DGWP=GWP−HotpixelBrightpixel−Hotpixel

(18) DNWP=NWP−HotpixelBrightpixel−Hotpixel

where GWP represents the water productivity of each pixel for each crop, 'Hot pixel' and 'Bright pixel' are the pixels with a local minimum and maximum average of 5 % of pixel water productivity for each crop, respectively.

DWPI is inspired by the water productivity score provided by Ref. [30], albeit with a different range. It is noteworthy that normalization is a common approach used to standardize various factors, making them dimensionless and comparable. For example, vegetation indices [43] and soil moisture content [44] are commonly formulated in a standardized manner. DWPI is a suitable for both local and regional assessments, allowing for comparisons at various scales without being influenced by specific units or environmental specifics. GWP and NWP layers were superimposed with the resize crop maps of the catchment to determine the contribution of each crop to overall productivity. Furthermore, given the potential impact of extreme values on the normalization and dimensionless representation of WP results, outliers for each crop were identified using a normal distribution. To better illustrate and clarify the methodological steps, the conceptual framework is depicted in Fig. 4.Fig (4) Framework of research.

Fig (4)

3 Results and discussion

3.1 Mapping cultivation types and crop type

The map of irrigated and rainfed areas is depicted in Fig. 5, generated using the random forest method. The accuracy of detection, based on supervised point access, reached 91 % for irrigated areas and 95 % for rainfed areas. Notably, agricultural land encompasses 34.7 % of the Honam catchment, with 16.6 % and 18.2 % dedicated irrigated and rainfed cultivation, respectively. The prevalence of rainfed lands can be attributed to sufficient precipitation level in our study area.Fig (5) Distribution of sample points versus irrigated and rainfed areas of Honam catchment.

Fig (5)

The spatial distribution revealed a concentration of irrigated regions in the lower elevations of the catchment, while rainfed regions were predominantly situated at higher altitudes. This observed pattern may be influenced by the availability of water resources, which tend to be more abundant in lower elevations. These findings align with the work of [45], who substantiated the relationship and correlation between elevation and irrigation zones. Variations in elevation and access to water resources are pivotal factors influencing the distribution of irrigated lands [46].

To accurately produce the annual crop map Fig. 7, we first mapped crop patterns for two distinct seasons: season 1 (fall/winter crops) and season 2 (spring/summer crops). There are two different seasonal cropping systems in our study area, making it impractical to rely on a single annual cropping map due to significant uncertainties. Therefore, the final annual cropping map was created by overlapping the two seasonal maps depicted in Fig. 6. This approach highlights that the annual cropping pattern is a composite of different seasonal maps. It is important to note that no fields in the study area are cropped in both seasons due to phenological and resources limitations.Fig (6) Seasonal cropping maps A) Fall/Winter cultivated crops, B) Spring/Summer cultivated crops.

Fig (6)

Fig (7) Annual cropping map.

Fig (7)

Fig. 7 illustrates the crop type distribution in Honam catchment. Irrigated chickpea and orchard crops were excluded from the analysis as they were relatively small. The map indicates that the area of rainfed barley is more extensive than that of other crops. However, field observations suggest that this area is a combination of both rainfed wheat and barley. Table .3 shows detail of each crop area in the catchment.Table (3) Area of each crop in Honam catchment.

Table (3)No	Crop	Area (ha)	No	crop	Area (ha)	
1	Alfalfa	214	7	Beans	79	
2	Barley Irrigated	253	8	Wheat irrigated	388	
3	Clover	158	9	Chickpea Rainfed	731	
4	Orchard	71	10	Wheat Rainfed	71	
5	Sugar beet	122	11	Barley Rainfed	1443	
6	Chickpea Irrigated	16	12	–		

Distinguishing between wheat and barley cultivations was difficult due to their similar phenology and planting schedules [47]. The study confirmed that high-resolution images are required to differentiate between cereals. Given the close alignment of planting, growing, and harvesting periods, as well as similar spectral reflection characteristics in wheat and barley, water productivity analysis for these crops was conducted by categorizing them collectively as “rainfed cereals.”

Table .4 shows the user's accuracy (UA) and producer's accuracy (PA) for each crop. According to the accuracy assessment metrics, beans had the highest UA, whereas Chickpea rainfed had the lowest significant UA. However, sugar beet had the lowest PA compared to other crops.Table (4) Accuracy assessment of crop type mapping.

Table (4)Crop	User Accuracy (UA)%	Producer Accuracy (PA)%	Sample Count	
Alfalfa	0.75	0.88	28	
Barely Irrigated	0.88	0.84	20	
Barely Rainfed	0.78	0.86	40	
Beans	1.00	1.00	10	
Chickpea Rainfed	0.46	0.88	37	
Clover	0.80	0.75	18	
Orchard	0.73	0.73	24	
Sugar beet	0.88	0.63	11	
Wheat irrigated	0.79	0.81	39	

Sugar beet is typically grown in various regions of Honam during the spring, but it is cultivated in small areas widely dispersed. As a result, obtaining sufficient ground samples of these classes for the network was challenging due to the difficulty of fieldwork [39]. Emphasized the importance of a sampling scheme while conducting the accuracy assessment. After calculating the Kappa statistic for the entire map, the kappa accuracy (KA) was found to be 73 %. The Kappa values were calssified as follows: <40 % as poor, 40–55 % as fair, 55–70 % as good, 70–85 % as very good, and >85 % as excellent [48]. The object-based approach yielded an average user accuracy index of approximately 79 %, suggesting that the provided cropping map is accurate. However, resampling data from a finer spatial resolution (30-m Landsat) to a coarser spatial resolution (250-m WaPOR) can introduce uncertainties. This resampling process can result in the loss of finer spatial details and obscure small-scale variations. Additionally, it can alter the spatial patterns and variability captured in the data, potentially impacting crop type mapping and water productivity results. Bilinear interpolation approach may also smooth out extremes, leading to errors. Consequnelty, we deterined the Kappa vaule for the resampled map. Accordingly, there was a Kappa value of 61 % for the resampled crop type map. Thus, despite a slight decrease in accuracy, the Kappa value obtained from the resampled cropping map, which range from 55 % to 70 %, indicate good accuracy of the results.

3.2 WaPOR water productivity

Fig. 8 displays NWP and GWP of the WaPOR dataset. The results have been divided into quartiles, which means that their values have been split into four equal parts, with each part containing 25 % of the total values. The values ranging from 0 to 25 % are in the first class, values ranging from 25 to 50 % are in the second class, and so on for the other two types.Fig (8) NWP and GWP of the Honam catchment.

Fig (8)

Overall based on raw data of WaPOR dataset, the upstream of the catchment exhibits superior water productivity compared to the downstream area based on both indices. This may be due to various factors, such as differences in field water management practices by farmers, variations in soil quality across different areas, reduced water volume downstream as it flows from upstream [49]. Emphasized that the decreased water availability may lead to reduced water productivity.

Also, Fig. 8 displays NWP being almost double GWP and a significant difference was observed. This highlights that soil evaporation (E) and interception (I) have the main role in reducing the values of water productivity. Therefore, a solution to improve water productivity in the region is to implement vapor shift [50]. To achieve this goal, management options such as increasing soil infiltration [51], enhancing water retention capacity, and leveling the lands should be considered [52].

Fig. 8 shows that some pixels with high Gross Water Productivity (GWP) values exhibit good water productivity (Q3 and Q4), whereas when using the Net Water Productivity (NWP), these same pixels are categorized as having a bad status (Q1 and Q2). When the Gross Water Productivity (GWP) is in the upper quartiles (Q3 and Q4) and the Net Water Productivity (NWP) is in the low quartiles (Q1 and Q2), the amount of total evapotranspiration (ETIa) was low. Therefore, NWP is high because maximum transpiration occurs when ETIa is low. In other words, the portion of transpiration is greater than that of evaporation (E). In addition, where the NWP proves to be suitable (Q3 and Q4) while the Gross Water Productivity GWP fails (Q1 and Q2), this means that the portion of transpiration is low as beneficial water consumption for that specific pixel.

Overall, these observations underscore the significance of considering the partitioning map of transpiration to evapotranspiration (T/ET) as an informative layer contingent in the particular type of cultivation for informed decision-making. Fig. 9 shows partitioning evapotranspiration (T/ET) for both irrigated and rainfed zones. Fig. 9 depicted T/ET ranges in (0.36–0.86) for irrigated areas, (0.39–0.80) for rain-fed areas. The results reported in this study align with the broad spectrum established by numerous prior global investigations [[53], [54], [55]].Fig (9) T/ET of irrigated and rain-fed areas in the Honam catchment.

Fig (9)

These results reveal no significant difference in the T/ET ratio between rain-fed and irrigated cultivation, and they oscillate within the same range. The increase in T/ET for rainfed cultivation is linked to a decline in factors such as soil moisture and rainfall events in water limited regions [56]. Under water-limited conditions, crops tend to increase their water uptake for transpiration to mitigate the impact of water restriction [57]. To further illuminate this matter, the standardized precipitation index (SPI) of the region is depicted in Fig. 10 as an indicator of droughts. This choice stems from the fact that SPI exhibits the strongest correlations between drought intensity and WP anomalies, which can be observed during or after drought events [58]. Specifically, SPI-12 was calculated based on CHIRPS precipitation dataset (1990–2021) to assess the drought status in 2020.Fig (10) SPI variation in the Honam catchment. The upper and lower red line start with wet, and drought; second, horizontal axis shows the monthly precipitation of the region.

Fig (10)

Fig. 10 shows that Honam experienced a very wet year in 2020, with precipitation levels reaching approximately 650 (mm). Therefore, the high rainfed productivity this year could be directly influenced by the significant amount of precipitation received. These findings are consistent with previous comparative studies conducted across different regions, which have consistently demonstrated the impacts of drought on rainfed systems in drylands [[59], [60], [61]].

3.3 Water productivity based on crop type

Fig. 11 illustrates NWP and GWP for various crops. The results indicate that water productivity for rainfed crops fluctuates more than for irrigated crops. Crops with a high degree of fluctuation suggest less satisfactory management of crops and resources (such as water). Conversely, areas with low fluctuation typically experience more effective management [62]. Irrigated crops are less affected by changes because farmers have more control over water availability [63,64].Fig (11) GWP and NWP for each crop.

Fig (11)

Fig. 11 illustrates that the difference between the indices for alfalfa, clover, and sugar beet is minimal. When both indices yield similar values, it suggests that the proportion of evaporation (E) as non-beneficial water consumption is lower in these crops compared to others. In the categorization of water-saving strategies, this result comes from deficit irrigation practices [65] or less frequent irrigation [66], which lead to reduced soil evaporation. Also, knowing that these crops command higher prices in the market, farmers might be more inclined to implement stricter management practices, such as precise irrigation and specialized fertilizers, to protect water productivity oscillation [67,68]. Although Fig. 11 shows the water productivity of each crop, but is not easier identification which crop of efficient water use practices. Therefore, it is necessary to consider the dimensionless method for each crop Fig. 12.Fig (12) Box plot of dimensionless GWP and NWP for each crop.

Fig (12)

Fig. 12 illustrates the dimensionless productivity indices (DGWP and DNWP) for crops that were cultivated in Honam. According to Ref. [69], the criterion for comparing the productivity of crops was set at 0.5. Therefore, any crop with an average value of less than 0.5 is considered to have poor water productivity and should be improved. Based on DWP indices, cereals and chickpeas rainfed are ranked as the most productive crops. The results show that the dimensionless indices change the ranking of crops, and consequently change water productivity interpretation and cereals rainfed shown better performance than irrigated crops.

3.4 Spatial mapping of DWP

Fig. 13 presents a regional map of DWP, which is classified into four clusters based on the WP values. The DWP values are classified as follows: (0–0.25) for low productivity, (0.25–0.5) for moderate productivity, (0.5–0.75) for good productivity, and (0.75–1) for excellent productivity. A regional map of DWP assists in pinpointing areas where water resources are effectively utilized and highlights locations where enhancements are feasible [64].Fig (13) Distribution of DWPI in the Honam plain.

Fig (13)

Fig. 13, reveals significant differences compared to Fig. 9, which shows raw water productivity. Several regions depicted as exhibiting low water productivity (0–0.25) in the first quartile of Fig. 8 are reclassified as demonstrating moderate water productivity when utilizing the DWP indices.

In regions with high water productivity, there is not always a corresponding high DWP. Improving water productivity in these areas might actually lead to increased water consumption. Conversely, some areas with low water productivity may exhibit higher DWP. In regions facing water scarcity, the narrative of water productivity becomes closely intertwined with the challenge of managing water consumption. It is important to recognize that this does not always translate to reduced water consumption. Therefore, solely focusing on water productivity increases paints an incomplete picture, especially in water-scarce regions [[70], [71], [72]].

While decision-making based solely on the DNWP seems appealing due to its consideration of transpiration alone, the DGWP, which takes into account both beneficial and non-beneficial water consumption, is more practical and advisable for the decision process. The DNWP can help DGWP as it is associated with which pixel zone, and each crop better uses water, but for the final decision it is better to decide with DGWP.

Using these maps, regions characterized by low water productivity can be identified, and by correlating them with spatial data such as soil type, irrigation management, and other relevant factors, efforts can be made to uncover the underlying causes. It can be employed to distinguish between farm plots exhibiting effective and inadequate management practices.

3.5 Limitations

Several limitations contributed to the uncertainties in our results. Firstly, some classes were not differentiated with the desired level of accuracy due to challenges such as pronounced temporal overlap between planting, growth, and harvesting phases. This results in spectral reflection similarities between crops like wheat and barley, leading to the inevitable combination of some classes and introducing uncertainties. Moreover, another constraint was the sample quantity for each crop in crop type mapping. Acquiring sufficient samples for each crop type is challenging due to the heterogeneous distribution of crops and the costs and logistics associated with data collection. Insufficient sample sizes may introduce inaccuracies in classification algorithms, potentially impacting the accuracy and reliability of crop identification and mapping. Additionally, resampling to coordinate the WaPOR and Landsat products introduces another source of error, albeit relatively low, in the results. This process can affect the finer spatial details and alter the spatial patterns and variability captured in the data, contributing to the overall uncertainty.

Future studies can reduce uncertainties and enhance the reliability of results by exploring alternative resampling algorithms or incorporating machine learning techniques that more effectively preserve spatial details during the resampling process.

4 Conclusions and remarks

This study aims to demonstrate how decision-makers can apply the WaPOR water productivity (WP) dataset. The results indicate that simply presenting a raw value may be misleading. Therefore, it is essential to consider the type of cultivation (irrigated and rainfed) and crop type mapping as a supporting layers to avoid erroneous conclusions that could mislead policymakers and decision-makers. Meanwhile, when comparing two crops in the same area, it is essential to first normalize WP to make it dimensionless. In fact, the dimensionless index is one of several indicators that can be used for the physical productivity of a crop compared to another ones. Moreover, in arid and semi-arid regions, efforts to enhance water productivity may not prioritize water consumption management. This could potentially lead to an unintended increase in water use. This suggests that solely pursuing high productivity, especially in arid and semi-arid regions, can exacerbate the risk of water insecurity in these areas. Datasets such as WaPOR have the potential to offer valuable insights to help decision-makers in the water resources sector of agriculture. However, in future versions, enhancing support layers can make this dataset even more useful.

Data availability

Data will be made available on request.

Funding

This research received no external funding.

CRediT authorship contribution statement

Shadman Veysi: Writing – original draft, Supervision, Project administration, Methodology, Formal analysis, Conceptualization. Eslam Galehban: Visualization, Software. Milad Nouri: Writing – review & editing, Resources, Investigation. Sina Mallah: Software, Data curation. Hamideh Nouri: Writing – review & editing, Methodology.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

We would like to express our sincere appreciation and gratitude to the Food and Agriculture Organization of the United Nations (FAO) for providing access to the WaPOR (Water Productivity Open-access Portal) dataset. The efforts of the FAO in collating and maintaining this dataset have been instrumental in advancing scientific understanding and promoting evidence-based decision-making in water management.
==== Refs
References

1 Schyns J.F. Hoekstra A.Y. Booij M.J. Review and classification of indicators of green water availability and scarcity Hydrol. Earth Syst. Sci. 19 2015 4581 4608 10.5194/hess-19-4581-2015
2 Molden D. Water for Food Water for Life: A Comprehensive Assessment of Water Management in Agriculture first ed. 2007 10.4324/9781849773799
3 Vallino E. Ridolfi L. Laio F. Measuring economic water scarcity in agriculture: a cross-country empirical investigation Environ. Sci. Pol. 114 2020 73 85 10.1016/j.envsci.2020.07.017
4 Rosa L. Chiarelli D.D. Rulli M.C. Dell'Angelo J. D'Odorico P. Global agricultural economic water scarcity Sci. Adv. 6 2020 eaaz6031
5 Nouri M. Homaee M. Contribution of soil moisture variations to high temperatures over different climatic regimes Soil Tillage Res. 213 2021 105115 10.1016/j.still.2021.105115
6 Mc Carthy U. Uysal I. Badia-Melis R. Mercier S. O'Donnell C. Ktenioudaki A. Global food security – issues, challenges and technological solutions Trends Food Sci. Technol. 77 2018 11 20 10.1016/j.tifs.2018.05.002
7 Cosgrove W.J. Loucks D.P. Water management: current and future challenges and research directions Water Resour. Res. 51 2015 4823 4839 10.1002/2014WR016869
8 Nazari B. Liaghat A. Akbari M.R. Keshavarz M. Irrigation water management in Iran: implications for water use efficiency improvement Agric. Water Manag. 208 2018 7 18 10.1016/j.agwat.2018.06.003
9 Whitcraft A.K. Becker-Reshef I. Justice C.O. Gifford L. Kavvada A. Jarvis I. No pixel left behind: toward integrating Earth observations for agriculture into the united Nations sustainable development goals framework Rem. Sens. Environ. 235 2019 111470 10.1016/j.rse.2019.111470
10 Caldera U. Breyer C. Strengthening the global water supply through a decarbonised global desalination sector and improved irrigation systems Energy 200 2020 117507 10.1016/j.energy.2020.117507
11 Molden D. Oweis T. Steduto P. Bindraban P. Hanjra M.A. Kijne J. Improving agricultural water productivity: between optimism and caution Agric. Water Manag. 97 2010 528 535 10.1016/j.agwat.2009.03.023
12 Drastig K. Singh R. Telesca F.-M. Carra S. Jordan J. A review of nutritional water productivity (NWP) in agriculture: why it is promoted and how it is assessed? Water 15 2023 10.3390/w15244278
13 Grafton R.Q. Williams J. Perry C.J. Molle F. Ringler C. Steduto P. Udall B. Wheeler S. Wang Y. Garrick D. The paradox of irrigation efficiency Science 361 2018 748 750 30139857
14 Wheeler S.A. Carmody E. Grafton R.Q. Kingsford R.T. Zuo A. The rebound effect on water extraction from subsidising irrigation infrastructure in Australia Resour. Conserv. Recycl. 159 2020 104755 10.1016/j.resconrec.2020.104755
15 Ren D. Xu X. Engel B. Huang Q. Xiong Y. Huo Z. Huang G. A comprehensive analysis of water productivity in natural vegetation and various crops coexistent agro-ecosystems Agric. Water Manag. 243 2021 106481 10.1016/j.agwat.2020.106481
16 Seckler D. Revisiting the “IWMI paradigm”: increasing the efficiency and productivity of water use Zeitschrift für Bewässerungswirtschaft 34 1999 85 94
17 T. Oweis, A. Hachum, J. Kijne. Water Harvesting and Supplementary Irrigation for Improved Water Use Efficiency in Dry Areas (1999). SWIM Paper 7, Colombo, Sri Lanka: International Water Management Institute.
18 Amarasinghe U.A. Mutuwatta L. Sakthivadivel R. Water Scarcity Variations within a Country: A Case Study of Sri Lanka 1999 International Water Management Institute Colombo, Sri Lanka Research Report 32
19 Ozcelik N. Rodríguez M. Lutter S. Sartal A. Indicating the wrong track? A critical appraisal of water productivity as an indicator to inform water efficiency policies Resour. Conserv. Recycl. 168 2021 105452 10.1016/j.resconrec.2021.105452
20 Steduto P. Hsiao T.C. Fereres E. On the conservative behavior of biomass water productivity Irrigat. Sci. 25 2007 189 207 10.1007/s00271-007-0064-1
21 Molden D. Murray-Rust H. Sakthivadivel R. Makin I. A water-productivity framework for understanding and action Water productivity in agriculture: Limits opportunities for improvement 1 2003 1 18
22 Nouri M. Drought assessment using gridded data sources in data-poor areas with different aridity conditions Water Resour. Manag. 2023 10.1007/s11269-023-03555-4
23 Veysi S. Nouri M. Jabbari A. Reference evapotranspiration estimation using reanalysis and WaPOR products in dryland Croplands Heliyon 10 2024 e26531 10.1016/j.heliyon.2024.e26531
24 Adamala S. An overview of big data applications in water resources engineering Mach. Learn. Res 2 2017 10 18
25 Veysi S. Naseri A.A. Hamzeh S. Relationship between field measurement of soil moisture in the effective depth of sugarcane root zone and extracted indices from spectral reflectance of optical/thermal bands of multispectral satellite images Journal of the Indian Society of Remote Sensing 48 2020 1035 1044 10.1007/s12524-020-01135-2
26 Barideh R. Veysi S. Ebrahimipak N. Davatgar N. The challenge of reference evapotranspiration between the WaPOR data set and geostatistical methods Irrigat. Drain. 71 2022 1268 1279 10.1002/ird.2738
27 Wulder M.A. White J.C. Goward S.N. Masek J.G. Irons J.R. Herold M. Cohen W.B. Loveland T.R. Woodcock C.E. Landsat continuity: issues and opportunities for land cover monitoring Rem. Sens. Environ. 112 2008 955 969 10.1016/j.rse.2007.07.004
28 Manfreda S. McCabe M.F. Miller P.E. Lucas R. Pajuelo Madrigal V. Mallinis G. Ben Dor E. Helman D. Estes L. Ciraolo G. Müllerová J. Tauro F. De Lima M.I. De Lima J.L.M.P. Maltese A. Frances F. Caylor K. Kohv M. Perks M. Ruiz-Pérez G. Su Z. Vico G. Toth B. On the use of unmanned aerial systems for environmental monitoring Rem. Sens. 10 2018 2018 10.3390/rs10040641
29 Weiss M. Jacob F. Duveiller G. Remote sensing for agricultural applications: a meta-review Rem. Sens. Environ. 236 2020 111402 10.1016/j.rse.2019.111402
30 Bastiaanssen W.G.M. Steduto P. The water productivity score (WPS) at global and regional level: methodology and first results from remote sensing measurements of wheat, rice and maize Sci. Total Environ. 575 2017 595 611 10.1016/j.scitotenv.2016.09.032 27712867
31 Leff B. Ramankutty N. Foley J.A. Geographic distribution of major crops across the world Global Biogeochem. Cycles 18 2004 10.1029/2003gb002108
32 Portmann F.T. Siebert S. Döll P. MIRCA2000—global monthly irrigated and rainfed crop areas around the year 2000: a new high‐resolution data set for agricultural and hydrological modeling Global Biogeochem. Cycles 24 2010 10.1029/2008gb003435
33 Pelgrum H. Miltenburg I. Cheema M. Klaasse A. Bastiaanssen W. ETLook a novel continental evapotranspiration algorithm. Remote Sensing and Hydrology Proceedings of a Symposium Held at Jackson Hole. Held in Wyoming, USA, September 2010 2010 IAHS Publ. 3XX, 2011 1087 2010
34 Bastiaanssen W.G.M. Cheema M.J.M. Immerzeel W.W. Miltenburg I.J. 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 2012 10.1029/2011WR010482
35 Mavroforakis M.E. Theodoridis S. A geometric approach to Support Vector Machine (SVM) classification IEEE Trans. Neural Network. 17 2006 671 682 10.1109/TNN.2006.873281
36 Peña-Barragán J.M. Ngugi M.K. Plant R.E. Six J. Object-based crop identification using multiple vegetation indices, textural features and crop phenology Rem. Sens. Environ. 115 2011 1301 1316 10.1016/j.rse.2011.01.009
37 Nitze I. Schulthess U. Asche H. Comparison of machine learning algorithms random forest, artificial neural network and support vector machine to maximum likelihood for supervised crop type classification Proceedings of the 4th GEOBIA vol. 79 2012 3540 Rio de Janeiro, Brazil
38 Belgiu M. Drăguţ L. Random forest in remote sensing: a review of applications and future directions ISPRS J. Photogrammetry Remote Sens. 114 2016 24 31 10.1016/j.isprsjprs.2016.01.011
39 Congalton R.G. A review of assessing the accuracy of classifications of remotely sensed data Rem. Sens. Environ. 37 1991 35 46 10.1016/0034-4257(91)90048-B
40 Safi A.R. Karimi P. Mul M. Chukalla A. de Fraiture C. Translating open-source remote sensing data to crop water productivity improvement actions Agric. Water Manag. 261 2022 107373 10.1016/j.agwat.2021.107373
41 Filippi P. Whelan B.M. Vervoort R.W. Bishop T.F.A. Identifying crop yield gaps with site- and season-specific data-driven models of yield potential Precis. Agric. 23 2022 578 601 10.1007/s11119-021-09850-7
42 Alauddin M. Amarasinghe U.A. Sharma B.R. Are there any ‘bright’ spots and ‘hot’ spots of rice water productivity in Bangladesh? A spatio-temporal analysis of district-level data Oxford Business Economics and Conference, St Hugh's College 2010 Oxford University Oxford UK 28 Oxford, U.K
43 Moesinger L. Zotta R.M. van der Schalie R. Scanlon T. de Jeu R. Dorigo W. Monitoring vegetation condition using microwave remote sensing: the standardized vegetation optical depth index (SVODI) Biogeosciences 19 2022 5107 5123 10.5194/bg-19-5107-2022
44 van Genuchten M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils Soil Sci. Soc. Am. J. 44 1980 892 898 10.2136/sssaj1980.03615995004400050002x
45 Jin N. Tao B. Ren W. Feng M. Sun R. He L. Zhuang W. Yu Q. Mapping irrigated and rainfed wheat areas using multi-temporal satellite data Rem. Sens. 8 2016 207
46 Mesgaran B. Azadi P. A national adaptation plan for water scarcity in Iran Working Paper 6, Stanford Iran 2040 Project 2018 Stanford University
47 Alzueta I. Abeledo L.G. Mignone C.M. Miralles D.J. Differences between wheat and barley in leaf and tillering coordination under contrasting nitrogen and sulfur conditions Eur. J. Agron. 41 2012 92 102
48 Monserud R.A. Leemans R. Comparing global vegetation maps with the Kappa statistic Ecol. Model. 62 1992 275 293
49 Falkenmark M. Lundqvist J. Towards water security: political determination and human adaptation crucial Nat. Resour. Forum 22 1998 37 51 10.1111/j.1477-8947.1998.tb00708.x
50 Rockström J. Water for food and nature in drought–prone tropics: vapour shift in rain–fed agriculture Phil. Trans. Roy. Soc. Lond. B Biol. Sci. 358 2003 1997 2009 14728794
51 Edreira J.I.R. Mourtzinis S. Conley S.P. Roth A.C. Ciampitti I.A. Licht M.A. Kandel H. Kyveryga P.M. Lindsey L.E. Mueller D.S. Assessing causes of yield gaps in agricultural areas with diversity in climate and soils Agric. For. Meteorol. 247 2017 170 180
52 Tokatlidis I. Koutroubas S. A review of maize hybrids' dependence on high plant populations and its implications for crop yield stability Field Crops Res. 88 2004 103 114
53 Ferreira M.I. Silvestre J. Conceição N. Malheiro A.C. Crop and stress coefficients in rainfed and deficit irrigation vineyards using sap flow techniques Irrigat. Sci. 30 2012 433 447
54 Valentín F. Sánchez J.M. Martínez-Moreno A. Intrigliolo D.S. Buesa I. López-Urrea R. Using on-the-ground surface energy balance to monitor vine water status and evapotranspiration under deficit irrigation and rainfed conditions Agric. Water Manag. 281 2023 108240
55 Kato T. Kimura R. Kamichika M. Estimation of evapotranspiration, transpiration ratio and water-use efficiency from a sparse canopy using a compartment model Agric. Water Manag. 65 2004 173 191
56 Shveytser V. Stoy P.C. Butterworth B. Wiesner S. Skaggs T.H. Murphy B. Wutzler T. El‐Madany T.S. Desai A.R. Evaporation and transpiration from multiple proximal forests and wetlands Water Resour. Res. 60 2024 e2022WR033757
57 Gu C. Tang Q. Zhu G. Ma J. Gu C. Zhang K. Sun S. Yu Q. Niu S. Discrepant responses between evapotranspiration-and transpiration-based ecosystem water use efficiency to interannual precipitation fluctuations Agric. For. Meteorol. 303 2021 108385
58 Pei F. Li X. Liu X. Lao C. Assessing the impacts of droughts on net primary productivity in China J. Environ. Manag. 114 2013 362 371
59 Peña-Gallardo M. Vicente-Serrano S.M. Domínguez-Castro F. Beguería S. The impact of drought on the productivity of two rainfed crops in Spain Nat. Hazards Earth Syst. Sci. 19 2019 1215 1234
60 Haro-Monteagudo D. Daccache A. Knox J. Exploring the utility of drought indicators to assess climate risks to agricultural productivity in a humid climate Nord. Hydrol 49 2018 539 551
61 Wang H. He B. Zhang Y. Huang L. Chen Z. Liu J. Response of ecosystem productivity to dry/wet conditions indicated by different drought indices Sci. Total Environ. 612 2018 347 357 28858745
62 Liu J. Williams J.R. Zehnder A.J. Yang H. GEPIC–modelling wheat yield and crop water productivity with high resolution on a global scale Agric. Syst. 94 2007 478 493
63 Molden D. Oweis T.Y. Pasquale S. Kijne J.W. Hanjra M.A. Bindraban P.S. Bouman B.A.M. Cook S. Erenstein O. Farahani H. Hachum A. Hoogeveen J. Mahoo H. Nangia V. Peden D. Sikka A. Silva P. Turral H. Upadhyaya A. Zwart S. Pathways for increasing agricultural water productivity Molden D. Water for Food, Water for Life: a Comprehensive Assessment of Water Management in Agriculture 310 2007 279
64 Zwart S.J. Bastiaanssen W.G.M. de Fraiture C. Molden D.J. A global benchmark map of water productivity for rainfed and irrigated wheat Agric. Water Manag. 97 2010 1617 1627 10.1016/j.agwat.2010.05.018
65 Pereira L.S. Oweis T. Zairi A. Irrigation management under water scarcity Agric. Water Manag. 57 2002 175 206
66 Wei Z. Paredes P. Liu Y. Chi W.W. Pereira L.S. Modelling transpiration, soil evaporation and yield prediction of soybean in North China Plain Agric. Water Manag. 147 2015 43 53
67 Araya A. Gowda P. Golden B. Foster A. Aguilar J. Currie R. Ciampitti I. Prasad P. Economic value and water productivity of major irrigated crops in the Ogallala aquifer region Agric. Water Manag. 214 2019 55 63
68 Nouri H. Stokvis B. Borujeni S.C. Galindo A. Brugnach M. Blatchford M. Alaghmand S. Hoekstra A. Reduce blue water scarcity and increase nutritional and economic water productivity through changing the cropping pattern in a catchment J. Hydrol. 588 2020 125086
69 Bastiaanssen W.G. Steduto P. The water productivity score (WPS) at global and regional level: methodology and first results from remote sensing measurements of wheat, rice and maize Sci. Total Environ. 575 2017 595 611 27712867
70 Bouman B. Tuong T.P. Field water management to save water and increase its productivity in irrigated lowland rice Agric. Water Manag. 49 2001 11 30
71 Passioura J. Increasing crop productivity when water is scarce—from breeding to field management Agric. Water Manag. 80 2006 176 196
72 Deng X.-P. Shan L. Zhang H. Turner N.C. Improving agricultural water use efficiency in arid and semiarid areas of China Agric. Water Manag. 80 2006 23 40 10.1016/j.agwat.2005.07.021
