Calibration and validation of solar radiation-based equations to estimate crop evapotranspiration in a semi-arid climate

The Abtew and Jensen-Haise solar radiation-based equations were used to estimate evapotranspiration, considering the limited climatic data in many locations. Both equations were proven to successfully predict the potential evapotranspiration (ETO) compared with the standard Penman–Monteith (PM) method in two Mediterranean countries. Calibration of the constant coefficient k of the Abtew equation showed substantial differences compared to recommended values (1.22 vs. 0.53), with the highest values observed during September (1.46). Validation of ETO measurements using calibrated Abtew equation against the PM method indicated a high correlation coefficient (r2 = 0.97, RMSE = 0.61). Further, evapotranspiration requirements, using the calibrated empirical equation, were calculated for olives (449 mm) and citrus (807 mm) showing a good agreement with recommended values for dry climate regions. Therefore, the tested equations could be safely used to predict frequencies and doses of irrigation in semi-arid climates, considering limited climatic data availability.


Introduction
The water consumption of a crop (i.e., Evapotranspiration-ET C -; soil evaporation-E, and plant transpiration-T) depends mainly on a daily phenomenon during which a transition occurs from the liquid phase of water to the vapour phase.It causes a reduction in plants' temperature and requires a certain amount of energy, which is provided by the available radiative and convective energy at the plant and soil evaporating surfaces.This available energy represents the crop water demand, ET C (Kittas 1990).In arid and semi-arid climates, agriculture is closely related to the rate of evapotranspiration and the probability of precipitation.For example, under arid conditions, 95% of the annual precipitation is estimated to be consumed by ET C (Melesse et al. 2009).Thus, considering the water and energy crises, the need for information on evapotranspiration has become more important today than a few decades ago.Frequently in agriculture practice, we estimate ET C by multiplying potential evapotranspiration (i.e., ET O ; the capacity of the atmosphere to remove water from a surface over a specified region) from a reference canopy such as the turf grass under a non-limiting water supply with a "crop coefficient value", which characterizes each crop and is related to its developmental stage (Nikolaou et al. 2022).However, experimental determination of ET O is only possible at very limited sites because of its complexity and the high cost of equipment required (Flores-Velazquez et al. 2022).
Today, ET O could be estimated based on real-time climatic data (i.e., wind speed, relative humidity, air temperature, solar radiation, and sunshine hours) recorded by automated agro-meteorological climatic stations (Shirmohammadi-Aliakbarkhani and Saberali 2020).The Penman-Monteith (PM) method, released by the Food and Agriculture Organization (FAO), represents an internationally recommended model since 1990 that could be used for evapotranspiration estimation (Allen et al. 2006;Achparaki et al. 2012;Wang et al. 2021).However, in many locations, the complete dataset of meteorological variables required for the PM method is not available (Bogawski and Bednorz 2014).Thus to overcome this problem, several scientists and researchers from around the world applied a variety of empirical equations to calculate evapotranspiration (such as the Hargraves, the Blanley Criddle, the Thornthwaite, and the Makkink), but none of them could be deemed flawless due to the vast variations in climatic conditions in different parts of the world (Yates and Strzepe 1994;Faruk Bin Poyen et al. 2016).However, the strong dependence of evaporation on the radiation energy term has been generally accepted (Islam and Rashidul Alam 2021;Flores-Velazquez et al. 2022).Evapotranspiration estimation based on solar radiation measurements eliminates the effect of the surface albedo and minimizes the contribution of the aerodynamic term.Xu and Singh (2000) found that under clear sky conditions, there is a great balance between evapotranspiration accuracy estimation, simplicity of the solar radiation method, and robustness.However, under low solar radiation and cloudiness, evapotranspiration underestimation was a common problem since only an amount between 70 and 75% of the daily ET C variance is explained by solar radiation alone (Melesse et al. 2009).In another case, for soillessbased greenhouse crops the accumulated solar-radiation method has been applied for a long time, matching the diurnal evapotranspiration fluctuation with solar radiation as a sustainable water-saving approach (Katsoulas et al. 2006;Nikolaou et al. 2017).However, the main disadvantage of using equations with very few input data is the performance efficiency under different climate systems as mentioned elsewhere (Aschale et al. 2022).
Among radiation-based methods, for warm climates, the Jensen-Haise method has a very good performance rating for ET O calculation (Shirmohammadi-Aliakbarkhani and Saberali 2020; Gharehbaghi and Kaya 2022).The Jensen-Haise was initially calibrated under semi-arid conditions using solar radiation and air temperature (Jensen and Haise 1963).The primary gain in Jensen-Haise accuracy comes from the inclusion of air temperature into the ET O equation.The argument to use temperature is that both components of evaporation in the PM equation are related to the air temperature (Yates and Strzepe 1994).Abtew is another simple method that could be applied for ET O calculations where the only available climatic data is solar radiation (Mengistu and Amente 2017;Islam and Rashidul Alam 2021).It was originally developed for warm and humid environments; therefore, recalibration under arid and semi-arid conditions is recommended to increase the model's accuracy ( Xu and Singh 2000;Samaras et al. 2014;Mengistu and Amente 2017).Considering the effectiveness of the Abtew equation, the k coefficient represents the relation between solar radiation and water consumption in the reference canopy used for the determination of the potential evapotranspiration.Initially, a constant k value of 0.53 was proposed.Working within greenhouses (Kittas 1990) showed that k values ranged from 0.44 to 0.52, increasing up to 0.67 in a semiarid environment.
Irrigation scheduling is not widely used in many parts of the world, despite the abundance of evapotranspiration estimation-based decision support systems that have been created over the past few decades (Giannakis et al. 2016;Taghvaeian et al. 2020).Therefore, the adaptability of simpler methods for use in predicting when and how much water is required for any particular irrigation scheme is worth investigating, especially in geographical areas where there is limited climate monitoring (Melesse et al. 2009).Therefore, in the present manuscript, we focus on two empirically based ET O equations that use limited and easily recorded climatic data (i.e., solar radiation and air temperature).Given that irrigation water is available, a wide variety of crops can be grown in the Mediterranean because of its temperate climate.For example, in arid and semi-arid regions, irrigated open-field crops are usually cultivated from March to September, in most cases under completely clear sky conditions.The total irrigation requirements for the main crops may rise from 350 (e.g., vegetables) to 2400 (e.g., colocasia) mm (i.e., 3500-24000 cubic meters of water per hectare) (Dalias et al. 2019;Nikolaou et al. 2020a, b).Because crops in these areas depend on sufficient supplies of highquality water, the availability of water has always been an issue (Sánchez-Molina et al. 2015;Nikolaou et al. 2020a, b).Particularly, olive crops (Olea europaea L.) are considered one of the most economically and ecologically important tree crops in the Mediterranean area (Sofoulaki et al. 2023).Olive is resistant to aridity; however, the higher frequency and severity of droughts in the future would result in an average increase of 18.5% of the irrigation demand over the Mediterranean (Fraga et al. 2021).On the other hand, there has been a tendency over the last decade toward high plant olive crop density (up to 2500 plants ha −1 as opposed to the current practice of 350 plants ha −1 ) as it has been reported to be a good strategy in terms of orchard productivity under semi-arid conditions (Egea et al. 2017).
Citrus is another crop widespread worldwide, with the countries around the Mediterranean basin constituting important producers (Sofoulaki et al. 2023).However, in semi-arid conditions, citrus evapotranspiration falls in the conventional range of 700-1300 (mm), with an average of 1000 (mm) during the irrigation period (Abou Ali et al. 2023).Thus, it is important to determine the citrus water requirements by designing simple, practical, and precise methods to optimize irrigation and to adopt precise irrigation scheduling and management techniques aiming for water savings (Puig-Sirera et al. 2021).
In view of the above, in this study, climatic data from two semi-arid Mediterranean regions, Cyprus and Italy, were utilized to compare the performances of the Abtew and Jensen-Haise radiation-based equations with the Penman-Monteith equation for calculating potential evapotranspiration (ET O ).A modified Abtew equation based on extraterrestrial solar radiation is also tested for citrus and olive crop ET C calculation as an ex-ante irrigation scheduling tool.

Experimental setup
For this study, three agro-meteorological stations were chosen, two of which were located in Cyprus and one in Italy, all of which were in semi-arid Mediterranean environments.
The selected study locations are characterized by the typical climates of the Mediterranean area, with mild winters (annual rainfall of about 320-340 mm between October and March), warm summers, and an average annual sunshine duration of about 3 332 h.The summer months are mostly dry, with relatively high global and horizontal beam radiation intensities.The Mediterranean climate means that evaporation rates are quite high but also reasonably constant from year to year.Historic data allows accurate predictions of evaporation losses, but such losses are considered inevitable (Cox 1999).
In each study location, an automatic agro-meteorological climatic station (MeteoSense 4.0; Netsens; Calenzano; Florence; Italy), with sensors mounted on a pole at 2 m height (Fig. 1).It consisted of a wind sensor (1-67 m s −1 , accuracy 5%; direction 0-360 • , accuracy 7 • ), a rain collector (tipping bucket, resolution 0.2 mm), a solar radiation sensor (0-1800 W m −2 , accuracy 5%), a thermo-hydrometer for air temperature (-25 to + 85 °C; accuracy 0.5 °C), and air humidity (0-100% RH; accuracy 3%).A solar panel is used as a source of power, integrated with a storage battery.Climatic data were forwarded in real-time to the Netsens LiveData platform (data cloud platform; IoT system) with the use of a SIM card.Climatic data were recorded at 30-s intervals and transmitted as average data every 5 min, while on the cloud platform, data can be downloaded on an hourly basis.

Penman-Monteith method
As a result of an expert consultation held in May 1990, the FAO Penman-Monteith method is recommended as the sole standard equation for the definition and computation of the reference evapotranspiration (Allen et al. 1998).
Where ET OPM is the estimated reference evapotranspiration value (mm d −1 ), R N is the net radiation at the crop surface (MJ m −2 day −1 ), G is the soil heat flux density (MJ m −2 day −1 ), T is the mean daily air temperature at 2 m height (°C), u 2 is the wind speed at 2 m height (m s −1 ), e s is the saturation vapour pressure (kPa), e a is the actual vapour pressure (kPa), e s -e a is the saturation vapor pressure deficit (kPa C), Δ is the slope of the vapour pressure curve (kPa C−1 ), γ is the psychrometric constant (kPa C−1 ).Jensen and Haise (1963) proposed an empirical equation for semi-arid and arid regions that estimated potential evapotranspiration based on solar radiation and air temperature as follows:

Jensen and Haise method
Where ET ΟJH is the estimated reference evapotranspiration value (mm d −1 ), T α the mean air temperature (°C d −1 ), R S is the mean solar radiation (MJ m −2 d −1 ).The air temperature indirectly introduced the contribution of the aerodynamic term into the estimate of ET O .The equation, which is referred to as a radiation-based model, has been derived from data collected in arid regions of the western part of the United States and should yield satisfactory results in areas with similar climates (Jensen and Haise 1963).

Abtew method
To account for cases where the only available climatic data is solar radiation, the Abtew equation is recommended (Abtew 1996): Where ET ΟΑ is the estimated reference evapotranspiration value (mm d −1 ), k represents a dimensionless coefficient, R S is the mean solar radiation (MJ m −2 d −1 ); λ is the latent heat vaporization (2.45 MJ kg −1 according to Allen et al. 2006).The Abtew equation was originally developed and used in warm and humid wetland environments (Samaras et al. 2014).Therefore, we proceeded with the calibration of the Abtew k coefficient for semi-arid environments by rewriting Eq. 3 as: For the calculation of potential evapotranspiration in mm d −1 , the solar radiation (R S ) must be converted to mm d −1 (2) ET OPM R S (Allen et al. 1998;Jaafar and Ahmad 2019;Mengistu and Amente 2017): Where ρ w is the density of water (1.000 kg m −3 ).

Estimating crop evapotranspiration (ETc)
In the single crop coefficient approach, the crop coefficient (K C ) integrates differences in crop transpiration rates and soil water evaporation between the crop and the grass reference surface.The Kc is the ratio of ETc to ETo; it represents an integration of the effects that distinguish a crop from the reference grass (Xiang et al. 2020).In the crop coefficient method, evapotranspiration (ET C ) is given below (Simbeye et al. 2023): Where ET C is the estimated crop evapotranspiration value (mm d −1 ), K C is a crop coefficient (dimensionless; varies with the growth stage of the crop), and ET O is a reference crop evapotranspiration (mm d −1 ).

Estimating ETc using a modified Abtew equation
The earth's solar radiation fluxes (R S ) can be estimated under clear sky conditions based on extraterrestrial (R A ) radiation (Appendix Table 4) and station elevation above sea level, following Allen, (1998) transformation equation (R S = 0.75 + 2 10-5 ζ) R A /λ).The importance of estimating R Ν based on R A could be applied in cases of total absence of current measurements of radiation.The modified Abtew equation can be rewritten as follows: Subsequently, crop evapotranspiration (ET C ) can be estimated by substituting ( 7) in (6) using the following formula: Where R A is the extraterrestrial radiation (Mj m −2 d −1 ), ζ is the station elevation (m), K C is the crop coefficient.

Experimental methodology
Considering climatic data recorded in stations A and B (Cyprus) and calculated potential evapotranspiration  2008) for olives (i.e., K C ini = 0.65, K C mid = 0.45, K C late = 0.65) and by Jamshidi et al. (2020) for citrus (i.e., K C = 0.71 to 0.96) for semi-arid regions.The olive season lasts in the Mediterranean regions from early March to November, and the lengths of the several crop development stages (L-ini, L-dev, L-mid, and L-late) are, respectively, 30, 90, 60, and 90 days.For citrus crops, the growing season typically starts in early February with flower bud induction, followed by flowering from mid-March to April.The length of crop development stages (L-ini, L-dev, L-mid, and L-late) is, respectively, 60, 90, 120, and 95 days (Allen et al. 1998;Fraga et al. 2021;Abou Ali et al. 2023).
The reliability of empirical equations for calculating ET O increases when the methods are calibrated for each crop in each region.This is necessary because the effective rooting depth and the permissible water deficit for each crop, as well as soil-water retention characteristics, are factors that have to be taken into consideration.In addition, the considerable advection of energy from unirrigated surroundings (oasis effects) affects the capacity of the atmosphere to remove water from a surface over a specified region of each location.

Statistical analysis
Climatic data were analyzed and comparisons of means were tested with ANOVA using a Statistical Package for the Social Sciences (IBM Corp. Release 2011.IBM SPSS Statistics for Windows, Version 20.0.Armonk, NY, USA: IBM Corp).Regression analysis was performed for the estimation of relationships between selected data.

Climatic conditions
Table 1 summarizes the mean values of air temperature, relative humidity, solar radiation, evapotranspiration (ET OPM ; potential evapotranspiration; Penman-Monteith FAO-56; Eq. 1), and 10-day accumulated precipitation for each agro-meteorological station during daylight hours.The mean daily solar radiation values were slightly higher over the study period for station A at 448 W m −2 (maximum 869 W m −2 ) compared to station B at 442 W m −2 (maximum 866 W m −2 ) or C at 411 W m −2 (maximum 899 W m −2 ) with an exception during the second 10-day interval measurement period in July where station C recorded higher solar radiation (Table 1 and Fig. 2A).This was probably due to a Sahara dust transport to Cyprus, reducing the amount of solar radiation reaching the soil surface.The minimum daily sunshine hours were 5 in March and increased to 12 in July.
Higher air temperature values of 3.4 C (station A) and 3.7 C (station B) were recorded in Cyprus compared with station C located in Italy.The maximum air temperature observed in July was 39.8 C (station A), 37.9 C (station B), 38.7 C (station C).The mean air temperature for a three-day representative period in July is illustrated in Fig. 2B.In July, during the second 10-day interval period, station C recorded the highest air temperature values.
In station A, the total recorded precipitation through the irrigation period was 53 mm.However, 43 mm of rain was recorded only on one day in June.In station B, the total precipitation was 3.7 mm, and in station C, it was 21.8 mm.Considering the number of rainy days in relation to the amount of precipitation, the effective rainfall (i.e., green water) during the irrigated period in the semi-arid Mediterranean region was considered negligible.
Over the study period, the mean ET OPM values in the study regions ranged from 1.1 to 6.1 mm (maximum value of 6.7; station A).Significantly higher potential evapotranspiration (ET OPM ) rates (p ≤ 0.05) were estimated for station A (12 and 18% increases compared to station B and station C respectively).The mean ET OPM was 4.4 mm (station A), 3.9 mm (station B), and 3.7 mm (station C). Figure 2C shows the diurnal ET OPM for a three-day period in July.It can be observed that in stations A and B, the ET OPM values followed a similar daily trend.However, higher values of ET OPM were recorded for station C during and after midday hours (Fig. 2D).
Figure 3 shows the mean wind velocities (Fig. 3A) and atmospheric pressure (Fig. 3B) during the measurement period in the three locations.It can be observed that the wind velocity was consistently lower at station B compared with stations A or C for the study period.These values show the expected differences caused by the crop surrounding the station.Considering the atmospheric pressure between Italy (station C) and Cyprus (stations A and B), higher values were observed on behalf of station C (Fig. 3B).

Comparison of the Abtew method with the Penman-Monteith
Considering the net solar radiation values recorded in stations A and B and reference evapotranspiration based on the Penman-Monteith equation (ET OPM ), we proceeded with the Abtew k coefficient calibration (Eq.4).The calculated mean k coefficient values (k = λ ΕΤ ΟPM / R S ) were in the range of 1 to 1.46 (dimensionless), with an overall mean value of 1.22 (± 0.01) for stations A and B (Table 2).Over the study period (March to September), on 10-day interval measurement calculations, the result shows variations in k coefficient values within the same station and between different stations.The higher mean k values were calculated for both stations during September (avg.1.46).
The mean daily variation of potential evapotranspiration calculated based on the Penman-Monteith method (eq. 1) and the Abtew method (eq.3, considering a new k coefficient value of 1.22) over a 10-day interval period (March to September), are shown in Fig. 4 (agro-meteorological station C, altitude 60 m, Italy).It can be observed that ET OPM and ET OA followed a similar trend over the study period.The slight over-and under-estimation of ET OA could be explained due to variations of the "constant k coefficient".Thus, using the calibrated k coefficient value (Cyprus) to measure ET OA against ET OPM , in another location (Italy) with similar conditions, we developed a linear fitting of the measured parameters (Fig. 5).The r 2 values (0.97) for many observations (n = 200) were statistically significant (P < 0.05), with a beta value of 0.98.The root mean square error (RMSE) was 0.61 mm d −1 .The mean daily ET OA was estimated at 3.75 (± 0.09) mm and ET OPM at 3.84 (± 0.12) mm, which statistically (p ≤ 0.05) are considered identical.

Comparison of the Jensen and Haise method with the Penman-Monteith
The mean daily variation of potential evapotranspiration estimated with the Penman-Monteith (ET OPM ; Eq. 1) and with the Jense and Haise method (ET OJH ; Eq. 2) based on climatic parameters derived from station C is shown in Fig. 6.It can be observed that the two methods applied for potential evapotranspiration estimation showed a similar trend, with a slight underestimation during the first two months.A significant linear regression between ET OJH and ET OPM for different locations is shown in (Fig. 7), suggesting a good correspondence.The mean daily ET OJH values were: 3.99 ± 0.11 mm (station A), 3.84 ± 0.10 mm (station B), and 3.48 ± 0.10 mm (station C), and accordingly, ET OPM values were 4.41 (± 0.10 mm-station A), 3.97 (± 0.09 mm-station B), 3.48 (± 0.10 mm-station C).The r 2 values for several observations (n = 200) were statistically significant (P < 0.05) accounting for 0.97 with a beta value of 0.98 for station A, 0.98 with a beta value of 0.97 for station B, and 0.97 with a beta value of 0.93 for station C.

Crop evapotranspiration requirements based on the modified Abtew equation
Using the single crop coefficient method (Eq.6) and the modified Abtew Eq. ( 7), the evapotranspiration requirements for olives and citrus were estimated for Cyprus (35.1264°N, 33.4299° E) on a ten-day interval basis (Table 3), considering historical extraterrestrial radiation (R A ; Appendix Table 4) (Allen et al. 1998).The evapotranspiration requirements were also estimated for the same crops following the single crop coefficient method and the reference evapotranspiration (Penman-Monteith; Equation; Table 1).The results are illustrated in Table 3.
The evapotranspiration requirements for olives and citrus, estimated using the modified Abtew equation with calibrated k coefficient and extraterrestrial radiation, closely matched the estimates obtained using the Penman-Monteith method (Table , 3) Furthermore, tresults on evapotranspiration for olive and citrus crops using a class A evaporation pan, along with the application of local crop coefficient values (K C values), are in line with the findings of the current study (Christou et al. 2017;Nikolaou et al. 2020a, b).

Discussion
The reliability of evapotranspiration models relying on limited climatic inputs is an important challenge in dry Mediterranean regions because the FAO-56 Penman-Monteith method  Fig. 6 Mean daily variation of potential evapotranspiration estimated with the Penman-Monteith (ET OPM ; green line) and with the Jense and Haise method (ET OJH ; blue line) based on climatic parameters derived from station C (agro-meteorological station, altitude 60 m, Italy) over a 10-day interval period starting in March till September cannot be applied in many situations due to the poor weather data collection facilities (Aschale et al. 2022).Thus, empirical equations of evapotranspiration based on easily measured parameters such as solar radiation that are adapted to the prevailing environmental conditions are considered useful tools, particularly in southern Europe, where irrigation scheduling is currently based on the experience of the growers (Incrocci et al. 2020).For example, incorporating these equations into irrigation schemes will enhance water productivity at the farm level.Furthermore, in the case of a large-scale irrigation policy, energy data could be retrieved based on satellite remote sensing or historical extraterrestrial radiation (Appendix Table 4).
Considering input climatic data and estimated potential evapotranspiration values (ET OPM ), it was observed that ET OPM in station B (agro-meteorological station, altitude 330 m, Cyprus) was lower than in station A (agro-meteorological station, altitude 165 m, Cyprus) and higher than in station C (agrometeorological station, altitude 60 m, Italy).In this context, many researchers have investigated the effect of microclimate on potential evapotranspiration.They concluded that solar radiation has the largest impact on the reference evapotranspiration values, followed by the air temperature, the saturation vapour pressure deficit, and lastly, the wind speed (Maček et al. 2018).Considering stations A and B (which are both located in  Cyprus), solar radiation and air temperature values were higher for station A, thereby resulting in higher evapotranspiration values.However, the higher wind speed that was observed at station C was not able to explain the lower evapotranspiration rate compared with station B. It is clear from the results that the effect of solar radiation was more intense on potential evapotranspiration than the air temperature values or wind speed, as may be expected.In addition, higher air humidity values in station C could also be associated with lower evapotranspiration rates.Considering the atmospheric pressure between Italy (station C) and Cyprus (stations A and B), higher values were observed on behalf of station C (Fig. 3).A distinct relationship between pressure and evaporative values is observed; thereby, at higher pressures, lower evaporation is expected.
It has been shown that the radiation-based methods tested have similar performance standards to the PM method in semiarid (type BSh) and hot Mediterranean (type Csa) climatic regions, and they can estimate reference evapotranspiration with high accuracy once the equation coefficients are adjusted to the local environmental conditions.Those k values of the calibrated Abtew equation obtained under our experimental conditions were much higher than those initially proposed and used by many researchers or even those found in the literature.However, according to Samaras et al. (2014), the k coefficient in Abrtew's model for evapotranspiration calculation is also related to the air Tmax of the region.Indeed, this is justified in our case by the lower values of the k coefficient for station B compared to station A, where the higher values of air Tmax were recorded, and probably explains the higher values of k found in our experiments.
Xu and Singh 2000, working with several radiation-based equations for determining evaporation, indicated that the original constant value of the Abtew equation (k = 0.53) agreed most closely with pan evaporation in Switzerland, compared with other models without recalibration (i.e., Hargreaves, Makkink, Priestley and Taylor and Turc).In contrast, Shirmohammadi-Aliakbarkhani and Saberali (2020) suggested that the Abtew equation based on a coefficient k = 0.53 might be unreliable for evapotranspiration estimation in arid regions of Iran.In West Africa, the calibrated Abtew method showed the best performance among nine tested ETo equations (Djaman et al. 2017).Similarly, (Samaras et al. 2014) reported that the Abtew model (calibrated and validated) showed the best overall performance to the data from all available climate stations under different Mediterranean climates in central Greece.In our case, with a regression slope of almost unity and a very high correlation coefficient between ET OPM and ET OA , the calibrated equation (k = 1.22) showed good performance and can be used for potential evapotranspiration estimation in semi-arid and hot Mediterranean climates, accounting for cases where solar radiation is the only available climatic parameter.In addition, the given modified Abtew equation using historical extraterrestrial solar radiation data (Appendix Table 4), can be a useful tool for decision-making when it comes to predicting water deliveries in advance.
The results also indicate that the Jensen and Haise model provides good estimates of ET O .Therefore, we may also consider this method for estimating potential evapotranspiration in semiarid and hot Mediterranean climate regions.Indeed, among the radiation-based methods, the Jensen − Haise method was the only one that exhibited consistent results in some regions (Shirmohammadi-Aliakbarkhani and Saberali 2020).Although other authors (Ahmadi and Javanbakht 2020) have suggested that the Jensen and Haise method tends to overestimate potential evapotranspiration in some cases.This can be explained, as cited by Yang et al. (2021), due to the differences in selected models, variables used to validate evaporation, evaluation criteria, or evaluation scales, different studies have reached different conclusions regarding the best radiation-based models.
Another point for consideration is the latent heat of vaporization (i.e., the amount of energy required to change liquid water into water vapor) in relation to the surrounding agro-meteorological station.The energy is provided by various sources, including the atmosphere above, the soil surrounding, and the inflowing water (Cox 1999).Comparing station A with station B, there was a great difference regarding the surrounding crop cultivation, thereby affecting the micro-environmental conditions.Station A was located in a field with aromatic plants (crop evapotranspiration requirements of 250-350 mm) compared with olive crops (crop evapotranspiration requirements of 450 mm).
Considering climatic data limitations and current results, a machine learning model based on the solar energy fluxes at the surface could be used as a simple and robust approach to controlling irrigation frequency.Given the daily or weekly amount of water needed by the crop and calculated with the abovedescribed methods, the frequency and duration of each irrigation shift can be automatically tuned, allowing the water to flow at root level while minimizing surface evaporation and stratum percolation.To achieve this aim, soil moisture sensors can be installed in the field to provide feedback on the water flow after the irrigation shift, to adapt the frequency and duration, and to ensure the proper amount of water is returned to the crop.

Conclusion
In semi-arid and hot Mediterranean climate regions, the Penman-Monteith method was compared to two empirical solar radiation methods (Abtew and Jensen and Haise) to compute daily reference evapotranspiration values.The findings demonstrated that the radiation-based methods were precise in estimating the potential evapotranspiration in the designated study areas.The Abtew method, in contrast to Jensen and Haise, is known for its simplicity when calculating reference evapotranspiration, provided that the empirical constant is calibrated.As a result, the Abtew method can be suggested as a novel water management approach for scheduling irrigation using smart irrigation controllers connected to a basic solar radiation sensor.Furthermore, by adjusting the Abtew equation to account for extraterrestrial radiation, effective irrigation requirements can be determined.Considering climatic data limitations in many locations, simple solar-radiation-based equations could be used for proper irrigation scheduling.

Appendix 1
Please see Table 4.

Fig. 2
Fig. 2 Daily solar radiation (A; mean hourly; W −2 ), daily air temperature (B; mean hourly; °C), and potential evapotranspiration based on the Penman-Monteith equation (C; mean hourly; mm) in a 3-representative day period in July; ET OPM in a 24-h period (D;

Fig. 3
Fig. 3 Mean daily wind velocities (A; Km h −1 ) and mean atmospheric pressure (B; hPa) over the study period; Red line (agrometeorologigal station A; altitude 165 m; Cyprus); Green line

Fig. 7
Fig. 7 Calculated potential evapotranspiration values (Jensen and Hasie equation; ET OJH ) plotted against potential evapotranspiration (Penman-Monteith equation; ET OPM ); Station A, agro-meteorological (ET OPM ; Penman-Monteith equation), we proceeded with calibration of the constant coefficient k of the Abtew equation.The model coefficient was then tested (ET OA ) and compared with the ET OPM in a different location and cropping conditions (station C; Italy).To check the suitability of the proposed modified Abtew equation (Eq.8),

Table 2
Ten-day interval

Table 3
Evapotranspiration requirements (mm) for olives and citrus estimated based on the modified Abtew equation and the FAO Penman-Monteith method R A, extraterrestrial radiation (MJ m −2 d −1 ); ET CA crop evapotranspiration (mm) based on the modified Abtew equation; ET CPM crop evapotranspiration (mm) based on the Penman-Monteith