# How to account for irrigation withdrawals in a watershed model

Elisabeth Brochet, Youen Grusson, Sabine Sauvage, Ludovic Lhuissier, Valérie Demarez

► To cite this version:

Elisabeth Brochet, Youen Grusson, Sabine Sauvage, Ludovic Lhuissier, Valérie Demarez. How to account for irrigation withdrawals in a watershed model. Hydrology and Earth System Sciences, 2024, 28 (1), pp.49-64. 10.5194/hess-28-49-2024. hal-04607383

HAL Id: hal-04607383

https://hal.science/hal-04607383

Submitted on 10 Jun 2024

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
# How to account for irrigation withdrawals in a watershed model

Elisabeth Brochet<fnref n="1" />, Youen Grusson<fnref n="1" />, Sabine Sauvage<fnref n="2" />, Ludovic Lhuissier<fnref n="3" />, and Valérie Demarez<fnref n="1" />

**Correspondence:** Sabine Sauvage (sabine.sauvage@univ-tlse3.fr) and Valérie Demarez (valerie.demarez@univ-tlse3.fr)

Received: 18 March 2023 – Discussion started: 24 April 2023

Revised: 28 July 2023 – Accepted: 10 October 2023 – Published: 3 January 2024

**Abstract.** In agricultural areas, the downstream flow can be highly influenced by human activities during low-flow periods, especially during dam releases and irrigation withdrawals. Irrigation is indeed the major use of freshwater in the world. This study aims at precisely taking these factors into account in a watershed model. The Soil and Water Assessment Tool (SWAT+) agro-hydrological model was chosen for its capacity to model crop dynamics and management. Two different crop models were compared in terms of their ability to estimate water needs and actual irrigation. The first crop model is based on air temperature as the main determining factor for growth, whereas the second relies on high-resolution data from the Sentinel-2 satellite to monitor plant growth. Both are applied at the plot scale in a watershed of 800 km<fnref n="2" /> that is characterized by irrigation withdrawals. Results show that including remote sensing data leads to more realistic modeled emergence dates for summer crops. However, both approaches have proven to be able to reproduce the evolution of daily irrigation withdrawals throughout the year. As a result, both approaches allowed us to simulate the downstream flow with a good daily accuracy, especially during low-flow periods.

## 1 Introduction

The water cycle is substantially affected by climate change. The scientific consortium World Weather Attribution has shown that the 2022 climatic and agricultural drought in Europe would have been about 5 times less likely to happen in

the 1900s, and the frequency of such drought will increase further in the coming years. This 2022 drought event has been triggered not only by heatwaves but also by a lack of soil water storage (Schumacher et al., 2022). Many regions suffer from water scarcity, inducing some conflicts between users. Water demands for drinking water and irrigation are, in addition, expected to increase due to population growth. Water management during low-flow periods will therefore be one of the major challenges in future years. Agriculture currently accounts for about 70 % of withdrawals in the world (UNESCO, 2015). In this context, it is important to better assess water needs for agriculture. Efficient tools are then needed to better know the cultivation practices, including irrigated areas, crop dynamics and irrigation practices. These parameters must be assessed in terms of the spatial extents relevant for water management, such as at the river basin or sub-basin scale.

Watershed models are already used by multiple actors on multiple timescales and to different spatial extents. Combined with climate change scenarios, they allow a better understanding of the upcoming challenges for the next decades. For instance, the Explore2070 project (Carroget et al., 2017) has shown with high confidence that the summer streamflow in France in the years 2050–2070 will be, on average, 30 % to 70 % lower than during the years 1990–2010. Combined with management scenarios, watershed models can support decision making for agricultural policies (Murgue et al., 2014; Allain et al., 2018). Furthermore, dam managers and local authorities can make use of short-term forecasts to optimize the low-flow management (E-tiage software: https://www.e-tiage.com/, last access: 9 October

<footnote n="1">¹Centre d'Etudes Spatiales de la Biosphère, Université de Toulouse, UT3/CNRS/IRD/CNES, 18 av. Edouard Belin, bpi 2801, 931401 Toulouse CEDEX, France</footnote>

<footnote n="2">²Centre de Recherche sur la Biodiversité et l'Environnement, ENSAT, CNRS/UT3/INPT, Av. de l'Agrobiopole, bpi 32607, 31326 Castanet Tolosan CEDEX, France</footnote>

<footnote n="3">³Compagnie d'Aménagement de Côteaux de Gascogne, Chemin de Lalette, CS 50449, 65004 Tarbes CEDEX, France</footnote>
# How to account for irrigation withdrawals in a watershed model

Elisabeth Brochet<fnref n="1" />, Youen Grusson<fnref n="1" />, Sabine Sauvage<fnref n="2" />, Ludovic Lhuissier<fnref n="3" />, and Valérie Demarez<fnref n="1" />

**Correspondence:** Sabine Sauvage (sabine.sauvage@univ-tlse3.fr) and Valérie Demarez (valerie.demarez@univ-tlse3.fr)

Received: 18 March 2023 – Discussion started: 24 April 2023

Revised: 28 July 2023 – Accepted: 10 October 2023 – Published: 3 January 2024

**Abstract.** In agricultural areas, the downstream flow can be highly influenced by human activities during low-flow periods, especially during dam releases and irrigation withdrawals. Irrigation is indeed the major use of freshwater in the world. This study aims at precisely taking these factors into account in a watershed model. The Soil and Water Assessment Tool (SWAT+) agro-hydrological model was chosen for its capacity to model crop dynamics and management. Two different crop models were compared in terms of their ability to estimate water needs and actual irrigation. The first crop model is based on air temperature as the main determining factor for growth, whereas the second relies on high-resolution data from the Sentinel-2 satellite to monitor plant growth. Both are applied at the plot scale in a watershed of 800 km<fnref n="2" /> that is characterized by irrigation withdrawals. Results show that including remote sensing data leads to more realistic modeled emergence dates for summer crops. However, both approaches have proven to be able to reproduce the evolution of daily irrigation withdrawals throughout the year. As a result, both approaches allowed us to simulate the downstream flow with a good daily accuracy, especially during low-flow periods.

## 1 Introduction

The water cycle is substantially affected by climate change. The scientific consortium World Weather Attribution has shown that the 2022 climatic and agricultural drought in Europe would have been about 5 times less likely to happen in

the 1900s, and the frequency of such drought will increase further in the coming years. This 2022 drought event has been triggered not only by heatwaves but also by a lack of soil water storage (Schumacher et al., 2022). Many regions suffer from water scarcity, inducing some conflicts between users. Water demands for drinking water and irrigation are, in addition, expected to increase due to population growth. Water management during low-flow periods will therefore be one of the major challenges in future years. Agriculture currently accounts for about 70 % of withdrawals in the world (UNESCO, 2015). In this context, it is important to better assess water needs for agriculture. Efficient tools are then needed to better know the cultivation practices, including irrigated areas, crop dynamics and irrigation practices. These parameters must be assessed in terms of the spatial extents relevant for water management, such as at the river basin or sub-basin scale.

Watershed models are already used by multiple actors on multiple timescales and to different spatial extents. Combined with climate change scenarios, they allow a better understanding of the upcoming challenges for the next decades. For instance, the Explore2070 project (Carroget et al., 2017) has shown with high confidence that the summer streamflow in France in the years 2050–2070 will be, on average, 30 % to 70 % lower than during the years 1990–2010. Combined with management scenarios, watershed models can support decision making for agricultural policies (Murgue et al., 2014; Allain et al., 2018). Furthermore, dam managers and local authorities can make use of short-term forecasts to optimize the low-flow management (E-tiage software: https://www.e-tiage.com/, last access: 9 October

<footnote n="1">¹Centre d'Etudes Spatiales de la Biosphère, Université de Toulouse, UT3/CNRS/IRD/CNES, 18 av. Edouard Belin, bpi 2801, 931401 Toulouse CEDEX, France</footnote>

<footnote n="2">²Centre de Recherche sur la Biodiversité et l'Environnement, ENSAT, CNRS/UT3/INPT, Av. de l'Agrobiopole, bpi 32607, 31326 Castanet Tolosan CEDEX, France</footnote>

<footnote n="3">³Compagnie d'Aménagement de Côteaux de Gascogne, Chemin de Lalette, CS 50449, 65004 Tarbes CEDEX, France</footnote>
2022; PREMHYCE platform https://webgr.inrae.fr/projets/projets-en-cours/onema-premhyce/, last access: 18 October 2022; Nicolle et al., 2020). Hence the need for continuous improvement of these models.

However, most studies at the watershed scale do not precisely include human activities, especially withdrawals for irrigation. For instance, most users of the SWAT (Soil and Water Assessment Tool) model do not explicitly include withdrawals in their studies (Fohrer et al., 2014; Boithias et al., 2014; Martin et al., 2016; Cakir et al., 2020). In such studies, dam releases and withdrawals are indirectly taken into account through the calibration performed to fit the downstream flow. As a result, the performance of low-flow simulation is quite low, as highlighted by Boithias et al. (2014).

Some modelers include cultivation practices through default crop management or through past years' statistics, but these models have shown limitations for the prediction of irrigation. Indeed, Leenhardt and Lemaire (2002) describe a case in which the sowing of summer crops was delayed due to the unusually wet spring, but the model did not account for this and predicted the maximum irrigation demand 1 month too early. In the same way, Senthilkumar et al. (2015) highlight the impact of cultivar earliness on their water needs. Models that better fit the actual crops dynamics and management are therefore needed. Remote sensing (RS) offers a real added value for spatial and temporal variability. Since 2016, thanks to the Copernicus program, optical satellite images (Sentinel-2) with high spatial (10 m) and temporal (5 d) resolutions are available anywhere in the world. These resolutions are appropriate to get information on crops dynamics and cultivation practices at the plot scale over a whole watershed.

Methods to retrieve the cultivation practices from RS data have been in development for more than 20 years (Guérif and Duke, 1998; Launay and Guerif, 2005) and have bloomed since short-revisit-time satellites were launched (Courault et al., 2010; Ferrant et al., 2014). The Sentinel constellation, offering a high spatial resolution and high revisit time, has made this task easier and more accurate (Pageot et al., 2020, for irrigated areas; Bazzi et al., 2021, for irrigation events; Rolle et al., 2022, for sowing dates).

For the purpose of water management, RS can also be combined with crop water requirements models. The FAO-56 method (Allen et al., 1998) is the most widely used land surface model in which the potential evapotranspiration (ET) is adjusted with a crop coefficient and a stress coefficient. The classical method uses standard values of crop coefficient, which implies standard growth conditions, but an increasing number of studies replace these with remote sensing data (Saadi et al., 2015; Etchanchu et al., 2017; Battude et al., 2017; Yousaf et al., 2021; Kharrou et al., 2021; Maguire et al., 2022) in order to account for the current-year conditions. In particular, Etchanchu et al. (2017) highlight the contribution of remote sensing data at the plot scale, making more robust calculations of LAI (leaf area index) and ET in

comparison to the classical method. In addition to ET calculation, these methods can be used to estimate irrigation through the soil water balance (Saadi et al., 2015; Battude et al., 2017; Kharrou et al., 2021; Maguire et al., 2022) or a comparison with a reference crop coefficient (Yousaf et al., 2021). Even if the retrieval of irrigation suffers from large uncertainties at the farm scale (Saadi et al., 2015, for instance), Olivera-Guerra et al. (2023) have shown that proper calibration can greatly reduce the uncertainties at the irrigation district scale.

Despite these recent advances in the modeling of crop dynamics and water needs, very few works investigate the effect of those cultivation practices on streamflow at high spatial and temporal scales. The MAELIA platform (Modelling of socio-Agro-Ecological system for Landscape Integrated Assessment, Murgue et al., 2014) appears as an exception since it combines the SWAT hydrological model with an agent-based model to perform dam operation and crop management. In this case, modelers have inferred a large set of decision rules either from surveys (Leenhardt et al., 2004; Maton et al., 2005, 2007; Murgue et al., 2014) or from a farm advisor (Clavel et al., 2011). The authors were able to check the relevance of the annual irrigation amount predicted by the model but not the timing of irrigation due to the lack of daily data. The present study intends to go further than the MAELIA project in terms of spatialization by substituting part of the decision rules with remote sensing estimates. Indeed, feeding the model with observations could improve its robustness while avoiding the complexity of an agent-based model.

Some studies already aimed to assess the benefits of optical remote sensing in agro-hydrological models without focusing on crop irrigation. They intended to improve the ET component of the water balance by forcing (Martin et al., 2016; Paul et al., 2021; Jin et al., 2022, all in SWAT model) or assimilating (Kumar et al., 2019, in Noah-MP model; Mohammadi Igder et al., 2022, in SWAT+ model) the RS LAI (remote sensing leaf area index) in the model and short-cutting the plant growth algorithm. Under temperate or arid climate types, ET indeed appears to be the main component of the water balance; hence, a small relative error in ET could lead to a larger relative error in the river discharge. However, even if these authors report an improvement in the yield or carbon fluxes, they report little modification to the streamflow, ranging from 2% to 7%. Explanations for these results could be that the modeled LAI, even if uncertain, is enough to compute ET or that the uncertainties in other hydrological or anthropogenic fluxes are larger than the one in ET. As Paul et al. (2021) and Jin et al. (2022) both performed calibrations with and without remote sensing data, this also raises the following question: can calibration correct the errors in ET values in a way that it hides the possible benefits of RS data? Moreover, as far as we know, no studies at the watershed scale used remote sensing data at the plot scale, even though Sentinel data would be suitable for that.
2022; PREMHYCE platform https://webgr.inrae.fr/projets/projets-en-cours/onema-premhyce/, last access: 18 October 2022; Nicolle et al., 2020). Hence the need for continuous improvement of these models.

However, most studies at the watershed scale do not precisely include human activities, especially withdrawals for irrigation. For instance, most users of the SWAT (Soil and Water Assessment Tool) model do not explicitly include withdrawals in their studies (Fohrer et al., 2014; Boithias et al., 2014; Martin et al., 2016; Cakir et al., 2020). In such studies, dam releases and withdrawals are indirectly taken into account through the calibration performed to fit the downstream flow. As a result, the performance of low-flow simulation is quite low, as highlighted by Boithias et al. (2014).

Some modelers include cultivation practices through default crop management or through past years' statistics, but these models have shown limitations for the prediction of irrigation. Indeed, Leenhardt and Lemaire (2002) describe a case in which the sowing of summer crops was delayed due to the unusually wet spring, but the model did not account for this and predicted the maximum irrigation demand 1 month too early. In the same way, Senthilkumar et al. (2015) highlight the impact of cultivar earliness on their water needs. Models that better fit the actual crops dynamics and management are therefore needed. Remote sensing (RS) offers a real added value for spatial and temporal variability. Since 2016, thanks to the Copernicus program, optical satellite images (Sentinel-2) with high spatial (10 m) and temporal (5 d) resolutions are available anywhere in the world. These resolutions are appropriate to get information on crops dynamics and cultivation practices at the plot scale over a whole watershed.

Methods to retrieve the cultivation practices from RS data have been in development for more than 20 years (Guérif and Duke, 1998; Launay and Guerif, 2005) and have bloomed since short-revisit-time satellites were launched (Courault et al., 2010; Ferrant et al., 2014). The Sentinel constellation, offering a high spatial resolution and high revisit time, has made this task easier and more accurate (Pageot et al., 2020, for irrigated areas; Bazzi et al., 2021, for irrigation events; Rolle et al., 2022, for sowing dates).

For the purpose of water management, RS can also be combined with crop water requirements models. The FAO-56 method (Allen et al., 1998) is the most widely used land surface model in which the potential evapotranspiration (ET) is adjusted with a crop coefficient and a stress coefficient. The classical method uses standard values of crop coefficient, which implies standard growth conditions, but an increasing number of studies replace these with remote sensing data (Saadi et al., 2015; Etchanchu et al., 2017; Battude et al., 2017; Yousaf et al., 2021; Kharrou et al., 2021; Maguire et al., 2022) in order to account for the current-year conditions. In particular, Etchanchu et al. (2017) highlight the contribution of remote sensing data at the plot scale, making more robust calculations of LAI (leaf area index) and ET in

comparison to the classical method. In addition to ET calculation, these methods can be used to estimate irrigation through the soil water balance (Saadi et al., 2015; Battude et al., 2017; Kharrou et al., 2021; Maguire et al., 2022) or a comparison with a reference crop coefficient (Yousaf et al., 2021). Even if the retrieval of irrigation suffers from large uncertainties at the farm scale (Saadi et al., 2015, for instance), Olivera-Guerra et al. (2023) have shown that proper calibration can greatly reduce the uncertainties at the irrigation district scale.

Despite these recent advances in the modeling of crop dynamics and water needs, very few works investigate the effect of those cultivation practices on streamflow at high spatial and temporal scales. The MAELIA platform (Modelling of socio-Agro-Ecological system for Landscape Integrated Assessment, Murgue et al., 2014) appears as an exception since it combines the SWAT hydrological model with an agent-based model to perform dam operation and crop management. In this case, modelers have inferred a large set of decision rules either from surveys (Leenhardt et al., 2004; Maton et al., 2005, 2007; Murgue et al., 2014) or from a farm advisor (Clavel et al., 2011). The authors were able to check the relevance of the annual irrigation amount predicted by the model but not the timing of irrigation due to the lack of daily data. The present study intends to go further than the MAELIA project in terms of spatialization by substituting part of the decision rules with remote sensing estimates. Indeed, feeding the model with observations could improve its robustness while avoiding the complexity of an agent-based model.

Some studies already aimed to assess the benefits of optical remote sensing in agro-hydrological models without focusing on crop irrigation. They intended to improve the ET component of the water balance by forcing (Martin et al., 2016; Paul et al., 2021; Jin et al., 2022, all in SWAT model) or assimilating (Kumar et al., 2019, in Noah-MP model; Mohammadi Igder et al., 2022, in SWAT+ model) the RS LAI (remote sensing leaf area index) in the model and short-cutting the plant growth algorithm. Under temperate or arid climate types, ET indeed appears to be the main component of the water balance; hence, a small relative error in ET could lead to a larger relative error in the river discharge. However, even if these authors report an improvement in the yield or carbon fluxes, they report little modification to the streamflow, ranging from 2% to 7%. Explanations for these results could be that the modeled LAI, even if uncertain, is enough to compute ET or that the uncertainties in other hydrological or anthropogenic fluxes are larger than the one in ET. As Paul et al. (2021) and Jin et al. (2022) both performed calibrations with and without remote sensing data, this also raises the following question: can calibration correct the errors in ET values in a way that it hides the possible benefits of RS data? Moreover, as far as we know, no studies at the watershed scale used remote sensing data at the plot scale, even though Sentinel data would be suitable for that.
The objective of this study is to take into account human activity during low-flow periods in the modeling process of an agricultural watershed. Unlike previous studies, our work is located on a watershed where irrigation has a huge impact on streamflow: in summer, the downstream flow can be 4 times lower than the upstream flow due to withdrawals. To reach this objective, crop models at the plot scale are used to predict crop water requirements and irrigation withdrawals, and the SWAT+ model is used to compute water stocks and fluxes. The native crop growth module of SWAT is compared to a new module that uses Sentinel-2 high-resolution data to compute ET. An innovative calibration method for hydrological fluxes is introduced to address several calibration issues: the short time depth of Sentinel-2 data, the fact that calibration can hide the possible benefits of RS data, and the highly influenced streamflow. Numerous data sources have been used to carry out the calibration of the irrigation trigger threshold, the estimation of irrigation sources, and the assessment of daily withdrawals. Finally, the relevance of a such model – combining the plot scale and watershed scale during the low-flow period – is assessed.

## 2 Materials and methods

### 2.1 Study area

The Gimone watershed (800 km²) has been selected for this study. It is located in southwestern France, in a temperate-climate zone. (Fig. 1). The average rainfall is 700 mm. The flow rate at the outlet ranges from 0.4 to 40 m³ s⁻¹, with an average of 3 m³ s⁻¹. Agricultural land is dominant (85%) in this watershed. Irrigated crops (mainly corn, soybean, orchards and vegetables) represent about 10% to 15% of the total area, and irrigation water is withdrawn from rivers and from individual reservoirs. Other crops, mainly winter wheat, winter rapeseed and sunflower, are not irrigated. Due to the high water demand and low water availability in summer, two dams were built upstream in the 1990s on the Gimone and Marcaoué rivers, with capacities, of respectively, 24 and 1.5 Mm³ (Fig. 1e). The biggest dam is filled with water diverted from the Neste river. Both dams are managed by the CACG (Compagnie d'Aménagement des Côteaux de Gascogne), and their purpose is to maintain a minimal flow rate in rivers during low-flow periods (DOE – Débit Objectif d'Etiage – minimum target flow of 0.4 m³ s⁻¹) and to limit the frequency of irrigation restrictions. As a result, the streamflow during irrigation periods is strongly influenced by dam releases and irrigation withdrawals and has little to do with hydrological processes.

### 2.2 Data overview

The topography of the watershed derives from the SRTM (Shuttle Radar Topography Mission, USGS, 2018) digital elevation model at ~30 m resolution (Fig. 1b). Daily weather

data – precipitation, air temperature, wind velocity, solar radiation and air humidity – were extracted from the SAFRAN product (Durand et al., 1993). It is an interpolation based on weather measurements and model reanalysis on a 8 km resolution grid. Soil properties, namely the depth, bulk density and clay and sand content, were found in the French soil map RRP (Référentiel Régional Pédologique – regional reference document for pedology, GisSol, 2014) at 1/250 000 resolution (Fig. 1c); the maximal root depths were set according to Rigou (2016); the available soil water content (AWC) was derived from Bruand et al. (2003) pedotransfer functions; the initial estimate of the saturated hydraulic conductivity comes from measurements. Land use data come from the French database of agricultural land use (RPG – Registre Parcellaire Graphique, IGN, 2021) for the years 2015 to 2021, completed by the OSO map (Inglada et al., 2017) for non-agricultural areas (Fig. 1d) (Inglada et al., 2018). Streamflow coming out of the two dams was provided by the CACG.

In addition to the classical data to set up the SWAT+ model, high-resolution remote sensing data were used. Since 2017, the two Sentinel-2 satellites provide images every 5 d, with resolutions ranging from 10 to 60 m depending on the wavelength (Drusch et al., 2012; ESA, 2012). The level-2A images from two Sentinel-2 tiles, 31TCJ and 30TYP (Fig. 1a), are used in this study. They are post-processed through the MAJA chain (Hagolle et al., 2015) that corrects atmospheric effects and produces a cloud mask. These clouds masks are used in the interpolation process to obtain gap-filled time series.

River discharge at the outlet was also provided by the CACG and has been used for calibration and validation. In order to get an overview of the origin of withdrawals, the PAR (Plan Annuel de Répartition – annual plan for water allocation) database was used. It gathers individual withdrawal permissions issued every year by the French administration and contains the authorized volumes in each type of water resource. The PKGC (Pratiques Culturales en Grandes Cultures – agricultural practices in field crops) database (Agreste, 2020), carried out every 5 years, consists of a large survey that covers all French departments, with information for about 28 000 agricultural plots. The dataset from 2017 has been used in this study, particularly the sowing dates and amounts of annual irrigation of 160 corn and soybean fields in southwestern France.

A large part of the irrigation pumps of the Gimone watershed are equipped with networked meters (Fig. 1e). These provide daily measurement of water withdrawn from the two main rivers, the same two that are regulated with dams and for which low flows are sustained during summer (Gimone and Marcaoué). Not all intake points are equipped for the years considered in the present study: 33% of the pumps were equipped in 2019, 62% were equipped in 2020, and 72% were equipped in 2021. For this reason, only 2020 and 2021 data were used and were corrected with a proportionality coefficient to convert those data into total withdrawals.
The objective of this study is to take into account human activity during low-flow periods in the modeling process of an agricultural watershed. Unlike previous studies, our work is located on a watershed where irrigation has a huge impact on streamflow: in summer, the downstream flow can be 4 times lower than the upstream flow due to withdrawals. To reach this objective, crop models at the plot scale are used to predict crop water requirements and irrigation withdrawals, and the SWAT+ model is used to compute water stocks and fluxes. The native crop growth module of SWAT is compared to a new module that uses Sentinel-2 high-resolution data to compute ET. An innovative calibration method for hydrological fluxes is introduced to address several calibration issues: the short time depth of Sentinel-2 data, the fact that calibration can hide the possible benefits of RS data, and the highly influenced streamflow. Numerous data sources have been used to carry out the calibration of the irrigation trigger threshold, the estimation of irrigation sources, and the assessment of daily withdrawals. Finally, the relevance of a such model – combining the plot scale and watershed scale during the low-flow period – is assessed.

## 2 Materials and methods

### 2.1 Study area

The Gimone watershed (800 km²) has been selected for this study. It is located in southwestern France, in a temperate-climate zone. (Fig. 1). The average rainfall is 700 mm. The flow rate at the outlet ranges from 0.4 to 40 m³ s⁻¹, with an average of 3 m³ s⁻¹. Agricultural land is dominant (85%) in this watershed. Irrigated crops (mainly corn, soybean, orchards and vegetables) represent about 10% to 15% of the total area, and irrigation water is withdrawn from rivers and from individual reservoirs. Other crops, mainly winter wheat, winter rapeseed and sunflower, are not irrigated. Due to the high water demand and low water availability in summer, two dams were built upstream in the 1990s on the Gimone and Marcaoué rivers, with capacities, of respectively, 24 and 1.5 Mm³ (Fig. 1e). The biggest dam is filled with water diverted from the Neste river. Both dams are managed by the CACG (Compagnie d'Aménagement des Côteaux de Gascogne), and their purpose is to maintain a minimal flow rate in rivers during low-flow periods (DOE – Débit Objectif d'Etiage – minimum target flow of 0.4 m³ s⁻¹) and to limit the frequency of irrigation restrictions. As a result, the streamflow during irrigation periods is strongly influenced by dam releases and irrigation withdrawals and has little to do with hydrological processes.

### 2.2 Data overview

The topography of the watershed derives from the SRTM (Shuttle Radar Topography Mission, USGS, 2018) digital elevation model at ~30 m resolution (Fig. 1b). Daily weather

data – precipitation, air temperature, wind velocity, solar radiation and air humidity – were extracted from the SAFRAN product (Durand et al., 1993). It is an interpolation based on weather measurements and model reanalysis on a 8 km resolution grid. Soil properties, namely the depth, bulk density and clay and sand content, were found in the French soil map RRP (Référentiel Régional Pédologique – regional reference document for pedology, GisSol, 2014) at 1/250 000 resolution (Fig. 1c); the maximal root depths were set according to Rigou (2016); the available soil water content (AWC) was derived from Bruand et al. (2003) pedotransfer functions; the initial estimate of the saturated hydraulic conductivity comes from measurements. Land use data come from the French database of agricultural land use (RPG – Registre Parcellaire Graphique, IGN, 2021) for the years 2015 to 2021, completed by the OSO map (Inglada et al., 2017) for non-agricultural areas (Fig. 1d) (Inglada et al., 2018). Streamflow coming out of the two dams was provided by the CACG.

In addition to the classical data to set up the SWAT+ model, high-resolution remote sensing data were used. Since 2017, the two Sentinel-2 satellites provide images every 5 d, with resolutions ranging from 10 to 60 m depending on the wavelength (Drusch et al., 2012; ESA, 2012). The level-2A images from two Sentinel-2 tiles, 31TCJ and 30TYP (Fig. 1a), are used in this study. They are post-processed through the MAJA chain (Hagolle et al., 2015) that corrects atmospheric effects and produces a cloud mask. These clouds masks are used in the interpolation process to obtain gap-filled time series.

River discharge at the outlet was also provided by the CACG and has been used for calibration and validation. In order to get an overview of the origin of withdrawals, the PAR (Plan Annuel de Répartition – annual plan for water allocation) database was used. It gathers individual withdrawal permissions issued every year by the French administration and contains the authorized volumes in each type of water resource. The PKGC (Pratiques Culturales en Grandes Cultures – agricultural practices in field crops) database (Agreste, 2020), carried out every 5 years, consists of a large survey that covers all French departments, with information for about 28 000 agricultural plots. The dataset from 2017 has been used in this study, particularly the sowing dates and amounts of annual irrigation of 160 corn and soybean fields in southwestern France.

A large part of the irrigation pumps of the Gimone watershed are equipped with networked meters (Fig. 1e). These provide daily measurement of water withdrawn from the two main rivers, the same two that are regulated with dams and for which low flows are sustained during summer (Gimone and Marcaoué). Not all intake points are equipped for the years considered in the present study: 33% of the pumps were equipped in 2019, 62% were equipped in 2020, and 72% were equipped in 2021. For this reason, only 2020 and 2021 data were used and were corrected with a proportionality coefficient to convert those data into total withdrawals.
<figure><img src="image_1.png" /><figcaption>**Figure 1.** Maps of the study area. (a) Location of the Gimone watershed in southwestern France, as well as Sentinel-2 footprints; (b) digital elevation model (1" resolution); (c) soil map (1/250 000 resolution); (d) crop pattern for year 2017 (plot scale); (e) hydrological network setup of the model.</figcaption></figure>

## 2.3 SWAT+ model

### 2.3.1 Global overview

The Soil and Water Assessment Tool (Arnold et al., 1993, 2012) is an agro-hydrological semi-distributed and process-based model. To set up this semi-distributed approach, the simulated watershed is divided into multiple sub-basins according to topography. Each sub-basin is divided into HRUs (hydrological response units), which are homogeneous areas in terms of slope, soil and land use. Calculation of the process-based water balance is then carried out at the HRU level on a daily basis (Neitsch et al., 2001). Daily rainfall is split between surface runoff and infiltration using the curve number method (Soil Conservation Service, 1972; Rallison and Miller, 1982). Soil water sustains evaporation and leads to subsurface runoff and percolation into aquifers. If the simulated water table of the shallow aquifer is high enough, ground water flow from the aquifer supports river flow. This contribution of aquifers to the river flow is called baseflow. Land use management (e.g., planting, irrigation, fertilization, harvest) occurs on the basis of user-defined decision rules. When irrigation is applied, the origin of withdrawals (river, aquifer, reservoir) must be specified so that the amount can be subtracted from the appropriate source. All fluxes reaching the rivers (surface and subsurface runoff, ground water flow) are aggregated at the sub-basin level and routed through the stream network to the next sub-basin downstream until reaching the watershed outlet.

### 2.3.2 Model setup

The 60.5.3 version of SWAT+ has been used in this study as the baseline model, to which features and options were added (see “Code availability” section).

The watershed has been divided into 149 sub-basins (Fig. 1). In this study, we forced HRUs to match individual fields in all agricultural land so that using remotely sensed indices at the plot scale would be consistent. This delineation resulted in about 26 000 HRUs. The model was fed with the true crop rotation from 2015 to 2021 according to the French database RPG.

## 2.4 Two approaches for crop growth and management

Two methods have been compared to simulate crops seasonal dynamics. The first one is the native crop growth module of SWAT, based on a simplified version of EPIC (Erosion Productivity Impact Calculator). The second one consists of observations of measured crop growth through the normalized difference vegetation index (NDVI), a remote-sensing-based index linked to vegetation development. The two approaches are hereafter called SWAT-O and SWAT-NDVI, respectively. The SWAT-O setup has also been run without any irrigation to assess the impact of withdrawals on the streamflow. Below is a short description of both methods.
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Maps of the study area. (a) Location of the Gimone watershed in southwestern France, as well as Sentinel-2 footprints; (b) digital elevation model (1" resolution); (c) soil map (1/250 000 resolution); (d) crop pattern for year 2017 (plot scale); (e) hydrological network setup of the model.</figcaption></figure>

## 2.3 SWAT+ model

### 2.3.1 Global overview

The Soil and Water Assessment Tool (Arnold et al., 1993, 2012) is an agro-hydrological semi-distributed and process-based model. To set up this semi-distributed approach, the simulated watershed is divided into multiple sub-basins according to topography. Each sub-basin is divided into HRUs (hydrological response units), which are homogeneous areas in terms of slope, soil and land use. Calculation of the process-based water balance is then carried out at the HRU level on a daily basis (Neitsch et al., 2001). Daily rainfall is split between surface runoff and infiltration using the curve number method (Soil Conservation Service, 1972; Rallison and Miller, 1982). Soil water sustains evaporation and leads to subsurface runoff and percolation into aquifers. If the simulated water table of the shallow aquifer is high enough, ground water flow from the aquifer supports river flow. This contribution of aquifers to the river flow is called baseflow. Land use management (e.g., planting, irrigation, fertilization, harvest) occurs on the basis of user-defined decision rules. When irrigation is applied, the origin of withdrawals (river, aquifer, reservoir) must be specified so that the amount can be subtracted from the appropriate source. All fluxes reaching the rivers (surface and subsurface runoff, ground water flow) are aggregated at the sub-basin level and routed through the stream network to the next sub-basin downstream until reaching the watershed outlet.

### 2.3.2 Model setup

The 60.5.3 version of SWAT+ has been used in this study as the baseline model, to which features and options were added (see “Code availability” section).

The watershed has been divided into 149 sub-basins (Fig. 1). In this study, we forced HRUs to match individual fields in all agricultural land so that using remotely sensed indices at the plot scale would be consistent. This delineation resulted in about 26 000 HRUs. The model was fed with the true crop rotation from 2015 to 2021 according to the French database RPG.

## 2.4 Two approaches for crop growth and management

Two methods have been compared to simulate crops seasonal dynamics. The first one is the native crop growth module of SWAT, based on a simplified version of EPIC (Erosion Productivity Impact Calculator). The second one consists of observations of measured crop growth through the normalized difference vegetation index (NDVI), a remote-sensing-based index linked to vegetation development. The two approaches are hereafter called SWAT-O and SWAT-NDVI, respectively. The SWAT-O setup has also been run without any irrigation to assess the impact of withdrawals on the streamflow. Below is a short description of both methods.
<table><caption>**Table 1.** Decision rules for corn emergence. The emergence occurs on the first day that all the conditions are met. YHU: total heat units of the year; SWC: soil water content; FC: field capacity.</caption><thead><tr><th>SWAT-O</th><th>SWAT-NDVI</th></tr></thead><tbody><tr><td>HU<sub>0</sub> &gt; 30 % YHU</td><td>20 April &lt; day &lt; 10 July</td></tr><tr><td>SWC &lt; 1.05 · FC</td><td>NDVI increasing<br>0.2 &lt; NDVI &lt; 0.5</td></tr></tbody></table>

## 2.4.1 SWAT-O: growth of plants and ET driven by heat units (EPIC)

Heat units (namely the average temperature of the day minus the base temperature) drive the growth of LAI (leaf area index) according to the following formulas (Barnard, 1948; Phillips, 1950):

$$ \text{fr}_{\text{LAlmx}} = \frac{\text{fr}_{\text{PHU}}}{\text{fr}_{\text{PHU}} + \exp(l_1 - l_2 \cdot \text{fr}_{\text{PHU}})}, \quad (1) $$

$$ \text{fr}_{\text{PHU}} = \frac{\sum_{i=1}^{d} \text{HU}_i}{\text{PHU}}, \quad (2) $$

where fr<sub>LAlmx</sub> is the fraction of maximum LAI, l<sub>1</sub> and l<sub>2</sub> are shape coefficients, HU<sub>i</sub> is the heat unit of day i (°C), and PHU is the potential heat unit required for plant maturity (°C). A more detailed description of this simplified EPIC module can be found in the SWAT theoretical documentation (Neitsch et al., 2001).

The beginning of growth is settled according to decision rules, where the most important parameter is the sum of heat units since 1 January (HU<sub>0</sub>). The emergence occurs on the first day that all the conditions are met. See an example for corn in Table 1.

Potential evapotranspiration (PET) of the crop depends on its LAI based on Penman-Monteith equation (Monteith, 1965):

$$ \lambda E_t = \frac{\Delta(H_{\text{net}} - G) + \gamma K_1 \cdot (0.622\lambda \cdot \frac{\rho_{\text{air}}}{P}) \cdot \frac{e_z^o - e_z}{r_a}}{\Delta + \gamma \cdot (1 + r_c/r_a)}, \quad (3) $$

$$ r_c = \frac{r_l}{0.5 \cdot \text{LAI}}, \quad (4) $$

where r<sub>c</sub> is the stomatal resistance (s·m<sup>-1</sup>), r<sub>l</sub> is the stomatal resistance (s·m<sup>-1</sup>) of a single leaf, E<sub>t</sub> is the maximum transpiration (mm·d<sup>-1</sup>), λ is the latent heat of vaporization (MJ·kg<sup>-1</sup>), Δ is the slope of the saturation vapor pressure-temperature curve (kPa·°C<sup>-1</sup>), H<sub>net</sub> is the net radiation (MJ·m<sup>-2</sup>·d<sup>-1</sup>), G is the heat flux density to the ground (MJ·m<sup>-2</sup>·d<sup>-1</sup>), γ is the psychometric constant (kPa·°C<sup>-1</sup>), ρ<sub>air</sub> is the air density (kg·m<sup>-3</sup>), P is the air pressure (kPa), K<sub>1</sub> is a dimension coefficient (K<sub>1</sub> = 86 400 s·d<sup>-1</sup>), e<sub>z</sub><sup>o</sup> is the saturation vapor pressure (kPa), e<sub>z</sub> is the vapor pressure of air (kPa), and r<sub>a</sub> is the aerodynamic resistance (s·m<sup>-1</sup>).

Actual evapotranspiration is adjusted from this PET through a stress coefficient calculated as follows:

$$ K_s = \begin{cases} \exp\left(5\left(\frac{\text{SWC}}{0.25 \cdot \text{Zr} \cdot \text{AWC}} - 1\right)\right) & \text{if } \frac{\text{SWC}}{0.25 \cdot \text{Zr} \cdot \text{AWC}} < 1, \\ 1 & \text{else.} \end{cases}, \quad (5) $$

where SWC is the soil water content (mm), and Zr is the root depth (mm).

## 2.4.2 SWAT-NDVI: growth of plants and ET driven by optical remote sensing

A normalized difference vegetation index (NDVI) is calculated from red and near-infrared (NIR) bands (fourth and eighth bands) of Sentinel-2 and is used as a proxy for vegetation growth:

$$ \text{NDVI} = \frac{\text{NIR} - \text{Red}}{\text{NIR} + \text{Red}}, \quad (6) $$

The bare-soil NDVI is lower than 0.2, whereas the crop NDVI ranges from 0.2 at the early stages to 0.7-0.9 at full development. Table 1 shows the rules used to detect crop emergence with NDVI.

Actual evapotranspiration (E<sub>t</sub>, mm) is calculated with the FAO-56 dual-coefficient method (Allen et al., 1998):

$$ E_t = (K_{cb} \cdot K_s + K_e) \cdot ET_0, \quad (7) $$

with ET<sub>0</sub> being the evaporative demand (mm), K<sub>cb</sub> being the crop basal coefficient, K<sub>e</sub> being the soil evaporative coefficient, and K<sub>s</sub> being the stress coefficient. K<sub>cb</sub> can be seen as a linear function of NDVI, the coefficients of which were fixed following Toureiro et al. (2017):

$$ K_{cb} = 1.464 \cdot \text{NDVI} - 0.253. \quad (8) $$

The K<sub>e</sub> factor depends on the crop coefficient K<sub>cb</sub>, the fraction of soil coverage F<sub>cover</sub>, and a reduction coefficient K<sub>r</sub>:

$$ K_e = \min(K_r \cdot (1.2 - K_{cb}); 1.2 \cdot (1 - F_{\text{cover}})), \quad (9) $$

with

$$ F_{\text{cover}} = 1.23 \cdot \text{NDVI} - 0.15, \quad (10) $$

and K<sub>r</sub> depends on the soil texture and SWC according to Merlin et al. (2011).

The stress coefficient is calculated as

$$ K_s = \frac{\text{SWC}}{(1 - p) \cdot \text{Zr} \cdot \text{AWC}}, \quad (11) $$

where p is the depletion coefficient that is usually equal to 0.55 in the FAO-56.
<table><caption>**Table 1.** Decision rules for corn emergence. The emergence occurs on the first day that all the conditions are met. YHU: total heat units of the year; SWC: soil water content; FC: field capacity.</caption><thead><tr><th>SWAT-O</th><th>SWAT-NDVI</th></tr></thead><tbody><tr><td>HU<sub>0</sub> &gt; 30 % YHU</td><td>20 April &lt; day &lt; 10 July</td></tr><tr><td>SWC &lt; 1.05 · FC</td><td>NDVI increasing<br>0.2 &lt; NDVI &lt; 0.5</td></tr></tbody></table>

## 2.4.1 SWAT-O: growth of plants and ET driven by heat units (EPIC)

Heat units (namely the average temperature of the day minus the base temperature) drive the growth of LAI (leaf area index) according to the following formulas (Barnard, 1948; Phillips, 1950):

$$ \text{fr}_{\text{LAlmx}} = \frac{\text{fr}_{\text{PHU}}}{\text{fr}_{\text{PHU}} + \exp(l_1 - l_2 \cdot \text{fr}_{\text{PHU}})}, \quad (1) $$

$$ \text{fr}_{\text{PHU}} = \frac{\sum_{i=1}^{d} \text{HU}_i}{\text{PHU}}, \quad (2) $$

where $fr_{LAlmx}$ is the fraction of maximum LAI, $l_1$ and $l_2$ are shape coefficients, $HU_i$ is the heat unit of day $i$ (°C), and PHU is the potential heat unit required for plant maturity (°C). A more detailed description of this simplified EPIC module can be found in the SWAT theoretical documentation (Neitsch et al., 2001).

The beginning of growth is settled according to decision rules, where the most important parameter is the sum of heat units since 1 January ($HU_0$). The emergence occurs on the first day that all the conditions are met. See an example for corn in Table 1.

Potential evapotranspiration (PET) of the crop depends on its LAI based on Penman-Monteith equation (Monteith, 1965):

$$ \lambda E_t = \frac{\Delta(H_{\text{net}} - G) + \gamma K_1 \cdot (0.622\lambda \cdot \frac{\rho_{\text{air}}}{P}) \cdot \frac{e_z^o - e_z}{r_a}}{\Delta + \gamma \cdot (1 + r_c/r_a)}, \quad (3) $$

$$ r_c = \frac{r_l}{0.5 \cdot \text{LAI}}, \quad (4) $$

where $r_c$ is the stomatal resistance (s·m⁻¹), $r_l$ is the stomatal resistance (s·m⁻¹) of a single leaf, $E_t$ is the maximum transpiration (mm·d⁻¹), $\lambda$ is the latent heat of vaporization (MJ·kg⁻¹), $\Delta$ is the slope of the saturation vapor pressure-temperature curve (kPa·°C⁻¹), $H_{net}$ is the net radiation (MJ·m⁻²·d⁻¹), $G$ is the heat flux density to the ground (MJ·m⁻²·d⁻¹), $\gamma$ is the psychometric constant (kPa·°C⁻¹), $\rho_{air}$ is the air density (kg·m⁻³), $P$ is the air pressure (kPa), $K_1$ is a dimension coefficient ($K_1 = 86400 \text{ s}^{-1}$), $e_z^o$ is the saturation vapor pressure (kPa), $e_z$ is the vapor pressure of air (kPa), and $r_a$ is the aerodynamic resistance (s·m⁻¹).

Actual evapotranspiration is adjusted from this PET through a stress coefficient calculated as follows:

$$ K_s = \begin{cases} \exp\left(5\left(\frac{\text{SWC}}{0.25 \cdot \text{Zr} \cdot \text{AWC}} - 1\right)\right) & \text{if } \frac{\text{SWC}}{0.25 \cdot \text{Zr} \cdot \text{AWC}} < 1, \\ 1 & \text{else.} \end{cases}, \quad (5) $$

where SWC is the soil water content (mm), and Zr is the root depth (mm).

## 2.4.2 SWAT-NDVI: growth of plants and ET driven by optical remote sensing

A normalized difference vegetation index (NDVI) is calculated from red and near-infrared (NIR) bands (fourth and eighth bands) of Sentinel-2 and is used as a proxy for vegetation growth:

$$ \text{NDVI} = \frac{\text{NIR} - \text{Red}}{\text{NIR} + \text{Red}}, \quad (6) $$

The bare-soil NDVI is lower than 0.2, whereas the crop NDVI ranges from 0.2 at the early stages to 0.7-0.9 at full development. Table 1 shows the rules used to detect crop emergence with NDVI.

Actual evapotranspiration ($E_t$, mm) is calculated with the FAO-56 dual-coefficient method (Allen et al., 1998):

$$ E_t = (K_{cb} \cdot K_s + K_e) \cdot ET_0, \quad (7) $$

with $ET_0$ being the evaporative demand (mm), $K_{cb}$ being the crop basal coefficient, $K_e$ being the soil evaporative coefficient, and $K_s$ being the stress coefficient. $K_{cb}$ can be seen as a linear function of NDVI, the coefficients of which were fixed following Toureiro et al. (2017):

$$ K_{cb} = 1.464 \cdot \text{NDVI} - 0.253. \quad (8) $$

The $K_e$ factor depends on the crop coefficient $K_{cb}$, the fraction of soil coverage $F_{cover}$, and a reduction coefficient $K_r$:

$$ K_e = \min(K_r \cdot (1.2 - K_{cb}); 1.2 \cdot (1 - F_{\text{cover}})), \quad (9) $$

with

$$ F_{\text{cover}} = 1.23 \cdot \text{NDVI} - 0.15, \quad (10) $$

and $K_r$ depends on the soil texture and SWC according to Merlin et al. (2011).

The stress coefficient is calculated as

$$ K_s = \frac{\text{SWC}}{(1 - p) \cdot \text{Zr} \cdot \text{AWC}}, \quad (11) $$

where $p$ is the depletion coefficient that is usually equal to 0.55 in the FAO-56.
## 2.5 Model for irrigation

Based on local water managers' knowledge, all fields with a slope lower than 10% and that are covered by corn, soybean, orchard or vegetables were considered to be irrigated. The actual source of irrigation withdrawals (river or reservoirs) is not known. The PAR database, however, provides information about the distribution of authorized withdrawals between sources. We assumed the real withdrawal distribution to be the same. It results in 58% of the water being pumped from the rivers, whereas 42% comes from small reservoirs (36% from connected reservoirs and 6% from disconnected ones). As small reservoirs (<10 ha) are not included in our setup, irrigation water was considered to come from the closest river in relation to each HRU. This should be kept in mind for the interpretation of results.

Rules are required to trigger the simulated irrigation. In order to implement them, we used several surveys conducted in southwestern France to identify farmers' behaviors (Leenhardt et al., 2004; Maton et al., 2005, 2007; Senthilkumar et al., 2015). Table 2 presents the implemented rules, which include soil water content, the rainfall of the day and the next day, and crop phenological stages.

The $p_{\text{trig}}$ parameter (Table 2) can vary depending on the crop, and it was calibrated for corn and soybean. The calibration method is similar to the one described by Olivera-Guerra et al. (2023). Calibration and validation data are the annual irrigation depth for 160 irrigated fields in the PKGC database. This led to a value of 0.57 for corn and soybean and 0.55 for silage corn.

## 2.6 Calibration of hydrological processes in SWAT

Streamflow in the Gimone River during summer is more impacted by human activities (release of water from dams, withdrawals) than by natural hydrological processes. In order to take this into account, calibration of natural hydrological parameters has been performed for months where anthropogenic impact is reduced: from November to May.

The calibration of parameters related to natural hydrology has been performed only for the SWAT-O setup, and calibrated values were also used for the SWAT-NDVI setup for several reasons. First, Sentinel-2 remote sensing data are only available from 2017. The SWAT-NDVI setup would have not allowed us to cover a sufficiently long period to perform calibration and validation over contrasted hydrological years. This choice has also been made to make a more straightforward comparison between the two crop models. Indeed, if calibration had been performed for both setups, the specificity of each crop model could have been hidden.

The years 2017 to 2021 have been used as the calibration period, and the years 2012 to 2016 have been as the validation period. The years 2010 and 2011 have been run as a warm-up period.

Sensitive parameters have been identified through a sensitivity analysis using a one-at-a-time procedure as described by Abbaspour (2015) and based on a previous study conducted in this watershed (Grusson et al., 2015; Cakir et al., 2020). This allowed us to adjust the range of each parameter before starting the calibration process. The parameters are listed in Table 3.

Calibration has been performed by randomly drawing 1000 sets of parameters following a uniform distribution within the tested range of each parameter. The model was then run for each of these 1000 sets over the 2010–2021 period. For each of the 1000 runs, several metrics have been calculated based on daily flow over the 5-year calibration period, excluding June to October. The metrics are the Nash–Sutcliffe efficiency (Nash and Sutcliffe, 1970) performed for the logarithm (NSElog) and the square root (NSEsqrt) of the discharge, as well as the Kling–Gupta efficiency performed for the square root of the discharge (KGEsqrt) (Gupta et al., 2009). The square root (sqrt) and logarithm (log) transforms, respectively, put the emphasis on middle-range and low flows (Pushpalatha et al., 2012; Santos et al., 2018). For each score independently, the top 1% parameter sets were analyzed using boxplots. The dispersion between those best sets allows us to highlight the equifinality issue during the calibration process.

The goodness of fit of the validation period was then evaluated with the same metrics and additionally with $R^2$ and $P_{\text{bias}}$. Eventually, only NSElog was used to select the parameter values because the calibration months (November to May) still contained low flows but were influenced very little by withdrawals. In order to assess the uncertainty due to calibration, the four best parameter sets according to NSElog were selected for the rest of this study.

# 3 Results and discussion

## 3.1 Calibration and hydrology outside the low-water period

Figure 2 shows the dispersion of the top 1% of sets of parameters according to each criterion independently (NSEsqrt, NSElog, KGEsqrt). Since the bounds of the y axis correspond to the range of each parameter, a low dispersion of values reveals a high sensitivity of the parameter and a low equifinality.

Accordingly, $\alpha_{\text{alpha\_bf}}$ and, to a lesser extent, $C_{n_2}$ are sensible parameters. A relative value of $C_{n_2}$ below 1 means that the runoff-over-infiltration ratio must be lower than the default values. A value of $\alpha_{\text{alpha\_bf}}$ around 0.005 to $0.01 \text{ d}^{-1}$ means that the transfer time from the aquifer to the river is about 3 to 6 months. Although no transfer time measurements are available to check the credibility of this value, the simulated base flow is consistent with the base flow index (BFI) of the river calculated with the Wallingford
## 2.5 Model for irrigation

Based on local water managers' knowledge, all fields with a slope lower than 10% and that are covered by corn, soybean, orchard or vegetables were considered to be irrigated. The actual source of irrigation withdrawals (river or reservoirs) is not known. The PAR database, however, provides information about the distribution of authorized withdrawals between sources. We assumed the real withdrawal distribution to be the same. It results in 58% of the water being pumped from the rivers, whereas 42% comes from small reservoirs (36% from connected reservoirs and 6% from disconnected ones). As small reservoirs (<10 ha) are not included in our setup, irrigation water was considered to come from the closest river in relation to each HRU. This should be kept in mind for the interpretation of results.

Rules are required to trigger the simulated irrigation. In order to implement them, we used several surveys conducted in southwestern France to identify farmers' behaviors (Leenhardt et al., 2004; Maton et al., 2005, 2007; Senthilkumar et al., 2015). Table 2 presents the implemented rules, which include soil water content, the rainfall of the day and the next day, and crop phenological stages.

The $p_{\text{trig}}$ parameter (Table 2) can vary depending on the crop, and it was calibrated for corn and soybean. The calibration method is similar to the one described by Olivera-Guerra et al. (2023). Calibration and validation data are the annual irrigation depth for 160 irrigated fields in the PKGC database. This led to a value of 0.57 for corn and soybean and 0.55 for silage corn.

## 2.6 Calibration of hydrological processes in SWAT

Streamflow in the Gimone River during summer is more impacted by human activities (release of water from dams, withdrawals) than by natural hydrological processes. In order to take this into account, calibration of natural hydrological parameters has been performed for months where anthropogenic impact is reduced: from November to May.

The calibration of parameters related to natural hydrology has been performed only for the SWAT-O setup, and calibrated values were also used for the SWAT-NDVI setup for several reasons. First, Sentinel-2 remote sensing data are only available from 2017. The SWAT-NDVI setup would have not allowed us to cover a sufficiently long period to perform calibration and validation over contrasted hydrological years. This choice has also been made to make a more straightforward comparison between the two crop models. Indeed, if calibration had been performed for both setups, the specificity of each crop model could have been hidden.

The years 2017 to 2021 have been used as the calibration period, and the years 2012 to 2016 have been as the validation period. The years 2010 and 2011 have been run as a warm-up period.

Sensitive parameters have been identified through a sensitivity analysis using a one-at-a-time procedure as described by Abbaspour (2015) and based on a previous study conducted in this watershed (Grusson et al., 2015; Cakir et al., 2020). This allowed us to adjust the range of each parameter before starting the calibration process. The parameters are listed in Table 3.

Calibration has been performed by randomly drawing 1000 sets of parameters following a uniform distribution within the tested range of each parameter. The model was then run for each of these 1000 sets over the 2010–2021 period. For each of the 1000 runs, several metrics have been calculated based on daily flow over the 5-year calibration period, excluding June to October. The metrics are the Nash–Sutcliffe efficiency (Nash and Sutcliffe, 1970) performed for the logarithm (NSElog) and the square root (NSEsqrt) of the discharge, as well as the Kling–Gupta efficiency performed for the square root of the discharge (KGEsqrt) (Gupta et al., 2009). The square root (sqrt) and logarithm (log) transforms, respectively, put the emphasis on middle-range and low flows (Pushpalatha et al., 2012; Santos et al., 2018). For each score independently, the top 1% parameter sets were analyzed using boxplots. The dispersion between those best sets allows us to highlight the equifinality issue during the calibration process.

The goodness of fit of the validation period was then evaluated with the same metrics and additionally with $R^2$ and $P_{\text{bias}}$. Eventually, only NSElog was used to select the parameter values because the calibration months (November to May) still contained low flows but were influenced very little by withdrawals. In order to assess the uncertainty due to calibration, the four best parameter sets according to NSElog were selected for the rest of this study.

# 3 Results and discussion

## 3.1 Calibration and hydrology outside the low-water period

Figure 2 shows the dispersion of the top 1% of sets of parameters according to each criterion independently (NSEsqrt, NSElog, KGEsqrt). Since the bounds of the y axis correspond to the range of each parameter, a low dispersion of values reveals a high sensitivity of the parameter and a low equifinality.

Accordingly, $\alpha_{\text{bf}}$ and, to a lesser extent, $C_n$ are sensible parameters. A relative value of $C_n$ below 1 means that the runoff-over-infiltration ratio must be lower than the default values. A value of $\alpha_{\text{bf}}$ around 0.005 to $0.01 \text{ d}^{-1}$ means that the transfer time from the aquifer to the river is about 3 to 6 months. Although no transfer time measurements are available to check the credibility of this value, the simulated base flow is consistent with the base flow index (BFI) of the river calculated with the Wallingford
<table><caption>**Table 2.** Decision rules for irrigation of corn. SWC: soil water content; frPHU: fraction of YHU.</caption><thead><tr><th>Type of rule</th><th>SWAT-O</th><th>SWAT-NDVI</th></tr></thead><tbody><tr><td>Water stress</td><td>SWC &lt; (1 - p<sub>trig</sub>) · Zr · AWC</td><td></td></tr><tr><td>Weather forecast</td><td>rain(t) &lt; 10 mm and rain(t + 1) &lt; 20 mm</td><td></td></tr><tr><td>Minimum time between two irrigations</td><td>10 d</td><td></td></tr><tr><td>Irrigation depth (starter)</td><td>If 0.1 &lt; frPHU &lt; 0.2: 15 mm</td><td></td></tr><tr><td>Irrigation depth (normal)</td><td>If 0.2 &lt; frPHU: 30 mm</td><td></td></tr><tr><td>Stop irrigation (senescence)</td><td>frPHU &lt; 0.9</td><td>If NDVI is decreasing: NDVI &gt; 0.95 · max (NDVI)</td></tr></tbody></table>

<table><caption>**Table 3.** Parameters selected for calibration, as well as the range within which they can vary.</caption><thead><tr><th>Name</th><th>Description</th><th>Range</th><th>Unit</th></tr></thead><tbody><tr><td>K<sub>sat</sub></td><td>Saturated hydraulic conductivity</td><td>0.4–1.5 K<sub>sat</sub></td><td>mm d<sup>-1</sup></td></tr><tr><td>Cn<sub>2</sub></td><td>Curve number for partition of rain between surface runoff and infiltration</td><td>0.9–1.1 Cn<sub>2</sub></td><td>–</td></tr><tr><td>perco</td><td>Ability of soil water to reach aquifer</td><td>0.5–1</td><td>–</td></tr><tr><td>latq_co</td><td>Adjustment coefficient for lateral subsurface flow</td><td>0–1</td><td>–</td></tr><tr><td>alpha_bf</td><td>Groundwater transfer time to the river</td><td>0.005–0.04</td><td>d<sup>-1</sup></td></tr></tbody></table>

<table><caption>**Table 4.** Four best parameter sets according to NSElog.</caption><thead><tr><th>K<sub>sat</sub> (rel.)</th><th>Cn<sub>2</sub> (rel.)</th><th>perco</th><th>latq_co</th><th>alpha_bf</th></tr></thead><tbody><tr><td>0.78</td><td>0.96</td><td>0.55</td><td>0.12</td><td>0.0057</td></tr><tr><td>0.75</td><td>0.97</td><td>0.53</td><td>0.20</td><td>0.0061</td></tr><tr><td>1.06</td><td>0.91</td><td>0.66</td><td>0.16</td><td>0.0071</td></tr><tr><td>0.65</td><td>0.92</td><td>0.86</td><td>0.36</td><td>0.0067</td></tr></tbody></table>

method (Hamon, 1963; ESPERE software from BRGM: https://www.brgm.fr/fr/logiciel/espere-estimation-pluie-efficace-recharge-selon-differentes-methodes; Lanini et al., 2020).

On the contrary, latq_co and Ksat values seem to be scattered over a wider range. However, further analyses have shown that these parameters are highly correlated because they are both used to calculate the sub-surface lateral flow: if latq_co is low then Ksat must be high and vice versa – hence the conclusion that when the value of one of these two coefficients is fixed then the other becomes very sensitive. Only the perco coefficient seems to have a lower sensitivity.

Due to the dispersion of the best parameters without large differences in the scores, it seemed appropriate to keep more than one set of parameters. This allows us to account for calibration uncertainty in future runs. The four best sets of parameters according to NSElog are therefore kept (Table 4). Table 5 shows the values of metrics during calibration (2017–2021) and validation (2012–2016) periods for the simulation using the first parameter set. It is worth noting that these scores are calculated only over the months of November to May.

The NSE value is only satisfactory (Moriasi et al., 2007), but since NSE mostly assesses the performance for high

flows, which are not the goal of this study, and knowing the limitation of discharge measurement for high flows, these values are seen as sufficient here. On the other hand, the NSElog values and the KGEsqrt values are very good in the calibration period and good in the validation period according to the classification of Moriasi et al. (2007, 2015). This indicates a good agreement for medium and low simulated flow rates compared to measurements. As for R² and Pbias, these values are considered to be good and very good, respectively.

The annual water balance averaged over 10 years (2012–2021) as simulated from the calibrated values is detailed in Table 6. Rainfall (around 700 mm per year) is mostly converted into evapotranspiration (around 580 mm per year) and flow (126 mm per year). In comparison, the measured annual flow is 116 mm per year. In this section of the study, the simulated streamflow is not yet reduced by withdrawals, which can explain most of the difference. This water balance is very similar to the one obtained by Boithias et al. (2014) in a nearby and similar watershed. By not including human influence on the streamflow, they reported a way lower performance in low-flow periods compared to in high-flow periods. This is the reason why our hydrological calibration is performed only during the November–May period, whereas the influenced flow (June–October) performance is investigated in the following sections. Table 6 also shows that two-thirds of the river flow comes from subsurface runoff. Baseflow accounts for about 14 % of the total flow.

The daily calibrated streamflow at the outlet is shown in Fig. 3. The four simulations with the four selected sets of parameters are plotted as a minimum–maximum colored range, along with a plain line for the median value.
<table><caption>**Table 2.** Decision rules for irrigation of corn. SWC: soil water content; frPHU: fraction of YHU.</caption><thead><tr><th>Type of rule</th><th>SWAT-O</th><th>SWAT-NDVI</th></tr></thead><tbody><tr><td>Water stress</td><td>SWC &lt; (1 - p<sub>trig</sub>) · Zr · AWC</td><td></td></tr><tr><td>Weather forecast</td><td>rain(t) &lt; 10 mm and rain(t + 1) &lt; 20 mm</td><td></td></tr><tr><td>Minimum time between two irrigations</td><td>10 d</td><td></td></tr><tr><td>Irrigation depth (starter)</td><td>If 0.1 &lt; frPHU &lt; 0.2: 15 mm</td><td></td></tr><tr><td>Irrigation depth (normal)</td><td>If 0.2 &lt; frPHU: 30 mm</td><td></td></tr><tr><td>Stop irrigation (senescence)</td><td>frPHU &lt; 0.9 If NDVI is decreasing: NDVI &gt; 0.95 · max (NDVI)</td><td></td></tr></tbody></table>

<table><caption>**Table 3.** Parameters selected for calibration, as well as the range within which they can vary.</caption><thead><tr><th>Name</th><th>Description</th><th>Range</th><th>Unit</th></tr></thead><tbody><tr><td>K<sub>sat</sub></td><td>Saturated hydraulic conductivity</td><td>0.4–1.5 K<sub>sat</sub></td><td>mm d<sup>-1</sup></td></tr><tr><td>Cn<sub>2</sub></td><td>Curve number for partition of rain between surface runoff and infiltration</td><td>0.9–1.1 Cn<sub>2</sub></td><td>–</td></tr><tr><td>perco</td><td>Ability of soil water to reach aquifer</td><td>0.5–1</td><td>–</td></tr><tr><td>latq_co</td><td>Adjustment coefficient for lateral subsurface flow</td><td>0–1</td><td>–</td></tr><tr><td>alpha_bf</td><td>Groundwater transfer time to the river</td><td>0.005–0.04</td><td>d<sup>-1</sup></td></tr></tbody></table>

<table><caption>**Table 4.** Four best parameter sets according to NSElog.</caption><thead><tr><th>K<sub>sat</sub> (rel.)</th><th>Cn<sub>2</sub> (rel.)</th><th>perco</th><th>latq_co</th><th>alpha_bf</th></tr></thead><tbody><tr><td>0.78</td><td>0.96</td><td>0.55</td><td>0.12</td><td>0.0057</td></tr><tr><td>0.75</td><td>0.97</td><td>0.53</td><td>0.20</td><td>0.0061</td></tr><tr><td>1.06</td><td>0.91</td><td>0.66</td><td>0.16</td><td>0.0071</td></tr><tr><td>0.65</td><td>0.92</td><td>0.86</td><td>0.36</td><td>0.0067</td></tr></tbody></table>

method (Hamon, 1963; ESPERE software from BRGM: https://www.brgm.fr/fr/logiciel/espere-estimation-pluie-efficace-recharge-selon-differentes-methodes; Lanini et al., 2020).

On the contrary, latq_co and Ksat values seem to be scattered over a wider range. However, further analyses have shown that these parameters are highly correlated because they are both used to calculate the sub-surface lateral flow: if latq_co is low then Ksat must be high and vice versa – hence the conclusion that when the value of one of these two coefficients is fixed then the other becomes very sensitive. Only the perco coefficient seems to have a lower sensitivity.

Due to the dispersion of the best parameters without large differences in the scores, it seemed appropriate to keep more than one set of parameters. This allows us to account for calibration uncertainty in future runs. The four best sets of parameters according to NSElog are therefore kept (Table 4). Table 5 shows the values of metrics during calibration (2017–2021) and validation (2012–2016) periods for the simulation using the first parameter set. It is worth noting that these scores are calculated only over the months of November to May.

The NSE value is only satisfactory (Moriasi et al., 2007), but since NSE mostly assesses the performance for high

flows, which are not the goal of this study, and knowing the limitation of discharge measurement for high flows, these values are seen as sufficient here. On the other hand, the NSElog values and the KGEsqrt values are very good in the calibration period and good in the validation period according to the classification of Moriasi et al. (2007, 2015). This indicates a good agreement for medium and low simulated flow rates compared to measurements. As for R² and Pbias, these values are considered to be good and very good, respectively.

The annual water balance averaged over 10 years (2012–2021) as simulated from the calibrated values is detailed in Table 6. Rainfall (around 700 mm per year) is mostly converted into evapotranspiration (around 580 mm per year) and flow (126 mm per year). In comparison, the measured annual flow is 116 mm per year. In this section of the study, the simulated streamflow is not yet reduced by withdrawals, which can explain most of the difference. This water balance is very similar to the one obtained by Boithias et al. (2014) in a nearby and similar watershed. By not including human influence on the streamflow, they reported a way lower performance in low-flow periods compared to in high-flow periods. This is the reason why our hydrological calibration is performed only during the November–May period, whereas the influenced flow (June–October) performance is investigated in the following sections. Table 6 also shows that two-thirds of the river flow comes from subsurface runoff. Baseflow accounts for about 14 % of the total flow.

The daily calibrated streamflow at the outlet is shown in Fig. 3. The four simulations with the four selected sets of parameters are plotted as a minimum–maximum colored range, along with a plain line for the median value.
<figure><img src="image_2.png" /><figcaption>**Figure 2.** Boxplot analysis of the parameter values of the 10 best simulations according to three different metrics.</figcaption></figure>

<table><caption>**Table 5.** Calibration and validation scores.</caption><thead><tr><th></th><th>NSEsqrt</th><th>NSElog</th><th>NSE</th><th>KGEsqrt</th><th>R<sup>2</sup></th><th>P<sub>bias</sub></th></tr></thead><tbody><tr><td>Calibration (2017–2021)</td><td>0.70</td><td>0.78</td><td>0.55</td><td>0.84</td><td>0.74</td><td>5.7%</td></tr><tr><td>Validation (2012–2016)</td><td>0.65</td><td>0.70</td><td>0.44</td><td>0.82</td><td>0.69</td><td>11.6%</td></tr></tbody></table>

## 3.2 Plant phenology

The emergence of crops follows the rules described in Sect. 2.4. As a result, the emergence dates have very different distributions depending on the method. In the SWAT-O setup, most of the crop growths are synchronized. For instance, in 2017, the model grows all corn during the last week of May. The reason is that the growth in SWAT-O is mostly determined by air temperature, which is nearly homogeneous over the watershed. On the other hand, the use of NDVI to detect the beginning of growth spreads the emergence dates over almost 2 months. Figure 4 shows this effect on the 700 corn fields in the Gimone watershed. To decide which distribution is the more realistic one, those dates have been compared to 150 corn fields from the PKGC database located within 80 km of the study area. The assumption is made that the dynamics of these plots and of the Gimone plots are similar due to the similarity in the altitudes, latitudes and climate. The PKGC database contains the sowing date, whereas the use of NDVI allows the detection of the phenological stages 6 to 8 of leaves (hereafter called emergence for more clarity). For consistency between NDVI and PKGC estimations, a fixed amount of 300 heat units has been added to the database of sowing dates to be converted into plausible emergence dates, 300 °C d⁻¹ being the average needed to reach the phenological stages 6 to 8 of leaves from sowing.

Comparison between the PKGC database and the models clearly shows that spread emergence dates (SWAT-NDVI) fit better to real dates. Surveys conducted in another sub-watershed in southwestern France in 2005 also showed a spread of sowing dates over 5 weeks (Maton et al., 2007). However, the distribution of emergence dates in SWAT-NDVI does not perfectly match the PKGC distribution either.

The major explanation might be the lack of cloud-free Sentinel-2 images during the months of April to June. In 2017, only four dates were usable: 6 April, 16 and 26 May, and 25 June. As images were not available for 40 d, the linear interpolation may lead to unrealistic NDVI time series. In addition, decision rules leads to some artifacts. The peak on 30 April is produced by a lot of plots that already fulfilled the condition of NDVI and that reached the condition of soil water content on the same day in the model. Even if the retrieved emergence dates in the SWAT-NDVI mode were closer to reality, the PKGC date could also not be a representative sample of the diversity of plots. As a consequence, the above comparison focuses on the range of emergence dates without a quantitative description of the distributions.

## 3.3 Evapotranspiration

Figure 5 shows the monthly simulated ET for the years 2017 to 2021. Both setups lead to quite similar ET, with a monthly average difference of 5 % and a daily average difference of 17 % at the watershed scale. On average, over the 5 years, ET simulated with SWAT-NDVI setup was slightly lower than that in the SWAT-O setup (by 1.5 %).

These low differences are in line with the findings of Martin et al. (2016), who also tested the influence of emergence dates retrieved with RS on the water balance. They reported a change in ET from 0 % to ±5 % depending on the months. Their biggest monthly change (around 10 %) occurred in one year where the actual growth of crops was far ahead of the theoretical growth calendar. In the present study, however, no year stands out in particular.
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Boxplot analysis of the parameter values of the 10 best simulations according to three different metrics.</figcaption></figure>

<table><caption>**Table 5.** Calibration and validation scores.</caption><thead><tr><th></th><th>NSEsqrt</th><th>NSElog</th><th>NSE</th><th>KGEsqrt</th><th>R<sup>2</sup></th><th>P<sub>bias</sub></th></tr></thead><tbody><tr><td>Calibration (2017–2021)</td><td>0.70</td><td>0.78</td><td>0.55</td><td>0.84</td><td>0.74</td><td>5.7%</td></tr><tr><td>Validation (2012–2016)</td><td>0.65</td><td>0.70</td><td>0.44</td><td>0.82</td><td>0.69</td><td>11.6%</td></tr></tbody></table>

## 3.2 Plant phenology

The emergence of crops follows the rules described in Sect. 2.4. As a result, the emergence dates have very different distributions depending on the method. In the SWAT-O setup, most of the crop growths are synchronized. For instance, in 2017, the model grows all corn during the last week of May. The reason is that the growth in SWAT-O is mostly determined by air temperature, which is nearly homogeneous over the watershed. On the other hand, the use of NDVI to detect the beginning of growth spreads the emergence dates over almost 2 months. Figure 4 shows this effect on the 700 corn fields in the Gimone watershed. To decide which distribution is the more realistic one, those dates have been compared to 150 corn fields from the PKGC database located within 80 km of the study area. The assumption is made that the dynamics of these plots and of the Gimone plots are similar due to the similarity in the altitudes, latitudes and climate. The PKGC database contains the sowing date, whereas the use of NDVI allows the detection of the phenological stages 6 to 8 of leaves (hereafter called emergence for more clarity). For consistency between NDVI and PKGC estimations, a fixed amount of 300 heat units has been added to the database of sowing dates to be converted into plausible emergence dates, 300 °C d⁻¹ being the average needed to reach the phenological stages 6 to 8 of leaves from sowing.

Comparison between the PKGC database and the models clearly shows that spread emergence dates (SWAT-NDVI) fit better to real dates. Surveys conducted in another sub-watershed in southwestern France in 2005 also showed a spread of sowing dates over 5 weeks (Maton et al., 2007). However, the distribution of emergence dates in SWAT-NDVI does not perfectly match the PKGC distribution either.

The major explanation might be the lack of cloud-free Sentinel-2 images during the months of April to June. In 2017, only four dates were usable: 6 April, 16 and 26 May, and 25 June. As images were not available for 40 d, the linear interpolation may lead to unrealistic NDVI time series. In addition, decision rules leads to some artifacts. The peak on 30 April is produced by a lot of plots that already fulfilled the condition of NDVI and that reached the condition of soil water content on the same day in the model. Even if the retrieved emergence dates in the SWAT-NDVI mode were closer to reality, the PKGC date could also not be a representative sample of the diversity of plots. As a consequence, the above comparison focuses on the range of emergence dates without a quantitative description of the distributions.

## 3.3 Evapotranspiration

Figure 5 shows the monthly simulated ET for the years 2017 to 2021. Both setups lead to quite similar ET, with a monthly average difference of 5 % and a daily average difference of 17 % at the watershed scale. On average, over the 5 years, ET simulated with SWAT-NDVI setup was slightly lower than that in the SWAT-O setup (by 1.5 %).

These low differences are in line with the findings of Martin et al. (2016), who also tested the influence of emergence dates retrieved with RS on the water balance. They reported a change in ET from 0 % to ±5 % depending on the months. Their biggest monthly change (around 10 %) occurred in one year where the actual growth of crops was far ahead of the theoretical growth calendar. In the present study, however, no year stands out in particular.
<table><caption>**Table 6.** Annual water balance from the calibrated SWAT model, averaged over the years 2012–2021, and streamflow separation.</caption><tbody><tr><td>In (mm)</td><td colspan="2">Rainfall input</td><td>695</td></tr><tr><td></td><td>from Neste channel</td><td></td><td>24</td></tr><tr><td rowspan="5">Out (mm)</td><td rowspan="5">Outlet flow</td><td>Surface runoff</td><td>8%</td></tr><tr><td>Sub-surface runoff</td><td>64%</td></tr><tr><td>Baseflow</td><td>14%</td></tr><tr><td>Neste channel</td><td>13%</td></tr><tr><td>Water diverted to Save river</td><td>7</td></tr><tr><td></td><td>Evapotranspiration</td><td>578</td></tr><tr><td></td><td>Deep percolation</td><td>2</td></tr><tr><td rowspan="2">Change in storage (mm)</td><td>Soil water content</td><td></td><td>6</td></tr><tr><td>Rivers + reservoirs + aquifers</td><td></td><td>&lt; 0.5</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 3.** Daily measured and simulated outlet streamflow. Gray areas highlight the periods removed from the calibration-validation process.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 4.** Distribution of corn emergence dates for the year 2017. Comparison between SWAT-O and SWAT-NDVI modes against the reconstructed emergence dates from PKGC database.</figcaption></figure>

## 3.4 Daily withdrawals

Figure 6 provides a daily comparison between the simulated water withdrawals and those that are measured by the networked water meters.

It can be seen that the timing and amounts of withdrawals are highly dependent on the year. The weather might therefore be the first determining factor for irrigation. The rainfall in June and July 2021 resulted in low evaporative demand and very low irrigation during this period. On the contrary, high temperatures and a lack of rain in July 2020 resulted in a high evaporative demand, followed by large amounts of irrigation.

Some differences can also be seen between simulated withdrawals in the SWAT-O and SWAT-NDVI modes. In particular, the peak of water demand in the SWAT-NDVI mode at the beginning of the irrigation season (June) is reduced in terms of intensity and spread in time compared to the SWAT-O mode. This is a direct effect of the spreading of emergence dates, which triggers a spread in crop management dates.

At the end of the 2020 irrigation season (September), the SWAT-NDVI mode produced a higher irrigation amount than the SWAT-O mode. This is due to a difference in crop phenology between the two setups: the NDVI showed that corn and soybean were not yet harvested and that their water needs were still high due to the high ET₀. On the contrary, in the SWAT-O mode, most crops were already harvested.

The remaining gaps between the simulations and measurements could be explained by several factors:

1. The actual irrigated fields are not known. Our choice to irrigate all corn, soybeans, vegetables and orchards in fields with slopes less than 10% is consistent but arbitrary. As several authors (Demarez et al., 2019; Pageot et al., 2020; Puy et al., 2022) highlight, irrigated areas are the main driver of water withdrawals. Indeed, a 10% error in areas logically leads to a 10% error in withdrawals.

2. This graph compares the withdrawals that occur only in the two main rivers. We assumed that it represents 56%
<table><caption>**Table 6.** Annual water balance from the calibrated SWAT model, averaged over the years 2012–2021, and streamflow separation.</caption><tbody><tr><td>In (mm)</td><td colspan="2">Rainfall input</td><td>695</td></tr><tr><td></td><td>from Neste channel</td><td></td><td>24</td></tr><tr><td rowspan="5">Out (mm)</td><td rowspan="5">Outlet flow</td><td>Surface runoff</td><td>8%</td></tr><tr><td>Sub-surface runoff</td><td>64%</td></tr><tr><td>Baseflow</td><td>14%</td></tr><tr><td>Neste channel</td><td>13%</td></tr><tr><td>Water diverted to Save river</td><td>7</td></tr><tr><td></td><td>Evapotranspiration</td><td>578</td></tr><tr><td></td><td>Deep percolation</td><td>2</td></tr><tr><td rowspan="2">Change in storage (mm)</td><td>Soil water content</td><td></td><td>6</td></tr><tr><td>Rivers + reservoirs + aquifers</td><td></td><td>&lt; 0.5</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Daily measured and simulated outlet streamflow. Gray areas highlight the periods removed from the calibration-validation process.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Distribution of corn emergence dates for the year 2017. Comparison between SWAT-O and SWAT-NDVI modes against the reconstructed emergence dates from PKGC database.</figcaption></figure>

## 3.4 Daily withdrawals

Figure 6 provides a daily comparison between the simulated water withdrawals and those that are measured by the networked water meters.

It can be seen that the timing and amounts of withdrawals are highly dependent on the year. The weather might therefore be the first determining factor for irrigation. The rainfall in June and July 2021 resulted in low evaporative demand and very low irrigation during this period. On the contrary, high temperatures and a lack of rain in July 2020 resulted in a high evaporative demand, followed by large amounts of irrigation.

Some differences can also be seen between simulated withdrawals in the SWAT-O and SWAT-NDVI modes. In particular, the peak of water demand in the SWAT-NDVI mode at the beginning of the irrigation season (June) is reduced in terms of intensity and spread in time compared to the SWAT-O mode. This is a direct effect of the spreading of emergence dates, which triggers a spread in crop management dates.

At the end of the 2020 irrigation season (September), the SWAT-NDVI mode produced a higher irrigation amount than the SWAT-O mode. This is due to a difference in crop phenology between the two setups: the NDVI showed that corn and soybean were not yet harvested and that their water needs were still high due to the high ET₀. On the contrary, in the SWAT-O mode, most crops were already harvested.

The remaining gaps between the simulations and measurements could be explained by several factors:

1. The actual irrigated fields are not known. Our choice to irrigate all corn, soybeans, vegetables and orchards in fields with slopes less than 10% is consistent but arbitrary. As several authors (Demarez et al., 2019; Pageot et al., 2020; Puy et al., 2022) highlight, irrigated areas are the main driver of water withdrawals. Indeed, a 10% error in areas logically leads to a 10% error in withdrawals.

2. This graph compares the withdrawals that occur only in the two main rivers. We assumed that it represents 56%
<figure><img src="image_5.png" /><figcaption>**Figure 5.** Monthly evapotranspiration in both setups.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 6.** Simulated and measured daily irrigation withdrawals in the two managed rivers (Gimone and Marcaoué) for the years 2020 and 2021.</figcaption></figure>

of the total withdrawals, but the actual distribution of withdrawal origin might be slightly different.

3. Intrinsic uncertainty in the networked meters is very low, but errors could originate from the fact that only 60% to 75% of the intake points are equipped and that we extrapolated these withdrawal data to the missing intake points.

4. Water restrictions are frequent at the end of the low-water season but are not taken into account in the model. The overestimation in September 2020 can indeed be explained by the drought decree on 29 August 2020, which resulted in a 50% reduction in withdrawals in the main rivers and the prohibition of irrigation from small tributaries. In contrast, 2021 was a wet year, and main rivers were been subjected to any restriction.

5. With April to June being a rainy period, Sentinel-2 images are often corrupted by clouds. During this period, summer crops are sown and begin to grow. For this reason,

the linear interpolation between cloud-free images fails to reproduce the convex shape of the NDVI and leads to an overestimation, which results in an overestimation of the ET. Therefore, irrigation can also be overestimated during this period.

**3.5 Effects of irrigation withdrawals on the river flow**

In addition to the SWAT-O and SWAT-NDVI setups, a third setup was run: SWAT-O with no irrigation. Each setup was run four times with different sets of parameters (see Calibration section). In the following, the error bars or uncertainty margins refer to the range of values obtain with those four runs.

For further quantitative analysis, we distinguish between two periods: on the one hand, the high-flow period from November to May, over which period the calibration of hydrological processes has been performed, and on the other hand, the administrative low-flow period from June to October.
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Monthly evapotranspiration in both setups.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Simulated and measured daily irrigation withdrawals in the two managed rivers (Gimone and Marcaoué) for the years 2020 and 2021.</figcaption></figure>

of the total withdrawals, but the actual distribution of withdrawal origin might be slightly different.

3. Intrinsic uncertainty in the networked meters is very low, but errors could originate from the fact that only 60% to 75% of the intake points are equipped and that we extrapolated these withdrawal data to the missing intake points.

4. Water restrictions are frequent at the end of the low-water season but are not taken into account in the model. The overestimation in September 2020 can indeed be explained by the drought decree on 29 August 2020, which resulted in a 50% reduction in withdrawals in the main rivers and the prohibition of irrigation from small tributaries. In contrast, 2021 was a wet year, and main rivers were been subjected to any restriction.

5. With April to June being a rainy period, Sentinel-2 images are often corrupted by clouds. During this period, summer crops are sown and begin to grow. For this reason,

the linear interpolation between cloud-free images fails to reproduce the convex shape of the NDVI and leads to an overestimation, which results in an overestimation of the ET. Therefore, irrigation can also be overestimated during this period.

## 3.5 Effects of irrigation withdrawals on the river flow

In addition to the SWAT-O and SWAT-NDVI setups, a third setup was run: SWAT-O with no irrigation. Each setup was run four times with different sets of parameters (see Calibration section). In the following, the error bars or uncertainty margins refer to the range of values obtain with those four runs.

For further quantitative analysis, we distinguish between two periods: on the one hand, the high-flow period from November to May, over which period the calibration of hydrological processes has been performed, and on the other hand, the administrative low-flow period from June to October.
ber, where the flow is highly influenced by dam releases and
irrigation withdrawals.

During the high-flow season, most of the metrics are not influenced by the introduction of the NDVI method into the model (Table 7-Up). The Pbias values rise from 5.7% to 11.2%, which is in the range of acceptable values (Moriasi et al., 2007). Since hydrological processes were calibrated without NDVI, a slight decrease in model performance could be expected after adding it.

On the other hand, Table 7 (bottom) shows the metrics calculated over the low-flow season. The scores over this period barely changed, but the significant increase of NSElog from about 0 to close to 0.4 suggests an improvement in the low-flow simulation using the SWAT-NDVI mode.

Figure 7 shows the simulated streamflow at the outlet for the period where Sentinel-2 data are available (2017 to 2021). The graph focuses on the irrigation period, from April to October, where the differences between setups are more likely to be seen, and this is the period of interest regarding the water demand and the strategies for water management.

The graph reflects as far as possible the uncertainties we are aware of. For simulation plots, plain lines are the median value, and the errors bars are minimum and maximum values obtained from the four runs with different parameters sets. In addition, there might be a ±10% error in the measured flow rate values according to the water manager. This uncertainty is included in the graph using a colored range.

The minimal flow rate that has to be maintained (DOE) is also included as a horizontal line.

Significant differences can be observed between the three setups. On the one hand, both models fail to simulate the flow rate at the end of the low-flow seasons for 2018, 2019 and 2020, for which heavy rains occurred in October. The reason for this might be that a lot of small reservoirs – not taken into consideration in this study – are close to being empty at the end of summer and are filled up by the first heavy-rain period (bar plots in Fig. 7), leading to less flow reaching the river.

On the other hand, the comparison with the no-irrigation mode clearly shows the importance of taking into account agricultural practices (through withdrawals here) that account for more than half of the streamflow during summer in such a watershed. It shows that the order of magnitude of simulated withdrawals is correct, which is a second validation of the retrieved withdrawals.

Finally, the improved timing of withdrawals in the SWAT-NDVI setup has some visible effects on the streamflow, mostly at the beginning of irrigation season (June, July). The flow is smoother in this mode because it is subject to fewer peaks in withdrawals. In general, the SWAT-NDVI mode produces a higher discharge than the SWAT-O mode. This is due to the lower ET in the SWAT-NDVI mode compared to in the SWAT-O mode (Fig. 5) in addition to the lower amount of withdrawals (Fig. 6).

## 4 Conclusions and outlooks

This study deals with the modeling of crop phenology and irrigation withdrawals in an agricultural zone and their impact on downstream flows. The SWAT+ agro-hydrological model was chosen because it allows a spatialized and detailed description of crop management. Two modeling approaches were compared for crop management. The first modeling approach (SWAT-O) uses heat units to compute the emergence and senescence dates of crops, their growth rate and LAI, and then their ET. The second approach (SWAT-NDVI) uses Sentinel-2 remote sensing data, through the NDVI vegetation index, to determine all these features.

Before tackling the crop modeling, the calibration of hydrological fluxes was addressed. Because of the very influenced streamflow during summer, calibration was performed during the remaining months, from November to May. This uncommon calibration method relies on the assumption that parameters are independent of the time of the year, and this hypothesis has proven to be plausible in the course of this study. A strength of this study is admitting to the uncertainties due to calibration, taking this into account in the analysis.

Then, focusing on crops, the SWAT-NDVI setup first allows us to better account for the spatial and temporal variability of the cultivated crops by spreading the emergence over time and over the watershed, which is more realistic. The modeled ET in this new setup follows the observed crop dynamics at the plot scale, even though the differences in ET seem to be barely significant at the watershed scale. The simulated daily withdrawals have been compared against daily irrigation data provided by the water manager. Both setups reproduce quite well the overall evolution of withdrawals throughout the year. However, the SWAT-NDVI mode shows the best performances at the beginning of the irrigation season without further calibration. This improvement can mainly be attributed to the spreading of crop management dates.

In our watershed, irrigation withdrawals and dam releases are the main determining factors for the river flow in summer. Since dam releases are known and irrigation withdrawals are simulated, the model is able to provide satisfactory daily streamflows in low-flow periods. The results suggest that the SWAT-NDVI mode allows for slightly better accuracy for very low flows, but this would need to be confirmed by other experiments.

To conclude, this study shows that it is possible to simulate a very influenced streamflow with modeling tools. Indeed, very few modelers could focus on daily flow during low-flow periods due to the difficulties in modeling both hydrology and anthropogenic processes. In this study, we assessed two methods for crop dynamics that allow us to retrieve withdrawals and flow dynamics during irrigation periods with an unusual accuracy.
ber, where the flow is highly influenced by dam releases and
irrigation withdrawals.

During the high-flow season, most of the metrics are not influenced by the introduction of the NDVI method into the model (Table 7-Up). The Pbias values rise from 5.7% to 11.2%, which is in the range of acceptable values (Moriasi et al., 2007). Since hydrological processes were calibrated without NDVI, a slight decrease in model performance could be expected after adding it.

On the other hand, Table 7 (bottom) shows the metrics cal-
culated over the low-flow season. The scores over this period
barely changed, but the significant increase of NSElog from
about 0 to close to 0.4 suggests an improvement in the low-
flow simulation using the SWAT-NDVI mode.

Figure 7 shows the simulated streamflow at the outlet for the period where Sentinel-2 data are available (2017 to 2021). The graph focuses on the irrigation period, from April to October, where the differences between setups are more likely to be seen, and this is the period of interest regarding the water demand and the strategies for water management.

The graph reflects as far as possible the uncertainties we
are aware of. For simulation plots, plain lines are the median
value, and the errors bars are minimum and maximum values
obtained from the four runs with different parameters sets. In
addition, there might be a ±10% error in the measured flow
rate values according to the water manager. This uncertainty
is included in the graph using a colored range.

The minimal flow rate that has to be maintained (DOE) is
also included as a horizontal line.

Significant differences can be observed between the three
setups. On the one hand, both models fail to simulate the flow
rate at the end of the low-flow seasons for 2018, 2019 and
2020, for which heavy rains occurred in October. The reason
for this might be that a lot of small reservoirs – not taken into
consideration in this study – are close to being empty at the
end of summer and are filled up by the first heavy-rain period
(bar plots in Fig. 7), leading to less flow reaching the river.

On the other hand, the comparison with the no-irrigation
mode clearly shows the importance of taking into account
agricultural practices (through withdrawals here) that ac-
count for more than half of the streamflow during summer
in such a watershed. It shows that the order of magnitude of
simulated withdrawals is correct, which is a second valida-
tion of the retrieved withdrawals.

Finally, the improved timing of withdrawals in the SWAT-
NDVI setup has some visible effects on the streamflow,
mostly at the beginning of irrigation season (June, July). The
flow is smoother in this mode because it is subject to fewer
peaks in withdrawals. In general, the SWAT-NDVI mode pro-
duces a higher discharge than the SWAT-O mode. This is due
to the lower ET in the SWAT-NDVI mode compared to in the
SWAT-O mode (Fig. 5) in addition to the lower amount of
withdrawals (Fig. 6).

4 Conclusions and outlooks

This study deals with the modeling of crop phenology and
irrigation withdrawals in an agricultural zone and their im-
pact on downstream flows. The SWAT+ agro-hydrological
model was chosen because it allows a spatialized and detailed
description of crop management. Two modeling approaches
were compared for crop management. The first modeling ap-
proach (SWAT-O) uses heat units to compute the emergence
and senescence dates of crops, their growth rate and LAI,
and then their ET. The second approach (SWAT-NDVI) uses
Sentinel-2 remote sensing data, through the NDVI vegetation
index, to determine all these features.

Before tackling the crop modeling, the calibration of hy-
drological fluxes was addressed. Because of the very influ-
enced streamflow during summer, calibration was performed
during the remaining months, from November to May. This
uncommon calibration method relies on the assumption that
parameters are independent of the time of the year, and this
hypothesis has proven to be plausible in the course of this
study. A strength of this study is admitting to the uncertain-
ties due to calibration, taking this it into account in the anal-
ysis.

Then, focusing on crops, the SWAT-NDVI setup first al-
lows us to better account for the spatial and temporal vari-
ability of the cultivated crops by spreading the emergence
over time and over the watershed, which is more realistic.
The modeled ET in this new setup follows the observed
crop dynamics at the plot scale, even though the differ-
ences in ET seem to be barely significant at the watershed
scale. The simulated daily withdrawals have been compared
against daily irrigation data provided by the water manager.
Both setups reproduce quite well the overall evolution of
withdrawals throughout the year. However, the SWAT-NDVI
mode shows the best performances at the beginning of the
irrigation season without further calibration. This improve-
ment can mainly be attributed to the spreading of crop man-
agement dates.

In our watershed, irrigation withdrawals and dam releases
are the main determining factors for the river flow in summer.
Since dam releases are known and irrigation withdrawals are
simulated, the model is able to provide satisfactory daily
streamflows in low-flow periods. The results suggest that the
SWAT-NDVI mode allows for slightly better accuracy for
very low flows, but this would need to be confirmed by other
experiments.

To conclude, this study shows that it is possible to simulate
a very influenced streamflow with modeling tools. Indeed,
very few modelers could focus on daily flow during low-flow
periods due to the difficulties in modeling both hydrology
and anthropogenic processes. In this study, we assessed two
methods for crop dynamics that allow us to retrieve with-
drawals and flow dynamics during irrigation periods with an
unusual accuracy.
<table><caption>**Table 7.** Scores for the SWAT-O and SWAT-NDVI setups for the years 2017 to 2021. At the top is the high-flow period from November to May. At the bottom is the low-flow period from June to October.</caption><thead><tr><th colspan="7">High flow: November to May</th></tr><tr><th></th><th>NSEsqrt</th><th>NSElog</th><th>NSE</th><th>KGEsqrt</th><th>R<sup>2</sup></th><th>P<sub>bias</sub></th></tr></thead><tbody><tr><th scope="row">SWAT-O</th><td>0.70</td><td>0.78</td><td>0.55</td><td>0.84</td><td>0.74</td><td>5.7 %</td></tr><tr><th scope="row">SWAT-NDVI</th><td>0.70</td><td>0.79</td><td>0.55</td><td>0.82</td><td>0.74</td><td>11.2 %</td></tr></tbody><tbody><tr><th colspan="7">Low flow: June to October</th></tr><tr><th></th><th>NSEsqrt</th><th>NSElog</th><th>NSE</th><th>KGEsqrt</th><th>R<sup>2</sup></th><th>P<sub>bias</sub></th></tr><tr><th scope="row">SWAT-O</th><td>0.63</td><td>-0.07</td><td>0.60</td><td>0.78</td><td>0.79</td><td>1.1 %</td></tr><tr><th scope="row">SWAT-NDVI</th><td>0.63</td><td>0.37</td><td>0.61</td><td>0.78</td><td>0.79</td><td>10.5 %</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 7.** Simulation against measured daily outlet streamflow during the June–October period for 5 consecutive years. The SWAT-O and SWAT-NDVI setups are compared with the theoretical simulation without irrigation withdrawals.</figcaption></figure>

It would of course be interesting to reproduce this experiment in different watersheds with different climate and agricultural practices. In areas where land use, crop dynamics and irrigated areas are not well known, remote sensing methods could be of great help. In this study, the hypothesis about

irrigated areas appeared to be sufficient. However, in other watersheds where fewer plots are irrigated, maps of irrigated areas could be necessary (Pageot et al., 2020). In addition, in order to obtain more accurate emergence dates, radar images
<table><caption>**Table 7.** Scores for the SWAT-O and SWAT-NDVI setups for the years 2017 to 2021. At the top is the high-flow period from November to May. At the bottom is the low-flow period from June to October.</caption><thead><tr><th colspan="7">High flow: November to May</th></tr><tr><th></th><th>NSEsqrt</th><th>NSElog</th><th>NSE</th><th>KGEsqrt</th><th>R<sup>2</sup></th><th>P<sub>bias</sub></th></tr></thead><tbody><tr><td>SWAT-O</td><td>0.70</td><td>0.78</td><td>0.55</td><td>0.84</td><td>0.74</td><td>5.7 %</td></tr><tr><td>SWAT-NDVI</td><td>0.70</td><td>0.79</td><td>0.55</td><td>0.82</td><td>0.74</td><td>11.2 %</td></tr></tbody><thead><tr><th colspan="7">Low flow: June to October</th></tr><tr><th></th><th>NSEsqrt</th><th>NSElog</th><th>NSE</th><th>KGEsqrt</th><th>R<sup>2</sup></th><th>P<sub>bias</sub></th></tr></thead><tbody><tr><td>SWAT-O</td><td>0.63</td><td>-0.07</td><td>0.60</td><td>0.78</td><td>0.79</td><td>1.1 %</td></tr><tr><td>SWAT-NDVI</td><td>0.63</td><td>0.37</td><td>0.61</td><td>0.78</td><td>0.79</td><td>10.5 %</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>**Figure 7.** Simulation against measured daily outlet streamflow during the June–October period for 5 consecutive years. The SWAT-O and SWAT-NDVI setups are compared with the theoretical simulation without irrigation withdrawals.</figcaption></figure>

It would of course be interesting to reproduce this experiment in different watersheds with different climate and agricultural practices. In areas where land use, crop dynamics and irrigated areas are not well known, remote sensing methods could be of great help. In this study, the hypothesis about

irrigated areas appeared to be sufficient. However, in other watersheds where fewer plots are irrigated, maps of irrigated areas could be necessary (Pageot et al., 2020). In addition, in order to obtain more accurate emergence dates, radar images
that are not affected by clouds could be used in combination
with optical images (Rolle et al., 2022).

In most watersheds, the lack of knowledge about the ori-
gin of withdrawals still remains a limitation of our approach.
Finally, a major unknown of many watersheds is the small
reservoirs that may have a huge impact on the flow rate
but about which very little information is available (Boisson
et al., 2022).

*Code availability.* The custom SWAT+ sources are available at https://github.com/ElisabethJustin/SWATplus-NDVI (Brochet, 2022).

*Data availability.* Most of the datasets used in this study are free and available online: https://theia.cnes.fr/atdistrib/rocket/#/search?collection=SENTINEL2 (ESA, 2012), https://doi.org/10.5066/F7PR7TFT (USGS, 2018), https://www.gissol.fr/donnees/liens-vers-les-referentiels-regionaux-pedologiques-5634 (GisSol, 2014), https://www.data.gouv.fr/fr/datasets/registre-parcellaire-graphique-rpg-contours-des-parcelles-et-ilots-culturaux-et-leur-groupe-de-cultures-majoritaire/ (IGN, 2021), https://doi.org/10.5281/zenodo.3613415 (Inglada et al., 2018) and https://agreste.agriculture.gouv.fr/agreste-web/disaron/Chd2009/detail/ (Agreste, 2020). Only two datasets are not freely available.

- The meteorological SAFRAN dataset is subject to restric-
tion and needs to be requested from the French mete-
orological agency: https://publitheque.meteo.fr/okapi/accueil/
okapiWebPubli/index.jsp (last access: 5 January 2022)

- Discharge and dam management data have been provided by the local water management company (CACG) and are subject to confidential restrictions (contact: https://www.cacg.fr/, last access: 5 December 2023).

*Author contributions.* EB, VD and SS designed the experiments; EB conducted the experiments; EB, VD, YG and SS analyzed the results; LL provided the data and gave expert assessment on the results. All the authors revised the paper.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

*Acknowledgements.* The authors would like to thank Sylvain Pu-
jol, Nicolas Barry and Damien Lilas from CACG for providing in-

formation about the Gimone watershed, as well as measurements
of streamflow and withdrawals. The authors also thank Vincent
Bustillo (Cesbio) for providing measurements of soil hydraulic con-
ductivity and Vincent Rivalland (Cesbio) for his help with the mod-
eling tools.

*Review statement.* This paper was edited by Ann van Griensven and reviewed by two anonymous referees.

References

Abbastpour, K. C.: SWAT-CUP: SWAT Calibration and Uncertainty Programs – A User Manual, Swiss Federal Institute of Aquatic Science and Technology, Eawag, Dubendorf, Switzerland, https://sndl.ucmerced.edu/files/San_Joaquin/Model_Work/SWAT_MercedRiver/SWATCUP/Usermanual_Swat_Cup_2012.pdf (last access: 10 October 2022), 2015.

Agreste: Enquête sur les pratiques culturales en grandes cultures et prairies 2017, https://agreste.agriculture.gouv.fr/agreste-web/disaron/Chd2009/detail/ (last access: 1 May 2021), 2020.

Allain, S., Ndong, G. O., Lardy, R., and Leenhardt, D.: 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, 2018.

Allen, R. G., Pereira, L. S., Raes, D., and Smith, M.: Crop Evapo-
transpiration: guidelines for computing crop water requirements,
FAO Irrigation and Drainage Paper, article no. 56, ISBN 92-5-
104219-5, 1998.

Arnold, J. G., Allen, P. M., and Bernhardt, G.: A comprehensive surface-groundwater flow model, J. Hydrol., 142, 47–69, https://doi.org/10.1016/0022-1694(93)90004-S, 1993.

Arnold, J. G., Moriasi, D. N., Gassman, P. W., Abbaspour, K. C., White, M. J., Srinivasan, R., Santhi, C., Harmel, R. D., Van Griensven, A., and Van Liew, M. W.: SWAT: Model use, calibration, and validation, T. ASABE, 55, 1491–1508, 2012.

Barnard, J. D.: Heat units as a measure of canning crop maturity, The Canner, 106, 28–29, 1948.

Battude, M., Al Bitar, A., Brut, A., Tallec, T., Huc, M., Cros, J., Weber, J.-J., Lhuissier, L., Simonneaux, V., and Demarez, V.: Modeling water needs and total irrigation depths of maize crop in the south west of France using high spatial and temporal resolution satellite imagery, Agr. Water Manage., 189, 123–136, https://doi.org/10.1016/j.agwat.2017.04.018, 2017.

Bazzi, H., Baghdadi, N., Fayad, I., Zribi, M., Demarez, V., Pageot, Y., and Belhouchette, H.: Detecting Irrigation Events Using Sentinel-1 Data, in: 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, 6355–6358, https://doi.org/10.1109/IGARSS47720.2021.9553587, 2021.

Boisson, A., Villesseche, D., Selles, A., Alazard, M., Chandra, S., Ferrant, S., and Maréchal, J.-C.: Long term monitoring of rainwater harvesting tanks: Is multi-years management possible in crystalline South Indian aquifers?, Hydrol. Process., 36, e14759, https://doi.org/10.1002/hyp.14759, 2022.

Boithias, L., Srinivasan, R., Sauvage, S., Macary, F., and Sánchez-Pérez, J. M.: Daily nitrate losses: Implication on long-term river quality in an intensive agricultural catchment of southwestern France, J. Environ. Qual., 43, 46–54, 2014.
that are not affected by clouds could be used in combination
with optical images (Rolle et al., 2022).

In most watersheds, the lack of knowledge about the ori-
gin of withdrawals still remains a limitation of our approach.
Finally, a major unknown of many watersheds is the small
reservoirs that may have a huge impact on the flow rate
but about which very little information is available (Boisson
et al., 2022).

*Code availability.* The custom SWAT+ sources are available at https://github.com/ElisabethJustin/SWATplus-NDVI (Brochet, 2022).

*Data availability.* Most of the datasets used in this study are free and available online: https://theia.cnes.fr/atdistrib/rocket/#/search?collection=SENTINEL2 (ESA, 2012), https://doi.org/10.5066/F7PR7TFT (USGS, 2018), https://www.gissol.fr/donnees/liens-vers-les-referentiels-regionaux-pedologiques-5634 (GisSol, 2014), https://www.data.gouv.fr/fr/datasets/registre-parcellaire-graphique-rpg-contours-des-parcelles-et-ilots-culturaux-et-leur-groupe-de-cultures-majoritaire/ (IGN, 2021), https://doi.org/10.5281/zenodo.3613415 (Inglada et al., 2018) and https://agreste.agriculture.gouv.fr/agreste-web/disaron/Chd2009/detail/ (Agreste, 2020). Only two datasets are not freely available.

- The meteorological SAFRAN dataset is subject to restric-
tion and needs to be requested from the French mete-
orological agency: https://publitheque.meteo.fr/okapi/accueil/
okapiWebPubli/index.jsp (last access: 5 January 2022)

- Discharge and dam management data have been provided by the local water management company (CACG) and are subject to confidential restrictions (contact: https://www.cacg.fr/, last access: 5 December 2023).

*Author contributions.* EB, VD and SS designed the experiments; EB conducted the experiments; EB, VD, YG and SS analyzed the results; LL provided the data and gave expert assessment on the results. All the authors revised the paper.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

*Acknowledgements.* The authors would like to thank Sylvain Pu-
jol, Nicolas Barry and Damien Lilas from CACG for providing in-

formation about the Gimone watershed, as well as measurements
of streamflow and withdrawals. The authors also thank Vincent
Bustillo (Cesbio) for providing measurements of soil hydraulic con-
ductivity and Vincent Rivalland (Cesbio) for his help with the mod-
eling tools.

*Review statement.* This paper was edited by Ann van Griensven and reviewed by two anonymous referees.

References

Abbastpour, K. C.: SWAT-CUP: SWAT Calibration and Uncertainty Programs – A User Manual, Swiss Federal Institute of Aquatic Science and Technology, Eawag, Dubendorf, Switzerland, https://sndl.ucmerced.edu/files/San_Joaquin/Model_Work/SWAT_MercedRiver/SWATCUP/Usermanual_Swat_Cup_2012.pdf (last access: 10 October 2022), 2015.

Agreste: Enquête sur les pratiques culturales en grandes cultures et prairies 2017, https://agreste.agriculture.gouv.fr/agreste-web/disaron/Chd2009/detail/ (last access: 1 May 2021), 2020.

Allain, S., Ndong, G. O., Lardy, R., and Leenhardt, D.: 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, 2018.

Allen, R. G., Pereira, L. S., Raes, D., and Smith, M.: Crop Evapo-
transpiration: guidelines for computing crop water requirements,
FAO Irrigation and Drainage Paper, article no. 56, ISBN 92-5-
104219-5, 1998.

Arnold, J. G., Allen, P. M., and Bernhardt, G.: A comprehensive surface-groundwater flow model, J. Hydrol., 142, 47–69, https://doi.org/10.1016/0022-1694(93)90004-S, 1993.

Arnold, J. G., Moriasi, D. N., Gassman, P. W., Abbaspour, K. C., White, M. J., Srinivasan, R., Santhi, C., Harmel, R. D., Van Griensven, A., and Van Liew, M. W.: SWAT: Model use, calibration, and validation, T. ASABE, 55, 1491–1508, 2012.

Barnard, J. D.: Heat units as a measure of canning crop maturity, The Canner, 106, 28–29, 1948.

Battude, M., Al Bitar, A., Brut, A., Tallec, T., Huc, M., Cros, J., Weber, J.-J., Lhuissier, L., Simonneaux, V., and Demarez, V.: Modeling water needs and total irrigation depths of maize crop in the south west of France using high spatial and temporal resolution satellite imagery, Agr. Water Manage., 189, 123–136, https://doi.org/10.1016/j.agwat.2017.04.018, 2017.

Bazzi, H., Baghdadi, N., Fayad, I., Zribi, M., Demarez, V., Pageot, Y., and Belhouchette, H.: Detecting Irrigation Events Using Sentinel-1 Data, in: 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, 6355–6358, https://doi.org/10.1109/IGARSS47720.2021.9553587, 2021.

Boisson, A., Villesseche, D., Selles, A., Alazard, M., Chandra, S., Ferrant, S., and Maréchal, J.-C.: Long term monitoring of rainwater harvesting tanks: Is multi-years management possible in crystalline South Indian aquifers?, Hydrol. Process., 36, e14759, https://doi.org/10.1002/hyp.14759, 2022.

Boithias, L., Srinivasan, R., Sauvage, S., Macary, F., and Sánchez-Pérez, J. M.: Daily nitrate losses: Implication on long-term river quality in an intensive agricultural catchment of southwestern France, J. Environ. Qual., 43, 46–54, 2014.
Brochet, E.: Modified SWAT+ NDVI, GitHub [code], https://github.com/ElisabethJustin/SWATplus-NDVI, 2022.

Bruand, A., Fernández, P. P., and Duval, O.: Use of class pedo-transfer functions based on texture and bulk density of clods to generate water retention curves, Soil Use Manage., 19, 232–242, https://doi.org/10.1111/j.1475-2743.2003.tb00309.x, 2003.

Cakir, R., Raimonet, M., Sauvage, S., Paredes-Arquiola, J., Grusson, Y., Roset, L., Meaurio, M., Navarro, E., Sevilla-Callejo, M., Lechuga-Crespo, J. L., Gomiz Pascual, J. J., Bodoque, J. M., and Sánchez-Pérez, J. M.: Hydrological Alteration Index as an Indicator of the Calibration Complexity of Water Quantity and Quality Modeling in the Context of Global Change, Water, 12, 115, https://doi.org/10.3390/w12010115, 2020.

Carroget, A., Perrin, C., Sauquet, E., Vidal, J.-P., Chazot, S., Chauveau, M., and Rouchy, N.: Explore 2070: quelle utilisation d'un exercice prospectif sur les impacts des changements climatiques à l'échelle nationale pour définir des stratégies d'adaptation, Sciences Eaux & Territoires, p. 4, https://doi.org/10.14758/SET-REVUE.2017.22.02, 2017.

Clavel, L., Soudais, J., Baudet, D., and Leenhardt, D.: Integrating expert knowledge and quantitative information for mapping cropping systems, Land Use Policy, 28, 57–65, https://doi.org/10.1016/j.landusepol.2010.05.001, 2011.

Courault, D., Hadria, R., Ruget, F., Olioso, A., Duchemin, B., Hagolle, O., and Dedieu, G.: Combined use of FORMOSAT-2 images with a crop model for biomass and water monitoring of permanent grassland in Mediterranean region, Hydrol. Earth Syst. Sci., 14, 1731–1744, https://doi.org/10.5194/hess-14-1731-2010, 2010.

Demarez, V., Helen, F., Marais-Sicre, C., and Baup, F.: In-Season Mapping of Irrigated Crops Using Landsat 8 and Sentinel-1 Time Series, Remote Sens., 11, 118, https://doi.org/10.3390/rs11020118, 2019.

Drusch, M., Del Bello, U., Carlier, S., Colin, O., Fernandez, V., Gascon, F., Hoersch, B., Isola, C., Laberinti, P., Martimort, P., Meygret, A., Spoto, F., Sy, O., Marchese, F., and Bargellini, P.: Sentinel-2: ESA's Optical High-Resolution Mission for GMES Operational Services, Remote Sens. Environ., 120, 25–36, https://doi.org/10.1016/j.rse.2011.11.026, 2012.

Durand, Y., Brun, E., Merindol, L., Guyomarc'h, G., Lesafre, B., and Martin, E.: A meteorological estimation of relevant parameters for snow models, Ann. Glaciol., 18, 65–71, https://doi.org/10.3189/S0260305500011277, 1993.

ESA: SENTINEL-2, ESA's Optical High-Resolution, Theia [data set], https://theia.cnes.fr/atdistrib/rocket/#/search?collection=SENTINEL2, 2012.

Etchanchu, J., Rivalland, V., Gascoin, S., Cros, J., Tallec, T., Brut, A., and Boulet, G.: Effects of high spatial and temporal resolution Earth observations on simulated hydrometeorological variables in a cropland (southwestern France), Hydrol. Earth Syst. Sci., 21, 5693–5708, https://doi.org/10.5194/hess-21-5693-2017, 2017.

Ferrant, S., Gascoin, S., Veloso, A., Salmon-Monviola, J., Claverie, M., Rivalland, V., Dedieu, G., Demarez, V., Ceschia, E., Probst, J.-L., Durand, P., and Bustillo, V.: Agro-hydrology and multi-temporal high-resolution remote sensing: toward an explicit spatial processes calibration, Hydrol. Earth Syst. Sci., 18, 5219–5237, https://doi.org/10.5194/hess-18-5219-2014, 2014.

Fohrer, N., Dietrich, A., Kolychalow, O., and Ulrich, U.: Assessment of the Environmental Fate of the Herbicides Flufenacet and

Metazachlor with the SWAT Model, J. Environ. Qual., 43, 75–85, https://doi.org/10.2134/jeq2011.0382, 2014.

GisSol: Référentiel Régional Pédologique, Gis-Sol [data set], https://www.gissol.fr/donnees/liens-vers-les-referentiels-regionaux-pedologiques-5634, 2014.

Grusson, Y., Sun, X., Gascoin, S., Sauvage, S., Raghavan, S., Antil, F., and Sáchez-Pérez, J.-M.: Assessing the capability of the SWAT model to simulate snow, snow melt and streamflow dynamics over an alpine watershed, J. Hydrol., 531, 574–588, https://doi.org/10.1016/j.jhydrol.2015.10.070, 2015.

Gupta, H. V., Kling, H., Yilmaz, K. K., and Martinez, G. F.: Decomposition of the mean squared error and NSE performance criteria: Implications for improving hydrological modelling, J. Hydrol., 377, 80–91, https://doi.org/10.1016/j.jhydrol.2009.08.003, 2009.

Guérif, M. and Duke, C.: Calibration of the SUCROS emergence and early growth module for sugar beet using optical remote sensing data assimilation, Eur. J. Agron., 9, 127–136, https://doi.org/10.1016/S1161-0301(98)00031-8, 1998.

Hagolle, O., Huc, M., Villa Pascual, D., and Dedieu, G.: A Multi-Temporal and Multi-Spectral Method to Estimate Aerosol Optical Thickness over Land, for the Atmospheric Correction of FormoSat-2, LandSat, VENUS and Sentinel-2 Images, Remote Sens., 7, 2668–2691, https://doi.org/10.3390/rs70302668, 2015.

Hamon, W. R.: Computation of Direct Runoff Amounts from Storm Rainfall, Int. Assoc. Sci. Hydrol. Publ., 63, 52–62, https://cir.nii.ac.jp/crid/1573387450184176896 (last access: 1 November 2022), 1963.

IGN: Registre Parcellaire graphique, data.gouv.fr [data set], https://www.data.gouv.fr/fr/datasets/registre-parcellaire-graphique-rpg-contours-des-parcelles-et-ilots-culturaux-et-leur-groupe-de-cultures-majoritaire/, 2021.

Inglada, J., Vincent, A., Arias, M., Tardy, B., Morin, D., and Rodes, I.: Operational High Resolution Land Cover Map Production at the Country Scale Using Satellite Image Time Series, Remote Sens., 9, 95, https://doi.org/10.3390/rs9010095, 2017.

Inglada, J., Vincent, A., and Thierion, V.: Theia OSO Land Cover Map, Zenodo [data set], https://doi.org/10.5281/zenodo.3613415, 2018.

Jin, X., Jin, Y., Fu, D., and Mao, X.: Modifying the SWAT Model to Simulate Eco-Hydrological Processes in an Arid Grassland Dominated Watershed, Front. Environ. Sci., 10, 939321, https://www.frontiersin.org/articles/10.3389/fenvs.2022.939321 (last access: 1 April 2021), 2022.

Kharrou, M. H., Simonneaux, V., Er-Raki, S., Le Page, M., Khabba, S., and Chehbouni, A.: Assessing Irrigation Water Use with Remote Sensing-Based Soil Water Balance at an Irrigation Scheme Level in a Semi-Arid Region of Morocco, Remote Sens., 13, 1133, https://doi.org/10.3390/rs13061133, 2021.

Kumar, S. V., Mocko, D. M., Wang, S., Peters-Lidard, C. D., and Borak, J.: Assimilation of Remotely Sensed Leaf Area Index into the Noah-MP Land Surface Model: Impacts on Water and Carbon Fluxes and States over the Continental United States, J. Hydrometeorol., 20, 1359–1377, https://doi.org/10.1175/JHM-D-18-0237.1, 2019.

Lanini, S., Caballero, Y., and Le Cointe, P.: ESPERE User Guide Version 2, Tech. rep., BRGM, https://www.brgm.fr/sites/default/files/documents/2020-11/logiciel-espere-user-guide-v2-en.pdf (last access: 1 February 2021), 2020.
Brochet, E.: Modified SWAT+ NDVI, GitHub [code], https://github.com/ElisabethJustin/SWATplus-NDVI, 2022.

Bruand, A., Fernández, P. P., and Duval, O.: Use of class pedo-transfer functions based on texture and bulk density of clods to generate water retention curves, Soil Use Manage., 19, 232–242, https://doi.org/10.1111/j.1475-2743.2003.tb00309.x, 2003.

Cakir, R., Raimonet, M., Sauvage, S., Paredes-Arquiola, J., Grusson, Y., Roset, L., Meaurio, M., Navarro, E., Sevilla-Callejo, M., Lechuga-Crespo, J. L., Gomiz Pascual, J. J., Bodoque, J. M., and Sánchez-Pérez, J. M.: Hydrological Alteration Index as an Indicator of the Calibration Complexity of Water Quantity and Quality Modeling in the Context of Global Change, Water, 12, 115, https://doi.org/10.3390/w12010115, 2020.

Carroget, A., Perrin, C., Sauquet, E., Vidal, J.-P., Chazot, S., Chauveau, M., and Rouchy, N.: Explore 2070: quelle utilisation d'un exercice prospectif sur les impacts des changements climatiques à l'échelle nationale pour définir des stratégies d'adaptation, Sciences Eaux & Territoires, p. 4, https://doi.org/10.14758/SET-REVUE.2017.22.02, 2017.

Clavel, L., Soudais, J., Baudet, D., and Leenhardt, D.: Integrating expert knowledge and quantitative information for mapping cropping systems, Land Use Policy, 28, 57–65, https://doi.org/10.1016/j.landusepol.2010.05.001, 2011.

Courault, D., Hadria, R., Ruget, F., Olioso, A., Duchemin, B., Hagolle, O., and Dedieu, G.: Combined use of FORMOSAT-2 images with a crop model for biomass and water monitoring of permanent grassland in Mediterranean region, Hydrol. Earth Syst. Sci., 14, 1731–1744, https://doi.org/10.5194/hess-14-1731-2010, 2010.

Demarez, V., Helen, F., Marais-Sicre, C., and Baup, F.: In-Season Mapping of Irrigated Crops Using Landsat 8 and Sentinel-1 Time Series, Remote Sens., 11, 118, https://doi.org/10.3390/rs11020118, 2019.

Drusch, M., Del Bello, U., Carlier, S., Colin, O., Fernandez, V., Gascon, F., Hoersch, B., Isola, C., Laberinti, P., Martimort, P., Meygret, A., Spoto, F., Sy, O., Marchese, F., and Bargellini, P.: Sentinel-2: ESA's Optical High-Resolution Mission for GMES Operational Services, Remote Sens. Environ., 120, 25–36, https://doi.org/10.1016/j.rse.2011.11.026, 2012.

Durand, Y., Brun, E., Merindol, L., Guyomarc'h, G., Lesafre, B., and Martin, E.: A meteorological estimation of relevant parameters for snow models, Ann. Glaciol., 18, 65–71, https://doi.org/10.3189/S0260305500011277, 1993.

ESA: SENTINEL-2, ESA's Optical High-Resolution, Theia [data set], https://theia.cnes.fr/atdistrib/rocket/#/search?collection=SENTINEL2, 2012.

Etchanchu, J., Rivalland, V., Gascoin, S., Cros, J., Tallec, T., Brut, A., and Boulet, G.: Effects of high spatial and temporal resolution Earth observations on simulated hydrometeorological variables in a cropland (southwestern France), Hydrol. Earth Syst. Sci., 21, 5693–5708, https://doi.org/10.5194/hess-21-5693-2017, 2017.

Ferrant, S., Gascoin, S., Veloso, A., Salmon-Monviola, J., Claverie, M., Rivalland, V., Dedieu, G., Demarez, V., Ceschia, E., Probst, J.-L., Durand, P., and Bustillo, V.: Agro-hydrology and multi-temporal high-resolution remote sensing: toward an explicit spatial processes calibration, Hydrol. Earth Syst. Sci., 18, 5219–5237, https://doi.org/10.5194/hess-18-5219-2014, 2014.

Fohrer, N., Dietrich, A., Kolychalow, O., and Ulrich, U.: Assessment of the Environmental Fate of the Herbicides Flufenacet and

Metazachlor with the SWAT Model, J. Environ. Qual., 43, 75–85, https://doi.org/10.2134/jeq2011.0382, 2014.

GisSol: Référentiel Régional Pédologique, Gis-Sol [data set], https://www.gissol.fr/donnees/liens-vers-les-referentiels-regionaux-pedologiques-5634, 2014.

Grusson, Y., Sun, X., Gascoin, S., Sauvage, S., Raghavan, S., Antil, F., and Sáchez-Pérez, J.-M.: Assessing the capability of the SWAT model to simulate snow, snow melt and streamflow dynamics over an alpine watershed, J. Hydrol., 531, 574–588, https://doi.org/10.1016/j.jhydrol.2015.10.070, 2015.

Gupta, H. V., Kling, H., Yilmaz, K. K., and Martinez, G. F.: Decomposition of the mean squared error and NSE performance criteria: Implications for improving hydrological modelling, J. Hydrol., 377, 80–91, https://doi.org/10.1016/j.jhydrol.2009.08.003, 2009.

Guérif, M. and Duke, C.: Calibration of the SUCROS emergence and early growth module for sugar beet using optical remote sensing data assimilation, Eur. J. Agron., 9, 127–136, https://doi.org/10.1016/S1161-0301(98)00031-8, 1998.

Hagolle, O., Huc, M., Villa Pascual, D., and Dedieu, G.: A Multi-Temporal and Multi-Spectral Method to Estimate Aerosol Optical Thickness over Land, for the Atmospheric Correction of FormoSat-2, LandSat, VENUS and Sentinel-2 Images, Remote Sens., 7, 2668–2691, https://doi.org/10.3390/rs70302668, 2015.

Hamon, W. R.: Computation of Direct Runoff Amounts from Storm Rainfall, Int. Assoc. Sci. Hydrol. Publ., 63, 52–62, https://cir.nii.ac.jp/crid/1573387450184176896 (last access: 1 November 2022), 1963.

IGN: Registre Parcellaire graphique, data.gouv.fr [data set], https://www.data.gouv.fr/fr/datasets/registre-parcellaire-graphique-rpg-contours-des-parcelles-et-ilots-culturaux-et-leur-groupe-de-cultures-majoritaire/, 2021.

Inglada, J., Vincent, A., Arias, M., Tardy, B., Morin, D., and Rodes, I.: Operational High Resolution Land Cover Map Production at the Country Scale Using Satellite Image Time Series, Remote Sens., 9, 95, https://doi.org/10.3390/rs9010095, 2017.

Inglada, J., Vincent, A., and Thierion, V.: Theia OSO Land Cover Map, Zenodo [data set], https://doi.org/10.5281/zenodo.3613415, 2018.

Jin, X., Jin, Y., Fu, D., and Mao, X.: Modifying the SWAT Model to Simulate Eco-Hydrological Processes in an Arid Grassland Dominated Watershed, Front. Environ. Sci., 10, 939321, https://www.frontiersin.org/articles/10.3389/fenvs.2022.939321 (last access: 1 April 2021), 2022.

Kharrou, M. H., Simonneaux, V., Er-Raki, S., Le Page, M., Khabba, S., and Chehbouni, A.: Assessing Irrigation Water Use with Remote Sensing-Based Soil Water Balance at an Irrigation Scheme Level in a Semi-Arid Region of Morocco, Remote Sens., 13, 1133, https://doi.org/10.3390/rs13061133, 2021.

Kumar, S. V., Mocko, D. M., Wang, S., Peters-Lidard, C. D., and Borak, J.: Assimilation of Remotely Sensed Leaf Area Index into the Noah-MP Land Surface Model: Impacts on Water and Carbon Fluxes and States over the Continental United States, J. Hydrometeorol., 20, 1359–1377, https://doi.org/10.1175/JHM-D-18-0237.1, 2019.

Lanini, S., Caballero, Y., and Le Cointe, P.: ESPERE User Guide Version 2, Tech. rep., BRGM, https://www.brgm.fr/sites/default/files/documents/2020-11/logiciel-espere-user-guide-v2-en.pdf (last access: 1 February 2021), 2020.
Launay, M. and Guerif, M.: Assimilating remote sensing data into a crop model to improve predictive performance for spatial applications, Agr. Ecosyst. Environ., 111, 321–339, https://doi.org/10.1016/j.agee.2005.06.005, 2005.

Leenhardt, D. and Lemaire, P.: Estimating the spatial and temporal distribution of sowing dates for regional water management, Agr. Water Manage., 55, 37–52, https://doi.org/10.1016/S0378-3774(01)00183-4, 2002.

Leenhardt, D., Trouvat, J. L., Gonzalès, G., Pérarnaud, V., Prats, S., and Bergez, J. E.: Estimating irrigation demand for water management on a regional scale: I. ADEAUMIS, a simulation platform based on bio-decisional modelling and spatial information, Agr. Water Manage., 68, 207–232, https://doi.org/10.1016/j.agwat.2004.04.004, 2004.

Maguire, M. S., Neale, C. M. U., Woldt, W. E., and Heeren, D. M.: Managing spatial irrigation using remote-sensing-based evapotranspiration and soil water adaptive control model, Agr. Water Manage., 272, 107838, https://doi.org/10.1016/j.agwat.2022.107838, 2022.

Martin, E., Gascoin, S., Grusson, Y., Murgue, C., Bardeau, M., Antil, F., Ferrant, S., Lardy, R., Le Moigne, P., Leenhardt, D., Rivalland, V., Sánchez Pérez, J.-M., Sauvage, S., and Therond, O.: On the Use of Hydrological 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, 2016.

Maton, L., Leenhardt, D., Goulard, M., and Bergez, J. E.: Assessing the irrigation strategies over a wide geographical area from structural data about farming systems, Agr. Syst., 86, 293–311, https://doi.org/10.1016/j.agsy.2004.09.010, 2005.

Maton, L., Bergez, J.-E., and Leenhardt, D.: Modelling the days which are agronomically suitable for sowing maize, Eur. J. Agron., 27, 123–129, https://doi.org/10.1016/j.eja.2007.02.007, 2007.

Merlin, O., Bitar, A. A., Rivalland, V., Béziat, P., Ceschia, E., and Dedieu, G.: An Analytical Model of Evaporation Efficiency for Unsaturated Soil Surfaces with an Arbitrary Thickness, J. Appl. Meteorol. Clim., 50, 457–471, https://doi.org/10.1175/2010JAMC2418.1, 2011.

Mohammadi Igder, O., Alizadeh, H., Mojaradi, B., and Bayat, M.: Multivariate assimilation of satellite-based leaf area index and ground-based river streamflow for hydrological modelling of irrigated watersheds using SWAT+, J. Hydrol., 610, 128012, https://doi.org/10.1016/j.jhydrol.2022.128012, 2022.

Monteith, J. L.: Evaporation and environment, Sym. Soc. Exp. Biol., 19, 205–234, https://repository.rothamsted.ac.uk/item/8v5v7/evaporation-and-environment (last access: 1 June 2021), 1965.

Moriasi, D. N., Arnold, J. G., Van Liew, M. W., Bingner, R. L., Harmel, R. D., and Veith, T. L.: Model evaluation guidelines for systematic quantification of accuracy in watershed simulations, T. ASABE, 50, 885–900, 2007.

Moriasi, D. N., Gitau, M. W., Pai, N., and Daggupati, P.: Hydrologic and water quality models: Performance measures and evaluation criteria, T. ASABE, 58, 1763–1785, 2015.

Murgue, C., Lardy, R., Vavasseur, M., Burger-Leenhardt, D., and Therond, O.: Fine spatio-temporal simulation of cropping and farming systems effects on irrigation withdrawal dynamics within a river basin, edited by: Ames, D. P., Quinn, N. W. T., and

Rizzoli, A. E., Proceedings of the 7th International Congress on Environmental Modelling and Software, San Diego, California, USA, 15–19 June, ISBN 978-88-9035-744-2, 2014.

Nash, J. E. and Sutcliffe, J. V.: River flow forecasting through conceptual models, J. Hydrol., 10, 282–290, https://doi.org/10.1016/0022-1694(70)90255-6, 1970.

Nicolle, P., Besson, F., Delaigue, O., Etchevers, P., François, D., Le Lay, M., Perrin, C., Rousset, F., Thiéry, D., Tilmant, F., Mangand, C., Leurent, T., and Jacob, É.: PREMHYCE: An operational tool for low-flow forecasting, Proc. IAHS, 383, 381–389, https://doi.org/10.5194/piahs-383-381-2020, 2020.

Neitsch, S. L., Arnold, J. G., Kiniry, J. R., and Williams, J. R.: SWAT: Soil and water assessment tool theoretical documentation, Temple, TX: USDA Agricultural Research Service, http://swat.tamu.edu/documentation/ (last access: 1 November 2022), 2001.

Olivera-Guerra, L.-E., Laluet, P., Altés, V., Ollivier, C., Pageot, Y., Paolini, G., Chavanon, E., Rivalland, V., Boulet, G., Villar, J.-M., and Merlin, O.: Modeling actual water use under different irrigation regimes at district scale: Application to the FAO-56 dual crop coefficient method, Agr. Water Manage., 278, 108119, https://doi.org/10.1016/j.agwat.2022.108119, 2023.

Pageot, Y., Baup, F., Inglada, J., Baghdadi, N., and Demarez, V.: Detection of Irrigated and Rainfed Crops in Temperate Areas Using Sentinel-1 and Sentinel-2 Time Series, Remote Sens., 12, 3044, https://doi.org/10.3390/rs12183044, 2020.

Paul, M., Rajib, A., Negahban-Azar, M., Shirmohammadi, A., and Srivastava, P.: Improved agricultural Water management in data-scarce semi-arid watersheds: Value of integrating remotely sensed leaf area index in hydrological modeling, Sci. Total Environ., 791, 148177, https://doi.org/10.1016/j.scitotenv.2021.148177, 2021.

Phillips, E. E.: Heat summation theory as applied to canning crops, The Canner, 27, 13–15, 1950.

Pushpalatha, R., Perrin, C., Moine, N. L., and Andréassian, V.: A review of efficiency criteria suitable for evaluating low-flow simulations, J. Hydrol., 420–421, 171–182, https://doi.org/10.1016/j.jhydrol.2011.11.055, 2012.

Puy, A., Sheikholeslami, R., Gupta, H. V., Hall, J. W., Lankford, B., Lo Piano, S., Meier, J., Pappenberger, F., Porporato, A., Vico, G., and Saltelli, A.: The delusive accuracy of global irrigation water withdrawal estimates, Nat. Commun., 13, 3183, https://doi.org/10.1038/s41467-022-30731-8, 2022.

Rallison, R. E. and Miller, N.: Past, present, and future SCS runoff procedure, Littleton, Colo., Water Resources Publications, 353–364, https://api.semanticscholar.org/CorpusID:133775552 (last access: 1 February 2021), 1982.

Rigou, L.: Typologie des sols agricoles du Gers – Rapport de présentation, Tech. rep., Atelier Sols, Urbanisme et Paysages, Angos, France, https://occitanie.chambre-agriculture.fr/agroenvironnement/agroecologie/guide-des-sols-de-midi-pyrenees/sols-du-gers/ (last access: 1 January 2021), 2016.

Rolle, M., Zribi, M., Tamea, S., and Claps, P.: Estimation of maize sowing dates from Sentinel 1&2 data, over South Piedmont, EGU22–10490, EGU General Assembly Conference Abstracts ADS, Vienna, Austria, 23–27 May 2022, https://doi.org/10.5194/egusphere-egu22-10490, 2022.
Launay, M. and Guerif, M.: Assimilating remote sensing data into a crop model to improve predictive performance for spatial applications, Agr. Ecosyst. Environ., 111, 321–339, https://doi.org/10.1016/j.agee.2005.06.005, 2005.

Leenhardt, D. and Lemaire, P.: Estimating the spatial and temporal distribution of sowing dates for regional water management, Agr. Water Manage., 55, 37–52, https://doi.org/10.1016/S0378-3774(01)00183-4, 2002.

Leenhardt, D., Trouvat, J. L., Gonzalès, G., Pérarnaud, V., Prats, S., and Bergez, J. E.: Estimating irrigation demand for water management on a regional scale: I. ADEAUMIS, a simulation platform based on bio-decisional modelling and spatial information, Agr. Water Manage., 68, 207–232, https://doi.org/10.1016/j.agwat.2004.04.004, 2004.

Maguire, M. S., Neale, C. M. U., Woldt, W. E., and Heeren, D. M.: Managing spatial irrigation using remote-sensing-based evapotranspiration and soil water adaptive control model, Agr. Water Manage., 272, 107838, https://doi.org/10.1016/j.agwat.2022.107838, 2022.

Martin, E., Gascoin, S., Grusson, Y., Murgue, C., Bardeau, M., Antil, F., Ferrant, S., Lardy, R., Le Moigne, P., Leenhardt, D., Rivalland, V., Sánchez Pérez, J.-M., Sauvage, S., and Therond, O.: On the Use of Hydrological 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, 2016.

Maton, L., Leenhardt, D., Goulard, M., and Bergez, J. E.: Assessing the irrigation strategies over a wide geographical area from structural data about farming systems, Agr. Syst., 86, 293–311, https://doi.org/10.1016/j.agsy.2004.09.010, 2005.

Maton, L., Bergez, J.-E., and Leenhardt, D.: Modelling the days which are agronomically suitable for sowing maize, Eur. J. Agron., 27, 123–129, https://doi.org/10.1016/j.eja.2007.02.007, 2007.

Merlin, O., Bitar, A. A., Rivalland, V., Béziat, P., Ceschia, E., and Dedieu, G.: An Analytical Model of Evaporation Efficiency for Unsaturated Soil Surfaces with an Arbitrary Thickness, J. Appl. Meteorol. Clim., 50, 457–471, https://doi.org/10.1175/2010JAMC2418.1, 2011.

Mohammadi Igder, O., Alizadeh, H., Mojaradi, B., and Bayat, M.: Multivariate assimilation of satellite-based leaf area index and ground-based river streamflow for hydrological modelling of irrigated watersheds using SWAT+, J. Hydrol., 610, 128012, https://doi.org/10.1016/j.jhydrol.2022.128012, 2022.

Monteith, J. L.: Evaporation and environment, Sym. Soc. Exp. Biol., 19, 205–234, https://repository.rothamsted.ac.uk/item/8v5v7/evaporation-and-environment (last access: 1 June 2021), 1965.

Moriasi, D. N., Arnold, J. G., Van Liew, M. W., Bingner, R. L., Harmel, R. D., and Veith, T. L.: Model evaluation guidelines for systematic quantification of accuracy in watershed simulations, T. ASABE, 50, 885–900, 2007.

Moriasi, D. N., Gitau, M. W., Pai, N., and Daggupati, P.: Hydrologic and water quality models: Performance measures and evaluation criteria, T. ASABE, 58, 1763–1785, 2015.

Murgue, C., Lardy, R., Vavasseur, M., Burger-Leenhardt, D., and Therond, O.: Fine spatio-temporal simulation of cropping and farming systems effects on irrigation withdrawal dynamics within a river basin, edited by: Ames, D. P., Quinn, N. W. T., and

Rizzoli, A. E., Proceedings of the 7th International Congress on Environmental Modelling and Software, San Diego, California, USA, 15–19 June, ISBN 978-88-9035-744-2, 2014.

Nash, J. E. and Sutcliffe, J. V.: River flow forecasting through conceptual models, J. Hydrol., 10, 282–290, https://doi.org/10.1016/0022-1694(70)90255-6, 1970.

Nicolle, P., Besson, F., Delaigue, O., Etchevers, P., François, D., Le Lay, M., Perrin, C., Rousset, F., Thiéry, D., Tilmant, F., Mangand, C., Leurent, T., and Jacob, É.: PREMHYCE: An operational tool for low-flow forecasting, Proc. IAHS, 383, 381–389, https://doi.org/10.5194/piahs-383-381-2020, 2020.

Neitsch, S. L., Arnold, J. G., Kiniry, J. R., and Williams, J. R.: SWAT: Soil and water assessment tool theoretical documentation, Temple, TX: USDA Agricultural Research Service, http://swat.tamu.edu/documentation/ (last access: 1 November 2022), 2001.

Olivera-Guerra, L.-E., Laluet, P., Altés, V., Ollivier, C., Pageot, Y., Paolini, G., Chavanon, E., Rivalland, V., Boulet, G., Villar, J.-M., and Merlin, O.: Modeling actual water use under different irrigation regimes at district scale: Application to the FAO-56 dual crop coefficient method, Agr. Water Manage., 278, 108119, https://doi.org/10.1016/j.agwat.2022.108119, 2023.

Pageot, Y., Baup, F., Inglada, J., Baghdadi, N., and Demarez, V.: Detection of Irrigated and Rainfed Crops in Temperate Areas Using Sentinel-1 and Sentinel-2 Time Series, Remote Sens., 12, 3044, https://doi.org/10.3390/rs12183044, 2020.

Paul, M., Rajib, A., Negahban-Azar, M., Shirmohammadi, A., and Srivastava, P.: Improved agricultural Water management in data-scarce semi-arid watersheds: Value of integrating remotely sensed leaf area index in hydrological modeling, Sci. Total Environ., 791, 148177, https://doi.org/10.1016/j.scitotenv.2021.148177, 2021.

Phillips, E. E.: Heat summation theory as applied to canning crops, The Canner, 27, 13–15, 1950.

Pushpalatha, R., Perrin, C., Moine, N. L., and Andréassian, V.: A review of efficiency criteria suitable for evaluating low-flow simulations, J. Hydrol., 420–421, 171–182, https://doi.org/10.1016/j.jhydrol.2011.11.055, 2012.

Puy, A., Sheikholeslami, R., Gupta, H. V., Hall, J. W., Lankford, B., Lo Piano, S., Meier, J., Pappenberger, F., Porporato, A., Vico, G., and Saltelli, A.: The delusive accuracy of global irrigation water withdrawal estimates, Nat. Commun., 13, 3183, https://doi.org/10.1038/s41467-022-30731-8, 2022.

Rallison, R. E. and Miller, N.: Past, present, and future SCS runoff procedure, Littleton, Colo., Water Resources Publications, 353–364, https://api.semanticscholar.org/CorpusID:133775552 (last access: 1 February 2021), 1982.

Rigou, L.: Typologie des sols agricoles du Gers – Rapport de présentation, Tech. rep., Atelier Sols, Urbanisme et Paysages, Angos, France, https://occitanie.chambre-agriculture.fr/agroenvironnement/agroecologie/guide-des-sols-de-midi-pyrenees/sols-du-gers/ (last access: 1 January 2021), 2016.

Rolle, M., Zribi, M., Tamea, S., and Claps, P.: Estimation of maize sowing dates from Sentinel 1&2 data, over South Piedmont, EGU22–10490, EGU General Assembly Conference Abstracts ADS, Vienna, Austria, 23–27 May 2022, https://doi.org/10.5194/egusphere-egu22-10490, 2022.
Saadi, S., Simonneaux, V., Boulet, G., Raimbault, B., Mougenot, B., Fanise, P., Ayari, H., and Lili-Chabaane, Z.: Monitoring Irrigation Consumption Using High Resolution NDVI Image Time Series: Calibration and Validation in the Kairouan Plain (Tunisia), Remote Sens., 7, 13005–13028, https://doi.org/10.3390/rs71013005, 2015.

Santos, L., Thirel, G., and Perrin, C.: Technical note: Pitfalls in using log-transformed flows within the KGE criterion, Hydrol. Earth Syst. Sci., 22, 4583–4591, https://doi.org/10.5194/hess-22-4583-2018, 2018.

Schumacher, D. L., Zachariah, M., and Otto, F.: High temperatures exacerbated by climate change made 2022 Northern Hemisphere droughts more likely, https://policycommons.net/artifacts/3174587/wce-nh-drought-scientific-report/3973082/ (last access: 1 June 2022), 2022.

Senthilkumar, K., Bergez, J.-E., and Leenhardt, D.: Can farmers use maize earliness choice and sowing dates to cope with future water scarcity? A modelling approach applied to south-western France, Agr. Water Manage., 152, 125–134, https://doi.org/10.1016/j.agwat.2015.01.004, 2015.

Soil Conservation Service: Section 4: Hydrology – Chapter 15: Travel time, time of concentration and lag, in: National engineering handbook, SCS USDA, https://directives.sc.egov.usda.gov/OpenNonWebContent.aspx?content=27002.wba (last access: 1 May 2021), 1972.

Toureiro, C., Serralheiro, R., Shahidian, S., and Sousa, A.: Irrigation management with remote sensing: Evaluating irrigation requirement for maize under Mediterranean climate condition, Agr. Water Manage., 184, 211–220, https://doi.org/10.1016/j.agwat.2016.02.010, 2017.

UNESCO: The United Nations world water development report 2015: water for a sustainable world, UNESCO Publishing, ISBN 978-92-3-100071-3, 2015.

USGS: Shuttle Radar Topography Mission 1 Arc-Second Global, USGS [data set], https://doi.org/10.5066/F7PR7TFT, 2018.

Yousaf, W., Awan, W. K., Kamran, M., Ahmad, S. R., Bodla, H. U., Riaz, M., Umar, M., and Chohan, K.: A paradigm of GIS and remote sensing for crop water deficit assessment in near real time to improve irrigation distribution plan, Agr. Water Manage., 243, 106443, https://doi.org/10.1016/j.agwat.2020.106443, 2021.
Saadi, S., Simonneaux, V., Boulet, G., Raimbault, B., Mougenot, B., Fanise, P., Ayari, H., and Lili-Chabaane, Z.: Monitoring Irrigation Consumption Using High Resolution NDVI Image Time Series: Calibration and Validation in the Kairouan Plain (Tunisia), Remote Sens., 7, 13005–13028, https://doi.org/10.3390/rs71013005, 2015.

Santos, L., Thirel, G., and Perrin, C.: Technical note: Pitfalls in using log-transformed flows within the KGE criterion, Hydrol. Earth Syst. Sci., 22, 4583–4591, https://doi.org/10.5194/hess-22-4583-2018, 2018.

Schumacher, D. L., Zachariah, M., and Otto, F.: High temperatures exacerbated by climate change made 2022 Northern Hemisphere droughts more likely, https://policycommons.net/artifacts/3174587/wce-nh-drought-scientific-report/3973082/ (last access: 1 June 2022), 2022.

Senthilkumar, K., Bergez, J.-E., and Leenhardt, D.: Can farmers use maize earliness choice and sowing dates to cope with future water scarcity? A modelling approach applied to south-western France, Agr. Water Manage., 152, 125–134, https://doi.org/10.1016/j.agwat.2015.01.004, 2015.

Soil Conservation Service: Section 4: Hydrology – Chapter 15: Travel time, time of concentration and lag, in: National engineering handbook, SCS USDA, https://directives.sc.egov.usda.gov/OpenNonWebContent.aspx?content=27002.wba (last access: 1 May 2021), 1972.

Toureiro, C., Serralheiro, R., Shahidian, S., and Sousa, A.: Irrigation management with remote sensing: Evaluating irrigation requirement for maize under Mediterranean climate condition, Agr. Water Manage., 184, 211–220, https://doi.org/10.1016/j.agwat.2016.02.010, 2017.

UNESCO: The United Nations world water development report 2015: water for a sustainable world, UNESCO Publishing, ISBN 978-92-3-100071-3, 2015.

USGS: Shuttle Radar Topography Mission 1 Arc-Second Global, USGS [data set], https://doi.org/10.5066/F7PR7TFT, 2018.

Yousaf, W., Awan, W. K., Kamran, M., Ahmad, S. R., Bodla, H. U., Riaz, M., Umar, M., and Chohan, K.: A paradigm of GIS and remote sensing for crop water deficit assessment in near real time to improve irrigation distribution plan, Agr. Water Manage., 243, 106443, https://doi.org/10.1016/j.agwat.2020.106443, 2021.