Precise and accurate delineation of flooding areas with synthetic aperture radar (SAR) and multi-spectral (MS) data is challenging because flooded areas are inherently heterogeneous as emergent vegetation (EV) and turbid water (TW) are common. We addressed these challenges by developing and applying a new stepwise sequence of unsupervised and supervised classification methods using both SAR and MS data. The MS and SAR signatures of land and water targets in the study area were evaluated prior to the classification to identify the land and water classes that could be delineated. The delineation based on a simple thresholding method provided a satisfactory estimate of the total flooded area but did not perform well on heterogeneous surface water. To deal with the heterogeneity and fragmentation of water patches, a new unsupervised classification approach based on a combination of thresholding and segmentation (CThS) was developed. Since sandy areas and emergent vegetation could not be classified by the SAR-based unsupervised methods, supervised random forest (RF) classification was applied to a time series of SAR and co-event MS data, both combined and separated. The new stepwise approach was tested for determining the flood extent of two events in Italy. The results showed that all the classification methods applied to MS data outperformed the ones applied to SAR data. Although the supervised RF classification may lead to better accuracies, the CThS (unsupervised) method achieved precision and accuracy comparable to the RF, making it more appropriate for rapid flood mapping due to its ease of implementation.@en

The ambition of H2020 OPERANDUM project is to develop and document Nature Based Solutions (NBS) to mitigate risks associated with hydro-meteorological (HM) hazards. NBS mitigate risks by reducing the vulnerability of a particular system. The aim of this work is to demonstrate the use of multisource remote sensing data in documenting the impact of extreme HM events to advance knowledge on vulnerability and exposure. In particular the focus is to document past impacts due to extreme events selected from a characterization of recent (3 0 years) HM events in 11 Open Air Laboratories (OALs) where co-design, co-development and deployment of NBS are taking place. The impacts were documented by applying a wide spectrum of satellite image data and other, close - range, remote sensing techniques. A better understanding of the consequences due to extreme HM events in a particular area (OALs) is essential to identify elements at risk and expected to provide a reference to evaluate the reduction of vulnerability and mitigation of risks past the completion of NBS.@en

Nature-based solutions (NBS) for hydro-meteorological risks (HMRs) reduction and management are becoming increasingly popular, but challenges such as the lack of well-recognised standard methodologies to evaluate their performance and upscale their implementation remain. We systematically evaluate the current state-of-the art on the models and tools that are utilised for the optimum allocation, design and efficiency evaluation of NBS for five HMRs (flooding, droughts, heatwaves, landslides, and storm surges and coastal erosion). We found that methods to assess the complex issue of NBS efficiency and cost-benefits analysis are still in the development stage and they have only been implemented through the methodologies developed for other purposes such as fluid dynamics models in micro and catchment scale contexts. Of the reviewed numerical models and tools MIKE-SHE, SWMM (for floods), ParFlow-TREES, ACRU, SIMGRO (for droughts), WRF, ENVI-met (for heatwaves), FUNWAVE-TVD, BROOK90 (for landslides), TELEMAC and ADCIRC (for storm surges) are more flexible to evaluate the performance and effectiveness of specific NBS such as wetlands, ponds, trees, parks, grass, green roof/walls, tree roots, vegetations, coral reefs, mangroves, sea grasses, oyster reefs, sea salt marshes, sandy beaches and dunes. We conclude that the models and tools that are capable of assessing the multiple benefits, particularly the performance and cost-effectiveness of NBS for HMR reduction and management are not readily available. Thus, our synthesis of modelling methods can facilitate their selection that can maximise opportunities and refute the current political hesitation of NBS deployment compared with grey solutions for HMR management but also for the provision of a wide range of social and economic co-benefits. However, there is still a need for bespoke modelling tools that can holistically assess the various components of NBS from an HMR reduction and management perspective. Such tools can facilitate impact assessment modelling under different NBS scenarios to build a solid evidence base for upscaling and replicating the implementation of NBS.

@en

A protected natural area in the Emilia Romagna region, Northern Italy is threatened by hydro-meteorological hazards, particularly sea storms. In the last 50 years the northern part of the Bellocchio Park (Sacca Bellocchio II Nature Reserve, Site code EUAPP0072 - Ferrara, Italy) was interested by an intensive urbanization (Lido di Spina) with the realization of infrastructures, e.g. roads and residential settlements. This land use change led to the construction of embankments and to the conversion of wetlands. These modifications, in combination to even more frequent storm surge events increased coastal erosion. In addition, inland flooding caused by storm surges acts with the reduction of the lagoon and the increase of soil salinity. As an example, the last event occurred in December 2020 eroded a large portion of the Bellocchio beach. Co-design, co-development and deployment of NBS solutions to reduce storm surge risk in the Bellocchio Park is one of the objectives of the H2020 project OPEn-air laboRAtories for Nature baseD solUtions to Manage environmental risks (OPERANDUM). BellocchioBellochio park is in fact one of the 10 Open Air Laboratories (OAL) where the evidence of mitigation of hydro-meteorological risk by NBS will be demonstrated by the combination of different models, approaches and data. During the co-design process in the Bellocchio park, potential deployment locations of sand dunes have been identified in collaboration with local authorities devoted to the management of the natural area and to the coast defense (CB and ARSTePC-RER) and an environmental engineering consultant assisting Arpae (IRIS sas). Field visits were devoted to the analysis of the environmental features, strengths and weaknesses of candidate sites. This work aims to explore the usefulness of the combined use of multisource remote sensing and modeling in decision making during the co-design process of a NBS. The impacts of the most intense extreme storm surge events in the last 30 years have been documented by delineating flooded areas along the coast using Synthetic Aperture Radar and Multispectral image data. Coastal erosion has been also described by means of change detection analysis and very high resolution multispectral EO data. This screening has given a picture of areas at the risk, i.e. the area most likely to be affected by storm-surge events. Auxiliary data like Digital Terrain Models has been assimilated in a dedicated model to produce flood maps under different scenarios, i.e. different locations and size of NBS and different intensities of storm surge. The integrated analysis was helpful in defining the priority sites, among the ones defined by the stakeholders and engineers, in term of effectiveness for storm surge risk reduction.@en

To bring to fruition the capability of nature-based solutions (NBS) in mitigating hydro-meteorological risks (HMRs) and facilitate their widespread uptake require a consolidated knowledge-base related to their monitoring methods, efficiency, functioning and the ecosystem services they provide. We attempt to fill this knowledge gap by reviewing and compiling the existing scientific literature on methods, including ground-based measurements (e.g. gauging stations, wireless sensor network) and remote sensing observations (e.g. from topographic LiDAR, multispectral and radar sensors) that have been used and/or can be relevant to monitor the performance of NBS against five HMRs: floods, droughts, heatwaves, landslides, and storm surges and coastal erosion. These can allow the mapping of the risks and impacts of the specific hydro-meteorological events. We found that the selection and application of monitoring methods mostly rely on the particular NBS being monitored, resource availability (e.g. time, budget, space) and type of HMRs. No standalone method currently exists that can allow monitoring the performance of NBS in its broadest view. However, equipments, tools and technologies developed for other purposes, such as for ground-based measurements and atmospheric observations, can be applied to accurately monitor the performance of NBS to mitigate HMRs. We also focused on the capabilities of passive and active remote sensing, pointing out their associated opportunities and difficulties for NBS monitoring application. We conclude that the advancement in airborne and satellite-based remote sensing technology has signified a leap in the systematic monitoring of NBS performance, as well as provided a robust way for the spatial and temporal comparison of NBS intervention versus its absence. This improved performance measurement can support the evaluation of existing uncertainty and scepticism in selecting NBS over the artificially built concrete structures or grey approaches by addressing the questions of performance precariousness. Remote sensing technical developments, however, take time to shift toward a state of operational readiness for monitoring the progress of NBS in place (e.g. green NBS growth rate, their changes and effectiveness through time). More research is required to develop a holistic approach, which could routinely and continually monitor the performance of NBS over a large scale of intervention. This performance evaluation could increase the ecological and socio-economic benefits of NBS, and also create high levels of their acceptance and confidence by overcoming potential scepticism of NBS implementations.

@en

Thermal InfraRed (TIR) image data at high temporal and spatial resolution are required to monitor the rapid development of crops during the growing season, taking into account the fragmentation of most agricultural landscapes. Moreover, integrating high-resolution satellite TIR data to calibrate hydrological models is a powerful information to efficiently monitor crop water use. Conversely, no single sensor meets these combined requirements in the TIR spectral region. Data fusion approaches offer an alternative to exploit observations from multiple sensors, providing image data to meet the combined requirements on spatial and temporal resolution. A novel spatio-temporal data fusion workflow based on a multi-sensor multi-resolution algorithm was developed and applied to generate TIR synthetic image data at high temporal and spatial resolution. The workflow includes two steps: in the first step, synthetic daily radiance images at Top of Atmosphere (TOA) and 30-m spatial resolution (at the ground) are generated using TIR radiometric data at TOA collected by the Moderate Resolution Imaging Spectroradiometer (MODIS) daily 1-km and Landsat 8/TIRS 16-day 30-m. This procedure is applied to two image pairs on different dates. The workflow yields an estimator to generate TIR TOA radiance data on any given date, provided a MODIS radiance image is available. The next step applies constrained unmixing of the 30 m (now considered as low-resolution) TIR images using the information about sub-pixel land-cover obtained from co-registered images at higher spatial resolution in the VNIR (Visible Near InfraRed) spectrum. In our case study, the L8/TIRS synthetic image data were unmixed to the Sentinel 2/MSI with 10 m × 10 m spatial resolution. Two geographically diverse experiments were carried out using the same procedure: one in The Netherlands to evaluate the procedure and other in Puglia (Italy) to generate a time series of the 10-m × 10-m TIR image data product. The validation experiment, where an actual TIRS image was applied as a reference, gave a RMSE value of 35.3 W/(m2 μm sr), which corresponds to a relative value of 8.5% against the TIRS reference values. The results confirm the feasibility of the proposed methodology, which yields a synthetic thermal band to integrate with the multi-spectral data provided by the S2/MSI at 10 m resolution. @en

This study on a large irrigation scheme in Morocco addressed a two-fold question: (a) does the local water management authority adapt water releases to atmospheric water demand ET₀-P? and (b) does crop actual evapotranspiration respond to inter- and intra-annual variability in water releases? We have evaluated the inter-annual variability of ET₀-P during the period 1992-2017 and compared its anomalies (i.e., deviations from average) with anomalies in annual water release. Overall, it appeared that anomalies in water release were consistent with anomalies in ET₀-P. The actual evapotranspiration (ETa) was estimated using a time series of multi-spectral satellite image data by applying the Surface Energy Balance (SEBAL) algorithm in a dry, wet, and reference year. We have determined the quartiles of the distribution of the ET₀-P time series to identify these three years. The dry year was 2015-2016, the wet year was 2014-2015, and the reference (median of ET₀-P) was 2013-2014. Finally, we compared seasonal and annual anomalies in ET₀-P, ETa and release, W_d of irrigation water. In the period 1992-2017, the relative anomalies in ET₀-P and W_d were similar only in 2000-2001 and 2015-2016. The analysis of anomalies in water inflow and stocks confirmed the response in increased W_d following wet years with higher inflow and replenishment of the reservoir. The response of crop water use to the irrigation water supply was evaluated by comparing anomalies in the ratio of actual to maximum ET, i.e., ET_a/ET_c with anomalies in W_d. As regards the Sidi Bennour, Faregh, and Gharbia districts, the relative anomalies in ET_a/ET_c and W_d were consistent, i.e., they had the same sign and comparable magnitude. Overall, the study shows that water delivery W_d responds to inter-annual variability in the pre-irrigation season water inflows into the reservoirs, rather than in ET₀-P. On the other hand, actual crop water use, i.e., ET_a/ET_c, does respond to inter- and intra-annual variability in W_d. This evidence suggests that there is scope for adaptive water management based on a flexible adaptation of water release to inter- and intra-annual variability in water demand.

@en

Small reservoirs play an important role in mining, industries, and agriculture, but storage levels or stage changes are very dynamic. Accurate and up-to-date maps of surface water storage and distribution are invaluable for informing decisions relating to water security, flood monitoring, and water resources management. Satellite remote sensing is an effective way of monitoring the dynamics of surface waterbodies over large areas. The European Space Agency (ESA) has recently launched constellations of Sentinel-1 (S1) and Sentinel-2 (S2) satellites carrying C-band synthetic aperture radar (SAR) and a multispectral imaging radiometer, respectively. The constellations improve global coverage of remotely sensed imagery and enable the development of near real-time operational products. This unprecedented data availability leads to an urgent need for the application of fully automatic, feasible, and accurate retrieval methods for mapping and monitoring waterbodies. The mapping of waterbodies can take advantage of the synthesis of SAR and multispectral remote sensing data in order to increase classification accuracy. This study compares automatic thresholding to machine learning, when applied to delineate waterbodies with diverse spectral and spatial characteristics. Automatic thresholding was applied to near-concurrent normalized difference water index (NDWI) (generated from S2 optical imagery) and VH backscatter features (generated from S1 SAR data). Machine learning was applied to a comprehensive set of features derived from S1 and S2 data. During our field surveys, we observed that the waterbodies visited had different sizes and varying levels of turbidity, sedimentation, and eutrophication. Five machine learning algorithms (MLAs), namely decision tree (DT), k-nearest neighbour (k-NN), random forest (RF), and two implementations of the support vector machine (SVM) were considered. Several experiments were carried out to better understand the complexities involved in mapping spectrally and spatially complex waterbodies. It was found that the combination of multispectral indices with SAR data is highly beneficial for classifying complex waterbodies and that the proposed thresholding approach classified waterbodies with an overall classification accuracy of 89.3%. However, the varying concentrations of suspended sediments (turbidity), dissolved particles, and aquatic plants negatively affected the classification accuracies of the proposed method, whereas the MLAs (SVM in particular) were less sensitive to such variations. The main disadvantage of using MLAs for operational waterbody mapping is the requirement for suitable training samples, representing both water and non-water land covers. The dynamic nature of reservoirs (many reservoirs are depleted at least once a year) makes the re-use of training data unfeasible. The study found that aggregating (combining) the thresholding results of two SAR and multispectral features, namely the S1 VH polarisation and the S2 NDWI, respectively, provided better overall accuracies than when thresholding was applied to any of the individual features considered. The accuracies of this dual thresholding technique were comparable to those of machine learning and may thus offer a viable solution for automatic mapping of waterbodies.

@en

Adaptation to climate change is a major challenge facing the agricultural sector worldwide. Olive (Olea europaea L.) is a global, high value crop currently cultivated in 28 countries worldwide. Global data to assess the vulnerability of the crop to climate variability are scarce, and in some notable cases, such the United Nations Food and Agricutlure Organization database (FAO, 2006), qualitative assessments rather than quantitative indicators are provided. The aim of this study is to demonstrate a new approach to help overcome these constraints toward a globally applicable method to assess the adaptability of olive cultivars. The adaptability of 11 cultivars, widely used in 11 countries worldwide, was studied using a new generic approach based on the evaluation of soil hydrological regime against cultivar-specific hydrological requirements. The approach requires local data, notably on soil hydrological properties, but it is easily transferable to other countries and regions. We applied an agrohydrological model in 60 soil units to determine hydrological indicators both in a reference (1961–1990) and a future (2021–2050) climate case. We compared indicators with cultivar-specific requirements to achieve the target yield; requirements were established using experimental yield response curves. We estimated the probability of adaptation, i.e., the probability that a given cultivar attains the target yield, and we used it to evaluate the cultivar potential distribution in the study area. At the locations where soil hydrological conditions were favorable, the probabilities of adaptation of the cultivars were high in both climate cases. The results show that the area with suitable conditions for the target yield (area of adaptability) decreased under future climate for all the cultivars, with higher reduction for Frantoio and Maiatica and smaller reduction for Itrana, Nocellara, Ascolana, and Kalamata. These cultivars are currently grown in Argentina, United States (US), Australia, France, Greece, and Italy. Our results indicate also that these cultivars require higher available soil water to attain the target yield, i.e., we may expect similar vulnerability in other parts of the world. Based on these findings, we provide some specific recommendations for enrichment of global databases and for further developments of our approach, to increase its potential for global application.

@en

Forest fires are a major source of ecosystem disturbance. Vegetation reacts to meteorological factors contributing to fire danger by reducing stomatal conductance, thus leading to an increase of canopy temperature. The latter can be detected by remote sensing measurements in the thermal infrared as a deviation of observed land surface temperature (LST) from climatological values, that is as an LST anomaly. A relationship is thus expected between LST anomalies and forest fires burned area and duration. These two characteristics are indeed controlled by a large variety of both static and dynamic factors related to topography, land cover, climate, weather (including those affecting LST) and anthropic activity. To investigate the predicting capability of remote sensing measurements, rather than constructing a comprehensive model, it would be relevant to determine whether anomalies of LST affect the probability distributions of burned area and fire duration. This research approached the outlined knowledge gap through the analysis of a dataset of forest fires in Campania (Italy) covering years 2003-2011 against estimates of LST anomaly. An LST climatology was first computed from time series of daily Aqua-MODIS LST data (product MYD11A1, collection 6) over the longest available sequence of complete annual datasets (2003-2017), through the Harmonic Analysis of Time Series (HANTS) algorithm. HANTS was also used to create individual annual models of LST data, to minimize the effect of varying observation geometry and cloud contamination on LST estimates while retaining its seasonal variation. LST anomalies where thus quantified as the difference between LST annual models and LST climatology. Fire data were intersected with LST anomaly maps to associate each fire with the LST anomaly value observed at its position on the day previous to the event. Further to this step, the closest probability distribution function describing burned area and fire duration were identified against a selection of parametric models through the maximization of the Anderson-Darling goodness-of-fit. Parameters of the identified distributions conditional to LST anomaly where then determined along their confidence intervals. Results show that in the study area log-transformed burned area is described by a normal distribution, whereas log-transformed fire duration is closer to a generalized extreme value (GEV) distribution. The parameters of these distributions conditional to LST anomaly show clear trends with increasing LST anomaly; significance of this observation was verified through a likelihood ratio test. This confirmed that LST anomaly is a covariate of both burned area and fire duration. As a consequence, it was observed that conditional probabilities of extreme events appear to increase with increasing positive deviations of LST from its climatology values. This confirms the stated hypothesis that LST anomalies affect forest fires burned area and duration and highlights the informative content of time series of LST with respect to fire danger.@en

Reliable and timely mapping of crop growth conditions and of their water resources is considered a prioritary application in light of the abrupt climate changes. With this view, the paper presents a novel approach that makes use of dense C-Band time-series for the timely estimation of crop growth stages and for the detection of changes in crop water conditions, especially related to precipitation and irrigation events. Aided by vegetation indexes extracted from Landsat and Sentinel-2 imagery, the proposed Sentinel-l centered method exploits both temporal patterns of crop growth and spatial patterns of water anomalies to enhance its classification robustness.@en

Remote sensing imagery to monitor global biophysical dynamics requires the availability of thermal infrared data at high temporal and spatial resolution because of the rapid development of crops during the growing season and the fragmentation of most agricultural landscapes. Conversely, no single sensor meets these combined requirements. Data fusion approaches offer an alternative to exploit observations from multiple sensors, providing data sets with better properties. A novel spatio-temporal data fusion model based on constrained algorithms denoted as multisensor multiresolution technique (MMT) was developed and applied to generate TIR synthetic image data at both temporal and spatial high resolution. Firstly, an adaptive radiance model is applied based on spectral unmixing analysis of. TIR radiance data at TOA (top of atmosphere) collected by MODIS daily 1-km and Landsat - TIRS 16-day sampled at 30-m resolution are used to generate synthetic daily radiance images at TOA at 30-m spatial resolution. The next step consists of unmixing the 30 m (now lower resolution) images using the information about their pixel land-cover composition from co-registered images at higher spatial resolution. In our case study, TIR synthesized data were unmixed to the Sentinel 2 MSI with 10 m resolution. The constrained unmixing preserves all the available radiometric information of the 30 m images and involves the optimization of the number of landcover classes and the size of the moving window for spatial unmixing. Results are still being evaluated, with particular attention for the quality of the data streams required to apply our approach.

@en

In the context of climate change strategies are needed towards sustainable agricultural production. The aim of this study is to identify crop adaptation options to face the expected changes in water availability by exploiting the existing intra-specific biodiversity of the tomato crop and accounting for irrigation management and the hydrological properties of soils. The biophysical dimension of crop adaptation is therefore addressed. A study is presented examining an irrigated district in southern Italy. Using as a climatic reference the period 1961–90 and as a future climate the period 2021–2050, a soil water availability indicator was determined by a soil water balance model, at optimal irrigation and at different deficit irrigation strategies, in 23 soil units. For five tomato cultivars, hydrological requirements were determined by means of yield response functions to soil water availability. Cultivar-specific hydrological requirements were evaluated against the soil water availability indicator to determine probabilities of adaptation of each cultivar. These cultivars are not currently being grown in the study area so their potential spatial distribution in the study area was estimated. For instance, with 60% of optimal irrigation, two cultivars were assessed as having probabilities of crop adaptation larger than 0.89, in 90% and 62% of the area. In the future climate, with limited water resources, a proper choice and combination of cultivars, irrigation strategies and soils would allow to maintain the current production system in a large part of the study area.

@en

The Caprivi basin in Namibia has been affected by severe flooding in recent years resulting in deaths, displacements and destruction of infrastructure. The negative consequences of these floods have emphasized the need for timely, accurate and objective information about the extent and location of affected areas. Due to the high temporal variability of flood events, Earth Observation (EO) data at high revisit frequency is preferred for accurate flood monitoring. Currently, EO data has either high temporal or coarse spatial resolution. Accurate methodologies for the estimation and monitoring of flooding extent using coarse spatial resolution optical image data are needed in order to capture spatial details in heterogeneous areas such as Caprivi. The objective of this work was the retrieval of the fractional abundance of water (γ_w) by applying a new spectral indices-based unmixing algorithm to Medium Resolution Imaging Spectrometer Full Resolution (MERIS FR) data using a minimum number of spectral bands. These images are technically similar to the OLCI image data acquired by the Sentinel-3 satellite, which are to be systematically provided in the near future. The normalized difference wetness index (NDWI) was applied to delineate the water surface and combined with normalized difference vegetation index (NDVI) to account for emergent vegetation within the water bodies. The challenge to map flooded areas by applying spectral unmixing is the estimation of spectral endmembers, i.e., pure spectra of land cover features. In our study, we developed and applied a new unmixing method based on the use of an ensemble of spectral endmembers to capture and take into account spectral variability within each endmember. In our case study, forty realizations of the spectral endmembers gave a stable frequency distribution of γ_w. Quality of the flood map derived from the Envisat MERIS (MERIS) data was assessed against high (30 m) spatial resolution Landsat Thematic Mapper (TM) images on two different dates (17 April 2008 and 22 May 2009) during which floods occurred. The findings show that both the spatial and the frequency distribution of the γ_w extracted from the MERIS data were in good agreement with the high-resolution TM retrievals. The use of conventional linear unmixing, instead, applied using the entire available spectra for each image, resulted in relatively large differences between TM and MERIS retrievals.

@en

Water deficit limiting yields is one of the negative aspects of climate change. However, this applies particularly when emphasis is on biomass production (e.g. for field crops), but not necessarily for plants where quality, not quantity is most relevant. For grapevine development, mild water stress occurring during specific phenological phases is an important factor when producing good quality wines. It induces the production of anthocyanins and aroma precursors and then could offer an opportunity to increase winegrower's income.

A multidisciplinary study was carried out in Campania region (Southern Italy), an area well known for high quality wine production. Growth of Aglianico grapevine cultivar, with a standard clone population on 1103 Paulsen rootstocks, was studied on two different types of soil: Calcisols and Cambisols occurring along a slope of 90 m length with 11% gradient.

The agro-hydrological model SWAP was calibrated and applied to estimate soil-plant water status during three consecutive seasons (2011–2013). Crop water stress index (CWSI), as estimated by the model, was related to leaf water potential, sugar content of grape bunches and wine quality (e.g. content of tannins). For both soils, the correlations between quality measurements and CWSI were high (e.g. − 0.97** with sugar; 0.895* with anthocyanins in the grape skins).

The model was also applied to explore effects of future climate conditions (2021–2051) obtained from statistical downscaling of Global Circulation Models (AOGCM) and to estimate the effect of the climate on CWSI and hence on grape quality. Effects of climate change on grape quality indicate: (i) a resilient behavior of Calcisol to produce high quality wine, (ii) a good potentiality for improving the quality wine in Cambisol.

The present study represents an example of multidisciplinary approach in which soil scientists, hydro-pedologists, crop modellers, plant physiologists and oenologists have integrated their knowledge and skills in order to deal with the complex interactions among different components of an agricultural system.

@en

Satellite observations of the terrestrial biosphere cover a period of time sufficiently extended to allow the calculation of a reliable climatology. The latter is particularly relevant for studies of vegetation response to climate variability. Observations from space of the land surface are hampered by clouds at shorter wavelength and affected by water in the atmosphere in the microwave range. Both polar orbiting and geostationary satellites have a revisit frequency high enough to allow for some redundancy relative to the processes being observed, so that time series where a fraction of observations are removed and the resulting gaps filled are still very useful to monitor land surface processes. Two examples illustrate this concept in two different spectral regions: Thermal Infrared (TIR) and observations of land surface temperature to study the thermal behavior of the land surface in response to weather and climate and 37 GHz observations of the polarization difference in brightness temperature to retrieve the fractional abundance of watersaturated soil. Three applications of time series of land surface temperature are presented: (a) monitoring of spectral thermal admittance of the land surface; (b) estimation and mapping of air temperature and (c) monitoring of thermal load to assess the risk of forest fires. Two methods were applied to identify and remove anomalous observations (outliers) and to fill the resulting gaps: Harmonic ANalysis of Time Series (HANTS) and theMultichannel Singular Spectrum Analysis (M-SSA). The HANTS algorithm has been widely used to reconstruct time series of Normalized Difference Vegetation Index (NDVI), Leaf Area Index (LAI), Land Surface Temperature (LST) as well as the polarization difference brightness temperature (PDBT) during the past 20 years to remove random noise or eliminate cloud/snow contamination. The M-SSA, an advanced methodology for time series analysis, was utilized to reconstruct gap-free LST time series using both the spatial and the temporal information content in the data set.

@en

The appraisal of crop water requirements (CWR) is crucial for the management of water resources, especially in arid and semi-arid regions where irrigation represents the largest consumer of water, such as the Doukkala area, western Morocco. Simple and (semi) empirical approaches have been applied to estimate CWR: the first one is called Kc-NDVI method, based on the correlation between the Normalized Difference Vegetation Index (NDVI) and the crop coefficient (Kc); the second one is the analytical approach based on the direct application of the Penman-Monteith equation with reflectance-based estimates of canopy biophysical variables, such as surface albedo (r), leaf area index (LAI) and crop height (hc). A time series of high spatial resolution RapidEye (REIS), SPOT4 (HRVIR1) and Landsat 8 (OLI) images acquired during the 2012/2013 agricultural season has been used to assess the spatial and temporal variability of crop evapotranspiration ETc and biophysical variables. The validation using the dual crop coefficient approach (Kcb) showed that the satellite-based estimates of daily ETc were in good agreement with ground-based ETc, i.e., R2 = 0.75 and RMSE = 0.79 versus R2 = 0.73 and RMSE = 0.89 for the Kc-NDVI, respective of the analytical approach. The assessment of irrigation performance in terms of adequacy between water requirements and allocations showed that CWR were much larger than allocated surface water for the entire area, with this difference being small at the beginning of the growing season. Even smaller differences were observed between surface water allocations and Irrigation Water Requirements (IWR) throughout the irrigation season. Finally, surface water allocations were rather close to Net Irrigation Water Requirements (NIWR).@en