Deficit irrigation is a management strategy to improve crop water productivity, especially in arid and semi-arid regions. Soil characteristics and weather parameters are among the factors affecting crop water productivity in water stress conditions. Due to spatial changes in soil characteristics and temporal and spatial variations in meteorological parameters, it can be expected that crop water productivity will also have temporal and spatial variations. In this study, by combining the Geographic Information System (GIS) with the grid weather generation tools from the Crop Growth Monitoring System (CGMS) and the plug-in version of the AquaCrop, a combined method was developed to investigate the spatial and temporal variation of crop yield, seasonal crop evapotranspiration, and water productivity of maize under various irrigation scenarios. The proposed model was implemented in a case study in the west of Iran. The study area was divided into 37 grid weather with 5 * 5 km and 19 soil units. By overlaying soil units and grid weathers, 94 homogeneous units were created. The model was executed for 94 homogeneous areas, using calibrated crop file of grain maize under four irrigation scenarios of 40, 60, 80, and 100% of potential irrigation requirement (S40, S60, S80, and S100, respectively) for 28 years (1988–2015) of weather data (10,528 runs). The results showed that by increasing water stress, the percentage of spatial and temporal variation of the studied parameters (crop yield, seasonal crop water requirement, and water productivity) would be increased. The percentage of spatial changes in crop yield and crop water productivity was more significant than temporal changes. The average of crop water productivity in the scenarios of S100, S80, S60, and S40 was determined as 1.5, 1.4, 1.2, and 0.5 kg m−3, respectively.
