Evapotranspiration and water availability response to climate change in the Middle East and North Africa

Quantifying the impact of climate change on evapotranspiration is necessary for devising accurate water and energy budgets in light of global warming. Nevertheless, in the Middle East and North Africa (MENA), little has been done to bridge this gap. This study, then, implements Penman and Budyko approaches to climatic data retrieved from the sixth phase of the Coupled Model Intercomparison Project (CMIP6) to assess evapotranspiration and water availability evolutions through the twenty-first century. Outcomes reveal that the MENA region is indeed vulnerable to a surge in temperature, which can increase evapotranspiration losses and decrease water availability. Under the shared socioeconomic pathway (SSP2-4.5), the potential evapotranspiration (PET) has been projected to increase throughout the MENA region by up to 0.37 mm per year during the middle of the twenty-first century (2021–2050) and by up to 0.51 mm per year during the end of the twenty-first century (2071–2100). Meanwhile, the actual evapotranspiration (AET) has been projected to increase by up to 0.3 (~0.2) mm per year before 2050 (2100). The trends in both projections (PET and AET) are exaggerated under SSP5-8.5. The analysis predicted a shortage of water availability (precipitation—AET), which is alarming for most MENA regions. Relative to the reference period (1981–2010), the decline in annual water availability would reach 26 (62) mm by 2100 under SSP2-4.5 (SSP5-8.5). The rise in temperatures appears to be the principal reason for MENA and water availability responses. This study’s outcomes can facilitate accurate and realistic predictions related to evapotranspiration and water availability, which are key elements in not only managing water resources but also in devising effective climate change mitigation and adaptation plans. Graphical abstract Graphical abstract


Introduction
Evapotranspiration refers to water movement from land to atmosphere through two processes: (1) evaporation from water surfaces and soil and (2) transpiration from plants. The evapotranspiration phenomenon is of high interest in many disciplines, including hydrology, agriculture, energy budget, and climate studies. It is a required component when identifying irrigation demands and runoff quantities (Li et al. 2020), when looking into an energy budget for the earth-atmosphere system (Diak et al. 2004), when improving the predictions of atmospheric circulation models (Pielke et al. 2003), and when undergoing climatic characterization (Wang et al. 2017). Accurate estimations of evapotranspiration enable farmers to apply suitable crop water requirements, which saves water, energy, and chemicals apart from preventing regional environmental contaminations. Erroneous estimations, by contrast, lead to poor irrigation practices (water shortages or waterlogging) that adversely affect crop yields (Rahimi-Moghaddam et al. 2018). The two terms mainly used to describe evapotranspiration are actual evapotranspiration (AET) and potential evapotranspiration (PET). The AET is computed under nonstandard (i.e., real) environmental and management conditions. Its measurement is complex, time-consuming, and costly, primarily because it depends on climatic parameters, plant properties, and environmental conditions (Goya and Harmsen 2014). Meanwhile, PET is calculated under standard conditions-in a uniform and completely vegetated area with a permanent water supply and shading.
Studying trends of evapotranspiration is challenging. No linear relationship can be generated between evapotranspiration and its controlling factors. Evapotranspiration estimation correlates with several influences, such as temperature and rainfall (Tabari and Talaee 2014), relative humidity (Gong et al. 2006), vapor pressure (Irmak et al. 2006), and land use (Eichelmann et al. 2018). The relevant literature, which has forecasted the trend in global evapotranspiration, displays some discrepancies (Anabalón and Sharma 2017;Dong and Dai 2016). For instance, while Wang et al. (2010) Mueller and Seneviratne (2014) reported that 14 CMIP5 models produced a systematic overestimation bias in global evapotranspiration projections of 0.17 mm per day on average. This bias underestimated evapotranspiration values in the Northern Hemisphere.
Accordingly, it is essential to evaluate evapotranspiration behavior under climate change on regional scales-particularly in the Middle East and North Africa (MENA) region, which is highly projected to climatic changes through the twenty-first century. These changes include a significant increase in temperatures and warm days, heat extremes, longevity of drought periods, and precipitation extremes in several regions (Ajjur and Riffi 2020;Lelieveld et al. 2016;Zittis 2017). Climatic changes can intensify evapotranspiration losses and regional water deficits (Erol and Randhir 2012). The MENA region consists predominantly of hyper-arid areas; therefore, several given methods for evapotranspiration estimation, such as Priestley and Taylor's (1972) and Makkink's (1957), are inapplicable (Goya and Harmsen 2014). Other methods, including Blaney and Criddle (1950), Hargreaves and Samni (1982), and Dalton (1802), may be unreliable (Li et al. 2016). A number of reliable methods, such as those that are physically based, demand extensive meteorological data that are not readily available in the MENA region, as those data provide neither sufficient coverage nor good quality (Zittis 2017). Li et al. (2016) evaluated the performance of six models for estimating PET in arid northwest China. They found that Penman (1948) and Shuttleworth (1993) are more reliable and accurate methods.
As such, this study used climatic data from an ensemble comprising 10 models in the CMIP6 to forecast the variations in PET, AET, and water availability across the MENA region throughout the twenty-first century. Two periods were considered for future projections: 2021-2050 (hereafter referred to as the mid-twenty-first century) and 2071-2100 (hereafter referred to as the end-twenty-first century). Two shared socioeconomic pathways (SSPs) were selected to simulate future projections from each model: SSP2-4.5 and SSP5-8.5. The study outcomes are key factors in managing water resources as well as in establishing climate change mitigation and adaptation plans.

Study area
According to the World Bank, the MENA region (12°N-40°N, 12°W-63°E) comprises 19 countries extending from Iran in the east to Morocco in the west (see Fig. 1). The MENA region has several common characteristics in terms of water resources and climate but is diverse in geography, natural resources, economic activities, and land use. Generally, MENA is one of the most water-stressed regions in the world. The average water availability in the MENA region is 1000 m 3 of freshwater per capita per year, whereas the global average exceeds 7000 m 3 (World Bank 2020). Anthropogenic stresses continue to pose pressure on water resources. During the last half-century, the population across MENA quadrupled, while urbanization nearly doubled (World Bank 2020). Most of the MENA region lies under hot desert class (BWh) as categorized by the Köppen climate classification (Peel et al. 2007). Based on an ensemble that comprises 10 CMIP6 models used in this study, MENA's mean near-surface air temperature (hereafter referred to as temperature) ranged between 5 and 31°C, while MENA's annual rainfall varied up to 231 mm from 1981 to 2010.

Data
Maximum and minimum temperature, precipitation, wind speed, radiation, and relative humidity simulations were obtained from 10 Atmosphere-Ocean General Circulation Models (AOGCMs) from CMIP6. Table 1 presents information about these models, including their institutions, countries, and horizontal resolution. The model data were retrieved from the archive of the earth system grid federation online system (https://esgf-node.llnl.gov/search/ cmip6/). Two SSPs were selected to simulate future projections from each model. These reflected the intermediate mitigation scenario (SSP2-4.5) and the very high greenhouse gas (GHG) emissions scenario (SSP5-8.5) as prescribed by the Intergovernmental Panel on Climate Change (IPCC 2014). The SSP2-4.5 (SSP5-8.5) considers a radiative forcing of 4.5 (8.5) Watts per square meter (W m −2 ) and predicts an increase of 1.8 ± 0.7 (3.7 ± 1.1)°C globally before the twenty-second century (IPCC 2014). Unless unavailable, the r1i1p1 ensemble member was selected for each model. The CMIP6 models have different resolutions;  therefore, each model dataset was remapped to a common grid of 1°× 1°using the bilinear interpolation method. Several studies (Hay et al. 2000;Wu and Huang 2016;Wu et al. 2020) reported biases in climate models. They suggested correcting the bias in original CMIP data to evade unreliable climate change assessments. Therefore, we compared monthly MENA-averaged precipitation and temperature AOGCMs simulations with observation-based global datasets. For precipitation, we referred to the Global Precipitation Climatology Centre (GPCC) V2018 datasets. The GPCC data are available at 1°× 1°grids from the National Oceanic and Atmospheric Administration's (NOAA) (https://psl.noaa.gov/data/gridded/data.gpcc.html#detail-last access March 18, 2021). We selected the GPCC precipitation datasets because their reliability have been proven in the Mediterranean region (Nastos 2011). For maximum and minimum temperature, we referred to the Climatic Research Unit -University of East Anglia datasets (CRU TS 4.05). The CRU TS 4.05 data are available at 0.5°× 0.5°grids from the https://sites.uea.ac.uk/cru/-last access March 18, 2021 (Harris et al. 2020). For ensuring the spatial resolution consistency, the CRU TS 4.05 datasets were bilinearly interpolated to the 1°× 1°grids similar to the CMIP6 simulations.

Methods
The first step was to validate and correct the bias in precipitation and maximum and minimum temperatures from CMIP6 simulations according to observational datasets. Next, the corrected data in addition to relative humidity, radiation, and wind speed simulations were applied to estimate PET, AET, and water availability for both the reference and future periods. The world meteorological organization (WMO 2017) recommends the most-recent 30-year period finishing in a year ending with zero as the climatological standard normal. Thus, the period between 1981 and 2010 was selected as the reference period. Finally, the response of evapotranspiration and water availability to climatic changes during the mid-twenty-first century was defined as the difference between evapotranspiration and water availability projections for 2021-2050 and 1981-2010. Similarly, the variations during the end-twentyfirst century were defined as the difference between evapotranspiration and water availability projections for 2071-2100 and 1981-2010. No temporal downscaling was applied in CMIP6 model projections, primarily because the projection period, 30 years, is short (Dezsi et al. 2018). The analyses mentioned above were performed using the Python and ArcGIS software.

Bias correction of CMIP6 simulations
This study used the "delta change" method (Hay et al. 2000) to correct the bias in precipitation and temperature AOGCMs simulations. The delta change method is a common way to correct the bias in GCM simulation data (Wu and Huang 2016;Wu et al. 2021;Wu et al. 2020;Xu et al. 2019). For precipitation, original CMIP6 simulations were multiplied by a correction factor as per Eq. 1: where x corri; j is the corrected precipitation at the j th grid point; x simulated i, j is the simulated precipitation at the j th grid point; and x observed i; j and x simulated i; j are the averages of the observed and simulated precipitation, respectively, at the j th grid point during the reference period . For maximum and minimum temperature, a correction factor was added to original CMIP6 simulations as per Eq. 2:

Potential evapotranspiration estimation
The study estimates the MENA PET using the Penman (1948) method. Penman (1948) is a widely used standard method to estimates PET and evaluate climate change impact on water resources in various locations and climatic conditions (Goya and Harmsen 2014;Li et al. 2013;Li et al. 2016;Wang et al. 2010;Wang et al. 2017;Wu et al. 2018). The Penman (1948) equation estimates PET as follows: where Δ is the slope of vapor pressure curve in kPa/°C; R n is the net radiation in MJ/m 2 .day; γ is the psychrometric constant (kPa/°C); U 2 is the wind speed at 2 m above the ground surface in m/s; RH mean is the mean relative humidity in %; e s is the mean saturation vapor pressure in kPa; and λ is the latent heat of vaporization of water in MJ/kg. The Rn is defined as the difference between incoming and outgoing radiation of both short and long wavelengths. The Rn can be calculated from other meteorological parameters like maximum air temperature T max (°C), minimum air temperature T min (°C), and RH mean . The psychrometric constant (γ) is valued at~0.665 × 10 −3 P, where P is the atmospheric pressure in KPa. The λ equals to 2.45 MJ/kg. Wind speed is simulated at a 10 m height in AOGCMs. Therefore, we used Eq. 4 to obtain the wind speed at 2 m above the surface: where u 2 is the wind speed at 2 m above the ground surface in m/s; u z is the simulated wind speed at z m above the ground surface in m/s; and z is the height of measurement above the ground surface in m. In other words, we multiplied U 10 by 0.748 to get U 2 . For additional information about PET calculation, the reader can refer to Doorenbos (1975); Goya and Harmsen (2014); and Penman (1948).

Actual evapotranspiration estimation
The annual AET was calculated using the method prescribed by Budyko and Miller (1974).
Their original framework is dependent on the long-term water balance principle, where AET is controlled by precipitation and PET. Several Budyko curves were developed later with different functional shapes, but the base hypothesis for all curves is identical: water limits evapotranspiration values in dry environments, while energy limits evapotranspiration values in wet environments (Baw-Puh 1981;Choudhury 1999;Donohue et al. 2012;Yang et al. 2008). The Budyko and Miller (1974) method is valuable for assessing climate change and surface condition impacts on water resources in scarce-data regions (Li et al. 2013;Williams et al. 2012;Yang et al. 2008;Zhang et al. 2001). This study applied PET and precipitation data to estimate the annual AET using the common Budyko Fu formula (Baw-Puh 1981) (Eq. 5): where P and PET are the annual precipitation and annual PET, respectively. Since precipitation data (projected from CMIP6) and PET data (calculated) have monthly time scales, all precipitations and PET outputs were aggregated into years by adding the values of their corresponding months. Parameter "ω" refers to the watershed characteristics parameter, which displays the capacity of the region to retain water for evapotranspiration based on its land use and cover and water-energy conditions. The "ω" parameter is linked with soil, topography, and vegetation properties in addition to rainfall and seasonality characteristics (Li et al. 2013;Williams et al. 2012;Zhang et al. 2001). It is also dimensionless and can be derived from Eq. 6 (Yang et al. 2008), as given in the following: where n is a dimensionless catchment-specific variable. Choudhury (1999) considered a default global value of 1.8 for "n." Other studies (Yang et al. 2009;Yang et al. 2008) reported that "n" values are mostly between 1.5 and 2.6. Li et al. (2013) proposed a range of n for 26 major global river basins (> 300,000 km 2 ) based on the long-term averaged annual vegetation coverage. Li et al. (2013) found that the basin-specific "n" performs better than the default global value of "n." In the Nile river basin, they set a value of "n" as 1.16 (ω = 1.88). As the purpose of this study is specifically to estimate the temporal evolution of AET, we adopted the Budyko parameters fitted in the Nile basin from Li et al. (2013) to represent the MENA region. This study also determined the aridity index (Ψ) (Eq. 7) and water availability (Eq. 8).
Water availability ¼ PÀAET 4 Results Figure 2 shows a comparison between monthly MENA-averaged precipitation and maximum and minimum temperature simulated by AOGCMs with observations during the reference period. Data are shown before bias correction (Fig. 2 a, b, and c) and after bias correction (Fig.  2d, e, and f). The AOGCM simulations, especially precipitation, have clearer discrepancy with the observations compared with the corrected simulations. After the bias correction, the improved multi-model ensemble median (MEM) precipitation agreed well with the GPCC observations. The correlation coefficient r = 0.92. For the maximum and minimum temperature, it is clear that the AOGCMs simulations had smaller discrepancy (compared with precipitation) with CRU TS observations during the reference period. After applying the bias correction, the correlation between the CMIP6 MEM and the CRU TS observations has improved to 0.98. In general, the bias-corrected precipitation and temperature simulations are able to represent the reference period since they are in good agreement with observational datasets. Therefore, the bias-corrected MEM simulations can predict the future changes in AET and water availability in the MENA region. Figure 3 illustrates the future variations of the CMIP6 MEM for mean temperature and precipitation in the MENA region through the twenty-first century. Future variations under SSP2-4.5 and SSP5-8.5 are presented as anomalies in°C for mean temperature and as percentages (%) for precipitation relative to the (1981-2010) average. Under SSP2-4.5, the monthly temperature was, relative to the reference period, projected to increase by 0.31°C by the mid-twenty-first century, and between 1.94 and 2.75°C by the endtwenty-first century. High values of temperature increases (during the end-twenty-first century) were noted in all MENA regions, except northern parts. Meanwhile, annual precipitation projections displayed substantial variations throughout the MENA region. When compared to 1981-2010, the annual precipitation changed from −4 to +32% during the mid-twenty-first century. This trend varied at the end-twenty-first century (see Figs. 2d and 3c). Values of annual precipitation increased (from 32 to 69% during the end-twenty-first century) in the North African countries, particularly south Egypt and Libya. High values of annual precipitation increased (up to 32% during the end-twenty-first century) across southern areas of the Arabian Peninsula, particularly in Yemen, Oman, and south Saudi Arabia.

Temperature and precipitation variations
Under SSP5-8.5, the most prominent climatic feature was a surge in temperature, which grew stronger at the end-twenty-first century. The MENA region would have been exposed to temperature increases between 3.56 and 5.19°C before 2100, with high values over the Arabian Gulf countries and south Algeria. Similar to SSP2-4.5, the CMIP6 ensemble projected that annual rainfall would increase in southeastern areas of the MENA region under SSP5-8.5. This increase would hit 105% in south Egypt and Libya.

Potential and actual evapotranspiration variations
The mean rates of PET throughout the MENA region during both the reference and future periods are displayed in Fig. 4. The PET is projected to rise uniformly across the MENA region. Under SSP2-4.5, the maximum value of annual PET increased from 340 mm for the reference period to 355 mm for the mid-twenty-first century and to 386 mm for the end-twenty-first century. Since the period between 1981 and 2021 (or between 2010 and 2050) is 40 years, the PET increased by up to 0.375 mm per year before 2050. The increasing rate of PET is higher through the end of the twenty-first century. Before 2100, the PET increased by up to 0.51 mm per year. High values of annual PET were noted in the Arabian Gulf region, while lower values of annual PET (~83 mm) were observed in northern Arabian Peninsula-such as Iraq and Iran. The PET losses increased significantly at the end-twenty-first century under SSP5-8.5. The analysis revealed that maximum annual PET is projected to increase by 90 mm during the endtwenty-first century (430-340 mm) relative to the reference period. This increase proposes a pronounced increased rate of PET at 1 mm year −1 .
By implementing the Budyko framework, annual AET values were derived across the MENA region under SSP2-4.5 and SSP5-8.5. Figure 5 illustrates the annual AET values during reference and future periods. Under SSP2-4.5, the annual AET varied up to 109 mm for the reference period, while up to 121 mm for the mid-twenty-first century and up to 126 mm for the end-twenty-first century. In other words, the AET is expected to rise by up to 0.3 (0.19) mm per year before 2050 (2100). High values of annual AET were noted in the north Arabian Gulf, while most other MENA areas had low values of annual AET (< 63 mm under SSP2-4.5) during the two future periods. Very low values of annual AET (< 32 mm under SSP2-4.5) were observed in south Egypt and Libya during the end-twenty-first century, which can be attributed to the precipitation increase projections in these regions (see Fig. 2).
The AET losses increased during the end-twenty-first century based on SSP5-8.5. The analysis showed that the maximum annual AET is projected to increase by 32 mm (141-109 mm) during the end-twenty-first century relative to the reference period. This increase proposes a significant growing rate of AET up to 0.35 mm per year.

Aridity index and water availability variations
Generally, the aridity index (PET/precipitation) is projected to see a slight increase in the twenty-first century (up to 0.15) throughout the MENA region. As portrayed in Fig. 6, the aridity index varied up to 0.48 for the reference period as well as up to 0.62 for the future Variations in the aridity index were not substantial between SSP2-4.5 and SSP5-8.5 during the mid-twenty-first century. Only a slight decrease was noticed under SSP5-8.5 in the Red Sea and south Egypt regions, relative to SSP2-4.5, during the end-twenty-first century (see Fig. 6 c and e). Under the SSP5-8.5, MENA aridity increases, and very low values of the aridity index (< 0.07) dominate most regions in the end-twenty-first century.
The water availability component revealed more spatial variability throughout the MENA region. Figure 7 shows that relative to the reference period, the MENA appeared vulnerable to water deficit with a significant shortage in all countries. Under SSP2-4.5, annual water availability is projected to decrease approximately 26 mm (102-76) through 2100 relative to the reference period in all MENA regions. In other words, annual water availability decreased at 0.29 mm per year under SSP2-4.5. Not many changes were observed between the middle of the twenty-first-century projections for annual water availability under SSP2-4.5 and SSP5-8.5, as Fig. 7 b and d show. This is ascribed to the projected increase in precipitation values in some areas of the MENA region, which increased water availability (precipitation-AET) in these areas. On the other hand, water availability changes were apparent under SSP5-8.5 throughout the end-twenty-first century. A comparison between SSP2-4.5 and SSP5-8.5 in Fig. 7 displayed that, except in east Yemen, the maximum annual MENA water availability, in 2100, decreased from 76 mm under SSP2-4.5 to 40 mm under SSP5-8.5. Considering that the variation period is 90 years, the difference occurred between 1981 and 2071 (or between 2010 and 2100). Therefore, on an annual scale, water availability decreased an extra 0.4 mm year −1 under SSP5-8.5, compared to SSP2-4.5.

Summary and discussion
The present study investigated the evolutions of evapotranspiration and water availability throughout the twenty-first century across the MENA region. The analysis under SSP2-4.5 revealed that PET and AET would increase by up to 0.37 mm per year during the mid-twentyfirst century . The PET rate would further increase by approximately 136% during the end-twenty-first century (2071-2100). The increasing trend in PET and AET over the MENA region is in line (to some extent) with that reported in the literature during past decades. Anabalón and Sharma (2017) noticed that the global PET (AET) increased at 1.09 (0.54) mm per year during 1980-2012. The trend reported by Anabalón and Sharma (2017) was clearly represented in water-limited regions like the MENA. Jung et al. (2010) and Wang et al. (2010) suggested higher values of global AET increments. Jung et al. (2010) found that global AET increased by 7.1 ± 1.0 mm per year during 1982-1997, whereas Wang et al. (2010 suggested an AET increasing rate of 15 mm per year during 1982-2002. In this study, the maximum decrease in MENA annual water availability would hit 26 mm in 2100 under SSP2-4.5, and the latter value was 62 mm under SSP5-8.5. These findings are consistent with those of McCabe and Wolock (2015), who demonstrated that the region with latitude band 0-40°N had experienced the most considerable decrease in water availability globally between 1905 and 2009. Another study in Iran (Tabari and Talaee 2014) found that arid regions experienced the greatest change in PET. This present analysis showed that global warming would keep pressurizing water deficits in the MENA region through the twenty-first century.
During the mid-twenty-first century, slight changes were noted between SSP2-4.5 and SSP5-8.5 in the heat, PET, AET, and aridity indices. On the contrary, the variations were larger (and clearer) in all terms (except the aridity index) during the end-twenty-first century between . This is compatible with the powerful atmospheric demand control on evapotranspiration. Sullivan et al. (2019) reported that the global trend in AET was 0.17-0.34 mm per year, which appeared greater in SSP5-8.5 relative to . This present analysis found that the AET trend was up to 0.17 mm per year, greater in SSP5-8.5 relative to SSP2-4.5 (see the differences in Fig. 5 c and e). It can be concluded, then, that arid environments are attributed to greater losses in AET than other climates.
Prior work had reported different parameters that affected PET estimations. Gong et al. (2006) found that the most sensitive variables to PET in the Yangtze River basin (China) were in the following order: relative humidity ≻ shortwave radiation ≻ air temperature ≻ wind speed. Irmak et al. (2006) discovered that daily PET in different climatic conditions (semiarid, Mediterranean, and humid) across the US appeared to be sensitive mainly to vapor pressure deficits followed by wind speed. Using the physically based method (Penman 1948) helped considering the most important variables that govern the evaporation process in the MENA region. As for AET and water availability, their sensitivities can be linked mainly to temperature rises across the whole MENA area. Through the twenty-first century, the MEM of 10 CMIP6 models projected temperature increases of up to 2.75 (5.2)°C under SSP2-4.5 (SSP5-8.5) while projecting variant trends in rainfall direction and magnitude across the MENA region. However, AET increased and water availability decreased relatively uniformly throughout all MENA areas. This suggests that temperature upsurge is the principal reason for evapotranspiration and water availability responses in arid areas. The latter finding is in line with Tabari and Talaee (2014) who reported that AET is primarily driven by temperature trends in arid regions across Iran. The situation may differ in humid regions where AET is influenced by different variables. For example, Anabalón and Sharma (2017) suggested that AET in the Amazon Basin and central Europe was affected largely by vegetation cover, and Sullivan et al. (2019) revealed that the AET's sensitivity was due primarily to the temperature trend over North America.
This work may carry some uncertainties. First, PET and AET are complex phenomena and may not be represented ideally using the Penman (1948) and Budyko and Miller (1974) formulas. The literature has suggested larger uncertainties in PET estimations. Although using the physically based reliable Penman (1948) method incorporates more variables and hence bridge some gaps, local calibration of these variables over the MENA might still be required. Regarding AET estimation, Wang et al. (2018) and Feng et al. (2019) concluded that the Budyko approach can better predict AET in non-humid areas (compared to humid areas). Hence, the uncertainty in AET estimation is smaller than that of PET. Second, the uncertainty surrounding GHG emission scenarios and the CMIP6 models should lead to subsequent uncertainty in future projections of evapotranspiration and water availability throughout the MENA region. Previous literature found larger uncertainties in CMIP models compared with emission scenarios (Wu et al. 2021;Xu et al. 2019). In this study, the original CMIP6 simulations were not able to accurately project the climatic changes across the MENA through 2100-especially rainfall. Rainfall is a major component in the Budyko framework, meaning that small inaccuracy in rainfall projections could lead to inaccuracy in AET estimations. Using an ensemble of the CMIP6 models as well as correcting the bias as proposed in this study helps bridge the gap. Furthermore, this study's aim was not to give ground truth estimations at local scales but to estimate the evapotranspiration and water availability evolutions through the twenty-first century.

Conclusion
Integrating the Penman (1948) and Budyko and Miller (1974) approaches with temperature, precipitation, wind speed, radiation, and relative humidity simulations in CMIP6 model projects changes in evapotranspiration and water availability (through 2100) across the MENA region. Comparing future periods (through the twenty-first century) with the historical (reference) period (1981-2010) offers guidance to policymakers and scientists for establishing climate change adaptation and mitigation measurements apart from managing water resources effectively. The temperature surge over the MENA region seemed to increase PET losses. The increase in PET amplified AET but declined water availability significantly during the twentyfirst century. This implies substantial changes in components of the hydrological cycle, which increased the imbalance between water supply and demand across the MENA region (Ajjur and Baalousha 2021). The study's outcomes are valuable, especially when ground data are unavailable. However, further regional work is warranted, especially in humid areas of the MENA region. This work should be followed by an assessment of how the projected changes in evapotranspiration can affect water resources, land cover, and soil moisture as well as runoff at regional scales. Such knowledge helps to increase the resilience of the MENA region's fragile systems against climatic changes.

Conflict of interest The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.