Soil Moisture Depletion Modelling Using a TDR Multi-Sensor System, GIS, Soil Analyzes, Precision Agriculture and Remote Sensing on Maize for Improved Irrigation-Fertilization Decisions †

: The effects of three drip irrigation (IR1: Farmer’s, IR2:Full (100%ET c ), IR3:Deﬁcit (80%ET c ) irrigation), and two fertilization (Ft1, Ft2) treatments were studied on maize yield and biomass by applying new agro-technologies (TDR—sensors for soil moisture (SM) measurements, Precision Agriculture, Remote Sensing—NDVI (Sentinel-2 satellite sensor), soil-hydraulic analyses and Geostatistical models, SM-rootzone modelling-2D-GIS mapping). A daily soil moisture depletion (SMDp) model was developed. The two-way-ANOVA statistical analysis results revealed that irrigation (IR3 = best) and fertilization treatments (Ft1 = best) signiﬁcantly affect yield and biomass. Deﬁcit irrigation and proper fertilization based on new agro-technologies for improved management decisions can result in substantial improvement on yield (+116.10%) and biomass (+119.71%) with less net water use ( − 7.49%) and reduced drainage water losses ( − 41.02%).


Introduction
Maize is a major irrigated crop worldwide that requires large quantities of water seasonally (400-800 m 3 ) [1]. Worldwide, the maize cultivated area and yield are 194.18 million ha and 1118.56 million metric tons (2019/2020) [2]. The agricultural sector accounts for 70% of global [1,3,4] and for 59% of Europe's [5] freshwater withdrawals. At present, many countries worldwide are experiencing a scarcity of fresh water [1,3,4] for potable and irrigation use. Global water demand is projected to increase by 55% between 2000 and 2050 and the climate change will have adverse impact on world water resources and food production with high degree of regional variability and scarcity [6], with the irrigation water amount being the main factor limiting crop production [1,4,7].
With the ever-increasing competition for finite water resources worldwide and the steadily rising demand for agricultural commodities, the call to improve the efficiency and productivity of water use for crop production, to ensure future food and crop products security and address the uncertainties associated with climate change, has never been more urgent [8]. The global challenge for the coming decades will be increasing crop production with less water using precision agriculture as a management strategy that helps farmers to improve crop production and water efficiency. The aim of the study was to determine the effects of three drip irrigation (IR1: Farmer's, IR2: Full [100%ET c ], IR3: Deficit [80%ET c ] irrigation), and two fertilization (Ft1, Ft2) treatments on maize yield and biomass by applying new agro-technologies.

Experimental Plot Design, Soil Sampling and Laboratory Soil and Hydraulic Analysis
The 2.90 ha field had a factorial split plot design with main factor the three drip irrigation treatments: (a) IR1: Farmer's, (b) IR2: Full (100%ET c ) and c) IR3: Deficit (80%ET c irrigation. The sub factor was two fertilization treatments: Ft1: N-P-K = 333.63-9.16-34.86 Kg ha −1 , and Ft2: N-P-K = 270.03-9.16-34.86 Kg ha −1 ). A GPS receiver was used to identify locations of soil samples that were collected at depth 0-40 cm and analyzed at the laboratory. The soil's pH was measured in a 1:2 soil/water extract with a glass electrode and a pH meter. Soil organic matter was analyzed by chemical oxidation with 1 mol L −1 K 2 Cr 2 O 7 and titration of the remaining reagent with 0.5 mol L −1 FeSO 4 [9]. Cation-exchange capacity (CEC) was analyzed by (i) saturation of cation exchange sites with Na by "equilibration" of the soil with pH 8.2, 60% ethanol solution of 0.4N NaOAc-0.1N NaCI; and (ii) extraction with 0.5N MgNO 3 . Total Na and CI were determined in the extracted solution [9]. The soil's nitrate and ammonium nitrogen were extracted with 0.5 mol L −1 CaCl 2 and estimated by distillation in the presence of MgO and Devarda's alloy, respectively. Available phosphorus P (Olsen method) was extracted with 0.5 mol L −1 NaHCO 3 and measured by spectroscopy [9]. Exchangeable potassium K forms were extracted with 1 mol L −1 CH 3 COONH 4 and measured with a flame Photometer. Field Capacity (FC) and wilting point (WP) were measured with the porous ceramic plate method with 1/3 Atm for FC and 15 Atm for WP [1]. Maize (Zea mays L., var. 1921 Hybrid) was seeded at the end of March and harvested at the end of September.

Remote Sensing Crop's NDVI, Evapotranspiration and Net Irrigation Requirement
Climatic data were obtained from a nearby meteorological station. The effective rainfall was calculated according to USDA-SCS (1970) [13]. The NDVI Vegetation Index [1] was calculated every week using remote sensing (RS) data (Sentinel-2 Satellite sensor) for studying spatial crop development and coefficients. The reference evapotranspiration was computed based on the F.A.O. Penman-Monteith method [1,7,8] and the net irrigation requirement (NIR) was calculated using a soil-water-crop-atmosphere model [1,7]. The crop evapotranspiration (ET c ) and actual evapotranspiration (ET a ) were computed using crop coefficients obtained from remote sensing-NDVI [1,7,8]. The moisture rootzone depletion was calculated using a daily soil moisture depletion (SMDp) model (1) [1,7,8]: where: Dr,(i) = root zone depletion at the end of day i (mm), Dr,(i − 1) (TDR) = TDR sensors measured rootzone soil-water content at the end of day i − 1 (mm), Pe = effective rainfall (mm), NIR(i) = net irrigation requirement (mm), GW(i) = groundwater contribution (capillary rise) from water table on day i (mm), ET c (i) = evapotranspiration on day i (mm), DP(i) = drainage water loss out of the root zone by deep percolation on day i (mm).

Study Area, Soil-Hydraulic Analysis and 2-D Moisture Maps Utilizing GIS, Precision Agriculture and Geostatistics
The experiment was carried out in a farm field located at Viotia Prefecture in Central Greece. The study area is characterized by a typical Mediterranean climate [1,4,7] with a cold winter, hot summer and low precipitation in spring and summer. Results found from sensors and analyses were used as input variables to delineate soil moisture profile digital 2D-GIS maps of the maize's root zone. SM is the major factor for crops' enhanced growth and production [1,7]. Spatial analysis revealed an excellent moisture distribution.

Daily Soil Moisture Depletion (SMDp) Model and NDVI Vegetation Index
The NDVI (Normalized Difference Vegetation Index) [1] was calculated using RS data (Sentinel-2 satellite sensor) for monitoring spatial crop development and coefficients for the ET c and SMDp model (Figure 1). The NDVI vegetation index quantifies crops' vegetation by measuring the difference between near-infrared (which vegetation strongly reflects) and red light (which vegetation absorbs). The daily TDR moisture measurements and 2D-GIS SM mapping, the accurate estimated ET c and monitoring of system inflows and estimated surface outflow and NDVI vegetation index mapping, showed a reliable daily Soil Moisture Depletion model [1,7] for the four stages of maize crop growth.

Study Area, Soil-Hydraulic Analysis and 2-D Moisture Maps Utilizing GIS, Precision Agriculture and Geostatistics
The experiment was carried out in a farm field located at Viotia Prefecture in Central Greece. The study area is characterized by a typical Mediterranean climate [1,4,7] with a cold winter, hot summer and low precipitation in spring and summer. Results found from sensors and analyses were used as input variables to delineate soil moisture profile digital 2D-GIS maps of the maize's root zone. SM is the major factor for crops' enhanced growth and production [1,7]. Spatial analysis revealed an excellent moisture distribution.

Daily Soil Moisture Depletion (SMDp) Model and NDVI Vegetation Index
The NDVI (Normalized Difference Vegetation Index) [1] was calculated using RS data (Sentinel-2 satellite sensor) for monitoring spatial crop development and coefficients for the ETc and SMDp model (Figure 1). The NDVI vegetation index quantifies crops' vegetation by measuring the difference between near-infrared (which vegetation strongly reflects) and red light (which vegetation absorbs). The daily TDR moisture measurements and 2D-GIS SM mapping, the accurate estimated ETc and monitoring of system inflows and estimated surface outflow and NDVI vegetation index mapping, showed a reliable daily Soil Moisture Depletion model [1,7] for the four stages of maize crop growth.

Statistical Analysis, Maize's Yield and Biomass Results
The two-way-ANOVA statistical analysis (p = 0.05) using IBM-SPSS (v.26) [1,7,14,15] revealed that the irrigation treatments (IR3: Deficit (80%ET c ) irrigation (best)) and the fertilization treatments (Ft1: (best)) significantly affect maize's yield and biomass. Whenever a prolonged shortage of soil moisture occurs during the most sensitive growth stages (flowering (L mid stage) and grain filling (L late stage)) ( Figure 1), this usually has as a result that the maize's crop growth is reduced because of less or excess water supply by the farmers and the yield is consequently lower [1,7,8]. Deficit irrigation and proper fertilization based on new agro-technologies for improved management decisions can result in substantial improvement in yield (+116.10%) and biomass (+119.71%) with less net water use (−7.49%) and reduced drainage losses (−41.02%).

Conclusions
Prolonged shortage of soil moisture often occurs during the most sensitive growth stages (flowering (L mid stage) and grain filling (L late stage)) of many crops. In practice, the usual outcome is that the maize's crop growth is reduced because of less (water stress) or excess (saturation) water supply by farmers who cannot estimate the correct soil moisture depletion, the MAD AD (Management Allowed Depletion) and the net water requirements of the crop, so yield is consequently lower. Proper fertilization and deficit irrigation based on new agro-technologies for improved management decisions can result in substantial improvement on yield (+116.10%) and biomass (+119.71%) with less net water use (−7.49%) and reduced drainage water losses (−41.02%) and sustainable management.