Predicting Maize Transpiration, Water Use and Productivity for Developing Improved Supplemental Irrigation Schedules in Western Uruguay to Cope with Climate Variability

Various maize irrigation treatments including full and deficit irrigation were used to calibrate and validate the soil water balance and irrigation scheduling model SIMDualKc at Paysandú, western Uruguay. The model adopts the dual crop coefficient approach to partition actual evapotranspiration (ET c act) into actual transpiration (T c act) and soil evaporation (E s). Low errors of estimation were obtained for simulating soil water content (Root mean square errors (RMSE) ď 0.014 cm 3 ¨ cm´3 with calibrated parameters, and RMSE ď 0.023 cm 3 ¨ cm´3 with default parameters). The ratio E s /T c act ranged from 26% to 33% and E s /ET c act varied from 20% to 25%, with higher values when the crop was stressed offering less soil coverage. Due to rainfall regime, runoff and deep percolation were quite large. The Stewarts phasic model was tested and used to predict maize yield from T c act with acceptable errors, in the range of those reported in literature. Water productivity values were high, ranging 1.39 to 2.17 kg¨m´3 and 1.75 to 2.55 kg¨m´3 when considering total water use and crop ET, respectively. Using a 22-year climatic data series, rainfed maize was assessed with poor results for nearly 40% of the years. Differently, alternative supplemental irrigation schedules assessed for the dry and very dry years have shown good results, particularly for mild deficit irrigation. Overall, results show appropriateness for using SIMDualKc to support the irrigation practice.


Introduction
Maize is a main summer crop in western Uruguay.In the last ten years, the cropped area increased by four times and that of maize by 2.5 [1].That increase was made at the expenses of pasture and grassland.Maize production in Uruguay is constrained by water availability [2], which, as for the Argentinian pampas, is subject to large variability related with the El Niño Southern Oscillation (ENSO) [3].Rainfed maize production largely varies due to a high variability of precipitation.Related adaptation measures include no-till systems [4,5] and supplemental irrigation since most of maize cropped area is rainfed [2].
In rainfed agriculture, mainly when soils have a small water storage capacity, it may not be possible to overcome impacts of the variability of rainfall as referred for maize [6].Thus, it is important to develop knowledge on maize water use in western Uruguay and to develop appropriate irrigation scheduling aimed at achieving the potential yields through efficient use of rainfall and irrigation water.A variety of factors influence irrigation management aimed at efficient use of water, improved water productivity and controlling environmental impacts [7][8][9][10].The selection of appropriate supplemental irrigation scheduling options not only require knowledge of the time dynamics of crops demand for water but also adequate prediction of yield responses to water [11].This justifies the use of models to develop strategies for supporting farmers' irrigation decision-making.
Numerous studies focused on the impacts of water stress on maize yields at various crop stages.Results of those studies (e.g., [12][13][14]) have shown that the most sensitive stages to water stress are emergence (VE stage), flowering which includes the tasseling and silking (VT and R1 stages), and early yield formation (R2 and R3 stages).Several approaches to assess the impacts of water deficits on yields have been developed.These include water-yield functions and crop growth and yield models.Applications of crop yield models to maize are numerous (e.g., [15][16][17][18]).However, these models are very demanding in terms of data and parameterization and may not be better than adopting a combination of a water balance simulation model with a water-yield function.Aiming at supporting irrigation scheduling, the Stewart's model [12,[19][20][21] is largely used (e.g., [13,14] for maize).In recent years, the model was applied to a variety of crops [22][23][24] and maize [25][26][27].Recently, that model was modified to relate yield with crop transpiration [28].
Numerous applications have been published on a variety of methods usable to optimize water and land allocation at basin and system levels [29][30][31], generally focusing multi-users and multi-crop systems.Differently, when focusing single crops, the search for more efficient water use is generally performed through simulation of alternative irrigation schedules which are compared through the respective impacts on total water use, yields and water productivity.With this purpose, both water balance simulation and crop-yield models may be used [6,32,33].This approach was tested for maize using the soil water balance simulation model SIMDualKc [34], which partitions ET into soil evaporation (E s ) and actual crop transpiration (T c act ).Knowing T c act , it is possible to adopt a water yield model based on impacts of transpiration deficits on yields to assess the performance of the studied irrigation schedules [28].Considering the need for an optimal conjunctive use of rainfall and irrigation water, this approach may be particularly useful given the large rainfall variability observed during the maize crop season in western Uruguay.
Based upon maize field experiments adopting various irrigation strategies developed at Paysandú, western Uruguay, the objectives of this study consist of: (1) assessing the various components of maize water use with support of the SIMDualKc model; (2) predicting transpiration, yields and water productivity relative to various irrigation management strategies; (3) developing alternative supplemental irrigation schedules to improve rainfall and irrigation water use; and (4) assessing the performance of simulated irrigation strategies.

Site Characterization
Maize experiments were performed at the Experimental Station "Dr.M.A. Cassinoni" (32 ˝22 1 S, 58 ˝41 W, and 50 m elevation), at Paysandú, western Uruguay.According to the Köppen-Geiger classification [35], climate is a Cfa, warm temperate climate, with humid and hot summers.The average annual temperature is 18.3 ˝C and the average annual precipitation is 1327 mm, however with high inter-annual variability.
Weather data were collected with an automatic meteorological station (Vantage Pro 2TM, Davis Instruments, Hayward, CA, USA) located near the experimental fields.Maximum and minimum air temperature ( ˝C), solar radiation (MJ¨m ´2¨d ´1), wind speed (m¨s ´1), air relative humidity (%), and precipitation (mm) were recorded daily; observations were performed above grass.These data were used to compute daily reference evapotranspiration (ET o ) with the FAO-PM method [20].Main climatic data for both seasons relative to October 2011 to February 2012 and October 2012 to February 2013 are presented in Figure 1.The soil is loamy in the first layer of 0.20 m and clay loamy thereafter.The main physical and hydraulic properties of the soil are presented in Table 1.Crop management was the same as recommended in the area to be adopted by local farmers.Plots cropped with maize were monitored to assess different irrigation schedules.The irrigation system consisted of pressure compensating in-line drippers spaced 0.20 m along the lateral; the operating pressure was 98 kPa; and the emitters discharge was 1.49 L•h −1 .Drip lines were located between every pair of crop rows, which were spaced by 0.70 m.Evaluations of the distribution uniformity (DU) of water applications were performed; an average DU = 90% was observed.Irrigation depths were measured with a flowmeter installed upstream of each plot.Soil water content (SWC) was measured with a calibrated neutron probe (503DR HYDROPROBE, InstroTek Inc., Martinez, CA, USA).Measurements were performed at each 0.10 m until a maximum of 0.80 m.However, considering the inappropriateness of using neutron probes in the surface soil layer [36], SWC was determined in the upper 0.10 m layer through soil sampling followed by oven drying.
The following irrigation strategies were adopted in the crop  The soil is loamy in the first layer of 0.20 m and clay loamy thereafter.The main physical and hydraulic properties of the soil are presented in Table 1.Crop management was the same as recommended in the area to be adopted by local farmers.Plots cropped with maize were monitored to assess different irrigation schedules.The irrigation system consisted of pressure compensating in-line drippers spaced 0.20 m along the lateral; the operating pressure was 98 kPa; and the emitters discharge was 1.49 L¨h ´1.Drip lines were located between every pair of crop rows, which were spaced by 0.70 m.Evaluations of the distribution uniformity (DU) of water applications were performed; an average DU = 90% was observed.Irrigation depths were measured with a flowmeter installed upstream of each plot.Soil water content (SWC) was measured with a calibrated neutron probe (503DR HYDROPROBE, InstroTek Inc., Martinez, CA, USA).Measurements were performed at each 0.10 m until a maximum of 0.80 m.However, considering the inappropriateness of using neutron probes in the surface soil layer [36], SWC was determined in the upper 0.10 m layer through soil sampling followed by oven drying.
The following irrigation strategies were adopted in the crop season of 2012-2013: -FI, full irrigation, aimed at minimizing water stress in all crop growth stages; Irrigations were scheduled using a simplified daily soil water balance applied to the effective root zone depth Z r = 0.75 m and using the observed weather and SWC data.Irrigation depths were set to refill SWC up to 90% of θ FC in the periods when water stress was not allowed and up to 60% of θ FC otherwise.Irrigation were applied whenever a depletion of 60% of the total available soil water (TAW, mm) was attained for the periods when aiming at inducing water stress.Otherwise, irrigations were applied when a depletion of 40% of TAW was reached.TAW was defined as TAW = (θ FC ´θWP ) 1000 Z r .No stress was allowed during emergence in order to assure a good crop establishment.
Water deficits during the targeted periods referred above were induced by withholding irrigation and preventing rainfall with a shelter that was moved to cover the plots when rainfall occurred.The experimental design consisted in completely random blocks with three replications of the above referred five irrigation strategies.Experimental plots were 3.5 m ˆ5 m, similar to those used by Echarte et al. [37] when assessing maize sowing densities for the region.The 1.5 m wide strips bordering the plots were also planted with maize to control advection effects, thus minimizing effects of plot sizes.It resulted that the continuous area cropped with maize totalized 352.5 m 2 .The irrigation schedules applied are described in Table 2.The maize hybrid DK 692 was selected as recommended by the Uruguayan National Cultivars Evaluation System adopting a plant density of 100,000 plants¨ha ´1, which is commonly adopted in the region following the studies by Echarte et al. [37].Densities larger than 80,000 plants¨ha ´1 are commonly used in Europe and in China [28,34,38] but vary with the variety used.The dates of each crop growth stages as proposed in FAO56 [20] and the cumulated growing degree days (CGDD) are presented in Table 3.The crop height (h, m) and the fraction of soil covered by the canopy (f c , dimensionless) are given in Table 4. f c along the crop seasons were visually estimated as the percentage of soil shaded by the crop near solar noon.Rooting depths (Z r , m) were measured in randomly distributed plants throughout the season, and this surveillance was performed collecting soil samples with an Edelman type probe (Eijkelkamp, NL, USA) to a depth of 1.0 m and visually checking for the existence of roots at each 0.10 m depth.The maximum root depth observed was 0.75 m, with most of the roots in the first 0.30 to 0.40 m of soil.
Biomass and grain yield were obtained when harvesting all experimental plots and standard deviations were computed; samples were oven dried to constant weight at 65 ˘5 ˝C and the yield was adjusted and measured at 13.5% grain moisture.

Modeling
The modeling strategy used in the present study consists of combining the soil water balance model SIMDualKc [34], which simulates water use and crop evapotranspiration, with the global (S1) and the multi-phasic (S2) water-yield models [12,19] as shown in the flow chart of Figure 2.
The SIMDualKc model performs a daily soil water balance at field scale, thus computing the water use components of a crop, particularly plant transpiration and soil evaporation using the dual crop coefficient approach (dual K c ) [20,39].The model is appropriate to generate irrigation scheduling alternatives and to support farm irrigation scheduling.The model performance has been assessed relative to its capability to simulate soil evaporation and plant transpiration (e.g., [38,40,41]).The adequateness of the dual K c approach has been independently assessed [42,43].
The model computes potential crop evapotranspiration (ET c , mm) as where ET o is reference evapotranspiration (mm), K cb is the basal crop coefficient (dimensionless) that characterizes crop transpiration and K e is the evaporation coefficient (dimensionless) characterizing soil evaporation.Therefore, the model computes separately crop transpiration T c = K cb ET o and soil evaporation E s = K e ET o .The terms and parameters in Equation ( 1) refer to potential crop evapotranspiration, which occurs when the crop is non-stressed and able to attain maximum yield.Otherwise, the actual ET (ET c act , mm) is computed as a function of the available soil water in the root zone.

Modeling
The modeling strategy used in the present study consists of combining the soil water balance model SIMDualKc [34], which simulates water use and crop evapotranspiration, with the global (S1) and the multi-phasic (S2) water-yield models [12,19] as shown in the flow chart of Figure 2.
The SIMDualKc model performs a daily soil water balance at field scale, thus computing the water use components of a crop, particularly plant transpiration and soil evaporation using the dual crop coefficient approach (dual Kc) [20,39].The model is appropriate to generate irrigation scheduling alternatives and to support farm irrigation scheduling.The model performance has been assessed relative to its capability to simulate soil evaporation and plant transpiration (e.g., [38,40,41]).The adequateness of the dual Kc approach has been independently assessed [42,43].
The model computes potential crop evapotranspiration (ETc, mm) as where ETo is reference evapotranspiration (mm), Kcb is the basal crop coefficient (dimensionless) that characterizes crop transpiration and Ke is the evaporation coefficient (dimensionless) characterizing soil evaporation.Therefore, the model computes separately crop transpiration Tc = Kcb ETo and soil evaporation Es = Ke ETo.The terms and parameters in Equation ( 1) refer to potential crop evapotranspiration, which occurs when the crop is non-stressed and able to attain maximum yield.
Otherwise, the actual ET (ETc act, mm) is computed as a function of the available soil water in the root zone.The total available soil water (TAW, mm) is defined as the soil water storage in the root zone between field capacity (θFC, cm 3 •cm −3 ) and the wilting point (θWP, cm 3 •cm −3 ).The readily available soil water (RAW, mm) is RAW = p TAW where p is the soil water depletion fraction for no stress.Water stress occurs when the available soil water content (ASW, mm) is below the threshold RAW.ASW is computed from the actual soil water content (SWC, cm 3 •cm −3 ) as ASW = (SWC − θWP) 1000 Zr.Water stress is expressed through the stress coefficient (Ks, 0-1) defined [20,39] as a function of the depletion in the effective root zone (Dr): with K 1 for Dr ≤ RAW.The total available soil water (TAW, mm) is defined as the soil water storage in the root zone between field capacity (θ FC , cm 3 ¨cm ´3) and the wilting point (θ WP , cm 3 ¨cm ´3).The readily available soil water (RAW, mm) is RAW = p TAW where p is the soil water depletion fraction for no stress.Water stress occurs when the available soil water content (ASW, mm) is below the threshold RAW.ASW is computed from the actual soil water content (SWC, cm 3 ¨cm ´3) as ASW = (SWC ´θWP ) 1000 Z r .Water stress is expressed through the stress coefficient (K s , 0-1) defined [20,39] as a function of the depletion in the effective root zone (D r ): with K s " 1 for D r ď RAW.ET c act equals ET c when K s = 1.Differently, ET c act < ET c when K s < 1.0, with ET c act given by Under these conditions, the actual plant transpiration is T c act = K s K cb ET o , thus depending upon K s .When stress affects crop growth and the coverage of the ground by the crop canopy K e also changes.K e is computed through a daily water balance of the evaporation soil layer, which is characterized by its depth (Z e , m), the total evaporable water (TEW, mm) and the readily evaporable water (REW, mm).TEW is the maximum depth of water that can be evaporated from the evaporation soil layer when it has been completely wetted, and REW is the depth of water that can be evaporated without water availability restrictions [20,39].The model computes the soil water balance in the root zone in terms of depletion at the end of every day [20,44]: where D r,i and D r,i´1 are the root zone depletion (mm) at the end of, respectively, Day i and Day i´1, P is precipitation (mm), RO is runoff (mm), I is net irrigation depth (mm), CR is capillary rise from a shallow groundwater table (mm), ET c act is the actual evapotranspiration (mm), and DP is deep percolation (mm), all terms referring to Day i. RO is computed using the curve number approach [45].CR was not considered in this study because the water table depth was below 10 m.DP fluxes were computed with the decay parametric equation proposed by Liu et al. [46], which relates the soil water stored after occurrence of a heavy rain or irrigation (W a , mm) with the time t (days) to attain field capacity.In that function, the parameter a D (mm) characterizes storage and b D refers to velocity of vertical drainage; both parameters depend upon the soil physical characteristics [46].Irrigation scheduling options available in SIMDualKc are based upon the concepts of depletion fraction for no stress (p), management allowed depletion (MAD), TAW and RAW.Those options include: (1) irrigation to prevent water stress, when MAD ď p; (2) deficit irrigation, when MAD > p, i.e., the available soil water may be depleted to a threshold level below RAW; (3) user selected irrigation depths and intervals between irrigation events; and (4) no irrigation [34].Furthermore, the computation of the seasonal net irrigation requirements (NIR, mm) as defined by Doorenbos and Pruitt [47] is also available.
The SIMDualKc model requires calibration when used with a crop different of those previous simulated, or, for the same crop, in a different environmental and managerial context as recently discussed [48,49].Calibration focuses on the crop parameters K cb and p, the soil evaporation parameters Z e , TEW and REW, and the RO and DP parameters.Otherwise, default parameters can be used because errors are small when using well selected default parameters [45,50]; thus, the parameters obtained for the calibration at Paysandú should be further used as default in other locations of western Uruguay.
In the present study, aiming at assessing irrigation scheduling impacts on yields, the SIMDualKc model was combined with the water-yield models S1 and S2 proposed by Stewart et al. [12].The model approach S1 assumes a linear relationship between the relative yield loss, RYL = 1 ´Ya /Y m , and the relative evapotranspiration deficit, RED = 1 ´ET c act /ET c .The modified approach [28] consists of adopting a relative transpiration deficit (RTD) instead of RED.This approach is justified because transpiration is the ET component directly responsible for yield formation [51] and, therefore, several crop growth models used T c act for the estimation of biomass and yield [21,52].Thus, the S1 model predicts actual yields as: where Y m is maximum (potential) yield (kg¨ha ´1), T d is the seasonal transpiration deficit (mm), with T d = T c ´Tc act and K y is the yield response factor (dimensionless).The Y a values observed at the experimental plots of Paysandú from 2010-2011 to 2012-2013 combined with the modeled T d were used to obtain the K y as previously proposed [28].The Y m values were obtained from the highest yields achieved for the non-stressed treatments and from yield data information collected from farmers in the study area, however adjusted to climate conditions using the "Wageningen method" [19].This method uses data on the seasonal shortwave radiation, saturation vapor pressure deficit, ET c and the duration of the crop cycle.Y m was 15.7 t¨ha ´1 for 2011-2012; 14.2 t¨ha ´1 in 2012-2013 for FI, DI FLO and DI VEG-REP .Because harvesting was anticipated for DI MAT and rainfed cases, it resulted a smaller Y m = 11.9 t¨ha ´1 (Table 3).The S2 model was also modified to consider T c act [28] and takes into consideration the time periods when water stress occurs, which is particularly important for maize due to its sensitivity to water stress, particularly at flowering and yield formation [12].The S2 model uses a three parameters function: where β v , β f and β m are the yield response factors (dimensionless) relative to the vegetative growth period (VE to V10 stages), the flowering (pollination) period, from first tassel until blister kernel (VT to R2 stages), and the maturation period, from milk stage until physiological maturity (R3 to R6 stages).T d,v , T d,f and T d,m are the transpiration deficits (mm) for the same crop periods.The parameters β v , β f and β m used in the present study were those proposed by Alves et al. [53], respectively, 1.2, 2.8 and 0.9, when the crop was conditioned by water deficits during the vegetative growth period, otherwise, β v = 2.1, β f = 7.9 and β m = 2.1 [53].The S2 model was also tested using all observed Y a data during 2010-2011, 2011-2012 and 2012-2013 as well as the model computed T d .

Calibration and Validation Procedures
Model calibration is the process of adjusting influential model parameters within their reasonable ranges aiming at achieving the best fit of observed SWC by the simulated SWC as discussed by Pereira et al. [49].Model calibration consisted therefore in searching the parameters relative to the crop (K cb , p), the soil evaporation layer (Z e , TEW, REW), runoff (CN) and deep percolation (a D and b D ) that lead to minimal deviations between the simulated and observed SWC values.Model validation consisted in evaluating the accuracy of the model when the calibrated parameters are used to simulate independent observed data sets.Calibration was performed using the DI FLO data set of 2011-2012, which was selected because its completeness of SWC data, and validation was performed for all five irrigation treatments of 2012-2013.
The calibration was performed through an iterative trial and error procedure of searching the best parameters values, first through visualizing the decrease of deviations between observed and simulated SWC values, later by observing when the root mean square errors were decreasing until nearly stabilizing.The trial and error procedure was applied by steps because the number of searched parameters is large.It was first applied to the K cb and p values, then to the runoff (RO), deep percolation (DP) and soil evaporation parameters, lately to the crop parameters again and finally considering all parameters.The search was performed considering the expected range of variation of the parameter values.
The initial values for K cb and p were those tabled by Allen et al. [20], for soil evaporation were those based on values proposed by Allen et al. [39], for RO were those based on CN values proposed by Allen et al. [45], and the initial DP parameters were those proposed by Liu et al. [46] for loamy soils.The water balance was initialized two days prior to sowing when SWC observations were performed.In 2011-2012, the initial soil water depletion in the surface layer (D e ) was null and the SWC in the underneath layers was 7% above field capacity.For 2012-2013, the initial D e was null and D r was also null for treatments DI MAT , DI VEG-REP and rainfed, while D e was 20% of TEW and D r was 15% of TAW for FI and DI FLO .Differences between treatments relate with soil water holding characteristics of the plots.
"Goodness-of-fit" indicators were used to assess the performance of SIMDualKc at calibration and validation.The adopted indicators were computed from the pairs of observed and predicted SWC values, respectively O i and P i (i = 1, 2, ..., n), whose means are, respectively, O and P. Following previous SIMDualKc applications [48,49] and Legates and McCabe Jr. [54], the indicators are: i The regression coefficient b 0 of a regression forced to the origin (FTO) relating O i and P i SWC values, which aim at recognizing how similar were the simulated and observed values, computed as b 0 " ii The determination coefficient of the ordinary least squares regression of the same variables aimed at assessing the fraction of the variance of observations that was explained by the model.
, . - iii The root mean square error (RMSE), which expresses the variance of the residual errors, computed as RMSE " which may vary between 0.0, when a perfect match would occur, and a positive value, which should be smaller than the mean of observations.iv The normalized RMSE (NRMSE), that is defined as the ratio between RMSE and the observations mean O, which expedites the comparison of its values for different variables, computed as v The average relative error (ARE), that expresses the estimation errors as a percentage of observation values ARE " 100 n vi The Nash and Sutcliff [55] modeling efficiency (EF), which expresses the relative magnitude of the mean square error (MSE = RMSE 2 ) when compared with the observed data variance [52]: EF = 1 is the target value and values close to 1.0 indicate that the model performance is very good since then the mean square error is much smaller than the variance of observations.Contrarily, when EF is null or negative, this means that there is no gain in using the model.
The same indicators were used when testing the water-yield parameters of Equations ( 5) and (6).

Generating and Assessing Alternative Supplemental Irrigation Scenarios
Using a 22-year series of weather variables provided by the Uruguayan Meteorological Institute (INUMET), a series of NIR values was obtained with SIMDualKc.Adopting an empirical frequency distribution for the NIR series, it was possible to characterize the climatic demand and selecting the years when the probabilities of NIR exceedance were 20% and 5%, which correspond to dry and very dry conditions [6].
For both the dry and very dry years, supplemental irrigation scheduling (SIS) alternatives were set with SIMDualKc, which was used to design various SIS, including full and deficit irrigation (DI).These SIS are different from the irrigation strategies used in the field trials because the latter aimed at assessing the impacts on yield of water stress imposed at selected crop growth stages while the SIS are aimed at searching alternative schedules leading to improved water use and productivity.The DI schedules were designed to control water stress during the most sensitive crop stages-germination/emergence, flowering/pollination, and maturation-therefore to minimize water stress impacts on yields.Simulations were performed considering the average CGDD as indicated in Table 3.The various SIS were built with fixed net irrigation depths (D = 30 mm) as proposed by García-Petillo [56], and ceasing irrigation 20 days before harvesting.A similar D was adopted by Martins et al. [57].The MAD irrigation thresholds were: The performances of the SIS and rainfed scenarios were assessed using various indicators: total water use (TWU), relative yield decrease (RYD) and water productivity (WP).The latter was computed [8] as the ratio between the predicted actual yield and the seasonal TWU computed as TWU " P e `GI `∆ASW where P e is effective precipitation (mm), i.e., the difference between total precipitation and runoff, GI is gross irrigation (mm), and ∆ASW is the variation of the available soil water (mm) between planting and harvesting, positive when ASW(harvesting) < ASW(planting).In addition, following Pereira et al. [8], the consumptive use WP (WP ET ), often called water use efficiency, which is the ratio between the actual yield and the actual crop evapotranspiration, was also computed.

Soil Water Balance Modeling and Model Parameterization
The calibration and validation of SIMDualKc through comparing predicted and observed SWC was the main modeling issue.Results relative to the calibration and two validation cases are shown in Figure 3 when using both calibrated and default parameters.
Analyzing Figure 3, it can be noticed that water stress occurred in the deficit irrigation case used for calibration (Figure 3a,b) and for the rainfed case used for validation (Figure 3e,f), respectively, from 10 to 30 January 2012 and from 13 January to 1 February 2013.Differently, no water stress occurred for the full irrigation case.
The "goodness-of-fit" indicators relative to all simulated cases using calibrated or default parameters are presented in Table 5. Regression coefficients for all treatments are very close to 1.0, therefore indicating that the predicted SWC values are statistically close to the observed ones.The values of the determination coefficients are generally close to 0.90, varying from 0.76 to 0.95, thus showing that most of the variance is explained by the model.Naturally, the indicators referring to the use of calibrated parameters are superior to those relative to using default parameters.The errors of estimate are small, with RMSE ď 0.014 cm 3 ¨cm ´3 when calibrated parameters were used, and ď0.023 cm 3 ¨cm ´3 when default parameters were used.The corresponding normalized values NRMSE are also quite small, not exceeding 0.07; similarly, ARE did not exceed 3.4% and 6.7%, respectively, when parameters used were calibrated or default.EF ranged from 0.71 to 0.87 when calibrated parameters were used, thus indicating that the mean square error was much smaller than the measured data variance.Moreover, lower but reasonably high EF values were also obtained when using default parameters.Overall, results in Figure 3 and Table 5 show that the model is appropriate for use in western Uruguay for applications aimed at supporting improved irrigation practices and management, including when adopting default parameters if these are well selected.
Water 2016, 8, 309 11 of 22 small (Table 6) because default parameters were well selected and resulted close to the calibrated ones, which is a main reason for the small errors obtained when using default parameters (Table 5).The calibrated potential Kcb mid = 1.05 (Table 6) equals that obtained by Rosa et al. [48] for Portugal.Slightly higher values, up to 1.15, are reported in other studies [20,28,39,57,58].Differently, a smaller Kcb mid was reported by Stricevic et al. [59].The value of the Kcb end depends upon crop management and its value is slightly lower than the one proposed by Allen et al. [20] because maize was harvest at low grain moisture.The p values are similar to those tabled by Allen et al. [20] except p at end season (0.75) because a large soil water depletion was intended by then.
The results above indicate that SIMDualKc is appropriate for further use in Uruguay and that the calibrated parameters are likely appropriate for being used as default ones in further applications; nevertheless, relative to soil evaporation and deep percolation, parameters may need to be adjusted considering the soil characteristics of the locations where the model will be applied.The default (initial) and calibrated model parameters-basal crop coefficients, depletion fractions for no stress, and parameters relative to soil evaporation, runoff and deep percolation-are presented in Table 6.It can be noted that differences between default and calibrated parameters are small (Table 6) because default parameters were well selected and resulted close to the calibrated ones, which is a main reason for the small errors obtained when using default parameters (Table 5).Notes: REW and TEW are the readily and total evaporable water, respectively; Z e is the depth of the soil evaporation layer; CN is the curve number; and a D and b D are the parameters of the deep percolation equation [46].* different values were obtained due to the spatial heterogeneity of the soil.
The calibrated potential K cb mid = 1.05 (Table 6) equals that obtained by Rosa et al. [48] for Portugal.Slightly higher values, up to 1.15, are reported in other studies [20,28,39,57,58].Differently, a smaller K cb mid was reported by Stricevic et al. [59].The value of the K cb end depends upon crop management and its value is slightly lower than the one proposed by Allen et al. [20] because maize was harvest at low grain moisture.The p values are similar to those tabled by Allen et al. [20] except p at end season (0.75) because a large soil water depletion was intended by then.
The results above indicate that SIMDualKc is appropriate for further use in Uruguay and that the calibrated parameters are likely appropriate for being used as default ones in further applications; nevertheless, relative to soil evaporation and deep percolation, parameters may need to be adjusted considering the soil characteristics of the locations where the model will be applied.

Crop Coefficients and ET Partitioning
As previously mentioned, the SIMDualKc model partitions ET into E s and T c act using the dual K c approach.Examples of results referring to the seasonal variation of the potential and actual basal crop coefficients, K cb and K cb act , and the evaporation coefficient K e are presented in Figure 4.In this figure, the K cb mid values obtained through model calibration were adjusted to the actual climate conditions when the average minimum relative humidity differed from 45% and the average wind speed was different from 2 m¨s ´1 [20].Figure 4 also includes irrigation and precipitation depths observed, which are depicted to ease perceiving the K e peaks, which essentially depend upon those wetting events.The K cb act and K cb curves are coincident when full irrigation (FI) was practiced in 2012-2013 (Figure 4a).The representation of wetting events help understanding why during various periods stress occurred in 2012-2013, when deficit irrigation was practiced, with K cb act < K cb in Figure 4b, which correspond to time periods when precipitation and irrigation were insufficient to satisfy the crop demand.Figure 4c shows that water stress of the rainfed crop occurred only during a brief period in the vegetative growth and maturation stages (K cb act < K cb ) because distribution of precipitation events along the maize season was favorable; thus, in 2012-2013, rainfed maize transpiration was quite high with relatively small impacts on yield.For all other treatments whose results are not shown, the K cb act curve lays below the K cb curve during the periods when water stress occurred.in the vegetative growth and maturation stages (Kcb act < Kcb) because distribution of precipitation events along the maize season was favorable; thus, in 2012-2013, rainfed maize transpiration was quite high with relatively small impacts on yield.For all other treatments whose results are not shown, the Kcb act curve lays below the Kcb curve during the periods when water stress occurred.The soil evaporation coefficient Ke presents numerous peaks in correspondence with the numerous soil wetting events as depicted in Figure 4. Peaks are larger during the initial crop growth stages, when the soil was not yet covered by the crop and more energy was available at the ground surface for evaporation.Peaks progressively decreased when the crop was growing, thus increasingly shadowing the ground.Ke peaks resulted smaller during the mid-season when the crop was fully developed and largely shadowed the ground, thus when fc values were maximal (Table 4).Ke peaks slightly increased again during the late season when the crop senesced and fc decreased.The partition of ETc act into soil evaporation and plant transpiration followed the dynamics of Kcb act and Ke discussed above.Results for Es and Tc act for all treatments are presented in Table 7 along with other soil water balance components.For 2012-2013, when various irrigation deficits were applied, it is evident that the ratio Es/Tc act increased with the decrease of transpiration when irrigation amounts were smaller and water deficits led to less crop development and smaller ground cover.Results show that Es was 22% of ETc act in 2011-2012 and ranged from 20% to 25% of ETc act in 2012-2013.These ratios are in agreement with those reported by other authors who used mini-and micro-lysimeters [60][61][62].Much lower Es/ETc act ratios were reported for a drip irrigated maize with straw mulch cover [57], which supports controlling Es/ETc act when adopting no-till systems [4,5].
Results obtained for the water balance terms in Table 7 evidence the importance of the non-consumptive terms, DP and RO, which relates with the high precipitation observed and, in case of runoff, to the high depths of rainfall per event.DP represents between 18% and 33% of the seasonal precipitation and RO represents 20% to 28% of the seasonal precipitation.In case of the rainfed treatment the non-beneficial water use components (DP and RO) represented 44% of the season precipitation.Thus, despite a large fraction of rainfall water was not used by the crop, a decrease in   The soil evaporation coefficient K e presents numerous peaks in correspondence with the numerous soil wetting events as depicted in Figure 4. Peaks are larger during the initial crop growth stages, when the soil was not yet covered by the crop and more energy was available at the ground surface for evaporation.Peaks progressively decreased when the crop was growing, thus increasingly shadowing the ground.K e peaks resulted smaller during the mid-season when the crop was fully developed and largely shadowed the ground, thus when f c values were maximal (Table 4).K e peaks slightly increased again during the late season when the crop senesced and f c decreased.
The partition of ET c act into soil evaporation and plant transpiration followed the dynamics of K cb act and K e discussed above.Results for E s and T c act for all treatments are presented in Table 7 along with other soil water balance components.For 2012-2013, when various irrigation deficits were applied, it is evident that the ratio E s /T c act increased with the decrease of transpiration when irrigation amounts were smaller and water deficits led to less crop development and smaller ground cover.Results show that E s was 22% of ET c act in 2011-2012 and ranged from 20% to 25% of ET c act in 2012-2013.These ratios are in agreement with those reported by other authors who used mini-and micro-lysimeters [60][61][62].Much lower E s /ET c act ratios were reported for a drip irrigated maize with straw mulch cover [57], which supports controlling E s /ET c act when adopting no-till systems [4,5].Results obtained for the water balance terms in Table 7 evidence the importance of the non-consumptive terms, DP and RO, which relates with the high precipitation observed and, in case of runoff, to the high depths of rainfall per event.DP represents between 18% and 33% of the seasonal precipitation and RO represents 20% to 28% of the seasonal precipitation.In case of the rainfed treatment the non-beneficial water use components (DP and RO) represented 44% of the season precipitation.Thus, despite a large fraction of rainfall water was not used by the crop, a decrease in T c act and yield (Table 8) was observed.This yield decrease 35% of the Y a observed for the FI treatment.These results show that achieving high yields is influenced by supplemental irrigation.The different irrigation treatments were assessed in terms of TWU and water productivity (Table 8).Results show that the highest WP value was for the year 2011-2012 (2.17 kg¨m ´3) where yields were highest and TWU was not high; the lowest WP values were obtained for the DI treatments because crop growth was impacted and yields were low (Table 8).The rainfed treatment had the lowest yield and TWU but WP is similar to that of full irrigation.WP values are similar to those reported in previous studies in Portugal [28].WP ET observed values are generally superior to those reported in literature [13,14,[63][64][65].

Water-Yield Relations and Yield Predictions
Pairs of actual yield-transpiration data obtained from a set of maize experiments at Paysandú were used to assess maize water yield relations.These data were available for three maize seasons and various irrigation treatments of 2010-2011, 2011-2012 and 2012-2013, some of them not analyzed above but simulated with SIMDualKc to estimate T c act following a previous study with maize [28].
The actual yield and T c act data were used to derive the yield response factor of the model S1 (Equation ( 5)) and K y = 1.42 was obtained.This value is slightly larger than that obtained by Paredes et al. [28] and is in the range of values reported by Stewart et al. [12] and Kresovic et al. [26].A slightly large value was reported by Howell et al. [66] and Payero et al. [14], while a smaller K y was referred by Popova and Pereira [67].Results in literature allow to consider that K y = 1.42 is likely appropriate for further use with the model S1.All yield data pairs were also used to assess the accuracy of S2 model (Equation (6)) using the β values referred before (Section 2.2).
The "goodness-of-fit" indicators relative to both S1 and S2 (Table 9) allow assuming the appropriateness of parameters used for both models.The regression coefficients are close to 1.0 for both models but the R 2 values are relatively low due to a large variability of observed yields.RMSE = 1.83 t¨ha ´1 was obtained for S1, which represents 18% of the average observed yield (Y a ); RMSE was smaller for the S2 model, representing 14% of Y a .Results for EF, 0.59 and 0.74, respectively, for S1 and S2, indicate that the mean square error is smaller than the observed data variance.Notes: b 0 is the regression coefficient forced to the origin; R 2 is the determination coefficient of the ordinary least squares regression; RMSE is the root mean square error; NRMSE is the normalized root mean square error; ARE is the average relative error; and EF is the modeling efficiency.
Table 9 also includes the "goodness-of-fit" indicators when yield predictions used default parameters.For S1 the K y value tabulated by Doorenbos and Kassam [19] leads acceptable but worse "goodness-of-fit" indicators, namely higher errors of estimation.For S2 the original β v , β f and β m parameters proposed by Stewart et al. [12] lead to much worse indicators than the selected ones and the S1 model (Table 9).
Various crop growth and yield models have been tested for maize yield predictions.Stöckle et al. [52] using the CropSyst model reported RMSE representing a NRMSE of 8% to 21% of the average observed yields.Applications of the AquaCrop model show a variety of results, e.g., Heng et al. [68] reported RMSE ranging 0.65 to 1.57 t¨ha ´1 with NRMSE of 7% to 43%, and Ahmadi et al. [16] reported RMSE of 0.7 and 1.77 t¨ha ´1 with NRMSE of 7% and 18%.Constantin et al. [17] with AqYield and STICs models, reported NRMSE of 15% and 18% respectively.Using CERES-Maize, Lin et al. referred NRMSE near 11% [18], and Ma et al. [69] reported NRMSE ranging 6% to 12% and 5% to 7% when using, respectively, the DSAAT-CERES and the RZWQ-CERES models.Ko et al. [70] reported RMSE of 0.5 and 0.87 t¨ha ´1 but representing NRMSE of 7% and 11% of the average observed yields.Monzon et al. [71] reported RMSE of 1.54 to 2.22 t¨ha ´1 using CropSyst and CERES-Maize.Our results with S2 model (RMSE = 1.47 t¨ha ´1 and NMRSE of 14.3%) are therefore in the range of values reported in literature, nevertheless, our approach is much less demanding than models quoted.
Overall, results above (Table 9) show that both models S1 and S2 are appropriate for predicting maize yields in further applications when using the parameters validated in this study.The S1 model should be applied when only seasonal T c and T c act are available while the S2 model can be used when transpiration data are available for the three crop growth stages of vegetation development, flowering-pollination, and maturation.

Assessing Supplemental Irrigation and Rainfed Scenarios under Water Scarcity
As referred in Section 2.4, the SIMDualKc model was used with a 22-year series of weather data to estimate maize net irrigation requirements.The resulting NIR series range from 0 to 736 mm, respectively for the wettest and driest years (Figure 5).The frequently high NIR, above 370 mm, indicate that rainfed maize may be often at risk and supplemental irrigation may be required.Thus, for the dry and very dry years, whose NIR values have the probability of non-exceedance of, respectively, 20% and 5% (Figure 5), full and deficit irrigation were assessed.Main climatic data characterizing both years are shown in Figure 6, where it is apparent that the main cause of climatic variability is precipitation while ET o variation is relatively small.Seasonal results of the soil water balance for the considered irrigation scheduling strategies (Section 2.4) and both dry and very dry scenarios are presented in Table 10.Results show that the seasonal irrigation depth in the dry year was 71% to 84% of Tc act, respectively, for the Mod and Full irrigation scenarios; in the very dry year, it corresponded to 81 to 86% of Tc act for the same scenarios.Adopting the improved schedules summarized in Table 10, results allow concluding that: (a) runoff may be reduced to nearly 1% of seasonal rainfall in the dry and very dry seasons; (b) deep   Seasonal results of the soil water balance for the considered irrigation scheduling strategies (Section 2.4) and both dry and very dry scenarios are presented in Table 10.Results show that the seasonal irrigation depth in the dry year was 71% to 84% of Tc act, respectively, for the Mod and Full irrigation scenarios; in the very dry year, it corresponded to 81 to 86% of Tc act for the same scenarios.Adopting the improved schedules summarized in Table 10, results allow concluding that: (a) runoff may be reduced to nearly 1% of seasonal rainfall in the dry and very dry seasons; (b) deep   Seasonal results of the soil water balance for the considered irrigation scheduling strategies (Section 2.4) and both dry and very dry scenarios are presented in Table 10.Results show that the seasonal irrigation depth in the dry year was 71% to 84% of T c act , respectively, for the Mod and Full irrigation scenarios; in the very dry year, it corresponded to 81 to 86% of T c act for the same scenarios.Adopting the improved schedules summarized in Table 10, results allow concluding that: (a) runoff may be reduced to nearly 1% of seasonal rainfall in the dry and very dry seasons; (b) deep percolation may be reduced to 2%-5% of the season water use in the dry season and less than 1% in the very dry one; however, high DP is expected to occur when very large daily rainfall depths occur; (c) small RYD are achievable but the need to cease irrigations 20 days prior to harvesting aimed at decreasing grain moisture may lead to late water deficits causing yield decreases, particularly when maize varieties are highly sensitive to water stress; (d) computed WP and WP ET are similar to the best observed (Table 8), with highest values referring to the Full scenario in very dry conditions, which relates with achieving the highest yield, but results for the Mild deficit scenario are quite similar; and (e) WP and WP ET are higher in the very dry year because TWU is smaller than for the dry year as a decrease of rainfall may be compensated by an increase of irrigation.Results show that under dryness and/or drought conditions rainfed maize is not feasible due to very high yield losses (>91%, Table 10).This study evidences that improved schedules can lead to the best conjunctive use of rainfall and irrigation but this is only achievable when irrigation decisions are supported by a water balance model like SIMDualKc.
As previously referred, rainfed maize is commonly practiced in western Uruguay.To assess related yield consequences, following the approach by Popova et al. [6], the 22-year climatic data series were simulated for rainfed conditions and results were analyzed in terms of TWU, RYD, WP and WP ET , with the actual yield estimated with the model S2 (Equation ( 6)).The respective empirical probability curves for the same 22 years are presented in Figure 7.
Water 2016, 8, 309 17 of 22 percolation may be reduced to 2%-5% of the season water use in the dry season and less than 1% in the very dry one; however, high DP is expected to occur when very large daily rainfall depths occur; (c) small RYD are achievable but the need to cease irrigations 20 days prior to harvesting aimed at decreasing grain moisture may lead to late season water deficits causing yield decreases, particularly when maize varieties are highly sensitive to water stress; (d) computed WP and WPET are similar to the best observed (Table 8), with highest values referring to the Full scenario in very dry conditions, which relates with achieving the highest yield, but results for the Mild deficit scenario are quite similar; and (e) WP and WPET are higher in the very dry year because TWU is smaller than for the dry year as a decrease of rainfall may be compensated by an increase of irrigation.Results show that under dryness and/or drought conditions rainfed maize is not feasible due to very high yield losses (>91%, Table 10).This study evidences that improved schedules can lead to the best conjunctive use of rainfall and irrigation but this is only achievable when irrigation decisions are supported by a water balance model like SIMDualKc.
As previously referred, rainfed maize is commonly practiced in western Uruguay.To assess related yield consequences, following the approach by Popova et al. [6], the 22-year climatic data series were simulated for rainfed conditions and results were analyzed in terms of TWU, RYD, WP and WPET, with the actual yield estimated with the model S2 (Equation ( 6)).The respective empirical probability curves for the same 22 years are presented in Figure 7.  TWU during the maize season ranged from 2283 to 10585 m 3 (Figure 7a) reflecting the high variability of precipitation.Actual yields varied enormously (Figure 7b), from 0 to 15.8 t•ha −1 ; under extremely dry conditions (cf. Figure 6b), when an extreme water deficiency occurs during flowering, pollination and yield formation, grain yield could not be produced.Relative yield decreases varied much, from just 4% in wet years to 100% in very dry years (Figure 7c) when water deficits do not allow pollinations and yield formation.It is likely that for nearly 40% of the years RYD is too much high and economic returns are insufficient; however, appropriate economic studies are required in future.TWU during the maize season ranged from 2283 to 10,585 m 3 (Figure 7a) reflecting the high variability of precipitation.Actual yields varied enormously (Figure 7b), from 0 to 15.8 t¨ha ´1; under extremely dry conditions (cf. Figure 6b), when an extreme water deficiency occurs during flowering, pollination and yield formation, grain yield could not be produced.Relative yield decreases varied much, from just 4% in wet years to 100% in very dry years (Figure 7c) when water deficits do not allow pollinations and yield formation.It is likely that for nearly 40% of the years RYD is too much high and economic returns are insufficient; however, appropriate economic studies are required in future.WP and WP ET also vary enormously (Figure 7d) due to the variability of both rainfall (and TWU, Figure 7a) and yield (Figure 7b): WP ranged from 0 to 2.1 kg¨m ´3 and WP ET varied from 0 to 2.6 kg¨m ´3.

Conclusions
Field data relative to several irrigation schedules, including irrigation at different crop stages, were used to calibrate and validate the water balance and irrigation scheduling model SIMDualKc by minimizing the errors of estimate of the soil water content throughout the crop seasons.The corresponding "goodness-of-fit" indicators were very good when calibrated parameters were used and also quite good when using well selected default parameters.Simulations identified well the periods when the crop was water stressed.Overall, results show that the model is appropriate for further use in Uruguay, namely to support improved irrigation scheduling practices and management.The calibrated parameters are likely appropriate for being used as default ones in further applications; nevertheless, the parameters relative to soil evaporation and deep percolation need to be adjusted considering the soil characteristics of the locations where the model will be applied.
Results have shown that the ratio E s /T c act increased with the decrease of transpiration when water deficits limited crop development and led to reduced ground cover.The ratio E s /ET c act behaved similarly.Results for this ratio agree well with literature.These ratios may be reduced if mulched direct planting is adopted.Non-consumptive water balance terms DP and RO were quite high due to the local rainfall regime.
Both water-yield models S1 and S2 were successfully parameterized and their test for various sets of data provided quite good "goodness-of-fit" indicators.The staged S2 model provided better predictions than the model S1 because it considers the effects of water stress during the critical crop stages.Their errors were in the range of those reported in literature for the application of more demanding crop growth and yield models.Therefore, predictions using those models may be used in the future to estimate yields required for water productivity estimations.WP results were quite high, with best values when higher yields were obtained.
Several irrigation scheduling scenarios were assessed for dry and very dry years.Respective results show that a mild deficit irrigation is likely appropriate if considering precise deficit irrigation thresholds.Modeling provided good solutions for the conjunctive use of rainfall and irrigation, with high control of RO and DP, thus leading to more efficient water use.Computed performance indicators, mainly WP and WP ET , were similar to those computed when observed data were used.Results for rainfed maize have shown a great range of yield and water productivity variation, which likely lead to non-feasible rainfed production in about 40% of the years.However, a more appropriate assessment of the feasibility of rainfed and deficit irrigation solutions requires using field economic data and developing an adequate economic analysis in addition to indicators used herein.

Figure 1 .
Figure 1.Daily weather data relative to the maize crop seasons of 2011-2012 (a,b) and 2012-2013 (c,d): maximum and minimum temperatures, T max and T min (a,c); precipitation and reference evapotranspiration (ET o ) (b,d).

Figure 2 .
Figure 2. Flow-chart for the simplified approach combining the soil water balance model SIMDualKc with the Stewart's water-yield model.

Figure 2 .
Figure 2. Flow-chart for the simplified approach combining the soil water balance model SIMDualKc with the Stewart's water-yield model.

i
Full irrigation (Full), aimed at preventing water stress, with MAD = p.ii Mild deficit irrigation (Mild): MAD = 1.20 p for the initial period, MAD = 1.30p for the crop development and the late season periods, and MAD = 1.10 p during mid-season, which includes flowering and yield formation.iii Moderate deficit irrigation (Mod) with MAD = 1.30p for the initial and crop development periods, MAD = 1.20 p for the mid-season period, and MAD = 1.40 p for the late-season, after grain filling until harvesting.iv Rainfed.

Figure 3 .
Figure 3. Simulated vs. observed seasonal variation of soil water content (SWC, cm 3 •cm −3 ) when using calibrated (left) and default parameters (right): a deficit irrigation case for 2011-2012 used for calibration (a,b); FI in 2012-2013 (c,d) and rainfed in 2012-2013 (e,f), both used for validation.(θFC and θWP represent SWC at field capacity and wilting point, and θp is the SWC threshold for no stress; error bars represent the standard deviation of the SWC measurements).Also depicted precipitation and irrigation.

Figure 3 .
Figure 3. Simulated vs. observed seasonal variation of soil water content (SWC, cm 3 ¨cm ´3) when using calibrated (left) and default parameters (right): a deficit irrigation case for 2011-2012 used for calibration (a,b); FI in 2012-2013 (c,d) and rainfed in 2012-2013 (e,f), both used for validation.(θ FC and θ WP represent SWC at field capacity and wilting point, and θ p is the SWC threshold for no stress; error bars represent the standard deviation of the SWC measurements).Also depicted precipitation and irrigation.

Figure 4 .
Figure 4. Seasonal variation of the coefficients Kcb, Kcb act, Ke, including net irrigation and precipitation for the following experimental conditions in 2012-2013: (a) FI; (b) DIFLO; and (c) rainfed.

Figure 4 .
Figure 4. Seasonal variation of the coefficients K cb , K cb act , K e , including net irrigation and precipitation for the following experimental conditions in 2012-2013: (a) FI; (b) DI FLO ; and (c) rainfed.

Figure 5 .Figure 6 .
Figure 5. irrigation requirements for maize in Paysandú with identification of the dry and very dry climatic demand conditions.

Figure 5 .Figure 6 .
Figure 5. Net irrigation requirements for maize in Paysandú with identification of the dry and very dry climatic demand conditions.

Figure 6 .
Figure 6.Daily precipitation and reference evapotranspiration (ET o ) during the maize crop seasons for: (a) dry years; and (b) very dry years.

Figure 7 .
Figure 7. Probability curves characterizing rainfed maize through 22 years data relative to: (a) total water use (TWU); (b) estimated actual yield; (c) relative yield decrease (RYD); and (d) total and consumptive use water productivity.

Figure 7 .
Figure 7. Probability curves characterizing rainfed maize through 22 years data relative to: (a) total water use (TWU); (b) estimated actual yield; (c) relative yield decrease (RYD); and (d) total and consumptive use water productivity.

Table 1 .
Main soil physical and hydraulic properties of the experimental site, Paysandú.
Note: θFC and θWP are respectively the soil water content at field capacity and wilting point.

Table 1 .
Main soil physical and hydraulic properties of the experimental site, Paysandú.

Layer Depth (m) Particle Size (%) Soil Water Content (cm 3 ¨cm ´3)
Note: θ FC and θ WP are respectively the soil water content at field capacity and wilting point.

Table 2 .
Irrigation dates and depths for all treatments and both season.

Table 3 .
Crop growth stages dates and accumulated growth degree days (CGDD) for each experimental season.

Table 4 .
Crop height (h) and fraction of ground covered by the crop (f c ) at the main crop growth stages.

Table 5 .
Indicators of "goodness-of-fit" of SIMDualKc model applied to all treatments when using calibrated and default parameters.
Notes: b 0 is the regression coefficient forced to the origin; R 2 is the determination coefficient of the ordinary least squares regression; RMSE is the root mean square error; NRMSE is the normalized root mean square error; ARE is the average relative error; and EF is the modeling efficiency.

Table 6 .
Standard and calibrated basal crop coefficients (K cb ), depletion fractions for no stress (p), and parameters characterizing the soil evaporation layer, deep percolation and runoff.

Table 7 .
Simulated water use components for drip irrigated maize, Paysandú.
Notes: P-precipitation; I-net irrigation; ∆ASW-seasonal variation of available soil water; DP-deep percolation; RO-runoff; E s -soil evaporation; T c act -actual crop transpiration; ET c act -actual crop evapotranspiration; * irrigation performed to allow good crop emergence and establishment.

Table 8 .
Observed total water use (TWU), yield and water productivity.: TWU-total water use; ET c act -crop evapotranspiration; WP-water productivity; WP ET -consumptive use water productivity; in brackets the standard deviation relative to three repetitions per irrigation treatment. Notes

Table 9 .
Indicators of "goodness of fit" relative to estimating maize yield using the S1 and S2 models.

Table 10 .
Simulations of the soil water balance of diverse irrigation schedules for the dry and very dry years.∆ASW-seasonal variation of the available soil water; DP-deep percolation; TWUtotal water use; Tc and Tc act-maximum and actual crop transpiration; RYD-relative yield decrease; -estimated actual yield; WP-water productivity; WPET-consumptive use water productivity.
Figure 5. Net irrigation requirements for maize in Paysandú with identification of the dry and very dry climatic demand conditions.

Table 10 .
Simulations of the soil water balance of diverse irrigation schedules for the dry and very dry years.
Notes: RO-runoff; ∆ASW-seasonal variation of the available soil water; DP-deep percolation; TWUtotal water use; Tc and Tc act-maximum and actual crop transpiration; RYD-relative yield decrease; -estimated actual yield; WP-water productivity; WPET-consumptive use water productivity.

Table 10 .
Simulations of the soil water balance of diverse irrigation schedules for the dry and very dry years.
Notes: RO-runoff; ∆ASW-seasonal variation of the available soil water; DP-deep percolation; TWU-total water use; T c and T c act -maximum and actual crop transpiration; RYD-relative yield decrease; Ŷa -estimated actual yield; WP-water productivity; WP ET -consumptive use water productivity.