Integrating Remote Sensing Techniques and Meteorological Data to Assess the Ideal Irrigation System Performance Scenarios for Improving Crop Productivity

: To increase agricultural productivity and ensure food security, it is important to understand the reasons for variations in irrigation over time. However, researchers often avoid investigating water productivity due to data availability challenges. This study aimed to assess the performance of the irrigation system for winter wheat crops using a high-resolution satellite, Sentinel 2 A/B, combined with meteorological data and Google Earth Engine (GEE)-based remote sensing techniques. The study area is located north of Erbil city in the Kurdistan region of Iraq (KRI) and consists of 143 farmer-owned center pivots. This study also aimed to analyze the spatiotemporal variation of key variables (Normalized Difference Moisture Index (NDMI), Normalized Difference Vegetation Index (NDVI), Precipitation (mm), Evapotranspiration (ETo), Crop evapotranspiration (ETc), and Irrigation (Hours), during the wheat-growing winter season in the drought year 2021 to understand the reasons for the variance in ﬁeld performance. The ﬁnding revealed that water usage ﬂuctuated signiﬁcantly across the seasons, while yield gradually increased from the 2021 winter season. In addition, the study revealed a notable correlation between soil moisture based on the (NDMI) and vegetation cover based on the (NDVI), and the increase in yield productivity and reduction in the yield gap, speciﬁcally during the middle of the growing season (March and April). Integrating remote sensing with meteorological data in supplementary irrigation systems can improve agriculture and water resource management by boosting yields, improving crop quality, decreasing water consumption, and minimizing environmental impacts. This innovative technique can potentially enhance food security and promote environmental sustainability.


Introduction
Water stress is a pressing issue that affects many regions of the world, particularly in arid and semi-arid areas where water resources are scarce.In recent years, governments and humanitarian organizations have attempted to address this issue by improving access to water for those living in water-stressed areas.However, as climate change and population growth continue to pressure water resources, the problem is expected to worsen in the coming years.In addition, water stress can impact agricultural production, which can have an effect on the economy and food security [1,2].Water stress can result in conflicts and mass migrations in certain cases.Moreover, the migration of people due to water stress can pose new challenges for governments and humanitarian organizations.Mitigating the challenges associated with water stress necessitates a collective endeavor from all stakeholders.Global collaboration in water management, in conjunction with efficacious policies and practices at the regional level, can facilitate the provision of secure and hygienic water access for all, and avert the societal, economic, and ecological repercussions stemming from water stress [3].Countries are being encouraged to implement innovative and sustainable measures such as rainwater harvesting and enhanced irrigation techniques [2].The Middle East has faced a severe water scarcity problem in recent decades due to climate change and inadequate management of water resources [4,5].Water scarcity considerably limits wheat production in Northern Iraq, especially in regions with low and erratic rainfall during the crop growth period.In such circumstances, supplemental irrigation can boost and sustain yields over time, particularly if applied during critical growth stages [6].
Supplementary irrigation has been identified as a plausible measure to address this quandary.This entails the provision of supplementary water by farmers to their crops to offset the shortfall in precipitation, especially during crucial growth phases, since wheat necessitates varying quantities of water at distinct stages of its growth.The application of supplementary irrigation at vital growth stages warrants the provision of adequate water supply to the crop during its peak water requirements, and thereby augments and preserves yields over a protracted time period, even in regions confronted with a dearth of or erratic precipitation [1,7].
The inaccessibility of sufficient water resources can pose a formidable obstacle to farmers, particularly during periods of drought [8].Shallow and deep tube wells are both employed as means to provide irrigation water.In 2007, the Ministry of Agriculture and Water Resources of the KRI initiated a program to assist farmers with contemporary irrigation systems, such as the center pivot, to enhance crop production after experiencing droughts on several occasions in previous years.Nevertheless, inadequate awareness regarding water scarcity and efficient water utilization has led some farmers to utilize the center pivot for prolonged durations and over-irrigate by extending the center pivot revolution time to amplify the water application depth.The assimilation of crop models and remotely sensed data via optimization algorithms has emerged as an efficacious and prospective technique for monitoring crop growth status and appraising crop yields, as it ameliorates specific deficiencies and amalgamates the benefits of individual methods [9,10].
In the KRI region, groundwater is abundant and generally of reasonable quality, but it has been heavily extracted for domestic, agricultural, and industrial purposes, resulting in a depletion of the groundwater table level.Mitigating the effects of climate change and rehabilitating nature requires more water management projects to ensure the efficient use of available water.Groundwater is the primary source of inaccessible surface water, and rainfall, groundwater, and rivers account for 51%, 48%, and 1% of crop production water, respectively [11].Despite an increasing number of wells and groundwater abstraction in the Erbil plain, the aquifer is being harmed, with an average drop in groundwater levels of 50 m [2].
Agriculture is the largest water-consuming sector globally, accounting for 78-90% of water use, followed by domestic and industrial use [12].In Iraq, the agricultural sector is the primary consumer of water resources.Irrigation in Iraq relies on three main sources of water: surface water, rainwater, and groundwater, and farmers use drip, sprinkler, and central pivot irrigation techniques [13].Wheat and barley are the dominant crops in Iraq, covering 73% of the total cultivated land [14].Wheat productivity under irrigation exceeds rainfed productivity by a factor of two to three, demonstrating the high productivity potential of wheat varieties under irrigation [2].Given the increasing scarcity of freshwater, it is critical to optimize water use, particularly in irrigated agriculture [15].Rainfall is the primary source of water for agriculture in Iraq, accounting for 51% of water supply during winter, while groundwater and rivers contribute 48% and 1%, respectively [16], as per the joint report of UNDP, USAID, and FAO (2019).Farmers in arid regions are confronted with a multitude of challenges, such as low agricultural productivity, frequent droughts, climate Water 2023, 15, 1605 3 of 32 variability, high soil erosion rates, and deforestation [17].Nevertheless, a comprehensive study of drought must be conducted for each area to develop a research framework [18,19].Understanding the spatial and temporal variations in crop growth is crucial for effective crop management and achieving food security.The combination of remote sensing data and crop growth models has proven to be a valuable approach for monitoring crop growth and estimating crop yields.
The impact of climate change, combined with reduced water discharge in the Tigris and Euphrates rivers, has led to more frequent droughts in Iraq, which have particularly affected agriculture, resulting in significant reductions in crop yields and vegetation cover [20].The Erbil province, KRI, is particularly vulnerable to severe drought, with reports indicating a reduction in precipitation and vegetation cover in 2000, 2008, 2012, and 2021 [21,22].Southern and western parts of the KRI have been the hardest hit by rainfall shortages, whereas the north and east have seen increased moisture levels.The FAO has predicted that wheat crop production in northern Iraq will be approximately 50% lower in 2021 than the previous year, according to KRI authorities [8,23].The use of secondary agricultural data in combination with remote sensing can provide more accurate estimates of irrigation performance [24].The NDVI is a useful parameter for crop monitoring due to its unambiguous response to irrigation interventions.Modern irrigation management strategies aim to enhance agricultural productivity while minimizing water consumption [25], which is crucial given the limited water supply in changing environmental conditions.The NDMI can be used to describe a crop's water stress level and is calculated as the ratio between the difference and the sum of the refracted radiation in the near-infrared and SWIR spectrums [26].Effective water resource management and enhancement of water productivity require a sound understanding and management of evapotranspiration [27].Monitoring significantly varying Normalized Difference Vegetation Index (NDVI) values in a region can provide valuable insight into potential pest or disease concerns and the presence of weeds in areas of a field [28].NDVI data can be used to display changes in vegetation cover and analyze patterns of drought occurrences.However, the performance of NDVI may be influenced by mistakes during the growing season and saturation effects on dense vegetation [29].Therefore, additional factors should be considered to improve the precision of the results [30].Various alternative vegetation and moisture indices are available to monitor crop growth, including the Enhanced Vegetation Index (EVI), Normalized Difference Water Index (NDWI), Leaf Area Index (LAI), Soil Adjusted Vegetation Index (SAVI), Thermal Vegetation Index (TVI), Crop Water Stress Index (CWSI), and Temperature Vegetation Dryness Index (TVDI).The most appropriate index selection depends on various factors, such as the type of crop, growth stage, and environmental conditions.NDVI have been widely used in semi-arid areas and have shown good performance in detecting vegetation changes [31][32][33].The Normalized Difference Moisture Index (NDMI) is a reliable indicator of vegetation water content and can be used to monitor water stress in crops or other vegetation, particularly in areas where water stress is a concern [34,35].NDMI is derived from the difference between the near-infrared and mid-infrared bands of satellite imagery.By monitoring NDMI values over time, changes in vegetation water content can be detected and the level of water stress experienced by crops or other vegetation can be assessed [15].
Inefficient water application occurs when the irrigation system is not appropriately managed, resulting in uneven water distribution across the field, reduced crop yields, and soil erosion.Such inefficiency happens when water evaporates or is lost through other means.Minimizing water loss requires proper irrigation scheduling, use of efficient irrigation systems, and managing the system appropriately based on remote sensing.Monitoring through the Sentinel satellite every five days and ETo, ETc data helps ensure proper irrigation scheduling based on remote sensing and GIS techniques.The capability to develop agriculture effectively with limited water resources is a crucial strategic objective for addressing future climate change and fulfilling Sustainable Development Goal 2 of the United Nations (SDG2).The primary objective of this research is to comprehend the impact of drought and compare wheat production between rainfed and center pivot systems during the study period.Additionally, it seeks to analyze the key factors that affect the spatial and temporal variation of NDVI and NDWI in wheat cultivated on center pivots.The study aims to investigate irrigation performance and water productivity at various scales to develop suitable water management strategies, especially considering decreasing water availability, rising threats from climate change, and growing population and food demand.This study provides a new strategy for agricultural resource management by providing consistent estimations of winter wheat water requirements and yield.This information can be used to optimize irrigation practices and improve crop management, particularly in arid and semi-arid regions where water availability is limited.

Study Area
The study area, comprising the Ain Kawa sub-district in the Erbil Governorate of northern Iraq, is illustrated in Figure 1 and spans an estimated area of 25,525 hectares.North Erbil was chosen as the case study region for a combination of reasons.Firstly, the wheat, which is a crucial crop that contributes significantly to the Kurdistan Region of Iraq's economy, is extensively grown in the region.Secondly, Erbil is a vital agricultural center, known as Iraq's breadbasket, and is responsible for a substantial portion of the country's food supply.Finally, the need to enhance and expand the existing water management sector in Iraq is pressing, and there are related challenges that require immediate attention.
irrigation scheduling based on remote sensing and GIS techniques.The capability to develop agriculture effectively with limited water resources is a crucial strategic objective for addressing future climate change and fulfilling Sustainable Development Goal 2 of the United Nations (SDG2).The primary objective of this research is to comprehend the impact of drought and compare wheat production between rainfed and center pivot systems during the study period.Additionally, it seeks to analyze the key factors that affect the spatial and temporal variation of NDVI and NDWI in wheat cultivated on center pivots.The study aims to investigate irrigation performance and water productivity at various scales to develop suitable water management strategies, especially considering decreasing water availability, rising threats from climate change, and growing population and food demand.This study provides a new strategy for agricultural resource management by providing consistent estimations of winter wheat water requirements and yield.This information can be used to optimize irrigation practices and improve crop management, particularly in arid and semi-arid regions where water availability is limited.

Study Area
The study area, comprising the Ain Kawa sub-district in the Erbil Governorate of northern Iraq, is illustrated in Figure 1 and spans an estimated area of 25,525 hectares.North Erbil was chosen as the case study region for a combination of reasons.Firstly, the wheat, which is a crucial crop that contributes significantly to the Kurdistan Region of Iraq's economy, is extensively grown in the region.Secondly, Erbil is a vital agricultural center, known as Iraq's breadbasket, and is responsible for a substantial portion of the country's food supply.Finally, the need to enhance and expand the existing water management sector in Iraq is pressing, and there are related challenges that require immediate attention.

Climate
The climatic conditions prevailing in Erbil, KRI, are predominantly continental, subtropical, and semi-arid, while a Mediterranean climate characterizes the mountainous parts.Rainfall in the mountainous areas usually occurs between December to February or Novem-  C at night [36,37].The period of precipitation in the region is typically between October to May, with the annual rainfall ranging between 100 and 200 mm.The Ministry of Agriculture and the Department of Meteorology in the KRI have validated the accuracy of the average precipitation recorded during this period.It is worth noting that farmers in the study region rely on rainfall to irrigate their crops, although they often supplement it with irrigation to enhance crop productivity.

Meteorological Data
Figure 2 presents the data on the mean daily maximum and minimum temperature, relative humidity, wind speed, sunshine, and hours, collected from the Ministry of Agriculture and Water Resources of Erbil for two stations during the year 2021.

Climate
The climatic conditions prevailing in Erbil, KRI, are predominantly continental, subtropical, and semi-arid, while a Mediterranean climate characterizes the mountainous parts.Rainfall in the mountainous areas usually occurs between December to February or November to April.The average daytime temperature during winter is approximately 16 degrees Celsius, which drops to around 2 °C at night, with a possibility of frost.Conversely, summers are characterized by high temperatures, with an average temperature exceeding 45 °C in July and August and falling to around 25 °C at night [36,37].The period of precipitation in the region is typically between October to May, with the annual rainfall ranging between 100 and 200 mm.The Ministry of Agriculture and the Department of Meteorology in the KRI have validated the accuracy of the average precipitation recorded during this period.It is worth noting that farmers in the study region rely on rainfall to irrigate their crops, although they often supplement it with irrigation to enhance crop productivity.

Meteorological Data
Figure 2 presents the data on the mean daily maximum and minimum temperature, relative humidity, wind speed, sunshine, and hours, collected from the Ministry of Agriculture and Water Resources of Erbil for two stations during the year 2021.

Sentinel 2 Satellite Imagery Acquisition
This study used open-access satellite data from the Sentinel 2 missions to calculate the vegetation and water indices as a proxy for factors.Sentinel-2 is a 10 m high-resolution, wide-swath, multispectral imaging mission that supports Copernicus Land Monitoring investigations, such as plant, soil, and water cover monitoring, and observation of both inland waterways and coastal regions [38].Thirteen spatial resolution spectral bands captured by Sentinel-2 represent TOA reflectance, including four bands at 10 m, six at 20 m,

Sentinel 2 Satellite Imagery Acquisition
This study used open-access satellite data from the Sentinel 2 missions to calculate the vegetation and water indices as a proxy for factors.Sentinel-2 is a 10 m high-resolution, wide-swath, multispectral imaging mission that supports Copernicus Land Monitoring investigations, such as plant, soil, and water cover monitoring, and observation of both inland waterways and coastal regions [38].Thirteen spatial resolution spectral bands captured by Sentinel-2 represent TOA reflectance, including four bands at 10 m, six at 20 m, and three at 60 m spatial resolution.In addition to data from the Sentinel 2 missions, detailed meteorological and secondary data from the literature were required, including corn winter wheat crops for monthly, seasonal periods in 2021, with season intervals ranging from January to May.The Copernicus Open Access Hub was used for this investigation to obtain Level-1C products of Sentinel-2 A/B satellite data.These products contain top-ofatmosphere (TOA) reflectance values for 13 spectral bands with spatial resolutions ranging from 10 to 20 m.The Sentinelsat Python API was used to obtain the data for the study area and time period of interest, which were subsequently preprocessed in Google Earth Engine (GEE) to acquire surface reflectance values.The preprocessing procedure involved resampling the data to a standardized spatial resolution of 10 m, masking out clouds and cloud shadows utilizing the SCL band, and applying a bi-directional reflectance distribution function (BRDF) correction to consider the impact of directional effects induced by surface roughness and slope.
In order to enable more in-depth analysis and visualization, the classified Sentinel 2 data underwent a process of division into smaller tiles.This was accomplished utilizing the ee.data.getTileUrlfunction in the GEE cloud-based platform, which allowed the extraction of individual tiles based on their geographical coordinates and zoom level.The resulting tiles were then saved in the GeoTIFF format, which can be easily imported into QGIS and other GIS software for further processing and analysis.The classified and divided Sentinel 2 data were subsequently subjected to a range of statistical analyses, including frequency distributions, cross-tabulations, and spatial autocorrelation analysis, which were performed using both GEE and QGIS software.These analyses were employed to quantify the spatial patterns and relationships between different meteorological variables in the study area.Additionally, it is imperative to comprehend the current status of the center pivot's performance for the purpose of determining its potential for improvement.As such, a manual selection process was undertaken, which involved identifying 143 center pivot farms in the vicinity for further investigation.
The biophysical factors were computed using the formula derived from Sentinel 2 Indices, which comprises two satellite sensors (S2A and S2B) that have generated products with a 5-day temporal resolution (at the equator) and a 10 m spatial grid-cell resolution since January 2021.The satellite products from the Sentinel 2 mission were obtained and analyzed using GEE.The QGIS (v 24.3) geographic information system program utilized the Semi-Automatic Classification Tool plugin, with the following options defined in the plugin's menu: (a) Band 4-RED and Band 8-NIR from Sentinel 2; (b) the study area coordinates; (c) the time search windows from 1 January to 21 May in 2021, as well as (Figure 3); and (d) an acceptable imagery cloud cover set to 100%.Figure 1 illustrates a Google Earth satellite image of a farm utilizing a center-pivot irrigation system in the study area.

Validation Fields
The irrigation hours and actual evapotranspiration estimates were validated using available measurements from 14 wheat fields (C1, C2, C3, …, C14) from January to May (Figure 4).Comparison has been made the weekly evolution of irrigation, observed actual ET, and the amount of irrigation and rainfall for the 14 wheat fields.A general coherence was observed, suggesting the good performance of irrigation hours in mitigating crop water stress through supplementary irrigation.The center pivot of the study area, character-

Validation Fields
The irrigation hours and actual evapotranspiration estimates were validated using available measurements from 14 wheat fields (C1, C2, C3, . . ., C14) from January to May (Figure 4).Comparison has been made the weekly evolution of irrigation, observed actual ET, and the amount of irrigation and rainfall for the 14 wheat fields.A general coherence was observed, suggesting the good performance of irrigation hours in mitigating crop water stress through supplementary irrigation.The center pivot of the study area, characterized by higher production and a regular irrigation system, had higher NDVI-based vegetation density in fields C1 to C14, while fields in other center pivots had lower density.

Validation Fields
The irrigation hours and actual evapotranspiration estimates were validated using available measurements from 14 wheat fields (C1, C2, C3, …, C14) from January to May (Figure 4).Comparison has been made the weekly evolution of irrigation, observed actual ET, and the amount of irrigation and rainfall for the 14 wheat fields.A general coherence was observed, suggesting the good performance of irrigation hours in mitigating crop water stress through supplementary irrigation.The center pivot of the study area, characterized by higher production and a regular irrigation system, had higher NDVI-based vegetation density in fields C1 to C14, while fields in other center pivots had lower density.

Spectral Indices 2.5.1. NDVI
The Near-Infrared (NIR) and Red bands were also applied to calculate the NDVI images [39].The NDVI index is calculated with the aid of the red (Red) and the Near-Infrared (NIR) bands of Sentinel-2 images, using Equation ( 1), as follows [39]: where B4 is RED, 664.5 nm, and B8 is NIR, 835.1 nm.
Theoretically, NDVI values ranged between −1.0 and +1.0.However, the typical range of NDVI gauged from vegetation and other earth surface materials is between approximately −0.1 (NIR less than Red) for non-vegetated surfaces and as high as 0.9 for dense vegetative cover.To better understand NDVI values, they are often classified into different ranges.Typically, values between 0.1 and 0.4 indicate low vegetation density or sparse vegetation, while values between 0.4 and 0.6 indicate moderate density or healthy vegetation.Values above 0.6 indicate high density or very healthy vegetation.These NDVI density classes can be beneficial when monitoring vegetation growth over time or comparing different areas.For instance, if an area consistently has NDVI values within the low-density range, it may indicate the need for irrigation or other interventions to promote plant growth.Conversely, if an area consistently exhibits high NDVI values, it may suggest a healthy and productive ecosystem [40].

NDMI
The NDMI normalizes the different moisture response bands between near-infrared (NIR) and shortwave infrared (SWIR) (Equation ( 2)).The linear correlation between the NIR/SWIR ratio and leaf relative water content was discovered by Hunt Jr and Rock [41].He calculated NDMI using the following equation: where NIR is the Near-Infrared band and SWIR is the Shortwave Infrared band.The values of these bands can be obtained from remote sensing data, such as satellite imagery.
The NDMI ranges from −1 to 1, where values closer to 1 indicate high moisture content in vegetation, and values closer to −1 indicate low moisture content.In general, vegetation with high moisture content reflects more NIR and absorbs more SWIR, resulting in a higher NDMI value.NDMI is often used in agriculture and environmental monitoring to assess vegetation health and drought conditions.It can be used to detect areas of vegetation stress, monitor changes in soil moisture, and predict crop yields [41].

Meteorological Indices Reference Evapotranspiration (ETo), and Crop Evapotranspiration (ETc)
The FAO Penman-Monteith equation is a close, simple representation of the physical and physiological factors governing the evapotranspiration process.Using the FAO Penman-Monteith definition for (ETo), the formula for ETo using the Penman-Monteith equation is: where: ETo = potential evapotranspiration (mm/day).delta = slope of the saturation vapor pressure-temperature curve (kPa/ • C).Rn = net radiation at the crop surface (MJ/m2/day).T = mean daily air temperature at 2 m height ( • C).U2 = wind speed at 2 m height (m/s).es = saturation vapor pressure (kPa).ea = actual vapor pressure (kPa).gamma = psychrometric constant (kPa/ • C) [42].Crop evapotranspiration (ETc) is the amount of water lost from a crop due to evaporation from the soil surface and transpiration from the crop itself.It is a measure of the amount of water required for a specific crop to achieve optimum growth and yield.Various factors, including the climate, soil characteristics, crop type, and stage of growth, influence ETc.The most common method used to estimate ETc is using reference evapotranspiration (ETo), which is the amount of water lost from a reference crop under standardized conditions.Once the ETo is calculated, crop coefficients are applied to adjust the ETo for the specific crop being grown [42].
The resulting value is the crop evapotranspiration (ETc) for that crop at that location and time.ETc is an important parameter in irrigation management, as it determines the amount of water that must be applied to a crop to maintain optimal growth and yield.Overirrigation can lead to waterlogging and the leaching of nutrients, while under-irrigation can reduce crop yield and quality.Therefore, an accurate estimation of ETc is crucial for the efficient and sustainable use of water resources in agriculture [42].The ETc equation is expressed as: where: ETc = crop evapotranspiration (mm/day).Kc = crop coefficient (dimensionless).

NDVI and NDMI
Based on remote sensing data, the NDMI values show more pronounced changes between February and March, potentially due to local farmers perceiving irrigation as unnecessary during January and February's cold weather.The higher variability in NDMI values indicates that it responds to changes in soil moisture content.Similarly, the greater vertical variability in biomass suggests that wheat biomass reacts more significantly to moisture stress during the initial growth stage (February).This implies that crops are more vulnerable to stress in March and April than in other months.A comparison of NDVI-and NDMI-designated fields indicates that irrigation's impact on crop yield varies geographically, with a substantial yield increase observed in areas with abundant light.To maximize yield, it is crucial to avoid water stress during the wheat's most susceptible seasons.However, there was a moderately negative correlation between NDVI and NDMI in 2021.Moreover, the majority of hotspot fields exhibited higher NDMI values, while bright spots showed early growth and greater biomass accumulation throughout the growing stages.This difference may be related to irrigation methods and soil moisture stress, particularly in March and April when the center pivot is essential for reducing temperature and relieving soil moisture stress.However, some farmers do not irrigate during January and February, which may lead to water stress.Remote sensing can detect the effects of water stress on plants with a lag time.
The NDVI has been widely used to investigate the correlation between spectral vegetation variability and growth rate changes.The findings of this investigation have demonstrated that the NDVI values fluctuated between 174.5 area/ha in January, the lowest area of dense NDVI, and 6028.3 ha in April, the highest area.The NDVI Value is a numerical index employed to assess vegetation health and growth.It indicates denser and healthier vegetation, with higher values ranging from −1 to +1.The NDVI value categorizes vegetation density into different classes such as Dense, Moderate, Sparse, and Open vegetation.The area of each category is measured in hectares (ha), a common unit in agriculture.Table 1 shows the NDVI Area/ha for Dense Vegetation Classes for different months, which include 1-Jan, 6-Jan, 11-Jan, 26-Jan, 15-Feb, 25 According to the data presented in Table 1, there was a noticeable change in the vegetation status of the Center pivot area from January to May.Specifically, the NDVI value of the Center pivot area decreased significantly in May due to the maturity stage of the crops and the increase in ETo, which resulted in a reduction in the agricultural land area.Moreover, the rise in temperature was detectable earlier, particularly during the drought season.Thus, NDVI with NDMI can help improve our understanding of how irrigation and climate change affect the yield and can be used as an early warning for Moisture stress.Moreover, this information can assist farmers and policymakers in making more informed management decisions.In order to achieve the greatest range of NDVI and NDMI values for wheat crops in 143 central pivot areas in north Erbil, we utilized Sentinel 2 images on coincident days.The dataset used for NDVI and NDMI was then split into two subsets: a training dataset and a validation dataset.These datasets were collected from 1 January 2021 to 21 May 2021, covering fourteen Center pivot areas throughout the phonological cycle to monitor various crops.The analysis of Tables 1 and 2 and Figures 5 and 6 revealed that the study region experienced drought episodes over the study period, with March and April 2021 being particularly affected.Table 1 shows the correlation between drought and reduced vegetation area/NDVI mean values, further supporting the importance of water availability for vegetation growth.It is also possible that other factors, such as temperature and nutrient availability, may contribute to the observed differences in vegetation indices between the two fields.Table 1 shows the correlation between drought and reduced vegetation area/NDVI mean values, further supporting the importance of water availability for vegetation growth.It is also possible that other factors, such as temperature and nutrient availability, may contribute to the observed differences in vegetation indices between the two fields.There is a clear difference in the vegetation indices (NDVI and NDMI) between the irrigated and rain-fed fields during the growing season months of March and April.The irrigated field shows a significant change in NDVI values and NDMI, which may indicate increased vegetation growth due to water availability.On the other hand, the rain-fed field experiences drought during these months, which is correlated with reduced vegetation area and lower NDVI mean values.It is worth noting that NDVI is a commonly used index for evaluating the health and growth of vegetation, as it quantifies the chlorophyll content in plant leaves.On the other hand, NDMI is an indicator of vegetation water content and can offer insights into soil moisture conditions [7].
Consequently, the observed dissimilarities in NDVI and NDMI values between irrigated and rain-fed fields imply that water availability is a crucial determinant of vegetation growth and productivity.Nonetheless, it was observed that NDMI values increased significantly during the months of March and April, especially in the irrigated field, when the crops were subjected to 72 to 96 h of irrigation per month.In contrast, the NDVI value showed very little vegetation, particularly in the non-irrigated section of the study area (see Figures 5 and 6).
These findings are consistent with those of a previous study by [43].The temporal fluctuations in NDVI and NDMI demonstrated that the drought resulted in a decline in the value and extent of spectral indices from normal to extreme levels, with the non-irrigated area being the driest.Conversely, the field where supplementary irrigation was applied during the specified period had high NDVI coverage.The NDVI values revealed that the non-irrigated area in the study region was affected by drought (see Figures 5 and 6).However, the severity and spatial extent of the drought varied.Low values indicated a dry season, whereas high values indicated a wet season [44].From January to May, the NDVI values exhibited temporal and spatial variations across the region (Figures 5 and 6), primarily attributed to differences in precipitation amount, frequency, and intensity [45].Moreover, meteorological factors such as precipitation, temperature, and relative humidity were crucial in influencing NDVI variability [46].The spatial variability of NDVI values was determined by several factors, including climate, soil, temperature, ETo, and ETc fluctuations (refer to Table 3).The NDMI is a remote sensing metric used to measure how much water is in vegetation.Higher values mean that the vegetation has more water.Such areas may be irrigated or have access to other sources of water.The NDMI values for certain months and areas.Larger values mean that plants are especially healthy because there is a lot of water.Specifically, Both the NDVI and the NDMI are commonly used as remote sensing indices to measure how healthy plants are and how much water they have.This information could be used to improve irrigation efficiency and determine how much water plants need based on satellite data.NDVI is typically used to estimate crop yield, vegetation cover, and photosynthetic activity.In this study, NDVI was used to find plants that were water stressed, which is a common problem in agriculture that uses water.Farmers can improve irrigation schedules, increase crop yields, and waste less water by keeping an eye on NDVI values.
The NDMI index is utilized for the estimation of soil moisture and water stress in crops, and it can aid farmers in adjusting their irrigation schedules to optimize water usage and improve crop yield.Changes in NDVI and NDMI values prior to and following irrigation can assist farmers in gauging the amount of water required to attain the desired vegetation growth or moisture content.The NDMI values were observed to be at their minimum during January and May, indicating reduced vegetation and lower NDVI values (Figures 5 and 6), which could be attributed to decreased rainfall during those months.However, during February, March, and April, NDMI values increased, which was positively correlated with the application of supplementary irrigation, suggesting that irrigation had a beneficial impact on vegetation growth during these months.The decline in NDMI values during the growth season could be attributable to the adverse effects of high temperatures on vegetation growth.The study area witnessed a significant decrease in annual rainfall averages, which could have impacted overall vegetation growth.The 14 fields received 445.2 mm and 579.6 mm of irrigation water, respectively.
During the 2021 monitoring season, irrigation amounts were measured in Fields C1 to C14 using a water meter, which recorded a total of 350 mm of water.One or two irrigation events were performed during the early growing stage to ensure the establishment of robust young plants.Subsequently, irrigation for Fields C1 to C14 was scheduled based on NDMI sensor measurements to irrigate when the depletion level within the rooting zone was ready to drop 60% of the total available below.An algorithm was executed with a crop file created based on Fields measurements and observations and dynamically modified during the cultivation season by adjusting weather parameters.Field-specific constraints, such as not allowing irrigation events to occur more often than once every four days, were integrated into the scheduling process.In cases where the generated irrigation schedule was not followed (e.g., due to an electrical failure in the pumping station or damage to the irrigation system), the algorithm was re-executed, and the irrigation schedule was adjusted accordingly.Despite these adjustments, consistent violations of the first depletion threshold occurred, as irrigation was applied either before the threshold was reached or after the water content had fallen below it.
In dryland areas, supplemental irrigation is often necessary to support crop growth, especially during the winter months from January to May.The region employs various irrigation technologies, including Center pivot irrigation and homogeneous irrigation systems.Monitoring vegetation variability using central pivots is feasible, and the NDMI is a useful input variable for irrigation prediction models, as it correlates with the NDVI, which is considered moderate in Center Pivot techniques (as demonstrated in Figures 6 and 7).It is crucial to validate modeling results to ensure their reliability.Figure 8 shows low variability, indicating that the modeling approaches used in this study were robust.The study utilized a combination of remote sensing techniques and meteorological data to evaluate the efficacy of irrigation systems, and the results indicated that this approach was reliable and robust.Remote sensing was employed to collect data on the water supply to crops and its impact on crop growth, while meteorological data provided information on weather patterns, including precipitation, temperature, humidity, ETo, and ETc, which affected crop growth and irrigation needs.The integration of these two data types enabled the development of a comprehensive model for assessing irrigation system performance.The study found that this approach was effective in providing accurate information about the system's performance, which could help optimize irrigation methods and enhance crop yields.
Notably, all methodologies employed in this study were linear based.Ideally, a well-designed irrigation system would ensure uniform water distribution, which would result in homogeneous water consumption patterns across fields.However, several factors, including deficient infrastructure, inadequate management practices, soil type, water quality, and fertilization, can lead to non-uniform water use across fields and zones.The amount of irrigation water applied in each event was measured through the use of rain gauges positioned above the crop canopy (see Figures 7 and 8).This methodology was proposed and implemented by María [47].
Water 2023, 15, x FOR PEER REVIEW 15 of 34 to crops and its impact on crop growth, while meteorological data provided information on weather patterns, including precipitation, temperature, humidity, ETo, and ETc, which affected crop growth and irrigation needs.The integration of these two data types enabled the development of a comprehensive model for assessing irrigation system performance.
The study found that this approach was effective in providing accurate information about the system's performance, which could help optimize irrigation methods and enhance crop yields.Notably, all methodologies employed in this study were linear based.Ideally, a welldesigned irrigation system would ensure uniform water distribution, which would result in homogeneous water consumption patterns across fields.However, several factors, including deficient infrastructure, inadequate management practices, soil type, water quality, and fertilization, can lead to non-uniform water use across fields and zones.The amount of irrigation water applied in each event was measured through the use of rain gauges positioned above the crop canopy (see Figures 7 and 8).This methodology was proposed and implemented by María [47].

Crop Water Use and Water Use Efficiency
During the period from January 2021 to May 2021, a total of 14 Center pivot (C) fields were subjected to rigorous monitoring to ascertain their irrigation schedules, crop growth, and soil moisture profiles, enabling the assessment of water productivity.The irrigation schedules were established by the farmers, who relied on weather forecasts to optimize their plans.In this regard, the irrigation plan was adjusted to compensate for the lack of spring precipitation, thereby ensuring a reliable supply of water for crop development in the long term.However, this approach was unable to account for the non-linear effect of water shortages during specific phonological stages on crop yields.Irrigation was carried out using the Center pivot technique, and the monthly amount of water was distributed evenly over two to three irrigation events, spaced two weeks apart.The first two rounds of watering were conducted in January, following a principle that was empirically established by the farmers as optimal for maintaining the moisture content of Vertisol, while minimizing water losses due to frequent irrigation.The remaining fields at each site served to verify the measurements obtained from the intensively monitored fields.Groundwater served as the primary source of irrigation water in the wheat fields.While

Crop Water Use and Water Use Efficiency
During the period from January 2021 to May 2021, a total of 14 Center pivot (C) fields were subjected to rigorous monitoring to ascertain their irrigation schedules, crop growth, and soil moisture profiles, enabling the assessment of water productivity.The irrigation schedules were established by the farmers, who relied on weather forecasts to optimize their plans.In this regard, the irrigation plan was adjusted to compensate for the lack of spring precipitation, thereby ensuring a reliable supply of water for crop development in the long term.However, this approach was unable to account for the non-linear effect of water shortages during specific phonological stages on crop yields.Irrigation was carried out using the Center pivot technique, and the monthly amount of water was distributed evenly over two to three irrigation events, spaced two weeks apart.The first two rounds of watering were conducted in January, following a principle that was empirically established by the farmers as optimal for maintaining the moisture content of Vertisol, while minimizing water losses due to frequent irrigation.The remaining fields at each site served to verify the measurements obtained from the intensively monitored fields.Groundwater served as the primary source of irrigation water in the wheat fields.While direct measurements are typically more accurate, transpiration and evapotranspiration are challenging to measure and must be derived from the Cropwat table).
In Table 3, the validation of irrigation measurements was conducted for 14 wheat fields (C1, C2, C3, . . ., C14) from January to May.The validation included a comparison of the weekly evolution of irrigation and observed actual ET, along with the amount of irrigation and rainfall for each of the 14 wheat fields.The validation fields exhibited a high level of coherence, indicating the effectiveness of the irrigation hours in mitigating crop water stress through supplementary irrigation.In April, the center pivot of C12, C13, and C14 with regular irrigation had higher vegetation density and yield, as evidenced by NDVI, whereas fields in other center pivots had lower density.Integrating remote sensing techniques and meteorological data can provide valuable insights into the performance of an ideal irrigation system, which can enhance crop productivity.The study results suggest that remote sensing and meteorological data integration can facilitate the identification of optimal irrigation schedules and strategies for specific crops and regions.By monitoring vegetation health using remote sensing techniques, farmers and irrigation managers can adjust irrigation schedules to match crop water requirements, thereby reducing water waste and improving crop productivity.More frequent irrigation may be necessary in hot and dry conditions, whereas irrigation may not be required for a certain period.

NDVI and NDMI Time Series of Wheat Crop
In the period from January to May 2021, some center-pivot fields were planted with a single crop type, while others had alternating crop types, resulting in variations in NDVI and NDMI time series behavior.To examine this variation, NDVI and NDMI were compared for the Wheat Crop.Specifically, NDVI and NDMI were identified from 143 fields where only Wheat Crop was planted throughout the study period.The fields, each with a diameter of approximately 800 m, were divided into 2500 pixels, and the mean value for each field was determined by averaging the NDVI and NDMI of these pixels.Crop statistics were derived from the NDVI and NDMI time series data for the wheat crop, including maximum (Max) and minimum (Min) values, amplitude variation (Amp), and standard deviation (Std), to establish differences in time series behavior (see Appendix B Tables A2 and A3).In March and April, NDVI and NDMI values fluctuated significantly more in irrigated lands than in non-irrigated lands, possibly due to the influence of ETo and temperature on vegetation growth.To assess the changes in vegetation, the changes in classes (1, 2, 3, and 4) for NDVI and Dense, Moderate, Sparse, and Open soil classes were calculated.The category with the greatest change relative to the other four was then identified as the dominant category.
To demonstrate changes over time, the relative dominance of each class category was computed for each period.The NDVI and NDMI indices revealed that the non-irrigated portions of the field experienced drought conditions in February, March, and April, as indicated by the low NDVI values.The severity and extent of the drought varied over time and space.The values for January to May displayed spatial and temporal variability across the region, which can be attributed to the prevailing temperature and humidity conditions.The fluctuation in NDVI values is influenced by weather-related parameters such as rainfall, temperature, and humidity levels.Temporal variations in climate, soil, and temperature cause changes in NDVI values.

Evolutions of Observed NDVI and NDMI
While Figures 7 and 8 focused on the relationship between NDMI-based soil moisture and NDVI-based vegetation health in Wheat Crop, we aim to assess the possibility of detecting irrigation events using Sentinel-2 data by analyzing their association with realworld data.As illustrated in Figure 8, the NDMI in the center pivot area that was irrigated was notably higher than that in the nearby bare soil area.
In the Center pivot field, soil moisture levels were high due to irrigation, while the other field showed significant variability in soil moisture levels before June 2021.The wheat crop was grown and harvested frequently from January to June 2021, with irrigation durations ranging from 12 to 96 h, as per farmer records.The field was planted with maize from August 2020 to December 2020, and irrigation maintained high soil moisture levels during that period.After the maize was picked, irrigation stopped until January 2021.This led to a sharp drop in the amount of water in the soil.However, soil moisture increased in mid-February 2021 after an 18 mm precipitation event from 11-15 February, which was captured by the Meteorological Station and resulted in an increase in both NDVI and NDMI (Figure 8).High temperatures are believed to impact vegetation strength and cause plant stress.The mean NDMI values increased in all center pivot irrigated fields, consistent with the trend of NDVI variation (see Figures 7 and 8).Although a lack of precipitation is typically the primary cause of drought, elevated potential evapotranspiration is associated with temperature and relative humidity.The evidence presented indicates that soil moisture plays a limiting role in the actual evapotranspiration process, especially during drought conditions, as it prevents excessive drying.Temperature and other factors that influence vegetation development and phenology also have an impact on actual evapotranspiration [48].Thus, NDVI is a reliable sign of drought, and it is important to use empirical NDVI-NDMI-based indices to improve the performance of irrigation systems and increase crop yields.The average NDVI, which reflects the greenness of vegetation, is strongly linked to seasonal rainfall, indicating the potential use of NDVI as a drought predictor.

Correlation Matrix between NDVI and NDMI
Figure 9 displays the average NDVI time series profile for 2021, which shows a significant correlation between NDVI and NDMI.The projected green cover increased from about 0.40 to 0.80 between February and April, likely due to irrigation, before dropping to 0.6 by the end of April as irrigation decreased in preparation for the harvest season.The average NDVI increased due to an increase in NDMI but then dropped to 0.20 in May after the maturation and harvesting phase.Sentinel2 was able to detect changes in vegetation or physiological density accurately.All monthly index values were significantly correlated with meteorological parameters from various locations with a 95% confidence level.NDVI increased at several locations in the 143 center pivot and rain-fed area due to increased rainfall, irrigation, and temperature.The combination of moisture content and temperature played a significant role in reducing NDVI values.Ongoing land and crop degradation and yield losses occurred due to a lack of vegetation cover and relative humidity caused by insufficient rainfall during the growing season (Table A1).Increased temperatures can have severe natural consequences, increasing ETo and ETc.This suggests that crop development was not uniformly distributed across the region, even in times of water shortage, and some farmers were able to achieve acceptable crop growth by effectively using water through supplementary irrigation and implementing good agricultural practices.Remote-sensing-based NDVI and NDMI can provide valuable insights into the spatial and temporal distribution of crop irrigation water requirements at the field level, which can be compared to in situ monitoring of observed irrigation water supplied at the irrigation district scale.Acceptable estimates of ET0 and ETc for wheat (73.7, 98.4, 125.1, 229.4,and 319.2 mm/month) and (42.7, 57.1, 93.8, 229.4, and 383.0 mm/month) were obtained under various weather and water management conditions with irrigation periods of 12, 24, 48, 72, and 96 h, respectively, according to farmer records for the months from January to May. Figure 9 provide further details on this.
The significant rise in temperature can be attributed to the absence of rainfall, which caused a lack of moisture and resulted in higher ETo and ETc.One of the consequences of the increased ETc in April and May in the study area is the elevated temperature and decreased NDVI (Figure 9).This indicates the adverse effect of high temperature on the vegetation growth environment, resulting in lower vegetation values and area (NDVI) in the region.Only a few Center pivot areas in the study saw an increase in NDMI values, which was reflected in a vegetation increase in NDVI.Due to its location, Iraq's climate has transformed into a semi-arid climate, with a heavy influence from the Mediterranean climate characterized by hot, dry summers and warm, wet winters.Rainfall, temperature, sunshine, and wind speed have a negative impact on wheat yield and chlorophyll content in semi-arid environments.The study suggests that the irrigation schedule followed by farmers compensated for the insufficient spring rainfall, resulting in sustainable crop growth.However, this approach did not consider the non-linear effects of water scarcity during specific phenological stages on crop yield.The newly proposed approach helps prevent harmful stress on wheat and ensures the proper distribution of water based on agronomic principles in a center pivot system.These results align with the findings of Polinova et al. (2019) [25].The significant rise in temperature can be attributed to the absence of rainfall, which caused a lack of moisture and resulted in higher ETo and ETc.One of the consequences of the increased ETc in April and May in the study area is the elevated temperature and decreased NDVI (Figure 9).This indicates the adverse effect of high temperature on the vegetation growth environment, resulting in lower vegetation values and area (NDVI) in the region.Only a few Center pivot areas in the study saw an increase in NDMI values, which was reflected in a vegetation increase in NDVI.Due to its location, Iraq's climate has transformed into a semi-arid climate, with a heavy influence from the Mediterranean climate The NDVI has been associated with numerous vegetative properties, including photosynthetic capacity.Sufficient irrigation must be provided for a larger wheat yield to encourage strong output intensity.Figures 7 and 8  water supply distribution, showcasing clear shifts from one month to another in both the applied and required irrigation water.The different shapes of the irrigation distribution reveal that farmers' use of irrigation water varies significantly from month to month, which may be influenced by both farmers' irrigation decision making and the small portion of crop water requirements met by rainfall, which was higher in the dry months of January and February 2021.During the research period (March), there was a slight increase in the average NDVI, although the increase was more noticeable.The average NDVI remained steady throughout the growing seasons (March to April).Even though 90% of the precipitation occurred during the growing seasons, annual precipitation varied greatly from 60.8 mm in January to 2.6 mm in May 2021 (Figures 9 and 10).These significant changes seemed to impact crop development, suggesting that sufficient irrigation water was available to support wheat growth throughout the growing season.This study defines water consumption as the amount of water that evaporates and seeps from an agricultural zone rather than the amount used for crop irrigation or diverted for that purpose.It was shown that by combining remote sensing and GIS techniques, a comprehensive evaluation of irrigation efficiency can be achieved for different irrigation systems [49].
during specific phenological stages on crop yield.The newly proposed approach helps prevent harmful stress on wheat and ensures the proper distribution of water based on agronomic principles in a center pivot system.These results align with the findings of Polinova et al. (2019) [25].
The NDVI has been associated with numerous vegetative properties, including photosynthetic capacity.Sufficient irrigation must be provided for a larger wheat yield to encourage strong output intensity.Figures 7 and 8 demonstrate the observed irrigation water supply distribution, showcasing clear shifts from one month to another in both the applied and required irrigation water.The different shapes of the irrigation distribution reveal that farmers' use of irrigation water varies significantly from month to month, which may be influenced by both farmers' irrigation decision making and the small portion of crop water requirements met by rainfall, which was higher in the dry months of January and February 2021.During the research period (March), there was a slight increase in the average NDVI, although the increase was more noticeable.The average NDVI remained steady throughout the growing seasons (March to April).Even though 90% of the precipitation occurred during the growing seasons, annual precipitation varied greatly from 60.8 mm in January to 2.6 mm in May 2021 (Figures 9 and 10).These significant changes seemed to impact crop development, suggesting that sufficient irrigation water was available to support wheat growth throughout the growing season.This study defines water consumption as the amount of water that evaporates and seeps from an agricultural zone rather than the amount used for crop irrigation or diverted for that purpose.It was shown that by combining remote sensing and GIS techniques, a comprehensive evaluation of irrigation efficiency can be achieved for different irrigation systems [49].

Crop Growing Period (NDVI)
Figure 11 displays the vegetation cover trend over a period of five months, indicating an upward trend in January, reaching its highest point in March and April, and declining to the lowest level at the end of May, which is consistent with the reported rates.The

Crop Growing Period (NDVI)
Figure 11 displays the vegetation cover trend over a period of five months, indicating an upward trend in January, reaching its highest point in March and April, and declining to the lowest level at the end of May, which is consistent with the reported rates.The fluctuations in NDVI values during these months are heavily influenced by the precipitation and irrigation water requirements, as evidenced by the strong correlation between NDVI and NDMI, which varied in unison with an R2 of 0.904 (Figure 11).Remote sensing provides a holistic perspective of ground conditions.This study utilized available data to explore the correlation between irrigation schedules, crop yield, and remote sensing data.The final stage of the research involved a comparison of the findings obtained through the analysis of remote sensing data.By considering all the environmental factors that impact crop growth, and by evaluating the performance of irrigation, the study aimed to determine the root causes of variations in NDVI and NDMI values [50].Nonetheless, it is worth mentioning that certain factors may not be pertinent or suitable in all situations [51,52].The study period covered five months from January to May, with one growing season per year for wheat (Figure 4).The variability of climatic influences on drought, temporal changes, and Water 2023, 15, 1605 20 of 32 their relationship with precipitation, irrigation hours, ETo, ETc, and temperature is evident (Figure 12).explore the correlation between irrigation schedules, crop yield, and remote sensing data.The final stage of the research involved a comparison of the findings obtained through the analysis of remote sensing data.By considering all the environmental factors that impact crop growth, and by evaluating the performance of irrigation, the study aimed to determine the root causes of variations in NDVI and NDMI values [50].Nonetheless, it is worth mentioning that certain factors may not be pertinent or suitable in all situations [51,52].The study period covered five months from January to May, with one growing season per year for wheat (Figure 4).The variability of climatic influences on drought, temporal changes, and their relationship with precipitation, irrigation hours, ETo, ETc, and temperature is evident (Figure 12).The NDVI is influenced by ETo, soil temperature, and irrigation, resulting in correlation values of 0.49, 0.55, and 0.76, respectively.Despite this, the growth rate of the crop is significantly affected by the NDMI.Although there is a strong relationship between NDMI and irrigation, with a correlation of approximately 0.91, it is not statistically significant, as shown in Figure 11 and Table A1.By analyzing 143 selected fields, it was found that the correlation between NDVI and NDMI increased the yield output when sufficient irrigation was provided in February, March, and April (Figure 11).Using EXLSTAT, the correlation coefficients for NDVI, NDMI, Irrigation, ETo, and ETc were calculated during the growing season spanning from January to May, an average of five months.These correlation coefficients are illustrated in Figures 12-14  The findings revealed a strong and positive association between ETc and both NDVI and NDMI.Similarly, there was a positive correlation observed between NDVI and NDMI.These correlations between remote sensing-derived spectral indices and meteorological parameters were statistically significant.After 2020, the planting of maize ceased in the field, and the surface was bare until September 2021.Both time series of NDVI and NDMI present fluctuations in vegetation, with peaks up to 1 and 0.5 for NDVI and NDMI, respectively, occurring during growing seasons.High values with a maximum NDVI of 1 and NDWI of 0.5 were observed in March 2021, which may have been caused by the presence of supplementary irrigation.The study also noticed through the trend of increasing the value of NDVI that the effect of supplementary irrigation during the months of the study was different, and this change mainly depended on the amount of precipitation falling for each month in the study area, and other environmental factors that also play a prominent role in the growth of the wheat crop (Figure 12).The results are consistent with multiple investigations [53][54][55].y = 0.7312x -0.1859 0.50 The NDVI is influenced by ETo, soil temperature, and irrigation, resulting in correlation values of 0.49, 0.55, and 0.76, respectively.Despite this, the growth rate of the crop is significantly affected by the NDMI.Although there is a strong relationship between NDMI and irrigation, with a correlation of approximately 0.91, it is not statistically significant, as shown in Figure 11 and Table A1.By analyzing 143 selected fields, it was found that the correlation between NDVI and NDMI increased the yield output when sufficient irrigation was provided in February, March, and April (Figure 11).Using EXLSTAT, the correlation coefficients for NDVI, NDMI, Irrigation, ETo, and ETc were calculated during the growing season spanning from January to May, an average of five months.These correlation coefficients are illustrated in Figures 12-14.
and NDWI of 0.5 were observed in March 2021, which may have been caused by the presence of supplementary irrigation.The study also noticed through the trend of increasing the value of NDVI that the effect of supplementary irrigation during the months of the study was different, and this change mainly depended on the amount of precipitation falling for each month in the study area, and other environmental factors that also play a prominent role in the growth of the wheat crop (Figure 12).The results are consistent with multiple investigations [53][54][55].Temperature is a key factor affecting the transpiration process of crops, as demonstrated by the significant fluctuation in ET0 values between January and May, corresponding to the growth and harvest stages of the crop (Figure 13).Additionally, the wheat fields studied through farmer interviews were located in 14 Center pivot fields.The fields had similar irrigation schedules, with an average of 48 irrigation hours distributed over two periods in January due to sufficient rainfall.However, in February, with a decrease in the precipitation rate and the application of fertilizers, the irrigation requirements of the fields increased to one hour, resulting in the need for three irrigation periods.PCA was utilized to determine the primary modes of variability in a dataset that examined the temporal changes in NDVI and NDMI in relation to meteorological parameters.NDVI and NDMI are remote sensing indices that are commonly used to monitor vegetation and moisture content, respectively (Figure 14).By examining the changes in these indices over time and their connection to temperature and precipitation, PCA identified the crucial factors influencing the changes in vegetation and moisture patterns.To The findings revealed a strong and positive association between ETc and both NDVI and NDMI.Similarly, there was a positive correlation observed between NDVI and NDMI.These correlations between remote sensing-derived spectral indices and meteorological parameters were statistically significant.After 2020, the planting of maize ceased in the field, and the surface was bare until September 2021.Both time series of NDVI and NDMI present fluctuations in vegetation, with peaks up to 1 and 0.5 for NDVI and NDMI, respectively, occurring during growing seasons.High values with a maximum NDVI of 1 and NDWI of 0.5 were observed in March 2021, which may have been caused by the presence of supplementary irrigation.The study also noticed through the trend of increasing the value of NDVI that the effect of supplementary irrigation during the months of the study was different, and this change mainly depended on the amount of precipitation falling for each month in the study area, and other environmental factors that also play a prominent role in the growth of the wheat crop (Figure 12).The results are consistent with multiple investigations [53][54][55].
Temperature is a key factor affecting the transpiration process of crops, as demonstrated by the significant fluctuation in ET0 values between January and May, corresponding to the growth and harvest stages of the crop (Figure 13).Additionally, the wheat fields studied through farmer interviews were located in 14 Center pivot fields.
The fields had similar irrigation schedules, with an average of 48 irrigation hours distributed over two periods in January due to sufficient rainfall.However, in February, with a decrease in the precipitation rate and the application of fertilizers, the irrigation requirements of the fields increased to one hour, resulting in the need for three irrigation periods.PCA was utilized to determine the primary modes of variability in a dataset that examined the temporal changes in NDVI and NDMI in relation to meteorological parameters.NDVI and NDMI are remote sensing indices that are commonly used to monitor vegetation and moisture content, respectively (Figure 14).By examining the changes in these indices over time and their connection to temperature and precipitation, PCA identified the crucial factors influencing the changes in vegetation and moisture patterns.To carry out the PCA analysis, a covariance matrix of the dataset was constructed, and the matrix was eigendecomposed to obtain the principal components, which are linear combinations of the original variables that represent unique patterns of variability in the data.Typically, the first few principal components account for most of the variance in the dataset and are used to summarize the dominant modes of variability.The PCA of the temporal pattern changes in NDVI and NDMI in relation to meteorological parameters yielded valuable insights into the underlying factors responsible for changes in vegetation and moisture patterns over time and can aid in the identification of critical variables for monitoring and predicting these changes.

Conclusions
The NDMI is a valuable tool for detecting moisture deficiencies in crops and identifying under-irrigated areas.This information can be used to divide fields into zones with different water needs and schedule precision irrigation events as needed.Crop Monitoring provides historical and current precipitation data, NDVI index graphs, and precipitation monitoring graphs.By analyzing seasonal weather and precipitation patterns, it is possible to plan precision irrigation strategies for different fields.In addition, Crop Monitoring offers a 5-day weather forecast for each field, which helps farmers decide on the need for watering activities to ensure proper soil moisture for crops.Precision irrigation is cost effective, providing sufficient water supply to crops in areas with limited rainfall.Crop Monitoring tracks changes in the NDVI for individual fields throughout the season.This helps farmers identify areas of weak and strong productivity across the field and create special maps for variable-rate applications of seeds and fertilizers.NDVI values can vary throughout the growing season, indicating water stress or waterlogging, and can be visualized through maps and graphs.Crop Monitoring also provides current and historical soil moisture data and NDVI index graphs, which help farmers track the correlation between rainfall and moisture levels in the field.NDMI values also vary throughout the growing season and can be visualized through maps and graphs, indicating water stress or waterlogging.A decrease in NDMI values suggests water stress, while abnormally high values could signal waterlogging.Visualization of NDMI through maps and graphs helps farmers detect problem areas in the field and save time and resources.Water supply is a critical factor for plant growth, along with sunlight, nutrients, and soil temperature.While fields in areas with frequent rainfall receive sufficient water for crop growth, additional watering is necessary to maximize yields in semi-arid regions [56][57][58].
For a crop to progress from germination to the reproductive growth stage, where it can begin to produce grain, it requires at least 100 mm of water.It is worth mentioning that, according to the local meteorological station data, the 2021 crop season had the least precipitation and was classified as a drought season.In this study, the focus was on quantifying the loss of water due to evaporation from soil and transpiration from plants, which is collectively known as evapotranspiration.The most commonly used method for calculating evapotranspiration is to determine the reference evapotranspiration (ETo), which is the amount of water that would be lost from a standardized grass surface under specific conditions [59].
The Penman-Monteith equation is often used to calculate ETo, taking into account weather variables such as temperature, humidity, wind speed, and solar radiation.To adjust ETo for different crop types, a crop coefficient (Kc) is used, which represents the ratio of the actual evapotranspiration of a crop to the reference evapotranspiration.The Kc value is dependent on the crop growth stage, canopy cover, and other factors, and can be determined from tables or estimated using empirical formulas.Crop evapotranspiration (ETc) is calculated by multiplying ETo with Kc, and it represents the amount of water required to meet the crop's water needs.ETc is calculated for each crop growth stage and added up over the entire growing season to estimate irrigation water requirements, which is the amount of water needed to supplement natural precipitation to meet the crop's water requirements.Irrigation water requirements can be calculated for different time periods, such as monthly, depending on the capacity of the irrigation system.This study engaged in collaboration with farmers and implemented a systematic approach to determine optimal irrigation water quantity and timing, aiming to satisfy crop water requirements, improve yield, reduce water usage, and address environmental concerns.The outcomes of this research align with findings from prior investigations.Furthermore, it was observed that the elevated temperatures in the growing season and spring led to heightened evaporation rates, necessitating longer irrigation periods [56][57][58][59].
This suggests a considerable need for additional irrigation during the drought year 2021.This study investigated the use of remote sensing techniques and meteorological data to assess the ideal irrigation system performance scenarios for improving crop productivity in the Erbil province of Iraq.The study period covered five months of the wheat growing season.The results demonstrated that vegetation coverage was significantly affected by climatic factors such as precipitation, irrigation hours, and temperature, particularly in the months of March, April, and May.The study also showed that the use of supplementary irrigation in rain-fed areas is essential to mitigate the impact of drought.The correlation between NDMI and NDVI was significant, likely due to the delayed effect of insufficient precipitation on vegetation.The Sentinel2 satellite provided accurate predictions of irrigation conditions, and the developed model was suitable for different time scales, ensuring the reliability and efficiency of water requirement predictions.Severe water stress was observed throughout the rain-fed area in March, April, and May, despite the widespread drought conditions in Iraq during that year.However, north Erbil's center pivot irrigation areas were unaffected by the lack of precipitation and water stress [56,57,60].
This study employed meteorological data and remote sensing-based indices to identify optimal irrigation schedules and water use anomalies.Including agricultural production and surface evaporation data could further improve the assessment of these factors during drought periods.Integrating remote sensing and evapotranspiration models can support irrigation managers in optimizing irrigation schedules and detecting areas with irrigation problems, leading to more efficient and sustainable irrigation practices.However, remote sensing techniques have limitations, such as inadequate spatial and temporal resolutions and inaccurate calibration.While continuous monitoring of the same field over multiple years is impractical due to the study area's crop rotation method, precision irrigation systems that depend on integrated remote sensing and meteorological data can significantly impact agriculture and water resource management.Such systems can provide up-todate information on crop health and moisture levels, which can be used to design and optimize precision irrigation systems that deliver water and nutrients where and when needed [56][57][58]60].This approach can increase crop yields and quality, reduce water usage, and mitigate adverse environmental impacts such as soil erosion, nutrient leaching, and water runoff.Adopting precision irrigation systems based on remote sensing can also help address water scarcity issues at both the global and local levels, given that agriculture accounts for most freshwater usage worldwide.The development and adoption of such systems represent a significant advancement in agriculture and water resource management, which has the potential to enhance food security and environmental sustainability.Despite some moderate research linking wheat crops with NDMI and NDVI, there is still little investigation of wheat crop monitoring based on these indices in the KRI.

Figure 1 .
Figure 1.Site map of the study area in Erbil, KRI.

Figure 1 .
Figure 1.Site map of the study area in Erbil, KRI.

Figure 2 .
Figure 2. Monthly precipitation, relative humidity, actual evaporation, maximum, minimum, and mean temperature of North Erbil and surrounding areas recorded for 2021 years.

Figure 2 .
Figure 2. Monthly precipitation, relative humidity, actual evaporation, maximum, minimum, and mean temperature of North Erbil and surrounding areas recorded for 2021 years.

Water 2023 , 34 Figure 3 .
Figure 3. Flowchart of the methodology adopted in this study.

Figure 3 .
Figure 3. Flowchart of the methodology adopted in this study.

Figure 3 .
Figure 3. Flowchart of the methodology adopted in this study.

Figure 5 .
Figure 5. Temporal Variation of the NDVI Value-Based Vegetation Density Classes of Wheat Field in 2021.Figure 5. Temporal Variation of the NDVI Value-Based Vegetation Density Classes of Wheat Field in 2021.

Figure 5 .
Figure 5. Temporal Variation of the NDVI Value-Based Vegetation Density Classes of Wheat Field in 2021.Figure 5. Temporal Variation of the NDVI Value-Based Vegetation Density Classes of Wheat Field in 2021.

Figure 6 .
Figure 6.Temporal Variation of the NDMI Value-Based Vegetation Density Classes of Wheat Field in 2021.

Figure 6 .
Figure 6.Temporal Variation of the NDMI Value-Based Vegetation Density Classes of Wheat Field in 2021.

Figure 7 .
Figure 7. Temporal Variation of the NDVI Value of Centre pivot fields in 2021.Figure 7. Temporal Variation of the NDVI Value of Centre pivot fields in 2021.

Figure 7 .
Figure 7. Temporal Variation of the NDVI Value of Centre pivot fields in 2021.Figure 7. Temporal Variation of the NDVI Value of Centre pivot fields in 2021.

Figure 8 .
Figure 8. Temporal Variation in the NDMI Value of Center pivot fields in 2021.

Figure 8 .
Figure 8. Temporal Variation in the NDMI Value of Center pivot fields in 2021.

Figure 10 .
Figure 10.Temporal Variation of the NDVI and NDMI Value Based Vegetation Density Classes of 143 Center pivot wheat Field in 2021.

Figure 12 .
Figure 12.Temporal variation of the NDVI and NDMI value-based vegetation density classes of Center pivot wheat field from 1 January to 21 May 2021.

Figure 13 . 34 Figure 13 .
Figure 13.NDMI Correlation with NDVI-wheat in 143 Center pivot wheat field products with a 5-day temporal resolution.

Figure 14 .
Figure 14.Principal Component Analysis of temporal Pattern changes of NDVI and NDMI with Meteorological parameters.

Figure 14 .
Figure 14.Principal Component Analysis of temporal Pattern changes of NDVI and NDMI with Meteorological parameters.

Table 1 .
Statistical indices of measured NDVI value, classes density, and area of each class.
Water 2023, 15, x FOR PEER REVIEW 11 of 34

Table 2 .
Statistical indices of measured NDMI Value, and area of each Classes.

Table 2 .
Statistical indices of measured NDMI Value, and area of each Classes.

Table 3 .
The fourteen fields were monitored during five months of irrigation, with irrigation hours used for each center pivot by month.

Irrigation/h/mm Location Longitude Latitude Yield Kg/ Hectare January February March April May Total/ Hours Mm/ Hours
NDMI values for the High Values category in the months of January through May.The values for this category are 0.28, 1.68, 3.72, 0.04, 65.96, 115.88, 239.64, 337.28, 877.52, 979.24, 5.12, 367.44, 37.7, 5.48, 2.52, and 1.84, respectively.High Values in the table indicate areas with abundant water in the vegetation, which signifies particularly healthy vegetation.
Temporal Variation of the NDVI and NDMI Value Based Vegetation Density Classes of 143 Center pivot wheat Field in 2021.