Net irrigation requirement under different climate scenarios using AquaCrop over Europe

. Global soil water availability is challenged by the effects of climate change and a growing population. On average, 70 % of freshwater extraction is attributed to agriculture, and the demand is increasing. In this study, the effects of climate change on the evolution of the irrigation water requirement to sustain current crop productivity are assessed by using the Food and Agriculture Organization (FAO) crop growth model AquaCrop version 6.1. The model is run at 0 . 5 ◦ lat × 0 . 5 ◦ long resolution over the European mainland, assuming a general C3-type of crop, and forced by climate input data from the Inter-Sectoral Impact Model Intercomparison Project phase three (ISIMIP3). First, the AquaCrop surface soil moisture (SSM) forced with two types of ISIMIP3 historical meteorological datasets is evaluated with satellite-based SSM estimates in two ways. When driven by ISIMIP3a reanalysis meteorology, daily simulated SSM values have an unbiased root mean square difference of 0.08 and 0.06 m 3 m − 3 , with SSM retrievals from the Soil Moisture Ocean Salinity (SMOS) and Soil Active Passive (SMAP)

Abstract. Global soil water availability is challenged by the effects of climate change and a growing population. On average, 70 % of freshwater extraction is attributed to agriculture, and the demand is increasing. In this study, the effects of climate change on the evolution of the irrigation water requirement to sustain current crop productivity are assessed by using the Food and Agriculture Organization (FAO) crop growth model AquaCrop version 6.1. The model is run at 0.5 • lat × 0.5 • long resolution over the European mainland, assuming a general C3-type of crop, and forced by climate input data from the Inter-Sectoral Impact Model Intercomparison Project phase three (ISIMIP3).
First, the AquaCrop surface soil moisture (SSM) forced with two types of ISIMIP3 historical meteorological datasets is evaluated with satellite-based SSM estimates in two ways. When driven by ISIMIP3a reanalysis meteorology, daily simulated SSM values have an unbiased root mean square difference of 0.08 and 0.06 m 3 m −3 , with SSM retrievals from the Soil Moisture Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) missions, respectively, for the years 2015-2016 (2016 is the end year of the reanalysis data). When forced with ISIMIP3b meteorology from five global climate models (GCMs) for the years 2015-2020, the historical simulated SSM climatology closely agrees with the satellite-based SSM climatologies.
Second, the evaluated AquaCrop model is run to quantify the future irrigation requirement, for an ensemble of five GCMs and three different emission scenarios. The simulated net irrigation requirement (I net ) of the three summer months for a near and far future climate period (2031-2060 and 2071-2100) is compared to the baseline period of 1985-2014 to assess changes in the mean and interannual variabil-ity of the irrigation demand. Averaged over the continent and the model ensemble, the far future I net is expected to increase by 22 mm per month (+30 %) under a high-emission scenario Shared Socioeconomic Pathway (SSP) 3-7.0. Central and southern Europe are the most impacted, with larger I net increases. The interannual variability in I net is likely to increase in northern and central Europe, whereas the variability is expected to decrease in southern regions. Under a high mitigation scenario (SSP1-2.6), the increase in I net will stabilize at around 13 mm per month towards the end of the century, and interannual variability will still increase but to a smaller extent. The results emphasize a large uncertainty in the I net projected by various GCMs.

Introduction
Global crop production has vastly increased over the past century, leading to the expansion of irrigated areas by almost 6 times and more pressure on the irrigation water demand (Siebert et al., 2015). With changing climatic conditions and a growing population, future water availability is expected to further decline, raising demands for more efficient irrigation systems Taylor et al., 2013) and a higher crop water productivity (Brauman et al., 2013). In this context, a range of modeling studies have tried to assess future impacts on agricultural water demands and possible actions, but this remains a difficult task due to high uncertainties in future climate and socioeconomic scenarios Haddeland et al., 2014;Wada et al., 2013).
Future meteorological variables are typically modeled by global climate models (GCMs) for different scenarios, usu-ally represented by the Representative Concentration Pathways (RCPs;van Vuuren et al., 2011). Some challenges associated with climate forcing data are the consistency and the representation of the uncertainty of the data. The Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) is an initiative to provide consistent, bias-corrected climate datasets for impact modeling Warszawski et al., 2014). The project is currently at its third simulation round (ISIMIP3) and provides reanalysis historical climate (ISIMIP3a) and GCM-driven historical and future climate (ISIMIP3b), following different emission scenarios and using various GCMs. Data from the previous simulation round (ISIMIP2) have already been used in several studies of historical and future water resources (e.g., Lange et al., 2020;Reinecke et al., 2021).
Based on such climate projections, it is possible to derive meteorological drought indicators, which are determined by precipitation (P ) and the atmospheric evaporation demand (ET 0 ). These meteorological droughts propagate into agricultural and hydrological droughts, characterized by a reduction in the soil water content and a reduction in streamflow. Over this century, droughts are expected to become more frequent in the Northern Hemisphere (Sheffield and Wood, 2008), in most parts of Europe (Spinoni et al., 2018;Grillakis, 2019), and especially in southern Europe Russo et al., 2013;Ruosteenoja et al., 2018). Common meteorological drought indices are directly associated to variations in P and ET 0 (Vicente-Serrano et al., 2015). The difference between these two fluxes (P -ET 0 ), also referred to as the climatic water balance, has served as proxy to investigate drying trends (Greve et al., 2014;Prȃvȃlie et al., 2019). For agriculture, P -ET 0 is also a major factor determining the need for additional water, i.e., for irrigation.
In Europe, rainfall fulfills the largest part of the crop water requirement (green water), but irrigation (blue water) becomes essential in most southern parts of the continent Liu and Yang, 2010;Siebert and Döll, 2010). For the past few decades, the yearly net irrigation requirement in Europe has been estimated between 53 to 1120 mm yr −1 in Denmark and Spain, respectively (Wriedt et al., 2009). The effectively applied amounts of irrigation could be much lower or higher but are unknown due to the lack of good observational data (Massari et al., 2021). Future global and regional irrigation trend assessments have commonly used hydrological models (e.g., WaterGAP; Döll and Siebert, 2002, in Döll, 2002, and Schaldach et al., 2012, agro-ecosystem models (Lund-Potsdam-Jena managed Land -LPJmL; Bondeau et al., 2007, in Fader et al., 2016, and Konzmann et al., 2012, or agro-ecological zone (AEZ) models (FAO-AEZ methodology applied in Fischer et al., 2007). The earliest global study addressing future irrigation requirement under climate change was performed by Döll (2002), using the WaterGAP model for two GCMs. The results indicate clear effects on the long-term average irrigation re-quirement, with an average global increase of ∼ 10 % by the 2070s under the Intergovernmental Panel on Climate Change (IPCC) IS92a scenario (Leggett et al., 1992). Similar increases were found later by Fischer et al. (2007), who also used two GCMs applied to an emission scenario from the IPCC Special Report on Emission Scenarios (SRES A2r; Nakicenovic et al., 2000;Riahi et al., 2007). In contrast, global decreases in irrigation water demand were simulated by Pfister et al. (2011) and Konzmann et al. (2012) for the end of the century. These studies only assessed one emission scenario, both from the IPCC SRES (Nakicenovic et al., 2000), namely the A1B and A2 scenarios, respectively. However, in Europe, all of these studies indicate clear increases in irrigation water requirement for most parts of the continent where irrigation is currently applied.
The outcomes of the different irrigation assessments can diverge quite significantly. Wada et al. (2013) provided an ensemble of seven general hydrological models (GHMs; including LPJmL and WaterGAP) and analyzed the sources of uncertainty on the final predictions. The results showed that the fraction of the variance due to the GCMs is larger than the fraction caused by the future emission scenarios, and that more than 50 % of the variance resulted from the GHMs. The experimental setup also plays a major role, as many parameters can influence the irrigation requirement. Global figures are highly different, depending on whether the expansion of irrigated areas is considered or not, which explains why Fischer et al. (2007) found increases in the average water requirement, whereas Pfister et al. (2011) and Konzmann et al. (2012) expected a global decrease. The conclusions also depend on (i) whether irrigation efficiencies are considered (i.e., including socioeconomic factors), (ii) the delineation of the growing season (a whole year; a fixed or flexible start), and (iii) the type of implementation of irrigation in the model (gross or net requirement, threshold to trigger irrigation, and amount of water applied; Telteu et al., 2021).
The irrigation requirement can also be estimated with crop models, which have the added benefit of estimating future trends in crop production and thereby provide useful information to farmers and decision-makers in their adaptation management strategies under climate change. Crop models mainly aim to present quantitative knowledge about the crop development and crop yield for a given crop with specific features and subject to given environmental conditions (Monteith, 1996). Crop modeling integrates physiological processes and the interactions between the crop and its environment. Several studies have shown the added value of upscaling field-scale crop models to a regional level (e.g., Balkovič et al., 2013;Boogaard et al., 2013;de Wit and van Diepen, 2007;Stöckle et al., 2014), allowing current and future crop yield and irrigation assessments. Liu and Yang (2010) used a geographic information system (GIS)-based version of the EPIC (Williams et al., 1989) crop model to spatially evaluate the crop consumptive water use, partitioning the precipitation input, and the irrigation requirement for the year 2000. Pfister et al. (2011) used CROPWAT (Smith, 1992 to compute the global increase in irrigation requirements to meet future food and biomass demands.  provided estimations of the potential irrigation water consumption with 10 GHMs (similar to Wada et al., 2013) and six global gridded crop models (GGCMs; developed within the Agricultural Model Intercomparison and Improvement Project (AgMIP) framework, Rosenzweig et al., 2014), of which three are upscaled site-based crop models. Global-scale crop modeling remains challenging, especially at coarser resolutions (e.g., 0.5 • × 0.5 • lat-long), where one grid cell may contain the information of many heterogeneous agricultural fields . In addition, field management practices (e.g., irrigation practices and fertilizer application) are even more challenging to integrate at regional and global levels. In this study, the spatial version of AquaCrop developed by de Roos et al. (2021b) will be used. AquaCrop Raes et al., 2009), set up as a field-scale model, was developed by the Food and Agriculture Organization of the United Nations (FAO) and is based on the soil water balance. Compared to other, more complex canopy-level models, AquaCrop stands out for its relatively few and intuitive input parameters ). AquaCrop has already been used in regional agricultural climate impact studies by Dale et al. (2017), where an open-source version of AquaCrop (AquaCrop-OS; Foster et al., 2017) was used to project crop yields for a high number of GCMs under different climate scenarios at a resolution of 2 • × 2 • . In this study, the impact of climate change on the future net irrigation requirement is assessed for different emission scenarios and GCMs, using the spatial version of AquaCrop (de Roos et al., 2021b) forced with ISIMIP3 meteorological data over the European continent for the first time. First, the model performance is evaluated by comparing historical spatial AquaCrop v6.1 simulations without any irrigation, forced with (i) reanalysis data from ISIMIP3a and (ii) GCM-based meteorological data from ISIMIP3b, against satellite-based surface soil moisture (SSM) observations. Next, AquaCrop v6.1 simulations are performed using an ensemble of five ISIMIP3b GCMs as forcing to provide estimates of changes in the net irrigation water requirement (I net ) during the summer months (June, July, and August) for two periods in the future (2031-2060 and 2071-2100). The focus is mainly on estimating water demand during the summer period and not on crop water productivity. The objective is to regionally quantify the mean and interannual variability in summer I net for a near and far future climate period and relate this to the current (baseline) I net and future changes in P -ET 0 following various climate scenarios. Compared to previous studies, the advantages are that the simulations are performed with (i) climate data from the latest generation of reanalyses and GCMs, (ii) the most recent set of future scenarios, and (iii) a crop model (AquaCrop), in which the dynamic interactions between water and vegetation are the main fo-cus and where irrigation and management practices can be included with more detail than in a land surface or hydrological model. Future I net projections could be used to inform on climate change adaptation strategies (e.g., climatesmart irrigation, crop type selection, and water conservation). The new AquaCrop-ISIMIP3 model setup can be run at any spatial domain and resolution, providing future opportunities for further climate analysis that also include other irrigation practices and management options.
2 Model and data

Model setup
The study domain focuses on the part of the European continent with latitudes (lat) ranging from 34.75 to 59.75 • N and longitudes (long) from −10.75 to 41.25 • E. The spatial and temporal resolutions of the model simulations are set to those of the ISIMIP3 input datasets, i.e., 0.5 • lat × 0.5 • long, and daily time steps. The same spatial AquaCrop (v6.1) model structure, as described by de Roos et al. (2021b), is used for this study, but adaptations are made to the spatial resolution, input datasets, and simulated periods. Simulations are performed from 1985 through 2100, either with or without considering irrigation, and with the respective associated croprelated parameters.

Model parameters
Climate impact assessments are subject to large uncertainties, which increase with longer temporal projections. Therefore, several assumptions are made in this study to limit the uncertainty of factors other than climate. We will present net irrigation requirement values that are independent of the irrigated area, period, infrastructure, and the exact crop type. First, simulations are performed over all pixels of the entire study domain (i.e., the main European continent), and the irrigation estimates for the entire hypothetically irrigated agricultural domain are normalized by area to make the results independent of the actual irrigated area. This avoids the need to include estimates of future hypothetical land use (Prestele et al., 2016) and the uncertain evolution of the extent of irrigated areas (Schaldach et al., 2012;Hurtt et al., 2020). Second, the spatial resolution of this study matches that of the ISIMIP input data resolution. In contrast to finescale agricultural studies, which usually assess actual irrigation under historical conditions, future climate projections are dependent on the resolution of the driving climate models (or downscaled output). Such studies mainly aim at estimating the irrigation requirement that is needed for crop root uptake, thereby omitting the part of irrigation that is lost to the atmosphere or retained on the soil surface or in the soil profile. Also, state-of-the-art global and continentalscale climate impact assessments are typically performed at the same resolution (e.g., Jägermeyr et al., 2021;Lange et al., 2020;Thiery et al., 2021). Third, each pixel is defined as a hypothetical homogeneous field in which the vegetation conditions are identical. For future projections, the use of a representative field crop is supported by the current lack of detailed year-and location-specific crop maps, and by the unpredictability of changes and developments in crop type and distribution. Finally, the uncertainty and high spatial and temporal variability in the start and end of the growing season (King et al., 2018;Menzel and Fabian, 1999;Schaldach et al., 2012) restricts the modeling possibilities. Some previous studies (e.g., Elliott et al., 2014;Fader et al., 2016;Fischer et al., 2007;Konzmann et al., 2012) have used dynamic growing seasons, but the choice has been made to avoid this additional level of uncertainty for this study. Therefore, only the summer months are considered to make the future requirement directly comparable to the baseline I net . On average, these are the months presenting the highest I net (Siebert and Döll, 2010) and are expected to remain important months for irrigation requirements, even if growing seasons might shift in the future.
Soil data are extracted from the ISIMIP3 soil input dataset that has been used in the AgMIP GGCM intercomparison (GGCMI; Rosenzweig et al., 2014). ISIMIP3 uses the Harmonized World Soil Database version 1.2 (HWSD1.2) aggregated to 0.5 • resolution. The soil dataset represents dominant soil types on croplands within each pixel. There are two soil layers implemented in AquaCrop, i.e., one topsoil layer of 0.30 m and an underlying layer of 1 m, both with the same ISIMIP (topsoil only) textural properties (clay, sand, and silt fractions) and gravel content but with different derived soil hydraulic parameters. More specifically, the volumetric soil water content at saturation, field capacity, and permanent wilting point (θ s , θ FC , and θ PWP ) and the saturated hydraulic conductivity (K sat ) are derived using depth-specific (topsoil and subsoil) pedotransfer functions described by De Lannoy et al. (2014). Because the crop rooting depth is set to 1 m and various bedrock maps indicate that the soil depth over Europe reaches below 1 m (Dirmeyer and Oki, 2002;Mahanama et al., 2015;Shangguan et al., 2017), no limitations to root development need to be considered . A total profile depth of 1.30 m is defined, but without the presence of a groundwater table or confining layers, the actual depth below the maximum rooting depth has no influence on the simulations. For the historical model evaluation with satellite retrievals, the choice is made to use a general C3-type of crop with a 1 m rooting depth to describe the vegetation component, similar to de Roos et al. (2021b). This choice is motivated by the coarse spatial resolution and the high uncertainty in crop modifications over time and follows the methodology of well-known hydrological and land surface models that also make use of general vegetation descriptions (e.g., Niu et al., 2011;Rodell et al., 2004). C3 crops are dominant in Europe (Monfreda et al., 2008;Still et al., 2003). A detailed description of the crop characteristics is given in Table 1 of de Roos et al. (2021b). For the model evaluation, no irrigation is activated, the soil fertility stress of 30 % is maintained (de Roos et al., 2021b), and the AquaCrop default record of mean annual CO 2 concentration observed at Mauna Loa (Hawaii, USA) is considered in the simulations. In contrast, the simulations with irrigation follow the yearly CO 2 concentrations of the emission scenarios from ISIMIP3 and assume near-optimal soil fertility.
For the determination of I net , a representative field crop is considered. The crop characteristics that determine crop transpiration and hence I net are listed in Table 1. The considered crop transpiration coefficient of 1.10 is a good indicative value of the basal crop coefficient for the mid-season for a large range of field crops (Allen et al., 1998). Moreover, it is assumed that, in the summer months (in which I net is determined), the crop has reached its maximum canopy cover and is prior to senescence. Since I net is determined by keeping the soil water content in the root zone above 50 % of the readily available water (RAW, which is 25 % of the total available water, TAW, for the representative field crop), water stress does not affect crop transpiration. Also, air temperature stress affecting crop transpiration will be small or absent in the summer months with the settings of the thresholds in Table 1. To be sure of a well-developed crop canopy during the 3 summer months, it is assumed, in the simulations with irrigation, that the crop germinated in early spring and that the natural crop senescence occurred in late autumn. Irrigated fields are assumed to be well managed. Hence, a nearoptimal soil fertility is defined in AquaCrop, corresponding to a potential achievable biomass production (without any other stress) of 80 % (compared to 70 % for the simulations without irrigation). Future elevated CO 2 concentrations are expected to increase biomass production by reducing crop transpiration and stimulating crop production (CO 2 fertilization effect; Vanuytrecht et al., 2012). This response can vary according to intrinsic crop characteristics or nutrient availability (Vanuytrecht et al., 2011). To avoid overexpression of this effect, the sink term in AquaCrop is lowered to 0 %.

Meteorological data
The AquaCrop model is run with both reanalysis (IS-MIP3a) and GCM-based (ISIMIP3b) meteorological input. The ISIMIP3a forcing data extend up to the end of 2016 and are based on the bias-corrected ECMWF Reanalysis data fifth generation (ERA5; Cucchi et al., 2020;Lange, 2019a). The GCM (ISIMIP3b) data start in 2015 and are derived from the following five different GCMs contributing to the Coupled Model Intercomparison Project phase 6 (CMIP6): GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, and UKESM1-0-LL (Lange, 2019b(Lange, , 2020. These future climate data are separated into different scenarios, which are based on the new scenario framework  , and y refers to the level of radiative forcing (in W m −2 ) in 2100 (RCP). There are three scenarios evaluated, i.e., SSP1-2.6 (low emissions thanks to strong mitigation), SSP3-7.0 (high emissions), and SSP5-8.5 (extreme emissions or unmitigated), for five GCMs, resulting in a total of 15 SSP-GCM scenarios. Under SSP3-7.0 and SSP5-8.5, a global warming of 2 • C will likely be exceeded by mid-century. AquaCrop requires minimum and maximum temperature, rainfall, and reference evapotranspiration (ET 0 ) on a daily basis. Meteorological variables extracted from ISIMIP3 are the daily maximum and minimum temperatures, total precipitation, near-surface relative humidity, near-surface wind speed (at a 10 m height), and the shortwave downwelling radiation. Daily ET 0 values are estimated with the FAO Penman-Monteith equation, according to the guidelines presented in the FAO "Irrigation and Drainage Paper No. 56" (Allen et al., 1998), with the available variables and ISIMIP elevation data (for the estimation of the atmospheric pressure).

Satellite-based evaluation data
To evaluate the performance of the regional AquaCrop simulations forced with ISIMIP input, the following two Lband microwave-based level 2 SSM products are used: (i) the SMUDP2 data product version 650 from the ESA Soil Moisture Ocean Salinity (SMOS) mission (Kerr et al., 2010), from 2011 onwards, and (ii) the SPL2SMP product version 7 from the NASA Soil Moisture Active Passive (SMAP) mission (Chan et al., 2016), from 2015 onwards. For both data sources, only the recommended quality retrievals are included. Additionally, retrievals for daily minimum temperatures below 4 • C are screened out to avoid retrievals of nearfrozen conditions. Both satellite products are projected on a 36 km Equal-Area Scalable Earth version 2 (EASEv2) grid, for SMOS data after reprojection, as in De Lannoy and Reichle (2016). It should be noted that SMOS data over Europe have been affected by radio frequency interference and are filtered out, especially in the early years after launch in 2010 (Oliva et al., 2012).

Simulations
There are three types of simulation experiments performed, and these are referred to as SIM1, SIM2, and SIM3, with the corresponding settings described in Table 2. SIM1 and SIM2 constitute the historical model evaluation against satellite retrieval products. For SIM1, reanalysis meteorological data (ISIMIP3a) are used as input, and simulated SSM is compared to satellite reference data at a daily resolution (shortterm variability). AquaCrop is run over the study area for the period from 1 January 2011 through 31 December 2016 with reanalysis data, i.e., until the end of the available reanalysis data. The second set of historical simulations (SIM2) is GCM-driven (ISIMIP3b) SSM simulations. The purpose of SIM2 is to determine whether the GCM-based forcing is reliable to use for future simulations, i.e., via an evaluation of multi-year average SSM (long-term distribution). For each GCM, AquaCrop is run with climate input data for the period 2011-2020. These input data gather historical simulated climate for 2011-2014 and scenario-based simulated climate for the period 2015-2020, only accounting for SSP5-8.5 (only small differences occur between the three SSPs for this time period). The meteorological time series of the two periods are stitched together to provide continuous AquaCrop forcing fields for 2011-2020. SIM1 and SIM2 have a spinup period of 4 years, and only output from 2015 onwards is used for evaluation, i.e., starting when both SMOS and SMAP data are available.
Once the model has been evaluated with the first two experiments, simulations of SIM3 are run with GCM-driven meteorological input (ISIMIP3b) for the baseline (historical reference period, 1985-2014) and into the future from 2021 through 2100. Irrigation is activated in AquaCrop, and the net irrigation water requirement I net for the 3 summer months is extracted from the simulations for the reference time window and two future time horizons (near future 2031-2060; far future 2071-2100). For the baseline simulation, the initial moisture conditions are set to field capacity, while the future periods have a spin up of at least 10 years (continuous simulation from 2021 through 2100). For SIM3, irrigation is introduced, using the net irrigation requirement option in AquaCrop, whereby a small amount of water (just covering the crop ET for that day) is injected into the root system on days when a certain fraction of the RAW is depleted . With this option, only the amount of water taken up by the roots is considered, where the wetting of the soil surface and interval and application amount specific to a particular irrigation method are not relevant. By selecting a threshold of 50 % RAW depletion, which is the average depletion in an optimal irrigation interval (Smith, 1992), crop water stress affecting the canopy development and transpiration of the representative field crop is avoided, and effective rainfall (the part stored in the root system up to field capacity) is still considered. All simulations performed in this research are uncoupled, i.e., feedback mechanisms from irrigation on atmospheric climate (e.g., Hirsch et al., 2017;Thiery et al., 2017 are neglected.

Skill metrics
To compare the spatial and temporal patterns of SSM from 0.5 • AquaCrop simulations with 36 km satellite data, nearest-neighbor sampling is used to spatially match simulated SSM with SMOS and SMAP retrievals. The output variable extracted from the AquaCrop simulation is the volumetric water content of the topsoil compartment, corresponding to the first 0.1 m of the soil (output variable WC01 in AquaCrop). After quality screening of the satellite data (see Sect. 2.4), about 0.9 ×10 6 and 1.9 ×10 6 usable observations are kept over the study domain (composed of 3882 pixels) for SMOS and SMAP, within the period April 2015-December 2020. The most widely used validation metrics for SSM estimates from large-scale model simulations and retrievals are the Pearson correlation coefficient (R), the bias, the root mean square difference (RMSD), and the unbiased RMSD (ubRMSD), which are calculated as follows: where x is the simulated SSM, y the reference observations, N the number of observation-simulation pairs, and (.) is the temporal mean. A minimum threshold of N = 100 reference data points in time are set per pixel for all analyses. Anomaly correlations are discussed in Appendix A. The aim of the historical evaluation is to assess the performance of AquaCrop to integrate ISIMIP3 meteorological forcings and to provide SSM estimates. By design, the model and satellite retrievals are biased due to model parameters related to the soil and the uniform vegetation type (generic C3 crop), vertical representativeness bias, etc. Therefore, bias-free metrics (R and ubRMSD) are essential to assess whether the main temporal variations in SSM are captured by the model forced with ISIMIP3 data.

Difference in evaluation for SIM1 and SIM2
Both the time series of historical SIM1 and SIM2 SSM are compared to satellite observations through the skill metrics described in Sect. 3.2.1 for the time period with available data for both SMOS and SMAP, i.e., from April 2015 through December 2016. SIM1 is a short-term evaluation, since daily SSM simulations are compared to satellite observations. All months of the year with available and qualitative satellite data were included in this first validation step. For the historical SIM2 SSM simulations, the multi-year average (long-term) results driven by the five different GCMs and the median SSM time series across the GCMs are evaluated. However, for each simulation year, only the period between the 1 March up to the 31 October is considered because only summer months will be considered for the sub-sequent analysis of future I net (Sect. 3.3). Climate models are developed to indicate changing climatic trends but do not present daily accurate data if they are not constrained by observational data. Therefore, the multi-year average (i.e., climatology) of SIM2 SSM is computed and then compared to the climatologies of satellite SSM during the observation period, using the same skill metrics presented in Sect. 3.2.1. Climatologies are calculated using a sliding window of 31 d with a minimum threshold of three data points of data within the window. The computation of the climatology is restrained to the availability of reference satellite data (i.e., SMAP data available from April 2015), as it is also the case in satellite data assimilation systems (e.g., SMAP Level 4 product; Reichle et al., 2019).

Future net irrigation I net requirement (SIM3)
This study focuses on the evaluation of the change in I net during the period for which the highest irrigation demand is expected in all parts of Europe, i.e., June, July, and August (Siebert and Döll, 2010). For the evaluation of the future irrigation water requirement, daily I net values (directly available from the model output) are first extracted from the SIM3 output of the 15 different SSP-GCM combinations. The results are expressed in millimeters per month by averaging the I net of the 3 summer months. The summer irrigation is then used for evaluation, following two approaches. First, the summer I net is averaged over the 30-year time window, allowing us to compare future (2031-2060 and 2071-2100) and baseline  average I net by computing the difference ( I net ). A statistical t test is carried out to define whether the difference of mean I net between the two periods is significant (p < 0.05). Second, interannual variation is assessed based on the I net range (RI net ), defined as the difference between maximum and minimum summer I net of the 30-year time window. Again the difference between future and baseline RI net is evaluated ( RI net ). I net values simulated for the different SSP-GCM combinations are analyzed individually. Additionally, the median results across the GCMs for each scenario are presented. A simple climate index (P -ET 0 ), computed for the 3 summer months, is used to identify where drying trends are potentially occurring, and how this is reflected in the irrigation requirement.  Table 3. The skill is generally better relative to SMAP SSM than relative to SMOS SSM. The expected errors of both missions are 0.04 m 3 m −3 when comparing the satellite data to in situ reference data (Entekhabi et al., 2014). Here, slightly higher ubRMSDs of 0.06 and 0.08 m 3 m −3 are obtained.
The spatial distribution of ubRMSD is presented in Fig. 1a. For further discussion, a partitioning of the study domain in various zones is shown in Appendix A (Fig. A1). Simulated SSM deviate more from SMOS retrievals in north and central-eastern Europe, whereas pixels located in southern regions (e.g., Spain) present a better model performance when comparing to SMOS. Central-eastern Europe presents on average a higher ubRMSD, stressing a lower performance in this region. Time series of SSM estimates at two locations are shown in Fig. 1b. The modeled SSM are close to satellite retrievals for the first pixel (left), and a mismatch is found between simulations and retrievals for the second pixel (right). For the latter, AquaCrop simulations are underestimating SSM during summer, and it can be noticed that SMOS and SMAP retrievals substantially diverge for this location.

Climatological evaluation (SIM2)
The SIM2 AquaCrop SSM for the period 2011-2020 is forced with ISIMIP3b GCM-driven meteorology. The modeled SSM is converted to a multi-year average climatology for the five GCMs and compared to climatologies of SMOS and SMAP SSM (2015-2020) for the months March through October. Spatially averaged temporal skill metrics are shown in Fig. 2.
All GCM-driven simulations are similarly biased compared to the satellite products. The larger dry bias with SMOS (on average −0.05 m 3 m −3 ) compared to SMAP observations (on average −0.02 m 3 m −3 ) agrees with the shortterm evaluation results of the reanalysis-driven simulations (Sect. 4.1.1). The evaluation of predicted SSM compared to satellite data results in spatially averaged mean ubRMSDs ranging between 0.02 and 0.04 m 3 m −3 , with the lowest values for the multi-model median SSM (Fig. 2a), is shown. Figure 3 presents the spatiotemporal skill metrics comparing the multi-model median SIM2 SSM climatology with the two satellite-based SSM climatologies. The GCM-driven SSM climatology remains close to satellite SSM climatologies in drier conditions, but there is a wet model bias (or dry satellite retrieval bias) in wetter conditions (Fig. 3). Correlations between simulated climatologies and satellite data are slightly lower when considering individual GCMs (no median) with ranges of 0.41-0.45 and 0.47-0.51 for SMOS and SMAP, respectively (not shown). By design, GCM climatologies are unbiased against the reanalysis climatology, indicating that GCM-driven projections are representative of the re- Table 3. Spatial mean (± spatial standard deviation) of R, RMSD, bias, and ubRMSD between SIM1 SSM estimates, SMOS, and SMAP, for April 2015 through December 2016. analysis climate. From the evaluation of SIM1 and SIM2, it can be concluded overall that AquaCrop forced by ISIMIP3 input demonstrates a reasonable performance in terms of spatiotemporal SSM pattern representation; we therefore assume that the model can be used to project I net changes across the study area.

Climate impact on mean I net
The change in summer I net is assessed by the difference ( I net ) between the mean I net of the future horizons (2031-2060; 2071-2100) and the baseline period . In Fig. 4, spatial box plots of I net are presented for five GCMs individually and for the median across the GCMs for each scenario. Based on Fig. 4, increases in I net are expected in the future for all scenarios, where the severity of the increase de-pends on the emission scenario. SSP1-2.6 presents a stabilization of I net towards the end of the century, which is in line with the evolution of CO 2 for this scenario, whereas the other scenarios show increases from 2031-2060 to 2071-2100. The differences between the GCMs within an SSP are considerable, and these disparities increase with rising emission scenarios. According to the first GCM (GFDL-ESM4), on average about 7 mm per month extra irrigation water will be required in the summer months by 2050 for SSP5-8.5, whereas for UKESM1-0-LL, more than 20 mm per month will be required by mid-century for the same emission scenario. Decreases in I net (box plot whiskers below 0; Fig. 4c) are only observed in a few southern coastal locations under the high and severe emission scenarios. In these historically warm and dry regions with insignificant rainfall in the summer months, the effect of stomatal closure of 5 % in response to CO 2 concentrations above 550 ppm (parts per million) is stronger than the increase in ET (less than 5 %). These neg- ative differences are statistically non-significant (except for GFDL-ESM4, but the total area subjected to decreases is negligible). Figure 5 presents the spatial distribution of I net , for the median across the GCMs. Regions where all GCMs present significant changes are stippled. Once the results are presented in terms of medians, no statistically significant decrease in I net is observed.
Under the low-emission scenario ( Fig. 5a and d), the whole continent will face a mild increase in summer I net by about 13 mm per month (+18 %) in the near and far future, and regions undergoing severe increases cannot be identified. Towards the end of the century, for high and extreme emissions, the most affected areas (where all GCMs agree on a significant change) are situated in the central and southern latitudes (Fig. 5e and f). For the end of the century, the spatial mean summer I net increases by 22 and 26 mm per month (+30 % and +35 %) for SSP3-7.0 and SSP5-8.5, respectively. The most eastern parts are, on average, presenting large I net for the far future (2071-2100), but according to GFDL-ESM4 alone (not shown), these changes are nonsignificant and therefore not stippled in Fig. 5e and f. All SSP-GCM combinations agree on the evolution of I net in the northern Alps, where the situation is likely to remain stable in terms of amounts of required irrigation water. Figure 6 shows the spatial relationship between the expected change in summer I net with reference to the baseline period. Areas with historically extreme (> 150 mm per month) or low (< 20 mm per month) I net will not see their future needs increase drastically, whereas regions with a relatively moderate to high baseline I net will face the strongest changes. Table 4 summarizes the baseline summer I net and I net (median and standard deviation across GCMs), for six selected countries and the Benelux region included in this study area. The difference between the I net for various scenarios is of the same order of magnitude as, and often smaller than, the variability introduced by the various GCMs. Note again that the presented numbers are expressed in millimeters per month, but only averaged over 3 summer months, and the results are purely based on climate projections that are integrated into AquaCrop, assuming a hypothetical C3 crop, near-optimal fertilization, and without accounting for the presence or quality of the irrigation network. Figure 6 also shows how the atmospheric conditions in the summer, i.e., (P -ET 0 ), are directly related to I net . The largest increases in I net correlate with strong decreases in P -ET 0 . The few locations showing a positive (P -ET 0 ) (black crosses in Fig. 6a, b, and d) are still subjected to a slight increase in irrigation requirement. The I net estimates obtained with AquaCrop provide additional information over the mere (P -ET 0 ) estimates because the soil-plant system has a memory and temporally integrates the past P -ET 0 and irrigation events. Since the crop and management parameters are constant for the entire study domain, the only factor affecting I net for a given climate (P and ET 0 ) is the buffering capacity of the root zone, i.e., soil characteristics. An analysis of the influence of soil characteristics showed, for instance, that sandy soils see their I net enlarge more rapidly compared to loamy soils. However, no clear conclusions could be drawn because the vast majority of Europe at the resolution of this study is dominated by a loamy soil texture.

Climate impact on the interannual variability in I net (RI net )
To assess the potential change in interannual variability in summer I net , the difference between the maximum and minimum summer I net within a 30-year time period (range of I net = RI net ) is evaluated. The future RI net values are assessed with reference to the baseline period, resulting in RI net for each scenario and GCM. Results are presented in Fig. 7, where expansions of RI net are indicated in red and reductions in blue.
For all SSPs, future RI net are likely to decrease in most of southern Europe, whereas the gap between the highest and lowest irrigation requirement in the 30-year time window is expected to grow in northern and central regions of Europe. Similar to I net (Fig. 5), Fig. 7 shows that changes   Table 4 summarizes the baseline RI net and changes in interannual variability for some selected countries in Europe. Figure 8 presents how the change in interannual variability ( RI net ) of the two future periods relates to the baseline RI net and to I net . Regions with severe increases in I net do not necessarily present the highest enlargements in RI net . The largest baseline RI net values correlate to lower RI net for the far future (SSP3-7.0 and SSP5-8.5; Fig. 8e and f), in combination with high values of I net (dark blue dots; Fig. 8e and f). In other words, the Mediterranean region, western France, and the region around the Black Sea, currently with a high interannual variability in irrigation requirements, will see their requirement significantly increase to more steady high irrigation requirement. Large RI net values follow in the Carpathian Mountains (central Europe) for SSP5-8.5 (Fig. 7f). According to the model, only a little irrigation was required in these mountainous regions during the baseline, whereas future requirements are projected to increase. In the future, I net peaks to larger values for several years, increasing RI net .
To obtain a better understanding of changes in the interannual variability in I net , time series at two different locations for one GCM are presented in Fig. 9. Figure 9a shows the evolution of summer I net for a pixel in central-western Europe, with a RI net of 100 mm per month (for IPSL-CM6A-LR, a randomly chosen GCM). During the baseline period, summer I net fluctuates between zero and about 35 mm per month, while at the end of the century, the maximum I net of the time window will reach 135 mm per month for SSP5-8.5 with the same minimum I net as for the baseline. For the second pixel in southern Europe (Fig. 9b), a stabilization of the yearly summer requirement is expected. Overall, more water will be required here, but summer I net will not vary much relative to the average requirement from 1 year to another. This second location results in a decrease in RI net of about 55 mm per month for the presented GCM under SSP5-8.5. Figure 5. Changes in summer I net ( I net ) (mm per month), median across five GCMs for the two future time horizons (rows) and the three scenarios (columns) with reference to the baseline period. The stippled areas represent pixels where all five GCMs present statistically significant changes (t test; p < 0.05). Figure 6. Scatterplots of I net relative to the baseline period for the two future periods (rows) and the three scenarios (columns). The coloring refers to the corresponding (P -ET 0 ). Increases in (P -ET 0 ) are represented by pink crosses. All values (I net and (P -ET 0 )) are medians across the GCMs (mm per month).

A new model setup for climate change impact assessment
The regional setup of the AquaCrop model using ISIMIP3 meteorological data has the potential to assess impacts of climate change on the irrigation requirement and possibly also on future crop production. First, the short-term evalu-ation proved that the model forced with reanalysis meteorology (ISIMIP3a) has an acceptable performance, i.e., the ubRMSD between SIM1 SSM simulations and satellite retrievals is 0.06 and 0.08 m 3 m −3 for SMAP and SMOS, respectively (Table 3). The lower model performance compared to SMOS SSM could be due to remaining radio frequency interference contamination (Oliva et al., 2012). It is important to note that the satellite target uncertainty is 0.04 m 3 m −3 over areas with less than 5 kg m −2 vegetation Table 4. Median across GCMs (± standard deviation of GCMs) of baseline summer I net , I net , baseline RI net , and RI net (mm per month), spatially averaged over the country, for six European countries and the Benelux region. The changes are presented for the two future horizons (2031-2060 and 2071-2100, columns) and for the different emission scenarios (line 1, 2, 3 of a cell corresponding to SSP1-2.6, SSP3-7.0, and SSP5-8.5, respectively).  Figure 7. Future changes in RI net ( RI net ) (mm per month) median of five GCMs for the two future horizons (rows) and the three scenarios (columns) with reference to the baseline period.  water (i.e., excluding dense vegetation; Entekhabi et al., 2014). Even though a conservative screening was used, this target may be exceeded at some times and locations. The findings of the evaluation against SMAP SSM are comparable to the results found by de Roos et al. (2021b). However, the latter study used Modern-Era Retrospective analysis for Research and Applications (version 2; MERRA-2) input and showed a slightly higher performance between the simulated SSM of the regional AquaCrop model and SMAP SSM. Furthermore, the difference in study domain and especially in resolution play a major role in explaining this difference (larger domain and soil characteristics aggregated to coarser pixels in this study).
The strong agreement between the SIM2 SSM climatologies obtained with GCM-driven input (ISIMIP3b) and the reference satellite SSM climatologies (Fig. 2) further confirmed that the historical GCM-driven input is also reliable. A larger bias was observed in wetter moisture conditions ( Fig. 3), possibly coming from the model itself or from biases in satellite retrievals. Overall, the provided simulated atmospheric data could represent the main variations in SSM for the past (2015-2020) and is therefore reliable to be used for climate change assessments.
The model simulations for the historical evaluation did not include any irrigation, and therefore, some mismatches in SSM could possibly be expected in areas which are currently irrigated. Earlier studies suggest that the contrast between satellite observations and model simulation could identify unmodeled processes (Brocca et al., 2018). In a separate analysis (not shown), this effect was assessed by evaluating the correlation values with regard to irrigation areas that are equipped for irrigation (AEI; similar to de Roos et al., 2021b). By using the FAO global maps of irrigated areas version 5  aggregated to ISIMIP resolution, pixels were divided into the following two groups: (1) less than 10 % of the area is equipped for irrigation and (2) more than 10 % of the pixel is equipped. The correlations between AquaCrop SIM1 SSM and satellite SSM for the irrigated pixels (AEI > 10 %) were nearly identical to correlations for the locations with an AEI < 10 % (0.65 versus 0.67 for SMOS and 0.51 for both classes of AEI for SMAP), therefore not revealing where irrigation was missing in the simulations. Even if irrigation could be captured by observation-based SSM (Kim et al., 2020), the low amount of reference satellite data in regions presenting high percentages of AEI compromised the evaluation in this study. The use of a conservative screening of satellite SSM retrievals for both SMOS and SMAP resulted in a significant amount of data loss, especially in densely vegetated areas which are generally masked out (Kim et al., 2020). Furthermore, irrigated areas are usually much smaller than the considered 0.5 • pixel, and a recent study stressed the low potential to detect irrigation at coarse resolutions (Dari et al., 2021). Nevertheless, the second time series presented in Fig. 1b showed an underestimation of SSM during summer by the model and large differences between SMOS and SMAP, both suggesting potential irrigation applications, as confirmed by a high AEI percentage.

Future mean and interannual variability in summer I net
The evolution of future summer I net with climate change is highly dependent on the scenario (SSP) but also on the GCM (Table 4; Fig. 4). Our results agree with several earlier droughts and irrigation projection assessments (Döll, 2002;Elliott et al., 2014;Konzmann et al., 2012;Pfister et al., 2011;Ruosteenoja et al., 2018;Satoh et al., 2021;Schaldach et al., 2012;Wada et al., 2013). Under high and extreme emission scenarios, the whole continent will be significantly impacted by the end of the century ( Fig. 5e and f), with the most drastic changes in central and upper southern latitudes of the study domain, confirmed by the high increases in meteorological (Spinoni et al., 2018) and soil water (Ruosteenoja et al., 2018) shortages in these regions. I net spatial patterns (Fig. 5) are comparable to the findings of Konzmann et al. (2012) and Wada et al. (2013) for the most affected areas where all GCMs present significant changes. Eastern Europe shows on average large positive I net values, but not all GCMs converge towards significant changes in this region. Moreover, the model evaluation showed a lower performance over this area (Fig. 1a). For these two reasons, the results may be less certain.
For the year 2050, Schaldach et al. (2012) estimated an average increase of 70 mm yr −1 (+15 % compared to the baseline of 470 mm yr −1 ) in the I net over the European continent under a high-emission scenario. For the same scenario, Fischer et al. (2007) predicted an increase of 53 mm yr −1 (+36 % compared to the baseline of 147 mm yr −1 ) over western Europe by 2080. These values are comparable to the results of this study, presenting a mean increase of 38 mm yr −1 (+17 % compared to the baseline of 221 mm yr −1 and integrated over the 3 summer months only) and 67 mm yr −1 (+30 %) for 2050 and 2090, respectively, and for a comparable emission scenario (SSP3-7.0). Also, the evolution of I net under SSP5-8.5 (+35 %) can be related to the findings in Wada et al. (2013), for which increases larger than 25 % are estimated over almost the entire European continent. Under this same scenario, the results show that France will face extreme increases in I net (+75 %), as confirmed by Fader et al. (2016) (+80 %). Smaller increases were found in other studies (e.g., Döll, 2002;Elliott et al., 2014); however, absolute values are hard to directly confront in the literature because of the differences in methodology. In the literature, I net is often assessed under the assumption of potential irrigation during the entire year or growing season (as opposed to the summer only in this study) and considering other factors such as irrigation efficiencies and strategies, varying crop types, and even population increase or economic growth ultimately impacting, for example, irrigation efficiencies. Furthermore, Elliott et al. (2014) demonstrated that irrigation estimates from GHMs deviate strongly from crop model predictions, potentially due to the differences in agro-hydrological processes between the two types of models. Also, Wada et al. (2013) proved that the largest part of uncertainty in future I net estimation is first due to the impact model and thereafter to climate uncertainty. ET 0 is a determinant factor for these kinds of studies, and its calculation procedures can have an important influence on the final results (Webber et al., 2016).
Atmospheric data alone could give an indication of the crop water requirement, as is done in meteorological drought assessments. However, the integration of P and ET 0 into a crop model with application of irrigation is more realistic to estimate I net because it benefits from the land system memory. It should be noted though that the wetness of the irrigated land area will in turn affect turbulent fluxes and thus atmospheric variables in general Thiery et al., 2017Keune et al., 2018). This feedback loop is not in-cluded in the presented simulations and needs to be carefully considered in future attempts to design climate-smart irrigation systems.
Whereas the focus of this study was on the irrigation requirement, a similar analysis can be performed in terms of agricultural productivity. An increase in I net is expected but, following increasing CO 2 concentrations, biomass production is also expected to increase (Schleussner et al., 2018;Vanuytrecht, 2020;Vanuytrecht et al., 2012). The yield water productivity (WP Y/ET ; i.e., the ratio between crop yield and the amount of water lost by evapotranspiration) will improve due to the rising CO 2 concentrations. Since crops can only fully profit of the CO 2 fertilization when soil fertility is high , an increase in WP Y/ET is likely to occur in irrigated fields that are generally well fertilized. In the absence of soil water and soil fertility stress, crop production might increase by about 25 % up to 45 % for an atmospheric CO 2 concentration of 550 ppm . Effects above this concentration remain more uncertain.

Future adaptations of irrigation infrastructure and management
Different practical future pathways can be considered starting from the current state of irrigation requirement. In regions where I net is currently low (low baseline I net ), there is typically no irrigation infrastructure available or needed to achieve a fairly high crop production. However, to maintain crop production in the future, large investments will be required to develop or extend the irrigation infrastructure . In regions with an existing water shortage and irrigation infrastructure, the focus will be on improving irrigation efficiencies, aiming to buffer the effects of climate change (Jägermeyr et al., 2016). Our study did not consider specific irrigation practices and efficiencies. The latter have been estimated at around 50 % in Europe (Fischer et al., 2007;Rohwer et al., 2007;Wada et al., 2013), meaning that the I net values presented in this study could roughly be doubled to obtain the gross requirement. Sprinkler irrigation remains the most widely used practice in Europe, but the share of drip irrigation is progressively increasing in southern countries such as Spain and Italy (Monaghan et al., 2013), aiming at improving irrigation efficiency. Furthermore, with a lower availability of freshwater, the introduction of other irrigation strategies, such as deficit irrigation, also gain importance. Deficit irrigation intends to maximize crop water productivity, therefore stabilizing crop yields through time (Geerts and Raes, 2009;Mushtaq and Moghaddasi, 2011).

Model uncertainty
Model uncertainty is an important factor influencing climate scenario analyses (Lehner et al., 2020). This starts with the high variability between climate scenarios that are input to the crop model simulations. The uncertainty of future climate was included by using meteorological input from three scenarios and five GCMs, resulting in 15 different SSP-GCM combinations. The process of using only a small fraction of the various existing GCMs has been criticized (McSweeney and Jones, 2016). However, previous drought and irrigation projections often used fewer than five GCMs or used more but for only one emission scenario. Additionally, the ISIMIP GCMs are carefully selected to represent the entire CMIP ensemble Warszawski et al., 2014). The AquaCrop model setup also adds uncertainty. First, the constantly evolving field practices in terms of, for example, crop type and cultivars, water management, and soil fertility management were not included in the model simulations. However, this aspect is almost impossible to include. Second, the model generalizations (generic C3-type of crop, unconstrained water availability, and constant small soil fertility stress for the whole domain) increase the uncertainty in the projections. It should be noted that actual area of irrigated land is not considered, and consequently, the expansion thereof is not simulated (estimated by, e.g., Schaldach et al., 2012). Nevertheless, the intention of this study is to limit the uncertainty in time and space (as described in Sect. 2.2) by assuming these generalizations, therefore aiming to present in a simple way the evolution of I net during summer months.

Conclusions
Large-scale AquaCrop simulations over Europe were performed using ISIMIP3 meteorological forcings at a spatial resolution of 0.5 • lat × 0.5 • long to assess future changes in net irrigation requirements I net . Because this is the first large-scale AquaCrop application with ISIMIP3 input, the model was first evaluated using satellite-based SSM. The reanalysis-driven (ISIMIP3a) simulated SSM have a mean spatial ubRMSD of 0.06 m 3 m −3 with SMAP retrievals, and thereby deviate slightly more than the assumed intrinsic error of the satellite retrieval error (0.04 m 3 m −3 ). The performance of AquaCrop compared with SMOS (ubRMSD = 0.08 m 3 m −3 ) is slightly lower than with SMAP, most likely because the SMOS sensor suffers more from radio frequency interference. When using GCM-driven (ISIMIP3b) meteorology as input, the multi-year average SSM of the simulations is comparable to that of reference satellite data (ubRMSD = 0.03 m 3 m −3 ), which reinforces the reliability of the ISIMIP3 climate data for future projections.
In the second part of this paper, the summer irrigation requirement of a near (2031-2060) and far (2071-2100) future horizon was simulated using five different GCMs and three emissions scenarios. We present net irrigation requirement values that are independent from the irrigated area, period, infrastructure, and the exact crop type. The mean and interannual variability in net irrigation requirement I net for the summer months were quantified for the two future climate horizons and compared to the baseline period . This evaluation showed that the effect of climate change on future I net depends on the emission scenario, but more strongly on the GCM. Under high and extreme emission scenarios (SSP3-7.0 and SSP5-8.5), almost the whole European continent will see an increase in summer I net , with on average 30 % and 35 % additional net irrigation water required in the far future relative to the baseline I net . Especially regions with a moderate baseline I net will experience strong increases in I net . All GCMs agree on significant increases in central to southern Europe, which is in line with meteorological and soil moisture drought projections for the same scenarios and previous irrigation demand projections.
The interannual variability in summer I net was quantified by the range between maximum and minimum I net within the 30 year climate periods, RI net . It was found that mild increases in I net result in larger gaps between maximum and minimum summer I net within a time window, corresponding to more extremes, and a high interannual variability (large RI net ). In the future, northern and central areas will face increased RI net , whereas southern Europe is likely to see the variability diminish, resulting in steady high I net . Under the strong mitigation scenario (SSP1-2.6), I net stabilizes towards the end of the century, consistent with the plateauing CO 2 concentrations in this scenario. The increase in variability is also reduced under this scenario. Overall, extra water will be required, but more production can be achieved under higher CO 2 concentrations. The exact effect of CO 2 fertilization remains uncertain, but it is expected that the yield, and especially yield water productivity, is likely to increase in the future in the absence of water and soil fertility stress. Our large-scale setup with AquaCrop is well suited to exploring the effect of climate scenarios on crop productivity in future research.
These results highlight the importance of climate change mitigation to keep future irrigation at reasonable levels, while it also stresses the high uncertainty of climate projections. This study aimed to demonstrate the effect of climate change on I net over Europe, without considering land use, crop types, and actual irrigated areas to avoid the inclusion of more uncertainty. Therefore, the results of this study should not be taken as predictions but as an indication of the potential consequences of climate change on the amount of and variability in I net for the summer months.

Appendix A: Comparison of SIM1 and SMOS SSM in terms of anomalies
To evaluate the short-term and interannual variability in the AquaCrop SSM in terms of temporal anomaly correlation (anomR), time series of anomalies are calculated by subtracting the climatology from simulation and satellite data for each daily time step. The climatology calculates the mean seasonal cycle as a long-term mean using a sliding window of 31 d, with a minimum threshold of three data points of data within the window.
Because anomR values can only be computed when several years of data are available, the AquaCrop SSM simulations forced with ISIMIP3a reanalysis data for the years 2011-2016 (SIM1) are evaluated against SMOS SSM retrievals only. Figure A1 shows the anomR over Europe with a spatial mean anomR of 0.44. Higher correlations are found in southwestern locations (anomR is often > 0.6) and lower performances occur in north and central-western Europe (anomR is generally < 0.4). Also shown in this figure is a partitioning of Europe in various zones for further discussion.
Author contributions. LB adapted the code to run the regional version of AquaCrop with ISIMIP data, prepared the input data, conducted the model evaluation, and performed all simulations and analyses. SdR provided the code of the regional version of AquaCrop (v6.1) and scientific guidance. GJMDL prioritized the main steps taken in the paper, provided supervision and scientific guidance throughout all research advances, and managed the HPC usage. WT provided scientific guidance through climate change impact assessments and ISIMIP. DR provided scientific guidance regarding the use and interpretation of AquaCrop, along with an appropriate methodology to assess the future irrigation requirement. LB wrote the paper, and all authors contributed.
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 in published maps and institutional affiliations.
Acknowledgements. This research is conducted as part of the H2020 project SHui, which stands for "Soil Hydrology research platform underpinning innovation to manage water scarcity in European and Chinese cropping systems". The resources and services used in this work were provided by the VSC (Flemish Supercomputer Center), funded by the Research Foundation -Flanders (FWO) and the Flemish government. The authors would like to acknowledge the ISIMIP, for providing climate input data used in this study. We would also like to thank Luke Grant, for his help with the ISIMIP data downloads and storage management. The authors appreciate the constructive reviews from the two anonymous reviewers.
Financial support. This research has been supported by SHui, which is funded by a European Union project (grant no. GA 773903), with additional support from the KU Leuven internal fund (grant no. C14/21/057).
Review statement. This paper was edited by Daniel Viviroli and reviewed by two anonymous referees.