Evaluation of SEBS and SEBAL algorithms for estimating wheat evapotranspiration (case study: central areas of Khuzestan province)

This study aimed to accurately estimate daily wheat evapotranspiration using two remote sensing algorithms, Surface Energy Balance System (SEBS) and Surface Energy Balance Algorithm for Land (SEBAL), in central Khuzestan province during 2019–2020. The results of two algorithms were compared with lysimeter (as a direct method), FAO-Penman–Monteith (FAO-PM), two temperature-based methods (Hargreaves-Samani and Blaney-Criddle), two radiation-based methods (Priestley–Taylor and Doorenbos–Pruitt), and two mass transfer-based methods (Mahringer and World Meteorology Organization) (as indirect methods). Coefficient of Determination (R2), Root-Mean-Square Error (RMSE), Percentage of Bias (PBIAS), Mean Bias Error, Mean Absolute Percentage Error, and Nash–Sutcliffe indicators used for comparing the results. According to the results, both SEBAL and SEBS algorithms showed the highest compatibility with lysimeter data (R2 = 0.92 and 0.96, RMSE = 2.15 and 1.53 mm/day, respectively). Comparing both algorithms with the FAO-PM method, resulted in RMSE and R2 of 2.42 mm/day and 0.87 for SEBS and 3.14 mm/day and 0.79 for SEBAL. The Hargreaves-Samani method (R2 = 0.72, RMSE = 16.4 mm/day) and (R2 = 0.8, RMSE = 10.4 mm/day) among temperature-based methods, Doorenbos–Pruitt (R2 = 0.71, RMSE = 3.33 mm/day) and (R2 = 0.79, RMSE = 2.63 mm/day) among radiation-based methods, and the Mahringer method (R2 = 0.6, RMSE = 6.8 mm/day mm/day) and (R2 = 0.68, RMSE = 5.51 mm/day) among mass transfer-based methods yielded better estimations than SEBAL and SEBS algorithms, respectively. Owing to the high accuracy of SEBAL and SEBS algorithms, in estimating the amount of evapotranspiration in the study area and close to the actual values in the field, using energy balance algorithms is recommended in Khuzestan province.


Introduction
Crop evapotranspiration (ET) plays a key role in determining crop water requirements and ultimately correctly designing irrigation systems (Djaman et al. 2015). ET is a major component of water balance, and its accurate estimation is of great importance in water use planning and optimization (Asadi et al. 2022). It is also necessary for irrigation and water resources management, increasing yield, and better crop management (Shamloo et al. 2021). ET is correctly estimated using different methods, the accuracy of which depends on the climatic conditions of the study areas (Racz et al. 2013). The actual ET can be measured directly using the water balance method, which is often a costly and complicated technique, but it is considered an effective tool for the validation and calibration of ET estimation models. Although indirect methods, such as mathematical models for ET estimation, can be employed easily, they are used when measurement methods are difficult (Obada et al. 2017). Based on the literature, conventional ET estimation methods are classified into three main groups, namely temperaturebased, radiation-based, and mass transfer-based methods (Xu and Singh. 2002). In the FAO-Penman-Monteith (FAO-PM) method, the weather parameters are considered related to these three categories (Obada et al. 2017). Climatic data are limited in most areas, and it is not the ability to use the FAO-PM method; thus, other methods are recommended for ET calculation (Djaman et al. 2015). Remote sensing makes it possible to obtain the amount of daily ET in different areas, with wide scales at the least possible time and with high economic benefit (Saboori et al. 2021).
Most RS-based Surface Energy Balance algorithms have been developed for crop ET determination in specific conditions in terms of land use and crop management, considering the accordance of empirical functions and parameters with these algorithms (Wolff et al. 2022). They are easily used in accurate agricultural management systems for better decision-making and higher yields (Shamloo et al. 2021).
To determine actual crop ET, many temperature-based, radiation-based, and mass transfer-based empirical methods were compared with FAO-PM and lysimeter methods in different climatic and spatial conditions. This indicates the importance of crop ET estimation, including cereals, to determine crop water requirements and the optimal use of water resources. In this context, Djaman et al. assessed 16 ET methods under coastal conditions in the Senegal River Delta, and the results indicated the good performance of the Mahringer and Terabert methods (Djaman et al. 2015). Lang et al. compared eight ET estimation methods with the FAO-PM method in southwest China. The results revealed that the performance of these methods was dependent on the regional climate type in each area. However, the best performance was reported for Makkink and Hargreaves-Samani methods from radiation and temperature-based methods, respectively (Lang et al. 2017). Zoratipour et al. studied the spatial and temporal evaluation of different methods for the prediction of ET in Khuzestan province. According to the results, the best ET estimations in this province were recorded for the Hargreaves-Samani among the temperature-based method, Doorenbos-Pruitt among the radiation-based method, and Mahringer and the World Meteorological Organization (WMO) among the mass transfer-based methods (Zoratipour et al. 2019). Empirical methods have limitations due to the need for measuring all effective parameters at more time and higher cost for the necessary equipment preparation and agricultural operations; also it is not easy to measure soil water balance data at depths. Given the limitations of various empirical methods, the possibility of errors in field measurements, and their non-allocation to large extents, the results can be feasibly extended with a combination of land and RS data as a Surface Energy Balance model. Accordingly, Lian et al. estimated (Tan et al. 2021). Khand et al. assessed wheat ET modeling based on SEBAL, METRIC, and SEBS algorithms and obtained a smaller RMSE (0.14 mm/day) for the METRIC model than those of SEBAL and SEBS algorithms (Khand et al. 2021). Shamloo et al. evaluated the SEBS algorithm to estimate maize ET and crop coefficient using Landsat 8 images in the Adana Mediterranean Area, Turkey. According to the results, the SEBAL estimated ET values mostly corresponded to the FAO-PM method with R 2 = 0.91 and RMSE = 1.14 mm/day. It was also highly correlated with Turc, Hargreaves, and Makkink methods (Shamloo et al. 2021). Asadi et al. compared actual wheat crop ET based on SEBAL algorithm using 12 Landsat 7 and Landsat 8 images during the crop development period in the Parsabad Moghan Plain, Northwestern Iran. The results demonstrated that the SEBAL algorithm (RSME = 0.633 mm/day and R 2 = 0.93) had the minimum error rate and maximum similarity compared to lysimeter data (Asadi and Valizadeh Kamran. 2022). Yang et al. investigated actual ET in different land uses based on the SEBAL algorithm and landsat 8 images in Ecotone, Northwestern China. The SEBAL algorithm had an appropriate estimation (RMSE = 0.9 mm/day and R 2 = 0.81), which was suitable for research on water resources, but it overestimated ET . In another study, Wei et al. evaluated daily ET estimation of rice fields, using the SEBAL algorithm in a subtropical region in Southern China. The results showed a high precision of the SEBAL algorithm, on a daily scale, with R 2 , NSE, and RMSE values of 0.85, 0.81, and 0.84 mm/day, which confirmed SEBAL application to logically allocate water resources in subtropical regions (Wei et al. 2022). Tariqul Islam et al. determined actual wheat evapotranspiration by the lysimeter method and SEBAL algorithm using Landsat 8 and Sentinel 2 satellite images. The results concluded the ET of wheat estimated by SEBAL had maximum similarity compared to lysimeter data. Also the average seasonal wheat ET was calculated to be 253 mm (Tariqul Islam et al. 2023). Liu et al. used the SEBS algorithm and Landsat 8 satellite images to estimate regional evapotranspiration in Beijing. The results showed that ET decreased with increasing land surface temperature (LST) and daily evapotranspiration ranged from 3.47 to 5.47 mm (Liu et al. 2023). Khoshnood et al. evaluated SEBS algorithm to investigate the rate of evapotranspiration change in different land cover use classes in the entire Urmia lake basin from 2016 to 2020. They stated that there is a high correlation between the results of the SEBS algorithm and the values measured (Khoshnood et al. 2023).
ET evaluation in every region and choosing the best methods tailored to each weather type contribute to the optimal use of water resources. With extensive crop cultivation, Khuzestan province is the central hub of agriculture, cereal production, and ultimately, its distributor to other provincial regions and Iran. The accurate ET estimation using RS can result in irrigation management, wheat water requirement supply, and optimal water use in a short temporal interval with a minimum economic cost. The climate of Khuzestan province is dry and hyperarid with annual evapotranspiration above 3000 mm and annual rainfall of 200 mm. Also, wheat is one of the most important strategic and widely consumed products (with extensive crop cultivation about 700-800 thousand hectares) in Iran, especially in Khuzestan province. Therefore, the evaluation of wheat evapotranspiration is necessary based on remote sensing algorithms and comparison with actual and empirical methods. Also, the evaluated results of the present algorithms in different researches are not the same and in many studies, the SEBAL algorithm has better results among the energy balance algorithms and in a number of others, the SEBS algorithm. In addition, in this research, remote sensing algorithms are investigated in comparison with various empirical methods categorized in the form of temperature, radiation, and mass transfer methods. Thus, it is essential to estimate the wheat crop ET based on SEBS and SEBAL algorithms compared to actual and empirical methods. Also, this topic has not been evaluated on the wheat crop in the study region, sofar. Therefore, this work intends to: • Estimating the wheat crop ET using SEBS and SEBAL algorithms. • Comparing the ET estimated using SEBS and SEBAL algorithms with the lysimeter method.

Study area
Khuzestan province is located at 47° 41′ to 50° 39′ E from the prime meridian and 29° 58′ to 33° 4′ N from the equator in the southwest of Iran, with an area of about 64,057 km 2 . The province shares borders with the Persian Gulf and Iraq in the south and west, respectively. In this study, daily weather data (2019-2020) were collected for the Shahid Modarres Basin in Ahvaz (Khuzestan province). The study area and the specifications of the studied station are represented in Fig. 1 and Table 1, respectively. The soil physicochemical and irrigation water properties are also listed in Tables 2 and 3, respectively.

Data calculation
In this research, weather data were obtained from the Emam Khomeini agro industry Synoptic Station to implement the SEBS and SEBAL algorithms and evaluate the obtained results in comparison with other methods. Minimum and maximum air temperatures (°C), dew point temperature (°C), average relative humidity (%), maximum and minimum relative humidity (%), wind speed at the height of 2 m (m/s), number and maximum hours of sunshine, precipitation, and air pressure (Pa) were the data used in this study. These data were prepared from the mentioned synoptic station in Khuzestan province according to satellite passing days. The results of SEBS and SEBAL algorithms were compared by lysimeter and conventional empirical methods for ET estimation, including two temperature-based, two radiation-based, two mass transfer-based, and FAO-PM methods. Due to little incomplete (missing) data, they were estimated using appropriate renewal techniques (regression relative to the situation of each adjacent stations) and were controlled qualitatively before use. The ET obtained from SEBS and SEBAL algorithms were compared with all the mentioned methods using statistical indices to introduce the best method. To this aim, six Landsat 8 OLI satellite images (2019-2020) obtained from the US Geological Survey were used during the wheat growth period. It was attempted to select images with good cloudless weather conditions to provide ET comparisons. The accurate dates of the images used in this study are presented in Table 4. To estimate actual ET, cloudless Landsat 8 satellite images were used during the wheat growth period    Ca (meq/l) 28 Mg (meq/l) 2.15 Na (meq/l) (November. 27 to May. 5 from 2019 to 2020). Since there are numerous errors in Landsat 8 images, the errors were eliminated by using two types of atmospheric and radiometric correction. The characteristics of the used images are shown in Table 5.

Surface Energy Balance Algorithm for Land (SEBAL)
The SEBAL algorithm provides a total balance of surface radiation and energy together with sensible heat flow and aerodynamic roughness, and ET is calculated as a component of energy per pixel. As mentioned, SEBAL includes an algorithm that solves complete energy balance (Eq. 1).
where λET is the latent heat flux (W / m 2 ), R n is the net radiations, G is the soil heat flux (W / m 2 ), and H is the sensible heat flux (W / m 2 ) (Bastiaanssen and Ali 2002).

Net solar radiation (Rn)
The Rn in each pixel was calculated using Eq. (2), where α is the surface albedo (dimensionless), RS↓ is the incoming short wavelength radiation flux, RL↓ is incoming wavelength radiation flux, RL↑ is the outgoing long wavelength radiation flux, and εo is the surface emissivity (a dimensionless quantity). (1) In the following, the relationships required are introduced for the calculation of the five factors mentioned above.

Soil heat flux (SHF)
Directly, it is not possible to measure soil heat flux (SHF) using RS, but associations between the G R n value and such factors as NDVI, surface temperature, and albedo have been reported in many investigations. In the present investigation, SHF was estimated with empirical Eq. (3) developed by Bastiaanssen 2000, The Normalized Difference Vegetation Index (NDVI) is showed rate of vegetation coverage and its condition. This index estimates the earth's surface reflectance, leaf area index, area under cultivation, and plant biomass growth intensity (Bastianssen and Chandrapala 2003). The NDVI index is calculated by Eq. (4), where NIR and R are, respectively, reflectances in the near-and the red infrared bands. This index ranges between − 1 and + 1. Densely vegetated lands have positive values of 0.3 to 0.8, whereas negative values belong to snowcovered regions.

Sensible heat flux
The sensible heat flux of air value is obtained from wind speed and the earth surface temperature using a unique internal calibration consisting of the difference between the earth surface temperature and the adjacent air temperature (dT). The formula developed by Bastiaanssen 2000 is calculated through Eq. (5), where ρ air is air density (kg/m 3 ), C p indicates specific heat of the air (J/kg/K), dT represents temperature difference T 1 and T 2 (K) between two heights Z 1 and Z 2 , and r ah is aerodynamic resistance to heat transport( sm −1 ).

Surface Energy Balance System (SEBS)
The Surface Energy Balance System (SEBS) was proposed by Su and Jacobs 2001, to estimate heat flow fluxes and evaporative fractions. Similar to the SEBAL algorithm, the SEBS algorithm is based on the energy balance equation (Eq. 1). Among the components of the energy balance equation, the calculation method for net solar radiation (R n ) in SEBAL is identical to that of SEBS. To avoid repetition, the equations mentioned above are not presented in this section. Here, the equations are present for calculating two other components of the energy balance equation in the SEBS algorithm.

Soil heat flux (SHF)
Soil heat flux (SHF) in the SEBS algorithm is determined from Eq. (6), where Γ c is the ratio of SHF to net radiation (R n ) for dense vegetation, which is considered equal to 0.05. Γ s is the ratio of SHF to R n for bare soil, which is taken equal to 0.315, and f c is the partial canopy coverage.

Sensible heat flux
In the SEBS algorithm, the sensible and latent heat flux are obtained from a similar theory. The definitions of planet boundary layer (PBL), atmospheric boundary layer (ABL), and atmospheric surface layer (ASL) are differentiated in this algorithm. ABL refers to the part of the atmosphere directly affected by the earth's surface reactions and forces at a timescale below one hour. ASL represents the 10% lower ABL in which tensions and turbulent fluxes change less than 10%. In the ASL layer, the simulation relationships for the average wind speed profile (u) are written as Eq. 7.
In these equations, u and u * represent wind speed and friction velocity (m/s), respectively, k is the von Karman constant equal to 0.4, z indicates the reference height (m), d 0 is the zero displacement height (m), z om denotes roughness height for momentum (m), Ψ m is the stability correction factor for atmospheric heat transfer, and L indicates the Monin-Obukhov length (m) (Su and Jacobs 2001).
To avoid the repeated introduction of more details on the applied equations, readers are referred to articles Bastiaanssen 2000; Bastiaanssen and Ali 2002 (for the SEBAL algorithm) and article Su and Jacobs 2001, (for the SEBS algorithm).
lysimeter Cylindrical lysimeters (with 0.5 m diameter and 1.2 m height) were established to directly measure ET and crop coefficients. After grass and wheat crop ET (ET c ) measurements, crop coefficients (K c ) were obtained at each growth stage using Eq. (8), where ET c and ET 0 represent the evapotranspiration of the crop and reference crop (mm/day), respectively, and K c is the crop coefficient (dimensionless).

Empirical methods of ET 0 estimation
In this research, the results of SEBS and SEBAL RS algorithms were evaluated using lysimeter and conventional empirical methods for ET 0 estimation to introduce the best method using statistical indices. The empirical methods used by each group are introduced in Tables 6, 7, 8. In Tables 6, 7, and 8, ET 0 is reference crop evapotranspiration (mm/day),T a ,T max ,T min , and T d are average, maximum, minimum, and dew point temperature, respectively (˚C),R a , extraterrestrial solar radiation (MJ/m 2 /day), a and b, empirical coefficients were presented by Doorenbos-Pruitt (1977), P, coefficient related to the length of the day, u, wind speed at the height of 2 m (m/s), n, actual sunshine duration (hr), φ, latitude (rad).RH , average relative humidity (%),R s , solar radiation (MJ/m 2 /day), h, elevation of sea level (m), K andα , empirical coefficients, Rn, is net solar radiation (MJ/ m 2 /day), e Smax ande Smin , maximum and minimum saturation vapor pressure, respectively (kPa), ∆ represents the slope of the vapor pressure curve (kPa/˚C),λ , latent heat of vaporization (MJ/kg), P, vapor pressure (kPa), Rnl and Rns, net solar radiation with long and short wavelength, respectively (MJ/m 2 /day), e x , ande a , saturated and actual vapor pressure, respectively (kPa).
If wind speed is measured at an height other than 2 m, it should be converted to speed at a 2 m height for use in the FAO-PM formula, the general equation of which is shown in the following: In this equation, U 2m and U Z are, respectively, wind speeds at height of 2 and Z m, and Z is the height measured at which wind speed.
Since this study used data from weather stations, where wind speed is measured at 10 m height, the wind speed was converted to a 2 m height using Eq. (10).

Statistical indicators
To validate the results, the ET values estimated by SEBAL and SEBS algorithms were compared with lysimeter, FAO-PM, and empirical methods through conventional statistical indicators. In this investigation, the best method for the study area was determined using R 2 , RMSE (mm/day), PBIAS (Lang et al. 2017), MBE, MAPE, and NS statistical indices (Ghaderi et al. 2020).   (1966) T, Tmin, Tmax, RH, u P i and O i are the values predicted with each RS algorithm and those obtained from comparative methods, respectively, P and O are the mean values predicted with each RS algorithm and those obtained from comparative methods, respectively, and n represents the total data.

Crop coefficient
The wheat crop coefficient was calculated at each cultivation stage by considering reference crop ET as ET 0 and wheat crop ET of lysimeter data (ET C ). According to Table 9, the highest crop coefficient (kc) belongs to the development stage on 29/12/2019 (42 days after cultivation), and the lowest level (0.58) was recorded for the final wheat growth period on 21/05/2020 (22 days after cultivation). The total actual ET was measured 460.1 mm during wheat cultivation for 144 days. The maximum and minimum water requirements were reported to be 231.23 and 19.47 mm/day during plant growth, respectively, in the Einkhosh Plain of Ilam in Iran (Ghaderi et al. 2020).
Also, the ET value of 203 mm was obtained for rainfed wheat for 155 days from November to April, in Khuzestan province (Porgholam Amiji et al., 2019). This difference can be attributed to cultivation conditions and time.

Calculation of wheat evapotranspiration using satellite image SV
To better illustrate different ET values during the growth period, the maximum and minimum values of wheat ET values are presented in Table 10. Accordingly, the minimum wheat ET values based on SEBAL and SEBS algorithms (0.28 and 0.93 mm/day) were recorded on December 29 and November 27, respectively. Moreover, the maximum wheat ET values based on SEBAL and SEBS algorithms (5 and 7.95 mm/day, respectively) occurred on April 03, 2020. Figures 2 and 3 depict the spatiotemporal changes of daily ET values in the whole study area for SEBAL and SEBS algorithms, respectively. As shown in the figures, the ET dispersion obtained from SEBAL images is slightly more from SEBS images. Reasons of this observation may be a high estimation of the net radiation flux, which is confirmed by Wang et al. 2017 andKhand et al. 2021.

Comparison of evapotranspiration estimated with lysimeter method
Tables 11 and 12 compare the ET values estimated by SEBAL and SEBS algorithms with those of the lysimeter method. According to Table 11, the maximum ET values estimated by SEBAL and SEBS algorithms (16.67 and 14.21 mm/day, respectively), similar to the lysimeter (11.50 mm/day), were recorded on 2020/05/21. The minimum ET values estimated by these two algorithms (1.29 and 1.22 mm/day, respectively), similar to the lysimeter (1.70 mm/day), were documented on 2019/12/29. Figure 4 displays the dispersion of ET values estimated by SEBAL and SEBS algorithms compared to the lysimeter method. As indicated by the assessments, SEBAL and SEBS algorithms correspond to the actual lysimeter method, and the results of energy balance algorithms can be generalized to the lysimeter method. These findings agree with those reported by Rawat et al. 2017, Asadi and Valizadeh Kamran. 2022, Tariqul Islam et al. 2023, Liu et al. 2023, and Khoshnood et al. 2023. Table 12 compares R 2 and RMSE values obtained for SEBAL (0.92 and 2.15 mm/day) and SEBS (0.96 and 1.53 mm/day) algorithms with the lysimeter method. Thus, the estimated and observed ET values are close to each other on the mentioned days. Hence, SEBS and SEBAL algorithms are relatively good alternatives to the lysimeter Page 9 of 15 137 Fig. 2 The process of spatial and temporal changes of actual wheat ET using the SEBAL algorithm

Comparison of evapotranspiration estimated and empirical methods
According to Table 13, the maximum ET values estimated from SEBS and SEBAL algorithms (16.67 and 14.21 mm/ day, respectively) and empirical methods (except for the Blaney-Criddle method) were recorded on 2020/05/21. The minimum ET values estimated from SEBS and SEBAL algorithms (1.67 and 1.22 mm/day, respectively) were similar to FAO-PM (1.80 mm/day) on 2019/12/29. Figures 5 and 6 illustrate the dispersion of ET values estimated by SEBAL and SEBS algorithms compared to empirical methods. According to the evaluations, SEBAL and SEBS algorithms (R 2 = 0.79 and 0.87, respectively) showed high compatibility with the FAO-PM method. Among temperature-based methods, the Hargreaves-Samani method showed good accuracy compared to SEBAL (R 2 = 0.72) and SEBS (R 2 = 0.80) data.  Also, in the temperature-based group, a weak estimation was observed for the Blaney-Criddle method relative to SEBAL and SEBS algorithms (R 2 = 0.31 and 0.39, respectively). The reasons may include limited parameters in this method and its unsuitability according to the climate of the study area. Among radiation-based methods, the Doorenbos-Pruitt method showed good outcomes with SEBAL (R 2 = 0.71) and SEBS (R 2 = 0.79) data. Thus, the Doorenbos-Pruitt method performs better in arid climates among radiation-based methods and can be used for ET estimation in arid and semiarid areas. Additionally, the Doorenbos-Pruitt method needs more input parameters than the other methods, which can account for its improved efficiency compared to SEBS and SEBAL. In the mass transfer-based group, both Mahringer and WMO methods presented less effective estimations versus SEBAL (R 2 = 0.60 and 0.68, respectively) and SEBS (R 2 = 0.55 and 0.64, respectively). This might be due to the unsuitability of related parameters in mass transfer-based methods with energy balance algorithms in the study area. It is noteworthy that a limitation of SEBAL and SEBS algorithms is that the presence of some empirical relationships during ET estimation may result in errors. These models also require a bright cloudless sky because even a thin cloud layer can reduce the estimated heat radiation and energy, leading to many errors. Therefore, better results were obtained with the Mahringer method among mass transfer-based methods. According to the results, SEBAL and SEBS algorithms mainly were compatible with the actual lysimeter method (R 2 = 0.92 and 0.96, respectively). The high accuracy of these two algorithms in ET estimation suggests their copious applicability for studying large extents. Table 14 compares the estimated data obtained from SEBAL and SEBS with ET data measured by the FAO-PM reference, Hargreaves-Samani, Blaney-Criddle, Priestley-Taylor, Doorenbos-Pruitt, Mahringer, and WMO methods. According to Table 14, the RMSE values for the FAO-PM method were obtained at 2.42 and 3.14 mm/day for SEBS and SEBAL algorithms, respectively. The scatter plots (Figs. 5 and 6) represent less dispersion and more correlation of SEBS than SEBAL compared to the other methods. The scatter plots show higher R 2 and lower RMSE values for SEBS (with a slight difference) than SEBAL, with more acceptable accuracy closer to actual values. The PBIAS index also reveals underestimations for both SEBAL and SEBS, with a more significant underestimation for SEBAL. This result can be attributed to the high estimation of the net radiation flux parameter, the main factor among the fluxes in the energy balance equation, which has caused a more Fig. 5 Scatter plot of actual ET estimated with SEBAL algorithm and empirical methods significant underestimation of SEBAL than SEBS. It is noteworthy that the heat flux is low in farmlands or lands with intensive and semi-intensive vegetation. The soil heat flux is also low in densely vegetated lands. This is because the soil under vegetation is not exposed to radiation energy in dense vegetation, and the flux reaches a medium value with decreasing vegetation intensity. SHFs decline during the growing season because of plant growth and increasing crop coverage. The heat flux will be high in barren and arid zones. Thus, the SEBAL algorithm shows some underestimation in comparison with direct and indirect methods, which can be attributed to a low SHF in the SEBAL model. Compared SEBAL (MBE = − 1.17,MAPE = 33%, and NS = 0.70) and SEBS (MBE = − 0.92,MAPE = 20%, and NS = 0.62) with empirical methods, the FAO-PM method showed the most compatibility with the data of these two algorithms among the other empirical methods. Based on the comparison of all tested methods with SEBAL (Table 15), the first, second, and third ranks belong to FAO-PM, Doorenbos-Pruitt, and Hargreaves-Samani methods, respectively, followed by   (Wang etal. 2017, Rawat et al. 2017, Khand et al. 2021, Wei et al. 2022, Liu et al. 2023, and Khoshnood et al. 2023 suggest the ability of energy balance algorithms, in particular SEBS, in crop actual ET estimation. Among temperature-based methods, the Hargreaves-Samani showed better performance, which agrees with Lang et al. 2017, andZoratipour et al. 2019. Among radiation-based methods, the Doorenbos-Pruitt performed better than the Priestley-Taylor, which is similar to its good performance in ET estimation reported by Wang et al. 2017, Zoratipour et al. 2019, and Shamloo et al. 2021. Among mass transfer-based methods, the Mahringer method presented a better estimation, as reported by Djaman et al. 2015 andZoratipour et al. 2019.

Conclusion
In this research, wheat daily ET values were estimated accurately using SEBAL and SEBS algorithms as well as Landsat 8 satellite images at six satellite passing dates from 2019 to 2020, and their results were compared with the lysimeter method. The results demonstrated that the ETc estimated with SEBS and SEBAL (R 2 = 0.96 and 0.92, RMSE = 1.53 and 2.42 mm/day, respectively) corresponded well to actual lysimeter data and can be compared with empirical methods. Therefore, SEBS and SEBAL results were compared with FAO-PM, Hargreaves-Samani, and Blaney-Criddle temperature-based methods, Priestley-Taylor and Doorenbos-Pruitt radiation-based methods, and Mahringer and WMO mass transfer-based methods, using statistical indices. Based on the findings, FAO-PM showed good compatibility with SEBS and SEBAL (R 2 = 0.87 and 0.79, RMSE = 2.42 and 3.14 mm/day, respectively) among the mentioned empirical methods. Altogether, the ET values obtained from SEBAL, SEBS, and the lysimeter method (26.9, 33.9, and 28.38 mm/day, respectively) were not highly different. In addition, the ET value estimated from SEBS is about 12% higher than that obtained from SEBAL. Thus, SEBS with more outstanding scores (R 2 = 0.96 and RMSE = 1.53 mm/ day) presented better results than SEBAL. Compared to the lysimeter method, SEBS and SEBAL algorithms are appropriate alternatives for use in the study area owing to their high compatibility in the absence of adequate data and no use of direct methods.