Dynamic Crop Models and Remote Sensing Irrigation Decision Support Systems: A Review of Water Stress Concepts for Improved Estimation of Water Requirements

: Novel technologies for estimating crop water needs include mainly remote sensing evapotranspiration estimates and decision support systems (DSS) for irrigation scheduling. This work provides several examples of these approaches, that have been adjusted and modiﬁed over the years to provide a better representation of the soil–plant–atmosphere continuum and overcome their limitations. Dynamic crop simulation models synthetize in a formal way the relevant knowledge on the causal relationships between agroecosystem components. Among these, plant–water–soil relationships, water stress and its e ﬀ ects on crop growth and development. Crop models can be categorized into (i) water-driven and (ii) radiation-driven, depending on the main variable governing crop growth. Water stress is calculated starting from (i) soil water content or (ii) transpiration deﬁcit. The stress a ﬀ ects relevant features of plant growth and development in a similar way in most models: leaf expansion is the most sensitive process and is usually not considered when planning irrigation, even though prolonged water stress during canopy development can consistently reduce light interception by leaves; stomatal closure reduces transpiration, directly a ﬀ ecting dry matter accumulation and therefore being of paramount importance for irrigation scheduling; senescence rate can also be increased by severe water stress. The mechanistic concepts of crop models can be used to improve existing simpler methods currently integrated in irrigation management DSS, provide continuous simulations of crop and water dynamics over time and set predictions of future plant–water interactions. Crop models can also be used as a platform for integrating information from various sources (e.g., with data assimilation) into process-based simulations.


Introduction and Scope
Agriculture has for long been the first consumer of freshwaters, with agricultural irrigation contributing for about 70% of global human freshwater withdrawal [? ]. At the present time, 57% of blue water (surface and groundwater) usage is unsustainable, meaning that it exceeds the available renewable blue water: around 70% of the unsustainable portion of the global blue water is due to the production of wheat, cotton, sugarcane and fodder [? ]. Considering also that irrigated agriculture covers only 20% of the world's cultivated lands, but is responsible for 40% of crop production [? ], it should be considered paramount to properly manage water resources for irrigated agriculture to obtain "more crop per drop", a term that dates back as early as the nineties, as pointed out by Giordano et al. [? ]. This means reasoning in terms of water productivity rather than merely on crop yield maximization,

•
Section ?? dealing with RS-related products used for direct or indirect estimation of crop water status, focusing in particular on approaches to model evapotranspiration and to support farmers' decision for irrigation scheduling. The section also highlights the limits of these methods and tracks the evolution of some of these examples over time, with reference to the implementation of additional WS indices/functions borrowed from crop models' formalism or conceptualization. • Section ?? aggregates dynamic crop models by similar conceptualization of the plant-water relationships, describing the main approaches used to simulate plant growth and development, the state variable used as input for WS calculations and the processes targeted by WS, in order to provide readers with useful reference for investigating approaches to model WS and schedule irrigation. Considerations raised by model comparison and perspectives of integrating RS with crop models are also addressed.

Modeling Evapotranspiration with Remote Sensing
The FAO (Food and Agricultural Organization) Irrigation and Drainage paper no. 56 [5] has become a standard reference for the estimation of crop water requirements based on evapotranspiration (ET) assessment, which can be easily adapted to RS data [21]. Crop water requirements (CWR) are defined in FAO Irrigation and Drainage paper no. 24 [22] as "the depth of water needed to meet the water loss through evapotranspiration (ETc) of a disease-free crop, growing in large fields under non-restricting soil conditions including soil water and fertility and achieving full production potential under the given growing environment". In this regard, the FAO-24 paper recommends considering three aspects of the agroecosystem to estimate CWR: (i) the effect of climate on the ET demand, (ii) the effect of crop characteristics on modulating ET and (iii) the effect of local conditions (e.g., soil characteristics) and agricultural practices (e.g., irrigation system) on determining soil water availability.
In this regard, some fundamental concepts and terminologies of the FAO-56 method are necessarily to be briefly recalled here, since they are widely used in RS, DSS and crop modeling communities. The FAO-56 method prescribes the use of the Penman-Monteith equation to estimate reference ET (ET0), with subsequent adjustments to obtain estimates of crop potential ET (ETp) and actual evapotranspiration (ETa). ET0 is the evapotranspiration of a reference hypothetical standard unstressed grass surface. ETp is the evapotranspiration of a specific crop, unstressed, and ETa is the

Modeling Evapotranspiration with Remote Sensing
The FAO (Food and Agricultural Organization) Irrigation and Drainage paper no. 56 [? ] has become a standard reference for the estimation of crop water requirements based on evapotranspiration (ET) assessment, which can be easily adapted to RS data [? ]. Crop water requirements (CWR) are defined in FAO Irrigation and Drainage paper no. 24 [? ] as "the depth of water needed to meet the water loss through evapotranspiration (ETc) of a disease-free crop, growing in large fields under non-restricting soil conditions including soil water and fertility and achieving full production potential under the given growing environment". In this regard, the FAO-24 paper recommends considering three aspects of the agroecosystem to estimate CWR: (i) the effect of climate on the ET demand, (ii) the effect of crop characteristics on modulating ET and (iii) the effect of local conditions (e.g., soil characteristics) and agricultural practices (e.g., irrigation system) on determining soil water availability.
In this regard, some fundamental concepts and terminologies of the FAO-56 method are necessarily to be briefly recalled here, since they are widely used in RS, DSS and crop modeling communities. The FAO-56 method prescribes the use of the Penman-Monteith equation to estimate reference ET (ET0), with subsequent adjustments to obtain estimates of crop potential ET (ETp) and actual evapotranspiration (ETa). ET0 is the evapotranspiration of a reference hypothetical standard unstressed grass surface. ETp is the evapotranspiration of a specific crop, unstressed, and ETa is the Remote Sens. 2020, 12, 3945 5 of ?? evapotranspiration of a crop under the real field conditions, accounting for WS. Equations (1) and (2) represent the FAO-56 method [? ] for calculating potential and actual ET, respectively.
where the 2 alternatives depend on the method used: the "single coefficient" approach (the first ones) consider the crop field as a whole, while the "dual coefficient" approach (the second ones) separates crop transpiration from soil evaporation. In Equation (1) and (2)  , of which we provide only a brief but necessary recap. The first two approaches use RS to estimate parameters for ET calculations following the FAO-56 approach.
Reflectance-based Kcb are estimated from VIs and used in the crop dual coefficient approach. The recent review of Pôças et al. [? ] on the relationships between crop coefficients and VIs highlighted that, due to the advances of this approach over the last 30 years, such relationships constitute simple and robust operational tools for CWR estimations. Most frequently, these are empirical linear or non-linear relationships between the normalized difference vegetation index (NDVI), the soil adjusted vegetation index (SAVI) and Kc or Kcb. The authors clarified that, due to their empirical nature, these relationships vary with experimental conditions (crop type, irrigation management and climatic conditions), and integration with other approaches (like soil water balance models) are advisable.
RS-based Penman-Monteith direct approaches use RS data to directly estimate the parameters of the Penman-Monteith ETp equation. The rationale behind this approach is that the crop coefficient is a proxy of the parameters describing the canopy development [? ? ]. D'Urso and Menenti [? ] originally developed the approach by defining analytical functions for the relationships between surface reflectance and vegetation parameters (LAI, surface albedo and crop height) used for the calculation of ETp according to the Penman-Monteith equation. Albedo is used to estimate the net shortwave radiation and LAI and crop height are used in the estimation of canopy resistance. The other weather-and climate-related parameters are derived from air temperature, relative humidity, wind speed and incoming radiation, as usual practice for the Penman-Monteith equation. The equation used in the Penman-Monteith direct method can be expressed as follows (Equation (3) [? ? ? ]): where λ, ∆, R nl , G, ρ a , c p , D E and γ are weather-related variables: λ is the latent heat of vaporization, ∆ the slope of vapor pressure, R nl the long-wave net radiation, G the soil heat flux, ρ a the atmospheric density, c p the specific heat capacity of air, D E the vapor pressure deficit and γ the psychrometric constant. R ns , r a , r s , are canopy-related variables derived from RS: R ns is the shortwave net radiation, r a the air and r s the surface resistance. These variables are derived from albedo, LAI and crop height through Equations (4)-(6), respectively: where R ts is the total incoming shortwave radiation, the surface albedo α is the spectrally integrated hemispherical solar reflectance, r l (the average daily value of single leaf resistance) and LAI are used to calculate the surface resistance, while crop height (h), z u and z h (measurement height of wind and humidity) and U (wind speed) are used to calculate the aerodynamic resistance. Optical data is used to estimate the required canopy biophysical variables, for example through the use of RS products of the Landsat constellation, SPOT4, ASTER [? ] or Sentinel 2 [? ]. RS-based surface energy balance (SEB) approaches estimate ET as a residual term of the surface energy balance equation, using thermal data to compute heat fluxes. The concept behind SEBs is that the net energy absorbed by the surface (Rn) should be balanced by the sensible, latent and ground heat fluxes [? ]. The balance can also be expressed to explicit the latent heat flux as in Equation (7): where LE is the latent heat flux (flux of energy from the Earth's surface to the atmosphere, that is associated with water evapotranspiration), Rn is the net radiation flux, G the soil heat flux and H the sensible heat flux. Rn consists of the sum of the net shortwave radiation flux (R ns , determined as in Equation (4)) and the net longwave radiation flux (difference between incident and emitted longwave radiation). Heat fluxes are determined using a wide variety of approaches, that span from empirical direct methods that use locally calibrated constants to simplify the calculations, through methods that merge empirical and mechanistically-based parameterization, up to more deterministic methods [? ]. We invite the readers to refer to the review of Courault et al.
[? ] for a deeper description of the methods used by SEB models to estimate energy fluxes, while some relevant features and differences among SEBs are worth to be mentioned here. A key step of all the SEB models is the estimation of land surface temperature (LST) from thermal imagery for the calculation of the sensible heat flux, while optical imagery in the visible and NIR spectra are used to estimate reflectance-derived canopy characteristics [? ]. The loss of sensible heat flux from the surface is determined by aerodynamic temperature, that is usually not directly measured but needs to be derived from RS-obtained radiative temperature [? ]. The way this is done defines the main difference between various SEB models, e.g., SEBAL (surface energy balance algorithm for land), TSEB (two source energy balance) and METRIC (mapping evapotranspiration at high resolution with internalized calibration). SEBAL works around the problem of the computation of the relationship between aerodynamic and radiative temperature through the identification of a "wet" and a "dry" pixel in the study area. At the wet surfaces, water vapor is released at a rate that roughly corresponds to 100% of the atmospheric demand, while at dry surfaces the latent heat flux is considered approximately null: this gradient permits the calculation of the sensible heat flux using the spatially-distributed radiometric temperature obtained from RS products [? ? ? ]. The TSEB model (two source energy balance) uses a two-source parameterization to account for the contribution of soil and canopy to the heat fluxes and reduce the error of single-source models when applied to sparsely vegetated surfaces [? ? ]. METRIC has its foundation on SEBAL, but calibrates internally the model on ground-measured ET0 to reduce the bias in the computation of the energy balance and provide consistency with the traditional methods for ET0 estimations used in agricultural sciences [? ]. Overall, Calera et al. [? ] provided several examples of successful use of the above-mentioned approaches for ET estimations, but also pointed out their limitations, which are mainly related to: uncertainty in the relationships of reflectance-derived VIs with Kcb or ET equation parameters (for the first two methods), uncertainty in the role of soil evaporation (for the direct ETp approach) and uncertainty in the representativeness of ET estimates over time in the SEB approach, due to the instantaneous nature of RS estimates of ET, corresponding to the timing of satellite overpasses. On the other hand, more accurate estimations of canopy biophysical variables (e.g., albedo, FCOVER and 7 of ?? LAI) and more stable relationships with crop coefficients are also more likely to be provided, given the increasing availability of physically-based approaches like radiative transfer models and hybrid methods that combine these models with nonparametric regression and machine learning algorithms [? ]. These methods do not present the limits of empirical relationships between VIs and biophysical variables, having a more general validity and being able to make full use of the information stored in the entire radiation spectrum, whereas VIs consider only a few spectral bands.

Decision Support Systems for Irrigation Scheduling
CWR estimation and irrigation planning can be carried out using decision support systems (DSSs), that simplify features of crop models and integrate them with communication tools to manage inputs and outputs. In practice, a DSS often takes the form of a software that aims at assisting farmers or technicians with decision making on a specific aspect of crop management [? ]. To acquire data, DSSs for irrigation scheduling may use either online data (mainly weather data, but also crop and soil data) or data from sensors providing real time information (e.g., from weather stations or soil moisture sensors) [? ]. A balance between simplicity and accuracy needs to be struck, since greater accuracy means that more data is required by the model. DSS come in a great variety of forms and complexity, but usually rely on estimates of ET or soil moisture, using water balance modules. The use of crop and soil simulation modules allow DSSs to provide continuous data on CWR and WS, in order to plan irrigation intervention beforehand and/or adjust irrigation amount based on real-time information. RS data can also be integrated into DSSs to improve predictions. RS permits to collect information on soil moisture, plant growth and plant water status, avoiding at the same time costly and time-consuming field work and the use of more expensive monitoring technologies, with the added benefit of providing spatialized data. This information can be used parameterize the DSS simulation modules or to provide real-time adjustments, in a similar way to more complex crop models (e.g., through data assimilation, discussed in Section ??). Recent advances in computer science have also led to the development of DSSs relying on machine learning approaches, that do not need to simulate the behavior of the agroecosystem, but rather find patterns on the collected data (e.g., weather and soil data) and suggest actions to be taken to schedule irrigation [? ]. Integrating field, proximal sensing and RS data with DSSs has the advantage of providing more accurate and continuous information on crop growth and development, soil water content, CWR and WS. Indeed, Cahn and Johnson [? ] pointed out several limitations of measurements, devices and methods used without a comprehensive look at the crop system. Among them we can summarize some that are particularly relevant for this work: (i) soil moisture sensors alone can be useful for determining irrigation timing, but are less useful to determine the irrigation amount, (ii) ET estimates alone may be impractical for scheduling irrigation if not backed up by dynamic assessment of the crop coefficients over time, (iii) soil evaporation calculations require knowledge of soil properties, irrigation methods and interval between irrigation and rainfall events, and (iv) information about the irrigation system (i.e., wetting area, irrigation efficiency, application rate, etc.) needs to be integrated into water requirements calculations.

Towards Greater Accuracy in Crop Water Stress and Water Requirements Estimations
Some of the approaches mentioned in previous subsections, like RS-based ET estimates or DSS that integrate RS data into simple water balances, have been adjusted and improved over the years to obtain more accurate estimates of WS and CWR. In general, the trend was to increase complexity that could span from the simple addition of CWR indices or thresholds to the integration of entire new modules for better predictions. This tendency can be considered as a logical development aimed at providing an answer to overcome the previously mentioned limitations. Here are some examples: proposed to integrate RS data with a simple radiation-based approach for the calculation of biomass, depending on fAPAR and on RUE (see Section ?? for definitions and details). The authors used a compact airborne spectrographic imager (CASI) mounted on an airborne platform to obtain hyperspectral RS data during the growing season. In addition, Landsat 5 and Landsat 7 cloud-free images were employed. Surface reflectance from these sensors (visible and NIR range) was combined and used to calculate 2 VIs: NDVI and MTVI2 (modified triangular vegetation index 2). fAPAR was estimated from MTVI2. The canopy structure dynamics model (CSDM) was used to simulate fAPAR dynamics over the season depending on temperature. CSDM was fitted to fAPAR estimated from RS, to derive fAPAR seasonal trend used for the seasonal estimation of biomass accumulation through the radiation driven Monteith approach. In addition, a crop stress index varying between 0 and 1 was added to act as a modifier of the biomass-radiation equation, nullifying the equation when the crop is completely stressed and leaving it unchanged when the crop is unstressed. This approach integrated RS data with simulation of crop cover dynamics, using the radiation-biomass relationship for estimating biomass and yield and a water stress index to limit crop growth, constituting an upgrade of methods based directly on RS data only.
[? ] used VIs obtained from RS data to estimate fAPAR and Kt, a transpiration coefficient similar to Kcb, in order to use these as inputs of simple radiation-driven and water-driven models, respectively (details on the two approaches are provided in Section ??). fAPAR and Kt were obtained from relationships with VIs, but the authors contemplated that other analytical approaches could be used to obtain the two biophysical parameters. NDVI obtained by Landsat 5, 7 and 8 were interpolated to obtain daily values to use in the estimation of fAPAR and Kt. Biomass accumulation was then simulated using the radiation-driven and water-driven models. • Olivera-Guerra et al.
[? ] developed a method to retrieve irrigation timing and amount at field scale from RS data. Landsat 7 and 8 images with <30% of cloud cover were used over 4 cropping seasons, with an average of 20 images per season. Optical and thermal data from these sensors were used to estimate LST and FCOVER. LST was estimated using the revised single channel algorithm [? ], using the thermal band of Landsat, the atmospheric water vapor content from Modis 6.0 and the spectral surface emissivity estimated using the NDVI threshold method proposed by Sobrino et al. [? ], that uses FCOVER to weight emissivity from soil and vegetation. Soil emissivity was acquired from ASTER-GED using bands 13 and 14. Vegetation cover was estimated linearly based on minimum and maximum NDVI, following Duchemin et al. [? ], and NDVI was obtained from Landsat red and infrared bands. A first guess of root zone soil moisture was obtained for each Landsat overpass date using LST and FCOVER, with a partitioning method based on the LST-FCOVER feature space (e.g., as in Jiang and Islam [? ]). The FAO-56 dual coefficient method water balance was applied from the last soil moisture guess, backwards, in a daily time step (recursive mode). The first time that soil moisture reached field capacity (FC) was considered as the day after irrigation was applied, and using the water balance in forward mode from the previous Landsat overpass it was possible to estimate the soil moisture before irrigation started. Since irrigation was assumed to bring soil moisture at FC, irrigation amount could be calculated. In this way, it was possible to reconstruct the soil water and irrigation dynamics (timing and amount) of the entire season of specific fields (by pixel aggregation), providing estimates of water consumption and of crop water requirements when field data are scarce or not available.
In general, in the context of RS and DSS technologies for CWR estimation, WS assessment and irrigation planning, there seems to be a trend on the part of the scientific community to improve these technologies at the cost of added complexity, in order to achieve better characterization of plant-water relationships. In Section ?? we present a categorization of the formalism used by more complex crop models to simulate CWR and WS. This categorization does not aim at providing an exhaustive description of the plant-water relationships of each model, for which we invite the readers to refer to previous studies of Saseendran et al. [? ] and Dhakar et al. [? ]. It is rather aimed at highlighting common and well-established features that can be useful when improving RS and DSS WS and CWR estimates. Crop models represent a synthesis of the current knowledge about crop processes [? ], i.e., an up to date storage of mathematically summarized soil-plant-atmosphere processes that should be considered relevant in the agricultural context. This should also prove useful for the implementation of deficit irrigation, where it is necessary to more precisely quantify also the effects of WS on yield production. As pointed out by Jovanovic et al. [? ], deficit irrigation is a complex practice that involves many variables, and optimizing modeling and DSS with holistic approaches are paramount to design the most appropriate strategies for specific crops, varieties and environments, striking a balance between water productivity, environmental safeguard and economic feasibility.

Crop Water Stress in Dynamic Crop Models
Complex crop models simulating growth and development of herbaceous and vegetables crops (e.g., maize, wheat, soybean, tomato, etc.) and WS were selected on the basis of the following characteristics. They had to be: (i) freely available, (ii) with an easily accessible documentation, (iii) including water balance and (iv) with options for setting irrigation scheduling rules. The selected models are: APSIM, AQUACROP, CROPSYST, DAISY, DSSAT, EPIC, RZWQM2, STICS. A list of the models, websites and main developer and maintainer is provided in Table ??. The formalism of these models will be reviewed in the next subsections by grouping the models according to common features on different aspects of WS and CWR topics. The models will be categorized on the basis of: (i) the functions used to calculate biomass growth and development; (ii) the main variable used to calculate WS coefficients and (iii) the influence of WS on vegetative and reproductive growth. Tracking back the rationale behind crop growth and development can provide opportunities for adjustments of DSSs and RS-related tool used for CWR estimations. As shown in the previous Section ??, modeling concepts based on plant processes and soil-water-plant dynamics have often been used to improve existing techniques. Some of these concepts and dynamics are represented in similar ways among different models and different crop modules within the same model. Overall, there is a broad consensus around relevant processes affected by water stress and stress thresholds.

Water-and Radiation-Driven Crop Growth and Crop Development
A definition of the terms "growth" and "development" is paramount to approach the description of plant-related processes. Following Basso and Ritchie [? ], we will consider crop growth as the plant biomass accumulation and crop development as the progression of the plant through growth stages (phenological development). The distinction between these two concepts is of primary importance for the adequate simulation of crop dynamics in different environments.
The selected models can be categorized in "water-driven" and "radiation-driven" models, depending on the mechanism chosen for the simulation of biomass growth.

Water-Driven Models
AQUACROP is a water-driven model developed by FAO. The concepts behind the model are extensively detailed in FAO-66 paper [? ], that updates the FAO-33 approach [? ] on water productivity and estimates ETa following the FAO-56 method [? ]. Model concepts, usage and parameterization are widely described a series of paper from 2009 [? ? ? ]. AQUACROP is water-driven, in that the biomass accumulation is directly related to crop transpiration and evapotranspiration, via a crop-dependent water-biomass productivity coefficient (Figure ??). Biomass water productivity has shown a remarkable conservative behavior, remaining nearly constant within species when normalized for the vapor pressure deficit (VPD), as demonstrated in the well-known studies by Tanner and Sinclair [? ] and Sinclair et al. [? ]. The approach considers that the biomass accumulated by plants in the same environment is directly proportional to the accumulated water volume transpired, being the stomata both the entry point of CO 2 for photosynthesis and the exit point of the absorbed water. The slope of the relationship constitutes the water use efficiency (WUE), i.e., aboveground biomass production per each unit of water transpired. Normalized WUE is considered to be fairly constant for different crops and variates, with major differences between C3 and C4 plants.   [88], that updates the FAO-33 approach [89] on water productivity and estimates ETa following the FAO-56 method [5]. Model concepts, usage and parameterization are widely described a series of paper from 2009 [90][91][92]. AQUACROP is waterdriven, in that the biomass accumulation is directly related to crop transpiration and evapotranspiration, via a crop-dependent water-biomass productivity coefficient (Figure 2). Biomass water productivity has shown a remarkable conservative behavior, remaining nearly constant within species when normalized for the vapor pressure deficit (VPD), as demonstrated in the well-known studies by Tanner and Sinclair [93] and Sinclair et al. [94]. The approach considers that the biomass accumulated by plants in the same environment is directly proportional to the accumulated water volume transpired, being the stomata both the entry point of CO2 for photosynthesis and the exit point of the absorbed water. The slope of the relationship constitutes the water use efficiency (WUE), i.e., aboveground biomass production per each unit of water transpired. Normalized WUE is considered to be fairly constant for different crops and variates, with major differences between C3 and C4 plants. Figure 2 provides an exemplification of this approach as presented in AQUACROP, where VPD is substituted by ET0 to account for environmental conditions. This modeling concept can give comparable results to that of other more complicated crop models in most situations [95,96], even though some authors questioned its validity under conditions of adequate water supply [97] where agronomic management can play an important role in increasing crop yield, the transpiration efficiency being equal [98]. AQUACROP is still maintained and currently at his 6th version and an open-source Matlab version has recently been released to permit wider use for research purposes by granting access to the source code [99]. This allows easier implementation of inverse modeling techniques and integrating with RS data [100].

Radiation-Driven Models
DAISY, DSSAT, EPIC, STICS and RZWQM2 are radiation-driven models, i.e., the biomass accumulation is directly (linearly) related to the photosynthetically active radiation (PAR, the radiation between 400 and 700 nm) absorbed or intercepted by canopy (APAR and IPAR, respectively). The slope of this relationship is the radiation use efficiency (RUE, also called light use  efficiency, LUE), as detailed in Rosati et al. [101]. The approach used in radiation-driven models is based on early work by Monteith [102,103] (Figure 3). DAISY is described in Abrahamsen and Hansen [104], Hansen [105] and Hansen et al. [106]. The DSSAT (Decision Support System for Agrotechnology Transfer) has been developed and maintained for a long time, currently collecting together two main crop modules: CERES for cereals and CROPGRO, originally for legumes [107,108]. EPIC was initially developed to study the relationship between soil erosion and crop productivity in the US, including simulations of economic related processes [109,110]. STICS is a model with a generic crop module able to adapt to different plants [111][112][113][114]. RZWQM2 (root zone water quality model 2, described in [115]), is an improved version of RZWQM [116] that aims at simulating mainly the fluxes of water, nutrients and chemicals in the root zone, implementing a DSSAT-derived crop module for the simulation of crop growth and development. .

Hybrid Models
APSIM [117,118] and CROPSYST [119,120] use a hybrid approach. These models calculate both radiation-and water-driven potential biomass accumulation. Following the law of limiting factors, the lower of the two is kept as the actual biomass accumulation. APSIM 7.1 is the current version to which we will refer, though the research team is working on an improved version called "APSIM Next Gen". In APSIM, when transpiration is not limiting, soil water demand (transpiration) is calculated backward from the biomass accumulation using the transpiration efficiency normalized by VPD (i.e., using the Tanner-Sinclair relationship in the opposite direction respect to AQUACROP).

Crop Development
The basic concept used to model crop development is common among all the selected models. Crop development is modeled as a function of thermal time (see the definition in Trudgill et al. [121]), which is expressed in growing degree days (GDD). This is calculated daily, as the difference between the actual temperature and a crop-specific base temperature under which development stops. AQUACROP, for example, simulates crop development using thermal time only, while other models like APSIM include also modifying factors to account for the effects of photoperiod and vernalization. In general, though, further additions to the basic thermal time concept vary between crop modules within the same model (e.g., in DSSAT) and among crop models. The development rate can be accelerated by WS, and WS may have peculiar effects on plant processes in specific growth stages [122], e.g., during flowering and grain filling.  ] use a hybrid approach. These models calculate both radiationand water-driven potential biomass accumulation. Following the law of limiting factors, the lower of the two is kept as the actual biomass accumulation. APSIM 7.1 is the current version to which we will refer, though the research team is working on an improved version called "APSIM Next Gen". In APSIM, when transpiration is not limiting, soil water demand (transpiration) is calculated backward from the biomass accumulation using the transpiration efficiency normalized by VPD (i.e., using the Tanner-Sinclair relationship in the opposite direction respect to AQUACROP).

Crop Development
The basic concept used to model crop development is common among all the selected models. Crop development is modeled as a function of thermal time (see the definition in Trudgill et al. [? ]), which is expressed in growing degree days (GDD). This is calculated daily, as the difference between the actual temperature and a crop-specific base temperature under which development stops. AQUACROP, for example, simulates crop development using thermal time only, while other models like APSIM include also modifying factors to account for the effects of photoperiod and vernalization. In general, though, further additions to the basic thermal time concept vary between crop modules within the same model (e.g., in DSSAT) and among crop models. The development rate can be accelerated by WS, and WS may have peculiar effects on plant processes in specific growth stages [? ], e.g., during flowering and grain filling.
Key state variables that express the canopy development in the selected models are LAI for all the models except AQUACROP. The latter uses the relative canopy cover (called CC in the model, equivalent to the FCOVER of RS), i.e., the percentage of ground surface covered by the canopy, due to its more DSS-oriented nature, since canopy cover is easy to estimate even just by visual estimates or by downward-looking photographs of the crop. This distinction is particularly relevant when the models are integrated with RS since FCOVER is usually easier to estimate than LAI and is not subject to saturation at high values. For example, Silvestro et al. [? ] assimilated FCOVER data estimated from HJ1A/B satellite into the AQUACROP model to provide wheat yield estimations at field and district level.

Soil Water Content and Transpiration Deficit as Inputs of Water Stress Functions
All the selected models can estimate ET0 using the Penman-Monteith formula, except AQUACROP that requires ET0 as an input variable. Within the different models, other less data hungry options are available (e.g., Priestley-Taylor [? ]), as well as more complicated ones (e.g., Shuttleworth-Wallace [? ]). All the models separate the evaporative fraction of the soil from the transpirative fraction of the plants. The movement of water in soil is simulated either by using solutions of the Richards equation (CROPSYST, DAISY, RZWQM2) or by the simpler cascade approach (AQUACROP, CROPSYST, DSSAT, EPIC, APSIM). The basis of the most frequently used method for separating evaporation and transpiration and for modeling downward flow and root water uptake in a layered soil (cascade approach) are described by Ritchie [? ? ]. As for the upward flow (ET-driven), separation of soil evaporation from plant transpiration is dependent on the canopy cover (e.g., through LAI) and is particularly relevant in earlier stages of crop growth, while when the canopy completely covers the ground transpiration contributes the most to the overall ET. As for the downward flow (gravity driven), in the cascade approach excess water (above field capacity) drains to deeper soil layer with a speed that depends on soil hydraulic conductivity.
WS is calculated on the basis of either (i) soil water content or (ii) transpiration deficit.

Soil Water Content
Soil water content for WS calculations is related to the total available water (TAW), where TAW is the difference between field capacity (FC) and plant wilting point (PWP). TAW is also known as plant available water (PAW), root available water or, simply, available water. The input state variable of WS functions can be expressed in relative or absolute values of TAW (e.g., ????). Other common ways to indicate the soil water limits of TAW are DUL (drain upper limit), that corresponds to field capacity, and LL (lower limit of plant water availability), that corresponds to the wilting point [? ]. This alternative nomenclature is common in models derived or inspired by DSSAT. The soil water limits are agronomic concepts used by most crop models to describe in a simplified way relevant soil water potential. However, Ritchie [? ? ] warned against the use of water potential to define field water limits to model water balance, especially when the potential is measured with laboratory pressure extraction equipment. The main concerns regard the disturbances related to the extraction of soil cores and the uncertainty around the water potential value of field capacity. These limits should therefore be measured in the field, or through pedotransfer functions [? ], that may be already included in the model code and operate automatically.  [88]). The state variable used to calculate the coefficients is the total available water (TAW) relative depletion ratio. The example refers to maize water stress (WS) calculated at ET0 of 5 mm. The first stress threshold (leaf expansion affected) is met at 14% of TAW depletion. Transpiration and senescence processes begin to be affected both at 69% TAW depletion. Stress threshold for senescence is usually met at equal or higher depletion ratio than for transpirative stress, depending on the crop.  ). The state variable used to calculate the coefficients is the total available water (TAW) relative depletion ratio. The example refers to maize water stress (WS) calculated at ET0 of 5 mm. The first stress threshold (leaf expansion affected) is met at 14% of TAW depletion. Transpiration and senescence processes begin to be affected both at 69% TAW depletion. Stress threshold for senescence is usually met at equal or higher depletion ratio than for transpirative stress, depending on the crop.   [88]). The state variable used to calculate the coefficients is the total available water (TAW) relative depletion ratio. The example refers to maize water stress (WS) calculated at ET0 of 5 mm. The first stress threshold (leaf expansion affected) is met at 14% of TAW depletion. Transpiration and senescence processes begin to be affected both at 69% TAW depletion. Stress threshold for senescence is usually met at equal or higher depletion ratio than for transpirative stress, depending on the crop.

Transpiration Deficit
Transpiration deficit is expressed as the ratio of potential root water uptake (RWUp) to potential transpiration (Tp), or as a supply/demand ratio (which is the same thing with a different name). Examples are shown in Figures ?? and ??b. RWUp is usually calculated from root density, soil resistance and soil water conductivity [? ], and constitutes the potential supply of water from the soil. Tp, obtained by separating the evaporative and transpirative fractions of ETp, represents the potential atmospheric demand of water. When the water supply is greater than the demand (RWUp/Tp > 1), there is no limitation in the transpiration process and the actual transpiration (Ta) is equal to Tp. When the supply is smaller than the demand (RWUp/Tp < 1), the potential water supply is absorbed and transpired completely (RWUp = Ta) and there is a transpirative stress (Ta < Tp).
Remote Sens. 2020, 12, x FOR PEER REVIEW 18 of 34 Figure 6. WS coefficients for vegetative growth initially developed for DSSAT (redrawn from Ritchie [51]). The stress is calculated on the ratio between potential root water uptake and potential transpiration. Transpiration and photosynthesis are affected when RWUp/Tp drops below 1, leaf expansion is affected when RWUp/Tp drops below 1.5, to account for the greater sensitivity of this process. ). The stress is calculated on the ratio between potential root water uptake and potential transpiration. Transpiration and photosynthesis are affected when RWUp/Tp drops below 1, leaf expansion is affected when RWUp/Tp drops below 1.5, to account for the greater sensitivity of this process.

WS Functions Used in Crop Models
To calculate WS functions, AQUACROP uses TAW depletion ratio (Figure ??) and STICS the absolute TAW value (Figure ??). CROPSYST, DAISY, DSSAT, EPIC and RZWQM2 use the transpiration deficit (an example taken from DSSAT is provided in Figure ??). APSIM uses either TAW ratio (Figure ??a) or transpiration deficit (Figure ??b) depending on the process involved. The functions calculate a water stress factor/coefficient/modifier that acts to modify the relative process rate of several target processes. The visual representation of the stress functions of these 4 models covers essentially all the approaches used in the selected models. The processes targeted by WS are discussed in Section ??.

Plant Processes Targeted by Water Stress
For the description of the processes targeted by WS, we will refer from now on to a generic cereal with determinate growth, for simplicity. In particular, for DSSAT we will refer to the DSSAT-CERES crop module and for APSIM to the "Plant" crop module, that replaces previous modules covering the relevant crops. Several functions are used to simulate the WS effects on different plant processes, but the number and type of the effects can differ between crops, especially between determinate and indeterminate crops. Almost all the selected models and crop modules, share similar concepts and formulations of the WS effects on (i) leaf expansion, (ii) transpiration/photosynthesis and (iii) senescence. These functions often come in the models' documentation or in reference papers with a visual representation of the stress functions on the basis of the state variables used for the calculations (as shown in Section ??). Other stress effects are sometimes simulated, but with a lower degree of agreement among different crops and models. These include modifications of the harvest index, acceleration/deceleration of phenology development, shift in photosynthates partitioning to plants' organs and restrictions to root growth. Table ?? sums up the relevant information on crop models' features and processes affected by WS, discussed in detail in the following Sections. Figure 6. WS coefficients for vegetative growth initially developed for DSSAT (redrawn from Ritchie [51]). The stress is calculated on the ratio between potential root water uptake and potential transpiration. Transpiration and photosynthesis are affected when RWUp/Tp drops below 1, leaf expansion is affected when RWUp/Tp drops below 1.5, to account for the greater sensitivity of this process.  Table 2. Summary of the relevant crop models' features discussed in the paper. R: radiation-driven; W: water driven; R-W: radiation-water hybrid approach (the minimum is kept); LAI: leaf area index; CC: canopy cover (usually called FCOVER in RS); avW: available water (above PWP); Tdef: transpiration deficit (Ta/Tp).
Crosses indicate that the model simulates the direct effects of water stress on the specific plant growth or development process.

WS Effect on Leaf Expansion
In the progress towards lower soil water content or higher transpiration deficit, the leaf expansion stress is always the first to be encountered, being the most sensitive process. Water stress factors target the equations that regulates canopy/leaf area expansion. From the physiological point of view, WS effect can be thought as a reduction in plant cells turgor and is rarely considered in irrigation planning, especially with RS applications (Section ??). However, some authors pointed out that prolonged mild water stress during the canopy growth period, even if not yet severe enough to inhibit transpiration and photosynthesis, might significantly reduce plant productive capabilities through the steadily reduction of the photosynthetic surface [? ? ]. Models that consider this effect of WS are: APSIM (Figure ??b), AQUACROP (Figure ??), CROPSYST (where leaf expansion completely stops when the Ta/Tp ratio drops below 0.9), DSSAT (Figure ??), EPIC (where the WS can reduce LAI expansion) and STICS (Figure ??). In all these cases leaf expansion is considered the process most sensitive to WS. DAISY and RZWQM2 do not simulate direct effect of WS on leaf expansion.

WS Effect on Transpiration/Photosynthesis
The WS effect on stomatal conductance affects biomass accumulation through reduced transpiration or reduced photosynthesis (due to reduced transpiration), depending on the formulation of each model. From the physiological point of view, water deficit causes stomata closure, reducing transpiration and CO 2 assimilation and, in turn, photosynthesis and biomass accumulation. Irrigation scheduling usually aims at avoiding completely stomatal closure (except for deficit irrigation), as highlighted by Jovanovic et al. [? ]. This stress is so important and limiting that it is simulated by all the selected models. APSIM (Figure ??b), CROPSYST, DAISY, DSSAT (Figure ??), EPIC and RZWQM2 use the Ta/Tp ratio to reduce the biomass accumulation (in RZWQM2 the relationship between WS coefficient and transpiration deficit is nonlinear [? ]. In DAISY, the biomass reduction is the only effect of WS considered. The reduction is based on soil water content in AQUACROP (Figure ??) and STICS (Figure ??), where the stress coefficient reduces transpiration and photosynthesis, respectively.

WS Effect on Senescence
WS can increase the senescence rate and also cause an earlier onset of senescence. The senescence process is affected by WS in most of the models and is the least sensitive process in respect to leaf expansion and transpiration reduction (see ????). Senescence is affected either through a direct increase in senescence rate or through a decrease in leaf area duration (LAD) and adjustment of the LAI curve. WS effect on senescence is modeled in APSIM, using WS factor for photosynthesis (Figure ??b), AQUACROP (Figure ??), CROPSYST, where senescence is determined by LAD and STICS (Figure ??). DAISY, DSSAT-CERES, EPIC and RZWQM2 do not directly model this effect.

WS Effect on Other Processes
WS can also influence harvest index (HI), phenology, partitioning and root growth, but these effects are not modeled consistently in all the models. HI is modified by WS in APSIM, AQUACROP, CROPSYST and EPIC, mainly during flowering and grain filling period. In APSIM, WS reduces HI depending on stress severity and a sensitivity factor. AQUACROP derives a HI modifier from the stress coefficients of leaf expansion (that increases HI) and of transpiration (that instead decreases HI). CROPSYST reduces HI using the lowest of water and N stress, with different sensitivity depending on stage. In EPIC, WS acts on HI depending on the stress range: with low WS (coefficient > 0.9) the HI increases, with higher WS (coefficient between 0.3 and 0.9) the HI decrease [? ]. In general, it is suggested that moderate stress can be beneficial (increase in HI, in some models), while more severe WS is always detrimental for production (decrease in HI).
WS effect on development (phenology) is modeled in APSIM (Figure ??a), CROPSYST, STICS and RZWQM2. Development is slowed down by WS in APSIM from emergence to flowering, while in CROPSYST and STICS it is hastened (due to an increase in crop temperature, that governs development) [? ]. In RZWQM2 moderate WS (Ta/Tp > 0.6) can hasten development, while severe stress (Ta/Tp < 0.6) slows it down.
A shift in partitioning can also be caused by WS in some of those models that simulate photosynthates allocation in different parts of the plant, but there is no consistent way of simulating it among the models. DSSAT-CERES shift allocation towards roots during WS, APSIM towards grain during grain filling period if there are not enough photosynthates, EPIC simulates root growth on the basis of RWU, RZWQM2 reduces partitioning to shoots under WS. Root growth can also be restricted by low soil water content, as in APSIM, DSSAT and STICS.

Comparison of Dynamic Crop Models
To choose the most suitable approach to model crop growth, CWR and WS is up to the modeler, depending on the purpose of the research or application, the desirable level of model complexity, data availability and user friendliness. The question of which model is "the best" has, therefore, virtually no univocal answer. Choosing a more complex and detailed model is not always a guarantee of getting better results than using a simpler, well calibrated model. However, research experience accumulated in scientific literature, comparing the performances of different models to field data in the presence of water limitations, can help to understand specific features of each model. For this reason, a few examples of model comparison will be discussed. Only the models previously discussed and reported in Table ?? will be considered.
Kloss et al.
[? ] compared among others DAISY and APSIM performances in simulating crop growth under deficit irrigation using field data from India, Malawi, France and Oman. The authors noticed that the yield reduction simulated by DAISY in cases of severe water stress was smaller and further from reality than for APSIM. The authors suggested that in the cases of severe water stress, especially when after a stress period full irrigation was restored, DAISY was unable to simulate irreversible damage to the crop. Even if photosynthesis was reduced during drought stress, when water availability was fully restored, both potential transpiration and photosynthesis where reached again at the same levels than before. On the other hand, APSIM also models WS effects on leaf expansion, phenology and partitioning, as previously shown in Section ??. However, the authors also stressed that the severe and alternating stress conditions where DAISY simulation fell short are unlikely to be met in non-experimental conditions.
Todorovic et al.
[? ] compared AQUACROP and CROPSYST performances in simulating sunflower growth under deficit irrigation over two cropping seasons in Southern Italy, i.e., in a Mediterranean climate. The experimental treatments were designed to provide deficit irrigation during the whole season, but also to differentiate it before and after flowering. AQUACROP performed similarly to CROPSYST, but required fewer inputs. CROPSYST simulated better the WUE under severe water stress. The authors suggested that the physiological submodel that simulates crop growth and partitioning in AQUACROP is more simplified than in CROPSYST, thus justifying the lower performances under severe water stress. Indeed, as explained in Section ??, CROPSYST simulates WS effect on phenology (hastening development under WS), while AQUACROP does not, and the effect on HI is also more detailed. However, the authors stressed that both models performed well under moderate stress conditions, and that their results might have been influenced by the manual calibration and the short duration of the experiment, suggesting that automatic calibration over multiple years may produce better results.
Abi Saab et al.
[? ] also compared AQUACROP and CROPSYST. The field experiment was on barley cultivated over 3 years, with treatments including rainfed conditions, full irrigation (maintaining full transpiration) and 50% of full irrigation. The calibration year was alternated to check the effect of calibration on the performance of the models on the validation datasets. In general, AQUACROP showed better performances, but the model was recalibrated for each non-optimal N supply conditions, while CROPSYST was calibrated once for each year (under optimal N supply conditions). This was due to the lack of a N balance module in AQUACROP. The authors suggested that CROPSYST should be a better choice in when nutrient and water stress are balanced, while when the yield reduction is due to water stress alone, AQUACROP may be preferable. Rötter et al. [? ] tested several models over 14 years of data for spring barley over multiple locations and various weather conditions. Among the models, CROPSYST, DAISY, DSSAT-CERES and STICS were included. The models were given to groups of users, to perform a "blind test": limited data was available for calibration (only phenological observations) and the modelers were allowed to calibrate only phenological parameters. Overall, yield estimates were unsatisfactory, with measurement error > 10-15 % respect to the observed yield. The uncertainty in the estimation was similar to that found by Palosuo et al. [? ] for wheat, even though the initial hypothesis of the authors was that the uncertainty was higher since barley had been less studied by the modeling community, thus having more imprecise formalism and parameters. The authors proposed that one of the reasons for the lower-than-expected error of barley simulations could have been the shorter cycle of spring barley and the absence of vernalization processes to be simulated, that instead were thought to have played a relevant role in increasing the uncertainty of winter wheat predictions in Palosuo et al. [? ]. Overall, the phenology showed considerable discrepancies between simulated and observed data. The authors suggested that the complexity of effects of water and nutrient limitation in delaying or hastening the phenological development was probably not caught in sufficient detail by all the selected models. At the same time, though, the models with more accurate estimations of phenology, they did not necessarily predict yield more accurately. As for grain yield, DAISY and CROPSYST systematically underestimated yield, while the greatest variance in model residuals (inaccuracy of the estimation) was found in DSSAT-CERES, while DAISY showed the strongest correlation between observed and simulated yield (lowest uncertainty).
Kimball et al.
[? ] compared ET estimates among 29 models or model versions given to several modelers' groups to simulate 8 years of maize monocropping in Iowa. Among the models, CROPSYST, RZWQM2, and various versions of APSIM, DSSAT-CERES and STICS were used. Differences between versions regarded mainly the ET estimation process (e.g., using different equations). The study was conducted in 4 phases: (i) a "blind phase" where only key input data on soil, weather and management and anthesis and maturity dates were provided; (ii) a "LAI phase" where LAI observations time series were provided, to see if improving simulations of leaf expansion also improved ET predictions; (iii) a phase where also observations of biomass, yield, ET and soil moisture of a typical year for rainfall was provided and (iv) a phase where all the observations were provided for all years to allow for detailed calibration. Overall, supplying LAI for all the years led only to slight improvements of ET estimates. With the additional information of phase 4, APSIM versions that used the Priestley-Taylor approach to estimate ET were the models that provided better evaporation estimates during early phases of crop growth. This highlighted the importance of providing information on crop variables to improve ET estimations and of providing good early-stage predictions of soil evaporation and leaf area development to obtain accurate measurements of seasonal ET. ETa during mid-season (when Ta contributes to a greater extent) was usually better predicted through the FAO-56 Penman-Monteith approach or the Shuttleworth-Wallace approach (RZWQM2 and a version of STICS). The authors noted little overlap among the models identified as best for the early-stage, evaporation-dominated period and for the mid-season, T-dominated period, suggesting that models should be improved by combining the best approaches to simulate evaporation with those to simulate T. Considering the performances on estimating both ET and agronomic variables (LAI, biomass and yield) after full calibration (phase 4), the median of all the models ranked first, followed by versions of APSIM and DSSAT-CERES (among the models selected for this review).

Integration of RS Data and Crop Models
DSS platforms to support irrigation scheduling often already include simple statistical or mechanistic models. In this framework, the application to DSSs of the concepts highlighted for dynamic models consists of a revision of the formalism of the DSS, or the addition of further simulation modules, e.g., water balance module or crop growth and development module, when this is considered appropriate. On the other hand, RS imagery requires subsequent elaborations to be used to estimate crop growth and development, WS and CWR, ranging between different degrees of complexity and directness. From the operational point of view, one of the most complicated and indirect, yet promising procedure, is the integration of RS data with dynamic crop models. RS data can be considered a natural partner of crop models, as testified by historical developments highlighted by Pinter et al. [? ] and Jones et al. [? ]. Kasamplis et al. [? ] stressed that one of the major limitations of crop models at field scale is the lack of spatial information of the conditions within the fields, but RS can provide this information to improve models' predictions. RS can be used to estimate state variables of crop models like FCOVER, LAI and fAPAR. One of the methods for integrating RS data with crop models that has received increased attention is data assimilation [? ]. Three different strategies can be used to integrate biophysical variables derived from RS observations into crop models: calibration, forcing and updating [? ]. As a rule of thumb, calibration adjusts the model parameters to minimize differences between historical observed and simulated data without considering the uncertainty contained in the observed data, forcing methods replace simulated with observed data having the same time-step (daily in most crop models) and the updating methods usually provide a real-time adjustment of model parameters considering the uncertainty related to both simulated and observed data [? ].
Crop models integrating RS data can be used to provide a continuous characterization of the agroecosystem over time and forecasts of its future state, which is something that RS alone struggles to do [? ]. In addition, modeling approaches permit to separate crop stresses and processes, while RS measures the actual growth conditions that includes the combined effects of growth reduction factors (e.g., water and nutrient stress together) [? ]. Data assimilation can also integrate simultaneously several types of data into the same modeling scheme, e.g., multisource RS data, as well as data from airborne, Unmanned Aerial Vehicles (UAVs) and ground proximal sensors platforms. It also allows the use of near real-time meteorological data and forecasts.
In this context, users can select models based on their needs and on the best approach that fits available data from RS or other sources. The modeled effects of WS on specific plant processes of interest can also be considered. For example, correct quantification of water stress effects on crop biomass and yield are paramount in cases of deficit irrigation, where the assessment of a stress/unstressed condition is not enough to maximize the efficiency of the irrigation water use, rather than maximizing yield. Or else, a representation of water stress effects on leaf expansion might be useful in environments characterized by prolonged moderate stress during early stages of crop growth. Equations and approaches from crop models presented in this work can provide examples to implement further RS and DSS-based estimations of WS and CWR. At the same time, an overview that can guide the users through the choice of the model best suited for their purposes is provided.

Conclusions
We have provided an overview of recent approaches based on remote sensing and decision support systems conceived for aiding in the management of crop water requirements and water stress, highlighting relevant steps undertaken over the years to improve these tools. A general trend of increasing complexity can be noticed, with the aim of introducing relevant features or adjustment to better fulfill irrigation management support. ET estimation from remote sensing imagery and advices from decision support systems offer great opportunities for scheduling irrigation, but are limited by an incomplete representation of the soil-plant-atmosphere dynamics. Dynamic crop models offer a better understanding of the processes related to water stress, but are usually more suitable for research purposes than for direct applications by farmers and agricultural technicians. However, freely available dynamic crop models with accessible documentation (such as APSIM, AQUACROP, CROPSYST, DAISY, DSSAT, EPIC, RZWQM2 and STICS) stand as a source of condensed knowledge on the processes governing plant water stress and stress effects on plant growth and development. To some degree, these models share concepts and formalism that can be classified together. The models simulate biomass accumulation from crop transpiration, intercepted radiation or using both the approaches and selecting the most limiting. Stress functions are based on soil water content or deficit transpiration. The main processes affected by water stress in similar ways among the models include: leaf expansion (the most sensitive), transpiration/photosynthesis (usually determining irrigation thresholds, since it heavily affects biomass accumulation) and senescence. The information provided, along with the references to relevant documents and papers regarding these concepts, may constitute a starting point for researchers seeking to improve RS and DSS approaches in the field of water stress and crop water requirements estimations. The best of these worlds can be combined and integrated in several ways, keeping in mind that the choice to be made should comply with the final goal to be pursued. Who is the end user? Which degree of operational complexity should be allowed? What is the temporal perspective of interest: predictions, real-time information, or both? Is a qualitative indication of water status enough, or is a quantification of water stress effects on plant processes more appropriate? Does irrigation timing and amount need to be accurately assessed? All are questions worth asking when choosing the most appropriate approach.

Acknowledgments:
The authors acknowledge the provision of documentation and in some cases the source code of the crop models mentioned in the manuscript by the model developers. The authors are thankful to the four anonymous reviewers who contributed with their careful work to improving the focus and quality of the paper.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1. Glossary of abbreviations and recurrent terms used in the paper.

RS Remote sensing
Technical and scientific discipline to acquire information from an object without making physical contact with it.

RUE Radiation use efficiency
Slope of the relationship between dry biomass accumulation and absorbed or intercepted PAR (dry matter produced per each unit of radiation).

RWUp
Root water uptake (potential) Amount of water supplied by the soil that can be absorbed by roots.

SAVI Soil adjusted vegetation index
Reflectance-based vegetation index originally developed by Huete [? ], includes adjustment for canopy background (soil).

SEB (approach)
Surface energy balance Model/approach to estimate energy fluxes between Earth's surface and atmosphere.

Ta
Actual transpiration Effective transpiration of a specific crop when stresses and restrictions are accounted for.

TAW Total available water
Amount of water stored in a specified soil layer between FC (also called DUL) and PWP (also called LL). It is the water that can be uptaken by the plant.

Tp Potential transpiration
Potential transpiration of a specific crop when transpiration and photosynthesis are not reduced by stresses and restrictions.

VI Vegetation Index
Index obtained from crop/cropping field reflectance that can be related to crop physiological parameters.

VPD Vapor pressure deficit
Difference between the amount of moisture present in the air and the amount of moisture air can hold when saturated. Calculated from temperature and humidity.

WS crop Water Stress
Suboptimal condition of plant water status, when the plant is not able to fully sustain all its processes due to water shortage. WS affects several plant processes, with different sensitivity.