2 Abstract Probable changes in mean and extreme precipitation in East Africa are estimated from global climate models (GCMs) prepared for the IPCC Fourth Assessment Report (AR4). Bayesian statistics are used to derive the relative weights assigned to each member in the multimodel ensemble. There is substantial evidence in support of a positive shift of the whole rainfall distribution in East Africa during the wet seasons. The models give indications for an increase in mean precipitation rates and intensity of high rainfall events, but less severe droughts. Upward precipitation trends are projected from early this (twenty-first) century. As in the observations, a statistically significant link between sea-surface temperature (SST) gradients in the tropical Indian Ocean and short-rains (October December) in East Africa is simulated in the GCMs. Furthermore, most models project a differential warming of the Indian Ocean during boreal autumn. This is favourable for an increase in the probability of positive Indian Ocean zonal mode (IOZM) events, which have been associated with anomalously strong short-rains in East Africa. On top of the general increase in rainfall in the tropics due to thermodynamic effects, a change in the structure of the Eastern Hemisphere Walker circulation is consistent with an increase in East Africa precipitation relative to other regions within the same latitudinal belt. A notable feature of this change is a weakening of the climatological subsidence over eastern Kenya. East Africa is shown to be a region in which a coherent projection of future precipitation change can be made, supported by physical arguments. Although the rate of change is still uncertain, almost all results point to a wetter climate with more intense wet seasons and less severe droughts. 1

3 1. Introduction The chemical composition of the atmosphere is changing due to human activities, triggering numerous studies aimed at understanding the sensitivity of the climate system to rising greenhouse gas concentrations (Houghton et al. 2001; Solomon et al. 2007, and references therein). Although the increase in greenhouse gases is relatively uniform around the globe, the response in a myriad of climate variables inevitably exhibits spatial inhomogeneities. For this reason, many studies have focused on possible impacts of climate change on selected climate variables in specific regions of the globe, especially those variables considered economically and socially significant (e.g., Gillett et al. 2004a; Karoly and Braganza 2005; Hegerl et al. 2006; van Ulden and van Oldenborgh 2006; van den Hurk et al. 2006; Beniston et al. 2007). Considerable effort has been devoted to investigating possible changes in mean climate, and increasingly also to changes in variability and extremes, in the Northern Hemisphere continental areas. This is related to the availability of expertise and reliable data in these areas. By comparison, many of the least developed countries, especially in Africa, suffer from both a lack of high-quality data and lack of research attention and capacity. Hence, there are far fewer climate analyses for these regions, especially in relation to variability and extremes. Impact analyses show that such variability and extremes disproportionally affect the poorest countries and the poorest people (e.g. African Development Bank (AfDB) et al. 2003). Trying to address that gap, Shongwe et al. (2009)(hereafter SHO09) presents an analysis of projected precipitation changes for southern Africa. In this paper, a similar approach is applied to East Africa, defined as the area lying between E, 10 S 4 N. The region experiences a semi-annual rainfall cycle with two major rainfall peaks in boreal spring [March May (MAM); also known as long-rains] and autumn [September December (SOND); 2

4 also known as the short-rains]. During the latter season (short-rains), atmospheric dynamics in the first month (September) differ from the rest of the season (October December). Also, teleconnections with large-scale features such as the El-Niño/Southern-Oscillation (ENSO) are different between September and OND (Mutai and Ward 2000). It is common practice therefore to omit September when referring to the short-rains in East Africa (e.g. Clark et al. 2003; Anyah and Semazzi 2007). We adopt the same convention in this paper. Some portion of the selected area has one rainfall maximum during boreal winter spring months (November April). Traditionally, climate change studies in which Africa features have tended to focus on likely precipitation changes during boreal summer and winter (Solomon et al. 2007). Much less is known about the transition seasons (i.e. boreal autumn and spring), which are the rainy seasons in East Africa. Furthermore, there has been a tendency to focus at large regions within Africa or in some cases on the Giorgi regions (Giorgi and Francisco 2000). While this approach is useful in diagnosing the climate change signal on rainfall as a first step, it tends to overlook the effect of local features such as East Africa s varied topography (e.g. mountain ranges and rift valleys) and large water bodies (e.g. Lake Victoria) on the large scale signal. Although most of these local controls are not yet adequately represented in low-resolution climate models such as those used here, the homogenous subregions identified by Indeje et al. (2000) are used to show that rainfall response to global warming is not uniform across the region. In recent years, East Africa has suffered frequent episodes of both excessive (e.g. Webster et al. 1999; Latif et al. 1999) and deficient rainfall (e.g. Hastenrath et al. 2007). In particular, the frequency of anomalously strong rainfall causing floods has increased. Our analysis of data from the International Disaster Database EM-DAT (http://www.em-dat.net/), reveals that there has also been an increase in the number of reported hydrometeorological disasters in the region, from an average of less than 3 events/year in the 1980s, to over 7 events/year 3

5 in the 1990s, and almost 10 events/year from 2000 to 2006, with a particular increase in floods: from an average of less than 1 event/year in the 1980s to 7 events/year between 2000 and In the period , these hydrometeorological disasters affected on average almost two million people per year. In addition, they have severe impacts on economic performance and poverty alleviation (e.g., Hellmuth et al. 2007). Furthermore, we note that many of the region s epidemics, which dominate the remainder of the reported disasters in the EM-DAT database, are also affected by climatic conditions. While part of the trend is due to better reporting, it also reflects a rising vulnerability to natural hazards, and potentially an underlying trend in climate variability and extremes. These rising risks are receiving increasing attention from policy makers, for example, in the Decision and Declaration on Climate Change and Development adopted by the African Union Head of State Summit in January 2007, and in new efforts by East African governments to better manage climate risks. In Kenya for instance, the Kenya Adaptation to Climate Change in Arid Lands (KACCAL) project aims to address the increasing risk of both floods and droughts (World Bank 2006). In this context, there is an obvious need for better analyses of the likely response of extreme climate events in this region to global warming, to inform disaster preparedness and development planning. Such demands provided the motivation for the present study. The primary aim is to assess how the intensity of seasonal precipitation extremes is likely to change in the region, against the backdrop of probable changes in mean precipitation. A secondary aim is to assess the spatial inhomogeneities in the model projections. We also present physical mechanisms that may explain the simulated changes in the precipitation probability distribution. 4

6 2. Data and Methods a. Model simulations and observations A brief description of the datasets used and the analysis methods applied are given below. Further details can be found in SHO09. The models used in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4; Solomon et al. 2007) form the major input to the analysis of changes in precipitation patterns. The output has been made available as part of the World Climate Research Program (WCRP) Coupled Model Intercomparison Project Phase 3 (CMIP3). We have chosen model projections driven by the intermediate SRES A1B scenario, a standard emission scenario currently roughly corresponding to observed CO 2 concentrations. In this scenario, there is no drastic reduction of CO 2 emissions up to 2100, when concentrations reach twice their pre-industrial levels. After 2100, experiments have been continued with constant CO 2 concentrations at this level. A subset of 12 coupled general circulation model (CGCM) simulations has been selected from the data available at the time of first submission (January 2008). The selection was based on a pattern comparison of the seasonal cycle in the simulations for the twentieth century (20c3m) with the Climate Research Unit (CRU TS2.1; New et al. 2000) gridded precipitation. The 12 selected CGCMs were among the six best-performing models (as measured by the mean of the monthly correlation coefficient and root-mean-squared difference with the analysis) over at least one of four domains of sub-saharan Africa: southern Africa, East Africa, North-East Africa and West Africa. We also demanded that enough data at high CO 2 concentrations was available to be able to compute return levels. Further details of the selection procedure can be found on the web site scenarios under technical details. The selected CGCMs, their spatial resolutions, ensemble sizes and references are given in 5

7 Table 1 of SHO09. The model simulations are linearly interpolated to a common lat/lon grid, corresponding roughly to T95 resolution. As in SHO09, two subsamples have been selected from the model simulations. The and periods define the present (20c3m) and future climate, respectively. In addition to monthly CMIP3 precipitation, the horizontal wind vector and pressure vertical velocity (ω) fields are used. These variables are used to estimate the strength and structure of the Eastern Hemisphere zonal (Walker) circulation and their projected changes. Three of the models used in the precipitation change analysis have been omitted in the analysis of the zonal circulation. The CGCM background Walker circulation is compared with that obtained from the European Centre for Medium-Range Weather Forecasts (ECMWF) 40- year Reanalysis (ERA40; Uppala et al. 2005). Observed twentieth century precipitation used in this study was taken from the CRU TS2.1 gridded station data. Unlike in SHO09, dearth of station observations in the Global Historical Climatology Network (GHCN) data precluded delineating homogenous rainfall zones. For this reason, the climate zones defined by Indeje et al. (2000) have been adopted in this study. These zones are shown in Fig. 1. The western zones of Indeje et al. (2000) are combined into one zone (Zone III; Fig. 1) to allow a larger spatial sampling. In most models used here, precipitation in these areas have a similar probability density functions (PDFs; not shown). The CRU TS2.1 grid points falling within each zone are averaged. All the zones except much of Tanzania (Zone IV) have a semi annual cycle. In Zone IV, the annual cycle of area-averaged precipitation showed a peak during austral summer autumn months (November April). The CMIP3 simulated precipitation is spatially averaged analogous to the CRU TS2.1 data. The coastal region, the central region of Kenya and the Lake Victoria The ECHO-G model has been omitted because horizontal wind and pressure velocity fields are not available on the IPCC-AR4 archive in both their 20c3m and SRESA1b runs. For the CSIRO Mk3.0, the vertical velocity fields are not available. For the UKMO HadCM3, SRESA1b horizontal wind and vertical velocity fields are not available. 6

8 area are small to be adequately resolved by some low-resolution CGCMs. These regions have therefore been omitted in the analysis. Monthly precipitation from the CRU TS2.1 gridded station and CMIP3 data are accumulated into seasonal totals for each rainfall season of East Africa. Mean precipitation rates (mm day 1 ) in a given season are calculated by dividing the seasonal totals by the number of days within that season. Climatologically wettest seasons are defined as those with the highest mean precipitation rates. On the other hand, driest seasons, characterising meteorological droughts, have the lowest mean precipitation rates. b. Extreme value analysis The peak over threshold or Generalised Pareto Distribution (GPD) is used in this study to represent the distribution of observed and simulated extreme seasonal precipitation rates. A description of GPD and its applications is found in Coles (2001). The quality of the GPD fit is then assessed using Anderson-Darling goodness-of-fit tests (Laio 2004). Details of how the GPD has been applied in our work including statistical tests is discussed by SHO09. Return levels are estimated from the fitted GPD. Return levels are frequently used in extreme precipitation studies in climate research (e.g., Kharin and Zwiers 2000; Meehl et al. 2005; Kharin et al. 2007) because of the simplicity of their interpretation. The return level (z p ) is the threshold likely to be exceeded in a given year with probability p, or the level likely to be exceeded once in every 1/p years. In this study 10-, 20-, 50- and 100-year return levels have been estimated. The return values are computed for each homogenous zone (Fig. 1) in both the OND and MAM seasons. Owing to the small sample size of threshold excesses, only results based on the least biased 10-year return levels are presented in this paper. These GPD quantiles (10-year return levels) are obtained from an interpolation rather than an extrapolation. Results obtained using 100-year return levels are available on the KNMI 7

9 Africa scenarios web-site (http://www.knmi.nl/africa scenarios/). When modelling dry extremes following the approach discussed by SHO09, a problem is encountered during extrapolation to much longer return periods because these extremes are by definition constrained by zero above (i.e. dry extremes have a finite upper bound). In practice almost all GPD fits used in this study have a negative shape parameter ξ and the resulting finite return level for infinite return times are in fact positive. The finite upper bound of dry extremes is therefore not a concern in this study. c. Multi-model ensembling Uncertainties are inherent in long-term climate model simulations. These can be attributed to natural climate variability, different model responses to a given forcing (such as increases in greenhouse gas concentration) and those associated with the emission scenarios used to force the climate models. In climate change studies devoted to precipitation extremes and their possible future changes, it has become common practice to use simulations from a range of climate models (e.g., Kharin et al. 2007). Research on how best to combine simulations from several climate models through objective weighting is ongoing (e.g. Giorgi and Mearns 2002, 2003; Tebaldi et al. 2004, 2005). Using Bayes theorem and making certain assumptions, it is possible to objectively assign weights to different climate models leading to a probability distribution of future climate change. These assumptions detailed by Tebaldi et al. (2005) relate to model bias, model independence and the similarity of physical mechanisms determining the unforced and future climate. This method has been adopted in the present study. We note here that equal weighting of the models yields similar mean results to the Bayesian weighing method used to obtain the results presented in Section 3. However, the credibility intervals obtained by the Bayesian method are smaller due to, amongst other things, down-weighting obvious outliers. 8

10 It should be noted that there is the possibility that this makes the method overconfident. For comprehensive details of this method, the reader is referred to Tebaldi et al. (2005). The Bayesian method requires that the input variables have Gaussian likelihoods. Mean precipitation taken over sufficiently long periods (100 or 150 years in this study) meet this requirement from the central limit theorem. To fulfil this requirement in the case of extreme quantiles, the first step involved estimating the GPD parameters from maximum likelihood. After checking for model quality and inverting the GPD to obtain the 10-year return levels, we then use the delta method described by Coles (2001) to estimate the uncertainty on the return levels. We verified that the 95% credibility interval obtained using this method is very similar to the interval computed with a non-parametric bootstrap method. More details of the methodology used to obtain the GPD quantiles can be found in Appendix A of SHO09. To summarise, in the Bayesian framework used, the ratio of the conditional posterior mean of future climate (ν) to control climate (µ) simulations is used to define the relative precipitation change ( P ) as P = 100( ν 1). (1) µ For mean precipitation rates and 10-year wettest events, P > 0 is indicative of an increase in their intensity. An increase in the severity of 10-year driest events is indicated by P < 0. As previously mentioned, these changes compare statistics computed using data for the (present) period with their counterparts in the (future) period. 3. Changes in mean and extreme precipitation The Bayesian method adopted in this study enables an estimation of the posterior distribution of the parameters of interest (Tebaldi et al. 2005). The uncertainties associated with 9

11 these parameters is then inferred from their posteriors. It should be kept in mind that these depend on the prior and the availability of GCM data, and as such only give an indication of the uncertainties in the return levels. In the following sub-sections, we present the 95% credibility interval of P as obtained from the Bayesian method for each zone and rainy season of East Africa. We present the results on maps to allow a concise visual interpretation and an easy assessment of the spatial pattern of the projected changes. We also give the (symmetrical) 95% credibility interval. When this interval excludes zero (i.e., no change) the percentage is shaded gray. When zero is outside the 99% credibility interval of P this is denoted by two asterisks. a. Short-rains season (October December; OND) Mean OND precipitation increases are simulated almost everywhere in East Africa (Fig. 2a). Over the semi-arid areas in northern Kenya (Zones I and II; Fig. 1), the western parts, which include Rwanda and Burundi (Zone III), mean precipitation is projected to increase by more than 10%. These increases achieve exclude zero from the 95% credibility interval almost everywhere, and zero is excluded from the 99% credibility interval to the north. Over much of Tanzania (Zone IV), the model projections provide evidence in support of an increase in precipitation rates during austral summer autumn months (November April). While most models show an increase in mean OND precipitation, the HadGEM1 model shows decreases in rainfall almost everywhere in the area. The convergence criterion used in the Bayesian weighting (Tebaldi et al. 2005) treats this model as an outlier and down-weighs it. A reduction of the severity of 10-year driest seasons are found over almost the entire East African region (Fig. 2b). The magnitudes of the simulated reduction in the severity of these OND dry extremes are comparable with those of the mean precipitation rates. Larger reductions ( 32%), are found in northern Kenya and Uganda (Zone II). Elsewhere in the 10

12 north, the projected lessening of 10-year droughts severity is not statistically significant. Austral summer autumn 10-year droughts are becoming less severe in much of Tanzania by about 14%. In common with changes in the mean precipitation rates, widespread increases in the intensity of 10-year wettest OND seasons are simulated (Fig. 2c). Significant increases exceeding 10% are found in Kenya, Uganda and the western parts. In northeastern Democratic Republic of Congo (DRC) and much of Tanzania (Zone IV), increases in the 10-year wettest austral summer seasons average 15%. b. Long-rains season (March May) Changes in mean precipitation during the long-rains season (MAM) have generally the same sign and magnitudes as those simulated for the short-rains (Fig. 3a). The Bayesian weighted model simulations project more than 15% increase in mean precipitation rates over much of East Africa. The upper credibility limit exceeds 25% over a large area to the north and west. However, in contrast to the short-rains, the models perform poorly in simulating the twentieth century climate during this season. During the model pre-selection process, relative to the OND season, correlation between the monthly CRU and model simulated MAM precipitation was found to be low while the RMSE scores were found to be higher (not shown). The difficulty in modelling precipitation during the long-rains has been found in seasonal climate predictability studies, and is due to a dominance of internal atmospheric variations (not provably connected to other components of the climate system; Mutai et al. 1998). The uncertainty associated with these patterns is higher than in the short-rains season. The models show a reduction in the severity of 10-year droughts almost everywhere except over eastern Kenya (Fig. 3b). In eastern Kenya, where changes in dry and wet extremes have opposite signs, an increase in interannual rainfall variance is implied. Changes 11

13 in 10-year wettest seasons are of the same sign as those of mean precipitation rates implying that floods are likely to become more intense (Fig. 3c). However, we place less emphasis on these projections owing to their lower reliability. Again, the HadGEM1 simulations that simulate a drying trend everywhere has been down-weighted by the convergence criterion. c. Time series analysis Time series plots for East Africa precipitation during the twentieth (observed and simulated) and twenty-first centuries (simulated) and in each rainy season are shown in Fig. 4. Unlike in Figs. 2 and 3, where the CGCMs are subjected to a Bayesian weighting, the models are equally weighted in these plots. Despite the wider inter-model dispersion in their twenty-first century projections [caused mainly by the outlying HadGEM1 (dry) and MIROC3.2 (wet) CGCMs], there is substantial evidence in support of an increase in the amount of rainfall. The sharp precipitation increase in East Africa emerges from the early part of the twentyfirst century and in most cases is a reversal of drier conditions experienced during much of last century. We note however that in some regions (e.g., Burundi, Rwanda and Tanzania; Figs. 4c,d), the observed trend over the last decades of the 20th century is opposite the long-term CGCMs simulated trend. This can be due to decadal variability that happens to be opposed to the long-term trend in the models, or due to CGCMs deficiencies in the small regions considered in our study. Further modelling work to distinguish between these possibilities is required. A trend toward predominantly positive precipitation anomalies is also present in the longrains twenty-first century time series (Fig. 4e f). However, the signal-to-noise ratio in these projections is low. Noteworthy, the results from the Bayesian weighted and unweighted simulations are broadly similar. By down-weighing outlying CGCMs, the Bayesian procedure offered the additional 12

14 value of reducing the model dispersion by about a factor two, particularly in the twentysecond century simulations when the model spread is largest. 4. Projected changes in large-scale forcing The robust climate change signal on East Africa precipitation presented above prompts us to endeavour to diagnose likely physical mechanisms and forcings. There is considerable consensus between observational, theoretical and modelling studies with regard to an increase in vertically-integrated atmospheric water vapor as the climate warms (Zveryaev and Chu 2003; Trenberth et al. 2005; Zveryaev and Allan 2005). Related to this is a robust projection of precipitation increase in the deep tropics, which has been detected in the current climate (Zhang et al. 2007). The zonal mean precipitation in the tropics increases because the 7%/ C increase in precipitable water, which follows the Clausius-Clapeyron relationship, is not completely compensated by a slow-down of the tropical circulation (Held and Soden 2006; Vecchi and Soden 2007). The increase in precipitation in the tropics is not zonally uniform, with East Africa precipitation projected to increase more than the zonal mean in the multi-model mean. These zonal asymmetries in tropical precipitation response to global warming suggest that other mechanisms also exert a significant influence. These partly stem from the atmospheric dynamic response, which is not zonally symmetric and leads to a horizontal redistribution of water vapor. To understand this we turn to the well understood interannual variability in the region. Our hypothesis is that the climate change signal will be communicated to the rainfall field partly through the same factors which control interannual rainfall variability in the present climate. At interannual timescales, Hastenrath et al. (1993) discussed atmospheric forcing 13

15 on the Indian Ocean hydrosphere. A number of papers documenting atmosphere-ocean coupling have appeared in literature (e.g. Webster et al. 1999; Saji et al. 1999; Goddard and Graham 1999). It has been found that atmospheric anomalies of pressure and wind typically precede sea-surface temperature (SST) anomalies (Hastenrath and Polzin 2005). Subsequent SST anomalies, in part forced by wind-stress anomalies, feed back to the atmosphere, further reinforcing atmospheric anomalies (Goddard and Graham 1999), which often result in precipitation anomalies. Against this background, modelled changes in SST and atmospheric circulation during the short-rains season are examined in this section. a. Projected OND Indian Ocean SST pattern and rainfall relationship Observational and modelling studies have linked East African high-frequency rainfall variability during boreal autumn to SST anomalies (Goddard and Graham 1999; Black et al. 2003). Earlier studies aimed at understanding the link between East Africa precipitation and equatorial Pacific Ocean SSTs associated with the El Niño-Southern Oscillation (ENSO; e.g., Hastenrath et al. 1993; Indeje et al. 2000). It was found that ENSO exerts some influence on East Africa short-rains such that rainfall is enhanced (suppressed) during warm (cold) events in the eastern equatorial Pacific Ocean. It has been shown that the Indian Ocean is one of the pathways through which the ENSO signal propagates into East Africa precipitation. Often, the west Indian Ocean warming lags that of the eastern equatorial Pacific by a few months (Klein et al. 1999). However, evidence has been presented supporting the existence of an Indian Ocean mode of SST variability independent of ENSO (Saji et al. 1999; Webster et al. 1999). In fact OND 1961, one of the wettest in the twentieth century, was not related to ENSO (figure not shown). The models used in this study do not give a robust signal with respect to changes in the mean state of ENSO (not shown). Based on these findings and the fact that changes in ENSO properties 14

16 are unlikely to exceed natural variability (van Oldenborgh et al. 2005), possible changes in Indian Ocean SST patterns are investigated in this paper. This is motivated by results from studies which found Indian Ocean SST anomalies to be the dominant factor controlling East Africa short-rains (Latif et al. 1999; Goddard and Graham 1999; Black et al. 2003; Clark et al. 2003; Behera et al. 2005; Washington and Preston 2006). In this study, the CGCM simulated East Africa short- and long-rains were subjected to a standard empirical orthogonal function (EOF) analysis (see e.g. van den Dool 2007). The first EOFs are a monopole patterns describing the largest source of interannual variability in this region, particularly during the short-rains. The corresponding principal component series are then regressed on simultaneous Indian Ocean SSTs from each model. The regression of the first principal components (PC 1) on Indian Ocean SSTs are displayed in the first and last columns of Fig. 5. Areas with significant regression at the 5% level (from a two-sided local t-test ) are shaded. A statistically significant positive link is found with SSTs in the western tropical Indian Ocean (WTIO), west of 80 E in most models. In some models (e.g. CSIRO-Mk3.0, ECHAM5/MPI-OM, GFDL-CM2.0, GFDL-CM2.1 and ECHO-G) the ocean area with significant positive association with precipitation extends into the Bay of Bengal. Negative regression coefficients are found in the southeastern Indian Ocean off the coast of Sumatra in almost all models. The figures demonstrate that wetter conditions in East Africa in the CGCMs during boreal autumn often occur when SSTs are warmer in the western Indian Ocean and cooler close to the Maritime Continent. These results are broadly similar to previous findings on Indian Ocean East Africa rainfall relationships in observational and modelling research work (Black et al. 2003; Clark et al. 2003). The major atmospheric response to the anomalous zonal SST gradients relevant for East Africa short-rains is a perturbed Indian Ocean Walker cell (Behera et al. 2005). Low-level easterly anomalies south and close to the equator are a prominent feature in this perturbed local Walker circulation. 15

17 These moisture-laden winds feed the diabatic heating-induced anomalous convection close to the warm SST anomaly, with the obvious consequence of enhanced precipitation. There is no similarly significant link between the long-rains (MAM) and simultaneous Indian Ocean SSTs in the models (last column of Fig. 5). Even in models which show some association (e.g. CNRM-CM3, MIROC3.2, ECHAM5/MPI OM, ECHO G and GFDL-CM2.0), the remarkable rainfall link to tropical Indian Ocean SST gradients found during the OND season is not present. We do not know of any documented external control (from other components of the climate system) for the long-rains interannual variability. For this reason, our discussion will focus on SST patterns likely to influence OND rainfall. Having established the existence of the Indian Ocean SST signal on rainfall, the model simulated SST differences ( minus ) are examined. Boreal autumn (OND) Indian Ocean SST differences are displayed in the middle column of Fig. 5. Statistically significant (at the 1% level) basin-wide warming is found. Although SSTs warm throughout the tropical Indian Ocean basin, zonal asymmetries in the warming are evident in most models. The western part of the basin is generally projected to warm more than the eastern. These results are in qualitative agreement with Vecchi and Soden (2007) who show a tendency towards upwelling (downwelling) along the eastern (western) equatorial Indian Ocean in the twenty-first century (their Fig. 15). A similar differential warming is observed in interannual variability. It stems from windevaporation-sst and wind-thermocline-sst feedbacks which cool the ocean surface in the south-east Indian Ocean (SEIO) from boreal summer to fall before giving way to basin wide warming from around December (Shinoda et al. 2004). Vecchi and Soden (2007) argue that atmospheric anomalies drive changes in tropical ocean dynamics, and hence the ocean thermal structure in the twenty-first century. The simulated differential warming and associated anomalous zonal SST gradients have implications for the east-west Indian SST dipole that occasionally peaks during boreal fall months (Behera et al. 2005). As demonstrated above, 16

18 enhanced rainfall in East Africa is often associated with anomalously warm (cold) SSTs in the WTIO (SEIO). We proceed by defining an index WTI expressed as average SST in the rectangular domain to the west (50 70 E, 10 N 10 S) and another index SEI for the box in the eastern extremity ( E, 0 15 S). The Indian Ocean zonal mode index (IOZMI), which represents zonal SST gradients across the near-equatorial Indian Ocean is then defined as IOZMI = WTI SEI. These boxes are indicated in the first two columns of Figure 5. Probability density functions (PDFs) of IOZMI are shown in Fig. 6. Almost all the models displayed in the figure show a clear shift in the PDFs towards a higher probability of positive IOZM phases. This enhanced probability translates to a higher probability of excessive shortrains in East Africa in the future climate. These findings are consistent with the results presented in Section 3a. On average, IOZM-shift contributes to 30% (95% confidence interval is [12,49]) of the overall precipitation change (results not shown). b. Changes in Eastern Hemisphere Walker circulation In this section, we investigate the circulation changes that give rise to the correlation between the OND Indian Ocean zonal SST gradient and simultaneous precipitation in East Africa. The main changes are in the Walker circulation, which are shown for much of the eastern Hemisphere stretching from 20 W 160 E. To define the Walker circulation, we employ the approach used by Chen (2005), and many others, to define atmospheric overturning circulations. In this approach, the first step involves decomposing the horizontal wind vector into its divergent and rotational components. A zonal (meridional) cell is then defined from the zonal (meridional) component of the divergent wind and the vertical velocity. The tropical zonal circulation plotted in Fig. 7a show three overturning cells: 1) A narrow cell featuring strong ascending motion to the west and subsidence to the east is found in 17

19 East Africa extending towards the western Indian Ocean. The descending branch of this cell coincides with the semi-arid areas in East Africa. The ascending branch correspond to the major source of diabatic heating over Congo. We refer to this as the East African Walker cell. 2) A shallower East Atlantic cell is found on the west flank of the East African Walker cell. 3) To the east, the East African cell is flanked by a broader Indian Ocean cell. In most models used here, the background zonal circulation is broadly similar to that found in the ERA40 reanalysis despite the differences in the analysis period considered ( and for the PCMDI and ERA40 data, respectively). A dominant feature in projected changes in atmospheric upward vertical velocity is a weakening of the ascending branch of the East Atlantic Walker cell over central Africa. A reduction in ω exceeding mbs 1 is found in most models in the middle to upper troposphere. Almost all the CGCMs, show positive differences in the omega velocities over the descending branch of the narrow East African Walker cell. Although not statistically significant in some of the models, these changes imply a weakening of the climatological subsidence over the eastern semi-arid regions. These results are in qualitative agreement with Vecchi and Soden (2007), who used different metrics to measure the intensity of the zonal circulation than are used here. The future behaviour of the eastern Hemisphere Walker circulation show a tendency towards positive IOZM-like state, consistent with the results presented in Sec. 4a suggesting that the IPCC-AR4 models have physical parameterizations of convective processes that respond realistically to the changes in forcings found here. Vecchi and Soden (2007) argue that atmospheric thermodynamics are the principal drivers behind the weakening of the Indo- Pacific Walker circulation. Investigating the cause of the weakening is not the objective of the present study. The aim of our analysis has been to diagnose physical factors underlying the modelled precipitation changes. 18

20 5. Discussion and Conclusions In this study, we investigated possible changes in the intensity of mean and extreme precipitation rates in East Africa from the CMIP3 multi-model dataset. Evidence in support of a future positive shift in the rainfall distribution under global warming has been presented for most models, with only a single model simulating a trend to less rain. Increases in both mean precipitation rates and the intensity of 10-year wettest events are simulated almost throughout the region, while dry extremes are becoming less severe. This is the case even in the semi-arid climate in northern Kenya and southern Somalia. The presence of a dipole between southern Africa (cf. SHO09) and East Africa in the precipitation response, a robust feature of interannual rainfall variability in the climate system (e.g. Goddard and Graham 1999), is particularly reassuring. A prime example of this feature is the austral summer season. Unprecedented wet conditions in East Africa (Kenya) in November preceded one of the worst droughts in much of southern Africa. The qualitative agreement with previous findings (Kharin and Zwiers 2000; Giorgi and Mearns 2002, 2003; Tebaldi et al. 2004) provides an additional reassurance. A robust feature across the CGCMs is an increase in atmospheric column integrated water vapor during the transition seasons (autumn and spring), particularly in the tropical region. Despite discrepancies in the magnitudes of the linear trends, all models show upward trends in precipitable water in the future climate (not shown). Although the increase in water vapor in a warmer atmosphere is offset to some extent by a slowdown of the tropical circulation (Vecchi and Soden 2007), this trend alone is favourable for an increase in the precipitation almost everywhere in the tropics by about 3% (Held and Soden 2006). Dynamical effects however force spatial inhomogeneities in the global warming-induced tropical rainfall increases. For instance, projections for rainfall in East Africa are higher than the zonal mean. This is in part caused by anomalous moisture flux convergence over East Africa (discussed 19

IGAD CLIMATE PREDICTION AND APPLICATIONS CENTRE CLIMATE WATCH FOR THE PERIOD JUNE -DECEMBER 2015 EL NIÑO EVOLUTION OVER EASTERN EQUATORIAL OCEAN REGION AND POTENTIAL IMPACTS OVER THE GREATER HORN OF FRICA

Queensland rainfall past, present and future Historically, Queensland has had a variable climate, and recent weather has reminded us of that fact. After experiencing the longest drought in recorded history,

1 of 10 3/09/2015 1:39 PM ENSO Wrap-Up Current state of the Pacific and Indian Ocean El Niño strengthens but a warm Indian Ocean Issued on 1 September 2015 Product Code IDCKGEWW00 The 2015 El Niño is now

SECTION 3 Making Sense of the New Climate Change Scenarios The speed with which the climate will change and the total amount of change projected depend on the amount of greenhouse gas emissions and the

Seasonal Climate Watch August to December 2016 Date: Jul 20, 2016 1. Advisory Current observations show a neutral ENSO (El-Niño Southern Oscillation) state. Despite the fact that during the previous month,

Climate Extremes Research: Recent Findings and New Direc8ons Kenneth Kunkel NOAA Cooperative Institute for Climate and Satellites North Carolina State University and National Climatic Data Center h#p://assessment.globalchange.gov

3072 JOURNAL OF CLIMATE Ocean Roles in the TBO Transitions of the Indian Australian Monsoon System JIN-YI YU ANDSHU-PING WENG Department of Earth System Science, University of California, Irvine, Irvine,

Ocean Circulation and Climate In addition to the atmospheric circulation heat is also transported to the poles by the ocean circulation. The ocean circulation is therefore an important part of the climate

Kathryn Sullivan, Ph.D, Acting Under Secretary of Commerce for Oceans and Atmosphere and NOAA Administrator Thomas R. Karl, L.H.D., Director,, and Chair of the Subcommittee on Global Change Research Jessica

THE CURIOUS CASE OF THE PLIOCENE CLIMATE Chris Brierley, Alexey Fedorov and Zhonghui Lui Outline Introduce the warm early Pliocene Recent Discoveries in the Tropics Reconstructing the early Pliocene SSTs

VOL. 16, NO. 11 JOURNAL OF CLIMATE 1JUNE 2003 The Asian Monsoon, the Tropospheric Biennial Oscillation, and the Indian Ocean Zonal Mode in the NCAR CSM* JOHANNES LOSCHNIGG International Pacific Research

Climate Change on the Prairie: A Basic Guide to Climate Change in the High Plains Region - UPDATE Global Climate Change Why does the climate change? The Earth s climate has changed throughout history and

What are the Experts Saying About Effects of Climate Change on Rainfall and Streamflow in the Southeast? by David H. Moreau Director, Water Resources Research Institute of the University of North Carolina

The Oceans Role in Seasonal and Longer Term Climate Why the recent cooling is likely just the start Although, I believe ultimately the sun is the primary driver for the changes to global climate, the oceans

Improving Hydrological Predictions Catherine Senior MOSAC, November 10th, 2011 How well do we simulate the water cycle? GPCP 10 years of Day 1 forecast Equatorial Variability on Synoptic scales (2-6 days)

The IPCC Special Report on Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation A changing climate leads to changes in extreme weather and climate events 2 How do changes

Summary for Policymakers: The Science of Climate Change - IPCC Working Group I Contents 1. Greenhouse gas concentrations have continued to increase 2. Anthropogenic aerosols tend to produce negative radiative

376 THE SIMULATION OF TROPICAL CONVECTIVE SYSTEMS William M. Frank and Charles Cohen Department of Meteorology The Pennsylvania State University University Park, PA, 16801 -U.S.A. ABSTRACT IN NUMERICAL

2138 JOURNAL OF CLIMATE Coupled Ocean Atmosphere Dynamical Processes in the Tropical Indian and Pacific Oceans and the TBO GERALD A. MEEHL, JULIE M. ARBLASTER, AND JOHANNES LOSCHNIGG* National Center for

Appendix. Global Circulation and El Niño The radiant energy from the Sun is not equally distributed over the surface of the Earth. Equatorial latitudes receive much more energy than polar latitudes, and

A changing climate leads to changes in extreme weather and climate events the focus of Chapter 3 Presented by: David R. Easterling Chapter 3:Changes in Climate Extremes & their Impacts on the Natural Physical

Reply to No evidence for iris Richard S. Lindzen +, Ming-Dah Chou *, and Arthur Y. Hou * March 2002 To appear in Bulletin of the American Meteorological Society +Department of Earth, Atmospheric, and Planetary

2: The Causes of Climatic Change Mick Kelly The character of the Earth's climate system is shaped by the general circulation of the atmosphere and oceans, the winds and currents that move heat and moisture

SPATIAL DISTRIBUTION OF NORTHERN HEMISPHERE WINTER TEMPERATURES OVER THE SOLAR CYCLE DURING THE LAST 130 YEARS Kalevi Mursula, Ville Maliniemi, Timo Asikainen ReSoLVE Centre of Excellence Department of

El Niño: How it works, how it affects the South Pacific William Kessler NOAA / Pacific Marine Environmental Laboratory Seattle, USA Outline of talk The normal situation in the tropical Pacific: a coupled

Will we ever be able to attribute individual weather events to anthropogenic climate change? Santiago de Compostela, July, 2003 Myles Allen Departments of Physics, University of Oxford myles.allen@physics.ox.ac.uk

The Weather of Pemba 3 years of SUN hours, winds & precipitation A unique weather system Pemba Island has a tropical monsoon climate with high precipitation figures all year round. Pemba Island has two

Diurnal and Semi-diurnal Variations of Rainfall in Southeast China Judy Huang and Johnny Chan Guy Carpenter Asia-Pacific Climate Impact Centre School of Energy and Environment City University of Hong Kong

LESSON 3: AFRICA S CLIMATE REGIONS Key Concepts You must know, or be able to do the following: Name, understand the characteristics and position of Africa s major climate regions Be able to link the African

Supplementary Material to A simple scaling approach to produce climate scenarios of local precipitation extremes for the Netherlands G. Lenderink and J. Attema Extreme precipitation during 26/27 th August

4 June 2015 GIEWS Update El Niño Southern Oscillation (ENSO) Possible impact on agricultural production during the second half of 2015 Highlights: In March 2015, the arrival of El Niño was officially declared

3) The movement of air in convective currents, such as Hadley cells, does not occur in straight lines because of the rotation of the earth. This is known as a) The Rotation Effect b) The Coriolis Effect

Near Real Time Blended Surface Winds I. Summary To enhance the spatial and temporal resolutions of surface wind, the remotely sensed retrievals are blended to the operational ECMWF wind analyses over the

1 Rwanda Rwanda is a landlocked country which lies within latitudes 1-3 S and longitudes 28-31 E and bordered by Uganda in the north and Tanzania in east while in the south and west are Burundi and the

Australian Temperature Variations - An Alternative View John McLean 11 October 2007 Appendix by Dr. Thomas Quirk added 17 October 2007 Australia's temperature since 1950 (see Figure 1) is usually described

Indo-French Workshop on Atmospheric Sciences 3-5 October 2013, New Delhi (Organised by MoES and CEFIPRA) Indian Ocean and Monsoon Satheesh C. Shenoi Indian National Center for Ocean Information Services

Winter 2014-2015: Illustrating the Strength of WCS Methods Amplified and persistent climate patterns in the winter of 2014-2015 brought unusual conditions to many parts of the Northern Hemisphere. Highly

IEAGHG Information Paper 2015-10; The Earth s Getting Hotter and So Does the Scientific Debate A recent study published in Nature Climate Change 1 suggests that the rate of climate change we're experiencing

Chapter Overview CHAPTER 6 Air-Sea Interaction The atmosphere and the ocean are one independent system. Earth has seasons because of the tilt on its axis. There are three major wind belts in each hemisphere.

1968-32 Conference on Teleconnections in the Atmosphere and Oceans 17-20 November 2008 Impact of the Indian and the North Pacific ocean on ENSO variability in a hybrid coupled model, YEH Sang-Wook Korea