# As simple as possible but not simpler?: the case of irrigation modeling at catchment scale in southwestern France

Myriam Soutif-Bellenger, Guillaume Thirel, Olivier Therond, Jean Villerd

► To cite this version:

Myriam Soutif-Bellenger, Guillaume Thirel, Olivier Therond, Jean Villerd. As simple as possible but not simpler?: the case of irrigation modeling at catchment scale in southwestern France. Irrigation Science, 2023, 10.1007/s00271-023-00846-x. hal-04004071

HAL Id: hal-04004071

https://hal.inrae.fr/hal-04004071

Submitted on 18 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# As simple as possible but not simpler?: the case of irrigation modeling at catchment scale in southwestern France

Myriam Soutif-Bellenger<fnref n="1" />,<fnref n="2" /> ID · Guillaume Thirel<fnref n="2" /> ID · Olivier Therond<fnref n="3" /> ID · Jean Villerd<fnref n="3" /> ID

Received: 8 April 2022 / Accepted: 3 January 2023 / Published online: 21 January 2023

© The Author(s) 2023

## Abstract

The estimation of irrigation water requirements (IWR) amount and timing is crucial for designing water management strategies at the regional scale. Irrigation requirements can be estimated with different types of models: very complex and detailed crop models, agent-based models, or simplified modeling approaches. Because simplified approaches are often preferred, in this study, we evaluate the consequences of using simplified approaches for IWR assessment at a catchment scale and the consequences of various modeling choices, providing information on the uncertainties. To this end, different simple modeling approaches based on the CropWat model are compared with an agent-based approach (MAELIA), which serves as a benchmark. To assess simulations in detail, partial variance is calculated for several indicators characterizing daily simulated irrigation. Our sensitivity analysis, applied over a sub-catchment of the Aveyron River (southwestern France), shows a high variability in simulations produced by CropWat between the modeling assumptions tested, principally explained by the rules for irrigation triggering and the quantification of daily irrigation. The analysis also shows that several simplified approaches are able to reproduce the irrigation simulated by the high-accuracy MAELIA model, but not necessarily corresponding to an optimal irrigation scheme. Hence, this study confirms the possibility of assessing daily irrigation with simplified approaches, but warns about high modeling uncertainties, reflecting uncertainty in effective irrigation practices. This uncertainty can be taken into account by water managers and modelers through the combination of a set of irrigation models.

## Introduction

Water withdrawals for irrigation water requirements (IWR) have a huge impact on low flows in water-stressed agricultural catchments (Martin et al. 2016). Their impact may increase with climate change as a consequence of the decrease in water resources in summer and the increase in crop water requirements (Wanders and Wada 2015), as is expected in France (Collet et al. 2015). Hence, IWR assessment is essential for managing water resources in water-stressed agricultural catchments. It might become essential to take into account irrigation withdrawals processes and their evolution while performing hydrological modeling

and considering future climate change. Moreover, IWR assessment may help to evaluate the relevance of adaptation strategies to climate change scenarios in association with hydrological projections. The IPCC (2014) recommended planning water management at the catchment scale. In France, a new planning and dialog tool was created, coined "PTGE" (Projets de Territoire pour la Gestion de l'Eau, i.e., Territory Project for Water Management; MTES and MAA 2019), to build water management strategies at the local scale. To support the design of such planning or strategies, many hydrological simulations have been coupled with estimations of irrigation withdrawals based on IWR assessment (Collet et al. 2013; Dehghanipour et al. 2020; Gorguner and Kavvas 2020; Kolokytha and Malamataris 2020; Wanders and Wada 2015). Likewise, several studies estimated crop water stress induced by limited water resources for irrigation (Collet et al. 2013; Elliott et al. 2014).

IWR in a catchment can be estimated from observation data. However, observation data are often lacking, leading to a preferential use of automatic irrigation algorithms for irrigation assessment (Uniyal and Dietrich 2019; Wriedt et al. 2009), i.e., algorithms based on soil–crop water balance

Communicated by Erko Stackebrandt.

<figure><img src="image_2.png" /></figure>

<footnote>Myriam Soutif-Bellenger
myriam.soutif-bellenger@inrae.fr</footnote>

<footnote n="1">¹ AgroParisTech, 75005 Paris, France</footnote>

<footnote n="2">² Université Paris-Saclay, INRAE, UR HYCAR,
92160 Antony, France</footnote>

<footnote n="3">³ Université de Lorraine, INRAE, LAE, 68000 Colmar, France</footnote>
models. Specifically, irrigation algorithms simulate auto-
matic triggering and nominal irrigation depths to quantify
irrigations. This enables, for example, long-term IWR pro-
jections in the context of climate change. Most often, irri-
gation is triggered according to decision rules based on soil
moisture deficit or crop water stress thresholds. Irrigation
amounts can be fixed (Bouras et al. 2019; Rouhi Rad et al.
2020), or can be calculated to fill the soil reservoir totally
(Collet et al. 2013; Hori et al. 2008) or partially (Funes et al.
2021; Smith et al. 2012). To estimate soil moisture or crop
water stress, a soil–crop water balance model is often used.
The development and the use of sophisticated crop models
have been largely investigated (Di Paola et al. 2016). The
spatio-temporal distribution of farmers’ practices strongly
determines irrigation dynamics and amounts (Bergez et al.
2012; Zaccaria et al. 2013). McInerney et al. (2018) explored
the impacts of different spatio-temporal distribution methods
of observed irrigation among hydrological response units
(HRUs) of a catchment and showed that the choice of dis-
tribution method might have an important impact on flows.
An approach was developed to estimate IWRs at the regional
scale based on a high-accuracy reproduction of farmers’
practices, taking into account spatio-temporal variability.
This approach is agent-based and is called the “MAELIA”
(Modelling of socio-Agro-Ecological system for Landscape
Integrated Assessment) platform (Allain et al. 2018; Martin
et al. 2016; Therond et al. 2014). In this last example, farm-
ers’ practices are represented on an individual basis, using
decision algorithms and taking into account field size and
spatial allocation, irrigation equipment, and working time
constraints. Nevertheless, the complexity of the spatio-tem-
poral distribution of farmers’ practices is often overlooked in
efforts aiming to develop operational methods of irrigation
quantification and to deal with the lack of data. However,
the spatio-temporal variability of irrigation practices might
be high and might have strong impacts on hydrology, and
simplifying assumptions might lead to biases that should
be evaluated.

Few previous studies performed a sensitivity analysis
of IWR assessment methods. Multsch et al. (2015) quanti-
fied wheat IWR with different potential evapotranspiration
formulas and crop coefficients sets and analyzed the vari-
ability obtained. Wada et al. (2013) and Wan et al. (2018)
developed a hydrological multi-model approach to estimate
global IWR. However, to our knowledge, the sensitivity of
IWR assessment to different automatic irrigation rules at
a regional scale has not been quantified. Moreover, these
comparisons between IWR simulations are usually based
on annual or monthly scales, while regulatory constraints
and objectives are usually defined on the basis of shorter
time steps (e.g., day(s)) and on small catchment scale as

is the case, for example, in France (Mazzega et al. 2014;
MEDDTL 2011).

The present study aims to fill these gaps in IWR assess-
ment by comparing the impact of different automatic irri-
gation rules with other sources of potential variability on
which many researchers focus their efforts, such as crop
evapotranspiration estimation, or root growth estimation. In
an operational perspective and for integrated hydrological
modeling, it is important to identify the strongest source
of uncertainty in IWR assessment. Two types of mod-
eling approaches involving different levels of simplifica-
tion of farmers’ practices are compared with the MAELIA
benchmark in a French southwestern water basin with a
significant water deficit. MAELIA has already been cali-
brated and validated on this study area in a previous work
(Martin et al. 2016; Murgue et al. 2016). Different versions
of these approaches are explored and compared with the
MAELIA outputs. To this end, indicators characterizing
irrigation hydrographs at the daily time step are used. The
impact of modeling hypotheses was quantified with variance
decomposition.

The following research questions are investigated in the
present study:

1. What is the sensitivity of daily IWR assessment to simplifying modeling assumptions at different spatial scales? How should this sensitivity be taken into account by hydrological modelers and water managers?

2. Is it possible to adequately reproduce simulations of a high-accuracy agent-based model with simplified approaches?

Study area and data

Study area

The study focuses on an 840-km² downstream portion of
the Aveyron River basin, a tributary of the Garonne River,
located in southwestern France (Fig. 1). This area is mostly
agricultural, with 58,500 ha of cropland, 22,000 ha of for-
ests and semi-natural areas, and only 3500 ha of artificial-
ized soils. Irrigated crops comprise 8000 ha. Maize is by far
the main irrigated crop, followed by wheat and orchards.
Murgue et al. (2016) estimated that mono-cropping of grain
maize in alluvial soils reached an annual average irrigation
amount of 255–305 mm during the 2003–2007 period in
the study area.

This area is highly water-stressed and agricultural with-
d withdrawal restrictions are common in summer. These restric-
tions are aimed at avoiding ecological degradation of the
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Study area location (a) and land use (b)</figcaption></figure>

aquatic environment and ensuring water is available for domestic and industrial use.

The climate is temperate with an annual mean air temperature of 13 °C and mean annual precipitation of approximately 750 mm. However, summers are dry and hot, when maize water needs are high, leading to a high irrigation dependency. Monthly precipitations are less than 50 mm in July and August, while monthly evapotranspiration is more than 115 mm during the same period.

## Data

The following data were available:

* Crop rotation sequences and irrigated surfaces from the French Land Parcel Identification System (LPIS) database created in 2006, based on European Union Common Agricultural Policy declarations of cultivated plots;

* Soil characteristics from the Soil Geographical Data Base of France (SGDBF) (INRA 2018);

* Daily data of past climate conditions (precipitation and air temperature) from the SAFRAN reanalysis (Vidal et al. 2010).

Data were previously completed and adjusted through a survey among local farmers and stakeholders (see Murgue et al. 2015).

The major irrigated crops in the study area are cereals (other than maize), maize (six cultivars from very early to very late), maize seeds, maize silage, rapeseed, peas, soybean, sunflower, orchards, and grassland.

# Modeling protocol

## Irrigation modeling at regional scale

In the present study, we compare two approaches of varying complexity:

* The MAELIA platform, combining the soil–crop water balance model AqYield (Constantin et al. 2015) at plot scale and a high-accuracy agent-based automatic irrigation modeling. This modeling approach is called “MAELIA” (M);

* The soil–crop water balance model CropWat (Smith 1992) combined with two simpler regional automatic irrigation approaches called “Conceptual” (C) and “Semi-plot” (S).

These different approaches are described below and summarized in Fig. 2. Details about the soil–crop water balance models, CropWat and AqYield, are given in Appendix “Soil–crop water balance models”. Their parameter calibration procedure is detailed in Appendix “Parameter calibration”. The beginning of the simulation period is January 2008 and the end is December 2014. The year 2007 is used for initialization.

## The MAELIA approach benchmark

MAELIA (Therond et al. 2014) is a high-accuracy platform modeling explicitly farmers’ practices at a daily time step in each plot (n = 15,224) individually for each farm (n = 1143) described in the LPIS database of the study area. It provides
**Fig. 2** The three modeling approaches used to assess irrigation water requirements (IWRs). The MAELIA (M) irrigation approach distinguishes 232 pedo-climatic zones and each individual plot and farm. The AqYield soil–crop water balance model is coupled to an agent-based model in MAELIA to assess IWR. A Conceptual approach (C) is used to simulate daily irrigation amounts on the basis of data aggregated on 200 pedo-climatic zones crossed with existing crops based on CropWat outputs. A Semi-plot approach (S) in which the simulation units of the C approach are surface-eqally divided into irrigation water turn groups irrigated successively. At simulation unit scale, the S approach leads to irrigation events on singular days followed by several days without irrigation represented by peaks in the figure. AqYield outputs are used to calibrate the CropWat parameters

<figure><img src="image_4.png" /></figure>

Complexity

a plot-scale modeling of irrigation practices. Technical operation durations are taken into account to constrain the number of farmers’ actions within a single day. Irrigation management strategies described through IF-THEN decision rules are defined for different combinations of crops, soils, irrigation equipment, farm types, etc. Each strategy defines a possible period for irrigation, the dose for each irrigation (mm), the conditions triggering irrigation (e.g., soil moisture and past and future rainfall thresholds, crop water stress), and the minimal temporal interval between irrigation applications in the same field (hereafter called “water turn”). Details regarding these decision rules are provided by Murgue et al. (2014). The irrigation strategies in our study area stem from a farm survey performed and described by Murgue et al. (2015).

The soil–crop model AqYield runs on each plot (Fig. 2), characterized by crop, technical management, soil properties, and climate. For the study area, 51 climatic zones (CZs) and 14 soils were differentiated. Crossing soils and CZs generate 232 pedo-climatic areas.

Unfortunately, irrigation data are very coarse in space, time span and time resolution, and their exactitude is highly questionable in France and particularly in the South-west

region, where this information is highly sensitive. The only existing information is irrigation declared by farmers to the French Water Authority at an annual time step. Moreover, actual irrigation withdrawals cannot be directly compared with IWR, because withdrawals can be limited by water availability and restrictions. A comparison between simulated irrigation withdrawals with MAELIA, taking into account water availability and restriction rules, and annual individual farmer reports to French water authorities showed a very good fit (Martin et al. 2016). The temporal distribution of irrigation over weeks was also assessed through local expert interviews (Murgue et al. 2016). Therefore, we consider in our work that the MAELIA irrigation demand estimation, without taking into account water availability and restriction rules, is probably the best assessment of daily IWR that one can afford for this region given the available data and it will represent our reference in this study.

## Conceptual approach

In many studies (Collet et al. 2015; Smith et al. 2012), regional irrigation modeling consists of bringing the amount of water (Irr, mm) to fill the soil reservoir, to maintain the
soil water deficit (Dr, mm) to a prefixed threshold ($\theta$, mm) for each simulation unit, as expressed in Eq. (1). In this configuration, irrigation is triggered when Dr exceeds $\theta$. A similar approach is to bring the amount of water needed to compensate for the lack of soil water to evapotranspirate at the crop maximal evapotranspiration (CET, mm day$^{-1}$) level, maintaining water at the readily available water (RAW, mm) level equivalent to a $\theta$ value (Funes et al. 2021):

$$ \begin{cases} \text{if } Dr(d) < \theta(d), Irr(d) = 0 \\ \text{else, } Irr(d) = Dr(d) - \theta(d) \end{cases} \quad (1) $$

In this simplified approach, individual plots are not explicitly represented. Simulation units to calculate daily irrigation are defined for homogeneous crop, soil, and climate types (Fig. 2). In this study, 200 pedo-climatic zones were defined, aggregating MAELIA spatial units characterized by the same climate and water content characteristics.

## Semi-plot approach

We developed an intermediate approach in this study, called the “Semi-plot” approach. It aims at reproducing a more realistic modeling of irrigation over plots, but in a more simplified way than the reference (MAELIA). In this approach, pedo-climatic zones are the same as in the Conceptual approach but the crop simulation units of the Conceptual approach are divided into irrigation water turn groups of equal area and are irrigated successively (Fig. 2). For example, for a 5-day water turn, a simulation unit would be divided in 5 groups, and the first group could be irrigated on days 1, 6, 11, etc., the second group could be irrigated on days 2, 7, 12, etc., and so on. The number of groups for each crop type was defined based on the number of water

turns simulated by the reference experiment. This approach simulates more explicitly the dynamics of irrigation at the plot scale as each simulation unit can receive a high amount of water on a particular day followed by several days without irrigation (Fig. 2).

## Indicators characterizing irrigation

While most studies only consider monthly or even annual IWR amounts, this study explores numerous indicators to characterize the temporal variability of withdrawals under different simulations. Indeed, water crises can occur in short time intervals, and a monthly estimation of IWR might not be sufficient to identify them (Mazzega et al. 2014). Indicators were selected to characterize seven features of irrigation:

* irrigation volumes ($I, I_m$),

* irrigation frequency ($N, N_m$),

* period of inter-annual mean highest irrigation (Max_regime, Max_regime_date),

* mean annual highest irrigation (Max_10, Max),

* irrigation daily variability ($Q$, dispersion, $Var2$),

* accuracy compared with reference irrigation simulation (KGE, KGE_10),

* crop water stress ($ET_S, ET_S_m$).

Details on these indicators are provided in Appendix “Indicator calculation”. In particular, the $Var2$ indicator was developed to compare narrow temporal variabilities of irrigation.

<caption>**Table 1** Mean, minimal, and maximal values of annual indicators (“Indicators characterizing irrigation”) for Conceptual, Semi-plot and reference experiments (as described in “Irrigation modeling at</caption>

regional scale” and “Sensitivity analysis experiments”) for lumped outputs at the scale of the study zone

<table><thead><tr><th rowspan="2">Indicator</th><th colspan="3">Conceptual</th><th colspan="3">Semi-plot</th><th rowspan="2">Reference</th></tr><tr><th>Mean</th><th>Min</th><th>Max</th><th>Mean</th><th>Min</th><th>Max</th></tr></thead><tbody><tr><td>I (m³)</td><td>1.40E+07</td><td>6.08E+06</td><td>2.48E+07</td><td>1.54E+07</td><td>4.83E+06</td><td>2.99E+07</td><td>1.41E+07</td></tr><tr><td>N (days)</td><td>120</td><td>79</td><td>188</td><td>148</td><td>83</td><td>269</td><td>154</td></tr><tr><td>Max_regime (m³)</td><td>1.05E+05</td><td>7.20E+04</td><td>1.47E+05</td><td>1.10E+05</td><td>5.40E+04</td><td>1.72E+05</td><td>1.16E+05</td></tr><tr><td>Max_regime_date (DOY)</td><td>197 (16 July)</td><td>139 (19 May)</td><td>201 (20 July)</td><td>194 (13 July)</td><td>140 (20 May)</td><td>230 (18 August)</td><td>201 (20 July)</td></tr><tr><td>Max (m³)</td><td>5.10E+05</td><td>2.81E+05</td><td>1.41E+06</td><td>4.21E+05</td><td>1.91E+05</td><td>8.08E+05</td><td>3.29E+05</td></tr><tr><td>Max_10 (m³)</td><td>3.10E+05</td><td>1.94E+05</td><td>4.31E+05</td><td>3.22E+05</td><td>1.46E+05</td><td>5.07E+05</td><td>1.47E+05</td></tr><tr><td>Q_dispersion (-)</td><td>3.7</td><td>2.4</td><td>6.4</td><td>4.3</td><td>2.2</td><td>13.7</td><td>5.1</td></tr><tr><td>Var2 (-)</td><td>0.0021</td><td>0.0019</td><td>0.0027</td><td>0.0010</td><td>0.0006</td><td>0.0016</td><td>0.0009</td></tr><tr><td>ET_S (-)</td><td>0.97</td><td>0.87</td><td>1.00</td><td>0.97</td><td>0.83</td><td>1.00</td><td>0.93</td></tr><tr><td>KGE (-)</td><td>0.57</td><td>-0.1</td><td>0.81</td><td>0.56</td><td>-0.26</td><td>0.91</td><td>-</td></tr><tr><td>KGE_10 (-)</td><td>0.70</td><td>0.17</td><td>0.95</td><td>0.62</td><td>-0.21</td><td>0.96</td><td>-</td></tr></tbody></table>
## Sensitivity analysis experiments

Balanced simulation plans were developed for both the Conceptual and Semi-plot approaches. They aim at studying the impact of modeling simplifications and at identifying sensitive parameters for each soil–crop water balance and regional irrigation modeling chain. The chosen variation factors are:

* Soil water deficit threshold θ. Many studies fix θ at a hypothetical RAW level (Neilsen et al. 2018; Rinaudo et al. 2013), corresponding to the value optimizing the amount of water brought to crops. Some studies fix various values of θ, e.g., at a value of 0 mm (Collet 2013), at a threshold of 50% of total available water (TAW, mm; Bouras et al. 2019), or at a threshold of 0.8×RAW + 0.2×TAW (Smith et al. 2012). Moreover, deficit irrigation amounts (θ > RAW) are also possible, particularly for crops resistant to water stress. Because

irrigation is calculated at a daily time step, contrary to Collet et al. (2013), who considered a 10-day time step, we consider that θ cannot be equal to 0 in this modeling configuration. We also include a value higher than RAW to explore a deficit irrigation hypothesis. Accordingly, we explore a range of values for θ from 0.25×RAW to RAW + 0.25×(TAW – RAW);

* Irrigation amounts (IA). For Semi-plot experiments, different quantification methods of irrigation in simulation units are tested: one similar to the Conceptual approach simulation units (Dr – θ), one with a fixed amount for each crop (FixMeanC), as done by Bouras et al. (2019) or Rouhi Rad et al. (2020), and one to fill the soil water reservoir (Dr), as done by Hori et al. (2008);

* Seasonal variation of depletion factor *p_var*. Simulations were made with variable *p* or constant *p* to evaluate the importance of taking into account this seasonal variability (see Appendix “The CropWat model”);

<figure><img src="image_5.png" /><figcaption>**Fig. 3** Daily regimes of irrigation obtained at the study zone scale from the reference experiments (purple) and from the ensemble, the mean and the 80% and 90% envelopes of all C and S experiments</figcaption></figure>

(blue), as described in “Irrigation modeling at regional scale” and “Sensitivity analysis experiments” (color figure online)

<figure><img src="image_6.png" /><figcaption>**Fig. 4** Partial variance explained by the factors of annual indicators. On the left, the Conceptual (C) experiments, on the right, the Semi-plot (S) experiments. A higher partial variance indicates that the variability of the indicator is more affected by this factor</figcaption></figure>
<table><caption>**Table 2** Characteristics of the “best experiments” compared to the reference (MAELIA) for Conceptual and Semi-plot approaches</caption><thead><tr><th colspan="2">Conceptual (96 experiments)</th><th colspan="3">Semi-plot (256 experiments)</th></tr><tr><th colspan="2"></th><th colspan="3">Number of experiments: 7</th></tr><tr><th></th><th>Factor</th><th>Modality</th><th>Number</th></tr></thead><tbody><tr><td rowspan="4">KGE > 0.9</td><td rowspan="4">Number of experiments: 0</td><td>θ</td><td>varCrop</td><td>7</td></tr><tr><td>IA</td><td>varCrop</td><td>7</td></tr><tr><td>IP</td><td>0</td><td>3</td></tr><tr><td>CET</td><td>1</td><td>4</td></tr><tr><td rowspan="10">KGE > 0.8</td><td rowspan="10">Number of experiments: 1</td><td colspan="3">Number of experiments: 58</td></tr><tr><th>Factor</th><th>Modality</th><th>Number</th></tr><tr><td>θ</td><td>varCrop</td><td>1</td></tr><tr><td>IP</td><td>0</td><td>1</td></tr><tr><td>CET</td><td>1</td><td>1</td></tr><tr><td>θ</td><td>0.25×RAW</td><td>8</td></tr><tr><td>IA</td><td>0.5×RAW</td><td>16</td></tr><tr><td>IP</td><td>0.75×RAW</td><td>10</td></tr><tr><td>CET</td><td>RAW</td><td>8</td></tr><tr><td>θ</td><td>varCrop</td><td>16</td></tr><tr><td>IA</td><td>Dr-θ</td><td>21</td></tr><tr><td>IP</td><td>FixMeanC</td><td>21</td></tr><tr><td>CET</td><td>varCrop</td><td>16</td></tr><tr><td>θ</td><td>0</td><td>28</td></tr><tr><td>IP</td><td>1</td><td>30</td></tr><tr><td>CET</td><td>0</td><td>19</td></tr><tr><td></td><td></td><td>1</td><td>39</td></tr><tr><td rowspan="10">KGE > 0.7</td><td rowspan="10">Number of experiments: 36</td><td colspan="3">Number of experiments: 100</td></tr><tr><th>Factor</th><th>Modality</th><th>Number</th></tr><tr><td>θ</td><td>0.5×RAW</td><td>8</td></tr><tr><td>IP</td><td>0.75×RAW</td><td>8</td></tr><tr><td>IP</td><td>RAW</td><td>8</td></tr><tr><td>CET</td><td>varCrop</td><td>16</td></tr><tr><td>IP</td><td>0</td><td>18</td></tr><tr><td>CET</td><td>1</td><td>18</td></tr><tr><td>IP</td><td>0</td><td>12</td></tr><tr><td>CET</td><td>1</td><td>24</td></tr><tr><td>θ</td><td>0.25×RAW</td><td>8</td></tr><tr><td>IA</td><td>0.5×RAW</td><td>17</td></tr><tr><td>IP</td><td>0.75×RAW</td><td>32</td></tr><tr><td>CET</td><td>RAW</td><td>20</td></tr><tr><td>IP</td><td>RAW + 0.25 (TAW-RAW)</td><td>7</td></tr><tr><td>θ</td><td>varCrop</td><td>16</td></tr><tr><td>IP</td><td>Dr-θ</td><td>41</td></tr><tr><td>CET</td><td>Dr</td><td>3</td></tr><tr><td>IP</td><td>FixMeanC</td><td>40</td></tr><tr><td>CET</td><td>varCrop</td><td>16</td></tr><tr><td>IP</td><td>0</td><td>48</td></tr><tr><td>CET</td><td>1</td><td>52</td></tr><tr><td>IP</td><td>0</td><td>45</td></tr><tr><td>CET</td><td>1</td><td>55</td></tr></tbody></table>

For three thresholds of performant KGE (best potential performance for KGE = 1), modalities of corresponding experiments are indicated, aimed at identifying modalities that clearly result in an improvement of KGE

* Root growth z_var. Simulations were made with increasing rooting depth between initial and maximal rooting depth or constant rooting depth equal to maximal rooting depth (see Appendix “The CropWat model”);

* Irrigation period (IP). MAELIA defines precise irrigation periods for each crop. However, in the context of

inter-annual climate variability increasing with climate change, it seems important to let the model calculate irrigation potential needs outside the usual irrigation periods. For this, we simulated experiments following irrigation periods for each crop defined in the reference
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Partial variance of monthly indicators related to different factors. For each month, partial variances are cumulated, and show the monthly variation in the contribution of each factor</figcaption></figure>

<table><caption>**Table 3** Mean, minimal, and maximal experiment values of mean CZ value for each indicator (“Indicators characterizing irrigation”) for Conceptual, Semi-plot and reference experiments (as described in “Irrigation modeling at regional scale” and “Sensitivity analysis experiments”)</caption><thead><tr><th rowspan="2">Indicator</th><th colspan="3">Conceptual</th><th colspan="3">Semi-plot</th><th rowspan="2">Reference</th></tr><tr><th>Mean</th><th>Min</th><th>Max</th><th>Mean</th><th>Min</th><th>Max</th></tr></thead><tbody><tr><td>I (m³)</td><td>2.98E+05</td><td>1.29E+05</td><td>5.28E+05</td><td>3.27E+05</td><td>1.03E+05</td><td>6.37E+05</td><td>3.01E+05</td></tr><tr><td>N (days)</td><td>75</td><td>47</td><td>127</td><td>85</td><td>33</td><td>165</td><td>76</td></tr><tr><td>Max_regime (m³)</td><td>4.53E+03</td><td>3.14E+03</td><td>6.28E+03</td><td>4.90E+03</td><td>2.39E+03</td><td>7.40E+03</td><td>4.86E+03</td></tr><tr><td>Max_regime_date (DOY)</td><td>198 (17 July)</td><td>143 (23 May)</td><td>218 (6 August)</td><td>196 (15 July)</td><td>142 (22 May)</td><td>222 (10 August)</td><td>197 (16 July)</td></tr><tr><td>Max_10 (m³)</td><td>6.71E+03</td><td>4.31E+03</td><td>9.43E+03</td><td>7.29E+03</td><td>3.27E+03</td><td>1.29E+04</td><td>6.20E+03</td></tr><tr><td>Max (m³)</td><td>1.12E+04</td><td>5.83E+03</td><td>3.01E+04</td><td>9.89E+03</td><td>4.00E+03</td><td>2.38E+04</td><td>8.80E+03</td></tr><tr><td>Q_dispersion (-)</td><td>1.47</td><td>1.09</td><td>1.87</td><td>2.21</td><td>1.26</td><td>3.68</td><td>1.45</td></tr><tr><td>Var2 (-)</td><td>0.0024</td><td>0.0021</td><td>0.0031</td><td>0.0018</td><td>0.0007</td><td>0.0031</td><td>0.0027</td></tr><tr><td>KGE_10</td><td>0.59</td><td>-0.03</td><td>0.86</td><td>0.47</td><td>-0.36</td><td>0.87</td><td>-</td></tr><tr><td>KGE</td><td>0.38</td><td>-0.39</td><td>0.64</td><td>0.28</td><td>-0.48</td><td>0.68</td><td>-</td></tr></tbody></table>

and compared them with experiments able to trigger irrigation at any time during the crop cycle;

* Crop maximal evapotranspiration calculation (CET) (see Appendix “The CropWat model”). The impact of the evapotranspiration calculation method has already been studied and quantified (Multsch et al. 2015). To identify the main sensitivity sources of regional irrigation uncertainty including evapotranspiration, we integrate experiments into our modeling scheme with a classic evaluation

of CET and with experiments considering CET equal to PET-PM.

The modalities for each factor are summarized in Appendix “Varying factors in C and S experiments”.

Some complementary experiments were conducted to evaluate the added value of defining different irrigation rules ($\theta$ and IA) between crops. For each crop, for each crossing of $p$-var, $z$-var, CET and IP hypotheses conditions, $\theta$ (for C experiments), and $\theta$ and IA (for S experiments), the

<figure><img src="image_8.png" /></figure>
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Comparison of indicator values at the CZ scale between the reference (MAELIA) and the best experiment. This experiment corresponds to a Semi-plot approach, with irrigation rules defined specifically for each crop, no irrigation period delimitation, and root growth as well as seasonal variation of *p* and *Kc* taken into account (modalities varCrop, 0, 1, 1, 1 for θ+IA, IP, z_var, p_var, CET, respectively)</figcaption></figure>

best-performing experiments in reproducing the reference experiment outputs at crop level were selected using the KGE-10 criterion (see “Indicators characterizing irrigation” and Appendix “Indicator calculation”). The resulting modality of θ and IA is called “varCrop.”

The variance decomposition procedure we used to estimate the sensitivity of indicators to variation factors is detailed in Appendix “Variance decomposition”.

## Results

In this section, we assess the irrigation obtained from the reference, Conceptual and Semi-plot experiments at the basin scale as well as the factors explaining the variability. Then, we analyze the best experiments in terms of KGE and assess the experiments at the CZ scale.

### Sensitivity of simulation to variation factors at the scale of the total study zone

#### General analysis of indicators

According to the results from the reference experiment, irrigation can start in April, is maximal in July (mean value of 6 million m³), and ends in October. The results of the reference, Conceptual (C) and Semi-plot (S) simulations are summarized in Table 1, presenting values of annual indicators aggregated for the total study zone. With higher mean values of *Q*, the S experiments show mean higher dephasing between daily values of high and low irrigation. However, with higher values of Var2, daily temporal variations of irrigation are globally higher in the C experiments. For almost all factors, the C results are within the envelope of the S results, which might be caused by the larger number of factors (IA) explored in the S experiments. Reference values of indicators are contained within the range of the C and S results except for the higher values reached by Max_10 and Var_2, which are not included within the range explored in the C experiments, showing that C experiments have a narrow higher daily variability than reference and S experiments (Appendix “Annual irrigation hydrographs”). The inclusion of the reference experiment in the CropWat-based experiments is confirmed in Fig. 3, showing that the reference regime of irrigation is contained within the simulations produced by our experiments.

### Distribution of the total variance in indicators between factors

The variance of annual indicators explained by different factors is plotted in Fig. 4. For the Semi-plot experiments, the six modalities of θ (0.25×RAW, ..., varCrop) and the four modalities of IA (Dr, ..., varCrop) were merged into one factor θ+IA of 16 modalities to avoid an unbalanced experimental plan, because the varCrop modality for θ was run only with the varCrop modality of IA.
<figure><img src="image_10.png" /><figcaption>**Fig. 7** Example of Kc (-) and p (-) curves for a set of parameters: crop coefficients (Kc<sub>ini</sub>, Kc<sub>mid</sub> and Kc<sub>end</sub>), depletion factor parameters (p<sub>ini</sub>, p<sub>mid</sub> and p<sub>end</sub>), and length of growth stages (L<sub>ini</sub>, L<sub>dev</sub>, L<sub>mid</sub> and L<sub>end</sub>)</figcaption></figure>

First, it is clear that the automatic irrigation parameters ($\theta$ for C experiments and $\theta$+IA for S experiments) override the effects of the other factors for most indicators. CET is also very impactful for both the Conceptual and Semi-plot experiments, particularly on indicators *I*, *N*, Max_regime_date, Var2, and ET_S, while IP is also very impactful on indicators *N* and *Q_dispersion*. Factors *p_var* and *z_var* are systematically not impactful, for each indicator of the Conceptual and Semi-plot experiments. Second, we note that the interaction effects might be high. The interaction between automatic irrigation parameters ($\theta$ and $\theta$+IA) and CET has a strong impact on the Max_regime_date and on KGE and KGE_10. Impacting factors on indicators are often the same between the C and S experiments, except for daily irrigation variability indicators, i.e., Var2, *Q_dispersion*. This difference between the C and S experiments is explained by the integration in the S experiments of several modalities of IA, leading to large variations between experiments for these indicators (Appendix “*Impact of IA modalities on Var2 and Q_dispersion*”).

The evaluation of the impact of each modality on each indicator is detailed in Appendix “*Impact of factor modalities on indicators*”. It is a necessary additional step after variance decomposition to select accurate modalities for modeling purposes and to identify the strengths and weaknesses of the chosen modalities. Notably, we naturally observed that higher $\theta$ leads to a decrease in annual irrigation and in irrigation maxima, and to a delay in irrigation. For the S experiments, the modality Dr of IA leads to higher amounts of irrigation.

## Characterization of best performing configurations

Table 2 summarizes modalities of experiments that reach high KGE values. First, we observe that the Semi-plot experiments clearly outperform Conceptual experiments in KGE performance in reproducing the reference experiment daily irrigation. Moreover, varCrop experiments are overrepresented, particularly among experiments with an excellent KGE (> 0.9). This observation is not surprising, because the selection of the modalities for automatic irrigation rules ($\theta$ and IA for S experiments) was based on their KGE_10 performance. However, some experiments with homogeneous automatic irrigation rules among crops also reach good KGE values. KGE does not reveal a strong difference between IP modalities among the performant experiments. However, among the performant experiments, those with modality 1 of CET are clearly more frequent than experiments with modality 0 of CET. This shows a clear and constant benefit of taking into account crop coefficients for IWR assessment, albeit the resulting variability is not as high as with other tested factors.

## A look at the monthly time step

In Fig. 5, partial variances of factors in monthly indicators are presented. The reactions of Semi-plot and Conceptual experiments with respect to the variation factors are relatively similar. The area under the curve for monthly variance of automatic irrigation parameters ($\theta$ and $\theta$+IA) shows that these parameters are the most impactful parameters on each monthly indicator. Their impact reaches a peak in July for each indicator. Consequently, the June–July period appears to be the period with the highest variance for *I_m* between experiments. However, some factors have a strong impact in other months. CET and IP have a strong impact in spring months until the beginning of summer and in autumn. Under the influence of these last two factors cumulated with automatic irrigation factors, the periods of the highest variance for *N_m* are spring and autumn, and autumn for ET_S_m. IP has a strong impact on *N_m* in winter and autumn months. Indeed, the irrigation periods defined for each crop often exclude those months, leading to large differences between experiments of modality 1 or 0 of IP factor. However, the IP factor is not as impactful on *I_m* during these periods. That can be explained by small amounts of irrigation applied in these periods, because cultivated surfaces are low, evapotranspiration is low, and rainfall is high. The variance of *N_m* is low at the end of autumn and winter because conditions for triggering irrigation are not reached and the modality 0 of IP is not sufficient to trigger irrigation in those extreme periods. The impact of CET is particularly
strong on ET_S in autumn. Indeed, in this season, the crop
coefficient is supposed to be low, but with the modality 0 of
CET, the crop coefficient is constantly equal to 1, which is
a high value.

Impact of spatial resolution on experiment outputs

For each experiment, each indicator was calculated at the CZ scale. Then the mean CZ value was computed for each experiment (Table 3). Experiments at the CZ scale reach similar values as the reference for indicators, similarly to what was observed at the lumped scale. Indeed, reference values are still contained within intervals explored by our CropWat experiments. However, we observe that KGE values at the CZ scale are overall lower than KGE for lumped outputs (Table 1), which can be explained by a higher irrigation sporadicity at the CZ scale. The best mean KGE at the CZ scale reaches 0.68, while a value of 0.91 reached at the lumped scale. There was a greater deterioration in the performance of KGE for Semi-plot experiments than for Conceptual experiments.

To complete this analysis, we drew maps show-
ing values of indicators for different CZs for the refer-
ence experiment and the experiment reaching the best
KGE_10 (Fig. 6), called the “best experiment” below.
This experiment corresponds to a Semi-plot approach,
with irrigation rules defined specifically for each crop,
no irrigation period delimitation, and root growth as
well as seasonal variation of *p* and *Kc* taken into account
(modalities varCrop, 0, 1, 1, 1 for θ+IA, IP, *z_var*, *p_var*,
CET, respectively). The spatial variability of indicators
between CZs is linked to climatic (temperature and pre-
cipitations), pedologic and agronomic spatial variabili-
ties. The spatial variability of the reference experiment
is almost perfectly reproduced by the best experiment for
I and Max_regime. Some CZs show different values for
Max_10. More important differences are observed in the
other indicators. For some indicators, we observe a bias
in comparison with the reference, but this bias is spatially
homogeneous. For example, Max_regime_date seems to
occur globally a few days later for the best experiment
than for the reference. The best experiment has more dif-
ficulties to fit reference values for *N*, *Q_dispersion*, and
Var_2. For these three indicators, the best experiment
seems to produce more spatially homogeneous results
than the reference experiment, showing that the simpli-
fied approaches tested here might face more difficulties
in reproducing spatial heterogeneity in daily irrigation
variability than the spatial heterogeneity in the other fac-
tors. Despite these differences, all spatialized indicator
values of the best experiment remain globally consistent
with those of the reference experiment.

Discussion

Impact of tested variation factors on irrigation
simulations

Our results showed that irrigation modeling choices have
an impact on irrigation modeling outputs, not only on
annual irrigation volumes, but also on the seasonal dis-
tribution of irrigation and high variations in irrigation in
short time periods. For example, the date of the maxi-
mal irrigation period (Max_regime_date) varies greatly
between 13 May and 10 August among our experiments.

Impact of modeling approaches

Tables 1 and 3 show that the ranges explored with both the
Conceptual and Semi-plot approaches are quite similar and
consistent with the reference experiment for most of indi-
cators at lumped and CZ scales. However, the Semi-plot
experiments are able to approach the results of the reference
for the Var2 indicator, which represents daily variability,
at the lumped scale, unlike Conceptual experiments. This
difference leads to lower values of KGE for the Conceptual
approach, because the daily variability of irrigation is differ-
ent from the daily variability simulated by MAELIA. Indeed,
the Conceptual approach can trigger irrigation in all simula-
tion units, which can lead to high irrigation peaks and dips,
which are not consistent with irrigation simulated by the
reference and Semi-plot approaches at the lumped scale. As
a consequence, taking into account water turns, which repro-
duce equipment availability constraints, might be decisive.
Nevertheless, several simulations of both the Conceptual and
Semi-plot approaches manage to approximate the reference
experiment, reaching high KGE_10 and KGE values at the
lumped scale: 37.5% of Conceptual experiments and 39% of
Semi-plot experiments reach a KGE higher than 0.7.

The ability to obtain good performances for both the
Semi-plot and Conceptual experiments in reproducing ref-
erence irrigation leads to the conclusion that calibration,
particularly of irrigation rules, can be more impactful than
the choice of the modeling approach among the approaches
tested. However, although adequately reproducing irrigation
at the lumped scale is possible, reproducing the daily irri-
gation simulated by the reference at the CZ scale is clearly
more difficult (see “*Impact of spatial resolution on experi-
ment outputs*”). This is easily explained by the sporadic
behavior of irrigation at the local scale.

However, we applied an approach using a fixed Kc curve
for CropWat experiments compared to an approach estimat-
ing crop growth based on a degree.day approach (Appendix
“Soil–crop water balance models”) for AqYield. Although
this difference did not seem to be very impactful for the
study period, the change of crop growth dynamic in a con-
text of climate warming might be very impactful in future. In
this context, it will become necessary to adjust Kc curves, or
to use models such as AqYield, to estimate irrigation needs.

Impact of variation factors

The impact of automatic irrigation factors (θ, IA) is very
strong among the factors tested. The impact of CET is also
important, confirming the value of taking into account
uncertainties in evapotranspiration estimation. Although the
impact of IP is not strong according to the KGE_10 indica-
tor, it has an important impact on the number of days of
irrigation. However, in a climate change context, restraining
or not restraining irrigation to specific periods might be
more impactful. Concerning the depletion factor and root
depth curves, the modalities evaluated in this work were to
take into account the variation in the parameters according
to crop growth on the one hand (modality 1), or to fix a
constant value on the other hand (modality 0). For example,
concerning the rooting depth, the maximal value could be
directly reached when the crop was sown (modality 0), or
progressively increased with crop growth for other experi-
ments (modality 1; Fig. 7). The maximal rooting depths
resulting from the crossing of crops and soils were main-
tained for all experiments. Consequently, our results show
that the dynamics of the evolution of these is not a key pro-
cess for irrigation assessment at catchment scale, and their
calibration should not be a priority. However, maximal root
depth might be a key factor of interest to evaluate in future
studies.

Beyond the statistical performance of the experiments, we
can question their agronomic relevance and robustness. For
example, extreme values of the modalities of factors inves-
tigated here can produce high KGE and KGE_10 values for
some experiments while they might also lead to unrealistic
simulations when combined with other modalities of other
factors. To reinforce the probability of modeling choices to
represent realistic irrigation over space and time, we advise
selecting realistic modalities of each variation factor repre-
sented. Following the same logic, our study also reveals the
importance of interactions between some variation factors.
Hence, with the modality 0 of IP but a high value of θ, irri-
gation can be triggered late enough in the year to reproduce
the reference scenario adequately.

There are numerous impacts related to the variation fac-
tors tested. Our study shows that evaluating and comparing
irrigation modeling based on a single indicator, for example,
annual irrigation, is not enough. Some modalities can lead
to a decrease in annual irrigation (I), but without changing
irrigation peaks in IWR (Max, Max_10, Max_regime_date)
and thereby without an impact on extreme values of hydro-
logic droughts.

Many studies approximate irrigation inputs and with-
drawals by calculating optimized values of irrigation, which
would correspond, for our CropWat simulations, to the use
of a θ equal to RAW (Funes et al. 2021; Neilsen et al. 2018;
Rinaudo et al. 2013). However, in our study site, more exper-
iments with a θ fixed at 0.5×RAW or 0.75×RAW were able
to reproduce the reference simulations very well. Our results
might be linked to the fact that farmers in our case study
tend to implement an over-irrigation strategy (Allain et al.
2018), as in other catchments (Battude 2017; Tan 2019). To
take into account this uncertainty, we would advise mod-
elers using automatic irrigation algorithms to use several
irrigation thresholds of θ to represent uncertainty linked to
farmers’ practices. Finally, defining irrigation parameters
adapted to each crop (varCrop) can lead to significantly
more accurate irrigation modeling. However, it leads to a
complexification of calibration while some experiments with
homogeneously calibrated irrigation for the different crops
were still able to reach high scores.

Potential consequences for hydrological modeling
and water management

Simulating dynamics of irrigation is of particular interest when considering the impact of irrigation during low-flow periods. Irrigation has two main impacts on hydrology: on the one hand withdrawals impacts, i.e., taking out water from the system, and on the other hand, irrigation rain impacts, bringing water to the system. Periods of irrigation, periods of maximal irrigation, and daily irrigation variations can change significantly between simulations obtained from different model configurations. A model bringing high amounts of irrigation in short periods could have an impact on hydrological modeling that is different from a model bringing low but regular amounts of irrigation.

Finally, similarly to Multsch et al. (2015), our results
showed that the evapotranspiration estimation method might
be an impactful variation factor between irrigation simula-
tions, even if this factor might be less important than the
irrigation rules. If an exhaustive coupling between crops and
hydrological modeling is intended, evapotranspiration esti-
mation might also be important as a direct input for hydro-
logical models.

Our simplified approaches (Conceptual and Semi-plot)
seem to be able to reproduce adequately the spatial vari-
ability of most indicators and should be compatible with
semi-distributed and distributed hydrological modeling.
However, we observed difficulties in reproducing irrigation
at a daily time step and at the local scale, showing the dif-
ficulty of mimicking farmers’ behavior regarding irrigation
at these scales.

<figure><img src="image_11.png" /></figure>
## Limitations to this work and other variation factors to explore

We can identify several limitations of our study. First, the benchmark irrigation data we used correspond to modeling outputs. As a consequence, these data are distinct from reality, even if they represent the best reference data existing in our study area and even if MAELIA showed a very good capacity to simulate irrigation withdrawals. Therefore, we could assume that some of our experiments, even if not identified as the best ones, could be more realistic than MAELIA simulations. This supports the need to keep several modeling hypotheses and even to keep modalities that were less performant to reproduce reference simulations if they are considered as realistic and robust.

Furthermore, this work was carried out in only one study zone. This choice is notably justified by the availability of MAELIA outputs and the complexity of the modeling protocol used, which cannot be easily generalized to other areas. Consequently, conclusions drawn in this area should be carefully used for other areas, particularly in very different agroclimatic zones and cropping systems.

The modeling approaches tested here aimed at studying different levels of modeling simplifications. However, more simplified models exist. For example, the impact of spatio-temporal aggregation of input data and the reduction in the number of simulated crops could also be explored in further studies. Moreover, the reference model, MAELIA, itself relies on several simplifications, and not all processes are described in detail. For example, the run-on of water and its impact on the redistribution of water between crop simulation units is not taken into account.

Climate data are deterministic in this study. However, for operational purposes (prediction, projection, generic characterization of irrigation distribution), climate inputs might result from hypotheses, simulations or estimations, which might also bring additional uncertainty. For example, Jie et al. (2022) considered precipitation and evapotranspiration statistical distributions as sources of variability to evaluate generic irrigation variability. Moreover, even deterministic precipitation amounts and spatio-temporal distribution carry unavoidable uncertainty. Comparing the relative impacts of climate modeling uncertainty and irrigation modeling uncertainty might be very informative. This is however not in the scope of this study.

Last but not least, coupling these simulations with hydrological modeling is necessary to confirm or contest the significance of the differences between irrigation simulations for water management and resource issues. Indeed, the differences between outputs might seem significant, but they may have a moderate impact on hydrology.

## Conclusion

This work described the methodology and the results of a sensitivity analysis of irrigation modeling at the local to river-basin scale. Two simplified modeling approaches (Conceptual: lumping irrigation simulation for homogeneous crop, soil, and climate conditions; Semi-plot: dividing simulation units into water turns groups) were compared with a more complex, agent-based, benchmark (MAELIA). For the two simplified approaches, the impacts of several modeling hypotheses regarding irrigation variation factors were analyzed. A sensitivity analysis based on variance decomposition was performed. The relative impacts of variation factors were measured based on several indicators of irrigation dynamics, with the objective of exploring the irrigation modeling effect beyond the simplistic annual sum of irrigation. This work highlighted that calibration of variation factors is more crucial than the choice of a given modeling approach. It showed the strong impact of irrigation-triggering rules and quantification of nominal irrigation amount parameters on regional irrigation assessment. It also confirmed that the definition of evapotranspiration and irrigation periods can have an important impact on irrigation modeling, a key issue for simulation under future climatic conditions. Several configurations of simpler approaches (Conceptual and Semi-plot) managed to reproduce adequately the simulations of the more complex approach (MAELIA). Experiments managing to reproduce adequately MAELIA were actually quite heterogeneous, showing a multiplicity of possible performing modeling configurations and the ability of the modalities tested to offset each other. Finally, this work enabled us to identify the following recommendations that might be followed for irrigation modeling: using multi-parameter simulations of irrigation; including different rules for triggering and quantifying irrigation; and evaluating irrigation with diverse indicators capturing its levels, frequency, and dynamics.

## Appendices

### Soil–crop water balance models

Soil water content mainly results from the balance between rain and irrigation inputs and evapotranspiration outputs. Evapotranspiration on cropland can be estimated with soil–crop water balance models. However, many of these models are complex, which makes them too computationally and data demanding for regional applications. On the other hand, semi-empiric crop coefficient methods
described by the FAO (Allen et al. 1998) are still largely used in research work and their performance and robustness, if well calibrated, have been demonstrated. The single crop coefficient (Kc) approach was implemented in the CropWat model and is used in this study.

## The CropWat model

In a single crop coefficient (Kc) approach, crop maximal evapotranspiration (CET, mm day⁻¹) is calculated at each time step (t) using Eq. (2) with the crop coefficient Kc (–) and potential evapotranspiration (PET, mm day⁻¹) estimated through the Penman–Monteith equation:

$$CET(t) = Kc(t) \times PET(t). \qquad (2)$$

In CropWat, three values of Kc are defined for each crop corresponding to the initial, mid-, and end stage of the crop cycle, linked by linear interpolation and associated with length of growth stages (Fig. 7). Water balance is calculated at a daily time step (d). Soil water availability for crops consists of a single bucket. Total available water (TAW, mm), describing the depth of the bucket, is calculated using Eq. (3), with ωfc the water content at field capacity (m³ m⁻³), ωwp the water content at wilting point (m³ m⁻³), and Zr the rooting depth (mm):

$$TAW(d) = (\omega_{fc} - \omega_{wp}) \times Z_r(d). \qquad (3)$$

The rooting depth is estimated at each time step by a linear interpolation between initial root depth and maximal root depth, the latter being reached at mid-stage. Readily available water (RAW) is calculated as follows in Eq. (4), with p (–) the depletion factor:

$$RAW(d) = p(d) \times TAW(d). \qquad (4)$$

In CropWat, p is represented by a curve with pini, pmid, pend defined for each crop associated with length of growth stages (Fig. 7).

Water level is estimated at a daily time step by the root zone depletion Dr (mm), i.e., the gap between TAW and soil water content. Soil water content is updated with daily rain, irrigation and evapotranspiration amounts. If Dr exceeds the RAW value, evapotranspiration is reduced because of water stress, leading to the calculation of actual evapotranspiration (AET, mm day⁻¹). AET is estimated with Ks (–), the water stress coefficient, as shown in Eq. (5). Ks is calculated through Eq. (6):

$$AET(d) = Ks(d) \times CET(d), \qquad (5)$$

$$Ks(d) = \begin{cases} \frac{(TAW(d) - Dr(d))}{(TAW(d) - RAW(d))} & \text{if Dr(d) > RAW(d)} \\ 1 & \text{if Dr(d) \le RAW(d)} \end{cases}. \qquad (6)$$

## The AqYield model

The MAELIA platform includes its own soil–crop water balance model, AqYield. Like CropWat, it is based on a Kc approach. The main differences are:

* Transpiration and evaporation are calculated separately. Transpiration takes water from the root zone, while evaporation takes water from the shallow soil horizon. Maximal crop transpiration (MT, mm day⁻¹) is calculated as follows in Eq. (7), with evaporation (E, mm day⁻¹):

$$MT(d) = Kc(d) \times (PET(d) - E(d)). \qquad (7)$$

* Developments of Kc and roots are represented by smooth functions depending on crop parameters and sum of degree day, with thresholds corresponding to flowering and maturity stages;

* Water stress impact on transpiration is a smooth function of the soil water amount, without any break between Dr ≤ RAW and Dr > RAW, and is influenced by the clay rate.

## Parameter calibration

MAELIA crop parameters were calibrated by experts of the AqYield model to fit the Aveyron basin context. To avoid bias linked to the crop parameters, CropWat parameters were estimated based on the MAELIA outputs. Hence, Kc curves were built on the basis of AqYield detailed outputs. However, Kc curves were kept fixed inter-annually unlike AqYield simulations. Moreover, AqYield detailed outputs were obtained for only one CZ due to constraints on data storage and computation time. As a consequence, Kc curves were built to correspond to AqYield evapotranspiration in one CZ and applied to the entire study area for CropWat experiments. A detailed explanation of CZ choice and Kc calibration methodology is given hereafter (“Choice of reference AqYield data for calibration”, “CropWat crop coefficients adjustment”).

Depletion factor p was adjusted locally with CET using the following FAO formula presented in Eq. (8):

$$p = pFAO + 0.04 \times (5 - CET), \qquad (8)$$
$p_{\text{ini}}$, $p_{\text{mid}}$, and $p_{\text{end}}$ values were calibrated to correspond to daily $p$ values calculated with the FAO formula. Minimal root depths were set to 30 cm for each crop and maximal root depths were taken from FAO report no. 56.

## Choice of reference AqYield data for calibration

AqYield detailed outputs were obtained in the CZ 2031 (Fig. 8).

The selection of this CZ was made with different criteria:

- each simulated crop is present;

- number of plots for every irrigated crop is high;

- number of different soils on which cultivated crops are present is high.

<figure><img src="image_12.png" /><figcaption>**Fig. 8** Location of the reference CZ used for CropWat crop coefficients calibration</figcaption></figure>

In CropWat, crop coefficients aim at calculating maximal crop evapotranspiration, cumulating evaporation, and transpiration, while in MAELIA, crop coefficients are designed to be proportional to maximal transpiration only. That is why we compiled for each plot the sum of maximal transpiration and evaporation simulated by AqYield between 2008 and 2014, considered as maximal evapotranspiration. Then, this maximal evapotranspiration was divided by potential evapotranspiration to obtain a crop coefficient curve comparable to the CropWat crop coefficient curve.

Finally, for each crop type, a mean daily inter-annual crop coefficient curve was calculated by the weighted mean of the crop coefficient curve of plots based on their surfaces (Fig. 9).

## CropWat crop coefficients adjustment

Each crop Kc curve was calibrated manually with the objective of reproducing adequately MAELIA outputs (Fig. 10). The intercrop and initial-stages Kc calibrated value ($Kc_{\text{ini}}$) is lower than the optimal value for this crop stage, but fits the crop Kc curve during the development stage. This choice was made to avoid overestimation of evapotranspiration during the transition between the initial and development stages, which would lead to an overestimation of soil reservoir depletion. This choice results in an underestimation of maximal evapotranspiration during the end of winter and spring, leading to a possible underestimation of water stress, but it does not have a strong impact on soil water depletion. Indeed, the water level is maintained near field capacity during this period (explaining high values of evapotranspiration in AqYield). Excess water in CropWat is not integrated in the soil water reservoir and is not converted into evapotranspiration, but is simply considered as lost water.

<figure><img src="image_13.png" /><figcaption>**Fig. 9** Inter-annual variability of crop coefficients for wheat for CZ 2031. Each colored curve corresponds to a mean crop coefficient curve weighted by the surfaces of plots for a specific year. The bold</figcaption></figure>

curve corresponds to the inter-annual mean of the colored curves (color figure online)
<figure><img src="image_14.png" /><figcaption>**Fig. 10** Annual crop coefficient curve calculated with AqYield (black line) and calibrated for CropWat (red line) for wheat (color figure online)</figcaption></figure>

Indicator calculation

<table><thead><tr><th>Indicator name</th><th>Unit</th><th>Definition</th><th>Scale of calculation</th><th>Calculation</th><th>Variables</th></tr></thead><tbody><tr><td>I</td><td>m<sup>3</sup></td><td>Mean annual irrigation</td><td>Lumped, CZ</td><td>I = <sup>N<sub>d</sub></sup>/<sub>N<sub>Y</sub></sub> I<sub>d</sub></td><td>N<sub>d</sub>: number of days (-)<br>d: day index (-)<br>I<sub>d</sub>: daily irrigation (m<sup>3</sup>)<br>N<sub>Y</sub>: number of years (-)</td></tr><tr><td>I_m</td><td>m<sup>3</sup></td><td>Mean monthly irrigation</td><td>Lumped</td><td>I_m = <sup>N<sub>d,m</sub></sup>/<sub>N<sub>Y</sub></sub> I<sub>d,m</sub></td><td>N<sub>d,m</sub>: number of days in month m (-)<br>d: day index (-)<br>I<sub>d,m</sub>: daily irrigation of day d in month m (m<sup>3</sup>)<br>N<sub>Y</sub>: number of years (-)</td></tr><tr><td>N</td><td>(-)</td><td>Mean annual number of days of irrigation</td><td>Lumped, CZ</td><td>N = <sup>N<sub>d,0</sub></sup>/<sub>N<sub>Y</sub></sub></td><td>N<sub>d,0</sub>: number of days with I<sub>d</sub> > 0 (I<sub>d</sub>: daily irrigation (m<sup>3</sup>)) (-)<br>N<sub>Y</sub>: number of years (-)</td></tr><tr><td>N_m</td><td>(-)</td><td>Mean monthly number of days of irrigation</td><td>Lumped</td><td>For each month m, N_m = <sup>N<sub>d,0,m</sub></sup>/<sub>N<sub>Y</sub></sub></td><td>N<sub>d,0,m</sub>: number of days of month m with I<sub>d,m</sub> > 0 (I<sub>d,m</sub>: daily irrigation of month m (m<sup>3</sup>)) (-)<br>N<sub>Y</sub>: number of years (-)</td></tr><tr><td>Max_ regime</td><td>m<sup>3</sup></td><td>Maximum of irrigation for 10-day rolling periods on mean annual regime curve</td><td>Lumped, CZ</td><td>Max<sub>regime</sub> = max(I<sub>10d</sub>)</td><td>I<sub>10d</sub>: vector of 10-day rolling mean annual regime of irrigation (m<sup>3</sup>,...,m<sup>3</sup>)<br>max(): maximal value (m<sup>3</sup>)</td></tr><tr><td>Max_ regime_ date</td><td>DOY</td><td>Day of year of maximum of irrigation for 10-day rolling periods on mean annual regime curve</td><td>Lumped, CZ</td><td>Max<sub>regime_date</sub> = DOY_max(I<sub>10d</sub>)</td><td>I<sub>10d</sub>: vector of 10-day rolling mean annual regime of irrigation (m<sup>3</sup>,...,m<sup>3</sup>)<br>DOY_max(): date of maximal value (DOY)</td></tr></tbody></table>
<table><thead><tr><th>Indicator name</th><th>Unit</th><th>Definition</th><th>Scale of calculation</th><th>Calculation</th><th>Variables</th></tr></thead><tbody><tr><td>Max</td><td>m³</td><td>Mean annual maximum of daily irrigation</td><td>Lumped, CZ</td><td>Max = <sup>Σ<sub>Y=1</sub><sup>N<sub>Y</sub></sup> max(I<sub>d,Y</sub>)</sup>/<sub>N<sub>Y</sub></sub></td><td>I<sub>d,Y</sub>: vector of daily irrigation for year Y (m³,...,m³)<br>max(): maximal value (m³)<br>N<sub>Y</sub>: number of years (–)</td></tr><tr><td>Max_10</td><td>m³</td><td>Mean annual maximum of 10-day mean rolling irrigation</td><td>Lumped, CZ</td><td>Max_10 = <sup>Σ<sub>Y=1</sub><sup>N<sub>Y</sub></sup> max(I<sub>10d,Y</sub>)</sup>/<sub>N<sub>Y</sub></sub></td><td>I<sub>10d,Y</sub>: vector of 10-day rolling mean of daily irrigation for year Y (m³,...,m³)<br>max(): maximal value (m³)<br>N<sub>Y</sub>: number of years (–)</td></tr><tr><td>Q_dispersion</td><td>(–)</td><td>(Q75–Q25)/Q50 of daily irrigation</td><td>Lumped, CZ</td><td>Q_dispersion = <sup>Q<sub>75</sub>(I<sub>d,d>0</sub>) – Q<sub>25</sub>(I<sub>d,d>0</sub>)</sup>/<sub>Q<sub>50</sub>(I<sub>d,d>0</sub>)</sub></td><td>I<sub>d,d>0</sub>: vector of daily irrigation I<sub>d</sub>, with I<sub>d</sub> > 0 (m³,...,m³)<br>Q<sub>d</sub>: quantile X% (m³)<br>N<sub>d</sub>: number of days (–)<br>d: day index (–)</td></tr><tr><td>Var2</td><td>(–)</td><td>Mean of absolute second derivative of daily irrigation (m³), divided by mean annual irrigation</td><td>Lumped, CZ</td><td>Var2 = <sup>Σ<sub>d=2</sub><sup>N<sub>d</sub>−1</sup> (I<sub>d+1</sub> – I<sub>d</sub>) – (I<sub>d</sub> – I<sub>d−1</sub>)</sup>/<sub>I</sub></td><td>I<sub>d</sub>: daily irrigation (m³)<br>I: see indicator I (m³)</td></tr><tr><td>KGE</td><td>(–)</td><td>Kling–Gupta efficiency comparing irrigation obtained from experiments with reference irrigation</td><td>Lumped</td><td>KGE = 1 – √<sup>(1 – r(I<sub>d</sub>))² + (1 – β(I<sub>d</sub>))² + (1 – α(I<sub>d</sub>))²</sup></td><td>r: the Pearson product–moment correlation coefficient between experiment values and reference values (–)<br>β: the ratio between the mean of the experiment values and the mean of the reference values (–)<br>α: the ratio between the standard deviation of the experiment values and the standard deviation of the reference values (–)</td></tr><tr><td>KGE_10</td><td>(–)</td><td>Kling–Gupta efficiency comparing 10-day rolling irrigation obtained from experiments with reference 10-day rolling irrigation</td><td>Lumped</td><td>KGE_10 = 1 – √<sup>(1 – r(I<sub>10d</sub>))² + (1 – β(I<sub>10d</sub>))² + (1 – α(I<sub>10d</sub>))²</sup></td><td>r: the Pearson product–moment correlation coefficient between experiment values and the reference values (–)<br>β: the ratio between the mean of the experiment values and the mean of the reference values (–)<br>α: the ratio between the standard deviation of the experiment values and the standard deviation of the reference values (–)<br>I<sub>10d</sub>: vector of 10-day rolling mean of daily irrigation for experiment and reference (m³,...,m³)</td></tr></tbody></table>
<table><thead><tr><th>Indicator name</th><th>Unit</th><th>Definition</th><th>Scale of calculation</th><th>Calculation</th><th>Variables</th></tr></thead><tbody><tr><td>ET_S</td><td>(-)</td><td>Sum of AET divided by sum of CET of cultivated surfaces during crop growth cycles</td><td>Lumped</td><td>ET_S = <sup>Σ<sub>d=1</sub><sup>N<sub>d</sub></sup> AET<sub>d</sub></sup>&frasl;<sub>Σ<sub>d=1</sub><sup>N<sub>d</sub></sup> CET<sub>d</sub></sub></td><td>N<sub>d</sub>: number of days (-)<br>d: day index(-)<br>AET<sub>d</sub>: daily actual evapotranspiration of all cultivated simulation units (between sowing date and harvesting date) (m<sup>3</sup>)<br>CET<sub>d</sub>: daily crop maximal evapotranspiration of all cultivated simulation units (between sowing date and harvesting date) (m<sup>3</sup>)</td></tr><tr><td>ET_S_m</td><td>(-)</td><td>ET_S calculated monthly</td><td>Lumped</td><td>ET_S = <sup>Σ<sub>d=1</sub><sup>N<sub>d,m</sub></sup> AET<sub>d,m</sub></sup>&frasl;<sub>Σ<sub>d=1</sub><sup>N<sub>d,m</sub></sup> CET<sub>d,m</sub></sub></td><td>N<sub>d,m</sub>: number of days in month m (-)<br>d: day index(-)<br>AET<sub>d,m</sub>: daily actual evapotranspiration of all cultivated simulation units (between sowing date and harvesting date) in month m (m<sup>3</sup>)<br>CET<sub>d,m</sub>: daily crop maximal evapotranspiration of all cultivated simulation units (between sowing date and harvesting date) in month m (m<sup>3</sup>)</td></tr></tbody></table>

## Varying factors in C and S experiments

<table><thead><tr><th>Factor name</th><th>Definition</th><th>Values</th></tr></thead><tbody><tr><td>θ</td><td>Deficit threshold for irrigation</td><td>0.25×RAW<br>0.5×RAW<br>0.75×RAW<br>RAW<br>RAW+0.25×(TAW - RAW)</td></tr><tr><td>IA<sup>a</sup></td><td>Irrigation amounts</td><td>Dr - θ: equivalent to C approach<br>Dr: irrigation amount is equal to Dr and soil is totally refilled<br>FixMeanC: irrigation amount is equal to a value defined for the crop based on MAELIA parameters</td></tr><tr><td>p_var</td><td>Depletion factor variation</td><td>1: depletion factor varies along time according to p<sub>init</sub>, p<sub>mid</sub>, p<sub>end</sub> (See Appendix “The CropWat model”)<br>0: depletion factor is constantly equal to p<sub>mid</sub></td></tr><tr><td>z_var</td><td>Root depth variation</td><td>1: root depth varies along time from minimal to maximal root depth<br>0: root depth is constantly equal to maximal root depth</td></tr><tr><td>IP</td><td>Irrigation period delimitation</td><td>1: irrigation period delimitation is injected for each crop and crops cannot be irrigated before and after this irrigation period delimitation<br>0: crops can be irrigated at any time during their crop cycle if θ is reached</td></tr><tr><td>CET</td><td>Crop modulation of potential evapotranspiration</td><td>1: CET is calculated as a modulation of Penman-Monteith PET<br>2: CET is considered equal to Penman-Monteith PET</td></tr></tbody></table>

<footnote><sup>a</sup>IA variations were simulated only for S experiments</footnote>
## Variance decomposition

Conceptual and Semi-plot approaches consist of balanced simulation plans, allowing for simple variance decomposition. Variance decomposition is used to estimate the sensitivity of IWR indicators to variation factors (θ, IP, etc.). For each indicator I and each factor F, we calculate partial variance $V_{I,F}$ with Eq. (9), with $N_{F,exp}$ the number of experiments for each modality of F, $N_F$ the number of modalities of F, $\overline{X}_i$ the mean value of F for modality I, $\overline{X}$ the mean value of F, and $N_{exp}$ the total number of experiments:

$$V_{I,F} = N_{F,exp} \frac{\sum_{i=1}^{N_F} (\overline{X}_i - \overline{X})^2}{N_{exp} - 1}. \qquad (9)$$

Moreover, sensitivity to first-order interactions is also calculated with Eq. (10). For a factor $F_1$ and a factor $F_2$, the sensitivity to their interaction is $V_{I,F_1,F_2}$, with $N_{F_1 \cap F_2,exp}$ the number of experiments for each crossing modality of $F_1$ and $F_2$:

$$V_{I,F_1,F_2} = N_{F_1 \cap F_2,exp} \frac{\sum_{j=1}^{N_{F_2}} \sum_{i=1}^{N_{F_1}} (\overline{X}_{i,j} - \overline{X}_{\cdot,j} - \overline{X}_{i,\cdot} + \overline{X}_{\cdot,\cdot})^2}{N_{exp} - 1}. \qquad (10)$$

Then partial variance can be divided by the total variance to get the contribution of a factor to total variability.

## Annual irrigation hydrographs

To complete our work based on inter-annual indicators, annual hydrographs are produced for dry-hot (2009) and wet-cold (2013) years globally over the basin. First, we observe the impact of inter-annual variability on the duration of crop cycles in MAELIA, which is not taken into account in other approaches: in a dry-hot year, MAELIA irrigation ends before the other experiments (Fig. 11), contrary to a cold-wet year (Fig. 12). We observe a higher difficulty to reproduce MAELIA irrigation in spring when irrigation is low. The same difficulty might be found at the CZ scale. We observe that the Conceptual experiments produce abrupt dips contrary to the Semi-plot experiments, explaining the lower KGE values and the higher Var2 values for the Conceptual experiments at the lumped scale. Irrigation during a wet-cold year seems to be more sporadic for all experiments, which can probably be explained by summer rain events limiting the irrigation needs during some short periods.

<figure><img src="image_15.png" /><figcaption>**Fig. 11** Comparison between CropWat experiments daily irrigation (salmon) and MAELIA reference irrigation (turquoise) for the year 2009 (dry) globally over the basin. C40 (a) corresponds to the Conceptual experiment yielding the best KGE value without varCrop modality; C0_0_1_1 (b) corresponds to the Conceptual experiment</figcaption></figure>

yielding the best KGE value with varCrop modality; S289 (c) corresponds to the Semi-plot experiment yielding the best KGE value without varCrop modality; S1_1_1_0 (d) corresponds to the Semi-plot experiment yielding the best KGE value with varCrop modality
<figure><img src="image_16.png" /><figcaption>**Fig. 12** Comparison between CropWat experiments daily irrigation (salmon) and MAELIA reference irrigation (turquoise) for the year 2013 (wet) globally over the basin. C40 (a) corresponds to the Conceptual experiment yielding the best KGE value without varCrop modality; C0_0_1_1 (b) corresponds to the Conceptual experiment</figcaption></figure>

yielding the best KGE value with varCrop modality; S289 (c) corresponds to the Semi-plot experiment yielding the best KGE value without varCrop modality; S1_1_1_0 (d) corresponds to the Semi-plot experiment yielding the best KGE value with varCrop modality

## Impact of IA modalities on Var2 and Q_dispersion

Figure 13 shows the high influence of IA variation factor on Var2 and Q_dispersion indicators.

## Impact of factor modalities on indicators

The impact of the different modalities of each factor on indicators is summarized in Fig. 14, compared with the MAELIA benchmark indicator values. The indicators for annual irrigation (I) and number of irrigated days (N) increase if θ

decreases. Indeed, if θ is higher, irrigation is triggered for higher values of Dr. As a consequence, irrigation is triggered less frequently in simulation units and later in the year. I and N are lower for modality 1 of IP than for modality 0, since for modality 1, irrigation can only be triggered during specific periods for each crop. I and N are lower for modality 1 of CET. Indeed, the Kc curves for modality 1 of CET result in globally lower evapotranspiration than evapotranspiration estimated from PET directly in modality 0. A higher evapotranspiration leads to an increase in instantaneous IWR, resulting in increased annual irrigation and number of

<figure><img src="image_17.png" /><figcaption>**Fig. 13** Impacts of automatic irrigation parameters on daily irrigation variability for Semi-plot experiments. IA modalities have a strong impact on Var2 (a) and Q_dispersion (b) values</figcaption></figure>

<figure><img src="image_18.png" /></figure>
<figure><img src="image_19.png" /><figcaption>**Fig. 14** Impact of the modalities of each factor on indicators compared with their reference values. The mean indicator value was calculated for each modality of each factor, then compared with the reference value (MAELIA). Differences between mean values for each modality and the reference value were then divided by the maximum absolute difference value for each indicator to get a relative variation between -1 and +1. Red: the mean modality value exceeds the reference value; blue: the mean modality value is below the reference</figcaption></figure>

value; yellow: the mean modality value is equal to the reference value. The color intensities are related to the variance of each factor separately, not to the total variance of all experiments. Consequently, some colors might be intense, but the total impact of this modality might remain relatively low compared to total variance. For the evaluation the partial effect of a variation factor, please refer to Fig. 4 (color figure online)

days of irrigation. For IA modalities, *I* is minimal and below the reference for Dr-θ, followed by varCrop near the reference value, by FixMeanC, and finally by Dr, both exceeding the reference value. It is clear that nominal amounts are higher for the Dr modality than for the Dr-θ modality, and this result shows that irrigation amounts brought with the Dr modality are globally higher than FixMeanC irrigation amounts. *N* is maximal for FixMeanC, followed by Dr-θ, by VarCrop, and by Dr. This observation can be linked to the explanation given for *I*: If nominal irrigation amounts are higher in simulation units, the frequency of irrigation is lower. Indeed, after an irrigation event in a simulation unit, if the irrigation amount was low, Dr after irrigation remains relatively high, and θ will be reached again after a shorter time than for a higher irrigation amount.

The maximum of irrigation for 10-day rolling periods on the mean annual regime curve (Max_regime) and its date of occurrence (Max_regime_date) are analyzed here. Max_regime_date is reached later for high values of θ. Indeed, with higher θ, irrigation is triggered later in simulation units, which results in a lag for the period of maximal irrigation. Max_regime decreases if θ increases. For lower values of θ, θ might be reached by more simulation units simultaneously, leading to higher values of irrigation during the maximum

irrigation period. Max_regime is higher than the reference value for the Dr modality of IA, and lower than the reference value for the other modalities, particularly the Dr-θ modality. As for annual irrigation, it is clear that nominal amounts are higher for the Dr modality than for the Dr-θ modality, leading to higher lumped irrigation amounts during the period of maximal irrigation. Again, as for annual irrigation, the difference in Max_regime between the Dr and FixMeanC modalities can be explained by higher irrigation amounts with the Dr modality than with the FixMeanC modality during the period of maximal irrigation. Finally, we observe that Max_regime is slightly higher for modality 1 of CET. During the annual maximal irrigation period, some major crops have a Kc value higher than 1.0, leading to higher evapotranspiration than modality 0 of CET, which might lead to an increase in IWR during this specific period. Max_regime_date is seen to occur later for modality 1 of CET. During the period preceding the maximal irrigation period, evapotranspiration is globally lower with modality 1 of CET than evapotranspiration estimated with modality 0 of CET, leading to a temporal dephasing of irrigation-triggering conditions.

The mean annual maximum of daily (Max) irrigation and 10-day (Max_10) irrigation indicators are both globally
higher than the reference. Max_10 and Max slightly decrease with higher values of θ, in a similar manner to the Max_ regime. For S experiments, Max_10 and Max are higher for the Dr modality than for the other three modalities.

Regarding daily variability, indicators Var2 and Qdispersion have opposite behavior, showing that these indicators measure different aspects of temporal variability. For C experiments, Var2 remains globally superior to the reference value, showing that C experiments have a narrow higher daily variability than reference and S experiments (Appendix “Annual irrigation hydrographs”). However, Qdispersion is globally lower than for the reference experiment. This result can be obtained for simulation showing unstable variations between successive days, for example, increasing and decreasing very frequently, but keeping lumped irrigation values in the same order of magnitude. On the contrary, a theoretical experiment with a regular increase in irrigation from day to day would lead to a low value of Var2, but a high value of Qdispersion.

Regarding crop water stress, the ET_S indicator decreases (meaning higher crop water stress) with higher modalities of θ. Obviously, with higher values of θ, Dr can be higher and as a consequence crop water stress too. However, we notice that only the modality RAW + 0.25 × (TAW – RAW), consisting of deficit irrigation, leads to lower ET_S (higher crop water stress) than the reference experiment. This can be explained by the difference in AET formulation between AqYield and CropWat. In AqYield, water stress is a smooth function of soil water deficit, with crop water starting for a zero deficit and progressively accelerating with an increasing deficit. In CropWat, there is a threshold effect since crop water stress begins when the deficit reaches RAW and is quickly high. Moreover, the calibration methodology of the Kc curve might also partly explain these differences (see Appendix “CropWat crop coefficients adjustment”). ET_S is higher (lower crop water stress) with modality 1 of CET. Indeed, in this case, the Kc curves result in globally lower evapotranspiration than evapotranspiration estimated from PET directly, and a lower evapotranspiration might lead to a lower crop water stress if conditions of irrigation are not triggered. ET_S is lower (higher crop water stress) for modality 1 of IP. Indeed, for this modality, irrigation cannot be triggered outside irrigation periods, leading to higher crop water stress. For IA impacts, ET_S is minimal for Dr-θ but still higher than the reference, followed by varCrop, FixMean_C, and then Dr: the higher the irrigation amounts for the same crop evapotranspiration, the lower the crop water stress.

**Acknowledgements** We are very grateful to Vazken Andreassian, Charles Perrin, and Gaëlle Tallec for their advice and for their re-reading of this manuscript, and to Renaud Misslin for his support on the MAELIA platform. We are also thankful to Météo-France

for providing climate data. Finally, we thank the French Agriculture Ministry and Sorbonne University for supporting the PhD of Myriam Soutif-Bellenger.

**Author contributions** MS-B designed the experimental protocol and the analysis strategy for the results. The complete protocol and analysis strategy for the results were amended and validated by GT, OT and JV. MS-B coded the Conceptual and Semi-plot experiments and performed the analysis. JV coded outputs of the MAELIA platform and supervised the use of MAELIA platform. GT supervised the Conceptual and Semi-plot experiments and analysis. OT suggested supplementary analysis. The article was written by MS-B and corrected and amended by GT and OT.

**Funding** This work was supported by Agence de l'eau Rhin-Meuse (REG-2020-00735) and CLIMAE metaprogram of INRAE (MP-P02051).

**Availability of data and materials** The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

## Declarations

**Conflict of interest** The authors have no relevant financial or non-financial interests to disclose.

**Ethical approval** Not applicable.

**Consent to participate** Not applicable.

**Consent to publish** Not applicable.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

Allain S, Ndong GO, Lardy R, Leenhardt D (2018) Integrated assessment of four strategies for solving water imbalance in an agricultural landscape. Agron Sustain Dev 38:60. https://doi.org/10.1007/s13593-018-0529-z

Allen RG, Pereira LS, Raes D, Smith M (1998) FAO irrigation and drainage paper No. 56. Rome Food Agric. Organ. U. N. 56, p 156

Battude M (2017) Estimation des rendements, des besoins et consommations en eau du maïs dans le Sud—Ouest de la France: apport de la télédétection à hautes résolutions spatiale et temporelle. Université Toulouse 3 Paul Sabatier, Toulouse

Bergez JE, Leenhardt D, Colomb B, Dury J, Carpani M, Casagrande M, Charron MH, Guillaume S, Therond O, Willaume M (2012) Computer-model tools for a better agricultural water management:
tackling managers' issues at different scales—a contribution from systemic agronomists. Comput Electron Agric 86:89–99. https://doi.org/10.1016/j.compag.2012.04.005

Bouras E, Jarlan L, Khabba S, Er-Raki S, Dezetter A, Sghir F, Tramblay Y (2019) Assessing the impact of global climate changes on irrigated wheat yields and water requirements in a semi-arid environment of Morocco. Sci Rep 9:19142. https://doi.org/10.1038/s41598-019-55251-2

Collet L (2013) Capacité à satisfaire la demande en eau sous contraintes climatique et anthropique sur un bassin méditerranéen (Thèse de doctorat). Université Montpellier 2, Montpellier

Collet L, Ruelland D, Borrell-Estupina V, Dezetter A, Servat E (2013) Integrated modelling to assess long-term water supply capacity of a meso-scale Mediterranean catchment. Sci Total Environ 461-462:528-540. https://doi.org/10.1016/j.scitotenv.2013.05.036

Collet L, Ruelland D, Estupina VB, Dezetter A, Servat E (2015) Water supply sustainability and adaptation strategies under anthropogenic and climatic changes of a meso-scale Mediterranean catchment. Sci Total Environ 536:589–602. https://doi.org/10.1016/j.scitotenv.2015.07.093

Constantin J, Willaume M, Murgue C, Lacroix B, Therond O (2015) The soil-crop models STICS and AqYield predict yield and soil water content for irrigated crops equally well with limited data. Agric for Meteorol 206:55–68. https://doi.org/10.1016/j.agrfo rmet.2015.02.011

Dehghanipour AH, Schoups G, Zahabioun B, Babazadeh H (2020) Meeting agricultural and environmental water demand in endorheic irrigated river basins: a simulation-optimization approach applied to the Urmia Lake basin in Iran. Agric Water Manag 241:106353. https://doi.org/10.1016/j.agwat.2020.106353

Di Paola A, Valentini R, Santini M (2016) An overview of available
crop growth and yield models for studies and assessments in agric-
culture: overview of crop models for agriculture. J Sci Food Agric
96:709–714. https://doi.org/10.1002/jsfa.7359

Elliott J, Deryng D, Müller C, Friel er K, Konzmann M, Gerten D, Glot-
ter M, Flörke M, Wada Y, Best N, Eisner S, Fekete BM, Folberth
C, Foster I, Gosling SN, Haddeland I, Khabarov N, Ludwig F,
Masaki Y, Olin S, Rosenzweig C, Ruane AC, Satoh Y, Schmid
E, Stacke T, Tang Q, Wisser D (2014) Constraints and potentials
of future irrigation water availability on agricultural production
under climate change. Proc Natl Acad Sci 111:3239–3244. https://
doi.org/10.1073/pnas.1222474110

Funes I, Savé R, de Herralde F, Biel C, Pla E, Pascual D, Zabalza
J, Cantos G, Borràs G, Vayreda J, Aranda X (2021) Modeling
impacts of climate change on the water needs and growing cycle
of crops in three Mediterranean basins. Agric Water Manag
249:106797. https://doi.org/10.1016/j.agwat.2021.106797

Gorguner M, Kavvas ML (2020) Modeling impacts of future climate change on reservoir storages and irrigation water demands in a Mediterranean basin. Sci Total Environ 748:141246. https://doi.org/10.1016/j.scitotenv.2020.141246

Hori T, Sugimoto T, Nakayama M, Ichikawa Y, Shiiba M (2008) Esti-
mation of field irrigation water demand based on lumped kin-
ematic wave model considering soil moisture balance. Phys Chem
Earth Parts ABC 33:376–381. https://doi.org/10.1016/j.pce.2008.
02.013

INRA (2018) Base de Données Géographique des Sols de France à 1/1 000 000 version 3.2.8.0, 10/09/1998. https://doi.org/10.15454/BPN57S

IPCC (2014) Summary for policymakers. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL (eds) Climate change 2014: impacts, adaptation, and vulnerability. Part A: global and sectoral aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate

Change. Cambridge University Press, Cambridge, United King-
dom and New York, NY, USA.

Jie F, Fei L, Li S, Hao K, Liu L, Peng Y (2022) Effects on net irriga-
tion water requirement of joint distribution of precipitation and
reference evapotranspiration. Agriculture 12:801. https://doi.org/
10.3390/agriculture12060801

Kolokytha E, Malamataris D (2020) Integrated water management approach for adaptation to climate change in highly water stressed basins. Water Resour Manag 34:1173–1197. https://doi.org/10.1007/s11269-020-02492-w

Martin E, Gascoin S, Grusson Y, Murgue C, Bardeau M, Anctil F, Fer-
rant S, Lardy R, Le Moigne P, Leenhardt D, Rivalland V, Sánchez
Pérez J-M, Sauvage S, Therond O (2016) On the use of hydrologi-
cal models and satellite data to study the water budget of river
basins affected by human activities: examples from the Garonne
Basin of France. Surv Geophys 37:223–247. https://doi.org/10.
1007/s10712-016-9366-2

Mazzega P, Therond O, Debril T, March H, Sibertin-Blanc C, Lardy R, Santana D (2014) Critical multi-level governance issues of integrated modelling: an example of low-water management in the Adour-Garonne basin (France). J Hydrol 519:2515–2526. https://doi.org/10.1016/j.jhydrol.2014.09.043

McInerney D, Thyer M, Kavetski D, Githubi F, Thayalakumaran T, Liu M, Kuczera G (2018) The importance of spatiotemporal variability in irrigation inputs for hydrological modeling of irrigated catchments. Water Resour Res 54:6792–6821. https://doi.org/10.1029/2017WR022049

MEDDTL (2011) Circulaire du 18 mai 2011 relative aux mesures exceptionnelles de limitation ou de suspension des usages de l'eau en période de sécheresse

MTES, MAA (2019) Instruction du Gouvernement du 7 mai 2019 rela-
tive au projet de territoire pour la gestion de l'eau 19

Multsch S, Exbrayat J-F, Kirby M, Viney NR, Frede H-G, Breuer L (2015) Reduction of predictive uncertainty in estimating irrigation water requirement through multi-model ensembles and ensemble averaging. Geosci Model Dev 8:1233–1244. https://doi.org/10.5194/gmd-8-1233-2015

Murgue C, Lardy R, Vavasseur M, Burger-Leenhardt D, Therond O (2014) Fine spatio-temporal simulation of cropping and farming systems effects on irrigation withdrawal dynamics within a river basin. In: Ames DP, Quinn NWT. Rizzoli AE (eds) Presented at the 7th int. congress on Env. Modelling and Software (iEMSs), San Diego, CA, USA, p 10

Murgue C, Therond O, Leenhardt D (2015) Toward integrated water and agricultural land management: participatory design of agricultural landscapes. Land Use Policy 45:52–63. https://doi.org/10.1016/j.landusepol.2015.01.011

Murgue C, Therond O, Leenhardt D (2016) Hybridizing local and generic information to model cropping system spatial distribution in an agricultural landscape. Land Use Policy 54:339–354. https://doi.org/10.1016/j.landusepol.2016.02.020

Neilsen D, Bakker M, Van der Gulik T, Smith S, Cannon A, Losso
I, Warwick Sears A (2018) Landscape based agricultural water
demand modeling—a tool for water management decision making
in British Columbia, Canada. Front Environ Sci 6:74. https://doi.
org/10.3389/fenvs.2018.00074

Rinaudo J-D, Maton L, Terrason I, Chazot S, Richard-Ferroudji A, Caballero Y (2013) Combining scenario workshops with modeling to assess future irrigation water demands. Agric Water Manag 130:103–112. https://doi.org/10.1016/j.agwat.2013.08.016

Rouhi Rad M, Araya A, Zambreski ZT (2020) Downside risk of aquifer depletion. Irrig Sci 38:577–591. https://doi.org/10.1007/s00271-020-00688-x

Smith M (1992) CROPWAT: a computer program for irrigation planning and management. Version 5.7. Food and Agriculture Organization of the United Nations, Rome, p 1992
Smith P, Calanca P, Fuhrer J (2012) A simple scheme for modeling irrigation water requirements at the regional scale applied to an Alpine River Catchment. Water 4:869–886. https://doi.org/10.3390/w4040869

Tan M (2019) Mismatch between crop water requirements and irrigation in Heihe River Basin, Northwestern China. Irrig Sci 37:95–103. https://doi.org/10.1007/s00271-018-0608-6

Therond O, Sibertin-Blanc C, Lardy R, Gaudou B, Balestrat M, Hong Y, Louail T, Nguyen VB, Panzoli D, Sanchez-Pérez J-M, Sauvage S, Taillandier P, Vavasseur M, Mazzega P (2014) Integrated modelling of social-ecological systems: the MAELIA high-resolution multi-agent platform to deal with water scarcity problems. In: Ames DP, Quinn NWT, Rizzoli AE (eds) Presented at the 7th int. congress on Env. Modelling and software (iEMSs), San Diego, CA, USA, p 8

Uniyal B, Dietrich J (2019) Modifying automatic irrigation in SWAT for plant water stress scheduling. Agric Water Manag 223:105714. https://doi.org/10.1016/j.agwat.2019.105714

Vidal J-P, Martin E, Franchistéguy L, Baillon M, Soubeyroux J-M (2010) A 50-year high-resolution atmospheric reanalysis over France with the Safran system. Int J Climatol 30:1627–1644. https://doi.org/10.1002/joc.2003

Wada Y, Wisser D, Eisner S, Flörke M, Gerten D, Haddeland I, Hanasaki N, Masaki Y, Portmann FT, Stacke T, Tessler Z, Schewe

J (2013) Multimodel projections and uncertainties of irrigation water demand under climate change: irrigation demand under climate change. Geophys Res Lett 40:4626–4632. https://doi.org/10.1002/grl.50686

Wan W, Zhao J, Li H-Y, Mishra A, Hejazi M, Lu H, Demissie Y, Wang H (2018) A holistic view of water management impacts on future droughts: a global multimodel analysis. J Geophys Res Atmos 123:5947–5972. https://doi.org/10.1029/2017JD027825

Wanders N, Wada Y (2015) Human and climate impacts on the 21st century hydrological drought. J Hydrol Drought Process Model Mitig 526:208–220. https://doi.org/10.1016/j.jhydrol.2014.10.047

Wriedt G, Van der Velde M, Aloe A, Bouraoui F (2009) Estimating irrigation water requirements in Europe. J Hydrol 373:527–544. https://doi.org/10.1016/j.jhydrol.2009.05.018

Zaccaria D, Lamaddalena N, Neale CMU, Merkley GP (2013) Simulation of peak-demand hydrographs in pressurized irrigation delivery systems using a deterministic–stochastic combined model. Part II: model applications. Irrig Sci 31:193–208. https://doi.org/10.1007/s00271-011-0308-y

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.