Satellite-based earth observations offer great opportunities to improve spatial model predictions by means of spatial-pattern-oriented model evaluations. In this study, observed spatial patterns of actual evapotranspiration (AET) are utilised for spatial model calibration tailored to target the pattern performance of the model. The proposed calibration framework combines temporally aggregated observed spatial patterns with a new spatial performance metric and a flexible spatial parameterisation scheme. The mesoscale hydrologic model (mHM) is used to simulate streamflow and AET and has been selected due to its soil parameter distribution approach based on pedo-transfer functions and the build in multi-scale parameter regionalisation. In addition two new spatial parameter distribution options have been incorporated in the model in order to increase the flexibility of root fraction coefficient and potential evapotranspiration correction parameterisations, based on soil type and vegetation density. These parameterisations are utilised as they are most relevant for simulated AET patterns from the hydrologic model. Due to the fundamental challenges encountered when evaluating spatial pattern performance using standard metrics, we developed a simple but highly discriminative spatial metric, i.e. one comprised of three easily interpretable components measuring co-location, variation and distribution of the spatial data. The study shows that with flexible spatial model parameterisation used in combination with the appropriate objective functions, the simulated spatial patterns of actual evapotranspiration become substantially more similar to the satellite-based estimates. Overall 26 parameters are identified for calibration through a sequential screening approach based on a combination of streamflow and spatial pattern metrics. The robustness of the calibrations is tested using an ensemble of nine calibrations based on different seed numbers using the shuffled complex evolution optimiser. The calibration results reveal a limited trade-off between streamflow dynamics and spatial patterns illustrating the benefit of combining separate observation types and objective functions. At the same time, the simulated spatial patterns of AET significantly improved when an objective function based on observed AET patterns and a novel spatial performance metric compared to traditional streamflow-only calibration were included. Since the overall water balance is usually a crucial goal in hydrologic modelling, spatial-pattern-oriented optimisation should always be accompanied by traditional discharge measurements. In such a multi-objective framework, the current study promotes the use of a novel bias-insensitive spatial pattern metric, which exploits the key information contained in the observed patterns while allowing the water balance to be informed by discharge observations.

In this study we evaluate a methodology for disaggregating land surface energy fluxes estimated with the Two-Source Energy Balance (TSEB)-based Dual-Temperature Difference (DTD) model which uses day and night polar orbiting satellite observations of land surface temperature (LST) as a remotely sensed input. The DTD model is run with MODIS input data at a spatial resolution of around 1 km while the disaggregation uses Landsat observations to produce fluxes at a nominal spatial resolution of 30 m. The higher-resolution modelled fluxes can be directly compared against eddy covariance (EC)-based flux tower measurements to ensure more accurate model validation and also provide a better visualization of the fluxes' spatial patterns in heterogeneous areas allowing for development of, for example, more efficient irrigation practices. The disaggregation technique is evaluated in an area covered by the Danish Hydrological Observatory (HOBE), in the west of the Jutland peninsula, and the modelled fluxes are compared against measurements from two flux towers: the first one in a heterogeneous agricultural landscape and the second one in a homogeneous conifer plantation. The results indicate that the coarse-resolution DTD fluxes disaggregated at Landsat scale have greatly improved accuracy as compared to high-resolution fluxes derived directly with Landsat data without the disaggregation. At the agricultural site the disaggregated fluxes display small bias and very high correlation (r ≈ 0.95) with EC-based measurements, while at the plantation site the results are encouraging but still with significant errors. In addition, we introduce a~modification to the DTD model by replacing the "parallel" configuration of the resistances to sensible heat exchange by the "series" configuration. The latter takes into account the in-canopy air temperature and substantially improves the accuracy of the DTD model.

This study evaluates three different metrics of water content of an herbaceous cover in a Mediterranean wooded grassland (dehesa) ecosystem. Fuel moisture content (FMC), equivalent water thickness (EWT) and canopy water content (CWC) were estimated from proximal sensing and MODIS satellite imagery. Dry matter (Dm) and leaf area index (LAI) connect the three metrics and were also analyzed. Metrics were derived from field sampling of grass cover within a 500 m MODIS pixel. Hand-held hyperspectral measurements and MODIS images were simultaneously acquired and predictive empirical models were parametrized. Two methods of estimating FMC and CWC using different field protocols were tested in order to evaluate the consistency of the metrics and the relationships with the predictive empirical models. In addition, radiative transfer models (RTM) were used to produce estimates of CWC and FMC, which were compared with the empirical ones. Results revealed that, for all metrics spatial variability was significantly lower than temporal. Thus we concluded that experimental design should prioritize sampling frequency rather than sample size. Dm variability was high which demonstrates that a constant annual Dm value should not be used to predict EWT from FMC as other previous studies did. Relative root mean square error (RRMSE) evaluated the performance of nine spectral indices to compute each variable. Visible Atmospherically Resistant Index (VARI) provided the lowest explicative power in all cases. For proximal sensing, Global Environment Monitoring Index (GEMI) showed higher statistical relationships both for FMC (RRMSE = 34.5 %) and EWT (RRMSE = 27.43 %) while Normalized Difference Infrared Index (NDII) and Global Vegetation Monitoring Index (GVMI) for CWC (RRMSE = 30.27 % and 31.58 % respectively). When MODIS data were used, results showed an increase in R2 and Enhanced Vegetation Index (EVI) as the best predictor for FMC (RRMSE = 33.81 %) and CWC (RRMSE = 27.56 %) and GEMI for EWT (RRMSE = 24.6 %). Differences in the viewing geometry of the platforms can explain these differences as the portion of vegetation observed by MODIS is larger than when using proximal sensing including the spectral response from scattered trees and its shadows. CWC was better predicted than the other two water content metrics, probably because CWC depends on LAI, that shows a notable seasonal variation in this ecosystem. Strong statistical relationship was found between empirical models using indices sensible to chlorophyll activity (NDVI or EVI which are not directly related to water content) due to the close relationship between LAI, water content and chlorophyll activity in grassland cover, which is not true for other types of vegetation such as forest or shrubs. The empirical methods tested outperformed FMC and CWC products based on radiative transfer model inversion.

Spatial patterns in long-term average evapotranspiration (ET) represent a unique source of information for evaluating the spatial pattern performance of distributed hydrological models on a river basin to continental scale. This kind of model evaluation is getting increased attention, acknowledging the shortcomings of traditional aggregated or timeseries-based evaluations. A variety of satellite remote sensing (RS)-based ET estimates exist, covering a range of methods and resolutions. There is, therefore, a need to evaluate these estimates, not only in terms of temporal performance and similarity, but also in terms of long-term spatial patterns. The current study evaluates four RS-ET estimates at moderate resolution with respect to spatial patterns in comparison to two alternative continental-scale gridded ET estimates (water-balance ET and Budyko). To increase comparability, an empirical correction factor between clear sky and all-weather ET, based on eddy covariance data, is derived, which could be suitable for simple corrections of clear sky estimates. Three RS-ET estimates (MODIS16, TSEB and PT-JPL) and the Budyko method generally display similar spatial patterns both across the European domain (mean SPAEF = 0.41, range 0.25–0.61) and within river basins (mean SPAEF range 0.19–0.38), although the pattern similarity within river basins varies significantly across basins. In contrast, the WB-ET and PML_V2 produced very different spatial patterns. The similarity between different methods ranging over different combinations of water, energy, vegetation and land surface temperature constraints suggests that robust spatial patterns of ET can be achieved by combining several methods.

Distributed hydrological models simulate states and fluxes of water and energy in the terrestrial hydrosphere at each cell. The predicted spatial patterns result from complex nonlinear relationships and feedbacks. Spatial patterns are often neglected during the modeling process, and therefore a spatial sensitivity analysis framework that highlights their importance is proposed. This study features a comprehensive analysis of spatial patterns of actual evapotranspiration (ET) and land surface temperature (LST), with the aim of quantifying the extent to which forcing data and model parameters drive these patterns. This framework is applied on a distributed model [MIKE Système Hydrologique Européen (MIKE SHE)] coupled to a land surface model [Shuttleworth and Wallace–Evapotranspiration (SW-ET)] of a catchment in Denmark. Twenty-two scenarios are defined, each having a simplified representation of a potential driver of spatial variability. A baseline model that incorporates full spatial detail is used to assess sensitivity. High sensitivity can be attested in scenarios where the simulated spatial patterns differ significantly from the baseline. The core novelty of this study is that the analysis is based on a set of innovative spatial performance metrics that enable a reliable spatial pattern comparison. Overall, LST is very sensitive to air temperature and wind speed whereas ET is rather driven by vegetation. Both are sensitive to groundwater coupling and precipitation. The conclusions may be limited to the selected catchment and to the applied modeling system, but the suggested framework is generically relevant for the modeling community. While the applied metrics focus on specific spatial information, they partly exhibit redundant information. Thus, a combination of metrics is the ideal approach to evaluate spatial patterns in models outputs.