Spatiotemporal Modeling of Soil Water Dynamics for Site-Specific Variable Rate Irrigation in Maize

: This study aimed to simulate dynamic irrigation management zones (MZs) in two maize fields for a variable rate hose reel fertigation machine (VRFM) with a four-section boom control. Soil moisture content was measured from nine and four soil moisture sensors in Field 1 (8.2 ha) and Field 2 (2.5 ha), respectively, on different dates during the 2022 crop season. Three and five MZs scenarios were simulated per irrigation and the theoretical maps were processed for implementation. The application maps fitted to the VRFM showed significant spatiotemporal variations in irrigation requirements. For instance, in Field 1, 3-MZ modelling showed that the areas requiring high (H), medium (M), and low (L)-level irrigation on 21 July were 1.60, 4.84, and 1.85 ha, respectively, even though the farmer applied uniform rate over the whole field. H-level sub-areas ranged between 1.22 ha (25 July) and 3.25 ha (7 July), showing a coefficient of variation (CV) of 43.32% for the three MZs, whereas H-level sub-areas for the five MZs varied from 0.41 ha (2 July) to 1.49 ha (7 July) with a CV value of 48.84%. High levels of within-field variability can be addressed using precise and dynamic irrigation MZs fitted to the irrigation technology used.


Introduction
Water shortages due to climate change are becoming more common and significantly affect agricultural production [1].Commonly used in arid farming regions and during times of low rainfall, irrigation is the process of adding water to the soil to increase crop growth and replant degraded soil [2].Irrigation depends on regional climatic conditions; however, as climate change begins to affect new regions that were not affected by drought in the past, the need for irrigation also arises in these regions.According to research and current practices, farmers' use of irrigation is increasing, even in northern EU countries.In countries where irrigation is widely used, irrigation technologies have improved substantially.Consequently, it is reasonable to assume that irrigation will receive greater attention in parts of the world where it was not a priority before.Precision irrigation is necessary because of the rising demand for water to combat climate concerns.To achieve higher agricultural yields and water savings, it is crucial to choose when and how much to irrigate.The four most popular methods of operation for irrigation scheduling, namely evapotranspiration and water balance (ET-WB), soil moisture status, plant water status, and crop models, along with their advantages and disadvantages, are introduced and contrasted to better understand their fundamental processes and principles [3].It is necessary to have knowledge of soil properties, including texture, water-holding capacity (WHC), wetness, and crop water requirements, at different growth stages [4].
With the advancement of contemporary sensors and sensing technology, a variety of technologies, such as wireless sensor networks (WSNs), remote sensing (RS), global navigation satellite system (GNSS), geographic information systems (GIS), LoRaWAN, and Bluetooth can now be used to collect field data for precision irrigation.WSNs have been widely used in agriculture for automated irrigation management because they are costeffective and have low energy consumption [5][6][7][8].In addition to wired solutions, wireless sensor technologies have become more important and widely used [9], particularly in weather forecasting, early warning systems, soil and crop monitoring, and automated and smart irrigation [10][11][12].Wireless sensor networks with cloud capabilities are promising because the end-user does not have to handle collecting and storing the data.
It has been well documented that parameters affecting the yield of any crop exhibit spatial and temporal variations [13][14][15][16].Despite awareness of this, the great majority of farmers today are known to apply a constant rate over the field while they vary the total amount of water from one irrigation to another during the irrigation season depending on the soil, plant, and weather conditions.At present, a limited number of farmers seem to determine irrigation needs based on climate data and some sensor data, if any.One of the most important factors affecting the adoption of precision or smart irrigation technologies is that the accuracy of the moisture sensor changes with time in the field, necessitating occasional calibration.Other factors that cause difficulties for the average end-user include cost, installation, data collection, and interpretation issues.Although the advantages of variable rate irrigation have been documented, the adoption rate of this technology has been low; however, efforts to develop the technology should continue for several reasons, including climate change, declining water resources, nutrient management regulations, and water management policies [17].
The development of variable-rate irrigation control systems has attracted considerable attention since the early 2000s.Some of these studies focused on center-pivot and linearly moving irrigation systems [18][19][20].The development of smart irrigation systems is of the utmost interest and requires the accurate determination of all parameters affecting irrigation [10,21].Smart systems require not only sensor-based monitoring but also the integration of Internet of Things (IoT) and cloud technologies [10,[22][23][24][25][26].Ref. [27] reviewed various water management techniques and the use of artificial intelligence (AI) and deep learning (DL), along with the IoT network and case studies.Although previous studies have covered new techniques to apply irrigation water more intelligently, there is still a need to improve the current applications by considering the dynamic variations in irrigation needs within fields more precisely.There have been studies focusing on variable rate irrigation (VRI) using hose reel machines, but they usually use the full-swath width as the lateral resolution [28][29][30][31][32].There appears to be a research gap regarding improvements in the lateral and longitudinal resolutions of VRI applications using hose reel machines.Furthermore, spatial and temporal changes during the irrigation period should be addressed in the management of VRI.
The general objective of this study was to simulate spatiotemporal variations in irrigation needs of maize crops based on dynamic irrigation management zones using a wireless soil moisture sensor network and GIS capabilities.This study compared the constant user rates of irrigation water on different irrigation dates with that of the sensorbased simulated variable irrigation rate.It assesses the right time and rate for irrigation, determines the accuracy of user rate applications, and analyses the dynamic irrigation management zones to be applied by VRFM.

Site and Irrigation System
The experimental maize fields were located in Bakırköy village of Karacabey, Turkey (40 • 08 ′ 22.61 ′′ N, 28 • 22 ′ 58.7 ′′ E).An 8.2 ha field named Field 1 (F1) and a 2.5 ha field named Field 2 (F2) were selected for the experiment (Figure 1).The soil texture in these fields was composed of 37.6% clay and 32.3% sand in Field 1, and 30.3% clay and 35.3% sand in Field 2. The type of the soil in this study area is alluvial [33].The fields have approximately rectangular shapes and are irrigated using a closed-loop pressurized wide-area drip irrigation system operated by an irrigation union in the district.The field experiment was conducted during the 2022 growing season.Maize seeds were planted on 16 May 2022 for Field 1 and on 11 May 2022 for Field 2, with the first drip irrigation performed on July 2 for both fields.

Site and Irrigation System
The experimental maize fields were located in Bakırköy village of Karacabey, Turkey (40°08′22.61″N, 28°22′58.7″E).An 8.2 ha field named Field 1 (F1) and a 2.5 ha field named Field 2 (F2) were selected for the experiment (Figure 1).The soil texture in these fields was composed of 37.6% clay and 32.3% sand in Field 1, and 30.3% clay and 35.3% sand in Field 2. The type of the soil in this study area is alluvial [33].The fields have approximately rectangular shapes and are irrigated using a closed-loop pressurized wide-area drip irrigation system operated by an irrigation union in the district.The field experiment was conducted during the 2022 growing season.Maize seeds were planted on 16 May 2022 for Field 1 and on 11 May 2022 for Field 2, with the first drip irrigation performed on July 2 for both fields.The drip irrigation system used during the experiment had drippers delivering 28.8 t h −1 ha −1 irrigation water at an operating pressure of 3 bar.The main and secondary plastic pipes were 200 mm in diameter each.Soil samples were collected randomly at 80 and 30 locations from F1 and F2, respectively, to determine the soil properties (pH, P, K, Mg, Ca, clay, sand, and silt) pertinent to the determination of irrigation needs.

Wireless Soil Moisture Sensor System
A cloud-based LoRaWAN wireless network was used with soil moisture sensors to collect soil moisture data continuously (Figure 2).Twelve and six MC sensors were installed horizontally at a depth of 300 mm for F1 and F2, respectively.The gateway was approximately 800 m away from the sensors with a clear view of sight.The quantity and placement of the Dragino LSE01 sensors were determined according to the relative management zones (MZs) developed for irrigation.These MZs were delineated based on the precise spatial data on predicted soil MC and clay.The sensors monitored the soil MC The drip irrigation system used during the experiment had drippers delivering 28.8 t h −1 ha −1 irrigation water at an operating pressure of 3 bar.The main and secondary plastic pipes were 200 mm in diameter each.Soil samples were collected randomly at 80 and 30 locations from F1 and F2, respectively, to determine the soil properties (pH, P, K, Mg, Ca, clay, sand, and silt) pertinent to the determination of irrigation needs.

Wireless Soil Moisture Sensor System
A cloud-based LoRaWAN wireless network was used with soil moisture sensors to collect soil moisture data continuously (Figure 2).Twelve and six MC sensors were installed horizontally at a depth of 300 mm for F1 and F2, respectively.The gateway was approximately 800 m away from the sensors with a clear view of sight.The quantity and placement of the Dragino LSE01 sensors were determined according to the relative management zones (MZs) developed for irrigation.These MZs were delineated based on the precise spatial data on predicted soil MC and clay.The sensors monitored the soil MC ((%) VWC), apparent electrical conductivity (ECa) (uS cm −1 ), and soil temperature ( • C).
The sensors enable MC measurements using frequency domain reflectometry (FDR); the accuracy of the MC sensors was reported to be ±3% by the manufacturer.
((%) VWC), apparent electrical conductivity (ECa) (uS cm −1 ), and soil temperature (°C).The sensors enable MC measurements using frequency domain reflectometry (FDR); the accuracy of the MC sensors was reported to be ±3% by the manufacturer.
The accuracy of the sensors was checked during the irrigation season once by collecting soil samples at sensor locations, resulting in approximately 3-5% error in soil moisture measurements.No further investigations were undertaken to check or improve the measurement accuracy of MC sensors.

Methods
The variable irrigation needs were determined by monitoring the water status using soil moisture content (MC) sensors by previous researchers [35,36] and the same method was used for the experimental fields in this study.Field capacity (FC), permanent wilting point (PWP), and volumetric weight (VW) were determined under laboratory conditions for each soil sample at each sensor location.Undisturbed soil samples were collected.VV was determined by drying soil samples in an oven for 24 h at 105 °C.The sand/kaolin box instrument (Eijkelkamp, Giesbeek, The Netherlands) was used to determine field capacity at 2.54 pF (0.34 bar).The permanent wilting point at a pressure of 4.2 pF (15.20 bar) was determined using a pressure membrane apparatus (Eijkelkamp, Giesbeek, The Netherlands).Available water (AW), management allowable depletion (MAD), and amount of irrigation water (irrigation need, IN) were determined based on laboratory analyses of the soil samples.IN was calculated based on the moisture content (MC) value of the soil to reach the FC values at a given location in the field.The following formulae were used to determine the aforementioned parameters: where depth is the plant root depth (300 mm in this study).From Equations (1), (2), and (5), the absolute values for FC, PWP, and MC in mm of water column were calculated for the root depth using the relative FC, PWP, and MC in percent and the measured volumetric weight of the soil.Available water is the difference between FC and PWP (Equation ( 3)).It represents the range of moisture that is theoretically available for plants, although at moisture levels nearing PWP, plants will start experiencing water stress.To prevent water stress, a management allowable depletion (MAD) is set, which is located between FC and PWP.In The accuracy of the sensors was checked during the irrigation season once by collecting soil samples at sensor locations, resulting in approximately 3-5% error in soil moisture measurements.No further investigations were undertaken to check or improve the measurement accuracy of MC sensors.

Methods
The variable irrigation needs were determined by monitoring the water status using soil moisture content (MC) sensors by previous researchers [35,36] and the same method was used for the experimental fields in this study.Field capacity (FC), permanent wilting point (PWP), and volumetric weight (VW) were determined under laboratory conditions for each soil sample at each sensor location.Undisturbed soil samples were collected.VV was determined by drying soil samples in an oven for 24 h at 105 • C. The sand/kaolin box instrument (Eijkelkamp, Giesbeek, The Netherlands) was used to determine field capacity at 2.54 pF (0.34 bar).The permanent wilting point at a pressure of 4.2 pF (15.20 bar) was determined using a pressure membrane apparatus (Eijkelkamp, Giesbeek, The Netherlands).Available water (AW), management allowable depletion (MAD), and amount of irrigation water (irrigation need, IN) were determined based on laboratory analyses of the soil samples.IN was calculated based on the moisture content (MC) value of the soil to reach the FC values at a given location in the field.The following formulae were used to determine the aforementioned parameters: where depth is the plant root depth (300 mm in this study).From Equations ( 1), (2), and (5), the absolute values for FC, PWP, and MC in mm of water column were calculated for the root depth using the relative FC, PWP, and MC in percent and the measured volumetric weight of the soil.Available water is the difference between FC and PWP (Equation ( 3)).It represents the range of moisture that is theoretically available for plants, although at moisture levels nearing PWP, plants will start experiencing water stress.To prevent water stress, a management allowable depletion (MAD) is set, which is located between FC and PWP.In Equation ( 4), MAD setting was considered as 50% of the available water (AW), according to [37].
The irrigation efficiency (EF) was taken (Equation ( 6)) as 90% for the drip irrigation system [38], to calculate sensor-based irrigation needs based on field capacity and moisture measurement.The irrigation time was based on MAD + 30% limit, which was taken to be 30% of the difference between the FC and MAD values.
The irrigations started before the sensors were ready to collect real-time MC data.After processing the sensor data, only eight irrigation applications were considered for the analyses.

Water Use Comparison between User-and Sensor-Based Constant Irrigation Rates
The farmer used inline linear drippers (Eurodrip-Eolos, Rivulis, Gvat, Israel) in the irrigations with an irrigation intensity of 28.8 t h −1 ha −1 at an operating pressure of 3 bar.The pressure at the end of each secondary line was measured using a handheld analog pressure gauge (Pakkens 0-6 bar) so that the flow rate was the same and the water distribution was uniform in each irrigation plot.For each irrigation, the farmer varied the irrigation time from 6 to 10 h per irrigation, depending on the job and employee status.The user irrigation rate was determined by the experience of the farmer, crop water demands given in the irrigation guide for the Bursa-Karacabey district, and evapotranspiration values provided by the farmer's weather station (Metos IMT280).
Wireless MC sensors were configured to deliver MC data every 20 min in this study.To confirm a sensor-based constant irrigation recommendation on a given irrigation day, the MC value measured just before irrigation was obtained from each sensor, and the average of all sensor readings was used as the sensor-based field-scale soil MC value.Equation ( 6) was used to determine the amount of water to be applied (irrigation need) to bring the moisture content of the field to the field capacity for a given irrigation date.The differences among the soil moisture sensor measurements for each irrigation date were calculated to reveal the range of soil moisture content within the field.Within-field soil moisture content variations were evaluated to assess the accuracy of uniform irrigation in the trial fields.A descriptive statistics analysis was performed in terms average, difference, standard deviation, and coefficient of variance (CV).

Variable Rate Irrigation Recommendations
The ordinary kriging (OK) interpolation method was implemented using ArcGIS ArcMap v10.8 (ESRI Redlands, Redlands, CA, USA) software to generate variable rate irrigation (VRI) maps.The irrigation needs (IN, mm da −1 , where 1.0 da = 0.1 ha) were mapped with three and five management zones (MZs) for each field to compare the effects of different resolutions for irrigation.The three and five MZs were selected to see the effects of the number of MZs on the variability in irrigation needs.
First, a theoretical application map was produced for the irrigation dates (2, 5, 7, 13, 17, 21, 25 July and 4 August) for F1 and F2, resulting in 32 maps (8 days × 2 different zone maps × 2 fields).The amount of water to be applied in each irrigation MZ was tabulated based on the irrigation rate and the total area of each irrigation management zone in each field.The user-applied uniform rates and variable-rate irrigation recommendations were then compared to meet the irrigation demands of the soil based on the postulated irrigation zones.
Then, the theoretical 3-and 5-zone irrigation need maps were further processed for VRI application to be applied using a variable rate hose reel fertigation system being developed as part of the project.The VRFM has a working width of 46 m, and it is proposed to have four sections, each to be individually controlled for VRI applications.Thus, using theoretical maps, applicable irrigation maps were obtained for cell sizes of 11.5 m × 1 m.Applicable maps were presented for the irrigation dates of 7, 17, and 25 July and 4 August to discuss the effect of management zones on the water distribution within the fields and were compared with theoretical IN maps.
Sentinel 2A satellite data were used to generate the Normalized Difference Vegetation Index (NDVI) maps corresponding to the nearest irrigation dates.The Sentinel-2A im-ages were downloaded from https://scihub.copernicus.eu/dhus/#/home(accessed on 14 September 2023) in bottom-of-atmosphere L2A format.A total of eight images from 1 July, 6 July, 13 July, 21 July, 26 July, and 5 August 2022, covering the study area, atmospherically corrected, and cloudless were used for analysis.For each S2 scene, the red (band 4), and NIR (band 8) spectral bands with a 10 m spatial resolution were selected for the NDVI calculation according to [39] formula.
Additionally, the texture (clay, silt, and sand) was mapped for F1 and F2 in an attempt to find potential cause-and-effects of NDVI and soil texture on the spatiotemporal variations in irrigation water needs across the fields studied.

User-Rate and Sensor-Based Constant Irrigation Comparisons
Nine out of twelve sensors were utilized for soil moisture measurements due to unexpected weak signals, battery, or installation problems in three of the sensors installed in F1; four of the six sensors were operable in F2 (Table 1).The average FC, PWP, and MAD values in F1 were approximately 139, 66, and 37 mm, respectively.The average volumetric weight of the soil samples was 1.45 (g cm −3 ) with a small coefficient of variation (CV) value (1.82%).Other parameters included in Table 1 had variations with CV values ranging from approximately 4.2 to 7.3%.In F2, based on measurements performed by the four soil moisture sensors, the FC was smaller and VW was higher compared with F1.The MAD values were smaller than those of F1.The coefficient of variation of the measured parameters in F2, except for PWP, ranged from 2.02 to 8.46%.The high CV value of PWP may be attributed to the fewer sensors used, and hence, to fewer data points to be averaged.The irrigation duration determined by the farmer varied from 7 to 10 h depending on the weather, crop conditions, and labor availability (Table 2).More water should have been applied throughout the season to align the soil MC with the FC.For instance, on July 7, the applied user rate was 23.04 t da −1 as a result of 8 h irrigation, which was 14.7% less than the recommended rate based on the average of the soil MC measurements.The overall average of user rates was 22.68 mm da −1 while the sensor-based application rate was calculated to be 29.83mm da −1 .The overall difference between the user and sensor-based rates was −22.76%, whereas the differences ranged from −4.6% to −35.9% on specific dates.This means that the soil received a smaller amount of water than it should have if the moisture sensor data would have been used to calculate the irrigation rate.It can be concluded that the irrigation duration could be adjusted more accurately by the farmer if moisture sensor data were used to align the available MC with the FC level.In F2, the ranges of irrigation rates for the user-and sensor-based irrigations were 17.28-28.80mm da −1 and 27.0-47.19mm da −1 , respectively (Table 2).On average, throughout the irrigation season, the user-based rate of irrigation water application was 45.09% less than the recommended IN.In this case, the difference was more significant than that for F1.Although the farmer intended to apply water every 3-4 days after 7 July, due to labor shortage, an irrigation was skipped between 7 July and 13 July.Automated irrigation systems can reduce the dependency on workers to accomplish time-dependent tasks.
It should be noted that the above discussion does not account for spatial variations in irrigation needs within a field.Even for a constant irrigation rate at the whole-field scale, the use of soil moisture sensors revealed opportunities for implementing irrigation at more accurate times and rates.In [40], it was concluded that the total profit was the same when sensor-based irrigation scheduling was compared with the user application, but the sensor-based application was found to be promising in terms of water and financial savings.
The variations in the measured MC values on specific irrigation dates are shown in Figure 3 for F1 and F2. Figure 3 also presents the average FC, PWP, and MAD values to visualize when irrigation had to be started and terminated.
The MC values were expected to be around the MAD line before the irrigation started and near the FC line after the irrigation was completed.The farmer always started irrigation at a soil moisture content slightly below the MAD + 30% value in F1 (Figure 3a).The irrigation timings were slightly after than the exact irrigation times.After irrigation, the soil MC was always lower than the FC.The farmer applied irrigation water that did not exceed the FC for the entire field.However, the timing and duration of the user rate can be improved by using the average soil moisture readings made by the sensors so that the soil MC values reach FC.
same when sensor-based irrigation scheduling was compared with the user application, but the sensor-based application was found to be promising in terms of water and financial savings.
The variations in the measured MC values on specific irrigation dates are shown in Figure 3 for F1 and F2. Figure 3 also presents the average FC, PWP, and MAD values to visualize when irrigation had to be started and terminated.The MC values were expected to be around the MAD line before the irrigation started and near the FC line after the irrigation was completed.The farmer always started irrigation at a soil moisture content slightly below the MAD + 30% value in F1 (Figure 3a).The irrigation timings were slightly after than the exact irrigation times.After irrigation, the At the early growth stages of maize plants, the penetration of the maize roots is shallower than in the later stages [3]; thus, maintaining the topsoil moisture should be more important to avoid water stress.The soil tended to generate cracks owing to the lack of moisture, especially during windy periods that are common in the region.Therefore, it may be recommended to start irrigation earlier than occurred because the measured MC values before irrigation were lower than the MAD + 30% values.MC values below the MAD + 30% value in Figure 3a,b suggest that plants might have experienced water stress during the early growth stages.To prevent this, frequent and sufficient irrigation could have been determined using the sensor-based data.
Figure 3b suggests that the timing of irrigation was not accurate for some irrigation dates in Field 2. The farmer started irrigating the field at a soil MC level lower than the MAD + 30% value in four of the irrigations.Furthermore, irrigation was terminated too early since the after-irrigation soil MC values were just above the MAD + 30% values and not reaching the targeted FC values.As demonstrated in Figure 3b, real-time data could be used to accurately determine the correct length and amount of irrigation, either done manually or automatically.
FC and after-irrigation MC values were compared on each irrigation day (Table 3).The MC of the soil was approximately 110 mm before irrigation and increased to 128.2-137.5 mm while the average FC was 138.9 mm for the field.The percent user difference in irrigation needs between the user and the average sensor readings varied from 1% to 7.6%, implying that the user should have applied slightly more water to bring the MC values to the FC level at the field scale.F1, field 1; F2, field 2; FC_Ave, average field capacity; MC_BI, before-irrigation moisture content; MC_AI, afterirrigation moisture content; Diff., absolute difference between the soil moisture content and the field capacity; %diff., difference between MC_AI and FC.
Table 3 shows the user applications compared with the after-irrigation MC-based measurements for F1 and F2.The user started irrigation when the soil MC value was between 86 mm and 97 mm in the field, depending on the irrigation date.The MC values following irrigations were between 107 and 126 mm.The applied water was less than the amount of water necessary to increase the soil MC values up to an average FC level of 128 mm in Field 2. Up to 17% more water should have been applied in Field 2 on 2 and 5 July, corresponding to the early growth stages of plants.
The spatial variation in the before-irrigation MC values for each irrigation day is shown in Table 4, where the first column shows the sensors installed in both fields.On average, the IN was between 30 and 35.5 mm before irrigation in F1.However, the IN was calculated to range from 35.94 to 47.19 mm for F2, which was a wider range than for F1.
The average irrigation need for each irrigation date was greater in F2 than in F1, probably because of the different ranges between the FC and PWP, and hence the MAD values in the two fields.The average MAD value of F1 was greater than that of F2, requiring more water delivery to irrigate the field up to field capacity.

Constant vs. Variable Rate Irrigation
As discussed in the previous subsection, irrigation scheduling based on soil moisture sensors has the potential to improve the timing and amount of irrigation water required.However, the above discussion was based on the average of all readings from the soil moisture sensors.As the sensors were in different parts of the field, it would be more useful to determine the spatial irrigation needs on a given irrigation day.Table 5 shows the difference between the average before-irrigation MC values and point-based MC values on each irrigation day.The data show that some parts of the field required more or less water than the previously calculated average IN.For instance, on 13 July, 11.5% more water would have been applied at S2, whereas 12.1% less water would have been applied at S10 if the average before-irrigation MC values had been used.Similarly, in F2, a maximum of 12.7% more and 11.5% less water would have been applied on 25 July at S8 and 4 August at S14, respectively, if the average of all MC sensors had been used to determine the irrigation need.The spatiotemporal changes in the soil MC values show that the irrigation needs vary within the fields; hence, appropriate precision application technologies should be developed to account for these variations (Figure 4).The variations in the soil MC in F1 during the monitoring period showed that the upper-left corner of the field always had lower moisture content values and the lower-left corner had the greatest measured MC values.Other parts of the field showed varying degrees of change in the level of water requirements.
Figures 4 and 5 show the three-and five-zone irrigation management zones for F1 and F2.These figures show the same patterns in irrigation needs with different precisions with low (L), medium-low (ML), medium (M), medium-high (MH), and high (H) application rates.
Soil Syst.2024, 8, x FOR PEER REVIEW 15 of 26 irrigation systems in use.Thus, it may be concluded that the maps obtained by the OK method are useful for revealing the spatiotemporal variabilities within the fields but need to be processed for use with the most appropriate irrigation technology available, for example, hose reel or pivot systems with, at minimum, section control mechanisms.The Irrigation maps for F2 in Figure 5 show that the areas requiring the least amount of water moved from the top-left corner (2 July and 5 July) to the bottom-left (7 July) and bottom-right (13 July) sides of the field.The locations of the areas needing the least amount of water continued to change for the remaining F2 during the irrigation season.The irrigation simulation maps in Figures 6 and 7 show different behaviors compared to F1, in that the dynamic changes in the spatial IN requirements can be prominent over time.
Table 6 shows the partitioning of the sub-section areas of the three-and five-zone irrigation maps for F1 and F2 derived from the maps in Figures 4 and 5, respectively.The analysis of five irrigation management zones suggests that irrigation water can be delivered more precisely to areas needing different irrigation rates, for example, M, ML, L, MH, and H levels.However, accurately applying five different application rates based on the management zones is unlikely considering the layout of the current drip irrigation systems in use.Thus, it may be concluded that the maps obtained by the OK method are useful for revealing the spatiotemporal variabilities within the fields but need to be processed for use with the most appropriate irrigation technology available, for example, hose reel or pivot systems with, at minimum, section control mechanisms.
The Irrigation maps for F2 in Figure 5 show that the areas requiring the least amount of water moved from the top-left corner (2 July and 5 July) to the bottom-left (7 July) and bottom-right (13 July) sides of the field.
The locations of the areas needing the least amount of water continued to change for the remaining F2 during the irrigation season.The irrigation simulation maps in Figures 6 and 7 show different behaviors compared to F1, in that the dynamic changes in the spatial IN requirements can be prominent over time.Table 6 shows the partitioning of the sub-section areas of the three-and five-zone irrigation maps for F1 and F2 derived from the maps in Figures 4 and 5, respectively.
The H level irrigation rate corresponded to areas from 0.89 ha to 32.6 ha with five-zone analysis in F1 whereas three-zone analysis resulted in 0.14 ha to 14.3 ha for H level rates on different irrigation days.Dividing the field into management zones shows that different parts of the field require different application rates to bring the soil MC to the FC.Five zones depicted the distribution of the rates more precisely for each irrigation date.In F2, the zone areas were smaller due to the smaller field size compared to that of F1.The CV values from about 45-65% both for F1 and F2 also show that the application rates for each zone need to be changed significantly over time.For instance, the CV level of 46.7% for M level area in the five-zone analysis was found in areas ranging from 3.0 ha (21 July) to 13.87 ha (4 August).The IN maps shown in Figures 5-8 are not readily applicable because of the lack of precise application technologies that can deliver irrigation water based on contour maps.The length of the lateral pipes in a drip irrigation system is designed to be the same length in each management zone in sub-plot areas of up to four hectares in the current practices in the Karacabey district, intersecting with the simulated IN maps.Furthermore, the dynamic variation in the irrigation management zone from one irrigation to the next requires an irrigation system that can respond accordingly.

VRI Maps for Hose Reel Irrigation Machine
Since an irrigation system in a fixed-size management zone is unlikely to respond accurately to spatiotemporal irrigation needs, a linearly moving hose reel irrigation system with section control was proposed as a solution in this study.For a four-section control machine, the IN maps were translated into grid maps so that the hose reel irrigation machine could achieve a lateral application resolution of 11.5 m (with a working width of 46 m) and a longitudinal resolution with sub-meters due to the low ground speeds (30-80 m h −1 ).Management zones with appropriate section widths that can be utilized by the VRFM are shown in Figure 6 for four selected irrigation days, both with three-and fivezone mapping.
Figures 6 and 7 show that the ideal irrigation management zones present in Figures 4 and 5, respectively, have to be changed due to the overlaps among the intersecting management zones when cells are formed with a width of 11.5 m with varying lengths.The lateral and longitudinal resolutions of the machine and the number of zones will be the initial determinants for the overall precision of the variable-rate irrigation application.

VRI Maps for Hose Reel Irrigation Machine
Since an irrigation system in a fixed-size management zone is unlikely to respond accurately to spatiotemporal irrigation needs, a linearly moving hose reel irrigation system with section control was proposed as a solution in this study.For a four-section control machine, the IN maps were translated into grid maps so that the hose reel irrigation machine could achieve a lateral application resolution of 11.5 m (with a working width of 46 m) and a longitudinal resolution with sub-meters due to the low ground speeds (30-80 m h −1 ).Management zones with appropriate section widths that can be utilized by the VRFM are shown in Figure 6 for four selected irrigation days, both with three-and five-zone mapping.
Figures 6 and 7 show that the ideal irrigation management zones present in Figures 4 and 5, respectively, have to be changed due to the overlaps among the intersecting management zones when cells are formed with a width of 11.5 m with varying lengths.The lateral and longitudinal resolutions of the machine and the number of zones will be the initial determinants for the overall precision of the variable-rate irrigation application.
Table 7 provides the updated sub-areas of each management zone after the maps were fitted to the widths of the sections of the VRFM.When rectangular cells with fixed widths and varying lengths were used, the total area of each management zone changed accordingly.For instance, three-zone H-level application areas for F1 ranged from 1.22 ha to 3.67 ha for different irrigation dates, which were 0.90 ha to 3.26 ha in Table 6.For F2, L-level areas with five-zone analysis ranged from 0.29 ha (4 August) to 0.98 ha (21 July) in Table 7, whereas these areas were 0.27 ha and 0.99 ha, respectively for 4 August and 21 July in Table 6.In terms of the variations among irrigation days, the CV values were high in the application maps fitted to the VRF machine, ranging from 25.77 to 67.4%.
The means of three-zone areas for H, M, and L levels were 2.22, 4.34, and 1.71 ha in Field 1 and 0.50, 1.24, 1.24, and 0.75 ha in Field 2. The large area differences among the zones show that the constant-rate applications introduce large errors in the application rates.Similar evaluations can be performed for a specific irrigation date.For instance, five-zone areas corresponding to 17 July in F1 range from 0.53 ha (H level) to 4.1 ha (M level) with the largest area, showing up to eight-fold difference in area between different zones.It is documented in the literature that precision agriculture has the potential to reduce the excess amount of inputs in production agriculture, along with an increase in profit and energy saving [41,42].Ref. [43] found large irrigation need differences within a study field, with a water requirement of up to 26% less in the management zones compared to user-rate applications and suggested that sensor-based variable-rate irrigation scheduling would prevent users from applying over-or under-irrigation.Neither variations in NDVI nor the soil texture seemed to explain the sharp spatial or temporal changes in irrigation needs in this study.From visual comparisons, Figures 8 and 9 could not reveal insight about the reasons for varying irrigation MZs within the field since there was no sign of correlation between the variations in MZs and NDVI or soil texture.Previous research has also concluded that NDVI data alone are not sufficient to make decisions in irrigation management [44].More data, such as spatial emergence rate, plant density, leaf area index, and pest (weed and disease) pressure on plant growth, may reveal the necessary information to explain the cause of varying irrigation management zones.
The application rates associated with the dynamic management zones simulated in this study were based on measurements using multiple number of soil moisture sensors installed in fields integrated with a cloud-based wireless system.A greater number of sensors may improve the accuracy of application maps, albeit at a higher cost.A limitation of this study was to utilize only the MC data.Some researchers have utilized areal measurements, particularly NDVI, to reduce the daily soil moisture variance [20,45] and/or to develop a smart system for linear irrigation systems, whereas others have focused on individual sprinkler control in continuously moving irrigation systems [46][47][48][49][50]. Artificial intelligence (AI) techniques have been used to improve the decision-support process in irrigation scheduling [51][52][53][54].
field since there was no sign of correlation between the variations in MZs and NDVI or soil texture.Previous research has also concluded that NDVI data alone are not sufficient to make decisions in irrigation management [44].More data, such as spatial emergence rate, plant density, leaf area index, and pest (weed and disease) pressure on plant growth, may reveal the necessary information to explain the cause of varying irrigation management zones.The application rates associated with the dynamic management zones simulated in this study were based on measurements using multiple number of soil moisture sensors installed in fields integrated with a cloud-based wireless system.A greater number of sensors may improve the accuracy of application maps, albeit at a higher cost.A limitation of this study was to utilize only the MC data.Some researchers have utilized areal measurements, particularly NDVI, to reduce the daily soil moisture variance [20,45] and/or to develop a smart system for linear irrigation systems, whereas others have focused on individual sprinkler control in continuously moving irrigation systems [46][47][48][49][50]. Artificial intelligence (AI) techniques have been used to improve the decision-support process in irrigation scheduling [51][52][53][54].

Conclusions
Determining the dynamic irrigation MZs, and hence the instantaneous irrigation needs in a given field is of utmost importance to optimize water use within the fields.In this study, the irrigation needs of two maize fields were determined based on spatiotemporal soil moisture content variations using LoRaWAN-based soil moisture sensor technology in an attempt to simulate the dynamic irrigation management zones to be used with a hose reel irrigation machine.Wireless soil MC monitoring technology enabled the determination of accurate timing and rates for irrigation to meet the field capacity requirement.Simulations showed that management zones changed temporally within the fields, requiring an update of the irrigation map for the next irrigation.It was demonstrated that resolutions of the theoretical (contour) irrigation maps could be converted to applicable grids or cells to be used by advanced variable rate irrigation systems.A section-controlled variable rate hose reel irrigation machine (VRIM) is currently under development to implement the solution postulated in this study.Future work should adopt data fusion of weather, electrical conductivity (EC), leaf area index (LAI), and satellite thermal and NDVI data in an attempt to improve irrigation management zones for a variable rate hose reel irrigation system.Future irrigation scheduling methods should also focus on the management of soil moisture based on an advanced understanding of its effects on crop growth, either by integrating existing irrigation scheduling methods or developing new models using advanced algorithms.

Figure 1 .
Figure 1.Study site (Field 1 and Field 2) with fields sub-plots for the drip irrigation.

Figure 1 .
Figure 1.Study site (Field 1 and Field 2) with fields sub-plots for the drip irrigation.

Figure 3 .
Figure 3.The soil moisture content (MC) before irrigation (BI) and after irrigation (AI) and their relationship with field capacity average (FC_Ave), permanent point average (PWP_Ave), and management allowable depletion (MAD) for Field 1 (a) and Field 2 (b).

Figure 3 .
Figure 3.The soil moisture content (MC) before irrigation (BI) and after irrigation (AI) and their relationship with field capacity average (FC_Ave), permanent point average (PWP_Ave), and management allowable depletion (MAD) for Field 1 (a) and Field 2 (b).

Figure 4 .
Figure 4. Three and five irrigation management zones maps for Field 1 for different irrigation dates in 2022: (a,b), 2 July; (c,d), 5 July; (e,f), 7 July; (g,h), 13 July; (i,j), 17 July; (k,l), 21 July; (m,n), 25 July; and (o,p), 4 August.The variations in the soil MC in F1 during the monitoring period showed that the upper-left corner of the field always had lower moisture content values and the lower-left corner had the greatest measured MC values.Other parts of the field showed varying degrees of change in the level of water requirements.Figures4 and 5show the three-and five-zone irrigation management zones for F1

Figure 8 .
Figure 8. NDVI for a given dates (a-g) and soil texture (sand, silt, and clay) (h-j) maps for Field 1.

Figure 8 .
Figure 8. NDVI for a given dates (a-g) and soil texture (sand, silt, and clay) (h-j) maps for Field 1.

Table 1 .
Laboratory results for soil samples taken at the sensor locations for Field 1 and Field 2.

Table 2 .
Irrigation needs to reach field capacity for each irrigation date based on average MC from all sensors for Field 1 and Field 2.

Table 3 .
Percent user difference between after-irrigation moisture content and field capacity (FC) for Field 1 and Field 2.

Table 4 .
Average irrigation needs (mm da −1 ) per irrigation and maximum differences among sensors for Field 1 and Field 2.

Table 5 .
The percent difference between the average sensor readings to point-based moisture content (MC) readings for each irrigation date for Field 1 and Field 2.

Table 6 .
Variation in area (ha) of each management zone over time for variable irrigation based on OK for Field 1 and Field 2.

Table 7 .
Area (ha) of each variable rate irrigation (VRI) management zone for VRFM for Field 1 and Field 2.