Seasonal Rainfall Variability in Ethiopia and Its Long-Term Link to Global Sea Surface Temperatures

: Investigating the influence of sea surface temperatures (SSTs) on seasonal rainfall is a crucial factor for managing Ethiopian water resources. For this purpose, SST and rainfall data were used to study a wide range of inhomogeneous areas in Ethiopia with uneven distribution of rainfall for both summer (1951–2015) and spring (1951–2000) seasons. Firstly, a preliminary subdivision of rainfall grid points into zones was applied depending on spatial homogeneity and seasonality of rainfall. This introduced new clusters, including nine zones for summer rainfall peak (July/August) and five zones for spring rainfall peak (April/May). Afterward, the time series for each zone was derived by calculating the rainfall averaged over grid points within the zone. Secondly, the oceanic regions that significantly correlated with the Ethiopian rainfall were identified through cross-correlations between rainfalls averaged over every homogeneous zone and the monthly averaged SST. For summer rainfall as a main rainy season, the results indicated that the Gulf of Guinea and southern Pacific Ocean had a significant influence on rainfall zones at a lag time of 5–6 and 6–7 months. Besides, for summer rainfall zones 8 and 9 at lag time 5–6 months, the common SST regions of the southern Pacific Ocean showed the opposite sense of positive and negative correlations. Thus, the difference in SSTs between the two regions was more strongly correlated ( r ≥ 0.46) with summer rainfall in both zones than others. For the spring season, the results indicated that SST of the northern Atlantic Ocean had a strong influence on spring rainfall zones (3 and 5) at a lag time 6–7 months, as indicated by a significant correlation ( r ≥ − 0.40). Therefore, this study suggests that SSTs of southern Pacific and northern Atlantic oceans can be used as effective inputs for prediction models of Ethiopian summer and spring rainfalls, respectively.


Introduction
Globally, rainfall variability across time and space affects all aspects of human activity, especially agricultural economies and social activities. In particular, rainfall is the most significant meteorological parameter in Ethiopia, as approximately 85% of the Ethiopian labor force is employed in rain-fed agriculture which highly depends on low or high amounts of rainfall availability vital for crop production [1]. Moreover, the Ethiopian Electric and Power Corporation (EEPCo) declared that the local hydropower contributed to more than 98% of the energy production in 2005, 2006, and 2007 [2].
The geographical location and topographic complexity of Ethiopia produce high rainfall variability in the region across space and time [3]. Spatial variations include the rainfall seasonal cycle, amount, onset, and cessation times, and length of growing season [3,4], but sometimes, rainfall can be temporally varied from days to decades in terms of the direction and magnitude of rainfall trends over regions and seasons [5,6]. Moreover, these spatiotemporal variations in rainfall are attributed to, in Ethiopia, the variable altitudes [3] and over the Pacific, Atlantic, and Indian oceans [7][8][9], the variable sea surface temperatures (SSTs), including the interannual and interseasonal variations in the strength of monsoon over the Arabian Peninsula [4,7,8]. Diro et al. in 2011 [2] described the highlands of Ethiopia to exhibit three cycles of seasonal rainfall: spring mid-rainy season (February-May), summer rainy season (June-September), and dry season (October-January), locally known as Belg, Kiremt, and Bega seasons, respectively. In addition, Korecha and Barnston [9] investigated the climatological factors influencing rainfall over the region: (i) migration of the intertropical convergence zone (ITCZ) meridionally; (ii) warm lows formed over the Arabian landmasses and Sahara; (iii) formation of sub-tropical high pressure over the Azores; (iv) flow of cross-equatorial moisture from the southern, central, and equatorial parts of the Indian Ocean, tropical Africa, and Atlantic Ocean, respectively; (v) flow of the upper-level tropical easterly jet over Ethiopia; and (vi) low-level Somali jet.
Different studies have shown that SST can provide crucial predictive information regarding hydrological variability in different regions of the globe [10][11][12]. For instance, most previous studies [1,2,4,[7][8][9][10]13,14] have focused only on addressing the linkage between summer rains (June-September) that contribute to about 65%-95% of the total annual rainfall in Ethiopia [15] and global SSTs. Generally, the majority of previous studies relied on SST as a main predictor of Ethiopian seasonal rainfall. For example, El Niño/Southern Oscillation (ENSO) as a large-scale phenomenon has a significant influence on the interannual variability of Ethiopian rains [13,16,17], while the effect of SSTs of Pacific Ocean, Indian Ocean, and Gulf of Guinea on the summer rainfall of Ethiopia was indicated in [2,18,19]. Meanwhile, low-level winds from the Atlantic and Indian oceans influenced rainfall of summer (Kiremt) season [7,8,14]. A few of these studies indicated some correlations between SSTs over Gulf of Guinea and southern Atlantic Ocean and summer rainfall over Ethiopia [2,4,7,8]. Moreover, associations between the spring rains (March-May) as a main rainfall season over the southern part of Ethiopia and SST are much less investigated. For example, two studies only [20,21] used SST anomalies from the Pacific, Indian, and Atlantic oceans to predict spring rainfall.
Generally, the majority of previous studies imply that the teleconnections between SST and Ethiopian rainfall are complicated on both spatial and temporal scales and still not well investigated. Therefore, the objective of this study is directed toward more understanding the characteristics of interannual variability of rainfall in the country and investigating clearly the influence of global SSTs on Ethiopian rainfall peaks at various regions and seasons (summer and spring), to reinforce the skill of rainfall predictions that would be valuable for operations of reservoirs [22], assessment and allocation of water resources [23,24], and for mitigation of disasters related to rainfall such as flooding and drought [25]. Also, better understanding of the rainfall variability plays a key role in several applications including hydrological analysis [26][27][28][29][30][31] and soil erosion risk assessment [32][33][34][35][36].
For these aforementioned purposes, based on the Climate Research Unit (CRU) database of rainfall and SST data over the period 1951-2015, it would be valuable to divide Ethiopia into homogenous rainfall zones due to high spatial variability of rainfall over the country. Different previous studies have relied on principal component analysis for zonation [37,38], but this technique delivered no significant results in case of complex variability of rainfall, as in Ethiopia [39]. For Ethiopia, there is no previous study for comprehensive regional zoning depending on spatial patterns of teleconnections for various seasons of rainfall. Only the summer rainfall seasons (June-September) were regionally classified in previous studies [1,2,13] based on interannual correlations of rainfall amounts and homogenous annual cycles of rainfall from rain gauges.
In this regard, this study considered two main seasons (summer and spring) of rainfall over Ethiopia to define the homogenous rainfall zones based on SST to rainfall teleconnections besides other factors mentioned in a previous study [13] such as standardized rainfall cross-correlations, annual cycles of rainfall, and geographical vicinity. Moreover, correlation analysis was applied to investigate the links between global SSTs in different oceanic regions and Ethiopian summer and spring rainfalls. SST is a main indicator because of its slowly varying rate of change and the high ocean-atmosphere coupling [2]. Finally, these significantly correlated SST areas with Ethiopian rainfall are suggested to sufficiently enhance the seasonal rainfall predictions over Ethiopia in further studies.

Study Area
The entire Ethiopian domain has an overall area of 1.13 × 10 6 km 2 and spans from latitude 3° to 15° N and longitude 33° to 48° E. The region is distinguished by highly complex topography in the northern and central highlands and by the lowland of the East African Rift Valley in the traversed northeast-southwest portion. Altitudes range between hundreds of meters below sea level in the northeastern regions to thousands of meters above sea level in the northern highlands ( Figure 1). Based on seasonal rainfall cycle, Ethiopia can be divided into three main regions: (i) northern and central western areas with one rainy season whose peak occurs in July/August; (ii) southern part with two seasons of short rainfall (September-November) and long rainfall (March-May); and (iii) eastern and central parts of the country with two rainy periods called spring season (February-May) and summer rainy season (June-September) [40]. The interannual rainfall in Ethiopia is highly variable during all seasons. Previous studies [4,7,8,13] observed that this variability is basically influenced by large-scale phenomena (e.g., ENSO), while other studies based on modeling suggested SSTs of Indian and Pacific oceans influence Ethiopian rainfall [18]. Investigating this large-scale relation between SST of different oceanic regions and rainfall over Ethiopia is defined as teleconnection. Definitely, rainfall is related historically to different climate parameters. However, SST is the prominent indicator due to its slow rate of change that helps for long-term investigations [2]. Therefore, variable teleconnections motivate more accurate regional-scale prediction models for seasonal rainfalls over the country.
The northwestern and central highlands of Ethiopia are considered main water resources not only for Ethiopia, but also for Sudan and Egypt where the Blue Nile basin contributes more than 50% of the Nile River flow [41]. Furthermore, at highlands (>1500 m), 90% of Ethiopian people reside and depend mainly on rain-fed agriculture, with a high percentage of labor forces (89%) involved in smallholder agriculture, while the remaining population lives at lowlands (<1500 m) surrounding the highlands [42]. Crop production in Ethiopia has two cycles: (i) long-cycle crops (e.g., sorghum) are cultivated during summer and spring seasons and account for half of the national production, and (ii) short-cycle crops (e.g., barley, wheat, teff) are grown during summer and spring rainfall seasons and are responsible for 40%-54% and 5%-10% of the Ethiopian crop production, respectively [43].
Given such a high dependence on rain-fed agriculture in Ethiopia, it becomes a vital matter to accurately identify rainfall trends and the influence of SSTs on seasonal rainfall for enhancing the seasonal prediction models within Ethiopia on a large scale, as these are expected to deliver various benefits for the development and achievement of environmental sectors planning to meet the food and water demand of its people.

Rainfall Data
In many Ethiopian regions, existing rain gauges do not provide timely or sufficient rainfall pattern data due to the scattering of observations, uneven allocation, and data gaps in the original dataset [44,45]. These common problems cause inhomogeneity in the climate data series, causing an abrupt change in the average values and in series trend [46]. Additionally, gridded rainfall data have been widely used for different hydro-climatological analyses and climate variability studies [10,[47][48][49][50][51].
This study relied particularly on the monthly rainfall data of about 381 rainfall grid points obtained from the CRU database (https://www.cru.uea.ac.uk/data). The data range covered 65 years of summer rainfall (1951-2015) and 50 years of spring rains (1951-2000) over Ethiopia. There are some previous studies [52][53][54] that also showed the availability of using CRU data.

SST Data
Global SST data (with resolution of 1° × 1°) used in this study covered the period (1950-2015), i.e., 1 year prior to rainfall data considering lag time. These data were provided by the Hadley Centre Global Sea Ice and SST (HadISST) version 1.1 by the British Atmospheric Data Center and described by Rayner et al. [55]. SST represents a key indicator due to its slow rate of change, length, and intensive ocean-atmospheric coupling [2]. Many previous studies have focused on identifying the link between SSTs at various oceanic regions and Ethiopian rainfall [1,2,9,10,21,56].

Zoning of Rainfall Grid Points
Dividing Ethiopia into homogenous rainfall zones was a crucial process due to the highly spatial variation in rainfall. Different zoning methods have been described in many studies [13,[37][38][39]57]. For this study though, a new division scheme was introduced based on the abovementioned rainfall data, as the standard classification of Ethiopia into three zones seemed inaccurate. The methodology followed herein was that of Gissila et al. [13]. Briefly, in the first procedure, cross-correlation was applied to identify the relation between monthly averaged rainfalls for various grids on each zone. In the second step, the mean correlation coefficient for rainfall zones was derived. Thus, grids showing similarity in annual rainfall cycles could be categorized. Subsequently, the boundaries of each zone were adjusted through the mean intergrid correlations between and within zones to ensure zonal homogeneity.

Standardization
Standardization was conducted for each zone to detect statistically considerable changes in rainfall, where the dependent variable of annual zone rainfall was regressed on the independent variable of time. Here, the data used covered 65 (1951-2015) and 50 (1951-2000) years for summer and spring seasons, respectively. Next, the temporal tendency was examined for statistical significance of the Pearson correlation coefficient (r) following Student's t-test as explained in detail by Hirsch et al. [58]. The statistical significance of the trend was evaluated (at 0.01 significance level) against the null hypothesis of no trend for the series of data.

Correlation Analysis
The cross-correlation between SST at each grid point and averaged rainfall over every homogenous zone was used to identify the potential coverabilities or teleconnections between global SSTs and averaged rainfall peaks of summer July/August (JA) and spring April/May (AM) over Ethiopia. Here, the Pearson product-moment correlation coefficient (r) was a crucial element for evaluating the teleconnection (remote linkage) between rainfall and SST at different significance levels. The value of (r) varied from +1 (perfectly positive correlation) to −1 (perfectly negative correlation), while zero indicated no link. Student's t-test was applied against the null hypothesis of no correlation for the assessment of statistical significance [59]. Estimation of the statistical significance of a variable in time series has been explained in different previous studies [58,60], whereas the significance of lag-1 correlation coefficient was firstly introduced by Anderson in 1942 [61], followed by Kendall and Stuart in 1948 [62]. Higher-order lags were used when the first-order correlation insufficiently depicted the serial dependence [63]. For this purpose, lagged (r) was calculated for lags from 0 to 9 months, i.e., calculating (r) until October/November and July/August of the previous year to summer and spring peaks, respectively, for identification of lags that show significant correlations [11,64]. The normal distribution of SST data followed the distribution explained in previous studies [11,65]. Table 1 shows the calculated correlation coefficients |r| corresponding to statistical significance levels. For this study, the concentration was on the significance at 1% level (two-tailed) with the calculated correlation coefficients of 0.32 and 0.36 for summer and spring rainfalls, respectively.  Figure 2 displays the country of Ethiopia being divided into two main subregions: (i) northern and central regions with nine zones (1-9) for summer rainfall and (ii) southern region with five zones (1-5) for spring rains. The average cross-correlations between and within zones were tested for various times for the optimal selection of a proper number of zones, as shown in Tables 2 and 3, with a reasonable homogeneity delineating a new classification for Ethiopia into nine summer and five spring zones. In case of high spatial and temporal variability of Ethiopian rainfall, it would be crucial to delineate boundaries and define the country into homogenous zones, as the optimal zonation would be valuable for hydrologic modeling, prediction, ecological, and climate classification, or for any other analysis [66,67]. For this purpose, different research studies have been undertaken. For example, Korecha and Sorteberg [68] indicated that the National Meteorological Agency of Ethiopia divides the region into eight homogenous zones depending on the main systems of atmosphericoceanic circulation and rain production. Moreover, Gissila et al. [13] grouped Ethiopia into four clusters based on seasonal rainfall cycles and homogeneity of interannual rainfalls at 24 rain gauges. Despite following the same methodology of the latter study, both this study and Diro et al.'s [57] produced different results. Here, the country was divided into nine summer peak (JA) zones and five spring peak (AM) zones based on CRU seasonal rainfall data covering Ethiopia (381 grid points), whereas Diro et al. [57] delineated five homogenous summer rainfall regions (June-September) and six spring zones (February-May) depending on the data from 33 rainfall stations and considering the interannual variation of seasonal rainfalls. Additionally, a number of recent studies carried out [21,52,69] obtained different zoning outcomes. Therefore, the classification of Ethiopia into homogenous zones could be different from one study to another, because of the high variability in large-scale forces producing rainfall over the country.    Figure 3 shows the mean monthly rainfall values calculated for each zone. Apparently, the northern and central Ethiopian regions displayed a regime of monomodal summer rainfall cycle (June-September) with rainfall peaks in July/August (hereafter JA), which is consistent with that reported by Diro et al. [2]. In contrast, the southern region of the country experienced bimodal rainfall cycles for spring or long rainfall (March-May) with rainfall peaks in April/May (hereafter AM) and another short rainfall cycle from (October-November).  Figure 4a shows a time series of standardized summer rainfall variability for each zone from 1951 to 2015 with a considerable variation between wet and dry years from one zone to another. The northern and central rainfall zones displayed a successive time period of 2-4 years of dry and wet years without any distinct trend, in addition to three records of drought years in 1965, 1984, and 2002. Moreover, the magnitude of the series gradually decreased for all zones, except for the central zones (7 and 8). Figure 4b indicates the time series of standardized spring rainfall variability for each zone from 1951 to 2000, where the eastern and southern zones indicated a continual below-average amount of rainfall from 1952 to 1956, but recovered in the late 1950s. Moreover, the drought years appeared in 1974, 1980, and 1984 with a prevalent drier condition from the beginning of the 1990s. Both figures also confirmed Ethiopia's significant interannual rainfall variability with the year 1984 being the major year of drought that covered the whole of Ethiopia [6], in which the standard deviations varied between 0.7-below-average in the southern zones and greater-than-2 below-average in the northern zones. Moreover, the country experienced abundant amounts of rainfall for 4 years (1951, 1961, 1967, and 1982) of the given 65 years for the study. These results were consistent with Seleshi and Zanke [6] for their analysis of the annual rainfall totals and annual rainy days over the country. The Ethiopian rainfall variability has been described in numerous studies, depending on the leading mechanisms of the ITCZ migration, influences of topography, and various interactions of hydroclimate system on regional scales [6,70].

Cross-Correlation between Summer Rainfall (July/August Peak) and Oceanic SSTs
Firstly, for the 65 years of summer rainfall , the correlation coefficients ±0.24, ±0.32, and ±0.40 corresponding to significance levels of 0.05, 0.01, and 0.001 were calculated using two-tailed t-test. As shown in Figure 5, the common summer rainfall zones having the same significantly correlated SST regions were statistically extracted and represented at significance levels. In particular, this study relied on the 0.01 significance level with r = ±0.32, showing a significantly strong correlation at a lag time of 5-8 months between the July/August rainfall and SST at the nine zones. Positive and negative correlations were represented in red and blue colors, respectively ( Figure 5). From the standpoint of forecasting, possible correlations of SST with rainfall at longer lag times are potentially effective. For instance, as shown in the left side of Figure 6, at 5-6 months lag period, Gulf of Guinea SST showed a significantly positive correlation with July/August rainfall in the central zones (6 and 9), while two regions in the Pacific Ocean exhibited the strongest teleconnection (positive and negative) between SST and rainfall at zones 8 and 9, due to the well-known Pacific events of ENSO. A recent investigation of Gleixner et al. [15] revealed that half of Ethiopian summer rainfall is affected by the variation of equatorial Pacific SST. Additionally, various previous studies indicated the teleconnection between SST of the Pacific and Ethiopian summer rainfall [1,4,9,13,19,71] and excess of rainfall during spring season. At 6-7 months lag period, the central summer rainfall zones 4 and 6 indicated a significantly positive correlation with both Gulf of Guinea and a small area over the southeastern part in the Pacific. Diro et al. [2] identified the effect of the Pacific SST on rainfall, besides the influence of the Gulf of Guinea, over different parts of Ethiopia. The strongest influence of southern Indian Ocean SST on Ethiopian rainfall due to Ethiopia's location was indicated at rainfall zones 5 and 7 (right side of Figure 6); however, the effect of the Indian Ocean SST was disregarded, as the study's objective was directed toward the remote influence (teleconnection) of oceanic SST on Ethiopian rainfall. For the northern summer zones 1 and 2, the northwest Pacific showed insignificant influence on rainfall (negative correlation) at a lag time of 7-8 months (Figure 7). In 2005, Segele and Lamb [4] indicated that upper-level features such as the African easterly jet (AEJ), and tropical easterly jet (TEJ) affect summer rainfall over Ethiopia based on their strength and location. Moreover, moisture flux from the Atlantic and Indian oceans and as a lower-level circulation affected by the highs of Mascarene and St. Helena affects summer rainfall significantly [40]. Generally, various previous studies [1,2,14,72,73] identified the teleconnection between SST of the Indian Ocean and Ethiopian rainfall, also they suggested that sea surface temperature anomalies (SSTA) of eastern part of Pacific Ocean may affect summer rainfall as a main rainy season. These positive SSTAs are associated with deficits in summer rainfall and excess of rainfall during spring season.

Cross-Correlation between Spring Rainfall (April/May Peak) and Oceanic SSTs
As shown in Figure 8, the common spring rainfall zones having the same significantly correlated SST regions were statistically extracted and represented at significance levels. In particular, this study relied on the 0.01 significance level with r = ±0.36, showing a significantly strong correlation at the lag period of 5-8 months between April/May rainfall and SST at five zones. From the standpoint of prediction, possible correlations of SST with rainfall at longer lag times are potentially effective. For instance, at lag time of 6-7 months, the northern Atlantic Ocean SST showed a significantly negative correlation with spring rainfall zones 3 and 5 representing southern and southeastern regions of Ethiopia (left side of Figure 9). Moreover, at 7-8 months, two regions at the southern Indian Ocean exhibited a strong positive influence on spring rainfall due to Ethiopia's location for the southeastern zones 4 and 5 (right side of Figure 9). The passage of ITCZ controls the shape of seasonal rainfall cycle. Annually, it oscillates from 15° N in July to 15° S in January and back again causing monomodal and bimodal rainfall patterns over northern and southern parts of Ethiopia, respectively [72]. In addition to the difference between the land surface and Indian Ocean in the heat capacity inducing ITCZ's meridional arm. The latter produces rainfall during February and March over the southwestern part of Ethiopia [40]. Additionally, the strength and location of upper-level feature of the subtropical westerly jet (STWJ) affects spring rainfall over the country [4]. Generally, SST and spring rainfall teleconnections are spatially and temporally complicated, where Diro et al. 's [20] was the only previous study to investigate the influence of SSTs of the Pacific, Atlantic, and Indian oceans on spring rains (March-May), depicting the main rainfall season over the southern part of Ethiopia. Therefore, the results in this study provide a more detailed understanding of such teleconnections that would be valuable for enhancing rainfall prediction system over Ethiopia.
Finally, to visualize the SST regions showing significant correlation with seasonal rainfall over Ethiopia, the correlation coefficients between rainfall time series and SSTs were derived for the aforementioned summer and spring zones with significant correlations (Figures 10 and 11). These regions could be used as inputs for a forecasting model for further studies at a lead time before the rainfall season. For the summer season, a significant correlation was indicated (r ≥ ±0.37) at a lag time of 5-6 and 6-7 months prior to the rainy season. In particular, summer zones (8 and 9) at 5-6 months lag period indicated opposite sense of positive and negative correlations, with the difference in SSTs enhancing the correlation. Thus, both zones (8 and 9) exhibited the strongest correlation (r ≥ 0.46) between summer rainfall and SST among all other summer zones ( Figure 10). On the other side, the spring rainfall zones (3 and 5) also showed significant correlation (r ≥ −0.40) at 6-7 months lag time, representing the influence of northern Atlantic Ocean SST, while the southern Indian Ocean SST exhibited the strongest correlation (r ≥ 0.5) at 7-8 months of lag time ( Figure 11).
These results suggest that SSTs of the southern Pacific and northern Atlantic oceans could be effective as inputs for prediction models of Ethiopian summer and spring rainfalls, respectively, at a lead time before these seasons.

Conclusions
Ethiopian rainfall is characterized by both seasonal and interannual spatial variabilities. Thus, understanding the teleconnection (remote link) between Ethiopian rainfall and oceanic SSTs is a key factor for establishing future rainfall prediction models and further managing the country's water resources. For this reason, this study was directed toward analyzing rainfall time series on both seasonal and annual scales, investigating their potential trends, and detecting their significance over Ethiopia, along with investigating the oceanic SST regions that significantly influence Ethiopian rainfall in order to detect the potential SST regions that could be valuable as input data for establishing forecasting models in further studies. Firstly, in this study, rainfall data for summer  and spring  seasons were analyzed for identification of the rainfall peaks of July/August and April/May for summer and spring seasons, respectively. Secondly, preliminary clustering or dividing of Ethiopia into homogenous rainfall zones was applied, depending on seasonality and spatial homogeneity of rainfall data, which introduced a new number of rainfall zones for Ethiopia with 14 homogenous rainfall zones: nine rainfall zones for summer peak July/August and five zones for spring peak April/May. Thirdly, cross-correlations between SST and rainfall peaks at the 0.01 significance level revealed oceanic regions that significantly influence rainfall over Ethiopia at different lag times. For summer rainfall (main rainy season), the influence of Gulf of Guinea (positive correlation) and two regions in the southern Pacific Ocean (positive and negative correlations) was significantly detected at a lag time of 5-6 months.
In addition, SST of northwestern Pacific showed a weak negative influence on summer northern rainfall zones (1 and 2) at a lag time 7-8 months, whereas SST of the Indian Ocean indicated a significantly negative effect on summer eastern and western zones (5 and 7) at a lag time of 6-7 months, being close to Ethiopia. For spring rainfall, the negative influence of the northern Atlantic Ocean SST was strongly detected on spring rainfall zones (3 and 5) at a lag time of 6-7 months, whereas two areas of the Indian Ocean showed a positive correlation between SSTs and spring rainfall zones (4 and 5).
Finally, correlation coefficients (r) between these significant SST oceanic regions and both summer and spring rainfalls were derived for detection of the SST regions that could be valuable as input data for prediction models at lead times before both seasons. For the summer season, at a lag time of 5-6 months, both central zones (8 and 9) exhibited the strongest positive correlation (r ≥ 0.46) between summer rainfall and southern Pacific Ocean SST, among all other summer zones due to the difference between positive and negative correlations at the two areas of the southern Pacific Ocean. Alternatively, a significantly negative correlation (r ≥ −0.40) was identified between spring rainfall zones (3 and 5) and the northern Atlantic Ocean SST at a lag time of 6-7 months before spring season.
Therefore, for later studies, it is recommended to use SST regions of the southern Pacific and northern Atlantic oceans as effective input data for forecasting models of Ethiopian summer and spring rainfall peaks at lead times of 5-6 and 6-7 months, respectively.