Near-term regional climate change in East Africa

In the coming few decades, projected increases in global temperature and humidity are generally expected to exacerbate human exposure to climate extremes (e.g., humid-heat and rainfall extremes). Despite the growing risk of humid-heat stress (measured by wet-bulb temperature), it has received less attention in East Africa, where arid and semi-arid climatic conditions prevail. Moreover, no consensus has yet been reached across models regarding future changes in rainfall over this region. Here, we screen Global Climate Models (GCMs) from CMIP5 and CMIP6 and use, for boundary conditions, simulations from only those GCMs that simulate successfully recent climatic trends. Based on these GCMs and Regional Climate Model (RCM) simulations, we project that annual mean temperature is likely to rise by 2 ℃ toward midcentury (2021–2050) at a faster rate than the global average (about 1.5 ℃), under the RCP8.5 and SSP5-8.5 scenarios, associated with more frequent and severe climate extremes. In particular, low-lying regions in East Africa will be vulnerable to severe heat stress, with an extreme wet-bulb temperature approaching or exceeding the US National Weather Service’s extreme danger threshold of 31 ℃. On the other hand, population centers in the highlands of Ethiopia will receive significantly more precipitation during the autumn season and will see more extreme rainfall events, with implications for flooding and agriculture. The robustness of these results across all GCM and RCM simulations, and for both of CMIP5 and CMIP6 frameworks (CMIP: Coupled Model Inter-comparison Project) supports the reliability of these future projections. Our simulations of near-term climate change impacts are designed to inform the development of sound adaptation strategies for the region.


Introduction
With global warming accelerating, climate extremes, such as heatwaves and heavy rainfall, have attracted much attention in past decades due to their impacts on natural and human systems (Meehl and Tebaldi 2004;Kharin et al. 2013;Sillmann et al. 2013;Fischer and Knutti 2015). The risks from climate change vary significantly depending on The original online version of this article was revised: In this article some text was changed under the section Funding and acknowledgement.
Elfatih A.B. Eltahir eltahir@mit.edu the main focus of this paper. These countries have recently experienced an upward trend as well as increasing variability in rainfall, leading to more frequent floods and droughts (Masih et al. 2014;Siam and Eltahir 2017;Gebrechorkos et al. 2019). For example, a heavy rainfall event in 2019 with severe flooding affected an estimated 2.8 million people across the region (Wainwright et al. 2021a). Without aggressive mitigation efforts, the increasing population in East Africa will face a severe climate crisis in the next few decades. Thus, it is necessary to produce reliable near-term climate change projections to help policymakers design effective and well-informed adaptation policies.
Simulating the East African climate is, however, a great challenge for Global Climate Models (GCMs) due to the region's complex topography and diverse climate conditions, ranging from semi-arid to tropical monsoon (Flato et al. 2013). Indeed, most state-of-the-art global climate models in the CMIP archive tend to exhibit significant systematic errors in their representation of the climate system in East Africa (e.g., Daron 2014;Ayugi et al. 2021;Wainwright et al. 2021b). These errors could be caused by several factors, such as coarse spatial resolution and inadequate representation of physical processes (Seneviratne et al. 2012;Flato et al. 2013;Im et al. 2017b;Finney et al. 2019;Finney et al. 2020;Wainwright et al. 2021b), deficiencies in simulating the African monsoon (McSweeney et al. 2015), and poor representation of ENSO and its teleconnection to Africa (Rowell 2013). Indeed, future climate projections using GCMs provide insufficient evidence for the change in rainfall in this region (Conway and Schipper 2011;Park et al. 2015;Daron 2014;Nangombe et al. 2018;Bichet et al. 2020). These substantial uncertainties make it difficult to  CIESIN, 2018). Blue dashed line in (a) indicates the Nile River. Red crosses in (a) denote the location of 16 weather stations in East Africa for the use of daily maximum temperature, daily mean dew point temperature, and surface pressure. Red crosses and circles in (b) indicate two stream flow gauge stations and five precipitation stations, respectively determine the impacts of climate change which are needed to inform and support the development of sound adaptation strategies.
In contrast, Regional Climate Models (RCMs), with finer resolution and adequate representation of physical processes, can mitigate the known problems in GCMs (Giorgi et al. 2009;Giorgi and Gutowski 2015). Thus, RCMs are appropriate tools for assessing regional climate change for a region with complex topography and land use distribution, like Africa (e.g., Laprise et al. 2013;Im et al. 2014;Alter et al. 2015;Giorgi and Gutowski 2015;Im and Eltahir 2018a). Recent studies through the Coordinated Regional Downscaling Experiment for Africa (CORDEX-AFRICA) have constructed regional climate change information for a pan-Africa domain. However, their horizontal resolution (50 km) still has limitations for use in climate impact studies.
So far, most climate change studies over East Africa using RCMs confine their attention to the end of the 21st century (2071-2100), or to a specific period in which the global mean temperature rises by 1.5 ℃ or 2 ℃ above pre-industrial levels (e.g., Dosio 2017;Déqué et al. 2017;Osima et al. 2018;Nikulin et al. 2018;Dunning et al. 2018). Moreover, only a few studies have attempted to examine the anthropogenic influences on weather extremes, such as heavy rainfall and heatwaves in the near-term future (Russo et al. 2016). Existing studies of heatwaves are primarily based on maximum temperature and do not consider the humidity effect (Engelbrecht et al. 2015;Russo et al. 2016;Dosio 2017), while heat stress impacts are significantly mediated by humidity (Pal and Eltahir 2016;Im et al. 2017a). This study considers CMIP5 GCMs, CMIP6 GCMs, and the MIT Regional Climate Model (MRCM; Pal and Eltahir 2016) to project near-term climate change and its impacts on heat stress and rainfall over East Africa.

Data collection
In-situ rainfall measurements at Addis Ababa, Debre Markos, Dessie, Gondar, and Nekemet from the Global Historical Climatology Network (GHCN; Peterson and Vose 1997) (Fig. 1b) are used to elucidate recent trends in rainfall in the Upper Blue Nile basin. In-situ runoff data at Diem and Atbara stations (Fig. 1b) are provided through personal communications with the ministries of water in the countries of the Nile basin. Various gridded precipitation datasets are obtained from the Climate Research Unit product (CRU; Harris et al. 2020), the Global Precipitation Climatology Project (GPCP; Adler et al. 2003), the Integrated Multi-Satellite Retrievals for the Global Precipitation Measurement mission (Huffman et al. 2019), the Tropical Rainfall Measuring Mission (TRMM; Huffman et al. 2007), and the Climate Hazards Group InfraRed Precipitation with Station (CHIRPS; Funk et al. 2015).
Dry-bulb temperature (T), dew point temperature, and surface pressure data are taken from the ERA5 reanalysis (Hersbach et al. 2020) at 3-hour intervals on 0.25° × 0.25° horizontal resolution for the period 1976-2019. Daily mean and maximum wet-bulb temperatures (TW and TWmax) (Davies-Jones 2008) are calculated using the ERA5 data. To quantify level of human exposure to extreme heat, the categories of heat stress, defined by the US National Weather Service (USNWS), are applied. Note that a wet-bulb temperature of 31 ℃, sustained for 6 h, is assumed for an extreme danger threshold of humid heat stress for various regions (Im et al. 2017a;Im et al. 2018;Choi et al. 2021).

Selection of gridded precipitation data
It is important to use accurate estimates of precipitation to assess the model performance and ultimately provide more reliable climate change projections. However, gridded precipitation datasets have limitations due to quality issues (Alaminie et al. 2021), especially in African countries where observational sites are spatially sparse (Sun et al. 2018).
To select a gridded precipitation data set with the best quality standards for East Africa, various precipitation estimates are compared to the in-situ observational data which provide relatively accurate measurements (Fig. 2). The decision on data selection is based primarily on its ability to represent recent climate trends. In modern records, precipitation at five stations in the Upper Blue Nile basin shows an overall increasing trend ( Fig. 2a), consistent with the finding of Siam and Eltahir (2017). In addition, runoff at Diem and Atbara stations supports these rainfall trends over the region (Fig. 2b). Contrary to the observed trends, drying trends are evident in two precipitation estimates, including GPCP and TRMM (Fig. 2c). By contrast, the CRU, GPM and CHIRPS data exhibit similar trends to the in-situ observations (Fig. 2c). The CHIRPS, which has a 0.05° resolution and exhibits trends consistent with those of the station observations, is mainly used to evaluate climate simulations rather than CRU (which has a coarser spatial resolution (0.5° × 0.5°)) and GPM (which is limited to a relatively short period (2001-2019)).

Selecting GCMs for downscaling over East Africa
As is well known from previous observational and modelling studies, the signal of anthropogenic climate change was detected in recent decades (e.g., IPCC 2013). In particular, Park et al., 2016;Giannini and Kaplan 2019). This suggests that GCMs, which reasonably well simulate not only climatic conditions but also recent climate-change trends, can provide reliable future climate projections.
the Sahel zone of Sudan and northern parts of Ethiopia have experienced notable climatic shift from dryer to wetter conditions since the 1980s (Fig. 2), primarily due to escalating greenhouse gas emissions (e.g., Dong and Sutton 2015; and CNRM-CM5) judged to be satisfactory in their representation of key aspects of the African climate (i.e., African monsoon circulation, climatological annual cycles of temperature and precipitation, the influence of teleconnection patterns, such as El Niño -Southern Oscillation and the Indian Ocean Dipole). Among the six models, we then identify that MRCM driven by three models (GFDL-CM3, HadGEM2-ES, and NorESM1-M) reproduce the observed climate trends (Figs. 3, S1, and S2). Meanwhile, MRCM forced by CMIP5 GCMs, including ACCESS1-0, CNRM-CM5, and MPI-ESM-MR, perform poorly ( Fig. 3) and are therefore eliminated from further analysis. As expected, climatological features are well captured by the three selected CMIP5 GCMs (Fig. S3).
However, most existing studies have placed emphasis on the ability of a GCM to describe the climatological mean features (Baek et al. 2013;McSweeney et al. 2015). In this study, among the many options in the CMIP5/CMIP6 data archive, a set of 6 GCMs (three from the CMIP5 archive and three from the CMIP6 archive) are carefully selected through a rigorous screening process. We mainly consider the model's ability in capturing not only average climatic conditions, but also recent climate trends in East Africa.
To avoid duplicating existing efforts, the selection process for CMIP5 GCMs relied on previous assessment performed by McSweeney et al. (2015). We, followed their suggestion, first select six CMIP5 GCMs (GFDL-CM3, Had-GEM2-ES, NorESM1-M, ACCESS1-0, MPI-ESM-MR, Project Phase 5 (CMIP5) GCMs is smaller than that for the MRCM simulations forced by CMIP6 GCMs. The horizontal resolution is 20 km in both longitudinal and latitudinal directions. As explained above, to specify the boundary conditions for MRCM, we carefully select three GCMs (GFDL-CM3, HadGEM2-ES, and NorESM1-M) from the CMIP5 archive and three GCMs (CMCC-ESM2, Had-GEM3-GC31-MM, and NorESM2-MM) from the CMIP6 archive (see Sect. 2.3). We perform two 30-year time slice simulations. The first one starts in 1975 and the second time slice experiment starts in 2020. For the MRCM simulations driven by CMIP5 GCMs, time-lagged ensemble members are produced with different starting dates (i.e., January 1st, 2nd, and 3rd in 1975, and January 1st, 2nd, and 3rd in 2020) ( Table 1).
To provide more reliable climate projections, it is required to eliminate systematic biases in the MRCM output. Following previous studies (Li et al. 2010;Choi et al. 2021), we apply the equidistant quantile-mapping bias correction procedure under the assumption that the bias does not change over time. Compared to the widely used quantile-mapping technique, this method explicitly considers changes in distribution between reference and projection periods. The ERA5 and CHIRPS data, which are regarded as relatively unbiased, are used to correct the bias in simulated monthly mean temperatures and precipitation (PR) over the entire domain, respectively. To correct the systematic bias apparent in model simulations at the station scale, we consider daily maximum dry-bulb temperature, dewpoint temperature, and surface pressure data derived from the GHCN at 16 stations (i.e., large population centers in Sudan and Ethiopia; Fig. 1) over East Africa. However, precipitation data over these stations are of limited use due to a large number of missing values. Thus, CHIRPS is applied to correct the precipitation bias at the 16 stations.
The six GCMs (CMIP5 GCMs: GFDL-CM3, HadGEM2-ES, and NorESM1-M; CMIP6 GCMs: CMCC-ESM2, Had-GEM3-GC31-MM, and NorESM2-MM) are finally selected in this study to specify the boundary conditions for MRCM. It is worth noting that the MRCM ensembles outperform most individual simulations due to the cancellation of some systematic biases.

Regional climate model and experimental setup
The MRCM is a three-dimensional hydrostatic numerical model, of which general structure and dynamical core are rooted in the International Center for Theoretical Physics (ICTP) Regional Climate Model version 3 (RegCM3; Pal et al. 2007). MRCM is notably coupled with the Integrated BIosphere Simulator land surface scheme (IBIS; Foley et al. 1996;Winter et al. 2009) and employs significant enhancements of model physics (Marcella 2012;Marcella and Eltahir 2012;Gianotti 2012;Gianotti andEltahir 2014a, 2014b). The latest version of MRCM well represents land surface processes and has sufficient skill to simulate the fine-scale structure of climatic variables over various domains, such as North America (Winter et al. 2009), West Africa (Im et al. 2014;Im and Eltahir 2018a), Southwest Asia (Pal and Eltahir 2016;Eltahir 2022a, 2022b), South Asia (Im et al. 2017a;Choi et al. 2021), and the Maritime Continent (Im and Eltahir 2018b).
The model domain cover East Africa on a Lambert conformal projection (Fig. 1a). Note that the domain for the MRCM simulations driven by Coupled Model Inter-comparison exhibit dry and wet biases in South Sudan and the highlands of Ethiopia, respectively. For dry-bulb temperature, the warm bias appears in the vicinity of South Sudan (possibly due to the lack of representation of swamps), while the cold bias is evident in the lowlands of Ethiopia and the northern part of the study area. Consistent with these cold biases, there is a significant underestimation of TW over Sudan and the lowlands of Ethiopia. Relative to MRCM/CMIP5, such structural biases are greatly reduced in MRCM/CMIP6, which could be due to the remarkable improvements in CMIP6 models (e.g., spatial resolution, physical processes, and biogeochemical cycles; Eyring et al. 2016). It is also noticeable that MRCM/CMIP6 overall shows better performance in representing the seasonal cycles of PR, T, and TW in East Africa, than GCMs and MRCM/CMIP5 (Fig. S3).

Evaluation of the historical simulations
The performance of MRCM/CMIP6 is superior overall compared to the GCMs and MRCM/CMIP5 in capturing the present climate conditions over East Africa (Fig. 4, S2, and S3). Figure 4 shows the spatial distribution of annual mean PR, T, and TW derived from the observations and MRCM ensembles. In the observational record, this region has diverse climate conditions from semi-arid in the north to tropical monsoon in the south (Fig. 4a). Also, a significant temperature gradient exists between cold regions in the highlands of Ethiopia and warm regions in the lowlands of southeastern Sudan (Fig. 4d, g). These observed features are well captured by MRCM/CMIP5 and MRCM/CMIP6, despite of noticeable biases in some regions. Compared to the precipitation observations, both the MRCM ensembles TWmax approaching or exceeding the USNWS extreme danger threshold of 31 ℃ by the middle of the 21st century under SSP5-8.5 scenario. Notably, Port Sudan, a city close to the Red sea, stands out due to its high TWmax values which are well above 31 ℃ and close to a human survivability threshold of 35 ℃ (Sherwood and Huber 2010;Pal and Eltahir 2016;Im et al. 2017a). The populous cities in Ethiopia are more likely to see a significant rise in TWmax, but the projected TWmax will be far below the USNWS extreme danger level, except for Gode. Though the projected risk of extreme heat stress will be generally lower in MRCM/CMIP6 than in MRCM/CMIP5, there is a robust inter-model consistency across different generations of models (Figs. 7 and S8).
The return period analysis further confirms the elevated risks of heat stress in a warmer world (Fig. S9). For the historical period, a recurrent period of extreme heat events (above the extreme danger threshold of 31 ℃) is close to or far longer than 100 years, depending on the given region. Though some discrepancy exists in different generations of CMIP experiments, significant decreases in the return period are consistently expected among different highemissions scenarios (RCP8.5 and SSP5-8.5). For instance, extreme heat events would occur once in less than 50 years over two major cities (Port Sudan and El Obeid) in Sudan in the near-term future. In particular, the most vulnerable city to heat stress is Port Sudan with a recurrent period of less than 5 years.

Projected changes in near-term precipitation
Consistent with the observed trend of rainfall in recent decades (Figs. 2 and 3), near-term climate projections exhibit an overall increase in precipitation over East Africa (Fig. 8). Recent observation suggests that annual mean precipitation over the region increased by about 9% between the periods 1976-1997 and 1998-2019. Based on MRCM/CMIP6, annual mean rainfall over East Africa is projected to increase by about 7% under the SSP5-8.5 scenario (Fig. 8 h). The consistency of the rainfall changes across MRCM/CMIP6 simulations is evident in Ethiopia. Similarly, MRCM/CMIP5 also projects an overall increase (more than 10%) in precipitation, especially in the Upper Blue Nile. On the other hand, there is limited inter-model agreement on the changes in precipitation over southwestern Sudan (Fig. 8 h). Nevertheless, overall increases in rainfall are expected in the southern part of Sudan in most of the MRCM simulations, except MRCM simulation driven by CMCC which presents a small and insignificant decrease in precipitation.
The seasonal cycle of rainfall in Ethiopia is likely to amplify by midcentury, consistently across CMIP5, CMIP6,

Projected changes in near-term temperature
MRCM ensemble projections (MRCM/CMIP5 and MRCM/ CMIP6) show significant warming in East Africa with a higher inter-model agreement. Figures 5 and 6 present the projected change in T and TW by the middle of the 21st century (2021-2050) under high-emissions scenarios (RCP8.5 and SSP5-8.5). These changes are, to a large degree, spatially uniform. The overall increases in T are evident across all simulations with a slightly greater increase in arid regions of Sudan (above 2.0 ℃; Fig. 5l, x). The change in T can lead to a substantial rise in TW (approximately 1.4 ℃ in MRCM/CMIP5 and 1.3 ℃ in MRCM/CMIP6). The more than 1.0 ℃ increase in TW is expected to increase the level of heat stress over the southeastern border of Sudan and the lowlands of Ethiopia where humid heat stress is already high in the current climate (Fig. 6). Although all MRCM simulations consistently predict the increase (more than 1.0 ℃) in TW in Ethiopia, the projected TW over a mountainous region (elevation higher than 1500 m) is much lower than the USNWS extreme danger threshold of 31 ℃ for the future climate.
A substantial increase in daily maximum temperature (Tmax) is likely to occur in the near future (Fig. S5), along with more frequent and severe temperature extremes in East Africa (Figs. S6 and S7). Fig. S5 shows the projected changes in the annual cycle of Tmax for major cities in the study region. MRCM/CMIP6 projected that the high values of Tmax approaching 45 ℃ mostly appear in central and eastern Sudan, including Port Sudan, Dongola, Atbara, Khartoum, Kassala, and Gedaref in the near future. By contrast, temperature extremes in Ethiopia will remain much below the levels of Sudan, despite the substantial increase in Tmax. Overall, these increases are not confined to a particular region and season (i.e., homogeneous warming), which implies that there would not be a significant change in the seasonality of extreme temperatures. Similar results are obtained with MRCM/CMIP5. It is important to note that the increases in Tmax would not necessarily result in a proportional increase in humid heat stress, since the high risk of heat stress is not just a function of temperature, but a combined consequence of temperature and humidity (Sherwood and Huber 2010;Pal and Eltahir 2016;Im et al. 2017a;Kang and Eltahir 2018).
The MRCM ensembles project that most regions at low elevations in East Africa are expected to be susceptible to humid heat extremes, as global warming progresses (Fig. 7). In the current climate, TWmax exceeding the extreme danger threshold of 31 ℃ is rarely observed. Based on MRCM/ CMIP6, the low elevation regions (Port Sudan, Atbara, Khartoum, Kassala, Gedaref, El Obeid, Nyala, and Gode) are expected to be exposed to frequent heatwaves with a  statistically insignificant. By contrast, Ethiopia may receive significantly more precipitation (p-value < 0.05) during the late summer and autumn seasons in the middle of the 21st century. In particular, all climate projections, regardless of the CMIP phase, expect significant increases in autumn rainfall over the highlands of Ethiopia. Note that the robustness of the wetting trends across all climate simulations (CMIP5 GCMs, CMIP6 GCMs, MRCM/CMIP5, and MRCM/CMIP6) suggests the reliability of future projections in Ethiopia. and MRCM simulations (Fig. 9). Figure 9 shows the projected changes in the annual cycle of PR under the high Green House Gas (GHG) emissions scenarios (RCP8.5 and SSP5-8.5). To identify the spatially heterogeneous effects of climate change, four sub-regions are considered where the precipitation response shows different behavior in climate projections (i.e., western Sudan, eastern Sudan, the highlands of Ethiopia, and the lowlands of Ethiopia). The MRCM ensembles generally project smaller precipitation changes in Sudan than GCM ensembles. Moreover, based on MRCM/CMIP6, future change in rainfall in this region is Large increases in precipitation with an increasing trend (p-value < 0.01) are projected to occur in the highlands of Ethiopia. The consistency between the GCM and MRCM simulations in this region is robust. It is important to note that this wetting tendency in autumn precipitation is continuously found in observations and projections, implying a century-long wetting trend that appears to be proportional to atmospheric GHG concentrations (Fig. S10).
The long-term (1976-2050) trend of precipitation helps to better discern the anthropogenic climate change signal. The evolution of seasonally (September to November) averaged precipitation over the sub-regions in East Africa is shown in Fig. 10. As expected from Fig. 9, MRCM simulations show smaller precipitation changes than their GCM counterparts. The MRCM ensembles exhibit an increase in precipitation by 16 to 22% (5 to 6%) in Sudan and by 5 to 8% (12 to 17%) in Ethiopia according to the RCP8.5 (SSP5-8.5) scenario. Analyses of the return period reveal the increased risk of heavy rainfall over the next few decades in East Africa. Fig. S11 shows the return periods of annual maximum precipitation in the most populous cities over the target region. Although the return periods of extreme precipitation do not dramatically change over most regions, it overall shortens, except for Dongola, Kassala, and Combolcha. In MRCM/ CMIP6 under the SSP5-8.5 scenario, a notable decrease Fig. 9 The 30-year mean seasonal cycle of rainfall for four subregions in East Africa derived from bias-corrected ensemble means of 25 CMIP5 GCM simulations (dashed lines in a-d; see Table S1), 29 CMIP6 GCM simulations (solid lines in a-d; see Table S1) Table S1), 29 CMIP6 GCM simulations (red lines in a-d; see Table S1) near-term climate projections, though still with some systematic biases and uncertainties.

Summary and conclusions
Based on the high-resolution MRCM simulations, near-term climate change would significantly affect the large population centers in Sudan and Ethiopia (i.e., East Africa) under very high emissions scenarios (RCP8.5 and SSP5-8.5). The projected regional climate change features a general warming of mean temperatures with a magnitude of about 2.0 ℃ across the region, which further results in more frequent and severe climate extremes. We project, however, considerably heterogeneous impacts of climate change across the regions. For example, the low plain region in East Africa (Port Sudan, Atbara, Khartoum, Kassala, Gedaref, El Obeid, Nyala, and Gode) is particularly vulnerable to heat extremes in the near future with a TWmax approaching or exceeding the USNWS extreme danger threshold of 31 ℃. Notably, the most vulnerable city to heat stress is Port Sudan, a city close to the Red sea.
In addition to the temperature increase, large population centers in East Africa are likely to experience small but significant increases in precipitation. Furthermore, all climate projections show a robust increase in autumn rainfall over the highlands of Ethiopia with more frequent and strong rainfall extremes. That is, the recurrence interval of extreme rainfall will decrease notably in Bahir Dar, Bole, and Awassa where a 1-in-50-year extreme rainfall event in the reference climate will become an approximately 1-in-25-year event in the near future under the RCP8.5/SSP5-8.5 scenario.
Our study indicates that East Africa would inevitably see exacerbated risks of extreme climate events, like dangerous heat and extreme rainfall, within a few decades without urgent climate change mitigation strategies. These heat extremes make the present semi-arid environments even harsher, contribute to food insecurity, and pose a significant public health challenge. Of particular note, most people in this region who are primarily engaged in agriculture will be vulnerable to these heatwaves. Moreover, enhanced rainfall extremes may cause severe flooding along the Blue Nile river with negative effects of agriculture. Our results identify specific regional hotspots of the near-term climate change and could help guide the development of a sound adaptation strategy.

Discussion
It has recently been suggested that the most aggressive scenarios, RCP8.5/SSP5-8.5, may become increasingly implausible, while approaching the end of the 21st century, assuming significant mitigation efforts (Hausfather and Peters 2020). However, these scenarios are possible, and can still help policy makers design sound adaptation policies, particularly for the near future (Schwalma et al. 2020).
Our results are not directly comparable with the previous regional climate studies that focus on the end of the 21st century or a specific period in which the global mean temperature rises by 1.5 ℃ or 2 ℃ above pre-industrial levels.
Although there are significant limitations to the comparison, our findings are broadly consistent with existing studies (Dosio 2017;Déqué et al. 2017;Osima et al. 2018;Nikulin et al. 2018;Dunning et al. 2018). For instance, based on CORDEX-Africa regional climate models, Déqué et al. (2017) investigated the consequences of a 2 ℃ global warming and showed that extreme heat waves become much more frequent. They also suggest that extreme precipitation could increase, although the change in rainfall is uncertain. In a similar way, Nikulin et al. (2018) and Osima et al. (2018) highlighted a pronounced summer warming and enhanced wet-season precipitation over East Africa in a warmer climate. However, contrary to our results, Mengistu et al. (2021) indicate a decrease in annual mean precipitation in the Ethiopian highlands by the end of this century, while most CORDEX-Africa models disagree on this drying trend (Dosio et al. 2019;Bichet et al. 2020).
This study naturally presents several important caveats. Although we found that CHIRPS data can capture the observed trend of rainfall in East Africa in recent decades by comparing to reliable in-situ station datasets (precipitation and runoff), caution must be taken in interpretation of the CHIRPS data as it contains some estimated values, not all true values. Complete comparison between the CMIP5 and CMIP6 scenarios is impossible because CMIP5 RCPs are different from the SSP scenarios, which contain regional forcings (land use and aerosols). For climate modeling, parameterized convection models are used in this study, which may be a source of uncertainty in the precipitation projections (Finney et a. 2019;Finney et al. 2020;Wainwright et al. 2021b). Also, regional climate projections unavoidably suffer from a significant body of uncertainties, including emission scenarios, model configurations, and regions of interest (Giorgi et al. 2009). Nevertheless, we believe that a well-validated MRCM can give reliable Data availability All data needed to evaluate the conclusions in the paper are present in the paper. Additional data related to this paper may be requested from the authors.
Code availability All the codes are programmed by NCAR Command Language. The codes are available and maintained by Yeon-Woo Choi (choiyw@mit.edu).

Conflicts of interest/Competing interests
The authors declare no conflict of interest.
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/.