Future Precipitation Projections over Central and Southern Africa and the Adjacent Indian Ocean: What Causes the Changes and the Uncertainty?

Futureprojectionsofprecipitationatregionalscalesarevitaltoinformclimatechangeadaptationactivities. Therefore, is it important to quantify projected changes and associated uncertainty, and understand model processes responsible. This paper addresses these challenges for southern Africa and the adjacent Indian Ocean focusing on the local wet season. Precipitation projections for the end of the twenty-ﬁrst century indicate a pronounced dipole pattern in the CMIP5 multimodel mean. The dipole indicates future wetting (drying) to the north (south) of the climatological axis of maximum rainfall, implying a northward shift of the ITCZ and south Indian Ocean convergence zone that is not consistent with a simple ‘‘wet get wetter’’ pattern. This pattern is most pronounced in early austral summer, suggesting a later and shorter wet season over much of southern Africa. Using a decomposition method we determine physical mechanisms underlying this dipole pattern of projected change, and the associated intermodel uncertainty. The projected dipole pattern is largely associated with the dynamical component of change indicative of shifts in the location of convection. Over the Indian Ocean, this apparent northward shift in the ITCZ may reﬂect the response to changes in the north–south SST gradient over the Indian Ocean, consistent with a ‘‘warmest get wetter’’ mechanism. Over land subtropical drying is relatively robust, particularly in the early wet season. This has contributions from dynamicalshifts in the locationof convection, which may be related to regional SST structures in the southern Indian Ocean, and the thermodynamic decline in relative humidity. Implications for understanding and potentially constraining uncertainty in projections are discussed.


Introduction
Africa is highly vulnerable to climate change, evident from the reliance on seasonal precipitation for agriculture, water supply, and energy generation over the majority of sub-Saharan Africa (Basher and Briceño 2006;Meadows 2006). The region exhibits relatively low adaptive capacity (Kusangaya et al. 2014), as the El Niño drought event of 2015/16 amply demonstrates (Archer et al. 2017;Baudoin et al. 2017). Therefore, future changes in rainfall over this region need to be identified and understood, such that stakeholders, from civil society to policymakers, can make informed decisions about future adaptation planning in key sectors (Collins et al. 2012;Knutti et al. 2010).
We focus here on southern Africa, where the rainfall climate (see Figs. S1 and S2 in the online supplemental material) is broadly characterized by a local summer wet season in austral summer, a pronounced rainfall gradient from the dry southwest to the humid north and northeast, and a peak wet season rainfall maximum extending from the continent into the adjacent southwest Indian Ocean, referred to as the south Indian Ocean convergence zone (SIOCZ; e.g., Lazenby et al. 2016), a preferred locus for tropical-temperate weather systems. For this region the multimodel mean projections from the CMIP3 and CMIP5 model ensembles show a degree of consistency in the dominant signals of change with drier austral summer conditions projected over the southwest of the region as well as a northeasterly shift in the SIOCZ and a later onset of the rainy seasons with indications of early cessation of rainfall (Niang et al. 2014;Shongwe et al. 2009).
For most of the tropics, at least at regional scales relevant to decision-making, considerable uncertainty is evident across the ensemble of global and regional models currently available in both sign and magnitude of future precipitation projections (Rowell 2012;Knutti and Sedlá cek 2013;McSweeney and Jones 2013). Central and southern Africa are no exception (see section 3) and there is a clear need for improved understanding due in part to a complex climatological setting in which regional climate drivers as well as remote influences affect the region (Christensen et al. 2013;IPCC 2007;Kusangaya et al. 2014). This presents challenges to climate adaptation policy, and the persistence of uncertainty in climate projections has led to the development of ''decision-making under climate uncertainty'' approaches in adaptation (e.g., Hallegatte et al. 2012). These approaches appreciate the cascade of uncertainty typically associated with climate risk assessment, resulting from climate projections, downscaling techniques, impact modeling, and so on, and emphasize appropriate methods for managing risk under such conditions including, for example, robust or low regrets options or adaptive management (e.g., Dessai et al. 2009).
There is therefore considerable interest in improving our understanding of physical mechanisms driving the particular patterns of projected model changes, so that we may determine and potentially improve the robustness and credibility of projections. The physical processes driving projected change are numerous and may vary among models. Projection uncertainty is a result of varying processes found within models including, among others, parameterization schemes, climate sensitivity, and regional patterns of SST changes.
Based on analyses of projections, various mechanisms have been proposed that link increased global temperatures and precipitation changes, notably the thermodynamic ''wet get wetter'' process of rainfall change Soden 2000, 2006;Allen et al. 2010;Christensen et al. 2013;Chou and Neelin 2004;Meehl et al. 2007;Chou et al. 2009;Seager et al. 2010). This operates at the largest scales (e.g., zonal means) and involves the increase in global specific humidity in a warmer atmosphere leading to an increase (decrease) in precipitation in the regions of mean moisture convergence (divergence). In the tropics this is likely to be offset by the weakening of the mean tropical overturning circulation associated with a reduction in convective mass flux in regions of present-day high ascent (Vecchi et al. 2006;Chadwick et al. 2013;Christensen et al. 2013;Ma and Xie 2013;DiNezio et al. 2013). Chadwick et al. (2013, hereafter C13), however, proposed that the wet-get-wetter mechanism alone does not adequately explain the global pattern of multimodel mean (MMM) projected rainfall change. The spatial correlation of future precipitation change DP and mean precipitation P globally is low, such that at regional and seasonal scales in the tropics other processes dominate. These processes are substantially related to changes in the spatial location of moisture convergence and hence convection. These include dynamic effects of regional gradients in near-surface temperature change over oceans, that is, the ''warmer get wetter'' process (Xie et al. 2010), land-sea temperature contrasts (Dong et al. 2009;Byrne and O'Gorman 2013), land surface processes (Pitman 2003), aerosol direct, indirect, and semidirect effects (Huang et al. 2007;Lohmann and Feichter 2005;Ackerman et al. 2000;Hansen et al. 1997), and changes in circulation (Shepherd 2014). The thermodynamic balance of the ''upped ante'' (Neelin et al. 2003) mechanisms additionally contributes, leading to results such as the so-called modified warmer-get-wetter mechanism (Huang et al. 2013). The modified warmer-get-wetter mechanism is described as the combination of SST changes (the warmer-get-wetter effect) modified by background climatological moisture and SSTs due to the nonlinear relationship between tropical convection and SSTs.
Generally, the processes operating over the ocean are better understood than those over land but in all cases differing representation of these processes across models is likely to drive projection uncertainty. Most previous analyses have focused on the MMM projected change quantities; however, Rowell et al. (2015) highlighted the importance of understanding the mechanisms of change within individual models and concluded that further investigation should be aimed at developing expert judgment of process-based mechanisms and their reliability of projections (Rowell et al. 2015).
In this context, the aims of this paper are 1) to determine the mechanisms of projected regional precipitation changes by decomposition into thermodynamic and dynamic components; 2) to quantify the contribution to total ensemble projection uncertainty, as represented by intermodel spread, associated with these mechanisms; and 3) to identify possible causes of uncertainty, and to draw inferences regarding the robustness and credibility of projected changes. The mechanism of change decomposition method of C13 is used to identify causes of precipitation change (aim 1) and associated uncertainty (aim 2) (reported in section 3a) and causes of uncertainty are inferred (aim 3) through analysis of the intermodel spread (reported in section 3b). This paper focuses on projected regional precipitation changes over southern Africa (SA) and the adjacent southwest Indian Ocean (SWIO) sector (08-308S, 108E-808E), where wet season rainfall is dominated by the SIOCZ (Cook 2000;Lazenby et al. 2016).

a. Data
We use output from simulations of the twentieth century and the twenty-first century [under the representative concentration pathway 8.5 (RCP8.5) emissions scenario] from 20 models [those used in Kent et al. (2015, (Meehl et al. 2007;Taylor et al. 2012) (see Table 1). We consider both the MMM of the 20 chosen CMIP5 models and the spread of model projections across the 20-model ensemble, to address aims 1 and 2 and aim 3, respectively. Lazenby et al. (2016) present the biases in the MMM of the larger full CMIP5 ensemble and we note that this is not significantly different from that of the MMM of the 20 CMIP5 models used here (see Figs. S1 and S2). Both MMMs represent well the dominant circulation features and the spatial structure of rainfall and its annual cycle. However, there is an overall positive precipitation bias compared to observations over the main climatological rainfall features over the Indian Ocean and central continent in austral summer, although not in the summer rainfall season over the Sahel zone (;158N).
Monthly data were extracted for key diagnostic variables to understand potential physical processes linked to precipitation changes over the SA-SWIO sector. The period of analysis for projected climate changes is 2071-2100 (in the RCP8.5 experiment) minus the historical period 1971-2000. Only the first ensemble member was utilized in creating the MMM. All model data were interpolated to a common grid of 1.58 3 1.58 to ensure uniformity. To determine the main contributors to both the MMM precipitation change (see section 3a) and associated uncertainty (i.e., intermodel spread; section 3b) over the study domain, the C13 decomposition of change methodology was applied to data from each model for each calendar month. We then focus specifically on the early and main SA-SWIO wet seasons October-December (OND) and December-February (DJF), respectively.
Projected changes in precipitation can be decomposed into the dominant mechanisms of change, the thermodynamic and dynamic components. A number of methods have been proposed for the purpose (e.g., Seager et al. 2010;Emori and Brown 2005). We utilize the method of C13, which is based on the assumption that in convective climate regimes mean precipitation P is equivalent to the vertical mass flux from the boundary layer to the free troposphere M multiplied by specific humidity in the boundary layer q : (1) As M (as defined here) is not directly available from most CMIP5 model outputs, in the decomposition a suitable surrogate M* is derived directly from model mean P and q:under C13 demonstrated M* to be a suitable replacement for actual column-integrated convective mass flux M int at the gridpoint level in global climate models. Therefore, the projected precipitation change DP can be expressed as where M* and q are the present day mean climatological values  of proxy mass flux and 2-m specific humidity respectively and DM* and Dq are the projected changes in those quantities over the period 2071-2100 compared to 1971-2000. In this formulation Dq and DM* represent the thermodynamic and dynamic components of change, respectively. The C13 method has the advantage that the dynamical term DM* can be further separated into DM weak * and DM shift * . The quantity DM weak * represents the tropics-wide weakening of the large-scale overturning circulation. It is derived from the climatological mean M* (DM weak * 5 2aM*), where a is a constant derived for each model, separately from the strong negative relationship observed between climatological M* and DM* across all grid cells in the tropics (assumed here to include all regions between 308N and 308S), that is, areas with higher M* experience greater declines in DM*, common to all models (see Figs. 3 and 6 from C13). The weakening in tropical circulation is due to the positive lapse-rate feedback response in which greater upper-level latent heating interacts with the mean circulation, leading to greater column warming in descent regions than ascent regions, thereby decelerating the tropical circulation (Ma et al. 2012). There is also a contribution to circulation weakening from the direct radiative effect of increased CO 2 concentrations (e.g., Bony et al. 2013).
Note that DM shift * 5 DM* 2 DM weak * ; that is, the deviation in DM* from that estimated directly from the regression of M* and DM* represents the effective ''shift'' at a given grid cell toward a greater or weaker convective mass flux.
The term Dq can also be further separated into two components: 1) Dq cc , the Clausius-Clapeyron change in surface q for the change in mean 2-m temperature at each location (expected under fixed relative humidity), expressed as Dq cc 5 (es 2 /es 1 )q 2 q, where es is the saturation vapor pressure and q is historical boundary layer specific humidity [the ratio expresses future (es 2 ) by historical (es 1 ) saturation vapor pressure], and 2) Dq rh , the residual of Dq 2 Dq cc , associated with changes in near-surface relative humidity.
On this basis DP for each model can be decomposed into its individual components as DP 5M*Dq cc 1 M*Dq rh 1 qDM weak * 1 qDM shift * 1 DqDM*, which for convenience in terminology can then be expressed as DP 5 DP T 1 DP RH 1 DP Weak 1 DP Shift 1 DP Cross , where d DP is the change in precipitation expressed per degree global warming, d DP T 2 M*Dq cc is the thermodynamic change due to Clausius-Clapeyron-driven increases in specific humidity, d DP RH 5 M*Dq rh is the change due to near surface relative humidity changes, d DP Weak 5 qDM weak * is the change due to the weakening tropical circulation, d DP Shift 5 qDM shift * is the change due to spatial shifts in the pattern of convective mass flux, and d DP Cross 5 DqDM* is the cross component of precipitation change, associated with interactions between the other components.
The total thermodynamic component consists of DP T 1 DP RH . For simplicity, we combine the total thermodynamic component (DP T 1 DP RH ) with the dynamic component associated with the weakening of the tropical circulation (DP Weak ) to create a new component named DP twrh [hereafter named the ''thermodynamic residual'' term, from Kent et al. (2015)]. We include DP Weak with the thermodynamic terms due to the strong anticorrelation between these components. This analysis of the mechanisms of change as applied to the ensemble MMM is presented in section 3a.

c. Analysis of uncertainty
To quantify model uncertainty in projected change we assess the spread in the precipitation change components over the 20-model member ensemble (section 3b). The dominant spatial pattern of variation in DP Shift between models is identified through EOF analysis, applied to the cross-model DP Shift fields (standardized) for two separate domains within our study region: the southern African continental land region (108-308S, 108-408E) and the Indian Ocean (08-308S, 408-808E) region, (see section 3c). Potential drivers of this uncertainty are then assessed through correlations of the DP Shift EOF component scores with diagnostic fields (SST and circulation indices) across the model ensemble and through composite analysis of diagnostic fields from models, sampling the upper and lower 25% of models from the EOF scores (see section 3c).
Note that all components of projected change in this analysis are normalized by the mean global surface temperature change DT global of each individual model, and therefore all quantities are expressed as per degree of global warming. This removes the uncertainty due to intermodel spread in climate sensitivity and makes results scalable to the magnitude of warming (assuming quasi-linearity of regional precipitation change with warming) and therefore tractable when applied to predefined warming levels for policymakers. Removing the effect on DP of uncertainty in model DT global related to model climate sensitivity (i.e., the high correlation between normalized and absolute precipitation changes; see Fig. 1) indicates a minor influence on our results.

Results and discussion
a. Changes in multimodel mean precipitation over southern Africa and the adjacent Indian Ocean First, we consider the annual cycle of zonally averaged MMM precipitation changes over our study domain for land and ocean regions (Figs. 2a and 2b, respectively). The most pronounced feature of projected changes in rainfall over both land and ocean is an opposing dipole of future wetter and drier conditions, oriented to the north and south, respectively, of the climatological axis of maximum rainfall. However, the drying-wetting dipole is not obvious in January-February (Fig. 2a), particularly over land. This DP wetter-drier dipole structure effectively straddles the ITCZ with wetting (feature A in Figs of the ITCZ axis. This indicates an effective northward shift in the ITCZ. This dipole pattern, however, has a strong seasonal cycle peaking in austral spring (September-December) over both land and ocean. The OND season is often overlooked in studies of SA-SWIO climate in favor of the main rainy season, DJF, but here we note a remarkably strong rainfall change signal in OND. The spatial patterns of MMM precipitation change over our study region for the key seasons OND and DJF (Figs. 3 and 4 ) clearly illustrate this dominant zonally oriented pattern of wetter (drier) conditions located to the north (south) of the mean ITCZ at ;108S, apparently connecting precipitation changes over the southern African continent with those in the equatorial and southwest Indian Ocean. The continental drying is most pronounced in OND (Fig. 2a, feature B) centered at ;158S over the southern African continent. This may indicate a later start and reduced length of the growing season rainfall with serious implications on the agricultural sector. In addition to the wetting-drying dipole there is evidence of additional wetting over 258-308S over eastern South Africa and across into the southwest Indian Ocean, especially in DJF.
The decomposition can aid in understanding the processes driving these changes. The thermodynamic component DP T (Figs. 3d and 4d) results in a wetting signal everywhere but whose magnitude is proportional to a rise in temperature through the Clausius-Clapeyron relation (i.e., it represents the wet-get-wetter process and maps substantially onto mean rainfall). It is substantially offset by an equivalent pattern of drying from the weakening of the tropical overturning circulation DP Weak (Figs. 3e and 4e), which is inversely proportional to M * and as such maps onto mean P. In addition, over land wetting from DP T is offset by drying from DP RH (Figs. 3f and 4f), due to reduced relative humidity, presumably resulting from moisture supply not keeping pace with increasing temperature. Together, these terms constitute the thermodynamic residual term DP twrh (Figs. 2e, 2f, 3g, and 4g), which drives a net wetting signal, peaking over the oceanic climatological ITCZ and humid land regions. Therefore, DP is partly composed of a thermodynamic wet-get-wetter process magnifying the mean ITCZ rainband. However, the thermodynamic DP RH response also drives up to ;50% of the drying over subtropical land regions during OND (feature C in Fig. 2e) centered on ;208S. This is broadly coincident with the maxima in DT (Fig. 2g, feature D), suggestive of a land-atmosphere positive feedback response in early wet season, which is overcome during the peak DJF wet season. As such over land in OND DP twrh is an important driver of the spatial pattern of DP (r 5 0.70, where r is the spatial correlation), substantially associated with DP RH.
However, over the study region as a whole, it is clear from Figs. 2-4 and Table 2 that the spatial pattern of DP in OND or DJF is not closely related to that of mean precipitation (P) (spatial correlations are not significant at the 95% confidence interval) or to the thermodynamic residual term DP twrh (r 5 0.46 and 0.30 for OND and DJF, respectively), such that the notable dipole features of DP are not driven by the wet-get-wetter process. Rather, DP most closely matches the dynamical component DP Shift (r 5 0.94 for OND and 0.96 for DJF) such that it is changes in the location of convection that explain the dipole. Note that DP Shift contributes most of the wetting-drying dipole over land and almost all of it over ocean, especially the drying signal south of the ITCZ. We postulate here that this projected northward shift of the ITCZ over the Indian Ocean may be related to projected changes in the SST structure through a ''warmest get wetter'' mechanism. The north-south gradient in SST over the Indian Ocean, which is broadly representative of the DP dipole [;(58N-208S)] in the present day of ;(2-4) K over the austral summerwinter seasons (Fig. 2h), is enhanced in the MMM projections by ;0.5 K (Figs. 2h, 3h, and 4h). This may lead to a northward displacement of convection toward the warmest oceanic waters, similar to that noted for the tropical Pacific (Widlansky et al. 2013).

b. Quantifying uncertainty in projected change in precipitation: Ensemble spread
Uncertainty across the multimodel ensemble remains a feature in projections of future precipitation across the tropics and is a major barrier to effective use of climate information in adaptation activities, notwithstanding approaches to ''decision making under climate uncertainty'' (e.g., Lempert and Collins 2007). Here, we are averaged over longitude bands indicative of (left) the SA continent (108-408E; land only) and (right) the Indian Ocean assess the contributions to total ensemble uncertainty from the various mechanisms of change in our decomposition. We assess the robustness of projected changes from the magnitude of intermodel spread for each component as represented in two forms: first as maps of the standard deviation of the ensemble at each grid cell (Figs. 5 and 6) and second as box-and-whisker plots of intermodel spread for each component averaged over specific regions of interest (Fig. 7). We average over areas where the MMM change signal shows mean future wetting or drying to help inform interpretation of the robustness of key signals emerging from MMM commonly used (e.g., Flato et al. 2013;IPCC 2013). Further, we consider land and ocean separately given the differing level of importance for adaptation actions and differing mechanisms driving changes. A number of broad signals emerge.
(i) Changes in precipitation averaged over the regions of wetting-drying from the MMM are more robust over land than ocean for both the future wetting and drying signals (i.e., total uncertainty in DP is lower for land compared to ocean; Fig. 7). This difference in the spread of DP is largely related to that in DP Shift (r 5 0.95 for OND and 0.97 for DJF for land regions), such that uncertainty in the dynamical drivers of projected change is dominant. We consider the contribution of intermodel differences in SST structures to uncertainty in DP Shift in section 3c. (ii) The MMM projected drying is more robust than wetting over both land and ocean and in OND, but differences are small in DJF. The future signal with the greatest robustness is the continental drying over (especially western) SA during the OND season (Figs. 5 and 7 and 6a), which raises the possibility that intermodel differences in background climatology may project onto the uncertainty in DP. However, only some of these hotspots of uncertainty correspond to locations of high intermodel spread in mean P (Figs. 5h and 6h), such that most of the spatial pattern in uncertainty in DP is unrelated to background climatology. The high uncertainty in DP often apparent in wetting-drying transition zones suggest that caution needs to be attached to interpretation of near-zero MMM DP change. (iv) Most of the total uncertainty in DP is contained in the DP Shift component, whose uncertainty is typically greater by a factor of 2-4 than that of the thermodynamic residual. Intermodel standard deviation within the DP Shift term is markedly higher over ocean than land, notably over the DP dipole region. While the MMM wetting-drying dipole over the Indian Ocean may result from a warmestget-wetter response to changes in the SST structure (Figs. 5h and 6h), high uncertainty in DP Shift suggests strong intermodel divergence in the form of the DSST patterns responsible (Chadwick 2016), which we explore in section 3c. Indeed, in many regions, notably the region of projected drying over the Indian Ocean at 108-158S, the model uncertainty in DP Shift is higher than that in DP. Note that DP twrh is especially robust over the ocean ITCZ where the thermodynamic response is closest to that of DP T . (v) Since the total uncertainty can be less than the sum of the individual components (e.g., Fig. 7) we can infer that the intermodel components are anticorrelated across models and offset each other, which acts to constrain total uncertainty in DP. For the drying over continental SA (during OND) and to a lesser extent over the Indian Ocean, uncertainty in the dynamical DP Shift and thermodynamic residual term DP twrh terms appear to offset each other, resulting in a total intermodel standard deviation of DP, which is much less than the sum of the component terms. (vi) The thermodynamic change due to Clausius-Clapeyron-driven increases in specific humidity (DP T ) exhibits the largest magnitude of all the decomposed components (Fig. 7); however, this is largely offset by the weakening of the tropical circulation component (DP Weak ) and relative humidity component (DP RH ) over land regions.

c. Understanding potential causes of uncertainty in the dynamic component of projected change
The dynamical DP Shift term provides the primary contribution to intermodel uncertainty in future projections of precipitation. Among the potential causes of this uncertainty in the location of tropical convection, a primary candidate is the varying patterns of model projected SST changes. We explore this using EOF analysis of the intermodel DP Shift patterns. Over the Indian Ocean domain during both OND and DJF the leading mode of intermodel DP Shift variability shows a loading pattern (Figs. 8a,b) that projects strongly onto the pattern of the future wetting-drying dipole in DP Shift (Figs. 3b and 4b), and indeed in DP explaining 37% and 31.5% of variability, respectively. These EOFs therefore represent well the strength of the MMM DP wettingdrying dipole, oriented broadly north-south, in individual models.
The component scores of this leading EOF correlate moderately with intermodel DSST over the Indian Ocean in both OND and DJF (Figs. 8c,d). In both cases the correlation is such that a stronger north-south wetting-drying dipole in individual models is associated with lower rates of SST warming in the Indian Ocean south of the equator, and hence an increased north-south SST gradient across the Indian Ocean. This response across models is consistent with the warmest-get-wetter mechanism inferred from the MMM changes in precipitation and SST (section 3a).  Moderate correlations also exist with DSST in the tropical Pacific broadly around the Niño-3.4 region. During OND there is an additional indication that the DP Shift EOF dipole is related to changes in the mean east-west SST gradient reminiscent of the Indian Ocean dipole (IOD) mode of interannual variability active in this season (Saji et al. 1999). Previous analysis has related changes in long-term mean precipitation to a shift toward a preference for the positive mode of the IOD in some coupled models (Shongwe et al. 2011), Over the domain centered on the SA landmass the intermodel DP Shift EOF analysis for both the OND and DJF seasons reveals a dominant EOF loading pattern oriented diagonally northwest-southeast across SA toward the southwest Indian Ocean (Figs. 9a,b). The peak loadings lie to the southwest of the position of the mean SIOCZ feature such that the EOF represents the magnitude of the MMM drying signal to the south of the mean rainfall maximum. The MMM projected pattern of DP Shift change across the region is more zonally oriented than the EOF SIOCZ-like orientation (Figs. 3b and 4b; see section 3a) such that the MMM change clearly masks considerable spread within the ensemble.
The EOF pattern of intermodel DP Shift appears associated with intermodel variability in SST, land heating and circulation and resulting moisture fluxes based on both correlations (Figs. 9c,d) and composites of diagnostic fields based on samples of the upper and lower 25% of models from the EOF component scores (Fig. 10). The EOF is rather weakly correlated to the intermodel structures of projected changes to SST across the southern Indian Ocean with an east-west dipole structure of correlation with the EOF component scores along ;(208-25) 8S (Figs. 9c,d) and corroborated in the composites (Fig. 10). For both the OND and DJF EOF, stronger future drying over SA in models is associated with weaker future warming in the Mozambique Channel-southwest Indian Ocean and stronger warming in the eastern subtropical Indian Ocean, as well as the eastern South Atlantic.
While the spatial patterns are similar, the correlations of EOF and SST are notably stronger in DJF than OND (Figs. 8c,d and 9c,d). The implication is that SSTs are most significant in austral summer than the transition season over the southern African region, particularly over the South Atlantic Ocean and central Indian Ocean. OND rainfall does not exhibit sufficient evidence of a significant link with changing SST gradients. Reason et al. (2006) state that the onset of the wet season over southern Africa appears to be associated with anomalous ridging of the South Atlantic high pressure, which requires further investigation and may shed more light on this research. Circulation features associated with the intermodel spread are identified by compositing projected changes in diagnostic fields in the top and bottom 25% of models from the EOF coefficients (Fig. 10). We find that in both seasons more intense versus weaker drying over SA in models is associated with cyclonic low-level circulation over the southwest Indian Ocean (Fig. 10) around a relatively weaker subtropical high, resulting in weaker Indian Ocean easterlies, and therefore an effective northeastward shift in the axis of the SIOCZ (Lazenby et al. 2016). In DJF (Fig. 10b) intense drying is also favored in models in which strong land heating occurs at ;208S. This may be a result of reduced soil moisture but the heating drives an apparent intensification of the continental low-pressure center (the ''Angolan low'' feature) and cyclonic low-level circulation. While this might be expected to favor rainfall over the continent, the composite of the drier models (not shown) indicates a marked anomalous southward shift of the low with implications for moisture transport. There is some indication of an intensified South Atlantic anticyclonic circulation in the models with more intense drying whose link with the SST changes requires further investigation.
These structures of intermodel differences in DSST, land temperature, and low-level circulation that may drive the leading pattern of intermodel spread in DP Shift over SA have some resonances with modes and structures of interannual variability in the region. Notable among these are the following. (i) The first is the regional response to ENSO characterized by the displacement in the SIOCZ associated with cyclonic lowlevel anomalies over the southwest Indian Ocean during El Niño events (Cook 2001), similar to that in Fig. 10. However, it is cautionary to note that the physical mechanisms and influence of local Indian Ocean SST (e.g., Goddard and Graham 1999) versus remote Pacific SSTs (e.g., Cook 2001;Ratnam et al. 2014) in driving this in the observed climate remain to be fully resolved. (ii) The second feature is the south Indian Ocean dipole (SIOD) mode of east-west SST gradient in the subtropical Indian Ocean (Behera and Yamagata 2001). In this case, although the form of the association of DP Shift over SA and DSST matches that of interannual variability in SA rainfall and the phase of the SIOD (Reason 2001), the circulation responses diverge (Fig. 10). (iii) Finally, there is the intensity of the Angolan low and associated circulation changes, known to be an important control on present day precipitation in observations (e.g., Manhique et al. 2011) and in models (Munday and Washington 2017). During wet (dry) years the strength of the Angola low tends to be stronger (weaker) with enhanced (suppressed) circulation bringing in moisture from the Atlantic and the Indian Ocean (Cook et al. 2004;Munday and Washington 2017). Differential rates of future warming over land versus ocean may be expected to intensify and displace the Angolan low and we see evidence of differential response across models strongly related to intermodel uncertainty in the structure of future dynamical precipitation responses. Of course, we would not expect the dynamical processes of intermodel spread in projected rainfall change to be consistent with all aspects of current interannual variability, given (i) the complex suite of interannual modes affecting the SA region and (ii) the highly variable and mixed ability of coupled models to represent the mean state and variability, both through local processes (Lazenby et al. 2016) and remote teleconnections (e.g., Rowell 2013). Nevertheless the link between processes of climate variability and projected change has proved to be a fruitful focus in recent years in explaining why models make the future changes they do, providing some basis for assessing potential credibility of projected change.

Summary and conclusions
Human influence is expected to drive considerable changes to the hydrological cycle, of particular concern in regions such as southern Africa that are currently vulnerable to climate. There is a need to quantify and understand the projections of future precipitation change and associated uncertainty to inform the possible use of climate information in adaptation planning (Hackenbruch et al. 2017). Our analysis of end-oftwenty-first-century projections from a sample of 20 models from the CMIP ensemble reveals a dominant dipole pattern of precipitation change over SA and the SWIO with a wetting (drying) response to the north (south) of the ITCZ, with therefore an effective northward shift in the ITCZ. Pronounced drying is exhibited over land in OND, implying a delay in the onset of the wet season, an important and relatively robust signal.
A decomposition framework is applied to separate the mechanisms of change in both the MMM and in individual models to shed light on the relative importance of the multiple, coincident, and often competing mechanisms that may drive future precipitation changes in a warmer world. We show that the dynamical component DP Shift representing spatial shifts in convection explains most of the dipole structure of DP, as evidenced by high spatial correlations no less than 0.94 over the study domain. The thermodynamic component of increased moisture is offset by the weakening of the tropical circulation leading to a relatively weak wet-get-wetter contribution. Considering the robustness of projected change, represented by the intermodel dispersal in response, we find that DP Shift holds the most uncertainty, whereas the thermodynamic component is typically more robust. Drying is more robust than wetting over the continent, notably in the early summer OND season.
For both the OND and DJF seasons MMM future precipitation projections show, at the broadest scale, a wetting-drying dipole across the continent and Indian Ocean. Over land the magnitude, extent, and robustness of drying is greater during in OND than DJF, with implication of a later onset in the wet season over southern Africa. Greater uncertainty during DJF is apparently associated with a stronger role of intermodel differences in the structure SST changes in the adjacent oceans.
Over the Indian Ocean sector, we relate the dominant dynamically driven DP Shift component to changes in SST structures. The wetting-drying dipole is associated with patterns of SST change, both in the MMM and across the model ensemble, that is consistent with a warmest-getwetter mechanism driving the northward shift in the location of convection and convergence. Much of the model uncertainty is therefore associated with intermodel differences in future SST patterns in the tropical Indian and Pacific Ocean and the atmospheric teleconnection response, such that better understanding of these projected changes and in the teleconnections linking these to rainfall may provide some basis for constraining uncertainty in projections. This inference, however, is particularly applicable to DJF as OND does not exhibit the same magnitude of significance over the region.
Changes in rainfall over land are likely more complex, first because over the continental interior of SA the contribution to mean drying of the thermodynamic reduction in relative humidity is roughly equal to that of the dynamical component. Reduced relatively humidity is likely itself to be related to the effects of enhanced land-sea temperature contrasts and physical mechanisms for suppression of convection through raised lifting condensation levels have been proposed (e.g., Fasullo 2010). Second, regarding the dynamical component of change, over the SA landmass we find intermodel uncertainty to be associated with low-level circulation patterns associated with zonal gradients in SST changes in the subtropical Indian Ocean and the intensity of continental heating and the thermal low pressure center. The land-sea heating contrast has been hypothesized by Bayr and Dommenget (2013) to be a driver of enhanced convergence over land with implications for convection at the continental scale, but our results indicate that the response can be complex in this case leading to a northward shift of the ITCZ/SIOCZ and wetting-drying at the regional scale. Over both the ocean and land examination of intermodel dispersal yields structures with less zonal uniformity of change than represented by the MMM response. The sensitivity of DP Shift across models to the pattern of SST structures and their resemblance to contemporary modes of variability suggests that further analysis of future SST patterns and assessment of model representation of present rainfall-SST teleconnections and mean state biases may be useful to inform our interpretation of the credibility of projected changes.
This paper provides new insight of projected precipitation mechanisms of change over southern Africa by identifying the dominant pattern of change over southern Africa (the wetting-drying dipole). Drivers of this projected change are evaluated and attempts are made to quantify the associated uncertainty. Regional circulation analysis is additionally performed to shed light on the patterns driving future dynamic precipitation changes in OND and DJF, which are in agreement with previous studies over other regions (e.g., Kent et al. 2015).