Revisiting the 1992 severe drought episode in South Africa: the role of El Niño in the anomalies of atmospheric circulation types in Africa south of the equator

During strong El Niño events, below-average rainfall is expected in large parts of southern Africa. The 1992 El Niño season was associated with one of the worst drought episodes in large parts of South Africa. Using reanalysis data set from NCEP-NCAR, this study examined circulation types (CTs) in Africa south of the equator that are statistically related to the El Niño signal in the southwest Indian Ocean and the implication of this relationship during the 1992 drought episode in South Africa. A statistically significant correlation was found between the above-average Nino 3.4 index and a CT that features widespread cyclonic activity in the tropical southwest Indian Ocean, coupled with a weaker state of the south Indian Ocean high-pressure. During the analysis period, it was found that the El Niño signal enhanced the amplitude of the aforementioned CT. The impacts of the El Niño signal on CTs in southern Africa, which could have contributed to the 1992 severe drought episode in South Africa, were reflected in (i) robust decrease in the frequency of occurrence of the austral summer climatology pattern of atmospheric circulation that favors southeasterly moisture fluxes, advected by the South Indian Ocean high-pressure; (ii) modulation of easterly moisture fluxes, advected by the South Atlantic Ocean high-pressure, ridging south of South Africa; (iii) and enhancement of the amplitude of CTs that both enhances subsidence over South Africa, and associated with the dominance of westerlies across the Agulhas current. Under the ssp585 scenario, the analyzed climate models suggested that the impact of radiative heating on the CT significantly related to El Niño might result in an anomalous increase in surface pressure at the eastern parts of South Africa.


Introduction
South Africa, situated at the descending branch of the Hadley circulation, is prone to drought (Baudoin et al. 2017). However, the 1992 drought episode in South Africa is characterized as one of the most severe in history (Glantz et al. 1997;Mason and Tyson 2000). Its occurrence during the 1991/1992 growing season significantly impacted agricultural activities and food security, resulting in the demand for increased efficiency in the management of drought events in South Africa (Bruwer 1993). Moreover, under greenhouse gas-induced climate change, dry spell is expected to increase in South Africa (IPCC 2013). Hence, using the 1992 severe drought episode as a reference, which equally is recorded as a strong El Niño year (Glantz et al. 1997), this study investigates how the interrelation between the El Niño Southern Oscillation (ENSO) and circulation types (CTs) in Africa south of the equator might have contributed to the drought episode and how the analysis might be beneficial for enhanced understanding of the synoptic-scale causes of drought events in South Africa.
Among the different teleconnection signals in the south Indian Ocean, ENSO plays a significant role in the interannual variability of rainfall in South Africa . ENSO has two phases-El Niño and La Niña-that can influence atmospheric circulations and patterns of precipitation in southern Africa (Jury et al. 1994;Reason et al. 2000). Strong El Niño years are usually accompanied by below-average rainfall in most regions in southern Africa and vice versa for La Niña years. The origin of ENSO is related to tropical Pacific sea surface temperature (SST) anomalies, and the energy released by the resulting moderations in convection in the tropical Indian Ocean and the Pacific Ocean induces Rossby waves over southern Africa (Hoell et al. 2015). According to Reason and Jagadheesha (2005), the regional fluxes of moisture in southern Africa can be modulated by the excited Rossby waves. Manatsa et al. (2015) noted that the relationship between southern Africa rainfall and ENSO is strongest in December-March (DJFM). The South Indian Ocean Convergence Zone (SICZ) is the major synoptic feature that modulates the austral summer (DJF) rainfall variability in southern Africa (Cook 2000;Ninomiya 2008). Cook (2000) reported that during strong El Niño events, the SICZ undergoes a northeastward shift. Thus, ENSO might be able to constrain the major regional rain-bearing atmospheric circulation features in southern Africa that result in the formation of the SICZ-such as the South Indian Ocean high-pressure and the Angola low. Hoell et al. (2015) reported that in southern Africa, vertical motions can be modified by ENSO events. Also, a weakening of the South Indian highpressure and the Angola low might be expected during strong El Niño events (Cook 2000;Reason and Jagadheesha 2005).
It is also necessary to note that given the remote influence of ENSO, the predictability of drought in South Africa is not always well correlated with El Niño (Reason and Jagadheesha 2005). For example, during the 1997/1998 El Niño event, rainfall reduction was projected in southern Africa but this failed to materialize (Lyon and Manson 2007). However, SST anomalies in the southwest Indian Ocean, Benguela current, the Greater Agulhas region, and variations in synoptic systems such as the South Indian Ocean high-pressure and the South Atlantic Ocean high-pressure influence moisture uptake in the adjacent oceans, moisture transport to South Africa, and so rainfall formation in the region (Walker 1990;Reason and Mulenga 1999;Vigaud et al. 2009;Driver 2014).
Several studies have used the Self-Organizing Maps (SOMs) to identify synoptic circulation over southern Africa (e.g., Lennard and Hegerl 2015;Engelbrecht et al. 2015, Engelbrecht and Landman 2016, Wolski et al. 2018Odoulami et al. 2020). The capability of the SOM to generate non-linear solutions and its tendency to treat the data set as a continuum make it a better choice for circulation typing compared to hard clustering algorithms (e.g., K-means) that allow the classified variable to be assigned only to one class. However, given the continuous and fuzzy nature of atmospheric circulation patterns, a classification scheme that allows overlapping of the classified variable can result in relatively better classifications (Gong and Richman 1995;Xu and Tian 2015). The rotated principal component analysis technique fits such classification criteria since it is fuzzy by the nature of the PCA loadings. When CTs are defined using obliquely rotated PCA applied to a climate variable represented in T-mode (i.e., variable is time series and observations are grid point), the classified CTs can exhibit physically interpretable dipoles (Compagnucci and Richman 2008) related to the probability of specific weather events in the study region (Richman 1981). The major difference between the SOM and obliquely rotated T-mode PCA is that SOM classification is based on model (see Xu and Tian (2015) and Hewitson and Crane (2002) for detailed explanation), while the rotated PCA technique finds vector positions that maximize hyperplane counts and simultaneously find subsets of highly related variables (Richman 1981(Richman , 1986Richman and Gong 1999).
In this study, CTs are defined using obliquely rotated T-mode PCA. Additionally, it is investigated if there is a modulation of the ENSO signal over the occurrence of the CTs. The physical interpretability of CTs that their amplitude or frequency are related to or influenced by ENSO will be examined by investigating if they reflect large-scale patterns of atmospheric circulation known a priori from existing works of literature that might be expected during ENSO events. During the 1992 El Niño season, the physical and statistical nature of the relationship between the selected synoptic situations and ENSO are investigated to understand how El Niño might constrain rain-bearing CTs in South Africa during the austral wet seasons in southern Africa. Thus, the strengths of this study are on the improved understanding of the mechanisms through which El Niño causes drought in South Africa in the light of anomalies in CTs in addition to highlighting how the statistical concept of circulation typing, achieved with obliquely rotated PCA, can be extended to the predictability of regional rainfall in southern Africa. The latter is of importance given the increased vulnerability of the hydroclimate of South Africa, and the need for simple approaches that can be applied by local weather agencies and water decision-makers in rainfall prediction in the study region.

Data and methodology
Sea level pressure (SLP); precipitation; specific humidity; and wind vectors at 850 hPa data sets are obtained from NCEP-NCAR reanalysis (Kalnay et al. 1996). The data sets are obtained at horizontal and temporal resolutions of 2.5° longitude and latitude, and daily from 1979 to 2020, respectively. The precipitation data set is obtained at a Gaussian grid type and interpolated to 0.75° longitude and latitude using First order conservative remapping (Jones 1999). Precipitation data set obtained from the Climate Prediction Center (CPC) (Xie et al. 2007) at a horizontal resolution of 0.5° longitude and latitude and daily temporal resolution from 1981 to 2010 is used to validate the NCEP precipitation data set in South Africa.
The spatial extent for the circulation typing is 5.25° E-55.25° E and 0° -50.25° S. The adjacent oceans are included to capture moisture advection processes. The classification steps with obliquely rotated T-mode PCA are the same as documented in Ibebuchi (2021a, b). Here it is applied to the daily z-score standardized SLP data set for the 1979-2020 period. The representation of the SLP data set in T-mode structure is based on the finding of Richman (1981Richman ( , 1986) that classification of a field that explains atmospheric circulation and represented in the T-mode structure results in a simplified time-series isolating subgroups of observations with a coherent spatial pattern.
In the classification process, a correlation matrix is used to relate the time series. Singular value decomposition is used in factorizing the matrix to obtain the PC scores, eigenvalues, and eigenvectors. The eigenvectors are multiplied with the square root of their corresponding eigenvalues that make them longer than a unit length, referred to as loadings (Richman and Lamb 1985). The number of components to retain is first based on the philosophy of North et al. (1982) which implies retaining components with separated eigenvalues. However, since the subsequent steps require rotation of the retained components, the recommendation of North et al. (1982) might not be so useful; however, it was used as a first guess. Preisendorfer et al. (1981) explained that even components with low eigenvalues might still uncover vital climatic patterns; thus, the ultimate decision of the number of components to retain is based on ensuring that the further added components have not been delineated by previous vectors (Richman 1981).
The retained components are rotated obliquely using Promax at a power of 2 (Richman 1986). The oblique rotation of the PCs relaxes the orthogonality constraint and maximizes the number of near-zero loadings so that each component clusters a unique number of days (Richman 1981(Richman , 1986. The absolute value of the loadings designates an important signal and threshold values within the range of 0.2-0.35 can be used to separate noise (i.e., loadings within the zero interval) from signal, depending on the effective sample size (Richman and Gong 1999). Given the large sample size, ±0.2 is used in this paper to cluster positive and negative loadings above the threshold for a given retained component. Thus, each retained component yields two clusters and the mean SLP of the days clustered under each class is the CT. This step increases the weather coherency and coherency in the time occurrence of the days clustered under each class. According to Gong and Richman (1995), another advantage that using rotated PCA offers despite the binarization with the (0.2) threshold is that it is inherently fuzzy-i.e., in this case, a day can have its loading above the 0.2 threshold in more than one retained component-thus, overlapping of the classified variables is possible. Given that atmospheric circulation is a continuum, this feature is quite desirable. Another implication of the continuum nature of atmospheric circulation patterns is that though some CT(s) might tend to exhibit some seasonality, CTs are not confined to occur at any specific season; thus, to avoid information loss on the time-frequency distribution of the CTs, the complete seasons in a year are used for the classification. However, the annual cycle of the selected CTs significantly modulated by ENSO is further analyzed to uncover their probability to be dominant at specific season(s).
Previous classifications of CTs in Africa south of the equator (Ibebuchi 2021a, b) show that for a given retained component, the dipoles (i.e., clusters of positive and negative loadings above the 0.2 threshold) can be related to the inverse of synoptic situations under different opposite seasons. For example, a synoptic situation associated with continental heating (upward vertical motion) dominates in austral summer (DJF) while its inverse exhibits subsidence and dominates in austral winter (JJA). This is also a further added value in using the complete seasons for the classification. Thus, the loadings (which designate the amplitude of the CTs) under each retained component, comprising of positive and negative dipoles, are related to the standardized Nino 3.4 index using correlation analysis. Nino 3.4 index is obtained from https:// psl. noaa. gov/ gcos_ wgsp/ Times eries/ Nino34/, and the correlation is done for the 1979 to 2020 period (excluding the years before the satellite era). A statistically significant correlation, at lag 0, based on the Kendall's-Tau-b test at a 95% confidence level (i.e., p value < 0.05), implies that ENSO is related to the component in question. Positive correlation coefficient (R) implies that above-average Nino 3.4 index (El Niño) is related to the positive (dipole) phase of the component and vice versa for a negative correlation coefficient. Also, the correlation analysis is done between the Nino 3.4 index and the annual frequency of occurrence of the CTs. To investigate if ENSO influences the amplitude of the CTs, lagged correlation of up to 3 months between ENSO and the loadings of the CTs during the DJFM season was done. For example, a 1-month lagged correlation implies correlating the seasonal NDJF mean of the Nino 3.4 index with the DJFM loadings of a given retained component for the 1979-2020 period.
For the 1979-2020 period, strong El Niño (La Niña) events for the DJF season are subjectively defined as when the index is > 0.7 (< − 0.7). ± 0.7 is used since generally, ESNO events ranked with it persisted into the JFM season, and DJFM is the period when the relationship between ENSO and southern African precipitation is the strongest. Years when the DJF index is within ± 0.3 were chosen as ENSO neutral years. Anomalies in the annual frequency of occurrence and persistence of the CTs modulated by ENSO were calculated in these years as the difference between their longterm mean and their mean observation during the selected strong ENSO and neutral ENSO years. The aim is to investigate if the ENSO signal impacts the aforementioned features of the CTs. Also, seasonal rainfall anomalies associated with the selected wet CTs modulated by ENSO are investigated for the 1979-2020 period, with more focus on the 1992 wet season.
Finally, changes in the pressure field of the CT that is statistically related to El Niño were investigated under the ssp585 scenario (2065-2100) relative to the historical experiment (1979-2014) using the two climate models pre-selected in Ibebuchi (2021a). The GCMs are MPI-ESM-LR and EC-EARTH -CC. Both GCMs simulated well, the mean shape of the CT related to El Niño. Statistical significance of the anomalies is done using the permutation test at a 95% confidence level. Figure 1 shows the seasonal climatology of precipitation in South Africa from the CPC and NCEP-NCAR data sets. Austral summer (DJF) is the wettest season and austral winter (JJA) is the driest season. During austral winter, the southwestern region of the Western Cape relatively receives rainfall by the activities of cold fronts. Southeastern coastal regions receive rainfall throughout the year. Generally, on the spatial scale, the northwestern regions are relatively arid due to proximity to the cold Benguela current while the northeastern regions are wettest due to proximity to the warm Agulhas current. From Fig. 2, December to March is relatively the wettest months. NCEP-NCAR overestimates rainfall during these months and this is evident in Fig. 1 at parts of the northeastern regions. In the austral winter months, NCEP-NCAR underestimates rainfall at southwestern regions of the Western Cape but overestimates rainfall at parts of the northeastern regions relative to the CPC data set. Generally, from Fig. 1 and Fig. 2, the annual cycle and spatial variability of precipitation from both data sets are The difference between the data sets is mostly in the magnitude of precipitation. The mean absolute error of precipitation averaged over South Africa between the two data sets for the analysis period is 0.89 mm/day. The variability of precipitation measured by the standard deviation is slightly higher in NCEP-NCAR (i.e., a bias of + 0.01 mm/day). The correlation coefficient of the daily values between the data Circulation types in Africa south of the equator. The CTs that are significantly related to the above-average (below-average) Nino 3.4 index are marked by the red (blue) frames sets is 0.7. NCEP-NCAR data set will be used for further analysis given that it provides the complete set of other climatic variables of interest in this work and also due to the longest available period of the data sets.

The relationship between ENSO and the classified circulation types in Africa south of the equator
By retaining 9 components, 18 CTs were classified and Fig. 3 shows the CTs. By visual comparison, the CTs are similar to the patterns of atmospheric circulation obtained with SOM on SLP (e.g., Wolski et al. 2018). Table 1 presents the correlation coefficient between the Nino 3.4 index, and the loadings (amplitude) of the retained components and the annual frequency of occurrence of the CTs clustered from the components, for the 1979-2020 period. For clarity, for the amplitude in Table 1, the (positive) correlation coefficient is written only for the phase directly related to the above-average Nino 3.4 index, implying that El Niño rather constrains the other phase. For the analysis period, at lag 0, a statistically significant relationship was found only between the Nino 3.4 index and the amplitude and frequency of Type 5 (CT5 + /CT5 − ). The above-average Nino 3.4 index (El Niño) is significantly related to CT5 − and the below-average Nino 3.4 index is significantly related to CT5 + . Both CTs that are statistically related to El Niño and La Niña are highlighted in Fig. 3 by the red and blue frames, respectively. A physical justification of the relationship between CT5 − and El Niño is that during El Niño the south Indian Ocean high-pressure is weakened (Cook 2000) and the tropical Indian Ocean (20°S -20°N) is anomalously warm (i.e., enhanced cyclonic activity) (Klein et al. 1999). From Fig. 3, these features are evident in CT5 − . The dominant period of CT5 − is DJFM (Appendix Fig. 11) when according to Manatsa et al. (2015) the relationship between ENSO and southern African precipitation is strongest. While CT5 + is not dominant during DJFM (Appendix Fig. 11), a strong dynamical linkage was observed between CT5 + and CT6 − that is a dominant wet DJFM pattern. CT6 − is associated with a stronger anti-cyclonic circulation at the western branch of the South Indian Ocean high-pressure, and widespread southeasterly fluxes in the southwest Indian Ocean (not shown) as observed during strong La Niña events. The dynamical linkage implied that during DJFM, the occurrence of CT5 + precedes CT6 − and this could be because both CTs feature strong anti-cyclonic circulation at the western branch of the South Indian Ocean high-pressure. From Fig. 4 and Fig. 5, under CT5 − , regardless of the weakening of southeast winds under the activities of the South Indian Ocean high-pressure, the ridging of the South Atlantic Ocean high-pressure, south of South Africa, drives abundant moisture from the warm southwest Indian Ocean into the eastern regions of southern Africa. This results in enhanced rainfall at the eastern coastal regions of South Africa (e.g., Ndarana et al. 2018). The western regions are dry since the moist southeast winds penetrate less into the western regions.
Under CT5 + , the wind at the eastern coastal region of South Africa is predominantly westerly. Thus, large parts of South Africa are dry since moisture advection by southeast winds is weakened. Also, anticyclonic activity dominates in large parts of the southwest Indian Ocean. However, south of South Africa, westerlies prevail, so that the southwestern regions of the Western Cape are wet under the activity of cold fronts. Table 2 shows that in a lag of up to 2-3 months, El Niño can enhance the amplitude of CT5 −.
From Table 1, even though the relationships are not statistically significant at a 95% confidence level, El Niño might be related to the amplitude of CT1 + and CT3 − while La Niña might exhibit some relationship with the amplitude of their inverse phases. It is also interesting that CT1 + is close to the overall mean pattern of atmospheric circulation in the study region (Appendix Fig. 10). According to Tyson and Preston-Whyte (2000), El Niño can account for about 30% of rainfall variability in southern Africa. Generally, CT1 + tends to be dominant during austral winter (Appendix Fig. 11) and overlaps with the active periods of all the classified CTs. Thus, a phenomenon such as ENSO that modulates it will have a significant influence on the overall patterns of atmospheric circulation in southern Africa. From Table 2, from a lag of up to 1 month, El Niño significantly enhances  Composites constrained (enhanced) by the El Niño signal are placed at the left (right) panels the amplitude of CT1 + . The physical justification of the relationship between CT1 + and El Niño can be because El Niño increases surface pressure over landmasses in southern Africa (Hoell et al. 2015) and also weakens easterly moisture fluxes. These features-i.e., the dominance of a high-pressure system over southern African landmasses and westerly moisture fluxes across the Agulhas current-are evident during CT1 + ( Fig. 3 and Fig. 5). As a result, CT1 + is relatively a dry CT (Fig. 4). Also, Table 1 indicates that during El Niño, the frequency of occurrence of CT1 + is likely to be enhanced. CT3 − is the same CT reported in Ibebuchi (2021a) to be significantly related to the negative phase of the Southern Annular Mode (SAM). It is associated with enhanced cyclonic activity and a band of westerlies at the mid-latitudes, and vice versa for its inverse phase (Fig. 5). While not statistically significant, Table 1 indicates that El Niño might strengthen/constrain CT3 − /CT3 + . According to Wang and Cai (2013), El Niño contributes to the negative SAM and La Niña contributes to the positive SAM. It is also interesting that CT3 + is close to the austral summer mean pattern of atmospheric circulation in southern Africa (Appendix Fig. 10 and Fig. 11). Thus, its modulation by a phenomenon such as El Niño can impact significantly easterly moisture fluxes in the southwest Indian Ocean. Since generally, CT3 + overlaps with other austral summer patterns, the stronger and persistent anti-cyclonic circulation at the western branch of the South Indian Ocean high-pressure associated with it is necessary to create  Fig. 4 and Fig. 5, Table 1 and Table 2, El Niño enhances (constrains) the amplitude and frequency of occurrence of dry (wet) CTs associated with westerly (easterly) fluxes over the southwest Indian Ocean (specifically the Agulhas current). However, the statistically significant relationship between El Niño and CT5 − from Table 1 might suggest that CT5 − plays the foundation role in how El Niño modulates the other CTs. Since CT5 − indicates a weaker South Indian Ocean highpressure and enhanced cyclonic activity in the southwest Indian Ocean that during strong El Niño events can favor weakening of onshore moisture fluxes by southeasterlies. When the moistening of the boundary layer in the southern African landmasses is weakened relative to the ambient atmosphere, then deep convection (vertical motion) can be relatively constrained. Table 3 and Table 4 show the anomalies in the frequency of occurrence and persistence of the CTs that were highlighted from the correlation analysis in Table 1, during strong ENSO events and ENSO neutral phase. It can be observed that generally, there is no definite pattern on how ENSO influences the frequency distribution of the CTs. This fits the idea that ENSO is not always well related to rainfall variability in southern Africa. However, from Table 3, there is the indication that strong El Niño events might enhance the frequency of occurrence of CT1 + . Also, generally, during strong El Niño (La Niña) events, the frequency of occurrence and persistence of CT3 + are likely to be constrained (enhanced). The persistence of CT5 − is equally likely to be enhanced (constrained) during strong El Niño (La Niña) events. It can be observed  It was also found that during the 1992 El Niño event, regardless of the increase in the persistence of CT5 − , its frequency of occurrence and circulation features seem to be constrained. Physically, this was examined in Fig. 6 based on anomalies in moisture fluxes at 850 hPa during the periods CT5 − occurred in 1992 for the JFM and MAM seasons, and the resultant effect on precipitation anomaly in South Africa in the 1992 JFM and MAM months. The anomaly in moisture flux is calculated as the difference between the JFM climatology of CT5 − from 1979-2020 and the JFM climatology of CT5 − during 1992 (the same holds for MAM). It can be seen that unlike the climatology of CT5 − in Fig. 5, in Fig. 6A, southeasterly fluxes in the southwest Indian Ocean tend to be southerly and southwest during the 1992 El Niño season. Thus, under CT5 − , onshore moisture advection by southeast winds onto the South Africa Plateau was shows that during these seasons, a significant rainfall deficit (i.e., < 50 mm/season) occurred in the eastern and northeastern parts of South Africa. The western regions were relatively less dry; this might be attributed to the northerly wind that is evident in parts of the western regions (c.f. Fig. 6A). Northerly winds are also evident during CT1 + and CT5 + (Fig. 5) that the 1992 El Niño event enhanced their frequency of occurrence. From Fig. 7, the wet patterns related to ENSO recorded deficit in JFM seasonal rainfall totals averaged over South Africa during the 1992 El Niño event.

The influence of radiative heating on the CT significantly related to El Niño under the ssp585 scenario
Using two GCMs that simulated well the mean shape of CT5 − (Fig. 8) that is significantly related to variations of ENSO, the influence of radiative heating on the spatial pressure field of CT5 − was investigated under the ssp585 scenario.
Though the changes as shown in Fig. 9 is subjected to intermodel uncertainty, there is a model agreement on (i) pressure  southeasterly moisture fluxes are advected onto the South Africa plateau and this will, in turn, decrease the rate at which the boundary layer is moistened relative to the ambient atmosphere, resulting to the possible weakening of vertical motion.

Discussion
According to Gray and Sheaffer (1991), teleconnections can alter regional atmospheric circulation patterns. In the regional context of southern Africa, studies by Cook (2000), Reason and Jagadheesha (2005), Hoell et al. (2015), and Pinault (2018) indicated that El Niño signal can modulate the regional fluxes of moisture advected onto the southern Africa plateau. This can be related both to the weakening of the Mascarene high, and the Rossby waves excited due to changes in tropical Indo-Pacific convection. Klein et al. (1999) also noted that anomalous warming (enhanced cyclonic activity) of the tropical Indian Ocean (20° S-20° N) are among the El Niño signals in the Indian Ocean.
Using the NCEP-NCAR reanalysis data set, this study examined CTs in Africa south of the equator that can be modulated by the ENSO signal. It should be equally noted that in previous classifications (e.g., Ibebuchi 2021a), the output of the CT classification in this work was reproduced using ERA5 reanalysis (Hersbach et al. 2020) which is closer to reality (e.g., Gleixner et al. 2020). Thus, both data sets represent well the underlying physics in atmospheric circulations in Africa south of the equator.
The synoptic features of CT5 − found in this study to be significantly related to El Niño are in line with the synoptic conditions during El Niño from the aforementioned works of literature. These include (i) weaker state of the south Indian Ocean high-pressure and (ii) enhanced cyclonic activity in the (tropical) south Indian Ocean. CT5 − appears to be the central CT through which El Niño controls other patterns of atmospheric circulation in southern Africa. During strong El Niño events, its amplitude and persistence might be expected to increase and this implies a further weakening of the Mascarene high and the associating fluxes of southeasterlies. These results in the enhancement of CT1 + which is the overall climatology of atmospheric circulation in southern Africa associated with anticyclonic activity over the southern Africa landmasses. The result is in line with the report by Hoell et al. (2015) that El Niño weakens vertical motion over southern Africa.
During the 1992 strong El Niño event, it was found that the El Niño signal could have caused a robust decrease in the frequency of occurrence of CT3 + which is a wet CT, and the austral summer climatology of atmospheric circulation in southern Africa that favors easterly moisture fluxes from the southwest Indian Ocean. Also, during the 1992 El Niño event, the anomalies in CT5 − revealed that southeast winds were adjusted to southerly and southwest winds. Thus, less moisture driven by the South Atlantic Ocean high-pressure penetrated the eastern regions of South Africa. This could have contributed to the anomalous dry conditions (mainly in the eastern parts) during the 1992 wet seasons.
Relative to 1992, stronger El Niño events were equally recorded in 2016, 1998, and 1983. According to Richard (2005), 1983 and 1992 are among the years with the most severe droughts at the 6-month time scale for the summer rainfall region of South Africa. During 1983, a similar scenario as in 1992 was observed which included a decrease (increase) in the frequency of occurrence and persistence of CT3 + (CT1 +). However, despite the strong El Niño signal during 1998, a robust decrease (increase) in the frequency of occurrence and persistence of CT3 + (CT1 +) was not observed, which could have contributed to why severe drought of a similar category as of the 1992 episode, failed to materialize (Lyon and Manson 2007). A similar argument can be inferred during the 2016 event. To this end, it should be noted that there have been other drought years in parts of South Africa (e.g., Western Cape) not completely linked to ENSO. Thus, ENSO is not the only climatic driver that can trigger hydrological extremes in South Africa. Variations of the SAM (e.g., Ibebuchi 2021a) and the Subtropical Indian Ocean dipole, coupled with the impact of anthropogenic climate change on atmospheric circulations, can equally be associated with drought episodes in different parts of South Africa.

Conclusions
Using obliquely rotated T-mode PCA applied to SLP, CTs were classified in Africa south of the equator. CT5 − associated with a weaker state of the south India Ocean high-pressure and enhanced cyclonic activity at the tropical southwest Indian Ocean was found in this paper to be statistically related to El Niño. El Niño enhances the amplitude of CT5 − and this implies weakening of atmospheric conditions associated with stronger anti-cyclonic circulation at the western branch of the Mascarene high, and enhancement of atmospheric conditions that favor westerly winds across the Agulhas current. During the 1992 El Niño event, drought conditions in South Africa might have amplified due to (i) robust decrease in the frequency of occurrence of CT3 + that favors southeasterly moisture fluxes, advected by the Mascarene high, onto southern Africa, during austral summer (ii) and the modulation of onshore moisture advection by southeast winds, under the activity of the South Atlantic Ocean high-pressure, ridging south of South Africa, under CT5 −. The latter scenario could be a result of the Rossby waves excited in the tropical Indo-Pacific.
Under the ssp585 scenario, from the analyzed climate models, the impact of radiative heating, on CT5 − suggests an increase in subsidence in the eastern parts of South Africa.

Figures 10, 11
Fig. 10 Probability of occurrence of the CTs in Fig. 3. The probabilities do not add up to 100% since the classification allows overlapping of the classified variables given the fuzzy/continuous nature of atmospheric circulation patterns. CT1 + , CT2 − , CT3 + , and CT4 − that are large classes are designated as the dominant or actual occurring states of the atmosphere ing the NCEP-NCAR reanalysis data sets and the CPC gridded precipitation data used in this work. I acknowledge the World Climate Research Programme's Working Group on Coupled Modelling, which is responsible for CMIP, and I thank the climate modeling groups for producing and making available the MPI-ESM and EC-EARTH CMIP6 simulations used in this work.
Funding Open Access funding enabled and organized by Projekt DEAL.
Author contribution Work was designed and executed by Chibuike Chiedozie Ibebuchi.
Code availability R studio was used for coding the methods as described in the methodology section. The codes are base packages in R (e.g., the PROMAX routine).