An Integrative Information Aqueduct to Close the Gaps between Satellite Observation of Water Cycle and Local Sustainable Management of Water Resources

: The past decades have seen rapid advancements in space-based monitoring of essential water cycle variables, providing products related to precipitation, evapotranspiration, and soil moisture, often at tens of kilometer scales. Whilst these data e ﬀ ectively characterize water cycle variability at regional to global scales, they are less suitable for sustainable management of local water resources, which needs detailed information to represent the spatial heterogeneity of soil and vegetation. The following questions are critical to e ﬀ ectively exploit information from remotely sensed and in situ Earth observations (EOs): How to downscale the global water cycle products to the local scale using multiple sources and scales of EO data? How to explore and apply the downscaled information at the management level for a better understanding of soil-water-vegetation-energy processes? How can such ﬁne-scale information be used to improve the management of soil and water resources? An integrative information ﬂow (i.e., iAqueduct theoretical framework) is developed to close the gaps between satellite water cycle products and local information necessary for sustainable management of water resources. The integrated iAqueduct framework aims to address the abovementioned scientiﬁc questions by combining medium-resolution (10 m–1 km) Copernicus satellite data with high-resolution (cm) unmanned aerial system (UAS) data, in situ observations, analytical- and physical-based models, as well as big-data analytics with machine learning algorithms. This paper provides a general overview of the iAqueduct theoretical framework and introduces some preliminary results.


Stakeholder Requirements and Potential Knowledge Gaps
To address the panoply of scaling issues described above and the multitude of user requirements (see Appendix A), iAqueduct will deploy six field observatories (five across Europe and one in Israel) for intensive studies: (i) The Twente site in the Netherlands, which serves as a core international site for SMOS/SMAP cal/val activities; (ii) the Zala catchment in Hungary, which has served as a study site for analyzing the performance of the European pedotransfer functions in deriving soil hydraulic maps; (iii) the Sde Yoav field in Israel, which is a well-documented site for soil investigation in Israel; (iv) the Alento River Hydrological Observatory in Italy, where intensive soil moisture and hydrological observations and modelling have been conducted; (v) the Corleto area in Italy, which has been used for detailed UAS research; and (vi) the Barranco del Carraixet area in Spain, a study site for dry land water management and interaction with stakeholders. The detailed description of each observatory and associated stakeholder requirements are given in Appendix A.
The analysis of stakeholder requirement identifies the need to support and facilitate the establishment of water management policies; addressing rapid climatic changes by involving researchers, water management authorities, companies, and farmers; and strongly supporting progress based on previous findings in each site/catchment and by dealing with local needs. In other words, the stakeholder calls for the translation of science and knowledge (about the response of hydrological cycle and water resources to climate change), into marketable tools, services, and/or products for the sustainable management of water resources. This is actually demanding the establishment of a science-policy-business-society interface to allow for continuous dialogues and interactions across different scales and levels, influencing stakeholders towards desirable behaviors [3].
To address the aforementioned stakeholder requirements, it requires the development of beyond the state-of-the-art approaches to derive local field-scale soil, vegetation, and water states and information (e.g., mainly precipitation, evapotranspiration, and profile soil moisture), using satellite, UAS, in situ observations, as well as modelling and big data analytics tools for water management under climate change. It is well-known that space-based EOs are highly effective to characterize water cycle variation at the regional to global scale (see Section 1.1) but are less so at the local and field scale Water 2020, 12, 1495 5 of 36 to provide more detailed information for the sustainable management of water resources. To this aspect, it is important to consider the heterogeneous characteristics of the soil and vegetation at these finer scales and to effectively bridge existing knowledge at different scales. We thus need to answer the following questions: -How to downscale the global water cycle products to the local scale using multiple sources and scales of EO data? -How to explore and apply the downscaled information at the management level for a better understanding of water-energy-soil-vegetation processes? -How can such fine-scaled information be used to improve the management of soil and water resources?
In the next section, we present the iAqueduct theoretical framework to address the above questions. Figure 1 describes the iAqueduct framework of methodologies and approaches. It includes six closely connected working blocks (WBs). WB1 deals with the scaling from global satellite water cycle products to field-scale water states, which includes both the surface and profile information on soil water states. Specifically, WB1 will advance the space-time characterization of soil moisture and evapotranspiration processes through the combined use of field, UAS, and satellite observations. In particular, the combined use of high-resolution soil characteristics and satellite data will increase our capabilities to describe soil moisture and evapotranspiration processes with high-level detail.

iAqueduct Framework
Water 2020, 12, x FOR PEER REVIEW 5 of 37 -How to downscale the global water cycle products to the local scale using multiple sources and scales of EO data? -How to explore and apply the downscaled information at the management level for a better understanding of water-energy-soil-vegetation processes? -How can such fine-scaled information be used to improve the management of soil and water resources?
In the next section, we present the iAqueduct theoretical framework to address the above questions. Figure 1 describes the iAqueduct framework of methodologies and approaches. It includes six closely connected working blocks (WBs). WB1 deals with the scaling from global satellite water cycle products to field-scale water states, which includes both the surface and profile information on soil water states. Specifically, WB1 will advance the space-time characterization of soil moisture and evapotranspiration processes through the combined use of field, UAS, and satellite observations. In particular, the combined use of high-resolution soil characteristics and satellite data will increase our capabilities to describe soil moisture and evapotranspiration processes with high-level detail. It has been demonstrated that soil hydraulic and thermal properties (SHP/STP) play a critical role in determining soil water and heat flow at field/plot scales, while such information is rarely available at such detailed scales [60,61]. WB2 will apply pedotransfer functions to derive local field specific SHP/STP properties for the modelling of soil water and heat dynamics at field-scale precision. It will bridge soil spectral information that can be obtained at a high resolution by satellites and UAS and the needed soil properties that are traditionally obtained at limited locations by in situ sample collections.

iAqueduct Framework
Using the information obtained from the two previous WBs, WB3 attempts to retrieve field-and grid-specific relationship functions between soil properties, soil moisture, and evapotranspiration. Such a relationship function is expected to advance the current hydrological modelling concepts, in which the actual evapotranspiration is parameterized on the availability of soil moisture using untested (linear) assumptions. Field-specific functions will be developed on the basis of downscaled It has been demonstrated that soil hydraulic and thermal properties (SHP/STP) play a critical role in determining soil water and heat flow at field/plot scales, while such information is rarely available at such detailed scales [60,61]. WB2 will apply pedotransfer functions to derive local field specific SHP/STP properties for the modelling of soil water and heat dynamics at field-scale precision. It will bridge soil spectral information that can be obtained at a high resolution by satellites and UAS and the needed soil properties that are traditionally obtained at limited locations by in situ sample collections.
Using the information obtained from the two previous WBs, WB3 attempts to retrieve field-and grid-specific relationship functions between soil properties, soil moisture, and evapotranspiration. Such a relationship function is expected to advance the current hydrological modelling concepts, in which the actual evapotranspiration is parameterized on the availability of soil moisture using untested (linear) assumptions. Field-specific functions will be developed on the basis of downscaled satellite observation of soil moisture and evapotranspiration, and their combined analysis with in situ measurements and UAS observations. WB4 is expected to advance ecohydrological modelling by intercomparing models with different levels of complexity, in terms of the soil-water-vegetation-atmosphere transfer processes involved. It sets to explore the advantages and disadvantages of these different models, while aiming at reducing the reliance on in situ observations for model parameterization and taking full advantage of UAS, airborne, and satellite observations. It focuses on the parameterization of the minimalist soil moisture models, coupled soil and plant models, and crop growth models.
WB5 will then demonstrate the benefits in closing water cycle gaps from the global to local scale, in terms of how to effectively handle spatiotemporal data (from in situ, UAS, and satellites), regarding ecohydrological model calibrations and accuracy evaluations of simulated spatial patterns of ecohydrological variables. Particularly, numerical experiments will be conducted for the calibration of a parsimonious-distributed ecohydrological daily model in ungauged basins using exclusively spatiotemporal information obtained from WB1, WB2, and WB3, to link the scales from plant to plot, sub-catchment, and catchment/basin. WB6 is about disseminating and communicating generated knowledge, data, and tools to water managers, companies, and farmers for actual sustainable water management of their responsible domains. Particularly, to address stakeholders' requirements, iAqueduct will develop an integrative information system (an open source iAqueduct toolbox), which will integrate models, soil parameters, forcing and field-scale observation, and gridded water states and fluxes to support the translation of science knowledge into water productivity information for the smart management of water resources. The goal is to develop potentially effective approaches connecting science to the society, thus influencing citizens towards desirable behavior in water management.
To address iAqueduct challenges, Table 1 lists the essential ecohydrological variables and parameters to be obtained or measured directly by means of various techniques from in situ, UAS, airborne, to satellite.

Downscaling of Satellite Water Cycle Products (WB1)
This WB will focus on the monitoring and downscaling of soil moisture data based on remotely sensed data. The aim is to enhance the level of accuracy and knowledge about soil moisture, in terms of its spatial distribution at surface layers and its vertical distribution at soil profiles (e.g., moving from the skin surface to root depth). Specifically, we aim at the spatial description of soil moisture and the prediction of soil moisture in the root zone. Soil moisture forms a natural link between precipitation, evapotranspiration, and runoff at different spatiotemporal scales. Its spatial and temporal patterns are influenced by several physical features that influence the structure of this pattern (Figure 2a). Moreover, soil moisture is measured through different systems and methodologies, but each of them provides information at specific temporal and spatial scales. In this context, the use of UAS may help to fill the existing gap between field observations and satellite data.
Water 2020, 12, x FOR PEER REVIEW 9 of 37 (3) Generation of high-resolution water cycle products of soil moisture, vegetation patterns, and vegetation stress (sub-meter spatial scale and daily interval). High-resolution maps will be provided with UAS equipped with thermal cameras, multispectral, and hyperspectral cameras. Such data will support the development of downscaling procedures, linking satellite to point measurements for calibration and validation at the selected field sites; (4) Characterization of the spatiotemporal distribution of soil moisture and evapotranspiration processes will be conducted after validation of the high-resolution imagery from UAS with outcomes of field measurements and outputs from ecohydrological models. The proper description of the controlling factors for the spatial variability of soil moisture is crucial to further advance the potential of downscaling methodologies; (5) Downscaling of the remote sensing data to the field scale (from the hectometer to plot scale) can be achieved by using a Bayesian approach exploiting the predicted variance and spatial correlation of the soil moisture process along with the ancillary data derived from UAS and WB2 activities on the physical characteristics of soil and vegetation. In particular, WB2 will support the development of new strategies aimed at the mapping of soil hydraulic and physical characteristics that will enhance the capabilities of soil moisture downscaling procedures (see, e.g., [11,42,43]).

Preliminary Results of Downscaling Surface Soil Moisture
This subsection will focus on surface soil moisture and present the preliminary results of downscaling satellite data with UAS measurements. Because soil moisture is determined by several physical features characterized by strong spatial gradients (e.g., terrain morphology and soil texture) and also temporal variability (e.g., vegetation patterns), these dynamics and features have to be taken into consideration in order to reach a reliable estimate of soil moisture using EO ( Figure 2a). As such, the use of combined technologies may help in describing the spatial patterns of land surface features closely related to soil moisture (or directly soil moisture itself), providing measurements over a range of scales moving from centimeters up to several meters, and thus enabling links to EO data from tens of meters to kilometers (Figures 2b and 3). It is to note from Figure 2b that it is not only about downscaling to the scale of interest but also upscaling. Thus, an evaluation of the physical consistency between different scales and corresponding used downscaling approaches and strategies is always needed [7,70]. Furthermore, ecohydrological model-based simulations at hyper-resolutions can reproduce the scale invariance property of soil moisture, which can be used to link scales from hundreds of meters in the field to tens of kilometers of satellite observations [59]. Big data analytics with machine learning (e.g., random forest, RF) can also effectively downscale satellite observations to UAS/in situ scales [58]. In the following, the RF-based soil moisture downscaling was demonstrated with

Spatial Downscaling Procedures
Before describing a number of procedures to downscale remotely sensed water cycle products to scales suitable for water management purposes, we first present a general framework for water budget closure for a basin by means of the mass conservation equation in the form of [6]: where S is the amount of water stored at the surface and subsurface, P GPCP is the GPCP (Global Precipitation Climatology Project) precipitation data; E SEBS is the SEBS-derived land evapotranspiration is a scaling factor to distribute the observed discharge to each pixel; P i,j , E i,j are GPCP precipitation and SEBS ET for pixel (i, j); and P, E are the mean GPCP precipitation and SEBS ET for the catchment area of interest, all expressed in water depth. It is understood that P GPCP , E SEBS can be replaced by any other remotely sensed similar data. In this mathematical form, Equation (1) can be applied for any catchment after integration. From a remote sensing point of view, ∂S/∂t can be obtained from the GRACE satellite observation of the change in terrestrial water storage (TWS) [62,63]. However, due to the rather coarse spatial resolution (e.g., 300 km), the GRACE TWS is of main utility for large river basins [64]. It is to note that, with earlier versions of GRACE data, the comparison of GRACE TWS with satellite/surface observation-based TWS shows the underestimation of seasonal cycles of TWS [62]. On the other hand, with the new development of local mascon solutions, GRACE data are no longer limited to large-basin hydrology, and are useful for groundwater monitoring [65][66][67]. For a much smaller catchment, we propose to utilize the generalized TOPMODEL concept [68,69] as follows to derive the total drainage: where Q s is the drainage at saturation, A is the area of the catchment, and γ is the spatial average of the soil topographic index; δ is the average of δ with D as the local soil moisture storage deficit (i.e., the difference between the maximum and the actual soil moisture storage), and m a scaling parameter describing the decrease of the subsurface transmissivity T with depth; n is a nondimensional scale parameter of the catchment; a, T 0 , β are the drainage area per unit contour for a specific location within the catchment, the subsurface transmissivity at saturation, and the local slope of the terrain, respectively. Similarly, the local subsurface transmissivity can be written as: and be related to the effective local recharge P i,j − E i,j in the form of: or by inserting Equation (3): Equations (8) and (9) are derived by assuming that the local water table is parallel to the local topography and that the steady state assumption for downslope discharge can be assumed as a power function. Equations (1)-(9) can then be used to link the forcing P i,j − E i,j and drainage Q o to the storage change ∂S ∂t . The obvious challenge in applying such a framework is to quantify the scaling parameters f (−), m, n from the observation scale (e.g., pixel scale) to the scale of management interest (i.e., field or a basin scale). Each of the scaling parameters mainly represents the land-atmospheric processes (e.g., f ), the vertical soil properties (e.g., m), and the lateral hydrological processes (e.g., n), and can be derived from satellite observations of precipitation, evapotranspiration, and soil moisture. Next, we present a number of procedures for downscaling individual variables: (1) Bayesian statistical bias correction of satellite data based on in situ observation. The calibration and validation of coarse-resolution satellite water cycle products at selected field sites with in situ observation is an integral part of this procedure (at the kilometer scale but corrected for spatio-temporal error, e.g., due to topography, soil texture, and climate, cf. those by [19] for precipitation; [26] for evapotranspiration; and [8] for soil moisture); (2) Development of downscaling methods based on Copernicus Sentinel data (from kilometer to hectometer scale). This procedure concerns evapotranspiration and soil moisture (by assuming the precipitation is homogeneous at the kilometer scale). Downscaling will be achieved by the combined use of optical, thermal, and radar data from Sentinel-1, 2, and 3; (3) Generation of high-resolution water cycle products of soil moisture, vegetation patterns, and vegetation stress (sub-meter spatial scale and daily interval). High-resolution maps will be provided with UAS equipped with thermal cameras, multispectral, and hyperspectral cameras. Such data will support the development of downscaling procedures, linking satellite to point measurements for calibration and validation at the selected field sites; (4) Characterization of the spatiotemporal distribution of soil moisture and evapotranspiration processes will be conducted after validation of the high-resolution imagery from UAS with outcomes of field measurements and outputs from ecohydrological models. The proper description of the controlling factors for the spatial variability of soil moisture is crucial to further advance the potential of downscaling methodologies; (5) Downscaling of the remote sensing data to the field scale (from the hectometer to plot scale) can be achieved by using a Bayesian approach exploiting the predicted variance and spatial correlation of the soil moisture process along with the ancillary data derived from UAS and WB2 activities on the physical characteristics of soil and vegetation. In particular, WB2 will support the development of new strategies aimed at the mapping of soil hydraulic and physical characteristics that will enhance the capabilities of soil moisture downscaling procedures (see, e.g., [11,42,43]).

Preliminary Results of Downscaling Surface Soil Moisture
This subsection will focus on surface soil moisture and present the preliminary results of downscaling satellite data with UAS measurements. Because soil moisture is determined by several physical features characterized by strong spatial gradients (e.g., terrain morphology and soil texture) and also temporal variability (e.g., vegetation patterns), these dynamics and features have to be taken into consideration in order to reach a reliable estimate of soil moisture using EO ( Figure 2a). As such, the use of combined technologies may help in describing the spatial patterns of land surface features closely related to soil moisture (or directly soil moisture itself), providing measurements over a range of scales moving from centimeters up to several meters, and thus enabling links to EO data from tens of meters to kilometers (Figures 2b and 3). It is to note from Figure 2b that it is not only about downscaling to the scale of interest but also upscaling. Thus, an evaluation of the physical consistency between different scales and corresponding used downscaling approaches and strategies is always needed [7,70]. Furthermore, ecohydrological model-based simulations at hyper-resolutions can reproduce the scale invariance property of soil moisture, which can be used to link scales from hundreds of meters in the field to tens of kilometers of satellite observations [59]. Big data analytics with machine learning (e.g., random forest, RF) can also effectively downscale satellite observations to UAS/in situ scales [58]. In the following, the RF-based soil moisture downscaling was demonstrated with preliminary results (Figure 3). The RF-based downscaling workflow is depicted in Figure 3a with four steps: Step 1 is to train and test the RF model with both predictors (i.e., land surface features) and soil moisture datasets (Sentinel-1) at a 1-km resolution. The relative importance of predictors ( Figure 3b) shows that the LST (land surface temperature), NDVI (normalized difference vegetation index), and DEM (digital elevation model) are the top three predictors; Step 2 is to train and test the RF model with only the three top predictors as identified in step 1; Step 3 is to apply the trained RF model (from step 2) with the UAS-derived surface features at 15 cm to predict high-resolution soil moisture at 15 cm; Step 4 is to compare the predicted high-resolution soil moisture with in situ measurements. Figure 3c shows the preliminary result of RF-based downscaling of Sentinel-1 soil moisture products at 1 km to 15 cm, taking land surface features derived from UAS (e.g., LST, NDVI and DEM) as predictors over the MFC2 sub-catchment of Alento catchment (hereafter as MFC2-Alento) (see Appendix A.3). Figure 3d shows the comparison between the downscaled soil moisture and in situ measurements. The downscaled soil moisture had 0.07 cm 3 /cm 3 unbiased root mean square error, and its Pearson correlation was 0.42 with the in situ measurements.
shows that the LST (land surface temperature), NDVI (normalized difference vegetation index), and DEM (digital elevation model) are the top three predictors; Step 2 is to train and test the RF model with only the three top predictors as identified in step 1; Step 3 is to apply the trained RF model (from step 2) with the UAS-derived surface features at 15cm to predict high-resolution soil moisture at 15cm; Step 4 is to compare the predicted high-resolution soil moisture with in situ measurements.   Figure 3d shows the comparison between the downscaled soil moisture and in situ measurements. The downscaled soil moisture had 0.07 cm 3 /cm 3 unbiased root mean square error, and its Pearson correlation was 0.42 with the in situ measurements.

From Surface Moisture Information to Profile Soil Moisture
While downscaling coarse-resolution remotely sensed water cycle products to a fine spatial resolution is achievable as described in the previous section, the remote sensing products typically refer to surface information that needs to be transferred to the depth, at least to the root zone, and be linked up with a physically consistent manner. We next describe how to derive profile soil moisture from surface soil moisture information: (1) Prediction of root-zone soil moisture (RZSM) with the SMAR-EnKF (soil moisture analytical relationship-ensemble Kalman filter) [16,54]. Such an approach derives RZSM based on the relative fluctuations of surface soil moisture (SSM) retrieved from the satellite or UAS. Figure 4 shows the workflow for this procedure with satellite data [71], while the same can be applied to UAS and/or downscaled data. Furthermore, given the physically based nature of the model, it will benefit from the information collected on the hydraulic characteristics of the soil (see WB2). Furthermore, the CDF (cumulative distribution function) depth scaling can be also used to derive RZSM from SSM [71]. The prediction of RZSM will help to derive useful information on dynamics of vegetation (e.g., via evapotranspiration). It is to note that both SSM and RZSM can be applied to determine the local soil moisture storage deficit, which is needed for estimating the total discharge (see Equation (6) in Section 3.1.1).
fluctuations of surface soil moisture (SSM) retrieved from the satellite or UAS. Figure 4 shows the workflow for this procedure with satellite data [71], while the same can be applied to UAS and/or downscaled data. Furthermore, given the physically based nature of the model, it will benefit from the information collected on the hydraulic characteristics of the soil (see WB2). Furthermore, the CDF (cumulative distribution function) depth scaling can be also used to derive RZSM from SSM [71].The prediction of RZSM will help to derive useful information on dynamics of vegetation (e.g., via evapotranspiration). It is to note that both SSM and RZSM can be applied to determine the local soil moisture storage deficit, which is needed for estimating the total discharge (see Equation (6) in Section 3.1.1). 2) Other than the above approach, the physical process-based model can be used to simulate SSM and RZSM, and to understand the mechanism behind the relationship of soil property, soil moisture, and evapotranspiration. The STEMMUS (Simultaneous Transfer of Energy, Momentum and Mass in Unsaturated Soil) -SCOPE (Soil Canopy Observation, Photochemistry, and Energy fluxes) numerical soil-water-atmosphere continuum model [72][73][74][75] can be applied to analyze the sensitivities of the predicted RZSM at sites with detailed observation of the soil hydro-thermal properties (soil hydraulic and thermal parameters) and states (profiles of soil moisture and soil temperature and surface radiation, sensible, and latent heat flux, precipitation and other meteorological forcing) [72,[76][77][78][79]. Figure 5 shows the physical processes considered in STEMMUS-SCOPE [80]. This modelling approach will provide high-resolution spatiotemporal patterns of RZSM that can be linked to the spatial distribution patterns of soil properties (see WB2) and evapotranspiration (see WB3) in different catchments. (2) Other than the above approach, the physical process-based model can be used to simulate SSM and RZSM, and to understand the mechanism behind the relationship of soil property, soil moisture, and evapotranspiration. The STEMMUS (Simultaneous Transfer of Energy, Momentum and Mass in Unsaturated Soil) -SCOPE (Soil Canopy Observation, Photochemistry, and Energy fluxes) numerical soil-water-atmosphere continuum model [72][73][74][75] can be applied to analyze the sensitivities of the predicted RZSM at sites with detailed observation of the soil hydro-thermal properties (soil hydraulic and thermal parameters) and states (profiles of soil moisture and soil temperature and surface radiation, sensible, and latent heat flux, precipitation and other meteorological forcing) [72,[76][77][78][79]. Figure 5 shows the physical processes considered in STEMMUS-SCOPE [80]. This modelling approach will provide high-resolution spatiotemporal patterns of RZSM that can be linked to the spatial distribution patterns of soil properties (see WB2) and evapotranspiration (see WB3) in different catchments.

Towards a Protocol for Field-Scale Data Collection
To establish the prediction model (or spectral transfer function) of soil properties, it needs to collect information on soil physical and hydraulic properties, terrain and environmental attributes (topographical, geological, pedological, and land-use/land-cover information together with hydro-meteorological datasets and soil physical and hydraulic properties), and topsoil spectral data [81].
In order to collect accurate and reliable spectral information from the field with physical meanings, the development of a protocol for field spectral measurements under a non-destructive scope is needed. Such a protocol will minimize soil disturbance, which will enhance the characterization of the soil surface hydraulic properties based on spectral data, because sampling soil to the laboratory may disturb the soil surface and hence its hydraulic properties may be compromised. In addition, the protocol to calibrate the field and airborne data and to measure soil temperature and emissivity should be developed as well [81]. Figure 6 shows the appearance of the disturbed and undisturbed soil surfaces at Afeka site, Israel.

Towards a Protocol for Field-Scale Data Collection
To establish the prediction model (or spectral transfer function) of soil properties, it needs to collect information on soil physical and hydraulic properties, terrain and environmental attributes (topographical, geological, pedological, and land-use/land-cover information together with hydrometeorological datasets and soil physical and hydraulic properties), and topsoil spectral data [81].
In order to collect accurate and reliable spectral information from the field with physical meanings, the development of a protocol for field spectral measurements under a non-destructive scope is needed. Such a protocol will minimize soil disturbance, which will enhance the characterization of the soil surface hydraulic properties based on spectral data, because sampling soil to the laboratory may disturb the soil surface and hence its hydraulic properties may be compromised. In addition, the protocol to calibrate the field and airborne data and to measure soil temperature and emissivity should be developed as well [81]. Figure 6 shows the appearance of the disturbed and undisturbed soil surfaces at Afeka site, Israel.   The main physical processes in the STEMMUS-SCOPE continuum model, integrating radiative transfer, vegetation photosynthesis, energy balance, root system dynamic, and soil moisture and soil temperature dynamic. The coupled model integrates vegetation photosynthesis and transfer of energy, mass, and momentum in the soil-vegetation system, via a simplified 1-D root growth model and a resistance scheme (from soil, through root zones and plants, to atmosphere) [80].

Towards a Protocol for Field-Scale Data Collection
To establish the prediction model (or spectral transfer function) of soil properties, it needs to collect information on soil physical and hydraulic properties, terrain and environmental attributes (topographical, geological, pedological, and land-use/land-cover information together with hydrometeorological datasets and soil physical and hydraulic properties), and topsoil spectral data [81].
In order to collect accurate and reliable spectral information from the field with physical meanings, the development of a protocol for field spectral measurements under a non-destructive scope is needed. Such a protocol will minimize soil disturbance, which will enhance the characterization of the soil surface hydraulic properties based on spectral data, because sampling soil to the laboratory may disturb the soil surface and hence its hydraulic properties may be compromised. In addition, the protocol to calibrate the field and airborne data and to measure soil temperature and emissivity should be developed as well [81]. Figure 6 shows the appearance of the disturbed and undisturbed soil surfaces at Afeka site, Israel.

Preliminary Results of Soil Spectroscopy and Hyperspectral Remote Sensing
The prediction model established from the field soil spectral library (SSL) can be upscaled to hyperspectral remote sensing. To achieve this, the soil spectral measurement performed in the field will be compared with the one acquired in the laboratory and from the remote sensing sensors (i.e., both airborne and spaceborne). We will focus on the spectral signature related to (a) soil properties, which are routinely measured and influence the soil hydrological processes, such as soil texture and organic matter content, and also (b) directly the soil hydraulic properties, such as soil water retention, hydraulic conductivity, and water infiltration in soil. In addition, the spectral signature in hyperspectral remote sensing (e.g., visible to near infrared) is expected to be extended to the thermal region.
As the first step, this WB has analyzed the relationship between the spectral information of the undisturbed soil surface and water infiltration into the soil. This soil hydraulic property is highly correlated with water runoff and soil erosion and therefore is important for the description of soil hydrological processes. This is done with a dataset containing 69 soil samples taken from the study areas of MFC2-Alento, Sde Yoav, and from an urban area with exposed soil in the neighborhood of Afeka in the city of Tel Aviv, Israel (see Appendix A.6).
This dataset contains infiltration rate measurements, laboratory spectral measurements, and field spectral measurements of undisturbed soil surfaces. For the measurements of the infiltration rate, we used the MiniDisk Infiltrometer (METER Group Inc., Pullman, WA, USA [83]), and for the reflectance measurements, we used an ASD Spectrometer.
For the field reflectance measurements, we connected the ASD spectrometer to SoilPRO ® [82], in order to obtain optimal spectral measurements in the field, and to neutralize atmospheric attenuation. SoilPRO ® (US patent number 10,473,580 B2) is an apparatus that can be connected to any portable spectrometer in order to extract the un-disturbed reflectance properties of soil field condition with a near laboratory quality. It consists of a large and lightweight closed chamber covering a wide surface area with a controllable illumination and a constant geometry [82]. The partial least squares regression (PLSR) models using the Scikitlearn package in Python [84] was then used to estimate the infiltration rate from the spectral measurements.
For every model, we adopted 5 components for estimation of the soil infiltration rate using soil reflectance in the 450-2400 nm spectral range. This analysis was applied to laboratory spectral measurements as well as to field spectral measurements to explore field sampling issues. Before the application of PLSR, the spectral data was pre-processed using the Savitsky-Golay derivative [85]. The Savitsky-Golay first derivative is a pre-processing method to calculate the variation of the measured reflectance in a given wavelength in relation to its neighbor bands. This pre-processing is a good alternative to enhance spectral properties/signals and reduces physical effects [85][86][87].
In Figure 7a, the result using the field-based model is presented, Figure 7b presents the result using the lab-based model, and Figure 7c presents a histogram of frequencies with the measured infiltration rate (cm/sec) in the different study areas. The PLSR model that was generated using the non-disturbed samples (at field) demonstrated that it is possible to use different soil types and still develop excellent models. The field-based model predicted the infiltration rate much better than the lab-based model. This is because with the application of laboratory protocols, we lost important information of the soil crust for the estimation of soil infiltration rate (see Figure 6).

Preliminary Results of Soil Spectroscopy and Hyperspectral Remote Sensing
The prediction model established from the field soil spectral library (SSL) can be upscaled to hyperspectral remote sensing. To achieve this, the soil spectral measurement performed in the field will be compared with the one acquired in the laboratory and from the remote sensing sensors (i.e., both airborne and spaceborne). We will focus on the spectral signature related to a) soil properties, which are routinely measured and influence the soil hydrological processes, such as soil texture and organic matter content, and also b) directly the soil hydraulic properties, such as soil water retention, hydraulic conductivity, and water infiltration in soil. In addition, the spectral signature in hyperspectral remote sensing (e.g., visible to near infrared) is expected to be extended to the thermal region.
As the first step, this WB has analyzed the relationship between the spectral information of the undisturbed soil surface and water infiltration into the soil. This soil hydraulic property is highly correlated with water runoff and soil erosion and therefore is important for the description of soil hydrological processes. This is done with a dataset containing 69 soil samples taken from the study areas of MFC2-Alento, Sde Yoav, and from an urban area with exposed soil in the neighborhood of Afeka in the city of Tel Aviv, Israel (see Appendix A.6).
This dataset contains infiltration rate measurements, laboratory spectral measurements, and field spectral measurements of undisturbed soil surfaces. For the measurements of the infiltration rate, we used the MiniDisk Infiltrometer (METER Group Inc., Pullman, WA, USA [83]), and for the reflectance measurements, we used an ASD Spectrometer.
For the field reflectance measurements, we connected the ASD spectrometer to SoilPRO ® [82], in order to obtain optimal spectral measurements in the field, and to neutralize atmospheric attenuation. SoilPRO ® (US patent number 10,473,580 B2) is an apparatus that can be connected to any portable spectrometer in order to extract the un-disturbed reflectance properties of soil field condition with a near laboratory quality. It consists of a large and lightweight closed chamber covering a wide surface area with a controllable illumination and a constant geometry [82]. The partial least squares regression (PLSR) models using the Scikitlearn package in Python [84] was then used to estimate the infiltration rate from the spectral measurements.
For every model, we adopted 5 components for estimation of the soil infiltration rate using soil reflectance in the 450-2400 nm spectral range. This analysis was applied to laboratory spectral measurements as well as to field spectral measurements to explore field sampling issues. Before the application of PLSR, the spectral data was pre-processed using the Savitsky-Golay derivative [85]. The Savitsky-Golay first derivative is a pre-processing method to calculate the variation of the measured reflectance in a given wavelength in relation to its neighbor bands. This pre-processing is a good alternative to enhance spectral properties/signals and reduces physical effects [85][86][87]. In Figure 7a, the result using the field-based model is presented, Figure 7b presents the result using the lab-based model, and Figure 7c presents a histogram of frequencies with the measured infiltration rate (cm/sec) in the different study areas. The PLSR model that was generated using the non-disturbed samples (at field) demonstrated that it is possible to use different soil types and still Given that the soil types in the study areas in question are very diverse, the nature of the samples varies. This implicates that field-spectral data has the potential to predict the infiltration rate using a generic approach for different soils using field-based spectral models (e.g., with SoilPRO ® ). Based on the current results, the next step will be to examine the adaptation of field-based spectral models to UAS and satellite platforms specifically in the study area of MFC2-Alento (see Appendix A.3).

Basic and Advanced Pedotransfer Functions
With the collected data in MFC2-Alento, the application and evaluation of already established pedotransfer functions (PTFs) [88][89][90] will be carried out to calculate soil hydraulic parameters. The 3-D Soil Hydraulic Database of Europe at a 250-m resolution [91] will be used as a baseline dataset. As such, this task will explore if and to what extent the predictive capability of these basic PTFs can be suitably improved, through in situ and remote measurements of spatial patterns of land cover, for the mapping of soil hydraulic and thermal properties [60].
Basic and advanced PTFs will enable soil hydraulic and thermal parameters (SHP/STP) to be estimated from spectral signatures and the knowledge of near-surface soil moisture dynamics. The world Soil Spectral Library [92], European Spectral Soil Library (LUCAS) [93], and some local SSL (e.g., the GEO-CRADLE Mediterranean Balkan SSL) will be used to generate global to local spectral-based models to assess soil properties. Spectral transfer functions (STFs) [94] will be derived to predict soil properties (SHP/STP) from high spatial-resolution EO data.
As the first step, we will validate STFs on the test set and samples of the MFC2-Alento catchment. We collected disturbed soil sample and undisturbed soil cores in the MFC2-Alento catchment at 20 locations, corresponding to the positions of the wireless sensor network end-devices (SoilNET) [95] ( Figure 8). We measured in the laboratory the soil particle-size distribution, oven-dry bulk density, soil organic carbon content, and the hydraulic properties of soil-water retention and hydraulic conductivity at the full suction range. Furthermore, we also acquired visible (VIS), hyperspectral, and thermal images with the UAS platform and conducted spectral analysis in the laboratory and in the field ( Table 2, Figures 8 and 9). These data will be used to relate soil spectral information with soil basic and hydrothermal properties.     Figure 8 shows the hydrographic basin of MFC2-Alento (detailed geographical location is referred to Appendix A.3), the wireless sensor network, the spatial coverage of RGB (red, green, blue) VIS, thermal, multispectral, and hyperspectral UAS imageries. Figure 9a shows the RGB image of MFC2-Alento, taken at noon time, Figure 9b the multispectral image taken at afternoon, and Figure 9c the hyperspectral image taken in the morning and the corresponding hyperspectral data cube and mean spectral curves of forest, grass, and soil features. The timing of UAS flights follow the best practice according to local conditions. Figure 9d shows the workflow on the combined use of different sources of data to produce soil texture information and corresponding soil hydro-thermal properties for the top 5 cm and root zone layers.   Figure 8 shows the hydrographic basin of MFC2-Alento (detailed geographical location is referred to Appendix A.3), the wireless sensor network, the spatial coverage of RGB (red, green, blue) VIS, thermal, multispectral, and hyperspectral UAS imageries. Figure 9a shows the RGB image of MFC2-Alento, taken at noon time, Figure 9b the multispectral image taken at afternoon, and Figure  9c the hyperspectral image taken in the morning and the corresponding hyperspectral data cube and mean spectral curves of forest, grass, and soil features. The timing of UAS flights follow the best practice according to local conditions. Figure 9d shows the workflow on the combined use of different sources of data to produce soil texture information and corresponding soil hydro-thermal properties for the top 5cm and root zone layers.

Linking Soil Properties, Soil Moisture, and Evapotranspiration (WB3)
With the downscaled water cycle products, as well as soil hydrothermal properties, this WB will retrieve field-and grid-specific relationship functions between soil properties, soil moisture, and evapotranspiration and their generalization. With the pilot study at MFC2-Alento, the field and grid-

Linking Soil Properties, Soil Moisture, and Evapotranspiration (WB3)
With the downscaled water cycle products, as well as soil hydrothermal properties, this WB will retrieve field-and grid-specific relationship functions between soil properties, soil moisture, and evapotranspiration and their generalization. With the pilot study at MFC2-Alento, the field and grid-specific relationship function will be firstly derived with hierarchical multi-scale EO data, i.e., from the local scale (by the wireless sensor network), to the field scale (by the cosmic-ray neutron probe), to the catchment scale (e.g., by UAS, Sentinel 1-2-3, etc.), and to the regional scale (e.g., by SMOS, SMAP, MODIS, etc.). Both upscaling and downscaling relationships will be developed to link soil properties, soil moisture, and evapotranspiration (see Figure 2). Tests will be done at other sites to validate the proposed procedures.
It is expected that different selected sites (i.e., with different climate zones) will have different specific relationship functions between soil properties, soil moisture, and evapotranspiration. However, as iAqueduct work is still in progress, we are not able to show the sensitivity of such relationship functions under different climate zones (and with the wide-enough variation range of soil properties) from the selected iAqueduct field sites. Nevertheless, in order to demonstrate this concept, we will use existing data from WaPOR (Water Productivity Open-access portal), FAO's (Food and Agriculture Organization) portal to monitor water productivity through open-access of remotely sensed derived data [97][98][99]. It is to note that the WaPOR data coverage is over Africa, which is different from the iAqueduct study sites (see Figure 1 and Appendices A.1-A.6).

Example: WaPOR Database
WaPOR includes the products of soil moisture and evapotranspiration over Africa and the Middle East between 2009 and 2018 at the spatial resolution of 250 m [98,99]. Figure 10 shows that the actual evapotranspiration (AET) increases generally with the soil moisture across all climate zones over Africa. Nevertheless, there are different characteristics within different climate zones.

Approach of iAqueduct
It is noted that the WaPOR relative soil moisture is derived using the trapezoid approach, which needs calibrations for different soils and vegetation covers, due to its assumption of constant meteorology and constant resistance (to heat transfer) [98,99]. The iAqueduct will use a more consistent way with the surface energy balance system (SEBS) equations [25], by which the relative evaporation (RE) can be calculated as [100]: where H is the sensible heat flux, , are the sensible heat flux at the wet limit (i.e., the surface is completely wet and the AET is not limited by soil moisture) and dry limit (i.e., there is no AET due to absence of soil moisture), respectively. As = ( − )/ ℎ ( − air density, − specific heat capacity of air, − land surface temperature, − air temperature, ℎ − aerodynamic resistance for heat transfer), if we assumed the climate (i.e., assuming constant) and land cover (i.e., assuming ℎ constant) do not change spatially over the area of interest, we can rewrite RE as below: where is modulated by SM; and _ , _ are land surface temperatures at the wet limit and dry limit. Equation (11) then describes the theoretical boundaries of the trapezoid approach (i.e., in the space of versus the albedo or vegetation index), when considering the uncertainties raised by assuming and ℎ constants can be addressed by calibrating the system under different surface conditions. Because can be easily measured at various spatial resolutions by remote sensors, For the hot (BWh) and cold (BWk) desert climates, one can see a slight positive relationship between soil moisture (SM) and actual evapotranspiration (AET), indicating the AET is limited by soil moisture for these arid zones, with a maximum rate of 1.5 mm d −1 . For the hot (BSh) and cold (BSk) semi-arid climates, the AET-SM relationship is generally positive with flat slopes alike for desert climates. The BSk is close to desert climates, as it is simultaneously limited by available energy (i.e., cold semi-arid climates) and partially by SM, while BSh has a higher AET as it is not limited by energy and has higher SMs than arid zones.
For the humid subtropical climates with dry winter (Cwa) or fully humid ones without a dry season (Cfa), it is clear that the AET is not limited by soil moisture and has no or a slight negative relationship between SM and AET, respectively. For the hot-summer Mediterranean climate (Csa) and temperate oceanic climates (Cfb), AET-SM is clearly negative, indicating the control of available energy on AET. While for the tropical wet savanna (Aw) and tropical dry savanna (As) climates, the tropical monsoon (Am), and tropical rainforest (Af) climates, the AET is generally high, but there is no relationship with SM.

Approach of iAqueduct
It is noted that the WaPOR relative soil moisture is derived using the trapezoid approach, which needs calibrations for different soils and vegetation covers, due to its assumption of constant meteorology and constant resistance (to heat transfer) [98,99]. The iAqueduct will use a more consistent way with the surface energy balance system (SEBS) equations [25], by which the relative evaporation (RE) can be calculated as [100]: where H is the sensible heat flux, H wet , H dry are the sensible heat flux at the wet limit (i.e., the surface is completely wet and the AET is not limited by soil moisture) and dry limit (i.e., there is no AET due to absence of soil moisture), respectively. As H = ρc P (T s − T a )/r ah (ρ− air density, c P −specific heat capacity of air, T s −land surface temperature, T a −air temperature, r ah −aerodynamic resistance for heat transfer), if we assumed the climate (i.e., assuming T a constant) and land cover (i.e., assuming r ah constant) do not change spatially over the area of interest, we can rewrite RE as below: where T s is modulated by SM; and T s_dry , T s_wet are land surface temperatures at the wet limit and dry limit. Equation (11) then describes the theoretical boundaries of the trapezoid approach (i.e., in the space of T s versus the albedo or vegetation index), when considering the uncertainties raised by assuming T a and r ah constants can be addressed by calibrating the system under different surface conditions. Because T s can be easily measured at various spatial resolutions by remote sensors, Equation (11) will relate ET and SM accordingly, which can be expressed as a sigmoid function by considering its dependence on soil texture [80,100]. Such a soil texture-dependent sigmoid function can be derived from multiscale data. It is expected with the availability of multiscale hierarchical information of soil moisture and evapotranspiration (e.g., from satellite, UAS, in situ measurements) for the selected study sites of iAqueduct (see Figure 1 and Appendix A), both upscaling and downscaling relationships can be developed to link soil properties, soil moisture, and evapotranspiration. Figure 11 shows an example of global evapotranspiration derived from MODIS satellite data in combination with global meteorological information, using the SEBS model [27]. This daily evapotranspiration dataset is at spatial resolution of 1 km. With Equation (11) and the ancillary data used to produce Figure 11, it is expected that the relationship function, as learned in iAqueduct, can be extended to the pan-European and global scale. Such a soil texture-dependent sigmoid function can be derived from multiscale data. It is expected with the availability of multiscale hierarchical information of soil moisture and evapotranspiration (e.g., from satellite, UAS, in situ measurements) for the selected study sites of iAqueduct (see Figure 1 and Appendix A), both upscaling and downscaling relationships can be developed to link soil properties, soil moisture, and evapotranspiration. Figure 11 shows an example of global evapotranspiration derived from MODIS satellite data in combination with global meteorological information, using the SEBS model [27]. This daily evapotranspiration dataset is at spatial resolution of 1km. With Equation (11) and the ancillary data used to produce Figure 11, it is expected that the relationship function, as learned in iAqueduct, can be extended to the pan-European and global scale.

Developing Plant-and Plot-Level Ecohydrological Models Using Remote Sensing Information (WB4)
Both ecologists and hydrologists have long recognized the important interactions between water availability and ecosystems [101]. Vegetation represents the main return path of terrestrial water to the atmosphere [102] and modifies the surface albedo and roughness [103], thereby the ecosystem species composition and structure exert strong controls on hydrologic fluxes at multiple scales [101]. These plant-water interactions have been captured by plant-and plot-level ecohydrological models, ranging in complexity from minimalist soil water balances to detailed descriptions of the water, Figure 11. An example of global evapotranspiration derived from the MODIS satellite data in combination with global meteorological information using the SEBS model [27].

Developing Plant-and Plot-Level Ecohydrological Models Using Remote Sensing Information (WB4)
Both ecologists and hydrologists have long recognized the important interactions between water availability and ecosystems [101]. Vegetation represents the main return path of terrestrial water to the atmosphere [102] and modifies the surface albedo and roughness [103], thereby the ecosystem species composition and structure exert strong controls on hydrologic fluxes at multiple scales [101]. These plant-water interactions have been captured by plant-and plot-level ecohydrological models, ranging in complexity from minimalist soil water balances to detailed descriptions of the water, energy, carbon, and momentum fluxes in the soil, canopy, and atmosphere (see [104,105] for reviews). WB4 will use information retrieved from WBs 1 through 3 to develop and parameterize ecohydrological models of increasing complexity at the plot and plant level and evaluate to what extent this information enhances descriptions of the soil-plant-atmosphere system to inform resource management.
The soil 'bucket filling' model (e.g., [106,107]) is the simplest description of interactions between soil water and vegetation. This approach mechanistically represents the soil water balance over the rooting zone, generally at a daily time scale, but the effects of vegetation are included via an empirical closure, linking directly soil water availability to losses via evapotranspiration [107][108][109]. 'Bucket filling' models have often been used to investigate the role of rainfall unpredictability, because, thanks to their simplicity, the probability density function of soil moisture [107] and plant water stress [110] can be obtained analytically. In some cases, this minimalist description has been linked to the plant carbon balance, to quantify the assimilation rate of carbon dioxide (CO 2 ) and its probability density function [111][112][113]. The minimalist soil moisture balance has been extended to include irrigation [114] and obtain analytically the probability density function of crop yield and irrigation requirements, thus directly linking to the management of water resources [115]. Once parameterized, these models provide a simple and practical output for WB6.
The 'big-leaf' approximation is a more realistic description of vegetation effects on the soil water balance. This approach couples lumped soil moisture dynamics to lumped vegetation models at daily or sub-daily time scales (see [116] for an example). 'Big-leaf' models up-scale mechanistic descriptions of leaf-level carbon and water fluxes to the whole canopy, assuming the entire canopy is subject to the same conditions and behaves in the same way [117,118]. 'Big-leaf' models allows also coupling the soil-vegetation water balance with the plant carbon balance, using various levels of complexity (see examples with increasing physiological detail in [105,113,[119][120][121]). The most mechanistic models include a canopy energy balance (e.g., [75,122]), describing explicitly the interactions between leaf activity and leaf temperature; or carbon and nitrogen dynamics, thus accounting for the cycling of critical resources also in the soil (see [123] for some examples). In addition, to enhance model realism while maintaining low parameter and computational requirements, vegetation and soil can be stratified in two or more discrete layers [123]. The more detailed models include energy, carbon, and possibly nitrogen flux information beyond the water flux provided by 'bucket-filling' models. The 'big-leaf' approximation can be a valuable alternative to 'bucket filling' models in WB6 for managers of agricultural systems, when information for their parameterization is available.
Vertically distributed models represent the highest level of detail, where water and energy balances are fully resolved and discretized over thin layers. Vertically distributed canopy models are often based on the same mechanistic principles informing big leaf models, but they explicitly describe water transport within the soil-plant-atmosphere continuum (see [124] for a review) and the interactions between leaves at different positions inside the canopy and the surrounding conditions. Similarly to the lumped models, vertically distributed models in the soil can focus on the soil water dynamics, incorporating a simplified description of the aboveground component [73,125,126], or considering the full above-and belowground details (e.g., [123,127,128]).
Increasing levels of detail can enhance the model performance but are inevitably associated to increasing parameter requirements and hence difficulties in parameterization when specific in situ observations are not available. In addition to hydrometeorological variables and soil moisture, information about vegetation structure, biomass, temperature, and activity can be retrieved from the iAqueduct platform at various temporal and spatial scales (Table 1), to develop simple to complex ecohydrological models. 'Bucket-filling' soil moisture balances require basic information on soil (soil porosity and hydraulic conductivity) and vegetation (active rooting depth, maximum evapotranspiration rate, and response to water stress). These ecohydrological parameters have been inferred successfully from point-, footprint-, and satellite-scale observations [129] and improve the description of soil moisture dynamics compared to reference constants [130]. Coupled soil and plant models, including a two-pool soil moisture model and vegetation dynamics relying on the radiation use efficiency approach, have been parameterized with satellite data [131]. Most recently, remote sensing has also contributed to the development of crop growth models (see [132,133] for a review).
The work within WB4 will follow these lines, exploiting downscaled data from WB1, data collected in WB2, and information derived from them in WB3. Particular attention will be devoted to the appropriate characterization of vegetation activity and its interactions with soil moisture dynamics. The ultimate goal is to explore the advantages and disadvantages of the different approaches to modelling soil-vegetation-atmosphere interactions when aiming at reducing the reliance on in situ observations and at taking full advantage of UAS, airborne, and satellite observations. The most promising models will be included in the iAqueduct toolbox and be available to improve distributed ecohydrological models in WB5.

Improving Distributed Catchment-Scale Ecohydrological Models Using Spatial Information (WB5)
Leaf-to plot-level plant activity affects and is affected by the spatial scale. The heterogeneity of soil moisture (vertical and horizontal) depends on the distribution of plant types [38,134,135]. At the same time, plants need water to survive, and thus, the composition and structure of plant communities are directly influenced by spatiotemporal patterns in water availability [136]. Piedallu et al. [137] demonstrated how an adequate spatial distribution of soil moisture can improve the performance of ecohydrological models by refining the prediction of the spatial distributions of species. It is obvious that the water cycle processes and their evolution in time and space are both a cause and consequence of vegetation [138].
Despite the pivotal role of plants on the water cycle, traditional hydrological models have neglected most plant-water interactions and vegetation dynamics [139,140]. Most traditional hydrological models include vegetation as a static parameter, which represents observed discharge fairly well at the catchment outlet [141]. As a consequence, there is still a lack of knowledge transfer, and the inclusion of vegetation in spatially distributed models remains challenging [35]. For instance, many recent large-scale models still incorporate vegetation as a constant or by using fixed phenology [142,143]. Fortunately, although including vegetation dynamics is not the norm yet, this trend is changing. During the last decade, new improved ecohydrological models have emerged, e.g., RHESSyS [144], SWIM [145], GEOTOP [146], LPJ-GUESS coupled with TOPMODEL [147], TETIS [136], STEMMUS-SCOPE [80], and STEMMUS-TeC [148].
The calibration of a hydrological model has traditionally relied on the temporal variation of the discharge at the catchment outlet. However, the discharge provides only limited insight on the spatial behavior of the catchment [149]. The development of distributed ecohydrological models and the availability of spatiotemporal data appear as a key alternative to overcome those limitations and can facilitate a spatial pattern-oriented model calibration [136]. In recent years, some researchers have used different sources of remotely sensed information from satellites to improve the implementation of ecohydrological models. For example, actual evapotranspiration was used by Immerzeel and Droogers  [159].
In this WB, iAqueduct will advance how to effectively handle the multiscale multi-source water cycle products, regarding model calibration, and how to evaluate the accuracy of the simulated spatial patterns of vegetation and water states, at sub-catchment and catchment scales. Numerical experiments will be conducted for the calibration of a parsimonious distributed ecohydrological daily model in ungauged basins using exclusively spatiotemporal information obtained from WB1, WB2, and WB3. As such, this WB links the scales from the plant to plot, subcatchment, and catchment/basin, respectively, with the representative size of 1-10 m 2 to 50-500 m 2 , 1-10 km 2 , and >100 km 2 , as derived by means of, for example, TDR observations, to cosmic ray/drone observations, drone/satellite, and satellite observations, respectively.

Summary
The main water challenge for many regions, including Europe, is to achieve the sustainable exploitation of natural and agricultural ecosystems, while enabling economic growth, under current and future climates. Addressing such challenge requires, among others, developing tools for sustainable integrative management of water resources, establishing networks and information sharing among existing research facilities/field labs and disciplines, and connecting science to society [3].
With the established framework of methodology and approaches, iAqueduct will enable the understanding of the space-time variability of EO data (e.g., regarding soil physical characteristics, soil moisture, and evapotranspiration fluxes) from the in situ/plot scale, to field and regional, and to global scales. As such, it will enable local-scale soil, vegetation, and water state information to be derived using satellite, UAS, in situ observations, ecohydrological modelling, and big data analytics tools, for sustainable integrative water management under climate change.
Several case studies are proposed at the local level, which cover a variety of climates, and hydrological and soil conditions, ensuring the generality of the methods developed. Meanwhile, iAqueduct establishes networks and information-sharing practices among existing research facilities/field labs, analytical methods, monitoring tools, and programs. Such knowledge, technology, and data hub complements available field sites with data products from satellite observations (e.g., bias corrected with in situ observation), and harmonizes protocols for the retrievals of soil properties and scaling between soil moisture and evapotranspiration for pan-European applications.
Furthermore, an end-to-end system will be developed to translate scientific data and knowledge into tailored water productivity information, establish a science-policy-business-society interface to allow for continuous dialogues and interactions across different scales and levels, influencing stakeholders towards desirable behaviors for sustainable water management. In the following section, we delineate the iAqueduct toolbox as an end-to-end system connecting science to society, and discuss the challenges to address science questions framing the scope of iAqueduct.

iAqueduct Toolbox
To address stakeholders' requirement (Section 2.1, Appendix A), WB6 will develop an integrative information flow to close the gaps between satellite observations of the water cycle and local sustainable management of water resources (iAqueduct). The generic iAqueduct open source toolbox will integrate models, soil parameters, forcing, and field-scale observation, and gridded water states and fluxes to support the translation of science knowledge into water productivity information for the sustainable management of water resources.
The functionality of iAqueduct toolbox centers on: (1) Monitoring soil, vegetation, and water state information of the selected field sites (Table 1); (2) assessing the long-term and extreme event impacts on local agricultural and natural ecosystems, based on different policy options and climate scenarios; and (3) informing stakeholders and the general public about the actual status of the water resources. The toolbox is to be designed to facilitate the consumption of multiscale multi-source EO data, uses of ecohydrological models, and dissemination of data products via a web-based dashboard. The architecture design and main components were prototyped and the technical details can be found in [160]. The toolbox will use inputs/outputs from WB1-WB5 and make emulated visualization for users and stakeholders, and also as a tool to engage them to help iteratively develop the scenario storylines (policy interventions) for improved scenario assessment.
As a start, scenarios can be worked out for each of the selected sites using a technique developed in the EC CORE-CLIMAX project [6], whereby the distribution of forcings can be derived from the 2018 summer European drought period (as an example) and by replacing the distribution with that of another site, mimicking potential future climate changes and impact to water resources. For example, the observed climate in the Twente region (Appendix A.1) during the drought of 2018 summer will be replaced by that of the Spanish (Appendix A.5) or Italian sites (Appendix A.3) and the spatiotemporal water situation in the Twente region be simulated. In collaboration with the water authority, potential management scenarios will be developed, and citizens will be invited to propose additional measures (e.g., water-saving measures) as a preparation for such a scenario, thus connecting science to the society more effectively and influencing citizens towards desirable behavior.

Challenges
There are challenges to answer the questions framing the scope of iAqueduct (Section 2.1), in order to address stakeholders' requirements: -How to downscale the global water cycle products to the local scale using multiple sources and scales of EO data?
WB1 is specifically designed to address this question, and is going to use multiple sources and scales of EO data, including in situ, UAS, medium-resolution satellites (e.g., Sentinel 1-2-3), and coarse-resolution satellites. Particularly, a number of procedures will be implemented to downscale remotely sensed water cycle products to scales suitable for water management purposes. The challenge here is to understand the physical consistency of the original satellite data and the downscaled one, when compared with in-situ measurements.
Taking soil moisture data as example, SMOS and SMAP are dedicated satellite missions for monitoring soil moisture globally at a 5-cm depth. However, Lv et al. [161,162] found that the soil moisture sensing depth of the L-Band sensor actually varies with the soil moisture itself. This means the direct comparison of in situ measured soil moisture (at 5 cm) to SMOS or SMAP observations may not be the optimal way to do the validation, which brings difficulties to the evaluation of the accuracy of downscaled soil moisture data. If we used the X-band sensor, the soil moisture sensing depth will be only limited to the top 1-2 cm [58], which renders the difficulties to placing soil moisture sensors in situ for a direct comparison with X-band data. These potential mismatches (between what we measure in situ and what is remotely sensed) highlight the need to understand physical processes behind different soil moisture products, which will guide the associated downscaling strategies and methodologies. Similar issues exist for other water cycle products [6,7,70].
-How to explore and apply the downscaled information at the management level for a better understanding of water-energy-soil-vegetation processes?
The EO data alone is not sufficient to understand water-energy-soil-vegetation processes, which requires synergetic use of soil-plant-atmosphere continuum (SPAC) models. There is a wide diversity of SPAC models, ranging from minimalist models to detailed descriptions of the water, energy, carbon, and momentum fluxes in the soil-plant-atmosphere system (see Section 3.4). The challenge is to identify the advantages and disadvantages of these models, while aiming at reducing the reliance on in situ observations for model parameterization and at taking full advantage of UAS, airborne, and satellite observations. Another dimension of this challenge is to provide model inputs at various spatial scales, from plant to plot, sub-catchment, and catchment/basin. As such, it is possible to evaluate the optimal model complexity (e.g., either minimalist or complex model) for different scales. The optical and thermal sensors (e.g., space-based or UAS-based), together with soil spectroscopy and associated pedotransfer functions, can provide the model-required soil hydro-thermal properties at a high resolution (of the order of centimeters). However, there are well-known constraints for sensors in VIS-NIR and TIR regions of the electromagnetic spectrum, the applicability of which depends on vegetation cover. For the purpose of deriving soil texture information and soil hydro-thermal properties, these sensors can only be applied over the non-vegetated soil surface. How to address this constraint is vital to answering the second iAqueduct question, which may be examined based on a global-scale approach to evaluate clean soil areas [163]. - How can such fine-scale information be used to improve the management of soil and water resources?
The aim of iAqueduct is to disseminate and communicate the generated knowledge and tools to water managers, companies, and farmers for actual sustainable water management. In order to be effective, stakeholders need to be engaged actively for the effective transfer of science knowledge and data into marketable tools, services, and/or products for addressing the actual needs for real life water management. This is actually demanding the establishment of an end-to-end system (iAqueduct toolbox), to enable a science-policy-business-society interface for continuous dialogues and interactions with stakeholders. The iAqueduct toolbox will build upon the existing open-source software system, MajiSys water information system [160], and will then be used for robust application (including machine learning algorithms) to the selected sites and also for use by stakeholders. The challenge here is how to facilitate such a dialogue to enable policy learning.
To sum up, iAqueduct will integrate various components from the global water cycle observation to local soil, vegetation, and water states in an open-source water information system, and test and demonstrate their utility at a set of carefully selected research sites for sustainable management of water resources. Furthermore, as iAqueduct links scales from plant to plot, field, catchment, and region, it is expected that the developed approach can be further upscaled to the pan-European scale.
weirs, and pumping stations. To meet this challenge, local and regional monitoring of the actual state of the water system and the anticipation of near future situations are needed. Agricultural water management requires, on the other hand, operational management of soil and water for adequate agricultural productions in the growth seasons at a field level. Other requirements are related to water quality management for nature conservation, including water treatments. A shortage of precipitation (i.e., precipitation minus potential evaporation) is used as a measure for water excess or shortage for the abovementioned tasks.
A regional soil moisture monitoring network has been installed since 2009 in the Twente region of The Netherlands [164], consisting of 20 stations continuously measuring soil moisture and soil temperature over an area of approximately 50 km × 40 km. The main objectives of Twente monitoring network are: (i) To investigate the sensitivity of active and passive microwave data to surface parameters, such as soil moisture, soil temperature, and vegetation cover; (ii) to run, calibrate, and validate new soil moisture retrieval algorithms; and (iii) to study new approaches to upscale soil moisture information from the point to large scale.
The precipitation data are available at the Royal Netherlands Meteorological Institute (KNMI). There are 15 KNMI stations measuring the precipitation in the area of the network since 1950 and the daily accumulated precipitation is available at http://www.knmi.nl/klimatologie/. According to the Koeppen Classification System, the climate in The Netherlands is a warm temperate humid climate (Cfb Climate). Other meteorological data are provided by KNMI at the Twente station located near Enschede. The precipitation is spread all over the year with an average of approximately 760 mm per year. The monthly average air temperature ranges between 3 • C in January to approximately 17 • C in July.
Detailed Land use map of Twente area are available as the Atlas of Overijssel by the Province of Overijssel (http://gisopenbaar.overijssel.nl/website/bodematlas/bodematlas.html). There are four main soil types in Twente: Sandy soils rich or poor of loam, loamy soils rich or poor of sand, man-made sandy thick earth soils, and peat soils covered by a layer of peat or sand. This information was retrieved from the soil maps (Bodemkaart van Nederland) by Stichting voor Bodemkartering (Wageningen) with a 1:50,000 scale and by Alterra, Wageningen UR (website www.bodemdata.nl). Groundwater monitoring data can be obtained from the Geological Survey of the Netherlands (TNO) (https://www.dinoloket.nl/en). Hydrological and water management information can be obtained at the water authority Vechtstromen (https://www.vechtstromen.nl/). The Twente SMST network serves as an SMOS and SMAP calibration and validation site [33].
Appendix A.2. Zala, Hungary (Cold, Humid Winter, Warm Summer) One of the main threats in the catchment of Zala river is the more frequently occurring extreme weather events at the catchment. Information on the impact of potential future climate changes to water resources and possible management scenarios to adapt to future extreme weather events will be shared with the General Directorate of Water Management, the farmer organizations (e.g., AGRYA, Agrion Top Kft.) and public bodies (e.g., Hungarian Chamber of Agriculture, Zala County, Zala County Office of Agricultural and Rural Development Agency). At the catchment, it is also important to analyse how the transport of fertilizers, pesticides, and herbicides will change due to the extreme weather events. Recently, the amount of nitrate and pesticide in the groundwater at the catchment is close to the threshold value of groundwater pollution, and at a few plots, even exceeds it.
The catchment of river Zala in western Hungary belongs to the watershed of Lake Balaton. The catchment area of the Zala River is 2622 km 2 , it is situated in Zala Hills. Mean discharge of Zala is 5.6 m 3 s −1 . The climate is moderately warm, moderately humid, and the number of sunshine hours per year ranges between 1800 and 2000 h. The mean annual temperature of the region is about 10 • C. The average amount of rainfall is between 600 and 700 mm year −1 .
In total, 37% of the total catchment area is arable land, which is much lower than the national average; 27% is forest, which exceeds the national average; 15% of the land is under grassland management; 5% is horticulture; 3% is pomiculture; 2% is viticulture; and 1% is reed management and fish farming. The "Kis-Balaton" nature conservation area, which is a wetland under protection of the Ramsar Convention habitat, is situated within the watershed of river Zala. The dominant soil types are Luvisols and Cambisols. Gleysols and Histosols occur in poorly drained valley bottoms [165,166].
Long-term data is available on water quantity and the water quality of Zala river, e.g., runoff, nutrients, dissolved oxygen, total dissolved solids, and water temperature. The depth of the groundwater level and daily rainfall are registered at 39 locations at the catchment of Zala river. Soil moisture is monitored with a TDR/MUX/mpts meter (Easy Test [167]) at two soil profiles till a 0.9-m depth at Vése and Keszthely since April 2018. Time series meteorological data for the sites is also available. Map of soil organic carbon content, texture, calcium-carbonate content, and pH [168] (http://dosoremi.hu/table.html) are available at a 100-m resolution at six soil depths up to 2 m. Soil water content at saturation, field capacity, wilting point (https://www.mta-taki.hu/en/kh124765), and saturated hydraulic conductivity have been mapped at a 100-m resolution at three soil depths up to 0.90 m [169]. For the mapping of soil properties, over 150 environmental parameters were collected to describe the topography, climate, parent material, state of vegetation, and land cover of the catchment. This study area partly belongs to the "Cilento and Vallo di Diano" National Park, the largest national park in Italy, and is included as a representative site within the UNESCO-HELP program. The Alento River Catchment is usually split in the Upper Alento, a hilly and mountain marginal area that suffered from severe land abandonment and subsequent land-use changes, and the Lower Alento, characterized by a flourishing economy especially because of tourism along the entire coastline. A system of barrages, the largest being the "Piano della Rocca" earthen dam, were built and is managed by the "Velia" Bureau of Reclamation to increase irrigated agriculture and the quality of livestock methods, hence reducing the gap between the two parts of the catchment.
However, as in most water-stressed zones of the Mediterranean belt, this area is experiencing an excessive demand for water partly because of the competition among different users, which could yield conflicts among them, especially during summer. Decision-makers and stakeholders are now concerned about future benefits and constraints deriving from the changes observed in land uses and climate seasonality, and are therefore interested in addressing the following main issues: (a) Predicting the storage capacity of the artificial reservoirs in view of projected climate and land-use changes so as to meet short-and medium-term water requirements from households, agriculture, tourism, and hydropower generation; (b) promoting the most effective demand-side adaptation options; and (c) identifying optimal land resource management to ensure adequate water availability to all sectors, reduce fire risk during the prolonged dry seasons, and, at the same time, alleviate natural hazards, such as flooding and soil erosion, during the wet season.
To meet these needs, the Alento River catchment is becoming a science-driven critical zone observatory (CZO), with a major aim of supporting the issues of rural environmental protection and sustainable management of natural resources. The "Alento" CZO not only relies on background geological, pedological, and hydrological studies carried out over the last decades but also benefits from a series of investigations currently underway in the Upper Alento catchment [95]. Since 2016, wireless sensor networks (WSNs) and cosmic-ray neutron probes (CRNPs) monitor soil moisture in two small sub-catchments, named MFC2 and GOR1, having different topographic, pedological, and land-use characteristics as well as slightly different weather conditions.
More hydrologically oriented investigations have been carried out in three sub-catchments (MFC1, MFC2, and GOR1) of UARC, each with drainage areas ranging from 6 to 10 hectares [43]. MFC1 is a small sub-catchment (area of about 5 ha) located near the village of Monteforte Cilento. This experimental area was subject to intense monitoring activities from 2006 to 2011. The MFC2 sub-catchment is situated near MFC1 and both represent typical farmland areas with olive orchards, vineyards, fruit trees and crops. A portion of MFC2 is planted with cherry and walnut trees for wood production only. Another sub-catchment, named GOR1, is instead located close to the rural village of Gorga and reflects a typical forested area with mixed chestnut and oak woods. Apart from their different land-use/land-cover features, the MFC2 and GOR1 test sites should also be viewed as representative of two different hydrogeological settings. On March 2016, the MFC2 and GOR1 sub-catchments were both instrumented with a wireless sensor network (WSN) and a cosmic-ray neutron probe (CRNP).
Each wireless sensor network (SoilNet, Forschungszentrum Jülich, Germany) comprises an array of 20 underground measuring nodes that are distributed in space to account for the local geomorphological and pedological features. The following sensors were embedded at the soil depths of 0.15 and 0.30 m of each node: (i) GS3 sensor (METER, Pullman, WA, USA) to measure simultaneously soil permittivity (converted to soil moisture using empirical calibration equations), soil temperature, and soil electrical conductivity, and (ii) MPS-6 sensor (METER, Pullman, WA, USA) to measure soil matric pressure head. In each of these sub-catchments, one cosmic-ray neutron probe (CRS2000/B, Hydroinnova LLC, Albuquerque, USA) complements the SoilNet sensor network by providing area-averaged soil moisture values over a footprint of approximately 7-14 ha. To our knowledge, these two cosmic-ray neutron probes are the first to be installed and operational since 2016 in a catchment of central/southern Italy. Streamflow gauging stations are operating at the outlet of both the MFC2 and GOR1 sub-catchments. One weather station in close proximity to each sub-catchment acquires rainfall, air temperature, relative humidity, wind speed, and net solar radiation (four-component net radiation sensors, Hukseflux Thermal Sensors, www.hukseflux.com) at hourly time-steps. Wind speed and air temperature are measured at about a 3-m height while solar radiation is measured at a 2-m height. At each weather station, three GS3 sensors and three MPS-6 sensors were installed at depths of 15, 30, and 45 cm. Fiumarella of Corleto belongs to the Basilicata region that is characterized by a significant diversity in terms of climatic conditions. For instance, the mean annual rainfall ranges between 400 and 2000 mm. Such variability reflects the regional hydrological patterns with areas affected by droughts and others that experience several floods and landslides. In this context, the study of river basin hydrology becomes critical from several points of view. The Fiumarella of Corleto is located in the water-rich part of the region that is crucial for the water supply of the region but also for the water supply of the Puglia region, which strongly relies on external resources for their agricultural and economical activities.
The experimental basin "Fiumarella of Corleto", located in Basilicata region (southern Italy), is a tributary of the Sauro river (Agri basin) and has an area of 32.5 km 2 . It is situated in a sub-humid climatic zone with a mean annual rainfall of approximately 720 mm and characterized by hot-humid summers and chilly to mild winters. The interest towards this basin is due to its peculiarities. In fact, the two slopes of the catchment have different land uses: The slope on the left is covered mostly by forests, the slope on the right is covered by agricultural land. In order to characterize with a high level of details the morphology of the two slopes, a DSM of the basin at high-resolution (1 × 1 m) was derived with a LiDAR.
Catchment pedology was investigated through field campaigns and laboratory measurements aimed at identifying the main soil and units of the basin by Santini et al. (1999) [170] and by Romano and Santini, 1997 [171] and Romano and Palladino (2002) [172]. These data were reported in the land cover map elaborated by Carriero et al. (2007) [173] to define the soil hydraulic properties of each unit.
Meteorological variables are monitored on both slopes with the aim to characterize differences between the two. Moreover, an additional rainfall station and streamflow gauge is placed at the basin outlet. Soil moisture is monitored on a transect of about 60 m with a sampling frequency of 1 h. The installed instrumental system consists of a TDR100 system connected to 22 probes located in 11 sampling sites at two different depths of 30 and 60 cm. The datalogger is a CR10X produced by Campbell Scientific that transmits the soil moisture values, elaborated by the TDR100, in real time via the GSM network.
Appendix A.5. Carraixet Creek, Spain (Semiarid, Steppe, Mediterranean) Barranco del Carraixet (or Carraixet Creek) is located in the east coast of Spain, has a catchment area of 314 km 2 , draining directly to the Mediterranean Sea, with a natural park in the upper part of the basin and with anthropogenic pressures in the middle and low basin. The human effect is quite important in this study site: The lowlands are characterized by alternation of the urban and industrial zones and agricultural fields, while the upper part is frequently affected by wildfires and it is a highly frequented leisure zone, subject to multiple pressures (hunters, several outdoor sports, owners, etc.). The climate is semiarid Mediterranean, with a mean annual precipitation of around 400 mm highly variable and potential evapotranspiration of 1100 mm. The hydrology is characterized by low or absent base flow, typical of Mediterranean ephemeral streams. Urban and irrigation water demands are supplied by the aquifer, mainly recharged by the upper catchment. The actual trend of the catchment is towards forest expansion in abandoned lands of the upper part and urbanization in the lower part. The main concern in Carraixet Creek is to improve forest management in order to increase aquifer recharge, increase the forest health, and to better control soil erosion. For this project, we will consider the upper and medium parts of the catchment, with an area of 250 km 2 . Within this area, there is one experimental watershed of 1 km 2 (with 3 meteorological stations, 1 cosmic ray, and 1 flowgauge) and 1 experimental forest plot heavily sensorized. At the catchment scale and operated by the Jucar Basin Water Authority, there is one additional flowgauge station and several raingauges and piezometric observations. Appendix A.6. Kibbutz Sde Yoav and Afeka, Israel (Arid, Dry Hot Summer) Kibbutz Sde Yoav is an agricultural settlement located in south-central Israel, between the cities of Ashkelon, Kyriat Gat, and Kyriat Malakhi. Like Sde Yoav, in Israel, there are several agricultural settlements that were created during the establishment of the state of Israel in order to ensure the food supply. In order to monitor the fields that sustain these settlements, farmers need chemical/physical analyses. However, traditional soil survey methods are expensive, time-consuming, and need high skilled professionals. Moreover, in order to represent these parameters spatially in these large agricultural fields, it is necessary to take several samples for a correct kriging methodology, and the measurements in question varies seasonally. Additionally, farmers cannot see the status of every point of interest in these fields with a common frequency.
Given the lack of rains and the dry climate of the region [174], water is a critical resource that is necessary to manage carefully. Remote sensing is a potential solution for this problem, because it can replace field chemical/physical measurements, and could monitor the infiltration rate in the fields of interest in a spatial scale.
Study sites for the estimation of water infiltration rate using surface spectra include not only Kibbutz Sde Yoav, and Afeka, Israel but also the Alento catchment. The idea is to create a dataset with infiltration rate measurements, with laboratory and field spectral measurements. The study sites from which we collected samples until now are the following: (i) Kibbutz Sde Yoav, Israel (30 Samples): Kibbutz Sde Yoav is an agricultural settlement located in south-central Israel, between the cities of Ashkelon, Kyriat Gat, and Kyriat Malakhi. According to a detailed map of the soils of Israel, the soil type of the study area of Kibbutz Sde Yoav is alluvial, and according to an updated version of the Koeppen climate classification [174], the climate of Sde Yoav is hot-semiarid (Bsh). In this study area, 30 samples were collected. (ii) Afeta, Tel Aviv, Israel (18 Samples): Afeka is a residential neighborhood located in the north of Tel Aviv. The soil type of the study area of Afeka is brown-red sandy soil, and the climate according the classification of Rubel and Kottek, 2010 [174] is hot-summer Mediterranean climate (Csa). From Afeka, we collected 18 samples. In Afeka, we only collected samples for the calibration of the model and to expand our dataset. We did not carry out UAS campaigns there because they are forbidden.
(iii) Alento, Italy (21 Samples): The Alento River Catchment is located in the Campania Region (Salerno Province, Italy). As in Afeka, the climate of Alento according the classification of Rubel and Kottek, 2010 is hot-summer Mediterranean climate (Csa). According to [175] in the book "Soils of Italy", Alento is located in an area characterized by three soil types: Cambisols, Leptosols, and Luvisols. In this book, Costantini and Dazzi remark that this area is characterized hills and mountains on limestones covered by volcanic ashes, including alluvial and coastal plains.