There are ample direct and proxy observational evidences that global climate is changing on a wide range of time scales. However, the changes since pre-industrial times are attributed to anthropogenic activities according to various IPCC reports. In particular, the use of fossil fuels led to increase in carbon dioxide and other greenhouse gases in the atmosphere that altered the global energy balance (Allan & Soden, 2008). However, the observed climate changes in response to this energy imbalance are not uniform spatially across the globe. For example, Kenya, from East Africa (EA), experienced rise in temperature on average by about 0.21o C per decade from the 1960s to 2006 (Hill, 1968). This rise in temperature over northern Kenya from October to February is coupled with a decrease in rainfall from long rains from March to May unlike that of southern Kenya (Stiebert et al., 2012). Consistent with this historical trend, climate projections for the northern Kenya have shown an increase in the incidence of drought, high temperature and water scarcity (Measham & Lumbasi, 2013; Muhati et al., 2018).
The spatial heterogeneity in climate change over EA is considerably high. The high rainfall variability and differences in trend over EA are attributed to the complex topography that allows varied local feedback to the radiative forcing (Yang et al., 2014; Omondi et al., 2014; Ntwali et al., 2016). In addition to the varied local response to climate change within EA, there are also changes in climate that are distinct and common to the whole region. For instance, the region receive mean total precipitation lower than the rest of the equatorial regions (Omondi et al., 2014). Moreover, the region regularly suffered from recurrent extreme weather conditions (Lyon & Dewitt, 2012; Liebmann et al., 2014) and has experienced wetting/drying trends in short (long) rains respectively (Lyon & Dewitt, 2012; Yang et al., 2014; Bahaga et al., 2015). Most models have also suggested that there is an increase in temperature and rainfall over EA under various climate change scenarios in the 21s century (Williams & Funk, 2011; Jacob et al., 2012) under various climate change scenarios. However, the increase in projected rainfall under the various scenarios is inconsistent with what is observed over EA. Moreover, the underlying reason for this difference in rainfall trend over EA between projection and observations is not yet fully understood and referred to as the EA climate paradox (Wainwright et al., 2019) and references therein).
Despite the limitations in the skills of models, GCM simulations forced by specified variations in GHGs are commonly used to understand future climate change. For example, using ensemble of observations, reanalysis and simulations data set, (Wainwright et al., 2019) have indicated the rainfall decline in the historical period is, to some extent, related to a later onset and earlier cessation of the long rains. Therefore, there is a concerted effort to improve and use climate models. Coupled Model Intercomparison Project Phase 5 (CMIP5) (Taylor et al., 2011; Krishnan et al., 2019) is one of these international efforts. CMIP5 is driven by historical forcing for current climate and Representative Concentration Pathway (RCP) for projection under different scenarios (Allen et al., 2014). Each Pathway embedded a set of assumptions that lead to four levels of radiative forcing of 2.6, 4.5, 6.0, and 8.5 W /m2, which is labeled as RCP2.6, RCP4.5, RCP6.0, and RCP8.5 respectively. The CMIP5 models are diverse in terms of model complexity, spatial and temporal resolutions. Yet, even the most complex and spatially highly resolved GCMs are not sufficient to resolve local scale processes. To overcome this limitation, dynamical downscaling of the GCMs to regional scale using RCM is preferred. As a result, Regional Climate Model (RCM) downscaling shows a considerable increase over the globe (Anyah & Semazzi, 2007) for simulations of both historical and projected scenarios. However, simulation of historical and projected scenarios using RCM is very few over EA (Brunswick et al., 2006; Segele et al., 2009; Diro et al., 2012). One of the few experiments is COordinated Regional Climate Downscaling EXperiment (CORDEX) which provides an opportunity of obtaining high-resolution RCM simulations for the historical and projection periods over the EA and the rest of the globe (Evans, 2011; Jacob et al., 2012; Gutowski et al., 2016). The CORDEX simulations of the historical ensemble and future climate as obtained from scenario projection have been used in the study of current and future climate variability and trend. Previous studies for example (Nguvava et al., 201; Mostafa et al., 2019) have employed CORDEX experiments to investigate current and projected climates under various scenarios over EA.
However, there is a significant model bias due to physical representation and parameterization in reproducing the observed rainfall which also varies from model to model and region to region over the globe (Allan & Soden, 2012; Kharin et al., 2007; Church et al., 2013). Various studies have been undertaken to assess the performance of models over EA. For example, (Yang et al., 2014) used five CMIP5 GCMs and identified that the outputs overestimate short rains and underestimate long rains. Models show weak performance in representing observed rainfall in the vicinity of the equator compared with the rest part of the region (Knutti & Sedláˇ, 2012; Woldemeskel et al., 2015. Ongoma et al., 2018; Ongoma et al., 2019) selected eight best performing GCMs from 22 CMIP5 coupled models over EA. (Ayugi et al., 2019) used Rossby Center of Atmospheric Models (RCA4) driven by some CMIP5 GCMS in the simulation of precipitation over the Greater Horn of Africa (GHA) from 1951 to 2005. The authors identified MIROC5, CSIRO, CM5A-MR, Max Plank Institute for Technology (MPI-ESM-LR), and EC-EARTH as best performing GCMs in reproducing observed rainfall over EA. All selected models captured the bimodal and unimodal rainfalls distribution pattern. (Yisehak et al., 2021) studied spatiotemporal characteristics of meteorological drought under changing climate over the northern parts of Ethiopia using five CMIP5 GCMs namely CCSM4, GGFDL-ESM2M, HadGEM2-ES, MIROC5, and MPI-ESM-MR. The author identified that the projected rainfall shows an increasing trend over the northern parts of Ethiopia. Projected drought events are more frequent in RCP4.5 than in RCP8.5 projections (Yisehak et al., 2021). Recently, (Ngoma et al., 2021) used EC-EARTH, IPSL-CMA-MR, MIROC5, CSIRO-MK3.6.0 and MPI-ESM-LR downscaled with RCA4 to assess current and future spatiotemporal precipitation variability and trends over Uganda. The projected rainfall has wet condition between April/May and October whereas reduction in wet condition in March over Uganda (Ngoma et al., 2021). Ongoma et al., 2018) also used CanESM2, CESM1-CAM5, CNRM-CM5, CSIRO-MK3.6-0, and MIROC5 CMIP5 GCM to evaluate change in the mean rainfall projected under RCP4.5 and RCP8.5 over EA. Assamnew & Gizaw (2020) have also assessed CORDEX simulations and they have indicated that MPI, MIROC, CCCma, IPSL, CSIRO, MOHC and MIROC are best performing GCMs when coupled to RCA4 and REMO RCMs for downscaling in most of the rainy season over EA.
Motivated by the on-going improvements in the performance of models over the years, the data have also been used to investigate various aspects of climate change and its impact on water, agriculture, and health sectors. For example, (Gebrechorkos et al., 2018) used regional climate projections for impact assessment in EA. (Lyon & Dewitt, 2012) have shown the long rains over EA exhibits an increasing trend over western Ethiopia and Kenya in contrast to a decreasing trend over Tanzania. (Ongoma et al., 2018) evaluated change in mean rainfall and temperature over EA based on scenario projections from five CMIP5 models. The authors identified increase in rainfall under both RCP4.5 and RCP8.5. Moreover, they have indicated that the increase in rainfall from October to December exceeds that of March to May season. Similarly, (Ongoma et al., 2018) have noted the increase in projected rainfall under RCP8.5 exceeds that of RCP4.5.
The majority of these studies have focused on model validation, drought, decadal variability and long term change in rainfall over identified homogeneous rainfall regimes in EA. For example, (Ongoma et al., 2018) have evaluated changes in mean rainfall and temperature based on CMIP5 models using nine homogeneous sub regions over EA and evaluated decadal variability and trend of annual mean rainfalls. However, there is limited number of studies on precipitation extremes which increased in frequency of occurrence and intensified in recent decades. Precipitation extremes can be quantified by the frequency analysis of rainfall series and precipitation heterogeneity indexes. To our knowledge, we are not aware of any study that employ precipitation heterogeneity index to understand rainfall extremes in EA. In this study, precipitation concentration index (PCI), defined in terms of monthly rainfall series to quantify temporal distribution of precipitation within a given season or year is used to investigate spatiotemporal rainfall characteristics (Oliver, 1980; De Luis et al., 2011). PCI represents the degree to which monthly precipitation is unevenly distributed throughout the season or year which then leads to rainfall extreme conditions such as drought. In addition, standardized precipitation index (SPI) is employed. Therefore, a comparative analysis of these extreme precipitation indices has scientific merit since they provide better understanding of recurrent and potential future drought events in the region. Moreover, the evaluation of rainfall trend and variability at each grid point allows identification of climate change hot spots which, to our knowledge, was not sufficiently covered in any of the previous studies. This is particularly important in view of complex topography that may not be captured by analysis over limited number of homogeneous rainfall regimes. Therefore, this study intends to address these problems and seeks answers to questions such as what is the skill of the mean of climate simulations for the historical period, hereafter referred to as historical Ensemble Mean (HEM), in capturing the observed seasonal precipitation heterogeneity over EA?; How projected rainfalls are changed relative to the HEM and observed (i.e. CRU) rainfalls over EA? In addition, comparison of rainfall trends from projections, observations and HEM during both short and long rains provides insight into a range of possible climate changes, and provides vital and appropriate information for the development of mitigation and adaptation options to climate change. The paper is organized such that Section 2 deals with methodology and Section 3 consists of results and discussion. Finally Section 4 presents conclusions.