Comment on hess-2022-125

Evapotranspiration (ET) is dynamic and influences water resource distribution. Sustainable management of water resources requires accurate estimations of the individual components that result in evapotranspiration, including the daily net radiation (DNR). Daily ET is more useful than the evaporative fraction (EF) provided by remote sensing ET models, and to account for daily variations, EF is usually combined with the DNR. DNR exhibits diurnal and spatiotemporal variations due to landscape heterogeneity. In the modified Two-Source Energy Balance (TSEB) approach by Zhuang and Wu, 2015, ecophysiological constraint functions of temperature and moisture of plants based on atmospheric moisture and vegetation indices were introduced, but the DNR was not spatially accounted for in the estimation of the daily ET. This research adopted a novel approach that accounts for spatiotemporal variations in estimated daily ET by incorporating the Bisht and Bras DNR model in the modified version of the TSEB model. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite imagery over the Tono irrigation watershed within the Upper East Region of Ghana and Southern Burkina Faso were used. We estimated the energy fluxes of latent and sensible heat as well as the net radiation and soil heat fluxes from the satellite images and compared our results with ground-based measurements from an eddy covariance (EC) station established by the West African Science Service Center on Climate Change and Adapted Land Use (WASCAL) within the watershed area. We noticed a similarity between our model estimated fluxes and ET with the ground-based EC station measurements. Eight different land use/cover types were identified in the study area, and each of these contributed significantly to the overall ET variations between the two study periods: December 2009 and December 2017. For instance, due to a higher leaf area index (LAI) for all vegetation types in December 2009 than in December 2017, the ET for December 2017 was higher than that for December 2009. We also noticed that the land use/cover types within the footprint area of the EC station were only six out of the eight. Generally, all the surface energy fluxes increased from December 2009 to December 2017. Mean ET varied from 3.576 to 4.486 (mm/d) for December 2009 while from 4.502 to 5.280 (mm/d) for December 2017 across the different land use/cover classes. Knowledge of the dynamics of evapotranspiration and adoption of cost-effective methods to estimate its individual components in an effective and efficient way is critical to water resources management. Our findings provide a tool for all water stakeholders within watersheds to manage water resources in an engaging and cost-effective way.

Multi-scale evaluation of global evapotranspiration products derived from remote sensing images: Accuracy and uncertainty

An efficient calibration with remotely sensed (RS) data is important for accurate predictions at ungauged catchments. This study investigates the advantages of streamflow-sensitive regionalization on calibration with RS evapotranspiration (ET). Regionalization experiments are performed at 28 catchments in Australia. The catchments are classified into three groups based on annual rainfall and runoff coefficients. Streamflow, RS ET, and a multi-objective RS ET-streamflow calibration are performed using the DiffeRential Evolution Adaptive Metropolis algorithm in each catchment. Simplified Australian Water Resource Assessment-Landscape model is calibrated for a selection of five parameters. Posterior probability distributions of parameters from three calibrations performed at donor catchments in each group are inspected to find the parameter for regionalization in the individual group. In group 1 of wetter catchments, regionalization of parameter FsoilEmax (soil evaporation scaling factor) helps to simplify the calibration without any deterioration in ET, soil moisture (SM) and streamflow predictions. Regionalization of parameter Beta (coefficient describing rate of hydraulic conductivity increase with water content) in group 2 assists to improve the streamflow predictions with no decrement in ET and SM predictions. However, regionalization is not able to provide satisfactory results in group 3. Group 3 includes low-yielding catchments, with average annual rainfall below 1000 mm/year and runoff coefficient less than 0.1, where traditional streamflow calibration also fails to produce accurate results. This study concludes that streamflow-sensitive regionalization is effective for improving the efficacy of RS ET calibration in wetter catchments.

One of the major production areas in Europe for Hardy Ornamental Nursery Stocks (HONS) is located around Pistoia (Tuscany, Italy). In this area, pot ornamental crops show low water use efficiency (WUE). Main reasons of poor irrigation efficiency are both the inaccurate irrigation scheduling, which is normally based on grower’s personal experience, and the choice to group different plant species together in the same irrigation plot in order to simplify irrigation and to reduce growing costs. An efficiency-oriented approach for the management of HONS requires the use of soil/substrate moisture sensors (SMS) and/or crop evapotranspiration (ET) model to regulate irrigation frequency. In order to compare different methods for irrigation scheduling in container beds hosting four different HONS species, two experiments were carried out in the spring-summer season of 2009 and 2010 at ‘Centro Speri-mentale per il Vivaismo’ (CESPEVI), Pistoia. In particular, we tested two irrigation control methods based on SMS or ET model as compared to conventional timer scheduling. The WET sensor was used for SMS control while the ET model consisted of a composite regression model of crop ET (as function of plant height and potential evapotranspiration) along with a specific algorithm implemented in the software of an irrigation controller prototype. Scheduling irrigation by means of the WET sensor (irrigation was activated when the matric potential of substrate was less than -5.0 kPa) or applying ET model reduced seasonal water use respectively by 40% and 45%, and water runoff by 74% and 81% in comparison with timer control strategy, without significant effects on plant ET, growth, and commercial value at the end of growing season. Both new irrigation scheduling are useful for growers, but for a large use of RZS method, a low cost RZS sensor is required. The use ET-model is less expensive, but it requires the knowledge of the Kc for the hundred of different species grown in the nurseries.

Evapotranspiration (ET) is a significant consumer of irrigation water and precipitation on cropland. Global and regional interest in the sustainable management of limited freshwater supplies to meet the rapidly increasing population and food demands has resulted in advanced scientific research on ET measurement, rapid water accounting, and irrigation schedules in the NENA region. The primary goal of this paper is to compare actual daily evapotranspiration (ET) collected by a remote sensing model and validated by Energy Balance (EB) flux tower field measurements. The flux tower was installed in a wheat field in Sids Agricultural Research Station in Beni Suef Governorate. Through the integration of Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Sentinel-2 data, a new remote sensing-based ET model is built on two parties: Thermal condition factor (TCF) and vegetation condition fraction (VCF). The remote sensing-based ET estimation model was evaluated using ET field measurements from the Energy Balance flux tower. The land use and land cover maps were created to assist the interpretation of remotely sensed ET data. Field data for five categories were collected to test the accuracy of the land use and cover maps: Water bodies (93 points), urban areas (252 points), trees (104 points), other field crops (227 points), and wheat (249 points), for a total of 925 ground points. The Google Earth Engine (GEE) imported sentinel-2 datasets and filtered the necessary dates and regions. From 1 October 2020 to 30 May 2021, sentinel-2 data were processed and transformed into the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), and Normalized Difference Built-up Index (NDBI), which were then combined. The composite layer data were classified using the Random Forest (RF) method on the GEE platform, and the results showed an overall accuracy of 91 percent. The validation factors revealed good indices when RS-based ET results were compared to ground-measured ET. The Root Mean Square Error (RMSE) was 0.84 mm/day. The ‘r’ and ‘d’ values indicated satisfactory results, where ‘r’ yielded a value of 0.785, which indicates that the correlation between predicted and reference results is robust. The analysis of d values revealed a high degree of correlation between predicted (RS-based ET) and reference results (measured ET). The d value was found to be 0.872. Between 21 November 2020 and 30 April 2021, RS-based accumulated ET was 418 mm/season, while ground-measured ET was 376 mm/season. The new RS-based ET model produced acceptable daily and seasonal results.

Development of a Benchmark Eddy Flux Evapotranspiration Dataset for Evaluation of Satellite-Driven Evapotranspiration Models Over the CONUS

Large sample datasets of in situ evapotranspiration (ET) measurements with well documented data provenance and quality assurance are critical for water management and many fields of earth science research. We present a post-processed ET oriented dataset at daily and monthly timesteps, from 161 stations, including 148 eddy covariance flux towers, that were chosen based on their data quality from nearly 350 stations across the contiguous United States. In addition to ET, the data includes energy and heat fluxes, meteorological measurements, and reference ET downloaded from gridMET for each flux station. Data processing techniques were conducted in a reproducible manner using open-source software. Most data initially came from the public AmeriFlux network, however, several different networks (e.g., the USDA-Agricultural Research Service) and university partners provided data that was not yet public. Initial half-hourly energy balance data were gap-filled and aggregated to daily frequency, and turbulent fluxes were corrected for energy balance closure error using the FLUXNET2015/ONEFlux energy balance ratio approach. Metadata, diagnostics of energy balance, and interactive graphs of time series data are included for each station. Although the dataset was developed primarily to benchmark satellite-based remote sensing ET models of the OpenET initiative, there are many other potential uses, such as validation for a range of regional hydrologic and atmospheric models.

<strong class="journal-contentHeaderColor">Abstract.</strong> Remotely sensed evapotranspiration (RS-ET) products have been widely adopted as additional constraints on hydrologic modeling to enhance the model predictability while reducing predictive uncertainty. However, vegetation parameters, responsible for key time/space variation in evapotranspiration (ET), are often calibrated without the use of suitable constraints. Remotely sensed leaf area index (RS-LAI) products are increasingly available and provide an opportunity to assess vegetation dynamics and improve calibration of associated parameters. The goal of this study is to assess the Soil and Water Assessment Tool (SWAT) predictive uncertainty in estimates of ET using streamflow and two remotely sensed products (i.e., RS-ET and RS-LAI). We explore how the application of RS-ET and RS-LAI products contributes to 1) reducing the parameter uncertainty; 2) improving the model capacity to predict the spatial distribution of ET and LAI at the sub-watershed level; and 3) assessing the model predictions of ET and LAI at the basic modeling unit (i.e., the hydrologic response unit [HRU]) under the assumption that ET and LAI are related in croplands. Our results suggest that most of the parameter sets with acceptable performances for two constraints (i.e., streamflow and RS-ET; 12 parameter sets) are also acceptable for three constraints (i.e., streamflow, RS-ET, and RS-LAI; 11 parameter sets) at the watershed level. This finding is likely because both the ET simulation algorithm and the RS-ET products consider LAI as an input variable. Relative to the watershed-level assessment, the number of parameter sets that satisfactorily characterize spatial patterns of ET and LAI at the sub-watershed level are reduced from 11 to 6. Among the 11 parameter sets acceptable for three constraints (i.e., streamflow, RS-ET and RS-LAI) at the sub-watershed level, two parameter sets appear to provide high spatial and temporal consistency between ET and LAI at the HRU level. These results suggested that use of multiple remotely sensed products as model constraints enables model evaluations at finer scales &ndash; thereby constraining acceptable parameter sets and accurately representing the spatial characteristics of hydrologic variables. As such, this study highlights the potential of remotely sensed data to increase the predictability and utility of hydrologic models.

A sensitivity and uncertainty analysis of the Penman‐Brutsaert evapotranspiration (ET) model was conducted using meteorological data collected from the humid South Florida region. Since net radiation and relative humidity were found to be correlated (ρ = −0.74) in this region, a method was developed to analyze the sensitivity with the correlation effect of these two independent variables. After conducting sensitivity analyses with and without the correlation effect, the results were evaluated using conditional probability density functions. Both theoretical and computed results show that the conditional probability with the correlation effect increases in proportion to the increasing absolute correlation coefficient of two independent variables. This finding suggests that the sensitivity analysis with the correlation effect represents more adequately a given set of meteorological data in South Florida than that without the correlation effect. After defining the random errors of the independent variables based on multiple measurements, the ET model error as well as the propagated error caused by erroneous independent variables are presented.

Model calibration and validation are challenging in poorly gauged basins. We developed and applied a new approach to calibrate hydrological models using distributed geospatial remote sensing data. The Soil and Water Assessment Tool (SWAT) model was calibrated using only twelve months of remote sensing data on actual evapotranspiration (ETa) geospatially distributed in the 37 sub-basins of the Lake Chad Basin in Africa. Global sensitivity analysis was conducted to identify influential model parameters by applying the Sequential Uncertainty Fitting Algorithm–version 2 (SUFI-2), included in the SWAT-Calibration and Uncertainty Program (SWAT-CUP). This procedure is designed to deal with spatially variable parameters and estimates either multiplicative or additive corrections applicable to the entire model domain, which limits the number of unknowns while preserving spatial variability. The sensitivity analysis led us to identify fifteen influential parameters, which were selected for calibration. The optimized parameters gave the best model performance on the basis of the high Nash–Sutcliffe Efficiency (NSE), Kling–Gupta Efficiency (KGE), and determination coefficient (R2). Four sets of remote sensing ETa data products were applied in model calibration, i.e., ETMonitor, GLEAM, SSEBop, and WaPOR. Overall, the new approach of using remote sensing ETa for a limited period of time was robust and gave a very good performance, with R2 &gt; 0.9, NSE &gt; 0.8, and KGE &gt; 0.75 applying to the SWAT ETa vs. the ETMonitor ETa and GLEAM ETa. The ETMonitor ETa was finally adopted for further model applications. The calibrated SWAT model was then validated during 2010–2015 against remote sensing data on total water storage change (TWSC) with acceptable performance, i.e., R2 = 0.57 and NSE = 0.55, and remote sensing soil moisture data with R2 and NSE greater than 0.85.

Evapotranspiration (ET) links the water and carbon cycles in the atmosphere, hydrosphere, and biosphere. In this study, we develop an ET modelling framework based on the idea that the transpiration and carbon uptake are closely coupled, as predicted by the ‘least-cost hypothesis’ that canopy conductance acclimates to environmental variations. According to eco-evolutionary optimality theory, which has been previously applied in monitoring and modelling land-surface processes, the total costs (per unit carbon fixed) for maintaining transpiration and carboxylation capacities should be minimized. We calculate gross primary production (GPP) assuming that the light- and Rubisco-limited rates of photosynthesis, described by the classical biochemical model of photosynthesis, are coordinated on an approximately weekly time scale. Transpiration (T) is then calculated via acclimated canopy conductance, with no need for plant type- or biome-specific parameters. ET is finally calculated from T using an empirical function of light, temperature, soil water content and foliage cover to predict the T/ET ratio at each site. The GPP estimates were well supported by (weekly) GPP data at 20 widely distributed eddy-covariance flux sites (228 site-years), with correlation coefficients (r) = 0.81 and root-mean-square error (RMSE) = 18.7 gC week−1 and Nash-Sutcliffe efficiency (NSE) = 0.61. Predicted ET was also well supported, with r =0.85, RMSE = 5.5 mm week–1 and NSE = 0.66. Estimated T/ET ratios (0.43–0.74) showed significant positive relationships to radiation, precipitation and green vegetation cover and negative relationships to temperature and modelled T (r = 0.84). Aspects of this framework could be improved, notably the estimation of T/ET. Nonetheless, we see the application of eco-evolutionary principles as a promising direction for water resources research, eliminating the uncertainty introduced by the need to specify multiple parameters, and leveraging the power of remotely sensed vegetation cover data as a key indicator of ecosystem function.

Abstract. We present a model to retrieve actual evapotranspiration (ET) from satellites' vegetation indices (Parameterization of Vegetation Indices for ET estimation model, or PaVI-E) for the eastern Mediterranean (EM) at a spatial resolution of 250 m. The model is based on the empirical relationship between satellites' vegetation indices (normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) from MODIS) and total annual ET (ETAnnual) estimated at 16 FLUXNET sites, representing a wide range of plant functional types and ETAnnual. Empirical relationships were first examined separately for (a) annual vegetation systems (i.e. croplands and grasslands) and (b) systems with combined annual and perennial vegetation (i.e. woodlands, forests, savannah and shrublands). Vegetation indices explained most of the variance in ETAnnual in those systems (71 % for annuals, and 88 % for combined annual and perennial systems), while adding land surface temperature data in a multiple-variable regression and a modified version of the Temperature and Greenness model did not result in better correlations (p &gt; 0.1). After establishing empirical relationships, PaVI-E was used to retrieve ETAnnual for the EM from 2000 to 2014. Models' estimates were highly correlated (R = 0.92, p &lt; 0.01) with ETAnnual calculated from water catchment balances along rainfall gradient of the EM. They were also comparable to the coarser-resolution ET products of the Land Surface Analysis Satellite Applications Facility (LSA-SAF MSG ETa, 3.1 km) and MODIS (MOD16, 1 km) at 148 EM basins with R of 0.75 and 0.77 and relative biases of 5.2 and −5.2 %, respectively (p &lt; 0.001 for both). In the absence of high-resolution (&lt; 1 km) ET models for the EM the proposed model is expected to contribute to the hydrological study of this region, assisting in water resource management, which is one of the most valuable resources of this region.

Abstract. We present an evapotranspiration (ET) model developed in the framework of the EUMETSAT "Satellite Application Facility" (SAF) on Land Surface Analysis (LSA). The model is a simplified Soil-Vegetation-Atmosphere Transfer (SVAT) scheme that uses as input a combination of remote sensed data and atmospheric model outputs. The inputs based on remote sensing are LSA-SAF products: the Albedo (AL), the Downwelling Surface Shortwave Flux (DSSF) and the Downwelling Surface Longwave Flux (DSLF). They are available with the spatial resolution of the MSG SEVIRI instrument. ET maps covering the whole MSG field of view are produced from the model every 30 min, in near-real-time, for all weather conditions. This paper presents the adopted methodology and a set of validation results. The model quality is evaluated in two ways. First, ET results are compared with ground observations (from CarboEurope and national weather services), for different land cover types, over a full vegetation cycle in the Northern Hemisphere in 2007. This validation shows that the model is able to reproduce the observed ET temporal evolution from the diurnal to annual time scales for the temperate climate zones: the mean bias is less than 0.02 mm h−1 and the root-mean square error is between 0.06 and 0.10 mm h−1. Then, ET model outputs are compared with those from the European Centre for Medium-Range Weather Forecasts (ECMWF) and the Global Land Data Assimilation System (GLDAS). From this comparison, a high spatial correlation is noted, between 80 to 90%, around midday. Nevertheless, some discrepancies are also observed and are due to the different input variables and parameterisations used.

Abstract. It remains unclear how the selection of a spatial domain affects the accuracy of evapotranspiration (ET) estimates from contextual remote sensing based surface energy balance (SEB) models, particularly at large spatial scales. We thus tested the effect of spatial domain on four widely implemented contextual remote sensing based SEB models: Surface Energy Balance Algorithm for Land (SEBAL), Mapping ET with Internalized Calibration (METRIC), Simplified Surface Energy Balance Index (S-SEBI), and Triangular ET models. We applied these models on 44 near cloud-free Moderate Resolution Imaging Spectroradiometer (MODIS) thermal images across all of India from 2004 to 2006. Four spatial domains were considered: all of India, agro-ecological regions, 300 km × 300 km grids, and 600 km × 600 km grids and we compared hourly ET estimates from the models against observed ET data at four Bowen Ratio sites in India. Model performance varied across all models and spatial domains. ET values of neighboring pixels across spatial domains formed sharp edges along the boundaries of agro-ecological regions, 300 km × 300 km grids, and 600 km × 600 km grids suggesting that all ET models are highly sensitive to the selection of spatial domain. No single spatial domain was found to be optimal for all models and hence potential uncertainties associated with the selection of spatial domain should be taken into consideration when implementing these models at regional scales. The results from this study provide guidance for future regional-scale implementation of ET models and potential approaches to overcome these challenges.

Core Ideas Increased bias in ET due to various data sources where generally utilized standard weather data poses similar biases. Increased understanding of the complex relationship among the various weather input and vegetation indices. Results showed local weather station data closely represented site condition indicating the difficulties and challenges in ET modeling using other commonly available and widely adopted data sources. ABSTRACTComplex interactions among meteorological data and vegetation indices are incompletely understood in relation to evapotranspiration (ET) calculations for larger spatial domains with higher spatial and temporal resolution. Objectives of this study were to evaluate contributions of inputs to uncertainty in ET calculations and to enhance understanding of interactions among weather data, vegetative indices, and resistances utilized in biophysical ET model. We evaluated individual and combined effects of weather variables and vegetation indices using BAITSSS (Backward‐Averaged Iterative Two‐Source Surface temperature and energy balance Solution). Local weather station (LWS) data at a lysimeter site were obtained for irrigated corn (Zea mays L.) during the growing season (May to September, 2016) at Bushland, Texas. Gridded meteorological data were obtained from North American Land Data Assimilation System (NLDAS) (∼ 12.5 km) and remotely‐sensed vegetation indices (Landsat 30 m). Standard weather station (SWS) data were obtained from a grass reference near lysimeter site. The r2 and RMSE of ET simulated using LWS data and measured vegetation indices were 0.90 and 0.85 mm for daily ET, and 0.90 and 0.10 mm, for hourly ET, compared to lysimeter ET (less than 4% cumulative error). However, r2 and RMSE were 0.74 and 1.64 mm for daily ET, and 0.81 and 0.14 mm for hourly ET using gridded data, with positive bias (∼ 25% from NLDAS data). Simulated ET from SWS data exhibited similar behavior to gridded data with increased ET up to 21%. Results quantify difficulties in ET modeling using commonly available and widely adopted data sources.